Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Competition-Based Model of Pheromone Component Ratio Detection in the Moth</title><author><name sortKey="Zavada, Andrei" sort="Zavada, Andrei" uniqKey="Zavada A" first="Andrei" last="Zavada">Andrei Zavada</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Buckley, Christopher L" sort="Buckley, Christopher L" uniqKey="Buckley C" first="Christopher L." last="Buckley">Christopher L. Buckley</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Martinez, Dominique" sort="Martinez, Dominique" uniqKey="Martinez D" first="Dominique" last="Martinez">Dominique Martinez</name><affiliation><nlm:aff id="aff2"><addr-line>CORTEX Team - LORIA, Vandoeuvre-lès-Nancy, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rospars, Jean Pierre" sort="Rospars, Jean Pierre" uniqKey="Rospars J" first="Jean-Pierre" last="Rospars">Jean-Pierre Rospars</name><affiliation><nlm:aff id="aff3"><addr-line>UMR 1272 Physiologie de l'insecte, INRA, Versailles, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Nowotny, Thomas" sort="Nowotny, Thomas" uniqKey="Nowotny T" first="Thomas" last="Nowotny">Thomas Nowotny</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">21373177</idno><idno type="pmc">3040183</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3040183</idno><idno type="RBID">PMC:3040183</idno><idno type="doi">10.1371/journal.pone.0016308</idno><date when="2011">2011</date><idno type="wicri:Area/Pmc/Corpus">000086</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000086</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Competition-Based Model of Pheromone Component Ratio Detection in the Moth</title><author><name sortKey="Zavada, Andrei" sort="Zavada, Andrei" uniqKey="Zavada A" first="Andrei" last="Zavada">Andrei Zavada</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Buckley, Christopher L" sort="Buckley, Christopher L" uniqKey="Buckley C" first="Christopher L." last="Buckley">Christopher L. Buckley</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author><author><name sortKey="Martinez, Dominique" sort="Martinez, Dominique" uniqKey="Martinez D" first="Dominique" last="Martinez">Dominique Martinez</name><affiliation><nlm:aff id="aff2"><addr-line>CORTEX Team - LORIA, Vandoeuvre-lès-Nancy, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rospars, Jean Pierre" sort="Rospars, Jean Pierre" uniqKey="Rospars J" first="Jean-Pierre" last="Rospars">Jean-Pierre Rospars</name><affiliation><nlm:aff id="aff3"><addr-line>UMR 1272 Physiologie de l'insecte, INRA, Versailles, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Nowotny, Thomas" sort="Nowotny, Thomas" uniqKey="Nowotny T" first="Thomas" last="Nowotny">Thomas Nowotny</name><affiliation><nlm:aff id="aff1"><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2011">2011</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>For some moth species, especially those closely interrelated and sympatric, recognizing a specific pheromone component concentration ratio is essential for males to successfully locate conspecific females. We propose and determine the properties of a minimalist competition-based feed-forward neuronal model capable of detecting a certain ratio of pheromone components independently of overall concentration. This model represents an elementary recognition unit for the ratio of binary mixtures which we propose is entirely contained in the macroglomerular complex (MGC) of the male moth. A set of such units, along with projection neurons (PNs), can provide the input to higher brain centres. We found that (1) accuracy is mainly achieved by maintaining a certain ratio of connection strengths between olfactory receptor neurons (ORN) and local neurons (LN), much less by properties of the interconnections between the competing LNs proper. An exception to this rule is that it is beneficial if connections between generalist LNs (i.e. excited by either pheromone component) and specialist LNs (i.e. excited by one component only) have the same strength as the reciprocal specialist to generalist connections. (2) successful ratio recognition is achieved using latency-to-first-spike in the LN populations which, in contrast to expectations with a population rate code, leads to a broadening of responses for higher overall concentrations consistent with experimental observations. (3) when longer durations of the competition between LNs were observed it did not lead to higher recognition accuracy.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Ljungberg, H" uniqKey="Ljungberg H">H Ljungberg</name></author><author><name sortKey="Anderson, P" uniqKey="Anderson P">P Anderson</name></author><author><name sortKey="Hansson, Bs" uniqKey="Hansson B">BS Hansson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anton, S" uniqKey="Anton S">S Anton</name></author><author><name sortKey="Hansson, Bs" uniqKey="Hansson B">BS Hansson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tumlinson, Jh" uniqKey="Tumlinson J">JH Tumlinson</name></author><author><name sortKey="Brennan, Mm" uniqKey="Brennan M">MM Brennan</name></author><author><name sortKey="Doolittle, Re" uniqKey="Doolittle R">RE Doolittle</name></author><author><name sortKey="Mitchell, Er" uniqKey="Mitchell E">ER Mitchell</name></author><author><name sortKey="Brabham, A" uniqKey="Brabham A">A Brabham</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="T Th, M" uniqKey="T Th M">M Tòth</name></author><author><name sortKey="Lofstedt, C" uniqKey="Lofstedt C">C Löfstedt</name></author><author><name sortKey="Blair, Bw" uniqKey="Blair B">BW Blair</name></author><author><name sortKey="Cabello, T" uniqKey="Cabello T">T Cabello</name></author><author><name sortKey="Farag, Ai" uniqKey="Farag A">AI Farag</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Mustaparta, H" uniqKey="Mustaparta H">H Mustaparta</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ishikawa, Y" uniqKey="Ishikawa Y">Y Ishikawa</name></author><author><name sortKey="Takanashi, T" uniqKey="Takanashi T">T Takanashi</name></author><author><name sortKey="Kim, Choong Gon" uniqKey="Kim C">Choong-gon Kim</name></author><author><name sortKey="Hoshizaki, S" uniqKey="Hoshizaki S">S Hoshizaki</name></author><author><name sortKey="Tatsuki, S" uniqKey="Tatsuki S">S Tatsuki</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Linster, C" uniqKey="Linster C">C Linster</name></author><author><name sortKey="Masson, C" uniqKey="Masson C">C Masson</name></author><author><name sortKey="Kerszberg, M" uniqKey="Kerszberg M">M Kerszberg</name></author><author><name sortKey="Personnaz, L" uniqKey="Personnaz L">L Personnaz</name></author><author><name sortKey="Dreyfus, G" uniqKey="Dreyfus G">G Dreyfus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hansson, Bs" uniqKey="Hansson B">BS Hansson</name></author><author><name sortKey="Anton, S" uniqKey="Anton S">S Anton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, W" uniqKey="Wu W">W Wu</name></author><author><name sortKey="Anton, S" uniqKey="Anton S">S Anton</name></author><author><name sortKey="Lofstedt, C" uniqKey="Lofstedt C">C Löfstedt</name></author><author><name sortKey="Hansson, Bs" uniqKey="Hansson B">BS Hansson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Linster, C" uniqKey="Linster C">C Linster</name></author><author><name sortKey="Kerszberg, M" uniqKey="Kerszberg M">M Kerszberg</name></author><author><name sortKey="Masson, C" uniqKey="Masson C">C Masson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Linster, C" uniqKey="Linster C">C Linster</name></author><author><name sortKey="Dreyfus, G" uniqKey="Dreyfus G">G Dreyfus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Heinbockel, T" uniqKey="Heinbockel T">T Heinbockel</name></author><author><name sortKey="Kloppenburg, P" uniqKey="Kloppenburg P">P Kloppenburg</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Laurent, G" uniqKey="Laurent G">G Laurent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stopfer, M" uniqKey="Stopfer M">M Stopfer</name></author><author><name sortKey="Bhagavan, S" uniqKey="Bhagavan S">S Bhagavan</name></author><author><name sortKey="Smith, Bh" uniqKey="Smith B">BH Smith</name></author><author><name sortKey="Laurent, G" uniqKey="Laurent G">G Laurent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Getz, Wm" uniqKey="Getz W">WM Getz</name></author><author><name sortKey="Lutz, A" uniqKey="Lutz A">A Lutz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bazhenov, M" uniqKey="Bazhenov M">M Bazhenov</name></author><author><name sortKey="Stopfer, M" uniqKey="Stopfer M">M Stopfer</name></author><author><name sortKey="Rabinovich, M" uniqKey="Rabinovich M">M Rabinovich</name></author><author><name sortKey="Huerta, R" uniqKey="Huerta R">R Huerta</name></author><author><name sortKey="Abarbanel, Hd" uniqKey="Abarbanel H">HD Abarbanel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bazhenov, M" uniqKey="Bazhenov M">M Bazhenov</name></author><author><name sortKey="Stopfer, M" uniqKey="Stopfer M">M Stopfer</name></author><author><name sortKey="Rabinovich, M" uniqKey="Rabinovich M">M Rabinovich</name></author><author><name sortKey="Abarbanel, Hd" uniqKey="Abarbanel H">HD Abarbanel</name></author><author><name sortKey="Sejnowski, Tj" uniqKey="Sejnowski T">TJ Sejnowski</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Laurent, G" uniqKey="Laurent G">G Laurent</name></author><author><name sortKey="Stopfer, M" uniqKey="Stopfer M">M Stopfer</name></author><author><name sortKey="Friedrich, Rw" uniqKey="Friedrich R">RW Friedrich</name></author><author><name sortKey="Rabinovich, Mi" uniqKey="Rabinovich M">MI Rabinovich</name></author><author><name sortKey="Volkovskii, A" uniqKey="Volkovskii A">A Volkovskii</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rabinovich, M" uniqKey="Rabinovich M">M Rabinovich</name></author><author><name sortKey="Volkovskii, A" uniqKey="Volkovskii A">A Volkovskii</name></author><author><name sortKey="Lecanda, P" uniqKey="Lecanda P">P Lecanda</name></author><author><name sortKey="Huerta, R" uniqKey="Huerta R">R Huerta</name></author><author><name sortKey="Abarbanel, Hd" uniqKey="Abarbanel H">HD Abarbanel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kwok, Yc" uniqKey="Kwok Y">YC Kwok</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hodgkin, Al" uniqKey="Hodgkin A">AL Hodgkin</name></author><author><name sortKey="Huxley, Af" uniqKey="Huxley A">AF Huxley</name></author><author><name sortKey="Katz, B" uniqKey="Katz B">B Katz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Traub, Rd" uniqKey="Traub R">RD Traub</name></author><author><name sortKey="Miles, R" uniqKey="Miles R">R Miles</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Destexhe, A" uniqKey="Destexhe A">A Destexhe</name></author><author><name sortKey="Mainen, Zf" uniqKey="Mainen Z">ZF Mainen</name></author><author><name sortKey="Sejnowski, Tj" uniqKey="Sejnowski T">TJ Sejnowski</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Destexhe, A" uniqKey="Destexhe A">A Destexhe</name></author><author><name sortKey="Mainen, Zf" uniqKey="Mainen Z">ZF Mainen</name></author><author><name sortKey="Sejnowski, Tj" uniqKey="Sejnowski T">TJ Sejnowski</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schneider, D" uniqKey="Schneider D">D Schneider</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kaissling, Ke" uniqKey="Kaissling K">KE Kaissling</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sanes, Jr" uniqKey="Sanes J">JR Sanes</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Matsumoto, Sg" uniqKey="Matsumoto S">SG Matsumoto</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Homberg, U" uniqKey="Homberg U">U Homberg</name></author><author><name sortKey="Montague, Ra" uniqKey="Montague R">RA Montague</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sch Fer, S" uniqKey="Sch Fer S">S Schäfer</name></author><author><name sortKey="Bicker, G" uniqKey="Bicker G">G Bicker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Waldrop, B" uniqKey="Waldrop B">B Waldrop</name></author><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Homberg, U" uniqKey="Homberg U">U Homberg</name></author><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Olsen, Sr" uniqKey="Olsen S">SR Olsen</name></author><author><name sortKey="Wilson, Ri" uniqKey="Wilson R">RI Wilson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Olsen, Sr" uniqKey="Olsen S">SR Olsen</name></author><author><name sortKey="Bhandawat, V" uniqKey="Bhandawat V">V Bhandawat</name></author><author><name sortKey="Wilson, Ri" uniqKey="Wilson R">RI Wilson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Av Ron, E" uniqKey="Av Ron E">E Av-Ron</name></author><author><name sortKey="Rospars, Jp" uniqKey="Rospars J">JP Rospars</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Av Ron, E" uniqKey="Av Ron E">E Av-Ron</name></author><author><name sortKey="Vibert, Jf" uniqKey="Vibert J">JF Vibert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Husch, A" uniqKey="Husch A">A Husch</name></author><author><name sortKey="Hess, S" uniqKey="Hess S">S Hess</name></author><author><name sortKey="Kloppenburg, P" uniqKey="Kloppenburg P">P Kloppenburg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Husch, A" uniqKey="Husch A">A Husch</name></author><author><name sortKey="Paehler, M" uniqKey="Paehler M">M Paehler</name></author><author><name sortKey="Fusca, D" uniqKey="Fusca D">D Fusca</name></author><author><name sortKey="Paeger, L" uniqKey="Paeger L">L Paeger</name></author><author><name sortKey="Kloppenburg, P" uniqKey="Kloppenburg P">P Kloppenburg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Husch, A" uniqKey="Husch A">A Husch</name></author><author><name sortKey="Paehler, M" uniqKey="Paehler M">M Paehler</name></author><author><name sortKey="Fusca, D" uniqKey="Fusca D">D Fusca</name></author><author><name sortKey="Paeger, L" uniqKey="Paeger L">L Paeger</name></author><author><name sortKey="Kloppenburg, P" uniqKey="Kloppenburg P">P Kloppenburg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Press, Wh" uniqKey="Press W">WH Press</name></author><author><name sortKey="Teukolsky, Sa" uniqKey="Teukolsky S">SA Teukolsky</name></author><author><name sortKey="Vetterling, Wt" uniqKey="Vetterling W">WT Vetterling</name></author><author><name sortKey="Flannery, Bp" uniqKey="Flannery B">BP Flannery</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vickers, Nj" uniqKey="Vickers N">NJ Vickers</name></author><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Mustaparta, H" uniqKey="Mustaparta H">H Mustaparta</name></author><author><name sortKey="Baker, Tc" uniqKey="Baker T">TC Baker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vickers, Nj" uniqKey="Vickers N">NJ Vickers</name></author><author><name sortKey="Christensen, Ta" uniqKey="Christensen T">TA Christensen</name></author><author><name sortKey="Hildebrand, Jg" uniqKey="Hildebrand J">JG Hildebrand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hansson, Bs" uniqKey="Hansson B">BS Hansson</name></author><author><name sortKey="Almaas, Tj" uniqKey="Almaas T">TJ Almaas</name></author><author><name sortKey="Anton, S" uniqKey="Anton S">S Anton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ditzen, M" uniqKey="Ditzen M">M Ditzen</name></author><author><name sortKey="Evers, Jf" uniqKey="Evers J">JF Evers</name></author><author><name sortKey="Galizia, Cg" uniqKey="Galizia C">CG Galizia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Todd, Jl" uniqKey="Todd J">JL Todd</name></author><author><name sortKey="Baker, Tc" uniqKey="Baker T">TC Baker</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">21373177</article-id><article-id pub-id-type="pmc">3040183</article-id><article-id pub-id-type="publisher-id">PONE-D-10-02316</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0016308</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Computational Biology</subject><subj-group><subject>Computational Neuroscience</subject><subj-group><subject>Sensory Systems</subject></subj-group></subj-group><subj-group><subject>Signaling Networks</subject></subj-group></subj-group><subj-group><subject>Neuroscience</subject><subj-group><subject>Sensory Systems</subject><subj-group><subject>Olfactory System</subject></subj-group></subj-group><subj-group><subject>Neural Networks</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Competition-Based Model of Pheromone Component Ratio Detection in the Moth</article-title><alt-title alt-title-type="running-head">Competition Model of Pheromone Ratio Detection</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Zavada</surname><given-names>Andrei</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Buckley</surname><given-names>Christopher L.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Martinez</surname><given-names>Dominique</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Rospars</surname><given-names>Jean-Pierre</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Nowotny</surname><given-names>Thomas</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Informatics, University of Sussex, Brighton, United Kingdom</addr-line></aff><aff id="aff2"><label>2</label><addr-line>CORTEX Team - LORIA, Vandoeuvre-lès-Nancy, France</addr-line></aff><aff id="aff3"><label>3</label><addr-line>UMR 1272 Physiologie de l'insecte, INRA, Versailles, France</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Brezina</surname><given-names>Vladimir</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1">Mount Sinai School of Medicine, United States of America</aff><author-notes><corresp id="cor1">* E-mail: <email>t.nowotny@sussex.ac.uk</email></corresp><fn fn-type="con"><p>Conceived and designed the experiments: AZ TN. Performed the experiments: AZ. Analyzed the data: AZ TN. Contributed reagents/materials/analysis tools: AZ TN. Wrote the paper: AZ TN. Designed software: AZ TN. Essential advice DM JR. Mathematical proofs: CLB.</p></fn></author-notes><pub-date pub-type="collection"><year>2011</year></pub-date><pub-date pub-type="epub"><day>16</day><month>2</month><year>2011</year></pub-date><pub-date pub-type="ecorrected"><day>22
</day><month>2</month><year>2011</year></pub-date><volume>6</volume><issue>2</issue><elocation-id>e16308</elocation-id><history><date date-type="received"><day>20</day><month>9</month><year>2010</year></date><date date-type="accepted"><day>22</day><month>12</month><year>2010</year></date></history><permissions><copyright-statement>Zavada et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</copyright-statement></permissions><abstract><p>For some moth species, especially those closely interrelated and sympatric, recognizing a specific pheromone component concentration ratio is essential for males to successfully locate conspecific females. We propose and determine the properties of a minimalist competition-based feed-forward neuronal model capable of detecting a certain ratio of pheromone components independently of overall concentration. This model represents an elementary recognition unit for the ratio of binary mixtures which we propose is entirely contained in the macroglomerular complex (MGC) of the male moth. A set of such units, along with projection neurons (PNs), can provide the input to higher brain centres. We found that (1) accuracy is mainly achieved by maintaining a certain ratio of connection strengths between olfactory receptor neurons (ORN) and local neurons (LN), much less by properties of the interconnections between the competing LNs proper. An exception to this rule is that it is beneficial if connections between generalist LNs (i.e. excited by either pheromone component) and specialist LNs (i.e. excited by one component only) have the same strength as the reciprocal specialist to generalist connections. (2) successful ratio recognition is achieved using latency-to-first-spike in the LN populations which, in contrast to expectations with a population rate code, leads to a broadening of responses for higher overall concentrations consistent with experimental observations. (3) when longer durations of the competition between LNs were observed it did not lead to higher recognition accuracy.</p></abstract><counts><page-count count="12"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>In moths, males are attracted to females by means of sexual pheromones, which are often composed of a blend of two or more components. Pheromone signalling serves as an interspecific barrier which, along with temporal and geographical separation of the populations, arises from and facilitates the process of speciation. In some species, <italic>e.g.</italic>, <italic>Spodoptera littoralis</italic>, the secondary component modulates the behavioural response to the primary component but is not strictly necessary to attract the male <xref ref-type="bibr" rid="pone.0016308-Ljungberg1">[1]</xref>, <xref ref-type="bibr" rid="pone.0016308-Anton1">[2]</xref>. Other species maintain a blend of two (<italic>e.g.</italic>, <italic>Manduca sexta</italic>: <xref ref-type="bibr" rid="pone.0016308-Tumlinson1">[3]</xref>) or more components, as in <italic>Agrotis segetum</italic> <xref ref-type="bibr" rid="pone.0016308-Tth1">[4]</xref>, <italic>Heliothis/Helicoverpa</italic> species <xref ref-type="bibr" rid="pone.0016308-Christensen1">[5]</xref>, <xref ref-type="bibr" rid="pone.0016308-Christensen2">[6]</xref> and <italic>Ostrinia</italic> species <xref ref-type="bibr" rid="pone.0016308-Ishikawa1">[7]</xref>. In the latter two cases, different species were shown to use subsets of the same group of chemical compounds in their pheromone blend, but in different combinations, roles and/or concentration ratios. In these cases, in order to be effective in the wild, relevant components must all be present in the pheromone plume in certain proportions, and recognising the correct ratio of pheromone component concentrations becomes critical for the mating success of the male.</p><p>Previous modelling work of the pheromone subsystem goes back to early 1990's. Linster <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Linster1">[8]</xref> were successful in replicating the various different response patterns observed experimentally in <italic>Manduca sexta</italic> projection neurons (PNs) by varying the number of excitatory and inhibitory local neurons (LNs) in an abstract model of the male moth Macroglomerular Complex (MGC) where all connections between pheromone-sensitive ORNs, LNs and PNs are made. They, however, refrained from identifying patterns which may indicate successful ratio discrimination, arguing that “it is not known whether the detection of a precise ratio is achieved at the level of the glomerulus or at higher olfactory centers”. Although whether the ratio detection occurs at the MGC or at a later stage (<italic>e.g.</italic>, in the mushroom bodies), is generally species-dependent <xref ref-type="bibr" rid="pone.0016308-Hansson1">[9]</xref>, there is evidence <xref ref-type="bibr" rid="pone.0016308-Wu1">[10]</xref> that, in <italic>Agrotis segetum</italic>, PNs relay both component-specific as well as blend-specific signals to higher brain centres. This suggests that at least partially, ratio recognition takes place within the LN-PN network in the MGC.</p><p>Extending their earlier work, Linster <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Linster2">[11]</xref> considered a binary blend of major and minor components and suggested that component ratio detection is signalled by the presence of oscillations rather than the overall firing rate in the PN population in the MGC. In another work they refined their models by introducing several physiological constraints <xref ref-type="bibr" rid="pone.0016308-Linster3">[12]</xref>. While oscillations have been observed in the MGC of <italic>M. sexta</italic> <xref ref-type="bibr" rid="pone.0016308-Heinbockel1">[13]</xref> as well as in the general olfactory subsystem of locusts <xref ref-type="bibr" rid="pone.0016308-Laurent1">[14]</xref>, recent work in our group indicates that oscillations do not seem to be present in the pheromone subsystem of <italic>Agrotis ipsilon</italic> (unpublished data; A.Chaffiol, pers. comm.). Similarly, while oscillations have been found in the locust, suppressing them with picrotoxin injections does not fully abolish odour recognition <xref ref-type="bibr" rid="pone.0016308-Stopfer1">[15]</xref>. It also remains open how the downstream brain centres would estimate the spectral power of the PN output if the presence or absence of oscillations was the main carrier of information.</p><p>Similarly to Linster and colleagues, Getz & Lutz <xref ref-type="bibr" rid="pone.0016308-Getz1">[16]</xref> used a detailed recurrent network model of LN and PN populations to explain the emergence of odour specific response patterns that allowed odour quality recognition for varying concentrations. Subsequent works of Bazhenov <xref ref-type="bibr" rid="pone.0016308-Bazhenov1">[17]</xref>, <xref ref-type="bibr" rid="pone.0016308-Bazhenov2">[18]</xref> put forward a model for the locust antennal lobe including both fast oscillations and slow patterning.</p><p>The slow patterning of odour responses and the spatio-temporal coding were incorporated into the theoretical framework of “winnerless competition” <xref ref-type="bibr" rid="pone.0016308-Laurent2">[19]</xref>, <xref ref-type="bibr" rid="pone.0016308-Rabinovich1">[20]</xref> to explain discrimination, and possibly learning, of a variety of odours in the general olfactory system. In this framework, the presence of synchronous oscillations is interpreted as a signal of the presence of odour input and slow patterning as an odour's distinguishing characteristic. More recently, Kwok <xref ref-type="bibr" rid="pone.0016308-Kwok1">[21]</xref> explored the behaviour of a model with direct ORN-PN and feedback LN-PN connections, and described the temporal evolution of the model state and related it to winnerless competition. Individual odours as well as odour blends in Kwok's model have distinct spatio-temporal representations. However, while the combinatorial explosion of spatio-temporal patterns is a distinctive advantage for coding the vast space of relevant odour inputs in the general olfactory system, it is difficult to identify a role for spatio-temporal patterns in the pheromone subsystem where the components are few and known beforehand, and only a certain fixed target ratio is biologically relevant.</p><sec id="s1a"><title>Competition prerequisites</title><p>In this paper we study minimalistic models of competition between LNs in the MGC and how they allow recognising a specific pheromone component ratio across a wide range of pheromone concentrations. The simplest model, “an elementary ratio recognition unit” shown in <xref ref-type="fig" rid="pone-0016308-g001">fig. 1</xref>, is an essentially feed-forward (<italic>i.e.</italic>, lacking feedback connections such as from PNs to LNs or LNs to ORNs) network that is capable of recognising a target ratio of a blend of two pheromone components. Over a range of one order of magnitude of ORN response frequency, the unit produces a discrete ‘yes/no’ signal, in the form of a steady spike train or quiescence, identifying the matching target component ratio determined by its intrinsic properties. If the ORN response is logarithmic in the concentration, this corresponds to an operating range of 7 orders of magnitude in pheromone concentration. The model uses conductance based, Hodgkin–Huxley type models for LNs and PNs <xref ref-type="bibr" rid="pone.0016308-Hodgkin1">[22]</xref>, <xref ref-type="bibr" rid="pone.0016308-Traub1">[23]</xref> and first order synapses <xref ref-type="bibr" rid="pone.0016308-Destexhe1">[24]</xref>, <xref ref-type="bibr" rid="pone.0016308-Destexhe2">[25]</xref>.</p><fig id="pone-0016308-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g001</object-id><label>Figure 1</label><caption><title>The structure of the elementary pheromone component ratio recognition unit in the MGC.</title><p>Two groups of ORNs (ORNa and ORNb), each tuned to a specific individual pheromone component, converge onto ipsilateral specialist LNs (LNa and LNb) and onto a generalist LN (LNc). All LNs are interconnected via inhibitory synapses. The response of the model is observed at the intermediary LNi and the PN.</p></caption><graphic xlink:href="pone.0016308.g001"></graphic></fig><p>While not biased towards a particular coding paradigm by design, we demonstrate below that decision-making in this model is entirely reduced to the competition between LNs in terms of latency to the first spike: the LN that spikes first becomes the effective ‘winner’ even though it may receive less excitation on average during the course of stimulation than other LNs. The noisy nature of the ORN signalling, instantiated by independent Poisson processes in our model, and the different initial conditions of the LNs can thus lead to a significant number of errors in recognising the correct ratio. Accordingly, the main factor affecting the accuracy and performance of this model is the ORN-LN convergence rate which controls the signal to noise ratio in the ORN signal.</p><p>We then extend our model in two different ways. First, we consider a population of the discussed elementary recognition units and then a population unit, where the individual LNs are replaced with more realistic populations. We find that the population models allow for a more gradual “rate” competition and perform better than the elementary unit. We also tested the different models within the framework of approximately equivalent rate models (see <xref ref-type="supplementary-material" rid="pone.0016308.s001">Text S1</xref>).</p></sec></sec><sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title><sec id="s2a"><title>The basic model topology and units</title><p>Our basic model of the male MGC for a species utilizing a blend of two pheromone components can be schematically represented as shown in <xref ref-type="fig" rid="pone-0016308-g001">fig. 1</xref>.</p><sec id="s2a1"><title>Olfactory Receptor Neurons</title><p>We consider two groups of ORNs, each composed of cells that are narrowly tuned to respond very specifically to one of the two components of the female sex pheromone. The number of ORNs of each type is, as we argue below, an important factor affecting the accuracy of ratio detection. The spiking response of ORNs in the presence of pheromone has been shown to follow Poisson statistics with stimulus dependent rate over a wide range of stimulus intensities <xref ref-type="bibr" rid="pone.0016308-Schneider1">[26]</xref>. The firing rate <italic>λ</italic> of the Poissonian ORN spike trains of a given ORN population is approximately proportional to the logarithm of the concentration of the pheromone component it responds to. Specifically, Kaissling <xref ref-type="bibr" rid="pone.0016308-Kaissling1">[27]</xref> gives a dose–response curve for <italic>Antheraea polyphemus</italic> where a flow of 10<sup>−8</sup> µg/sec of (E,Z)-6,11-hexadecadienyl acetate causes ORNs to fire at approximately 10 Hz. With the pheromone concentrations up to 10<sup>−2</sup> µg/sec, the unadapted (<italic>i.e.</italic>, not previously exposed to the pheromone) ORNs remain sensitive and increase the firing rate to 300 Hz in a fairly linear fashion. The range of firing rates of the ORNs in our study is 10–200 Hz. In terms of pheromone concentrations, this approximately corresponds to 10<sup>−8</sup>–10<sup>−3</sup> µg/sec <xref ref-type="bibr" rid="pone.0016308-Kaissling1">[27]</xref>.</p><p>Our ORNs have a resting potential of −60 mV and emit 80 mV ‘durationless’ (<italic>i.e.</italic>, lasting for the current value of the integration time step) impulses.</p><p>As the summation of individual Poisson processes is still a Poisson process, further optimisation has been applied in which a group of ORNs is emulated by a single “compound” Poisson model neuron. The compound ORN has frequency <italic>λ</italic> = <italic>λ</italic><sub>s</sub> <italic>N</italic><sub>ORN</sub>, where <italic>λ</italic><sub>s</sub> is the frequency of a single ORN and <italic>N</italic><sub>ORN</sub> is the size of the represented ORN population. The high spiking rate generated by these compound ORNs results in a significant probability of multiple spikes occurring within a single time step Δ<italic>t</italic>. To account for this we modified the standard first order synapse model by Destexhe et al <xref ref-type="bibr" rid="pone.0016308-Destexhe2">[25]</xref>:<disp-formula><graphic xlink:href="pone.0016308.e001"></graphic><label>(1)</label></disp-formula><disp-formula><graphic xlink:href="pone.0016308.e002"></graphic><label>(2)</label></disp-formula>where <italic>α</italic> and <italic>β</italic> are the rates of neurotransmitter release and decay, respectively, and <italic>T</italic> is a step function which equals unity for <italic>t</italic><sub>release</sub> after each pre-synaptic spike to simulate the finite duration of transmitter release following a pre-synaptic spike. <italic>E</italic><sub>pre</sub> and <italic>E</italic><sub>post</sub> are the pre- and postsynaptic membrane potentials and <italic>E</italic><sub>rev</sub> the reversal potential of the synapse. The release of transmitter due to spikes within a single time step and across time steps within <italic>t</italic><sub>release</sub> of each other is summed according to<disp-formula><graphic xlink:href="pone.0016308.e003"></graphic><label>(3)</label></disp-formula>where <italic>T<sub>i</sub></italic> = <italic>q</italic> in the interval [<italic>t<sub>i</sub></italic>, <italic>t<sub>i</sub></italic>+<italic>t</italic><sub>release</sub>] if <italic>q</italic> spikes occurred at time <italic>t<sub>i</sub></italic> and the sum is over all relevant time steps <italic>t<sub>i</sub></italic> in [<italic>t</italic>−<italic>t</italic><sub>release</sub>, <italic>t</italic>]. <italic>α</italic> is determined such that for a single spike the resulting maximal synaptic activation <italic>S</italic> is the same in both the ORN population and compound ORN descriptions. This is the case for <inline-formula><inline-graphic xlink:href="pone.0016308.e004.jpg" mimetype="image"></inline-graphic></inline-formula> = 0.556 <italic>α</italic>.</p></sec><sec id="s2a2"><title>Convergence rate of ORNs to LNs</title><p>In male <italic>Manduca sexta</italic> antenna, there are 42,848±3,374 <italic>sensilla trichodea</italic>, each bearing on average two pheromone-sensitive ORNs <xref ref-type="bibr" rid="pone.0016308-Sanes1">[28]</xref>, with a total of about 85,700 ORNs. These project to the MGC to form synapses with LNs, of which Matsumoto & Hildebrand <xref ref-type="bibr" rid="pone.0016308-Matsumoto1">[29]</xref> distinguish several types. The authors identified a class of LNs only found in males and responding to pheromone substances as ‘Type II’ LNs and estimated the proportion of Type II LNs, in a sample of over 1000 cells stained, to be about 25% (namely, (86+68)/(260+86+68)). In a later work Homberg <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Homberg1">[30]</xref> confirmed LNs of <italic>M. sexta</italic> to have their somata grouped in three clusters on the periphery of the AL with a total cell count of 991 in the lateral cell group, 16, in the anterior, and 218 in the medial cell group (of these, only the latter is sexually dimorphic having 186 cells in females). The total number of LN somata in all cell groups is thus 1225, which includes LNs both pheromone-sensitive as well as those involved in the general olfactory system. Ultimately, from the observation that Matsumoto & Hildebrand <xref ref-type="bibr" rid="pone.0016308-Matsumoto1">[29]</xref> “have not detected any pattern in the locations of the cell bodies of the various types of neurons within the lateral and medial clusters,” the closest gross estimate for the number of pheromone-sensitive LNs is thus 25% of 1225, or some 300 cells.</p><p>Given the fact that ORN axons arborise in the glomeruli, and hence, form synapses with more than a single LN in the MGC, the actual convergence ratio can be larger than 85,700/300, or upwards of 285∶1.</p></sec><sec id="s2a3"><title>Local Neurons and the competition principle</title><p>We consider two ORN populations, ORNa amd ORNb, which are each selectively responsive to one of two pheromone components, ORNa to component A and ORNb to component B. LNs receive excitatory inputs from one of the two ORN populations (LNa and LNb, sensitive to pheromone components A and B, respectively) or from both (LNc). According to their innervation, LNa and LNb are termed ‘specialist’ LNs, and LNc is a ‘generalist’ LN. Both response types of LNs have been found experimentally <xref ref-type="bibr" rid="pone.0016308-Christensen3">[31]</xref>, <xref ref-type="bibr" rid="pone.0016308-Christensen1">[5]</xref>. All LNs are interconnected through inhibitory synapses, in agreement with the findings of numerous studies <xref ref-type="bibr" rid="pone.0016308-Christensen3">[31]</xref>, <xref ref-type="bibr" rid="pone.0016308-Schfer1">[32]</xref>, <xref ref-type="bibr" rid="pone.0016308-Waldrop1">[33]</xref>, <xref ref-type="bibr" rid="pone.0016308-Homberg2">[34]</xref>, <xref ref-type="bibr" rid="pone.0016308-Olsen1">[35]</xref> showing that connections between LNs in the glomeruli are inhibitory (GABAergic), though see <xref ref-type="bibr" rid="pone.0016308-Olsen2">[36]</xref>. When a pheromone blend is presented which has a component ratio sufficiently close to the target value, the ORNs will activate their LN targets and the competition between LNs, mediated by the local inhibitory connections will, given appropriate inter-LN synaptic strengths, cause the generalist LN to activate most and inhibit the specialist LNs. This in turn will inhibit the LNi, disinhibiting the PN and causing it to spike. Conversely, when the component ratio is deviating from the target value, one of the specialist LNs will suppress the other specialist LN as well as the generalist LN and prevent them from spiking. The LNi will continue spiking and there will be no response in the PN. In this way the system implements a winner-takes-all scenario where one of the specialist LNs ‘wins’ if its own component dominates the blend and the generalist LN wins if the blend has the correct target ratio. The generalist LNc is the exit point for the ‘ratio detected’ signal (to be picked up by a projection neuron and delivered further downstream to higher brain centres).</p></sec><sec id="s2a4"><title>Projection Neurons</title><p>The spiking response elicited by the stimulation at LNc will be dependent on the intensity of the stimulation, <italic>i.e.</italic>, the pheromone component concentration. Since our concern is ratio recognition <italic>irrespective of</italic> component concentration, we utilise the mechanism of disinhibition in the intermediate neuron LNi and the PN, in order to achieve a discrete ‘yes/no’ response as model output. This disinhibition pathway also allows to adhere to the biological principle that typically only one transmitter is used in any given neuron — the competing LNs are inhibitory whereas the PN is excitatory. The latter two are also spontaneously active (implemented <italic>via</italic> the injection of a constant current of 0.1 nA in the model). LN-LNi and LNi-PN synapses are inhibitory. One previous work in which this same mechanism has been studied in an AL model at the PN stage, is by Av-Ron and colleagues <xref ref-type="bibr" rid="pone.0016308-AvRon1">[37]</xref>, <xref ref-type="bibr" rid="pone.0016308-AvRon2">[38]</xref>, who also substantiate a biological relevance of disinhibition occurring in PNs.</p><p>All LNs and the PN in the model are described by Hodgkin–Huxley model neurons <xref ref-type="bibr" rid="pone.0016308-Hodgkin1">[22]</xref>, <xref ref-type="bibr" rid="pone.0016308-Traub1">[23]</xref> connected by generic synapses following Destexhe <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Destexhe1">[24]</xref> amended as explained above; for details and equations, see <xref ref-type="supplementary-material" rid="pone.0016308.s001">Text S1</xref>. While there have been recent advances on Ca<sup>2+</sup> currents in the LNs of the cockroach <xref ref-type="bibr" rid="pone.0016308-Husch1">[39]</xref>, <xref ref-type="bibr" rid="pone.0016308-Husch2">[40]</xref>, <xref ref-type="bibr" rid="pone.0016308-Husch3">[41]</xref> no complete and detailed voltage clamp data set or existing specific Hodgkin-Huxley models for moth AL LNs has been accessible to us. We, therefore, chose the popular, excitability type 1 Hodgkin-Huxley model neuron of Traub and Miles <xref ref-type="bibr" rid="pone.0016308-Hodgkin1">[22]</xref>, <xref ref-type="bibr" rid="pone.0016308-Traub1">[23]</xref> for the LNs and the PN. This neuron model allows a wide range of gradual changes in spike frequency which seemed appropriate for our study. The model has also already been used extensively in previous models of an insect AL <xref ref-type="bibr" rid="pone.0016308-Bazhenov1">[17]</xref>, <xref ref-type="bibr" rid="pone.0016308-Bazhenov2">[18]</xref> and seems to be well suited for allowing the competition dynamics that we seek to evaluate.</p></sec></sec><sec id="s2b"><title>Simulation of the Hodgkin–Huxley neurons</title><p>All ordinary differential equations for Hodgkin–Huxley neurons and corresponding rate models were integrated using a 6-5-order Runge–Kutta algorithm adapted from Press <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Press1">[42]</xref>, with adaptive time step Δ<italic>t</italic> confined to 10<sup>−6</sup>–0.5 ms and accuracy goal <italic>ε</italic> = 10<sup>−12</sup>.</p><p>The code of our simulations comes in two parts, a general simulation engine “cnrun”, available at <ext-link ext-link-type="uri" xlink:href="http://sourceforge.net/projects/cnrun">sourceforge.net/projects/cnrun</ext-link> or <ext-link ext-link-type="uri" xlink:href="http://johnhommer.com/academic/cnrun">johnhommer.com/academic/cnrun</ext-link>, and model descriptions in NeuroML, available at <ext-link ext-link-type="uri" xlink:href="http://johnhommer.com/academic/ratiocoding">johnhommer.com/academic/ratiocoding</ext-link>.</p><sec id="s2b1"><title>Rate-based reduction of the Hodgkin–Huxley model</title><p>Inherent in the competition approach is the issue that the timing of phasic events (spikes) may not reflect the true input ratio due to the stochasticity of the ORNs, and a resulting non-vanishing probability of locking the state of competing LNs into a ‘wrong state.’ To assess the relevance of spikes we studied the ratio recognition capability of an otherwise identical model built with a rate approximation of the Hodgkin–Huxley neurons (for the derivation of this model from the conductance-based equations, see <xref ref-type="supplementary-material" rid="pone.0016308.s001">Text S1</xref>). In the rate model, the LNc spiking rate averaged over the duration of each stimulus presentation replaces the spike density function in the response matrix (see below).</p></sec></sec><sec id="s2c"><title>Experimental protocol and trial cost function</title><p>We tested the model systematically with 100 different component concentration combinations, using firing rates <italic>λ</italic><sub>a</sub>, <italic>λ</italic><sub>b</sub> for the two populations of ORNs in increments of <italic>λ<sub>i</sub></italic> = 1.3 <italic>λ<sub>i</sub></italic><sub>−1</sub> for <italic>i</italic> = 0…9, starting with <italic>λ</italic><sub>0</sub> = 0.01 ms<sup>−1</sup> which allows us to cover a range of 10 to 100 Hz for ORN firing, a realistic range. Each stimulus frequency pair was presented for 250 ms, with 250 ms between presentations which was inspired by the ability of moths to resolve pulsed odor trains at this order of frequency.</p><p>Calculating the spike density function (SDF, defined as <italic>f</italic>(<italic>t</italic>) = ∑exp−(<italic>t</italic>−<italic>t<sub>i</sub></italic>)<sup>2</sup>/σ<sup>2</sup>, where <italic>t<sub>i</sub></italic> is the time of the <italic>i</italic>th spike and <italic>σ</italic> = 400 ms.) for the PN at the midpoint of each stimulus presentation, we obtained a 10×10 response matrix <bold><italic>r</italic></bold> for all combinations of input rates <italic>λ<sub>i</sub></italic>. A trial cost function <italic>s</italic>(<bold><italic>r</italic></bold>) is then defined as the “convolution” of the response matrix with a target 10×10 matrix <bold><italic><underline>r</underline></italic></bold>: <italic>s</italic>(<bold><italic>r</italic></bold>) = −Σ<italic><sub>i</sub></italic><sub>,<italic>j</italic></sub> <bold><italic>r</italic></bold><italic><sub>ij</sub></italic><bold><italic><underline>r</underline></italic></bold><italic><sub>ij</sub></italic> (<xref ref-type="fig" rid="pone-0016308-g002">fig. 2</xref>). A smaller cost function indicates a “better” response profile. Using this multiplicative approach, rather than, for example, a Euclidean distance of profiles, avoids penalising upward deviations from the target profile on the diagonal. The target profile <bold><italic><underline>r</underline></italic></bold> was chosen heuristically (see <xref ref-type="fig" rid="pone-0016308-g002">fig 2</xref> and its caption). We validated the choice of the target profile by evaluating the cost function for a number of observed response profiles and verifying that the ranking according to the cost function coincided with a manual assessment of the quality of the response profiles.</p><fig id="pone-0016308-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g002</object-id><label>Figure 2</label><caption><title>The target profile matrix for the conductance-based setup, with the x axis referring to λ<sub>a</sub><sup>i</sup> = λ<sub>0</sub>×1.3<sup>i</sup>, i = 0,…,9, the y axis to λ<sub>b</sub><sup>j</sup> = λ<sub>0</sub>×1.3<sup>j</sup>, j = 0,…,9 and λ<sub>0</sub> = 0.01 ms<sup>−1</sup>.</title><p>Individual weights are <inline-formula><inline-graphic xlink:href="pone.0016308.e005.jpg" mimetype="image"></inline-graphic></inline-formula> for the target ratios <italic>R</italic> of 1∶1 (panel <bold>A</bold>), 1∶3 (<bold>B</bold>) and 1∶9 (<bold>C</bold>), where <italic>C<sub>i</sub></italic> = 2×1.3<italic><sup>i</sup></italic> and <italic>a</italic> = 18, <italic>b</italic> = 1.25 and <italic>c</italic> = 0.3. The functional form and parameters of this target response profile were chosen heuristically to ensure that a) highest firing rates on the diagonal are favoured, b) the “punishment” of off-diagonal responses increases gradually with the distance from the diagonal for ratios close to the diagonal and c) responses far off the diagonal are equally and strongly “punished” in the cost function.</p></caption><graphic xlink:href="pone.0016308.g002"></graphic></fig><p>We optimised the performance of the model with respect to the cost function by adjusting the synaptic strength on the following connections: ORNs to specialist LNs; ORNs to generalist LN; and on all interconnections between LNs. The optimal set of parameters was determined by simulated annealing <xref ref-type="bibr" rid="pone.0016308-Press1">[42]</xref>.</p><p>We simulated the model behaviour with target component ratios of 1∶1 for ORN-LN convergence rates of 200, 500 and 1000, for both conductance and rate based versions. For target component ratios 1∶3 and 1∶9, we used only the conductance-based version with an ORN-LN convergence rate of 1000.</p><sec id="s2c1"><title>Decision latency</title><p>Besides recognition accuracy assessed by the cost function detailed above, another measure of model performance is the decision latency. We define the latter as the time between stimulus onset and the first spike that is part of a stable, unbroken spike train, <italic>i.e.</italic>, the time from stimulus onset until the system reaches a steady state, irrespective of whether the response occurs to the target stimulus ratio (correctly) or not. We calculate decision latencies on the level of LNs as well as of the PN.</p></sec></sec><sec id="s2d"><title>Cell population models</title><p>In the minimal model we use only one LN per competing group (loosely identified as a ‘glomerulus’). As already mentioned in the introduction, in this model, the competition between LNs is masked by the effect of the first spike. In our rate-based approximation, such phasic phenomena are absent, and as we will see, the recognition of pheromone component ratios is more reliable. However, as insightful as the rate-based approach may be from a theoretical standpoint, a full conductance-based, spiking model is more convincing in order to elicit and assess the real potential of competition based ratio recognition systems.</p><p>We, therefore, built two more realistic models based on our minimal unit: (a) ‘<italic>stacked model</italic>’ where ratio recognition units are independent from each other, only connecting at the level of LNi, which now receives convergent input from ten generalist LNs; and (b) ‘<italic>grouped model</italic>’, where all LNs of the same type are organised into groups (‘glomeruli’) with all-to-all inter-LN connections between glomeruli. In both cases, respective connection strengths were downscaled by a factor of 10, while the amount of ORN stimulation to the LNs was kept the same as in the single-unit experiments.</p></sec></sec><sec id="s3"><title>Results</title><sec id="s3a"><title>Minimalist model</title><p>We initially simulated the minimal model depicted in <xref ref-type="fig" rid="pone-0016308-g001">fig. 1</xref> with all combinations of 10 different ORN firing rates in the two ORN populations and observed the response of the PN. In these simulations we used synapses with neurotransmitter release rate constant <italic>α</italic> = (20 ms)<sup>−1</sup> and decay rate <italic>β</italic> = (50 ms)<sup>−1</sup>. After adjusting the ORN to LNa,b,c connections appropriately (<italic>g</italic><sub>ORN-LNsp</sub> roughly twice as strong as g<sub>ORN-LNgen</sub>), the target ratio of 1∶1 in the input rates elicits a reliable response in the PN neuron for a wide range of parameters for inter-LN synapses. However, we quite frequently also observed PN activation for non-target input ratios which could not be avoided by adjusting the parameters of the model by hand.</p><p>In an attempt to find an optimal set of parameters, we then used simulated annealing, with a cost function that rewarded PN activation for a 1∶1 target ratio and punished PN activation for all other input ratios (see <xref ref-type="sec" rid="s2">Methods</xref>). To minimise non-deterministic effects of the Poissonian input signals, the cost function was evaluated by averaging over three trials for each input concentration ratio. All synaptic connection strengths except the output disinhibition pathway were subject to change in this procedure (<italic>g</italic><sub>syn</sub> of ORN-LNsp, ORN-LNgen, LNsp-LNgen, LNgen-LNsp, and LNsp-LNsp). Note that because of the symmetric pheromone target ratio of 1∶1, we fix <italic>g</italic><sub>ORNa-LNa</sub> = <italic>g</italic><sub>ORNb-LNb</sub> = <italic>g</italic><sub>ORN-LNsp</sub> and <italic>g</italic><sub>ORNa-LNc</sub> = <italic>g</italic><sub>ORNb-LNc</sub> = <italic>g</italic><sub>ORN-LNgen</sub>.</p><p>In the following, inter-LN connection strengths are abbreviated as, respectively, <italic>g</italic><sub>LNi</sub> = <italic>g</italic><sub>LNsp-LNgen</sub>, <italic>g</italic><sub>LNo</sub> = <italic>g</italic><sub>LNgen-LNsp</sub> and <italic>g</italic><sub>LNt</sub> = <italic>g</italic><sub>LNsp-LNsp</sub> (for <italic>inbound</italic>, <italic>outbound</italic> and <italic>transverse</italic>) and we assume symmetry of connection strengths between <italic>a</italic> and <italic>b</italic> populations again based on the 1∶1 symmetric target ratio. We tested three models with ORN-LN convergence ratios 200∶1, 500∶1 and 1000∶1 and both the full Hodgkin–Huxley model and an equivalent rate model (see <xref ref-type="sec" rid="s2">Methods</xref> and <xref ref-type="supplementary-material" rid="pone.0016308.s001">Text S1</xref>) were used. For each model up to 40,000 iterations of the simplex annealing have been completed. The algorithm was terminated when the initial simplex of parameter vectors has shrunk to 0.01% of its initial size. The resulting best parameters found are shown in <xref ref-type="table" rid="pone-0016308-t001">Table 1</xref>.</p><table-wrap id="pone-0016308-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.t001</object-id><label>Table 1</label><caption><title>Resulting parameters of the simplex annealing of the 1∶1 ratio recognition, for the various ORN-LN convergence rates.</title></caption><alternatives><graphic id="pone-0016308-t001-1" xlink:href="pone.0016308.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td><td colspan="3" align="left" rowspan="1">Conductance-based</td><td colspan="3" align="left" rowspan="1">Rate-based</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">200</td><td align="left" rowspan="1" colspan="1">500</td><td align="left" rowspan="1" colspan="1">1000</td><td align="left" rowspan="1" colspan="1">200</td><td align="left" rowspan="1" colspan="1">500</td><td align="left" rowspan="1" colspan="1">1000</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>g</italic><sub>ORN-LNsp</sub> [µS]</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">0.094</td><td align="left" rowspan="1" colspan="1">0.103</td><td align="left" rowspan="1" colspan="1">9.47·10<sup>−3</sup></td><td align="left" rowspan="1" colspan="1">0.086</td><td align="left" rowspan="1" colspan="1">0.034</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>g</italic><sub>ORN-LNgen</sub> [µS]</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">0.072</td><td align="left" rowspan="1" colspan="1">0.074</td><td align="left" rowspan="1" colspan="1">6.06·10<sup>−3</sup></td><td align="left" rowspan="1" colspan="1">0.045</td><td align="left" rowspan="1" colspan="1">0.021</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>g</italic><sub>LNi</sub> [µS]</td><td align="left" rowspan="1" colspan="1">0.185</td><td align="left" rowspan="1" colspan="1">0.490</td><td align="left" rowspan="1" colspan="1">0.499</td><td align="left" rowspan="1" colspan="1">0.146</td><td align="left" rowspan="1" colspan="1">0.153</td><td align="left" rowspan="1" colspan="1">0.098</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>g</italic><sub>LNt</sub> [µS]</td><td align="left" rowspan="1" colspan="1">0.185</td><td align="left" rowspan="1" colspan="1">0.500</td><td align="left" rowspan="1" colspan="1">0.389</td><td align="left" rowspan="1" colspan="1">0.081</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.187</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>g</italic><sub>LNo</sub> [µS]</td><td align="left" rowspan="1" colspan="1">0.184</td><td align="left" rowspan="1" colspan="1">0.347</td><td align="left" rowspan="1" colspan="1">0.386</td><td align="left" rowspan="1" colspan="1">0.076</td><td align="left" rowspan="1" colspan="1">0.181</td><td align="left" rowspan="1" colspan="1">0.099</td></tr><tr><td align="left" rowspan="1" colspan="1">Cost function</td><td align="left" rowspan="1" colspan="1">−1092</td><td align="left" rowspan="1" colspan="1">−1081</td><td align="left" rowspan="1" colspan="1">−1110</td><td align="left" rowspan="1" colspan="1">−2.65</td><td align="left" rowspan="1" colspan="1">−3.20</td><td align="left" rowspan="1" colspan="1">−3.69</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>The synaptic strengths are in µS.</p></fn></table-wrap-foot></table-wrap><p>The model performance with respect to parameter changes as revealed in the course of annealing, is illustrated in <xref ref-type="fig" rid="pone-0016308-g003">fig 3</xref>. Low cost function parameter sets (good solutions) are coloured green to blue whereas worse solutions are in red colours (see colour bar and not the logarithmic scale). There are some very clear features in both the projection to the parameter space of ORN-LN connection strengths (panel A) and of the LN-LN connection strengths (panel B).</p><fig id="pone-0016308-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g003</object-id><label>Figure 3</label><caption><title>The cloud of data points (with cost function <0, <italic>n</italic> = 159,180) obtained in the simplex annealing of a basic conductance-based model with an ORN∶LN convergence rate of 1000.</title><p>The colour of data points corresponds to the cost function value observed (see color bar). In each subfigure, the point with the best cost function is marked with a cross.</p></caption><graphic xlink:href="pone.0016308.g003"></graphic></fig><sec id="s3a1"><title>ORN-LN connection strengths</title><p>Obviously, there is a trade-off between having a reliable PN response for 1∶1 ratios (large enough <italic>g</italic><sub>ORN-LNgen</sub>) and spurious PN activation for close by ratios (too large <italic>g</italic><sub>ORN-LNgen</sub>). Naively, we expected <italic>g</italic><sub>ORN-LNsp</sub> = <italic>2 g</italic><sub>ORN-LNgen</sub>. The best performing parameter set in the simulated annealing appears to have a somewhat larger <italic>g</italic><sub>ORN-LNgen</sub>, with the ratio being approximately 1.7–1.8 ORNsp∶ORNgen (<xref ref-type="table" rid="pone-0016308-t001">Table 1</xref>, <xref ref-type="fig" rid="pone-0016308-g003">Fig. 3A</xref>). This may surprise on first sight but it is perfectly reasonable as at equality of ORNa and ORNb firing, the generalist neuron has to have <italic>stronger</italic> input than the specialists, not only equal input. This bias must be enough for reliable LNgen responses with equal input but weak enough to have LNsp take over for the first off-diagonal terms (1∶1.3 and 1.3∶1 ratios). The region of successful parameter combinations of <italic>g</italic><sub>ORN-LNsp</sub> and <italic>g</italic><sub>ORN-LNgen</sub> is very well constrained indicating that the correct ratio between these conductances is very important for the correct recognition of pheromone component concentration ratio.</p></sec><sec id="s3a2"><title>Inter-LN connection strengths</title><p>In contrast to the ORN-LN connections, our results do not indicate any clear rule for the inter-LN connection strengths other than that they have to be of approximately the same order of magnitude and not too weak overall. In this case, the “cloud” of appropriate synaptic strengths is more wide-spread (<xref ref-type="fig" rid="pone-0016308-g003">fig. 3B</xref>), although still well-delineated, with successful systems having an approximate ratio of 1∶1∶1 (LNi∶LNt∶LNo). Notably, there is a “side-arm” of parameter combinations for small <italic>g</italic><sub>LNo</sub> and <italic>g</italic><sub>LNi</sub> with strongly varying values for the strength of <italic>g</italic><sub>LNt</sub>, indicating that the value of <italic>g</italic><sub>LNt</sub> is not essential for success in this area of parameter space. At close view, the cost function appears to be considerably scattered in repeated trials with the same or very close parameters due to the random nature of the Poissonian input, i.e., points with high cost function (red) can be very close to points with very low cost (green/blue) depending on the instantiations of the Poisson processes describing the ORNs. Even though we minimized the variability of the estimated cost function by taking averages over three trials for each parameter set, some of this uncertainty remains.</p><p>We note that this result indicates that while the existence of inter-LN competition is essential for ratio recognition, the exact details of it do not seem to matter. This is an interesting results at it seems to imply that no costly detailed genetic specification of inter-LN connectivity is necessary to achieve a functional MGC.</p></sec><sec id="s3a3"><title>The balance of the inter-LN synaptic strengths</title><p>To elucidate the effect of inter-LN connections further, we obtained simulation results from a combination of discrete inter-LN <italic>g</italic> values set apart by a factor of 1.5 from the optimal solution as shown in <xref ref-type="fig" rid="pone-0016308-g004">fig. 4</xref>. The most detrimental change is caused by breaking the equality of <italic>g</italic><sub>LNi</sub> and its reciprocal <italic>g</italic><sub>LNo</sub>, resulting in the loss of ratio selectivity towards higher pheromone component concentrations if <italic>g</italic><sub>LNi</sub><<italic>g</italic><sub>LNo</sub> and partial loss of response if <italic>g</italic><sub>LNi</sub>><italic>g</italic><sub>LNo</sub> <italic>and g</italic><sub>LNi</sub>><italic>g</italic><sub>LNt</sub>. Most interestingly, given <italic>g</italic><sub>LNi</sub> = <italic>g</italic><sub>LNo</sub>, the value of <italic>g</italic><sub>LNt</sub> does not seem to affect the ratio recognition across the entire concentration range, as can be seen in the fact that the case where g<sub>LNi</sub> and g<sub>LNo</sub> are elevated (“<bold>o_o</bold>”) and the case where g<sub>LNt</sub> is elevated (“<bold>_o_</bold>”) are virtually identical despite having different <italic>g</italic><sub>LNt</sub>. This agrees well with the observation of the extension of the cloud of good solutions in <xref ref-type="fig" rid="pone-0016308-g003">fig. 3</xref> towards many different <italic>g</italic><sub>LNt</sub> values. The rate-based model (<xref ref-type="fig" rid="pone-0016308-g004">fig. 4B</xref>) behaved similarly with respect to these observations.</p><fig id="pone-0016308-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g004</object-id><label>Figure 4</label><caption><title>The balance of inter-LN connection strengths affects the ratio recognition accuracy.</title><p>A) Response profiles in the conductance-based setup. B) Response profiles in the rate-based setup. The three symbols in each panel indicate, in order, the synaptic strengths on LNsp-LNgen, LNsp-LNsp and LNgen-LNsp connections, where “<bold>_</bold>” is the value for that connection found in the annealing optimisation (<xref ref-type="table" rid="pone-0016308-t001">Table 1</xref>), and “<bold>o</bold>” is that value times 1.3 (A) or 1.5 (B). The <italic>x</italic>- and <italic>y</italic>-axes represent the concentration of individual components in the blend as the log of the Poisson firing rate of ORNs; the colour code of the squares shows the SDF (spike-density function) (A) and the average spiking rate (B), see respective colour bars. ORN-LN convergence rate used was 1000.</p></caption><graphic xlink:href="pone.0016308.g004"></graphic></fig><p>In a repeat of the above tests with baseline inter-LN conductances all set to 0.18 µS (rather than to the optimal values found in the annealing procedure (<xref ref-type="table" rid="pone-0016308-t002">Table 2</xref>)), this baseline case still had the best cost function of all combinations studied (data not shown).</p><table-wrap id="pone-0016308-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.t002</object-id><label>Table 2</label><caption><title>Decision latencies measured at the PN, correlated (Pearson's product-moment) with the cost function, in a grouped arrangement of 10 LNs per glomerulus.</title></caption><alternatives><graphic id="pone-0016308-t002-2" xlink:href="pone.0016308.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td colspan="2" align="left" rowspan="1">Case</td><td align="left" rowspan="1" colspan="1">Cost Function</td><td colspan="2" align="left" rowspan="1">PN</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>α</italic> [ms<sup>−1</sup>]</td><td align="left" rowspan="1" colspan="1"><italic>β</italic> [ms<sup>−1</sup>]</td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Latency [ms]</td><td align="left" rowspan="1" colspan="1">corr. coeff.</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">−1145±51.8</td><td align="left" rowspan="1" colspan="1">137.4±4.19</td><td align="left" rowspan="1" colspan="1">−0.18*</td></tr><tr><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">−1096±47.1</td><td align="left" rowspan="1" colspan="1">131.2±4.20</td><td align="left" rowspan="1" colspan="1">−0.27†</td></tr><tr><td align="left" rowspan="1" colspan="1">0.060</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">−1182±57.3</td><td align="left" rowspan="1" colspan="1">139.1±3.97</td><td align="left" rowspan="1" colspan="1">−0.37†</td></tr><tr><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.021</td><td align="left" rowspan="1" colspan="1">−1180±42.5</td><td align="left" rowspan="1" colspan="1">141.5±4.55</td><td align="left" rowspan="1" colspan="1">−0.32†</td></tr><tr><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">−1084±48.3</td><td align="left" rowspan="1" colspan="1">131.4±3.94</td><td align="left" rowspan="1" colspan="1">−0.29†</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Measurements were taken for different combinations of neurotransmitter release rise rate <italic>α</italic> and decay rate <italic>β</italic> (shown in ms<sup>−1</sup>) on inter-LN synapses. Higher cost function means better response. <italic>n</italic> = 120 for the case of baseline <italic>α</italic>−<italic>β</italic> combination, <italic>n</italic> = 100 for all other cases; values ±SD; (*) marks significance at <italic>p</italic> = 0.05 level, (†), at <italic>p</italic> = 0.01 level, n.s. means the <italic>p</italic> value is greater than 0.05.</p></fn></table-wrap-foot></table-wrap><p>We observed a surprising number of false positive responses, in particular towards higher total concentrations where the system responded even though the correct input ratio was not present. As we are convinced that the solutions found are optimal within the given architecture and constraints on parameters, we must assume that there are reasons why in principle a more specific response cannot be achieved. In order to uncover these underlying reasons, we undertook a closer analysis of the cases in which the system responded incorrectly to its input.</p></sec></sec><sec id="s3b"><title>Estimation of the probability of errors in ratio detection</title><sec id="s3b1"><title>“Spurious spikes”</title><p>Close analysis of the raw data revealed that the spurious activations for non-target ratio inputs are due to the excitation from both ORN populations producing the first spike in the generalist LN (LNc) before the excitation from the more active population (<italic>e.g.</italic>, ORNb) can activate its specialist LN (LNb). This occurs even though the rate of the ORNb population is higher than the rate of the ORNa population and thus the input to LNb should be stronger than the combined ORNa+ORNb input to LNc. The “spurious spike” in LNc will then cause a surge of inhibition onto the specialist LNs. Due to the time scales of the synapses, and the minimal value for the inter-LN conductances that allows winner-take-all dynamics, the outcome of the competition between LNs turns out to be determined largely by the first spike, <italic>i.e.</italic>, it is a competition in the latency to the first spike. Only in rare cases can the specialist LN overcome a spurious initial activation of LNc. As a result, LNs are locked into the wrong state for a significant part of or the entire time of stimulus presentation, and a wrong signal is sent downstream.</p><p>The incidence of spurious spikes depends on the properties of the Poisson spike trains arriving from the two ORN populations, in particular the firing rates <italic>λ</italic><sub>ORNa</sub> and <italic>λ</italic><sub>ORNb</sub> of the ORNs of each type and the number of ORNs converging onto each LN, as well as on the input integration in the LNs. The higher the overall population input rate to each LN, the lower is the variability of the summed input and hence, the probability of spurious spikes.</p><p>To focus on the LN latency based on input and cellular properties, we removed all inter-LN connections and used a simplified protocol including only the 1∶1 stimulus ratios and the ratio immediately adjacent to 1∶1, i.e. 1∶1.3, and estimated numerically the incidence of LN activations for several ORN-LN convergence rates, in the low, medium and high portions of the stimulus intensity, as shown in <xref ref-type="fig" rid="pone-0016308-g005">fig. 5</xref>. The lines descending as a function of the ratio <italic>g</italic><sub>ORN-LNsp</sub>/<italic>g</italic><sub>ORN-LNgen</sub> are the probability of LNc responses to the target 1∶1 ratio (true positives for 1∶1) and the ascending lines are the responses of LNa/LNb to the non-target ratios (true positives for 1∶1.3). There are three separate lines for three different exemplary overall concentrations, low (red), medium (green) and high (blue). Optimally, we would expect to find <italic>g</italic><sub>ORN-LNsp</sub>/<italic>g</italic><sub>ORN-LNgen</sub> such that true positives in both cases are (close to) 100%. However, there is no region in the plots where both percentages are high, not even for a single concentration value let alone for all three concentrations simultaneously. This convinced us that the model with single LNs as functional units, due to it being reduced to a latency-to-first-spike system, can on principle not perform to a better accuracy. We, therefore, extended our models to population models with the rationale that populations of neurons would be able to respond more gradually such that the latency to first spike does not dominate the competition. Alternatively, a population of the single-cell units could lead to a more informed response due to a second summation and thresholding process at the LNi, similar to boosting in machine learning.</p><fig id="pone-0016308-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g005</object-id><label>Figure 5</label><caption><title>Percentage of LNc responses as a function of the ratio <italic>g</italic><sub>ORN-LNsp</sub>/<italic>g</italic><sub>ORN-LNgen</sub>.</title><p><bold>A–D</bold>, ORN∶LN convergence rates of 100, 200, 500 and 1000. The descending lines in each panel represent the percentage of responses to the 1∶1 ratio where LNc spikes first (true positives for 1∶1); ascending lines, the percentage of responses where an LNsp (specifically, LNb) fired first (true positives for non-1∶1 ratios immediately adjacent to the target 1∶1 ratio). The colours correspond to low (red), medium (green) and high (blue) pheromone blend concentrations. Note that there is no value of <italic>g</italic><sub>ORN-LNsp</sub>/<italic>g</italic><sub>ORN-LNgen</sub> where both true positives are high simultaneously.</p></caption><graphic xlink:href="pone.0016308.g005"></graphic></fig></sec></sec><sec id="s3c"><title>Extended models</title><p>We have taken a total of 100 repeated measurements of the performance in the two types of extended models using two arrangements of the elementary ratio recognition unit: stacked (multiple copies of the basic unit circuit converge onto LNi) and grouped (each LN is replaced by a population of like copies). The best outcomes of this test are shown in <xref ref-type="fig" rid="pone-0016308-g006">fig. 6</xref>.</p><fig id="pone-0016308-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g006</object-id><label>Figure 6</label><caption><title>Response profiles for the population models.</title><p>The grey scale plots show the spike density function of the PN in response to the different input frequencies in the stacked arrangement (A) and the grouped arrangement (B). These are the best outcomes of a total of 120 trials.<xref ref-type="fig" rid="pone-0016308-g007">Figure 7</xref>. Per-stimulus pair decision latencies (in ms) observed at the individual LNs, averaged over 120 trials (<bold>A</bold>, stacked arrangement, <bold>B</bold>, grouped arrangement).</p></caption><graphic xlink:href="pone.0016308.g006"></graphic></fig><p>In order to identify the character of the competition process we introduced and analysed the measure of decision latency, defined as the time until the system reaches a form of steady state in response to a rectangular onset of a constant odour stimulus (see methods for details). Specifically, where a neuron did not produce a spike train, or produced one that did not last to the end of the stimulus presentation, no data point was taken. Similarly, when spike trains persisted at LNs from more than a single LN group during the entire stimulus presentation (thus producing multiple effective winners), such data points were also omitted. Many (four to six) incidents of such multiple winners systematically appeared in every trial in the stacked case, which is not surprising as the constituent elementary units are completely isolated. In the grouped case, multiple winners were less frequent.</p><p>The results of the latency measurements are illustrated in <xref ref-type="fig" rid="pone-0016308-g007">figure 7</xref>. As expected, the difficult off-diagonal ratios lead to the longest latencies followed by the 1∶1 ratio at small concentrations. It is interesting that the latency depends quite conspicuously on concentration, unlike the final outcome that, predominantly, does not.</p><fig id="pone-0016308-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g007</object-id><label>Figure 7</label><caption><title>Decision latencies of LNs.</title><p>The grey scale plots show the decision latency (in ms) per stimulus pair observed at the individual LNs, averaged over 120 trials. A) stacked arrangement. B) grouped arrangement.</p></caption><graphic xlink:href="pone.0016308.g007"></graphic></fig><p>To elucidate the relationship between the decision latency and the correct ratio recognition, we calculated the correlation between latency and correct recognition. We observed that in the grouped case, there is a significant <italic>negative</italic> correlation between latency and correctness, whereas we did not observe significant correlations in the stacked case (<xref ref-type="fig" rid="pone-0016308-g008">fig. 8</xref>). This came as a surprise because the intuition had been that a longer deliberation time (longer decision latency) should yield a better response. However, in this case the longer latency is caused by the problem being harder (the ratio being closer to 1∶1) which has the opposite effect on the response accuracy. Our simulation results indicate that in the trade-off between the two factors the second one appears to prevail.</p><fig id="pone-0016308-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0016308.g008</object-id><label>Figure 8</label><caption><title>Correlation of average per-trial PN latency and the cost function in the extended models.</title><p>(<bold>A</bold>, stacked arrangement, <bold>B</bold>, grouped arrangement).</p></caption><graphic xlink:href="pone.0016308.g008"></graphic></fig><p>Listed in <xref ref-type="table" rid="pone-0016308-t002">table 2</xref> are additional latency-accuracy correlations we have computed supporting this point of view. In the combination of higher <italic>α</italic> and lower <italic>β</italic> (which clearly means more stimulation overall), we see both a higher average cost function and stronger latency-cost function link. The former simply suggests our initial values for LN <italic>α</italic> and <italic>β</italic> were a little off the optimal mark.</p><p>The fact that decision latency for harder problems appears to be longer seems on first sight to contradict the results of <xref ref-type="bibr" rid="pone.0016308-Ditzen1">[46]</xref> who found in bees that the behavioural decision latency for recognizing binary mixtures did not depend on the difficulty of the problem. However, the latencies observed in this study (∼700 ms) are almost an order of magnitude larger than the latencies in our model (∼100 ms) and observed during upwind flight in moths (∼150 ms, see <xref ref-type="bibr" rid="pone.0016308-Todd1">[47]</xref>). We conclude that pheromone blend recognition in moths and recognition of binary mixtures in the general olfactory system of bees may not be directly comparable.</p></sec></sec><sec id="s4"><title>Discussion</title><p>The simple model we propose represents an elementary ratio recognition unit, consisting of three LNs innervated by two groups of ORNs and connecting to a blend-specific PN through a disinhibition pathway. Of these units, in the moth Antennal Lobe there could be many, with variable connection strengths and hence, variable sensitivity and target ratio.</p><p>The pheromone communication in <italic>Heliothis virescens</italic> and an allied species <italic>Helicoverpa zea</italic> is an intricate case where a principal pheromone component, (Z)-11-hexadecenal (Z11–16:AL) is shared by both species, while two other compounds, (Z)-9-tetradecenal (Z9–14:AL) and (Z)-9-hexadecenal (Z9–16:AL) are, respectively, secondary for <italic>H. virescens</italic> and <italic>H. zea</italic>. In field studies, Vickers <italic>et al</italic> <xref ref-type="bibr" rid="pone.0016308-Vickers1">[43]</xref> reported fairly loose accuracy of component ratio recognition in <italic>H. virescens</italic> and <italic>H. zea</italic>. The authors distinguished four gradual behaviours, the most basic of which (taking flight) is enacted by a much broader pheromone component ratio than the ultimate behaviour (contacting the female). This lends support to our conception of multiple elementary units variously adapted, or ‘tuned,’ to detect the target pheromone component ratio. In this scenario, progressively greater populations of such units would be recruited as the strength of the signal increases and/or the component ratio approaches the target one.</p><p>Pheromone signalling has been shown <xref ref-type="bibr" rid="pone.0016308-Vickers2">[44]</xref> to sometimes involve an ‘antagonist’ component which disrupts the behaviour normally evoked by the ‘agonist’ pheromone components in the right proportion. Hansson et al <xref ref-type="bibr" rid="pone.0016308-Hansson2">[45]</xref>, who studied the axonal arborisations of ORNs sensitive to the two agonist pheromone components in <italic>H. virescens</italic>, (Z)-9-tetradecenal (Z9–14:AL) and (Z)-11-hexadecenal (Z11–16:AL), found that they “do not display a clear-cut morphological separation into different glomeruli in the MGC” (denoted <italic>a</italic> and <italic>b</italic>) but rather innervate both, whereas the antagonist component's ORNs all terminate in a single glomerulus <italic>c</italic>. In this arrangement, the agonist and antagonist paths are independent from each other; the decision to trigger female-seeking behaviour is deferred to higher brain centres, where an agonist signal enables it unless overridden by the antagonist signal. In this case, our model would only cover the function of the agonist part.</p><p>Of the constraints on the synaptic strengths, we found that the ORN-LN connection strengths are much better constrained than the inter-LN connections, as evident from our cost function optimisation. Within the inter-LN connections it only seems to matter that g<sub>LNi</sub> must be close to g<sub>LNo</sub>. If this requirement is met, the inter-specialist LN connection strength (LNt) may vary to a large extent, both upward of the LNi and LNo absolute value and downward. One probable interpretation of this result is that the competition between specialist LNs is not instrumental to the success of ratio recognition because the readout of the result is entirely based on the generalist LN population. If the two specialist populations fail to compete properly against each other but the competition of either with the generalist LNs is intact, the recognition of pheromone component ratios is likely almost unimpaired.</p><p>Our original intuition was that the MGC would achieve ratio recognition by straightforward competition in terms of the firing rate of populations of LNs which could have equivalently been expressed in a rate model. By considering the more detailed Hodgkin-Huxley model we found that the competition was realised in terms of latency to first spike rather than rate. One consequence of this latency coding in the system was a broadening of responses for higher overall pheromone concentrations which seems unintuitive from a pure rate perspective, because higher overall rates should have improved the signal to noise ratio. However, this broadening is observed experimentally (S. Anton, personal communication) supporting our results. Consequently, we feel that to advance the understanding of biological systems it is necessary to appreciate the way the account of a biological phenomenon may change depending on the type of model description.</p></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0016308.s001"><label>Text S1</label><caption><p><bold>The inference of the rate-based model of the Hodgkin–Huxley neuron.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0016308.s001.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>We are grateful to Sylvia Anton for helpful discussions.</p></ack><fn-group><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was funded by the Biotechnology and Biological Sciences Research Council (grant number BB/F005113/1) within the French-British ANR BBSRC SysBio 006 01 “Pherosys” initiative. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group><ref-list><title>References</title><ref id="pone.0016308-Ljungberg1"><label>1</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Ljungberg</surname><given-names>H</given-names></name><name><surname>Anderson</surname><given-names>P</given-names></name><name><surname>Hansson</surname><given-names>BS</given-names></name></person-group> <year>1993</year> <article-title>Physiology and morphology of pheromone-specific sensilla on the antennae of male and female Spodoptera littoralis (Lepidoptera: Noctuidae).</article-title> <source>J Insect Physiol</source> <volume>39</volume> <fpage>253</fpage> <lpage>260</lpage> </mixed-citation></ref><ref id="pone.0016308-Anton1"><label>2</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Anton</surname><given-names>S</given-names></name><name><surname>Hansson</surname><given-names>BS</given-names></name></person-group> <year>1994</year> <article-title>Central processing of sex pheromone, host odour, and oviposition deterrent information by interneurons in the antennal lobe of female Spodoptera littoralis (Lepidoptera: Noctuidae).</article-title> <source>J Comp Neurol</source> <volume>350</volume> <fpage>199</fpage> <lpage>214</lpage> <pub-id pub-id-type="pmid">7884038</pub-id></mixed-citation></ref><ref id="pone.0016308-Tumlinson1"><label>3</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Tumlinson</surname><given-names>JH</given-names></name><name><surname>Brennan</surname><given-names>MM</given-names></name><name><surname>Doolittle</surname><given-names>RE</given-names></name><name><surname>Mitchell</surname><given-names>ER</given-names></name><name><surname>Brabham</surname><given-names>A</given-names></name><etal></etal></person-group> <year>1989</year> <article-title>Identification of a pheromone blend attractive to Manduca sexta (L.) males in a wind tunnel.</article-title> <source>Arch Insect Biochem Physiol</source> <volume>10</volume> <fpage>255</fpage> <lpage>271</lpage> </mixed-citation></ref><ref id="pone.0016308-Tth1"><label>4</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Tòth</surname><given-names>M</given-names></name><name><surname>Löfstedt</surname><given-names>C</given-names></name><name><surname>Blair</surname><given-names>BW</given-names></name><name><surname>Cabello</surname><given-names>T</given-names></name><name><surname>Farag</surname><given-names>AI</given-names></name><etal></etal></person-group> <year>1992</year> <article-title>Attraction of male turnip moths Agrotis segetum (Lepidoptera: Noctuidae) to sex pheromone components and their mixtures at 11 sites in Europe, Asia, and Africa.</article-title> <source>J Chem Ecol</source> <volume>18</volume> <fpage>1337</fpage> <lpage>1347</lpage> </mixed-citation></ref><ref id="pone.0016308-Christensen1"><label>5</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Mustaparta</surname><given-names>H</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1989</year> <article-title>Discrimination of sex-pheromone blends in the olfactory system of the moth.</article-title> <source>Chem Senses</source> <volume>14</volume> <fpage>463</fpage> <lpage>477</lpage> </mixed-citation></ref><ref id="pone.0016308-Christensen2"><label>6</label><mixed-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1994</year> <article-title>Neuroethology of sexual attraction and inhibition in heliothine moths.</article-title> <person-group person-group-type="editor"><name><surname>Schildberger</surname><given-names>K</given-names></name><name><surname>Elsner</surname><given-names>N</given-names></name></person-group> <source>Neural basis of behavioural adaptations—Progress in zoology</source> <publisher-loc>Stuttgart</publisher-loc> <publisher-name>Gustav Fischer Verlag</publisher-name> <fpage>37</fpage> <lpage>46</lpage> </mixed-citation></ref><ref id="pone.0016308-Ishikawa1"><label>7</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Ishikawa</surname><given-names>Y</given-names></name><name><surname>Takanashi</surname><given-names>T</given-names></name><name><surname>Kim</surname><given-names>Choong-gon</given-names></name><name><surname>Hoshizaki</surname><given-names>S</given-names></name><name><surname>Tatsuki</surname><given-names>S</given-names></name><etal></etal></person-group> <year>1999</year> <article-title>Ostrinia spp. in Japan: their host plants and sex pheromones 1999.</article-title> <source>Entomol Exp Appl</source> <volume>91</volume> <fpage>237</fpage> <lpage>244</lpage> </mixed-citation></ref><ref id="pone.0016308-Linster1"><label>8</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Linster</surname><given-names>C</given-names></name><name><surname>Masson</surname><given-names>C</given-names></name><name><surname>Kerszberg</surname><given-names>M</given-names></name><name><surname>Personnaz</surname><given-names>L</given-names></name><name><surname>Dreyfus</surname><given-names>G</given-names></name></person-group> <year>1993</year> <article-title>Computational diversity in a formal model of the insect olfactory macroglomerulus.</article-title> <source>Neural Comput</source> <volume>5</volume> <fpage>228</fpage> <lpage>241</lpage> </mixed-citation></ref><ref id="pone.0016308-Hansson1"><label>9</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Hansson</surname><given-names>BS</given-names></name><name><surname>Anton</surname><given-names>S</given-names></name></person-group> <year>2000</year> <article-title>Function and morphology of the antennal lobe: new developments.</article-title> <source>Annu Rev Entomol</source> <volume>45</volume> <fpage>203</fpage> <lpage>231</lpage> <pub-id pub-id-type="pmid">10761576</pub-id></mixed-citation></ref><ref id="pone.0016308-Wu1"><label>10</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wu</surname><given-names>W</given-names></name><name><surname>Anton</surname><given-names>S</given-names></name><name><surname>Löfstedt</surname><given-names>C</given-names></name><name><surname>Hansson</surname><given-names>BS</given-names></name></person-group> <year>1996</year> <article-title>Discrimination among pheromone component blends by interneurons in male antennal lobes of two populations of the turnip moth, Agrotis segetum.</article-title> <source>Proc Natl Acad Sci U S A</source> <volume>93</volume> <fpage>8022</fpage> <lpage>8027</lpage> <pub-id pub-id-type="pmid">8755596</pub-id></mixed-citation></ref><ref id="pone.0016308-Linster2"><label>11</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Linster</surname><given-names>C</given-names></name><name><surname>Kerszberg</surname><given-names>M</given-names></name><name><surname>Masson</surname><given-names>C</given-names></name></person-group> <year>1994</year> <article-title>How neurons may compute: the case of insect sexual pheromone discrimination.</article-title> <source>J Comput Neurosci</source> <volume>1</volume> <fpage>231</fpage> <lpage>238</lpage> <pub-id pub-id-type="pmid">8792232</pub-id></mixed-citation></ref><ref id="pone.0016308-Linster3"><label>12</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Linster</surname><given-names>C</given-names></name><name><surname>Dreyfus</surname><given-names>G</given-names></name></person-group> <year>1996</year> <article-title>A model for pheromone discrimination in the insect antennal lobe: Investigation of the role of neuronal response pattern complexity.</article-title> <source>Chem Senses</source> <volume>21</volume> <fpage>19</fpage> <lpage>27</lpage> <pub-id pub-id-type="pmid">8646488</pub-id></mixed-citation></ref><ref id="pone.0016308-Heinbockel1"><label>13</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Heinbockel</surname><given-names>T</given-names></name><name><surname>Kloppenburg</surname><given-names>P</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1998</year> <article-title>Pheromone-evoked potentials and oscillations in the antennal lobes of the sphinx moth Manduca sexta.</article-title> <source>J Comp Physiol [A]</source> <volume>182</volume> <fpage>703</fpage> <lpage>714</lpage> </mixed-citation></ref><ref id="pone.0016308-Laurent1"><label>14</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Laurent</surname><given-names>G</given-names></name></person-group> <year>1996</year> <article-title>Dynamical representation of odors by oscillating and evolving neural assemblies.</article-title> <source>Trends Neurosci</source> <volume>19</volume> <fpage>489</fpage> <lpage>496</lpage> <pub-id pub-id-type="pmid">8931275</pub-id></mixed-citation></ref><ref id="pone.0016308-Stopfer1"><label>15</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Stopfer</surname><given-names>M</given-names></name><name><surname>Bhagavan</surname><given-names>S</given-names></name><name><surname>Smith</surname><given-names>BH</given-names></name><name><surname>Laurent</surname><given-names>G</given-names></name></person-group> <year>1997</year> <article-title>Impaired odour discrimination on desynchronization of odour-encoding neural assemblies.</article-title> <source>Nature</source> <volume>390</volume> <fpage>70</fpage> <lpage>74</lpage> <pub-id pub-id-type="pmid">9363891</pub-id></mixed-citation></ref><ref id="pone.0016308-Getz1"><label>16</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Getz</surname><given-names>WM</given-names></name><name><surname>Lutz</surname><given-names>A</given-names></name></person-group> <year>1999</year> <article-title>A neural network model of general olfactory coding in the insect antennal lobe.</article-title> <source>Chem Senses</source> <volume>24</volume> <fpage>351</fpage> <lpage>372</lpage> <pub-id pub-id-type="pmid">10480672</pub-id></mixed-citation></ref><ref id="pone.0016308-Bazhenov1"><label>17</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Bazhenov</surname><given-names>M</given-names></name><name><surname>Stopfer</surname><given-names>M</given-names></name><name><surname>Rabinovich</surname><given-names>M</given-names></name><name><surname>Huerta</surname><given-names>R</given-names></name><name><surname>Abarbanel</surname><given-names>HD</given-names></name><etal></etal></person-group> <year>2001</year> <article-title>Model of transient oscillatory synchronization in the locust antennal lobe.</article-title> <source>Neuron</source> <volume>30</volume> <fpage>553</fpage> <lpage>567</lpage> <pub-id pub-id-type="pmid">11395014</pub-id></mixed-citation></ref><ref id="pone.0016308-Bazhenov2"><label>18</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Bazhenov</surname><given-names>M</given-names></name><name><surname>Stopfer</surname><given-names>M</given-names></name><name><surname>Rabinovich</surname><given-names>M</given-names></name><name><surname>Abarbanel</surname><given-names>HD</given-names></name><name><surname>Sejnowski</surname><given-names>TJ</given-names></name><etal></etal></person-group> <year>2001</year> <article-title>Model of cellular and network mechanisms for odor-evoked temporal patterning in the locust antennal lobe.</article-title> <source>Neuron</source> <volume>30</volume> <fpage>569</fpage> <lpage>81</lpage> <pub-id pub-id-type="pmid">11395015</pub-id></mixed-citation></ref><ref id="pone.0016308-Laurent2"><label>19</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Laurent</surname><given-names>G</given-names></name><name><surname>Stopfer</surname><given-names>M</given-names></name><name><surname>Friedrich</surname><given-names>RW</given-names></name><name><surname>Rabinovich</surname><given-names>MI</given-names></name><name><surname>Volkovskii</surname><given-names>A</given-names></name><etal></etal></person-group> <year>2001</year> <article-title>Odor encoding as an active, dynamical process: experiments, computation, and theory.</article-title> <source>Annu Rev Neurosci</source> <volume>24</volume> <fpage>263</fpage> <lpage>297</lpage> <pub-id pub-id-type="pmid">11283312</pub-id></mixed-citation></ref><ref id="pone.0016308-Rabinovich1"><label>20</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Rabinovich</surname><given-names>M</given-names></name><name><surname>Volkovskii</surname><given-names>A</given-names></name><name><surname>Lecanda</surname><given-names>P</given-names></name><name><surname>Huerta</surname><given-names>R</given-names></name><name><surname>Abarbanel</surname><given-names>HD</given-names></name><etal></etal></person-group> <year>2001</year> <article-title>Dynamical encoding by networks of competing neuron groups: winnerless competition.</article-title> <source>Phys Rev Lett</source> <volume>87</volume> <fpage>068102</fpage> <pub-id pub-id-type="pmid">11497865</pub-id></mixed-citation></ref><ref id="pone.0016308-Kwok1"><label>21</label><mixed-citation publication-type="other"> <person-group person-group-type="author"><name><surname>Kwok</surname><given-names>YC</given-names></name></person-group> <year>2007</year> <article-title>Encoding of odour blends in the moth antennal lobe.</article-title> <comment>PhD thesis, Department of Engineering University of Leicester</comment> </mixed-citation></ref><ref id="pone.0016308-Hodgkin1"><label>22</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Hodgkin</surname><given-names>AL</given-names></name><name><surname>Huxley</surname><given-names>AF</given-names></name><name><surname>Katz</surname><given-names>B</given-names></name></person-group> <year>1952</year> <article-title>Measurement of current-voltage relations in the membrane of the giant axon of Loligo.</article-title> <source>J Physiol</source> <volume>116</volume> <fpage>424</fpage> <lpage>448</lpage> <pub-id pub-id-type="pmid">14946712</pub-id></mixed-citation></ref><ref id="pone.0016308-Traub1"><label>23</label><mixed-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Traub</surname><given-names>RD</given-names></name><name><surname>Miles</surname><given-names>R</given-names></name></person-group> <year>1991</year> <source>Neuronal networks of the hippocampus</source> <publisher-name>Cambridge University Press</publisher-name> <size units="page"></size> </mixed-citation></ref><ref id="pone.0016308-Destexhe1"><label>24</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Destexhe</surname><given-names>A</given-names></name><name><surname>Mainen</surname><given-names>ZF</given-names></name><name><surname>Sejnowski</surname><given-names>TJ</given-names></name></person-group> <year>1994</year> <article-title>Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism.</article-title> <source>J Comp Neurosci</source> <volume>1</volume> <fpage>195</fpage> <lpage>230</lpage> </mixed-citation></ref><ref id="pone.0016308-Destexhe2"><label>25</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Destexhe</surname><given-names>A</given-names></name><name><surname>Mainen</surname><given-names>ZF</given-names></name><name><surname>Sejnowski</surname><given-names>TJ</given-names></name></person-group> <year>1994</year> <article-title>An efficient method for computing synaptic conductances based on a kinetic model of receptor binding.</article-title> <source>Neural Comput</source> <volume>6</volume> <fpage>14</fpage> <lpage>18</lpage> </mixed-citation></ref><ref id="pone.0016308-Schneider1"><label>26</label><mixed-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Schneider</surname><given-names>D</given-names></name></person-group> <year>1970</year> <article-title>Olfactory receptors for the sexual attractant (bombykol) the silk moth.</article-title> <person-group person-group-type="editor"><name><surname>Schmitt</surname><given-names>FO</given-names></name></person-group> <source>The Neurosciences Second Study Program</source> <publisher-name>Rockefeller Univ Press</publisher-name> <fpage>511</fpage> <lpage>518</lpage> </mixed-citation></ref><ref id="pone.0016308-Kaissling1"><label>27</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kaissling</surname><given-names>KE</given-names></name></person-group> <year>1996</year> <article-title>Peripheral mechanisms of pheromone reception in moths.</article-title> <source>Chem Senses</source> <volume>21</volume> <fpage>257</fpage> <lpage>268</lpage> <pub-id pub-id-type="pmid">8670704</pub-id></mixed-citation></ref><ref id="pone.0016308-Sanes1"><label>28</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Sanes</surname><given-names>JR</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1976</year> <article-title>Structure and development of antennae in a moth, Manduca sexta.</article-title> <source>Dev Biol</source> <volume>51</volume> <fpage>280</fpage> <lpage>299</lpage> <pub-id pub-id-type="pmid">955260</pub-id></mixed-citation></ref><ref id="pone.0016308-Matsumoto1"><label>29</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Matsumoto</surname><given-names>SG</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1981</year> <article-title>Olfactory mechanisms in the moth Manduca sexta: Response characteristics and morphology of central neurons in the antennal lobes.</article-title> <source>Proc Royal Soc London Series B Biol Sci</source> <volume>213</volume> <fpage>249</fpage> <lpage>277</lpage> </mixed-citation></ref><ref id="pone.0016308-Homberg1"><label>30</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Homberg</surname><given-names>U</given-names></name><name><surname>Montague</surname><given-names>RA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1988</year> <article-title>Anatomy of antenno-cerebral pathways in the brain of the sphinx moth Manduca sexta.</article-title> <source>Cell Tissue Res</source> <volume>254</volume> <fpage>255</fpage> <lpage>281</lpage> <pub-id pub-id-type="pmid">3197087</pub-id></mixed-citation></ref><ref id="pone.0016308-Christensen3"><label>31</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1987</year> <article-title>Male-specific, sex pheromone-selective projection neurons in the antennal lobes of the moth Manduca sexta.</article-title> <source>J Comp Physiol [A]</source> <volume>160</volume> <fpage>553</fpage> <lpage>569</lpage> </mixed-citation></ref><ref id="pone.0016308-Schfer1"><label>32</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Schäfer</surname><given-names>S</given-names></name><name><surname>Bicker</surname><given-names>G</given-names></name></person-group> <year>1986</year> <article-title>Distribution of GABA-like immunoreactivity in the brain of the honeybee.</article-title> <source>J Comp Neurol</source> <volume>246</volume> <fpage>287</fpage> <lpage>300</lpage> <pub-id pub-id-type="pmid">3700720</pub-id></mixed-citation></ref><ref id="pone.0016308-Waldrop1"><label>33</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Waldrop</surname><given-names>B</given-names></name><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1987</year> <article-title>GABA-mediated synaptic inhibition of projection neurons in the antennal lobes of the sphinx moth, Manduca sexta.</article-title> <source>J Comp Physiol [A]</source> <volume>161</volume> <fpage>23</fpage> <lpage>32</lpage> </mixed-citation></ref><ref id="pone.0016308-Homberg2"><label>34</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Homberg</surname><given-names>U</given-names></name><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1989</year> <article-title>Structure and function of the deutocerebrum in insects.</article-title> <source>Annu Rev Entomol</source> <volume>34</volume> <fpage>477</fpage> <lpage>501</lpage> <pub-id pub-id-type="pmid">2648971</pub-id></mixed-citation></ref><ref id="pone.0016308-Olsen1"><label>35</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Olsen</surname><given-names>SR</given-names></name><name><surname>Wilson</surname><given-names>RI</given-names></name></person-group> <year>2008</year> <article-title>Lateral presynaptic inhibition mediates gain control in an olfactory circuit.</article-title> <source>Nature</source> <volume>452</volume> <fpage>956</fpage> <lpage>960</lpage> <pub-id pub-id-type="pmid">18344978</pub-id></mixed-citation></ref><ref id="pone.0016308-Olsen2"><label>36</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Olsen</surname><given-names>SR</given-names></name><name><surname>Bhandawat</surname><given-names>V</given-names></name><name><surname>Wilson</surname><given-names>RI</given-names></name></person-group> <year>2007</year> <article-title>Excitatory interactions between olfactory processing channels in the Drosophila antennal lobe.</article-title> <source>Neuron</source> <volume>54</volume> <fpage>89</fpage> <lpage>103</lpage> <pub-id pub-id-type="pmid">17408580</pub-id></mixed-citation></ref><ref id="pone.0016308-AvRon1"><label>37</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Av-Ron</surname><given-names>E</given-names></name><name><surname>Rospars</surname><given-names>JP</given-names></name></person-group> <year>1995</year> <article-title>Modeling insect olfactory neuron signaling by a network utilizing disinhibition.</article-title> <source>BioSystems</source> <volume>36</volume> <fpage>101</fpage> <lpage>108</lpage> <pub-id pub-id-type="pmid">8573691</pub-id></mixed-citation></ref><ref id="pone.0016308-AvRon2"><label>38</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Av-Ron</surname><given-names>E</given-names></name><name><surname>Vibert</surname><given-names>JF</given-names></name></person-group> <year>1996</year> <article-title>A model for temporal and intensity coding in insect olfaction by a network of inhibitory neurons.</article-title> <source>Biosystems</source> <volume>39</volume> <fpage>241</fpage> <lpage>250</lpage> <pub-id pub-id-type="pmid">8894125</pub-id></mixed-citation></ref><ref id="pone.0016308-Husch1"><label>39</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Husch</surname><given-names>A</given-names></name><name><surname>Hess</surname><given-names>S</given-names></name><name><surname>Kloppenburg</surname><given-names>P</given-names></name></person-group> <year>2008</year> <article-title>Functional parameters of voltage-activated Ca<sup>2+</sup> currents from olfactory interneurons in the antennal lobe of Periplaneta americana.</article-title> <source>J Neurophysiol</source> <volume>99</volume> <fpage>320</fpage> <lpage>332</lpage> <pub-id pub-id-type="pmid">18003883</pub-id></mixed-citation></ref><ref id="pone.0016308-Husch2"><label>40</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Husch</surname><given-names>A</given-names></name><name><surname>Paehler</surname><given-names>M</given-names></name><name><surname>Fusca</surname><given-names>D</given-names></name><name><surname>Paeger</surname><given-names>L</given-names></name><name><surname>Kloppenburg</surname><given-names>P</given-names></name></person-group> <year>2009</year> <article-title>Distinct electrophysiological properties in subtypes of nonspiking olfactory local interneurons correlate with their cell type-specific Ca<sup>2+</sup> current profiles.</article-title> <source>J Neurophysio.l</source> <volume>102</volume> <fpage>2834</fpage> <lpage>2845</lpage> </mixed-citation></ref><ref id="pone.0016308-Husch3"><label>41</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Husch</surname><given-names>A</given-names></name><name><surname>Paehler</surname><given-names>M</given-names></name><name><surname>Fusca</surname><given-names>D</given-names></name><name><surname>Paeger</surname><given-names>L</given-names></name><name><surname>Kloppenburg</surname><given-names>P</given-names></name></person-group> <year>2009</year> <article-title>Calcium current diversity in physiologically different local interneuron types of the antennal lobe.</article-title> <source>J Neurosci</source> <volume>29</volume> <fpage>716</fpage> <lpage>726</lpage> <pub-id pub-id-type="pmid">19158298</pub-id></mixed-citation></ref><ref id="pone.0016308-Press1"><label>42</label><mixed-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Press</surname><given-names>WH</given-names></name><name><surname>Teukolsky</surname><given-names>SA</given-names></name><name><surname>Vetterling</surname><given-names>WT</given-names></name><name><surname>Flannery</surname><given-names>BP</given-names></name></person-group> <year>1992</year> <source>Numerical recipes in C (2nd ed.): The art of scientific computing</source> <publisher-name>Cambridge University Press</publisher-name> <size units="page"></size> </mixed-citation></ref><ref id="pone.0016308-Vickers1"><label>43</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Vickers</surname><given-names>NJ</given-names></name><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Mustaparta</surname><given-names>H</given-names></name><name><surname>Baker</surname><given-names>TC</given-names></name></person-group> <year>1991</year> <article-title>Chemical communication in heliothine moths: III. Flight behavior of male Helicoverpa zea and Heliothis virescens in response to varying ratios of intra- and interspecific sex Pheromone components 1991.</article-title> <source>J Compar Physiol A</source> <volume>169</volume> <fpage>275</fpage> <lpage>280</lpage> </mixed-citation></ref><ref id="pone.0016308-Vickers2"><label>44</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Vickers</surname><given-names>NJ</given-names></name><name><surname>Christensen</surname><given-names>TA</given-names></name><name><surname>Hildebrand</surname><given-names>JG</given-names></name></person-group> <year>1998</year> <article-title>Combinatorial odor discrimination in the brain: attractive and antagonist odor blends are represented in distinct combinations of uniquely identifiable glomeruli.</article-title> <source>J Comp Neurol</source> <volume>400</volume> <fpage>35</fpage> <lpage>56</lpage> <pub-id pub-id-type="pmid">9762865</pub-id></mixed-citation></ref><ref id="pone.0016308-Hansson2"><label>45</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Hansson</surname><given-names>BS</given-names></name><name><surname>Almaas</surname><given-names>TJ</given-names></name><name><surname>Anton</surname><given-names>S</given-names></name></person-group> <year>1995</year> <article-title>Chemical communication in heliothine moths: V. Antennal lobe projection patterns of pheromone-detecting olfactory receptorneurons in the male Heliothis virescens (Lepidoptera: Noctuidae).</article-title> <source>J Compar Physiol A</source> <volume>177</volume> <fpage>535</fpage> <lpage>543</lpage> </mixed-citation></ref><ref id="pone.0016308-Ditzen1"><label>46</label><mixed-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Ditzen</surname><given-names>M</given-names></name><name><surname>Evers</surname><given-names>JF</given-names></name><name><surname>Galizia</surname><given-names>CG</given-names></name></person-group> <year>2003</year> <article-title>Odor similarity does not influence the time needed for odor processing.</article-title> <source>Chemical Senses</source> <volume>28</volume> <fpage>781</fpage> <lpage>789</lpage> <pub-id pub-id-type="pmid">14654446</pub-id></mixed-citation></ref><ref id="pone.0016308-Todd1"><label>47</label><mixed-citation publication-type="other"> <person-group person-group-type="author"><name><surname>Todd</surname><given-names>JL</given-names></name><name><surname>Baker</surname><given-names>TC</given-names></name></person-group> <year>1999</year> <article-title>Function of Peripheral Olfactory Organs.</article-title> <person-group person-group-type="editor"><name><surname>Hansson</surname><given-names>BS</given-names></name></person-group> <source>Insect Olfaction</source> <fpage>67</fpage> <lpage>96</lpage> <comment>Springer, Heidelberg, Berlin</comment> </mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Lorraine/explor/InforLorV4/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000086 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000086 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Lorraine
|area= InforLorV4
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |