diff --git a/plos2015.bst b/plos2015.bst new file mode 100644 index 0000000..f696b8c --- /dev/null +++ b/plos2015.bst @@ -0,0 +1,1883 @@ +%% +%% +%% +%% +%% This `plos2015.bst' file is intended for use in PLOS submissions to +%% correctly compile the bibliography according to the PLOS manuscript +%% guidelines updated January 2015. +%% +%% +%% This is derived from the `vancouver.bst' bibliographic style file (for +%% LaTeX/BibTeX), generated with the docstrip utility and modified manually +%% to meet the +%% ``Uniform Requirements for Manuscripts Submitted to Biomedical Journals'' +%% as published in N Engl J Med 1997;336:309-315. +%% (also known as the Vancouver style) +%% This specification may be found on the web page of the +%% International Committe of Medical Journal Editors: +%% +%% http://www.icmje.org +%% +%%------------------------------------------------------------------- +%% +%% Copyright 2004 Folkert van der Beek +%% +%% This work may be distributed and/or modified under the +%% conditions of the LaTeX Project Public License, either version 1.3 +%% of this license or (at your option) any later version. +%% The latest version of this license is in +%% http://www.latex-project.org/lppl.txt +%% and version 1.3 or later is part of all distributions of LaTeX +%% version 2005/12/01 or later. +%% +%% This work has the LPPL maintenance status `maintained'. +%% +%% The Current Maintainer of this work is Folkert van der Beek. +%% +%% Complaints, suggestions and comments may be sent to +%% +%% Folkert van der Beek +%% +%%------------------------------------------------------------------- +%% +%% This bibliography style file is intended for texts in ENGLISH +%% This is a numerical citation style, and as such is standard LaTeX. +%% It requires no extra package to interface to the main text. +%% The form of the \bibitem entries is +%% \bibitem{key}... +%% Usage of \cite is as follows: +%% \cite{key} ==>> [#] +%% \cite[chap. 2]{key} ==>> [#, chap. 2] +%% where # is a number determined by the ordering in the reference list. +%% The order in the reference list is that by which the works were originally +%% cited in the text, or that in the database. + % +%% To change the reference numbering system from [1] to 1, +%% put the following code in the preamble: +%% \makeatletter % Reference list option change +%% \renewcommand\@biblabel[1]{#1} % from [1] to 1 +%% \makeatother % +%% +%%--------------------------------------------------------------------- + +%% List of all possible fields +ENTRY + { address + assignee % for patents + author + booktitle % for articles in books + chapter % for incollection, esp. internet documents + cartographer % for maps + day + edition + editor + howpublished + institution % for technical reports + inventor % for patents + journal + key + month + note + number + organization + pages + part + publisher + school + series + title + type + volume + word + year + eprint % urlbst + doi % urlbst + url % urlbst + lastchecked % urlbst + updated % urlbst + } + {} + { label } +%% Declaration of integer variables +INTEGERS { output.state before.all mid.sentence after.sentence after.block } +STRINGS { urlintro doiprefix } % urlbst... +INTEGERS { hrefform addeprints adddoiresolver } +% Following constants may be adjusted by hand, if desired +FUNCTION {init.config.constants} +{ + "Available from: " 'urlintro := % prefix before URL + "doi:" 'doiprefix := % text prefix printed before DOI ref + #0 'addeprints := % 0=no eprints; 1=include eprints + #0 'adddoiresolver := % 0=no DOI resolver; 1=include it + #0 'hrefform := % 0=no crossrefs; 1=hypertex xrefs; 2=hyperref refs +} +INTEGERS { + bracket.state + outside.brackets + open.brackets + within.brackets + close.brackets +} +% ...urlbst to here +FUNCTION {init.state.consts} +{ #0 'outside.brackets := % urlbst + #1 'open.brackets := + #2 'within.brackets := + #3 'close.brackets := + + #0 'before.all := + #1 'mid.sentence := + #2 'after.sentence := + #3 'after.block := +} +%% Declaration of string variables +STRINGS { s t} + +% urlbst +FUNCTION {output.nonnull.original} +{ 's := + output.state mid.sentence = + { ". " * write$ } + { output.state after.block = + { add.period$ write$ + newline$ + "\newblock " write$ + } + { output.state before.all = + 'write$ + { add.period$ " " * write$ } + if$ + } + if$ + mid.sentence 'output.state := + } + if$ + s +} + +% urlbst... +FUNCTION {output.nonnull} +{ % Save the thing we've been asked to output + 's := + % If the bracket-state is close.brackets, then add a close-bracket to + % what is currently at the top of the stack, and set bracket.state + % to outside.brackets + bracket.state close.brackets = + { "]" * + outside.brackets 'bracket.state := + } + 'skip$ + if$ + bracket.state outside.brackets = + { % We're outside all brackets -- this is the normal situation. + % Write out what's currently at the top of the stack, using the + % original output.nonnull function. + s + output.nonnull.original + } + { % Still in brackets. Add open-bracket or (continuation) comma, add the + % new text (in s) to the top of the stack, and move to the close-brackets + % state, ready for next time (unless inbrackets resets it). If we come + % into this branch, then output.state is carefully undisturbed. + bracket.state open.brackets = + { " [" * } + { ", " * } % bracket.state will be within.brackets + if$ + s * + close.brackets 'bracket.state := + } + if$ +} + +% Call this function just before adding something which should be presented in +% brackets. bracket.state is handled specially within output.nonnull. +FUNCTION {inbrackets} +{ bracket.state close.brackets = + { within.brackets 'bracket.state := } % reset the state: not open nor closed + { open.brackets 'bracket.state := } + if$ +} + +FUNCTION {format.lastchecked} +{ lastchecked empty$ + { "" } + { updated empty$ + { inbrackets "cited " lastchecked * } + { inbrackets "updated " updated * "; cited " * lastchecked * } + if$ + } + if$ +} +% ...urlbst to here + +FUNCTION {output} +{ duplicate$ empty$ + 'pop$ + 'output.nonnull + if$ +} + +FUNCTION {output.check} +{ 't := + duplicate$ empty$ + { pop$ "empty " t * " in " * cite$ * warning$ } + 'output.nonnull + if$ +} + +FUNCTION {fin.entry} +{ + bracket.state close.brackets = % urlbst + { "]" * } + 'skip$ + if$ + add.period$ + write$ + newline$ +} + +FUNCTION {new.block} +{ output.state before.all = + 'skip$ + { after.block 'output.state := } + if$ +} + +FUNCTION {new.sentence} +{ output.state after.block = + 'skip$ + { output.state before.all = + 'skip$ + { after.sentence 'output.state := } + if$ + } + if$ +} + +FUNCTION {add.blank} +{ " " * before.all 'output.state := +} + +FUNCTION {no.blank.or.punct} +{ "" * before.all 'output.state := +} + +FUNCTION {add.semicolon} +{ + ";" * + no.blank.or.punct +} + +FUNCTION {date.block} +{ + "." * + no.blank.or.punct +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% LOGICAL `NOT', `AND', AND `OR' % +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Logical 'not': +% If the first element on the stack is A then this function +% does the following: +% push { #0 } +% push { #1 } +% So now the first 3 elements of the stack are +% { #1 } { #0 } A +% The first 3 are popped and subjected to 'if': +% If A > 0 then { #0 } is executed, else { #1 } is executed: +% if A > 0 +% then 0 +% else 1 +% So consider integers as logicals, where 1 = true and 0 = false, +% then this does +% (if A then false else true) +% which is a logical 'not'. + +FUNCTION {not} +{ { #0 } + { #1 } + if$ +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Logical 'and': +% If the first 2 elements on the stack are A B +% then this function does the following: +% push 'skip$ +% push { pop$ #0 } +% So now first 4 elements are +% { pop$ #0 } 'skip$ A B +% The first 3 are popped and subjected to 'if' (B is on top of +% the stack): +% If A > 0 then 'skip$ is executed, else { pop$ #0 } is executed: +% if A > 0 +% then (B stays on top of stack) +% else (B is popped and #0 is pushed) +% So consider integers as logicals, where 1 = true and 0 = false, +% then this does +% (if A then B else false) +% which is a logical 'and'. + +FUNCTION {and} +{ 'skip$ + { pop$ #0 } + if$ +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Logical 'or': +% If the first 2 elements on the stack are A B +% then this function does the following: +% push { pop$ #1 } +% push 'skip$ +% So now first 4 elements are +% 'skip$ { pop$ #1 } A B +% The first 3 are popped and subjected to 'if' (B is on top of +% the stack): +% If A > 0 then { pop$ #1 } is executed, else 'skip$ is executed: +% if A > 0 +% then (B is popped and #1 is pushed) +% else (B stays on top of stack) +% So consider integers as logicals, where 1 = true and 0 = false, +% then this does +% (if A then true else B) +% which is a logical 'or'. + +FUNCTION {or} +{ { pop$ #1 } + 'skip$ + if$ +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% GENERAL PURPOSE FUNCTIONS FOR FORMATTING % +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% issues warning if field is empty +% call with +% "field" field warning.if.empty +% Note that the first field must be between quotes +% because it is the fieldname for use in the warning message. +% + +FUNCTION {warning.if.empty} +{ empty$ + { "No " swap$ * " in " * cite$ * warning$ } + { pop$ } + if$ +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + % + % encloses string in pre- and postfix string + % call with + % prefix postfix S enclose.check + % delivers empty string if S empty + % +FUNCTION {enclose.check} +{ duplicate$ empty$ + { pop$ pop$ pop$ + "" + } + { swap$ * * } + if$ +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% +% emphasizes top of stack +% call with +% string" emphasize.check +% + +FUNCTION {emphasize.check} +{ "\Bem{" swap$ + "}" swap$ + enclose.check +} + + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + % + % brackets top of stack + % call with + % "string" bracket.check + % +FUNCTION {bracket.check} +{ "[" swap$ + "]" swap$ + enclose.check +} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + % + % parenthesizes top of stack + % call with + % "string" parenthesize + % +FUNCTION {parenthesize.check} +{ "(" swap$ + ")" swap$ + enclose.check +} + +STRINGS {z} + +FUNCTION {remove.dots} +{ 'z := % expects string on top of the stack, pops the string and assigns it to variable z + "" % push empty string + { z empty$ not } % returns 0 if variable z is empty + { z #1 #1 substring$ % push the first character of variable z + z #2 global.max$ substring$ 'z := % assigns the 2nd to last character of variable z to variable z + duplicate$ "\" = % pushes 1 if the last character is "\", otherwise 0 + { * % concatenates the last 2 literals + z #1 #1 substring$ % push the first character of variable z + z #2 global.max$ substring$ 'z := % assigns the 2nd to last character of variable z to variable z + * % concatenates the last 2 literals, i.e. every character, even a dot, following a "\" will be printed + } + { duplicate$ "." = % pushes 1 if the last character is ".", otherwise 0 + 'pop$ % pushes the pop$ function + { * } % concatenates the last 2 literals + if$ % pops the last character if it is a dot, otherwise concatenates it with the string on top of the stack + } + if$ + } + while$ +} + +INTEGERS {l} +FUNCTION{string.length} +{ + #1 'l := + { duplicate$ duplicate$ #1 l substring$ = not } + { l #1 + 'l := } + while$ + pop$ l +} + +STRINGS {replace find text} +INTEGERS {find_length} +FUNCTION {find.replace} +{ + 'replace := + 'find := + 'text := + find string.length 'find_length := + "" + { text empty$ not } + { text #1 find_length substring$ find = + { + replace * + text #1 find_length + global.max$ substring$ 'text := + } + { text #1 #1 substring$ * + text #2 global.max$ substring$ 'text := + } + if$ + } + while$ +} + +FUNCTION {new.block.checka} +{ empty$ + 'skip$ + 'new.block + if$ +} + +FUNCTION {new.block.checkb} +{ empty$ + swap$ empty$ + and + 'skip$ + 'new.block + if$ +} + +FUNCTION {new.sentence.checka} +{ empty$ + 'skip$ + 'new.sentence + if$ +} + +FUNCTION {new.sentence.checkb} +{ empty$ + swap$ empty$ + and + 'skip$ + 'new.sentence + if$ +} + +FUNCTION {field.or.null} +{ duplicate$ empty$ + { pop$ "" } + 'skip$ + if$ +} + +FUNCTION {emphasize} +{ skip$ } + +FUNCTION {tie.or.space.prefix} +{ duplicate$ text.length$ #3 < + { "~" } + { " " } + if$ + swap$ +} + +FUNCTION {capitalize} +{ "u" change.case$ "t" change.case$ } + +FUNCTION {space.word} +{ " " swap$ * " " * } + + % Here are the language-specific definitions for explicit words. + % Each function has a name bbl.xxx where xxx is the English word. + % The language selected here is ENGLISH + +FUNCTION {bbl.and} +{ "and"} + +FUNCTION {bbl.etal} +{ "et~al." } + +FUNCTION {bbl.editors} +{ "editors" } + +FUNCTION {bbl.editor} +{ "editor" } + +FUNCTION {bbl.cartographers} +{ "cartographers" } + +FUNCTION {bbl.cartographer} +{ "cartographer" } + +FUNCTION {bbl.inventors} +{ "inventors" } + +FUNCTION {bbl.inventor} +{ "inventor" } + +FUNCTION {bbl.assignees} +{ "assignees" } + +FUNCTION {bbl.assignee} +{ "assignee" } + +FUNCTION {bbl.edby} +{ "edited by" } + +FUNCTION {bbl.edition} +{ "ed." } + +FUNCTION {bbl.volume} +{ "vol." } + +FUNCTION {bbl.of} +{ "of" } + +FUNCTION {bbl.number} +{ "no." } + +FUNCTION {bbl.nr} +{ "no." } + +FUNCTION {bbl.in} +{ "in" } + +FUNCTION {bbl.pages} +{ " p." } + +FUNCTION {bbl.page} +{ " p." } + +FUNCTION {bbl.chapter} +{ "chap." } + +FUNCTION {bbl.techrep} +{ "Tech. Rep." } + +FUNCTION {bbl.mthesis} +{ "Master's thesis" } + +FUNCTION {bbl.phdthesis} +{ "Ph.D. thesis" } + +FUNCTION {bbl.first} +{ "1st" } + +FUNCTION {bbl.second} +{ "2nd" } + +FUNCTION {bbl.third} +{ "3rd" } + +FUNCTION {bbl.fourth} +{ "4th" } + +FUNCTION {bbl.fifth} +{ "5th" } + +FUNCTION {bbl.st} +{ "st" } + +FUNCTION {bbl.nd} +{ "nd" } + +FUNCTION {bbl.rd} +{ "rd" } + +FUNCTION {bbl.th} +{ "th" } + +MACRO {jan} {"Jan."} + +MACRO {feb} {"Feb."} + +MACRO {mar} {"Mar."} + +MACRO {apr} {"Apr."} + +MACRO {may} {"May"} + +MACRO {jun} {"Jun."} + +MACRO {jul} {"Jul."} + +MACRO {aug} {"Aug."} + +MACRO {sep} {"Sep."} + +MACRO {oct} {"Oct."} + +MACRO {nov} {"Nov."} + +MACRO {dec} {"Dec."} + +FUNCTION {eng.ord} +{ duplicate$ "1" swap$ * + #-2 #1 substring$ "1" = + { bbl.th * } + { duplicate$ #-1 #1 substring$ + duplicate$ "1" = + { pop$ bbl.st * } + { duplicate$ "2" = + { pop$ bbl.nd * } + { "3" = + { bbl.rd * } + { bbl.th * } + if$ + } + if$ + } + if$ + } + if$ +} + +FUNCTION {bibinfo.check} +{ swap$ + duplicate$ missing$ + { + pop$ pop$ + "" + } + { duplicate$ empty$ + { + swap$ pop$ + } + { swap$ + pop$ + } + if$ + } + if$ +} + +FUNCTION {bibinfo.warn} +{ swap$ + duplicate$ missing$ + { + swap$ "missing " swap$ * " in " * cite$ * warning$ pop$ + "" + } + { duplicate$ empty$ + { + swap$ "empty " swap$ * " in " * cite$ * warning$ + } + { swap$ + pop$ + } + if$ + } + if$ +} + +STRINGS { bibinfo} +INTEGERS { nameptr namesleft numnames } + +FUNCTION {format.names} +{ 'bibinfo := + duplicate$ empty$ 'skip$ { + "." ". " find.replace 's := + "" 't := + #1 'nameptr := + s num.names$ 'numnames := + numnames 'namesleft := + { namesleft #0 > } + { s nameptr + "{vv~}{ll}{ f{}}{ jj}" + format.name$ + remove.dots + bibinfo bibinfo.check + 't := + nameptr #1 > + { + nameptr #6 + #1 + = + numnames #6 + > and + { "others" 't := + #1 'namesleft := } + 'skip$ + if$ + namesleft #1 > + { ", " * t * } + { + "," * + s nameptr "{ll}" format.name$ duplicate$ "others" = + { 't := } + { pop$ } + if$ + t "others" = + { + " " * bbl.etal * + } + { " " * t * } + if$ + } + if$ + } + 't + if$ + nameptr #1 + 'nameptr := + namesleft #1 - 'namesleft := + } + while$ + } if$ +} + +FUNCTION {format.names.org} +{ 'bibinfo := + duplicate$ empty$ 'skip$ { + 's := + "" 't := + #1 'nameptr := + s num.names$ 'numnames := + numnames 'namesleft := + { namesleft #0 > } + { s nameptr + "{ff~}{vv~}{ll}" + format.name$ + bibinfo bibinfo.check + 't := + nameptr #1 > + { + namesleft #1 > + { "; " * t * } + { + ";" * + s nameptr "{ll}" format.name$ duplicate$ "others" = + { 't := } + { pop$ } + if$ + t "others" = + { + " " * bbl.etal * + } + { " " * t * } + if$ + } + if$ + } + 't + if$ + nameptr #1 + 'nameptr := + namesleft #1 - 'namesleft := + } + while$ + } if$ +} + +FUNCTION {format.names.ed} +{ + format.names +} + +FUNCTION {format.authors} +{ + author "author" format.names + %%"." " " "author" find.replace format.names +} + +FUNCTION {format.organizations} +{ organization "organization" format.names.org +} + +FUNCTION {get.bbl.editor} +{ editor num.names$ #1 > 'bbl.editors 'bbl.editor if$ } + +FUNCTION {get.bbl.cartographer} +{ cartographer num.names$ #1 > 'bbl.cartographers 'bbl.cartographer if$ } + +FUNCTION {get.bbl.inventor} +{ inventor num.names$ #1 > 'bbl.inventors 'bbl.inventor if$ } + +FUNCTION {get.bbl.assignee} +{ assignee num.names$ #1 > 'bbl.assignees 'bbl.assignee if$ } + +FUNCTION {format.editors} +{ editor "editor" format.names duplicate$ empty$ 'skip$ + { + "," * + " " * + get.bbl.editor + * + } + if$ +} + +FUNCTION {format.assignees} +{ assignee "assignee" format.names.org duplicate$ empty$ 'skip$ + { + "," * + " " * + get.bbl.assignee + * + } + if$ +} + +FUNCTION {format.cartographers} +{ cartographer "cartographer" format.names duplicate$ empty$ 'skip$ + { + "," * + " " * + get.bbl.cartographer + * + } + if$ +} + +FUNCTION {format.inventors} +{ inventor "inventor" format.names duplicate$ empty$ 'skip$ + { + "," * + " " * + get.bbl.inventor + * + } + if$ +} + +FUNCTION {format.note} +{ + note empty$ + { "" } + { note #1 #1 substring$ + duplicate$ "{" = + 'skip$ + { output.state mid.sentence = + { "l" } + { "u" } + if$ + change.case$ + } + if$ + note #2 global.max$ substring$ * "note" bibinfo.check + } + if$ +} + +FUNCTION {format.title} +{ title +%%duplicate$ empty$ 'skip$ +%% { "t" change.case$ } +%%if$ + "title" bibinfo.check +} + +FUNCTION {format.type} +{ type empty$ + 'skip$ + { inbrackets type } + %%{ add.blank "[" type * "]" * } + if$ +} + + +FUNCTION {output.bibitem} +{ outside.brackets 'bracket.state := % urlbst + newline$ + "\bibitem{" write$ + cite$ write$ + "}" write$ + newline$ + "" + before.all 'output.state := +} + +FUNCTION {n.dashify} +{ + 't := + "" + { t empty$ not } + { t #1 #1 substring$ "-" = + { t #1 #2 substring$ "--" = not + { "--" * + t #2 global.max$ substring$ 't := + } + { { t #1 #1 substring$ "-" = } + { "-" * + t #2 global.max$ substring$ 't := + } + while$ + } + if$ + } + { t #1 #1 substring$ * + t #2 global.max$ substring$ 't := + } + if$ + } + while$ +} + +FUNCTION {word.in} +{ bbl.in capitalize + ":" * + " " * } + +FUNCTION {format.journal.date} +{ + "" + duplicate$ empty$ + year "year" bibinfo.check duplicate$ empty$ + { + swap$ 'skip$ + { "there's a month but no year in " cite$ * warning$ } + if$ + * + } + { swap$ 'skip$ + { + " " * swap$ + } + if$ + * + remove.dots + } + if$ + duplicate$ empty$ + 'skip$ + { + before.all 'output.state := + after.sentence 'output.state := + } + if$ +} + +FUNCTION {format.date} +{ + no.blank.or.punct + ";" + duplicate$ empty$ + year "year" bibinfo.check duplicate$ empty$ + { swap$ 'skip$ + { "there's a month but no year in " cite$ * warning$ } + if$ + * + } + { swap$ 'skip$ + { + swap$ + " " * swap$ + } + if$ + * + } + if$ +} + +FUNCTION {format.btitle} +{ title "title" bibinfo.check + duplicate$ empty$ 'skip$ + { + } + if$ +} + +FUNCTION {either.or.check} +{ empty$ + 'pop$ + { "can't use both " swap$ * " fields in " * cite$ * warning$ } + if$ +} + +FUNCTION {format.bvolume} +{ volume empty$ + { "" } + { bbl.volume volume tie.or.space.prefix + "volume" bibinfo.check * * + series "series" bibinfo.check + duplicate$ empty$ 'pop$ + { swap$ bbl.of space.word * swap$ + emphasize * } + if$ + "volume and number" number either.or.check + } + if$ +} + +FUNCTION {format.number.series} +{ volume empty$ + { number empty$ + { series field.or.null } + { series empty$ + { number "number" bibinfo.check } + { output.state mid.sentence = + { bbl.number } + { bbl.number capitalize } + if$ + number tie.or.space.prefix "number" bibinfo.check * * + bbl.in space.word * + series "series" bibinfo.check * + } + if$ + } + if$ + } + { "" } + if$ +} + +FUNCTION {is.num} +{ chr.to.int$ + duplicate$ "0" chr.to.int$ < not + swap$ "9" chr.to.int$ > not and +} + +FUNCTION {extract.num} +{ duplicate$ 't := + "" 's := + { t empty$ not } + { t #1 #1 substring$ + t #2 global.max$ substring$ 't := + duplicate$ is.num + { s swap$ * 's := } + { pop$ "" 't := } + if$ + } + while$ + s empty$ + 'skip$ + { pop$ s } + if$ +} + +FUNCTION {convert.edition} +{ extract.num "l" change.case$ 's := + s "first" = s "1" = or + { bbl.first 't := } + { s "second" = s "2" = or + { bbl.second 't := } + { s "third" = s "3" = or + { bbl.third 't := } + { s "fourth" = s "4" = or + { bbl.fourth 't := } + { s "fifth" = s "5" = or + { bbl.fifth 't := } + { s #1 #1 substring$ is.num + { s eng.ord 't := } + { edition 't := } + if$ + } + if$ + } + if$ + } + if$ + } + if$ + } + if$ + t +} + +FUNCTION {format.edition} +{ edition duplicate$ empty$ 'skip$ + { + convert.edition + output.state mid.sentence = + { "l" } + { "t" } + if$ change.case$ + "edition" bibinfo.check + " " * bbl.edition * + } + if$ +} +INTEGERS { multiresult } +FUNCTION {multi.page.check} +{ 't := + #0 'multiresult := + { multiresult not + t empty$ not + and + } + { t #1 #1 substring$ + duplicate$ "-" = + swap$ duplicate$ "," = + swap$ "+" = + or or + { #1 'multiresult := } + { t #2 global.max$ substring$ 't := } + if$ + } + while$ + multiresult +} + +FUNCTION {format.pages} +{ pages duplicate$ empty$ 'skip$ + { duplicate$ multi.page.check + { + bbl.pages swap$ + n.dashify + } + { + bbl.page swap$ + } + if$ + tie.or.space.prefix + "pages" bibinfo.check + * * + } + if$ +} + +FUNCTION {format.journal.pages} +{ pages duplicate$ empty$ 'pop$ + { swap$ duplicate$ empty$ + { pop$ pop$ format.pages } + { + ":" * + swap$ + n.dashify + "pages" bibinfo.check + * + } + if$ + } + if$ +} + +FUNCTION {format.vol.num} +{ volume field.or.null + duplicate$ empty$ 'skip$ + { + "volume" bibinfo.check + } + if$ + number "number" bibinfo.check duplicate$ empty$ 'skip$ + { + swap$ duplicate$ empty$ + { "there's a number but no volume in " cite$ * warning$ } + 'skip$ + if$ + swap$ + "(" swap$ * ")" * + } + if$ * +} + +FUNCTION {format.vol.num.pages} +{ volume field.or.null + duplicate$ empty$ 'skip$ + { + "volume" bibinfo.check + } + if$ + number "number" bibinfo.check duplicate$ empty$ 'skip$ + { + swap$ duplicate$ empty$ + { "there's a number but no volume in " cite$ * warning$ } + 'skip$ + if$ + swap$ + "(" swap$ * ")" * + } + if$ * + format.journal.pages +} + +FUNCTION {format.chapter.pages} +{ chapter empty$ + 'format.pages + { type empty$ + { bbl.chapter } + { type "l" change.case$ + "type" bibinfo.check + } + if$ + chapter tie.or.space.prefix + "chapter" bibinfo.check + * * + pages empty$ + 'skip$ + { ", " * format.pages * } + if$ + } + if$ + } + + FUNCTION {format.booktitle} + { + booktitle "booktitle" bibinfo.check + } + + FUNCTION {format.in.ed.booktitle} + { format.booktitle duplicate$ empty$ 'skip$ + { + editor "editor" format.names.ed duplicate$ empty$ 'pop$ + { + "," * + " " * + get.bbl.editor + ". " * + * swap$ + * } + if$ + word.in swap$ * + } + if$ + } + + FUNCTION {format.in.ed.title} + { format.title duplicate$ empty$ 'skip$ + { + editor "editor" format.names.ed duplicate$ empty$ 'pop$ + { + "," * + " " * + get.bbl.editor + ". " * + * swap$ + * } + if$ + word.in swap$ * + } + if$ + } + + FUNCTION {empty.misc.check} + { author empty$ title empty$ howpublished empty$ + month empty$ year empty$ note empty$ + and and and and and + { "all relevant fields are empty in " cite$ * warning$ } + 'skip$ + if$ + } + FUNCTION {format.thesis.type} + { type duplicate$ empty$ + 'pop$ + { swap$ pop$ + "t" change.case$ "type" bibinfo.check + } + if$ +} +FUNCTION {format.tr.number} +{ + number "number" bibinfo.check + %%type duplicate$ empty$ + %%{ pop$ bbl.techrep } + %%'skip$ + %%if$ + %%"type" bibinfo.check + %%swap$ duplicate$ empty$ + %%{ pop$ "t" change.case$ } + %%{ tie.or.space.prefix * * } + %%if$ +} + +FUNCTION {format.org.or.pub} +{ 't := + "" + address empty$ t empty$ and + 'skip$ + { + address "address" bibinfo.check * + t empty$ + 'skip$ + { address empty$ + 'skip$ + { ": " * } + if$ + t * + } + if$ + } + if$ +} + +FUNCTION {format.publisher.address} +{ publisher "publisher" bibinfo.warn format.org.or.pub +} + +FUNCTION {format.organization.address} +{ organization "organization" bibinfo.check format.org.or.pub +} + +FUNCTION {format.institution.address} +{ institution "institution" bibinfo.check format.org.or.pub +} + + +% urlbst... +% Functions for making hypertext links. +% In all cases, the stack has (link-text href-url) +% +% make 'null' specials +FUNCTION {make.href.null} +{ + pop$ +} +% make hypertex specials +FUNCTION {make.href.hypertex} +{ + "\special {html: }" * swap$ * + "\special {html:}" * +} +% make hyperref specials +FUNCTION {make.href.hyperref} +{ + "\href {" swap$ * "} {" * swap$ * "}" * +} +FUNCTION {make.href} +{ hrefform #2 = + 'make.href.hyperref % hrefform = 2 + { hrefform #1 = + 'make.href.hypertex % hrefform = 1 + 'make.href.null % hrefform = 0 (or anything else) + if$ + } + if$ +} + +FUNCTION {format.url} +{ url empty$ + { "" } + { hrefform #1 = + { % special case -- add HyperTeX specials + urlintro "\url{" url * "}" * url make.href.hypertex * } + { urlintro "\url{" * url * "}" * } + if$ + } + if$ +} + +%FUNCTION {format.eprint} +%{ eprint empty$ +% { "" } +% { eprintprefix eprint * eprinturl eprint * make.href } +% if$ +%} + +FUNCTION {format.doi} +{ doi empty$ + { "" } + { + new.block + "doi:{" doi * "}" * + } + if$ +} + +% Output a URL. We can't use the more normal idiom (something like +% `format.url output'), because the `inbrackets' within +% format.lastchecked applies to everything between calls to `output', +% so that `format.url format.lastchecked * output' ends up with both +% the URL and the lastchecked in brackets. +FUNCTION {output.url} +{ url empty$ + 'skip$ + { new.block + format.url output + format.lastchecked output + } + if$ +} + +%FUNCTION {output.web.refs} +%{ +% new.block +% output.url +% addeprints eprint empty$ not and +% { format.eprint output.nonnull } +% 'skip$ +% if$ +% adddoiresolver doi empty$ not and +% { format.doi output.nonnull } +% 'skip$ +% if$ +% addeprints +% { eprint empty$ +% 'skip$ +% { format.eprint output.nonnull } +% if$ +% } +% 'skip$ +% if$ +%} + + + +% Webpage entry type. +% Title and url fields required; +% author, note, year, month, and lastchecked fields optional +STRINGS {database} +FUNCTION {webpage} +{ output.bibitem + author empty$ + { editor empty$ + 'skip$ % author and editor both optional + { format.editors output.nonnull } + if$ + } + { editor empty$ + { format.authors output.nonnull } + { "can't use both author and editor fields in " cite$ * warning$ } + if$ + } + if$ +% author empty$ +% 'skip$ +% { format.authors output.nonnull } +% if$ + new.block + format.title "title" output.check + journal empty$ + { + format.type "type" + publisher empty$ + 'skip$ + { format.publisher.address output } + if$ + "database on the Internet" 'database := + type database = + { format.journal.date "year" output.check } + { format.date "year" output.check } + if$ + lastchecked empty$ + 'skip$ + { format.lastchecked output } + if$ + new.block + part empty$ + 'skip$ + { part output } + if$ + pages empty$ + 'skip$ + { pages bracket.check output } + if$ + } + { journal + remove.dots + "journal" bibinfo.check + "journal" output.check + format.type "type" + format.journal.date "year" output.check + lastchecked empty$ + 'skip$ + { format.lastchecked output + ";" no.blank.or.punct output + } + if$ + no.blank.or.punct format.vol.num output + pages empty$ + 'skip$ + { ":" no.blank.or.punct output + no.blank.or.punct pages bracket.check output + } + if$ + new.block + } + if$ + format.url "url" output.check + new.block + note output + fin.entry +} +% ...urlbst to here + +FUNCTION {misc} +{ output.bibitem + format.authors "author" output.check +% format.editors "author and editor" output.check + format.title "title" output.check + type missing$ + { skip$ } + { format.type "type" } + %%{ inbrackets type output } + if$ + new.block +% format.publisher.address output + format.date "year" output.check + new.block +% format.note output + new.block + howpublished new.block.checka + howpublished "howpublished" bibinfo.check output + output.url % urlbst + fin.entry + empty.misc.check +} + +FUNCTION {article} +{ output.bibitem + format.authors "author" output.check + organization empty$ + 'skip$ + { author empty$ + { + format.organizations "organization" output.check + } + { + "; " * + no.blank.or.punct + format.organizations "organization" output.check + } + if$ + } + if$ + new.block + format.title "title" output.check + new.block + journal + remove.dots + "journal" bibinfo.check + "journal" output.check + format.journal.date "year" output.check + add.semicolon + format.vol.num.pages output + new.block + format.doi output + new.block +% format.note output +% output.url % urlbst + fin.entry +} + +FUNCTION {book} +{ output.bibitem + author empty$ + { editor empty$ + { format.organizations "organization" output.check } + { format.editors "author and editor" output.check } + if$ + } + { format.authors output.nonnull + "author and editor" editor either.or.check + } + if$ + new.block + format.btitle "title" output.check + format.bvolume output + new.block + format.edition output + new.sentence + author empty$ not + editor empty$ not + and + { format.editors "author and editor" output.check } + 'skip$ + if$ + format.number.series output + format.publisher.address output + format.date "year" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {booklet} +{ misc } + +FUNCTION {dictionary} +{ output.bibitem + format.booktitle "booktitle" output.check + format.bvolume output + new.block + format.edition output + new.sentence + format.publisher.address output + format.date "year" output.check + format.btitle "title" output.check + add.semicolon +% add.blank + format.pages "pages" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {inbook} +{ output.bibitem + format.authors "author" output.check + new.block + chapter "chapter" output.check + new.block + format.in.ed.title "title" output.check + format.bvolume output + format.edition output + new.sentence + format.number.series output + format.publisher.address output + format.date "year" output.check + date.block +% add.blank + format.pages "pages" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {incollection} +{ output.bibitem + format.authors "author" output.check + new.block + format.title "title" output.check + new.block + format.in.ed.booktitle "booktitle" output.check + format.bvolume output + format.edition output + new.sentence + format.number.series output + format.publisher.address output + format.date "year" output.check + date.block +% add.blank + format.pages "pages" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {inproceedings} +{ output.bibitem + format.authors "author" output.check + new.block + format.title "title" output.check + new.block + format.in.ed.booktitle "booktitle" output.check + format.bvolume output + new.sentence + format.number.series output + publisher empty$ + { format.organization.address output } + { organization "organization" bibinfo.check output + format.publisher.address output + } + if$ + format.date "year" output.check + date.block +% add.blank + format.pages "pages" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {conference} +{inproceedings} + +FUNCTION {manual} +{misc} + +FUNCTION {phdthesis} +{ output.bibitem + format.authors "author" output.check + new.block + format.btitle + "title" output.check + format.type "type" bibinfo.check output + new.block + school "school" bibinfo.warn output + address "address" bibinfo.check output + format.date "year" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {mastersthesis} +{phdthesis} + +FUNCTION {proceedings} +{ output.bibitem + editor empty$ + { organization "organization" bibinfo.check output + } + { format.editors output.nonnull } + if$ + new.block + format.btitle "title" output.check + format.bvolume output + editor empty$ + { publisher empty$ + 'skip$ + { + new.sentence + format.number.series output + format.publisher.address output + } + if$ + } + { publisher empty$ + { + new.sentence + format.organization.address output } + { + new.sentence + organization "organization" bibinfo.check output + format.publisher.address output + } + if$ + } + if$ + format.date "year" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {techreport} +{ output.bibitem + format.authors "author" output.check + new.block + format.title + "title" output.check + new.block + format.institution.address output + format.date "year" output.check + format.tr.number output.nonnull + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {map} +{ output.bibitem + format.cartographers "cartographer" output.check + new.block + format.title + "title" output.check + format.type "type" + new.block + format.publisher.address output + format.date "year" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {patent} +{ output.bibitem + format.inventors "inventor" output.check + "; " * + no.blank.or.punct + format.assignees "assignee" output.check + new.block + format.title + "title" output.check + new.block + format.tr.number output.nonnull + format.date "year" output.check + new.block +% format.note output + output.url % urlbst + fin.entry +} + +FUNCTION {unpublished} +{ output.bibitem + format.authors "author" output.check + new.block + format.title "title" output.check + format.date output + new.block +% format.note "note" output.check + output.url % urlbst + fin.entry +} + +FUNCTION {default.type} { misc } +READ +STRINGS { longest.label } +INTEGERS { number.label longest.label.width } +FUNCTION {initialize.longest.label} +{ "" 'longest.label := + #1 'number.label := + #0 'longest.label.width := +} +FUNCTION {longest.label.pass} +{ number.label int.to.str$ 'label := + number.label #1 + 'number.label := + label width$ longest.label.width > + { label 'longest.label := + label width$ 'longest.label.width := + } + 'skip$ + if$ +} +EXECUTE {initialize.longest.label} +ITERATE {longest.label.pass} +FUNCTION {begin.bib} +{ preamble$ empty$ + 'skip$ + { preamble$ write$ newline$ } + if$ + "\begin{thebibliography}{" longest.label * "}" * + write$ newline$ +} +EXECUTE {begin.bib} +EXECUTE {init.config.constants} +EXECUTE {init.state.consts} +ITERATE {call.type$} +FUNCTION {end.bib} +{ newline$ + "\end{thebibliography}" write$ newline$ +} +EXECUTE {end.bib} +%% End of customized bst file +%% +%% End of file `vancouver.bst'. diff --git a/susceptibility1.tex b/susceptibility1.tex index d084228..9bd6a95 100644 --- a/susceptibility1.tex +++ b/susceptibility1.tex @@ -1,86 +1,175 @@ -% -%% Copyright 2007-2020 Elsevier Ltd -%% -%% This file is part of the 'Elsarticle Bundle'. -%% --------------------------------------------- -%% -%% It may be distributed under the conditions of the LaTeX Project Public -%% License, either version 1.2 of this license or (at your option) any -%% later version. The latest version of this license is in -%% http://www.latex-project.org/lppl.txt -%% and version 1.2 or later is part of all distributions of LaTeX -%% version 1999/12/01 or later. -%% -%% The list of all files belonging to the 'Elsarticle Bundle' is -%% given in the file `manifest.txt'. -%% - -%% Template article for Elsevier's document class `elsarticle' -%% with numbered style bibliographic references -%% SP 2008/03/01 -%% -%% -%% -%% $Id: elsarticle-template-num.tex 190 2020-11-23 11:12:32Z rishi $ -%% -%% -%% Use the option review to obtain double line spacing -%% \documentclass[authoryear,preprint,review,12pt]{elsarticle} - -%% Use the options 1p,twocolumn; 3p; 3p,twocolumn; 5p; or 5p,twocolumn -%% for a journal layout: -%% \documentclass[final,1p,times]{elsarticle} -%% \documentclass[final,1p,times,twocolumn]{elsarticle} -%% \documentclass[final,3p,times]{elsarticle} -%% \documentclass[final,3p,times,twocolumn]{elsarticle} -%% \documentclass[final,5p,times]{elsarticle} -%% \documentclass[final,5p,times,twocolumn]{elsarticle} - -% \documentclass[preprint,12pt]{elsarticle} -%\documentclass[final,5p,times,twocolumn]{elsarticle} -%\documentclass[final, 10pt, 5p, twocolumn, authoryear]{elsarticle} - -\documentclass[preprint, 10pt, 3p]{elsarticle} -\usepackage{natbib} -%\documentclass[preprint, 10pt, 3p]{elsarticle} -%\usepackage{natbib} -\setcitestyle{super,comma,sort&compress} -%\setcitestyle{super,comma,sort&compress} - -%%%%% units %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\usepackage[mediumspace,mediumqspace, Gray,amssymb]{SIunits} % \ohm, \micro - -%%%%% layout %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%\usepackage[left=20mm,right=20mm,top=20mm,bottom=20mm]{geometry} -%\usepackage[left=26mm,right=14mm,top=20mm,bottom=20mm]{geometry} -%\usepackage[left=26mm,right=14mm,top=20mm,bottom=20mm]{geometry} -%\usepackage[hmargin=35mm, vmargin=25mm]{geometry} -%\setcounter{secnumdepth}{-1} - -%\renewcommand*{\chaptermarkformat}{\chapapp~\thechapter\\} - - - - -%%%%% graphics %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\usepackage{xcolor} -\usepackage{graphicx} -\usepackage{ifthen} -\usepackage{booktabs} +% Template for PLoS +% Version 3.6 Aug 2022 +% +% % % % % % % % % % % % % % % % % % % % % % +% +% -- IMPORTANT NOTE +% +% This template contains comments intended +% to minimize problems and delays during our production +% process. Please follow the template instructions +% whenever possible. +% +% % % % % % % % % % % % % % % % % % % % % % % +% +% Once your paper is accepted for publication, +% PLEASE REMOVE ALL TRACKED CHANGES in this file +% and leave only the final text of your manuscript. +% PLOS recommends the use of latexdiff to track changes during review, as this will help to maintain a clean tex file. +% Visit https://www.ctan.org/pkg/latexdiff?lang=en for info or contact us at latex@plos.org. +% +% +% There are no restrictions on package use within the LaTeX files except that no packages listed in the template may be deleted. +% +% Please do not include colors or graphics in the text. +% +% The manuscript LaTeX source should be contained within a single file (do not use \input, \externaldocument, or similar commands). +% +% % % % % % % % % % % % % % % % % % % % % % % +% +% -- FIGURES AND TABLES +% +% Please include tables/figure captions directly after the paragraph where they are first cited in the text. +% +% DO NOT INCLUDE GRAPHICS IN YOUR MANUSCRIPT +% - Figures should be uploaded separately from your manuscript file. +% - Figures generated using LaTeX should be extracted and removed from the PDF before submission. +% - Figures containing multiple panels/subfigures must be combined into one image file before submission. +% For figure citations, please use "Fig" instead of "Figure". +% See http://journals.plos.org/plosone/s/figures for PLOS figure guidelines. +% +% Tables should be cell-based and may not contain: +% - spacing/line breaks within cells to alter layout or alignment +% - do not nest tabular environments (no tabular environments within tabular environments) +% - no graphics or colored text (cell background color/shading OK) +% See http://journals.plos.org/plosone/s/tables for table guidelines. +% +% For tables that exceed the width of the text column, use the adjustwidth environment as illustrated in the example table in text below. +% +% % % % % % % % % % % % % % % % % % % % % % % % +% +% -- EQUATIONS, MATH SYMBOLS, SUBSCRIPTS, AND SUPERSCRIPTS +% +% IMPORTANT +% Below are a few tips to help format your equations and other special characters according to our specifications. For more tips to help reduce the possibility of formatting errors during conversion, please see our LaTeX guidelines at http://journals.plos.org/plosone/s/latex +% +% For inline equations, please be sure to include all portions of an equation in the math environment. For example, x$^2$ is incorrect; this should be formatted as $x^2$ (or $\mathrm{x}^2$ if the romanized font is desired). +% +% Do not include text that is not math in the math environment. For example, CO2 should be written as CO\textsubscript{2} instead of CO$_2$. +% +% Please add line breaks to long display equations when possible in order to fit size of the column. +% +% For inline equations, please do not include punctuation (commas, etc) within the math environment unless this is part of the equation. +% +% When adding superscript or subscripts outside of brackets/braces, please group using {}. For example, change "[U(D,E,\gamma)]^2" to "{[U(D,E,\gamma)]}^2". +% +% Do not use \cal for caligraphic font. Instead, use \mathcal{} +% +% % % % % % % % % % % % % % % % % % % % % % % % +% +% Please contact latex@plos.org with any questions. +% +% % % % % % % % % % % % % % % % % % % % % % % % +\documentclass[10pt,letterpaper]{article} +\usepackage[top=0.85in,left=2.75in,footskip=0.75in]{geometry} + +% amsmath and amssymb packages, useful for mathematical formulas and symbols +\usepackage{amsmath,amssymb} + +% Use adjustwidth environment to exceed column width (see example table in text) +\usepackage{changepage} + +% textcomp package and marvosym package for additional characters +\usepackage{textcomp,marvosym} + +% cite package, to clean up citations in the main text. Do not remove. +\usepackage{cite} + +% Use nameref to cite supporting information files (see Supporting Information section for more info) +\usepackage{nameref,hyperref} + +% line numbers +\usepackage[right]{lineno} + +% ligatures disabled +\usepackage[nopatch=eqnum]{microtype} +\DisableLigatures[f]{encoding = *, family = * } + +% color can be used to apply background shading to table cells only +\usepackage[table]{xcolor} + +% array package and thick rules for tables +\usepackage{array} + +% create "+" rule type for thick vertical lines +\newcolumntype{+}{!{\vrule width 2pt}} +% create \thickcline for thick horizontal lines of variable length +\newlength\savedwidth +\newcommand\thickcline[1]{% + \noalign{\global\savedwidth\arrayrulewidth\global\arrayrulewidth 2pt}% + \cline{#1}% + \noalign{\vskip\arrayrulewidth}% + \noalign{\global\arrayrulewidth\savedwidth}% +} +% \thickhline command for thick horizontal lines that span the table +\newcommand\thickhline{\noalign{\global\savedwidth\arrayrulewidth\global\arrayrulewidth 2pt}% +\hline +\noalign{\global\arrayrulewidth\savedwidth}} + +% Remove comment for double spacing +%\usepackage{setspace} +%\doublespacing + +% Text layout +\raggedright +\setlength{\parindent}{0.5cm} +\textwidth 5.25in +\textheight 8.75in + +% Bold the 'Figure #' in the caption and separate it from the title/caption with a period +% Captions will be left justified +\usepackage[aboveskip=1pt,labelfont=bf,labelsep=period,justification=raggedright,singlelinecheck=off]{caption} +\renewcommand{\figurename}{Fig} + +% Use the PLoS provided BiBTeX style +\bibliographystyle{plos2015} + +% Remove brackets from numbering in List of References +\makeatletter +\renewcommand{\@biblabel}[1]{\quad#1.} +\makeatother + +% Header and Footer with logo +\usepackage{lastpage,fancyhdr,graphicx} +\usepackage{epstopdf} +%\pagestyle{myheadings} +\pagestyle{fancy} +\fancyhf{} +%\setlength{\headheight}{27.023pt} +%\lhead{\includegraphics[width=2.0in]{PLOS-submission.eps}} +\rfoot{\thepage/\pageref{LastPage}} +\renewcommand{\headrulewidth}{0pt} +\renewcommand{\footrule}{\hrule height 2pt \vspace{2mm}} +\fancyheadoffset[L]{2.25in} +\fancyfootoffset[L]{2.25in} +\lfoot{\today} + +\usepackage{lipsum} +\usepackage[mediumspace,mediumqspace,Gray,amssymb]{SIunits} +\usepackage{ifthen} +%% BEGIN MACROS SECTION -%%%%% math %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\usepackage{textcomp} -\usepackage{array} -\usepackage{amsmath} -\usepackage{amssymb} -\usepackage{graphicx} -\usepackage{float} +%%%%% species names %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +\newcommand{\Lepto}{\emph{Apteronotus leptorhynchus}} +\newcommand{\lepto}{\emph{A. leptorhynchus}} +\newcommand{\Eigen}{\emph{Eigenmannia virescens}} +\newcommand{\eigen}{\emph{E. virescens}} %%%%% notes %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newcommand{\note}[2][]{\textbf{[#1: #2]}} -\newcommand{\notesr}[1]{\note[SR]{#1}} +\newcommand{\noteab}[1]{\note[AB]{#1}} \newcommand{\notejb}[1]{\note[JB]{#1}} \newcommand{\notejg}[1]{\note[JG]{#1}} \newcommand{\notems}[1]{\note[MS]{#1}} @@ -89,21 +178,12 @@ %%%%% equation references %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \renewcommand{\eqref}[1]{(\ref{#1})} \newcommand{\eref}[1]{\ref{#1}} -\newcommand{\eqn}{Eq.} -\newcommand{\Eqn}{Eq.} -\newcommand{\eqns}{Eqs.} -\newcommand{\Eqns}{Eqs.} -\newcommand{\eqnref}[1]{\eqn~\eqref{#1}} -\newcommand{\enref}[1]{\eqn~\eref{#1}} +\newcommand{\eqn}{Eq} +\newcommand{\Eqn}{Eq} +\newcommand{\eqns}{Eqs} +\newcommand{\Eqns}{Eqs} \newcommand{\Eqnref}[1]{\Eqn~\eqref{#1}} -\newcommand{\eqnsref}[1]{\eqns~\eqref{#1}} \newcommand{\Eqnsref}[1]{\Eqns~\eqref{#1}} - -%%%%% hyperef %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\usepackage{xcolor} -\usepackage[breaklinks=true,colorlinks=true,citecolor=blue!30!black,urlcolor=blue!30!black,linkcolor=blue!30!black]{hyperref} -%\usepackage[breaklinks=true,colorlinks=true,citecolor=blue!30!black,urlcolor=blue!30!black,linkcolor=blue!30!black]{hyperref} - %%%%% figures %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % references to panels of a figure within the caption: @@ -115,9 +195,9 @@ \newcommand{\fref}[1]{\textup{\ref{#1}}} \newcommand{\subfref}[2]{\textup{\ref{#1}}\,\panel{#2}} % references to figures in normal text: -\newcommand{\fig}{Fig.} +\newcommand{\fig}{Fig} \newcommand{\Fig}{Figure} -\newcommand{\figs}{Figs.} +\newcommand{\figs}{Figs} \newcommand{\Figs}{Figures} \newcommand{\figref}[1]{\fig~\fref{#1}} \newcommand{\Figref}[1]{\Fig~\fref{#1}} @@ -128,13 +208,41 @@ \newcommand{\subfigsref}[2]{\figs~\subfref{#1}{#2}} \newcommand{\Subfigsref}[2]{\Figs~\subfref{#1}{#2}} % references to figures within brackets: -\newcommand{\figb}{Fig.} -\newcommand{\figsb}{Figs.} +\newcommand{\figb}{Fig} +\newcommand{\figsb}{Figs} \newcommand{\figrefb}[1]{\figb~\fref{#1}} \newcommand{\figsrefb}[1]{\figsb~\fref{#1}} \newcommand{\subfigrefb}[2]{\figb~\subfref{#1}{#2}} \newcommand{\subfigsrefb}[2]{\figsb~\subfref{#1}{#2}} +%%%%% changes %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +\newcommand{\remove}[2][]{% + \protect\ifthenelse{\value{showtrackchanges}<2}% + {\protect\ifthenelse{\value{showtrackchanges}=0}% + {}% 0 + {\ifthenelse{\equal{#1}{}}{}{\textit{\textbf{[#1:]}} }\st{#2}}}% 1 + {\protect\ifthenelse{\value{showtrackchanges}=2}% + {\protect\marginpar{\protect\rule[-0.3ex]{0.5mm}{2ex}}\ifthenelse{\equal{#1}{}}{}{\textit{\textbf{[#1:]}} }\st{#2}}% 2 + {\protect\marginpar{\small\textbf{deleted\ifthenelse{\equal{#1}{}}{}{ [#1]}:} #2}}}}% 3 + +\newcommand{\add}[2][]{% + \protect\ifthenelse{\value{showtrackchanges}<2}% + {\protect\ifthenelse{\value{showtrackchanges}=0}% + {#2}% 0 + {\textcolor{red}{\ifthenelse{\equal{#1}{}}{}{\textit{[#1:]} }#2}}}% 1 + {\protect\ifthenelse{\value{showtrackchanges}=2}% + {{\protect\marginpar{\protect\rule[-0.3ex]{0.5mm}{2ex}}\bfseries \ifthenelse{\equal{#1}{}}{}{\textit{[#1:]} }#2}}% 2 + {\protect\marginpar{\small\textbf{inserted\ifthenelse{\equal{#1}{}}{}{ [#1]}}}{\bfseries #2}}}}% 3 + +\newcommand{\change}[2][]{% + \protect\ifthenelse{\value{showtrackchanges}<2}% + {\protect\ifthenelse{\value{showtrackchanges}=0}% + {#2}% 0 + {\textcolor{red}{\ifthenelse{\equal{#1}{}}{}{\textit{[#1:]} }#2}}}% 1 + {\protect\ifthenelse{\value{showtrackchanges}=2}% + {\protect\marginpar{\protect\rule[-0.3ex]{0.5mm}{2ex}}{\bfseries \ifthenelse{\equal{#1}{}}{}{\textit{[#1:]} }#2}}% 2 + {\protect\marginpar{\small\textbf{changed\ifthenelse{\equal{#1}{}}{}{ [#1]}}}{\bfseries #2}}}}% 3 + %%%%% new commands %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newcommand{\fstim}{\ensuremath{f_{stim}}} \newcommand{\feod}{\ensuremath{f_{EOD}}} @@ -146,9 +254,6 @@ \newcommand{\feodintro}{\ensuremath{\rm{EOD}_{1}}} \newcommand{\ffstimintro}{\ensuremath{f_{2}}} \newcommand{\ffeodintro}{\ensuremath{f_{1}}} -%intro - -%\newcommand{\signalcomp}{$s$} \newcommand{\baseval}{134} \newcommand{\bone}{$\Delta f_{1}$} @@ -158,7 +263,6 @@ \newcommand{\ff}{$f_{1}$--$f_{2}$} \newcommand{\auc}{\rm{AUC}} - \newcommand{\n}{\ensuremath{N}} \newcommand{\bvary}{\ensuremath{\Delta f_{1}}} \newcommand{\avary}{\ensuremath{A(\Delta f_{1})}} @@ -179,14 +283,6 @@ \newcommand{\ffall}{$f_{EOD}$\&$f_{1}$\&$f_{2}$} \newcommand{\ffvary}{$f_{EOD}$\&$f_{1}$}%sum \newcommand{\ffstable}{$f_{EOD}$\&$f_{2}$}%sum -%\newcommand{\colstableone}{dark blue}%sum -%\newcommand{\colstabletwo}{blue}%sum -%\newcommand{\colvaryone}{dark red}%sum -%\newcommand{\colvarytwo}{red}%sum -%\newcommand{\colstableoneb}{Dark blue}%sum -%\newcommand{\colstabletwob}{Blue}%sum -%\newcommand{\colvaryoneb}{Dark red}%sum -%\newcommand{\colvarytwob}{Red}%sum \newcommand{\colstableone}{blue}%sum \newcommand{\colstabletwo}{cyan}%sum \newcommand{\colvaryone}{brown}%sum @@ -201,7 +297,6 @@ \newcommand{\rocf}{98}%sum %\Delta \newcommand{\roci}{36}%sum %\Delta - %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Beat combinations \newcommand{\boneabs}{\ensuremath{|\Delta f_{1}|}}%sum @@ -222,8 +317,25 @@ \newcommand{\bdiffehalf}{$\bdiff{}=f_{EOD}/2$}%diff of both %beat frequencies -%%%%%%%%%%%%%%%%%%%%%% -% Frequency combinations +%%%%% tables %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% references to tables: +\newcommand{\tref}[1]{\textup{\ref{#1}}} +% references to tables in normal text: +\newcommand{\tab}{Tab.} +\newcommand{\Tab}{Table} +\newcommand{\tabs}{Tabs.} +\newcommand{\Tabs}{Tables} +\newcommand{\tabref}[1]{\tab~\tref{#1}} +\newcommand{\Tabref}[1]{\Tab~\tref{#1}} +\newcommand{\tabsref}[1]{\tabs~\tref{#1}} +\newcommand{\Tabsref}[1]{\Tabs~\tref{#1}} +% references to tables within the bracketed text: +\newcommand{\tabb}{Tab.} +\newcommand{\tabsb}{Tab.} +\newcommand{\tabrefb}[1]{\tabb~\tref{#1}} +\newcommand{\tabsrefb}[1]{\tabsb~\tref{#1}} + +%%%%%%%%%%%%%%%%%%%%%%% Frequency combinations \newcommand{\fsum}{\ensuremath{f_{1} + f_{2}}}%sum \newcommand{\fdiff}{\ensuremath{|f_{1}-f_{2}|}}%diff of @@ -241,7 +353,6 @@ \newcommand{\fdiffe}{$\fdiff{}=f_{EOD}$}%diff of both \newcommand{\fdiffehalf}{$\fdiff{}=f_{EOD}/2$}%diff of both - \newcommand{\fctwo}{\ensuremath{f_{\rm{Female}}}}%sum \newcommand{\fcone}{\ensuremath{f_{\rm{Intruder}}}}%sum \newcommand{\bctwo}{\ensuremath{\Delta \fctwo{}}}%sum @@ -275,24 +386,13 @@ \newcommand{\bcdiffe}{\ensuremath{\bcdif{}f=f_{\rm{EOD}}}}%diff of both \newcommand{\bcdiffehalf}{\ensuremath{\bcdiff{}=f_{\rm{EOD}}/2}}%diff of both - -%\newcommand{\spnr}{$\langle spikes_{burst}\rangle$}%diff of \newcommand{\burstcorr}{\ensuremath{{Corrected}}} \newcommand{\cvbasecorr}{CV\ensuremath{_{BaseCorrected}}} \newcommand{\cv}{CV\ensuremath{_{Base}}}%\cvbasecorr{} \newcommand{\nli}{NLI\ensuremath{(\fbase{})}}%Fr$_{Burst}$ \newcommand{\nlicorr}{NLI\ensuremath{(\fbasecorr{}})}%Fr$_{Burst}$ -\newcommand{\rostra}{\textit{Apteronotus rostratus}} -\newcommand{\lepto}{\textit{Apteronotus leptorhynchus}} -\newcommand{\eigen}{\textit{Eigenmannia virescens}} \newcommand{\suscept}{$|\chi_{2}|$} \newcommand{\frcolor}{pink lines} -%\newcommand{\rec}{\ensuremath{f_{Receiver}}} -%\newcommand{\intruder}{\ensuremath{f_{Intruder}}} -%\newcommand{\female}{\ensuremath{f_{Female}}} -%\newcommand{\rif}{\ensuremath{\rec{} & \intruder{} & \female{}}} -%\newcommand{\ri}{\ensuremath{\rec{} & \intruder{}}}%sum -%\newcommand{\rf}{\ensuremath{\rec{} & \female{}}}%sum \newcommand{\rec}{\ensuremath{\rm{R}}}%{\ensuremath{con_{R}}} \newcommand{\rif}{\ensuremath{\rm{RIF}}} @@ -302,85 +402,75 @@ \newcommand{\wofemale}{\rm{ROC}\ensuremath{\rm{_{NoFemale}}}}%sum \newcommand{\dwithfemale}{\rm{\ensuremath{\auc_{Female}}}}%sum CV\ensuremath{_{BaseCorrected}} \newcommand{\dwofemale}{\rm{\ensuremath{\auc_{NoFemale}}}}%sum -%\newcommand{\withfemale}{\rm{ROC}\ensuremath{_{Female}}}%sum \textit{ -%\newcommand{\wofemale}{\rm{ROC}\ensuremath{_{NoFemale}}}%sum -%\newcommand{\dwithfemale}{\auc\ensuremath{_{Female}}}%sum CV\ensuremath{_{BaseCorrected}} -%\newcommand{\dwofemale}{\auc\ensuremath{_{NoFemale}}}%sum \newcommand{\fp}{\ensuremath{\rm{FP}}}%sum \newcommand{\cd}{\ensuremath{\rm{CD}}}%sum -%\dwithfemale - -%%%%% tables %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% references to tables: -\newcommand{\tref}[1]{\textup{\ref{#1}}} -% references to tables in normal text: -\newcommand{\tab}{Tab.} -\newcommand{\Tab}{Table} -\newcommand{\tabs}{Tabs.} -\newcommand{\Tabs}{Tables} -\newcommand{\tabref}[1]{\tab~\tref{#1}} -\newcommand{\Tabref}[1]{\Tab~\tref{#1}} -\newcommand{\tabsref}[1]{\tabs~\tref{#1}} -\newcommand{\Tabsref}[1]{\Tabs~\tref{#1}} -% references to tables within the bracketed text: -\newcommand{\tabb}{Tab.} -\newcommand{\tabsb}{Tab.} -\newcommand{\tabrefb}[1]{\tabb~\tref{#1}} -\newcommand{\tabsrefb}[1]{\tabsb~\tref{#1}} - -\usepackage{xr} - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%% notes -%% The line no packages add line numbers. Start line numbering with -%% \begin{linenumbers}, end it with \end{linenumbers}. Or switch it on -%% for the whole article with \linenumbers. -\usepackage{lineno} -%\linenumbers -%\setlength{\linenumbersep}{3pt} - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% caption packages -%\usepackage{graphicx} -\usepackage[format=plain,labelfont=bf]{caption} -%\usepackage{capt-of} -\captionsetup[figure2]{labelformat=empty} -%\usepackage[onehalfspacing]{setspace} - -\journal{iScience} -\setcounter{secnumdepth}{0} +%% END MACROS SECTION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{document} - -%\maketitle -\begin{frontmatter} - -\title{Second-order susceptibility in electrosensory primary afferents in a three-fish setting} - - -\author[1]{Alexandra Barayeu} %\corref{fnd1}} -% \ead{alexandra.rudnaya@uni-tuebingen.de} -\author[4,5]{Maria Schlungbaum} -\author[4,5]{Benjamin Lindner} -\author[1,2,3]{Jan Benda} -\author[1]{Jan Grewe\corref{cor1}} -\ead{jan.grewe@uni-tuebingen.de} -% \ead[url]{home page} - -\cortext[cor1]{Corresponding author} -\affiliation[1]{organization={Neuroethology, Institute for Neurobiology, Eberhard Karls University}, - city={T\"ubingen}, postcode={72076}, country={Germany}} -\affiliation[4]{organization={Bernstein Center for Computational Neuroscience T\"ubingen}, - city={T\"ubingen}, postcode={72076}, country={Germany}} -\affiliation[5]{organization={Werner Reichardt Centre for Integrative Neuroscience}, - city={T\"ubingen}, postcode={72076}, country={Germany}} -\affiliation[2]{organization={Bernstein Center for Computational Neuroscience Berlin}, city={Berlin}, country={Germany}} -\affiliation[3]{organization={Physics Department of Humboldt University Berlin},city={Berlin}, country={Germany}} +\vspace*{0.2in} +% Title must be 250 characters or less. +\begin{flushleft} + {\Large + \textbf\newline{Second-order susceptibility in electrosensory primary afferents in a three-fish setting} + } + \newline + % Insert author names, affiliations and corresponding author email (do not include titles, positions, or degrees). + \\ + Alexandra Barayeu\textsuperscript{1}, + Maria Schlungbaum\textsuperscript{2,3}, + Benjamin Lindner\textsubscript{2,3}, + Jan Benda\textsuperscript{1, 4} + Jan Grewe\textsuperscript{1, *} + \\ + \bigskip + \textbf{1} Institute for Neurobiology, Eberhardt Karls Universit\"at T\"ubingen, 72076 T\"ubingen, Germany\\ + \textbf{2} Bernstein Center for Computational Neuroscience Berlin, Berlin, Germany\\ + \textbf{3} Department of Physics, Humboldt University Berlin, Berlin, Germany\\ + \textbf{4} Bernstein Center for Computational Neuroscience Tübingen, Tübingen, 72076, Germany\\ + \bigskip + + % Insert additional author notes using the symbols described below. Insert symbol callouts after author names as necessary. + % + % Remove or comment out the author notes below if they aren't used. + % + \paragraph{Short title:} + \paragraph{Corresponding author:}Jan Grewe, E-mail: jan.grewe@uni-tuebingen.de + \paragraph{Conflict of interest:}The authors declare no conflict of interest. + \paragraph{Author Contributions:} All authors designed the study and discussed the results. AB performed the data analyses and modelling, AB and JG drafted the paper, all authors discussed and revised the manuscript. + + \paragraph{Keywords:} population coding $|$ conduction delay $|$ heterogeneity $|$ electric fish $|$ mutual information + + % Additional Equal Contribution Note + % Also use this double-dagger symbol for special authorship notes, such as senior authorship. + % \ddag These authors also contributed equally to this work. + + % Current address notes + % \textcurrency Current Address: Dept/Program/Center, Institution Name, City, State, Country % change symbol to "\textcurrency a" if more than one current address note + % \textcurrency b Insert second current address + % \textcurrency c Insert third current address + + % Deceased author note + % \dag Deceased + + % Group/Consortium Author Note + % \textpilcrow Membership list can be found in the Acknowledgments section. + + % Use the asterisk to denote corresponding authorship and provide email address in note below. + * jan.grewe@uni-tuebingen.de + + \end{flushleft} + % Please keep the abstract below 300 words + \section*{Abstract} + \lipsum[1-1] + % Please keep the Author Summary between 150 and 200 words + % Use first person. PLOS ONE authors please skip this step. + % Author Summary not valid for PLOS ONE submissions. + \section*{Author summary} + \lipsum[1-1] + \linenumbers %Nonlinearities contribute to the encoding of the full behaviorally relevant signal range in primary electrosensory afferents. @@ -395,105 +485,98 @@ %Nonlinearities facilitate the encoding of a wide range of beat frequencies and amplitudes in primary electrosensory afferents - -\begin{abstract} -In this work, the influence of nonlinearities on stimulus encoding in the primary sensory afferents of weakly electric fish of the species \lepto{} was investigated. These fish produce an electric organ discharge (EOD) with a fish-specific frequency. When the EOD of one fish interferes with the EOD of another fish, it results in a signal with a periodic amplitude modulation, called beat. The beat provides information about the sex and size of the encountered conspecific and is the basis for communication. The beat frequency is predicted as the difference between the EOD frequencies and the beat amplitude corresponds to the size of the smaller EOD field. Primary sensory afferents, the P-units, phase-lock to the EOD and encode beats with changes in their firing rate. In this work, the nonlinearities of primary electrosensory afferents, the P-units of weakly electric fish of the species \lepto{} and \eigen{} were addressed. Nonlinearities were characterized as the second-order susceptibility of P-units, in a setting where at least three fish were present. The nonlinear responses of P-units were especially strong in regular firing P-units. White noise stimulation was confirmed as a method to retrieve the socond-order suscepitbility in P-units.% with bursting being identified as a factor enhancing nonlinear interactions. -\end{abstract} - -\end{frontmatter} - -\section{Introduction} +\section*{Introduction} %with nonlinearities being observed in all sensory modalities -Nonlinear processes are fundamental in neuronal information processing. On the systemic level: deciding to take one or another action is a nonlinear process. On a finer scale, neurons are inherently nonlinear: whether an action potential is elicited depends on the membrane potential to exceed a certain threshold\citealp{Adelson1985, Brincat2004, Chacron2000, Chacron2001, Nelson1997, Gussin2007, Middleton2007, Longtin2008}. In conjunction with neuronal noise, nonlinear mechanism facilitate the encoding of weak stimuli via stochastic resonance\citealp{Wiesenfeld1995, Stocks2000,Neiman2011fish}. We can find nonlinearities in many sensory systems such as rectification in the transduction machinery of inner hair cells \citealp{Peterson2019}, signal rectification in electroreceptor cells \citealp{Chacron2000, Chacron2001}, or in complex cells of the visual system \citealp{Adelson1985}. In the auditory or the active electric sense, for example, nonlinear processes are needed to extract envelopes, i.e. amplitude modulations of a carrier signal\citealp{Joris2004, Barayeu2023} called beats. Beats are common stimuli in different sensory modalities enabling rhythm and pitch perception in human hearing \citealp{Roeber1834, Plomp1967, Joris2004, Grahn2012} and providing context for electrocommunication in weakly electric fish \citealp{Engler2001, Hupe2008, Henninger2018, Benda2020}. +Nonlinear processes are fundamental in neuronal information processing. On the systemic level: deciding to take one or another action is a nonlinear process. On a finer scale, neurons are inherently nonlinear: whether an action potential is elicited depends on the membrane potential to exceed a certain threshold\cite{Adelson1985, Brincat2004, Chacron2000, Chacron2001, Nelson1997, Gussin2007, Middleton2007, Longtin2008}. In conjunction with neuronal noise, nonlinear mechanism facilitate the encoding of weak stimuli via stochastic resonance\cite{Wiesenfeld1995, Stocks2000,Neiman2011fish}. We can find nonlinearities in many sensory systems such as rectification in the transduction machinery of inner hair cells \cite{Peterson2019}, signal rectification in electroreceptor cells \cite{Chacron2000, Chacron2001}, or in complex cells of the visual system \cite{Adelson1985}. In the auditory or the active electric sense, for example, nonlinear processes are needed to extract envelopes, i.e. amplitude modulations of a carrier signal\cite{Joris2004, Barayeu2023} called beats. Beats are common stimuli in different sensory modalities enabling rhythm and pitch perception in human hearing \cite{Roeber1834, Plomp1967, Joris2004, Grahn2012} and providing context for electrocommunication in weakly electric fish \cite{Engler2001, Hupe2008, Henninger2018, Benda2020}. -While the sensory periphery can often be well described by linear models, this is not true for many upstream neurons. Rather, nonlinear processes are implemented to extract special stimulus features\citealp{Gabbiani1996}. In active electrosensation, the self-generated electric field (electric organ discharge, EOD)\citealp{Salazar2013} that is quasi sinusoidal in wavetype electric fish and acts as the carrier signal that is amplitude modulated in the context of communication\citealp{Benda2013, Fotowat2013, Walz2014, Henninger2018} as well as object detection and navigation\citealp{Fotowat2013, Nelson1999}. In social contexts, the interference of the EODs of two interacting animals result in a characteristic periodic amplitude modulation, the so-called beat. The beat amplitude is defined by the smaller EOD amplitude, its frequency is defined as the difference between the two EOD frequencies ($\Delta f = f-\feod{}$, valid for $f < feod{}/2$)\citealp{Barayeu2023}. Cutaneous electroreceptor organs that are distributed over the bodies of these fish \citealp{Carr1982} are tuned to the own field\citealp{Hopkins1976,Viancour1979}. P-type electroreceptor afferents (P-units) innervate these organs via ribbon synapses\citealp{Szabo1965, Wachtel1966} and project to the hindbrain where they trifurcate and synapse onto pyramidal cells in the electrosensory lateral line lobe (ELL)\citealp{Krahe2014}. The P-units ot the gymnotiform electric fish \lepto{} encode such amplitude modulations (AMs) by modulation of their firing rate\citealp{Gabbiani1996}. They fire probabilistically but phase-locked to the own EOD and the skipping of cycles leads to their characteristic multimodal interspike interval distribution. Even though the extraction of the AM itself requires a nonlinearity\citealp{Middleton2006,Stamper2012Envelope,Savard2011} encoding the time-course of the AM is linear over a wide range\citealp{Xu1996,Benda2005,Gussin2007,Grewe2017,Savard2011}. In the context of social signalling among three fish we observe an AM of the AM, also referred to as second-order envelope or just social envelope\citealp{Middleton2006, Savard2011, Stamper2012Envelope}. Encoding this again requires nonlinearities\citealp{Middleton2006} and it was shown that a subpopulation of P-units are sensitive to envelopes\citealp{Savard2011} and exhibit nonlinearities e.g. when driven by strong stimuli\citealp{Nelson1997,Chacron2004}. +While the sensory periphery can often be well described by linear models, this is not true for many upstream neurons. Rather, nonlinear processes are implemented to extract special stimulus features\cite{Gabbiani1996}. In active electrosensation, the self-generated electric field (electric organ discharge, EOD)\cite{Salazar2013} that is quasi sinusoidal in wavetype electric fish and acts as the carrier signal that is amplitude modulated in the context of communication\cite{Benda2013, Fotowat2013, Walz2014, Henninger2018} as well as object detection and navigation\cite{Fotowat2013, Nelson1999}. In social contexts, the interference of the EODs of two interacting animals result in a characteristic periodic amplitude modulation, the so-called beat. The beat amplitude is defined by the smaller EOD amplitude, its frequency is defined as the difference between the two EOD frequencies ($\Delta f = f-\feod{}$, valid for $f < feod{}/2$)\cite{Barayeu2023}. Cutaneous electroreceptor organs that are distributed over the bodies of these fish \cite{Carr1982} are tuned to the own field\cite{Hopkins1976,Viancour1979}. P-type electroreceptor afferents (P-units) innervate these organs via ribbon synapses\cite{Szabo1965, Wachtel1966} and project to the hindbrain where they trifurcate and synapse onto pyramidal cells in the electrosensory lateral line lobe (ELL)\cite{Krahe2014}. The P-units ot the gymnotiform electric fish \lepto{} encode such amplitude modulations (AMs) by modulation of their firing rate\cite{Gabbiani1996}. They fire probabilistically but phase-locked to the own EOD and the skipping of cycles leads to their characteristic multimodal interspike interval distribution. Even though the extraction of the AM itself requires a nonlinearity\cite{Middleton2006,Stamper2012Envelope,Savard2011} encoding the time-course of the AM is linear over a wide range\cite{Xu1996,Benda2005,Gussin2007,Grewe2017,Savard2011}. In the context of social signalling among three fish we observe an AM of the AM, also referred to as second-order envelope or just social envelope\cite{Middleton2006, Savard2011, Stamper2012Envelope}. Encoding this again requires nonlinearities\cite{Middleton2006} and it was shown that a subpopulation of P-units are sensitive to envelopes\cite{Savard2011} and exhibit nonlinearities e.g. when driven by strong stimuli\cite{Nelson1997,Chacron2004}. -\begin{figure*}[h!] +\begin{figure*}[!ht] \includegraphics[width=\columnwidth]{motivation} - \caption{\label{motivation} Nonlinearity in an electrophysiologically recorded P-unit of \lepto{} in a three-fish setting. Receiver with EOD frequency $\feod{} =664$\,Hz encounters fish with EOD frequencies $f_{1}=631$\,Hz and $f_{2}=797$\,Hz. Both encountered fish lead to a beat contrast of 10\,\%. Top: Interference of the receiver EOD with the EODs of other fish. Second row: Spike trains of the P-unit. Third row: Firing rate, retrieved as the convolution of the spike trains with a Gaussian kernel ($\sigma = 1$\,ms). Bottom row: Power spectrum of the firing rate. \figitem{A} Baseline condition: Only the receiver is present. The baseline firing rate \fbase{} dominates the power spectrum of the firing rate. \figitem{B} The receiver and the fish with EOD frequency $f_{1}=631$\,Hz are present. \figitem{C} The receiver and the fish with EOD frequency $f_{2}=797$\,Hz are present. \figitem{D} All three fish with the EOD frequencies \feod{}, $f_{1}$ and $f_{2}$ are present. Nonlinear peaks occur at the sum and difference of the two beat frequencies in the power spectrum of the firing rate. + \caption{\label{fig:motivation} Nonlinearity in an electrophysiologically recorded P-unit of \lepto{} in a three-fish setting. Receiver with EOD frequency $\feod{} =664$\,Hz encounters fish with EOD frequencies $f_{1}=631$\,Hz and $f_{2}=797$\,Hz. Both encountered fish lead to a beat contrast of 10\,\%. Top: Interference of the receiver EOD with the EODs of other fish. Second row: Spike trains of the P-unit. Third row: Firing rate, retrieved as the convolution of the spike trains with a Gaussian kernel ($\sigma = 1$\,ms). Bottom row: Power spectrum of the firing rate. \figitem{A} Baseline condition: Only the receiver is present. The baseline firing rate \fbase{} dominates the power spectrum of the firing rate. \figitem{B} The receiver and the fish with EOD frequency $f_{1}=631$\,Hz are present. \figitem{C} The receiver and the fish with EOD frequency $f_{2}=797$\,Hz are present. \figitem{D} All three fish with the EOD frequencies \feod{}, $f_{1}$ and $f_{2}$ are present. Nonlinear peaks occur at the sum and difference of the two beat frequencies in the power spectrum of the firing rate. } \end{figure*} -Recent theoretical and modelling work on leaky integrate-and-fire (LIF) model neurons revealed specific pattern of nonlinear interactions when driven with weak sinewaves\citealp{Voronenko2017}. This situation is reminiscent of the real fish's situation in which three electric fish interact. The active electrosensory system is then exposed to the two AMs arising from the interactions with either foreign fish's field. Previous recordings in the natural habitat showed interactions of three animals in which two animals interacted closely (strong signals) and were interrupted by an intruding animal (weak signal). The intruder was detected at spatial distances that led to an extremely faint intruder signal while the strong signal of the other animal is present\citealp{Henninger2018}. When the receiver fish with EOD frequency \feod{} is alone, the P-unit fire action potentials at a spontaneous baseline rate. Accordingly, a peak at \fbase{} is present in the power spectrum of the neuronal response (\subfigrefb{motivation}{A}, bottom).\notejg{more baseline trials, fake it, if needed}. Phase-locking to the own field also leads to a representation of \feod{} in their firing rate (see \figref{cells_suscept}\panel{B})\citealp{Sinz2020}. The beat frequency when two fish interact is also represented in the respective responses and accordingly, there is a peak in power spectrum at $\Delta f_1$ or $\Delta f_2$ (\subfigrefb{motivation}{B, C}, respectively). When three fish encounter, all their waveforms interfere with both beat frequencies being present in the superimposed signal (\subfigrefb{motivation}{D}, top). The power spectrum of the neuronal activity contains both previously seen peaks, but also nonlinear interaction peaks at the sum frequency \bsum{} and the difference frequency \bdiff{} (\subfigrefb{motivation}{D}, bottom). The latter is known as the social envelope\citealp{Stamper2012Envelope, Savard2011}. The neuron shown here, clearly encodes the envelope. Whether P-units in general encode envelopes has been subject of controversy, some works do not consider P-units as envelope encoders\citealp{Middleton2006}, while others identify some P-units as envelope encoders\citealp{Savard2011}. +Recent theoretical and modelling work on leaky integrate-and-fire (LIF) model neurons revealed specific pattern of nonlinear interactions when driven with weak sinewaves\cite{Voronenko2017}. This situation is reminiscent of the real fish's situation in which three electric fish interact. The active electrosensory system is then exposed to the two AMs arising from the interactions with either foreign fish's field. Previous recordings in the natural habitat showed interactions of three animals in which two animals interacted closely (strong signals) and were interrupted by an intruding animal (weak signal). The intruder was detected at spatial distances that led to an extremely faint intruder signal while the strong signal of the other animal is present\cite{Henninger2018}. When the receiver fish with EOD frequency \feod{} is alone, the P-unit fire action potentials at a spontaneous baseline rate. Accordingly, a peak at \fbase{} is present in the power spectrum of the neuronal response (\subfigrefb{fig:motivation}{A}, bottom).\notejg{more baseline trials, fake it, if needed}. Phase-locking to the own field also leads to a representation of \feod{} in their firing rate (see \figref{fig:cells_suscept}\panel{B})\cite{Sinz2020}. The beat frequency when two fish interact is also represented in the respective responses and accordingly, there is a peak in power spectrum at $\Delta f_1$ or $\Delta f_2$ (\subfigrefb{fig:motivation}{B, C}, respectively). When three fish encounter, all their waveforms interfere with both beat frequencies being present in the superimposed signal (\subfigrefb{fig:motivation}{D}, top). The power spectrum of the neuronal activity contains both previously seen peaks, but also nonlinear interaction peaks at the sum frequency \bsum{} and the difference frequency \bdiff{} (\subfigrefb{fig:motivation}{D}, bottom). The latter is known as the social envelope\cite{Stamper2012Envelope, Savard2011}. The neuron shown here, clearly encodes the envelope. Whether P-units in general encode envelopes has been subject of controversy, some works do not consider P-units as envelope encoders\cite{Middleton2006}, while others identify some P-units as envelope encoders\cite{Savard2011}. -The P-unit responses can be partially explained by simple linear filters. The linear relation of the cellular response and the stimulus can be captured by the first-order susceptibility (or transfer function). As in Volterra series, higher-order terms describe the nonlinear interactions. We quantify the nonlinearity of P-unit encoding by estimating the second-order susceptibility from white-noise responses\citealp{Voronenko2017,Neiman2011fish,Nikias1993}. We address whether the theory still holds for neurons that do not encode the foreign signals directly but respond to the AM of a carrier. For this study we re-use a large set of recordings of afferents of both the active (P-units) and the passive electrosensory (ampullary cells) to work out which cellular features determine which cells show nonlinear encoding. +The P-unit responses can be partially explained by simple linear filters. The linear relation of the cellular response and the stimulus can be captured by the first-order susceptibility (or transfer function). As in Volterra series, higher-order terms describe the nonlinear interactions. We quantify the nonlinearity of P-unit encoding by estimating the second-order susceptibility from white-noise responses\cite{Voronenko2017,Neiman2011fish,Nikias1993}. We address whether the theory still holds for neurons that do not encode the foreign signals directly but respond to the AM of a carrier. For this study we re-use a large set of recordings of afferents of both the active (P-units) and the passive electrosensory (ampullary cells) to work out which cellular features determine which cells show nonlinear encoding. -\section{Results} -Theoretical work shows that leaky-integrate-and-fire (LIF) model neurons show a distinct pattern of nonlinear stimulus encoding when the model is driven by two cosine signals. In the context of the weakly electric fish, such a setting is part of the animal's everyday life as the sinusoidal electric-organ discharges (EODs) of neighboring animals interfere with the own field and each lead to sinusoidal amplitude modulations (AMs) that are called beats and envelopes \citealp{Middleton2006, Savard2011,Stamper2012Envelope}. The p-type electroreceptor afferents of the tuberous electrosensory system, i.e. the P-units, respond to such AMs of the underlying EOD carrier and their time-dependent firing rate carries information about the stimulus' time-course. P-units are heterogeneous in their baseline firing properties \citealp{Grewe2017, Hladnik2023} and differ with respect to their baseline firing rates, the sensitivity and their noisiness. We here explore the nonlinear mechanism in different cells that exhibit distinctly different levels of noise, i.e. differences in the coefficient of variation (CV) of the interspike intervals (ISI). Low-CV P-units have a less noisy, more regular, firing pattern whereas high-CV P-units show a less regular firing pattern in their baseline activity. +\section*{Results} +Theoretical work shows that leaky-integrate-and-fire (LIF) model neurons show a distinct pattern of nonlinear stimulus encoding when the model is driven by two cosine signals. In the context of the weakly electric fish, such a setting is part of the animal's everyday life as the sinusoidal electric-organ discharges (EODs) of neighboring animals interfere with the own field and each lead to sinusoidal amplitude modulations (AMs) that are called beats and envelopes \cite{Middleton2006, Savard2011,Stamper2012Envelope}. The p-type electroreceptor afferents of the tuberous electrosensory system, i.e. the P-units, respond to such AMs of the underlying EOD carrier and their time-dependent firing rate carries information about the stimulus' time-course. P-units are heterogeneous in their baseline firing properties \cite{Grewe2017, Hladnik2023} and differ with respect to their baseline firing rates, the sensitivity and their noisiness. We here explore the nonlinear mechanism in different cells that exhibit distinctly different levels of noise, i.e. differences in the coefficient of variation (CV) of the interspike intervals (ISI). Low-CV P-units have a less noisy, more regular, firing pattern whereas high-CV P-units show a less regular firing pattern in their baseline activity. -\subsection{Low-CV P-units exhibit nonlinear interactions} %frequency combinations withappearing when the input frequencies are related to \fbase{} are -Second-order susceptibility is expected to be especially pronounced for low-CV cells \citealp{Voronenko2017}. P-units fire action potentials probabilistically phase-locked to the self-generated EOD. Skipping of EOD cycles leads to the characteristic multimodal ISI distribution with maxima at integer multiples of the EOD period (\subfigrefb{fig:cells_suscept}{A}). In this example the ISI distribution has a CV of 0.2 which can be considered low among P-units\citealp{Hladnik2023}. Spectral analysis of the baseline activity shows two major peaks, the first is located at the baseline firing rate (\fbase), the second is located at the discharge frequency of the electric organ (\feod{}) and is flanked by two smaller peaks at $\feod \pm \fbase{}$ (\subfigref{fig:cells_suscept}{B}). +\subsection*{Low-CV P-units exhibit nonlinear interactions} %frequency combinations withappearing when the input frequencies are related to \fbase{} are +Second-order susceptibility is expected to be especially pronounced for low-CV cells \cite{Voronenko2017}. P-units fire action potentials probabilistically phase-locked to the self-generated EOD. Skipping of EOD cycles leads to the characteristic multimodal ISI distribution with maxima at integer multiples of the EOD period (\subfigrefb{fig:cells_suscept}{A}). In this example the ISI distribution has a CV of 0.2 which can be considered low among P-units\cite{Hladnik2023}. Spectral analysis of the baseline activity shows two major peaks, the first is located at the baseline firing rate (\fbase), the second is located at the discharge frequency of the electric organ (\feod{}) and is flanked by two smaller peaks at $\feod \pm \fbase{}$ (\subfigref{fig:cells_suscept}{B}). -\begin{figure*}[!h] -\includegraphics{cells_suscept} -\caption{\label{fig:cells_suscept} \notejg{dashed lines still a little faint} Estimation of linear and nonlinear stimulus encoding in a low-CV P-unit. Calculated based on the first frozen noise repeat. \textbf{A} Interspike interval (ISI) distribution of the cell's baseline activity, i.e. the cell is driven only by the unperturbed own field. The CV of the ISIs is a dimensionless measure quantifying the response regularity. Zero CV would indicate perfect regularity. \textbf{B} Power-spectrum of the baseline response. \textbf{C} Random amplitude modulation stimulus (top) and responses (spike raster in the lower traces) of the same P-unit. The stimulus contrast reflects the strength of the AM. \textbf{D} Transfer function (first-order susceptibility, \eqnref{linearencoding_methods}) of the responses to 10\% (light purple) and 20\% contrast (dark purple) RAM stimulation. \textbf{E, F} Absolute value of the second-order susceptibility, \eqnref{eq:susceptibility}, for the low, and high stimulus contrasts. Pink angles -- edges of the structure when \fone, \ftwo{} or \fsum{} are equal to \fbase{}. \textbf{G} Projected diagonals, calculated as the mean of the anti-diagonals of the matrices in \panel{E, F}. Gray dots mark \fbase{}. Horizontal dashed lines: medians of the projected diagonals.} +\begin{figure*}[!ht] +\includegraphics[width=\columnwidth]{cells_suscept} +\caption{\label{fig:cells_suscept} Estimation of linear and nonlinear stimulus encoding in a low-CV P-unit. Calculated based on the first frozen noise repeat. \figitem{A} Interspike interval (ISI) distribution of the cell's baseline activity, i.e. the cell is driven only by the unperturbed own field. The CV of the ISIs is a dimensionless measure quantifying the response regularity. Zero CV would indicate perfect regularity. \figitem{B} Power-spectrum of the baseline response. \figitem{C} Random amplitude modulation stimulus (top) and responses (spike raster in the lower traces) of the same P-unit. The stimulus contrast reflects the strength of the AM. \figitem{D} Transfer function (first-order susceptibility, \Eqnref{linearencoding_methods}) of the responses to 10\% (light purple) and 20\% contrast (dark purple) RAM stimulation. \figitem{E, F} Absolute value of the second-order susceptibility, \Eqnref{eq:susceptibility}, for the low, and high stimulus contrasts. Pink angles -- edges of the structure when \fone, \ftwo{} or \fsum{} are equal to \fbase{}. \figitem{G} Projected diagonals, calculated as the mean of the anti-diagonals of the matrices in \panel{E, F}. Gray dots mark \fbase{}. Horizontal dashed lines: medians of the projected diagonals.} \end{figure*} -Noise stimuli, the random amplitude modulations (RAM, \subfigref{fig:cells_suscept}{C}, top trace, red line) of the EOD, are commonly used to characterize the stimulus driven responses of P-units by means of the transfer function (first-order susceptibility), the spike-triggered average, or the stimulus-response coherence. Here, we additionally estimate the second-order susceptibility to capture nonlinear encoding. Depending on the stimulus intensity which is given as the contrast, i.e. the amplitude of the noise stimulus in relation to the EOD amplitude (see methods), the spikes align more or less clearly to AMs of the stimulus (light and dark purple for low and high contrast stimuli, \subfigrefb{fig:cells_suscept}{C}). The linear encoding (see \eqnref{linearencoding_methods}) is similar for the two RAM contrasts in this low-CV P-unit (\subfigrefb{fig:cells_suscept}{D}). First-order susceptibility is low for low frequencies, peaks in the range below 100\,Hz and then falls off again. +Noise stimuli, the random amplitude modulations (RAM, \subfigref{fig:cells_suscept}{C}, top trace, red line) of the EOD, are commonly used to characterize the stimulus driven responses of P-units by means of the transfer function (first-order susceptibility), the spike-triggered average, or the stimulus-response coherence. Here, we additionally estimate the second-order susceptibility to capture nonlinear encoding. Depending on the stimulus intensity which is given as the contrast, i.e. the amplitude of the noise stimulus in relation to the EOD amplitude (see methods), the spikes align more or less clearly to AMs of the stimulus (light and dark purple for low and high contrast stimuli, \subfigrefb{fig:cells_suscept}{C}). The linear encoding (see \Eqnref{linearencoding_methods}) is similar for the two RAM contrasts in this low-CV P-unit (\subfigrefb{fig:cells_suscept}{D}). First-order susceptibility is low for low frequencies, peaks in the range below 100\,Hz and then falls off again. -Theory predicts a pattern of nonlinear susceptibility for the interaction of two cosine stimuli when \fone{}, \ftwo{} or \fsum{} are equal to \fbase{} (pink triangle in \subfigsref{fig:cells_suscept}{E, F}). To estimate the second-order susceptibility from the P-unit responses to RAM stimuli, the noise stimulus was set in relation to the corresponding neuronal response in the Fourier domain, resulting in a matrix in which the nonlinearity at the sum frequency \fsum{} in the firing rate is depicted for two noise frequencies \fone{} and \ftwo{} (\eqnref{eq:susceptibility}, \subfigrefb{fig:cells_suscept}{E--F}). The decomposition of the RAM stimulus into the frequency components allows us to approximate the beat frequencies $\Delta f$, occurring in case of pure sine-wave stimulation (\subfigrefb{motivation}{D}). Thus, the nonlinearity accessed with the RAM stimulation at \fsum{} (\subfigrefb{fig:cells_suscept}{E}) is similar to the nonlinearity appearing during pure sine-wave stimulation at \bsum{} (orange peak, \subfigrefb{motivation}{D}). A band of elevated nonlinearity appears for the low RAM contrast when \fsumb{} (yellow anti-diagonal between pink edges, \subfigrefb{fig:cells_suscept}{E}) but vanishes for the stronger stimulus (20\% contrast, \subfigref{fig:cells_suscept}{F}). Further, the overall level of nonlinearity changes with the stimulus strength. To compare the structural changes in the matrices we calculated the mean of the anti-diagonals, resulting in the projected diagonal (\subfigrefb{fig:cells_suscept}{G}). For a low RAM contrast the second-order susceptibility at \fbase{} is slightly enhanced (\subfigrefb{fig:cells_suscept}{G}, gray dot on light purple line). For the higher RAM contrast, however, this is not visible and the overall second-order susceptibility is reduced (\subfigrefb{fig:cells_suscept}{G}, compare light and dark purple lines). +Theory predicts a pattern of nonlinear susceptibility for the interaction of two cosine stimuli when \fone{}, \ftwo{} or \fsum{} are equal to \fbase{} (pink triangle in \subfigsref{fig:cells_suscept}{E, F}). To estimate the second-order susceptibility from the P-unit responses to RAM stimuli, the noise stimulus was set in relation to the corresponding neuronal response in the Fourier domain, resulting in a matrix in which the nonlinearity at the sum frequency \fsum{} in the firing rate is depicted for two noise frequencies \fone{} and \ftwo{} (\Eqnref{eq:susceptibility}, \subfigrefb{fig:cells_suscept}{E--F}). The decomposition of the RAM stimulus into the frequency components allows us to approximate the beat frequencies $\Delta f$, occurring in case of pure sine-wave stimulation (\subfigrefb{fig:motivation}{D}). Thus, the nonlinearity accessed with the RAM stimulation at \fsum{} (\subfigrefb{fig:cells_suscept}{E}) is similar to the nonlinearity appearing during pure sine-wave stimulation at \bsum{} (orange peak, \subfigrefb{fig:motivation}{D}). A band of elevated nonlinearity appears for the low RAM contrast when \fsumb{} (yellow anti-diagonal between pink edges, \subfigrefb{fig:cells_suscept}{E}) but vanishes for the stronger stimulus (20\% contrast, \subfigref{fig:cells_suscept}{F}). Further, the overall level of nonlinearity changes with the stimulus strength. To compare the structural changes in the matrices we calculated the mean of the anti-diagonals, resulting in the projected diagonal (\subfigrefb{fig:cells_suscept}{G}). For a low RAM contrast the second-order susceptibility at \fbase{} is slightly enhanced (\subfigrefb{fig:cells_suscept}{G}, gray dot on light purple line). For the higher RAM contrast, however, this is not visible and the overall second-order susceptibility is reduced (\subfigrefb{fig:cells_suscept}{G}, compare light and dark purple lines). -\subsection{High-CV P-units do not exhibit increased nonlinear interactions at \fsum} -CVs in P-units can range up to 1.5 \citealp{Grewe2017, Hladnik2023}. \Figref{fig:cells_suscept_high_CV} shows the same analysis for an example higher-CV P-unit. Similar to the low-CV cell, high-CV P-units fire at multiples of the EOD period (\subfigrefb{fig:cells_suscept_high_CV}{A}). In contrast to low-CV P-unit, however, the higher CV characterizes the noisier baseline firing pattern and the peak at \fbase{} is less pronounced in the power spectrum of the baseline activity (\subfigrefb{fig:cells_suscept_high_CV}{B}). High-CV P-units do not exhibit a clear nonlinear structure related to \fbase{} neither in the second-order susceptibility matrices (\subfigrefb{fig:cells_suscept_high_CV}{E--F}), nor in the projected diagonals (\subfigrefb{fig:cells_suscept_high_CV}{G}). The overall level of nonlinearity, however shows the same dependence on the stimulus contrast. It is much reduced for high-contrast stimuli that drive the neuron much stronger (\subfigrefb{fig:cells_suscept_high_CV}{F}). +\subsection*{High-CV P-units do not exhibit increased nonlinear interactions at \fsum} +\notejg{it is not PLoS style to actively work with the supplementary figures... We may want to reduce this section here, move the text to the supporting information} +CVs in P-units can range up to 1.5 \cite{Grewe2017, Hladnik2023}. \nameref*{S1:highcvpunit} shows the same analysis for an example higher-CV P-unit. Similar to the low-CV cell, high-CV P-units fire at multiples of the EOD period (\subfigrefb{fig:cells_suscept_high_CV}{A}). In contrast to low-CV P-unit, however, the higher CV characterizes the noisier baseline firing pattern and the peak at \fbase{} is less pronounced in the power spectrum of the baseline activity (\subfigrefb{fig:cells_suscept_high_CV}{B}). High-CV P-units do not exhibit a clear nonlinear structure related to \fbase{} neither in the second-order susceptibility matrices (\subfigrefb{fig:cells_suscept_high_CV}{E--F}), nor in the projected diagonals (\subfigrefb{fig:cells_suscept_high_CV}{G}). The overall level of nonlinearity, however shows the same dependence on the stimulus contrast. It is much reduced for high-contrast stimuli that drive the neuron much stronger (\subfigrefb{fig:cells_suscept_high_CV}{F}). -\subsection{Ampullary afferents exhibit strong nonlinear interactions} -Irrespective of the CV, neither cell shows the complete proposed structure of nonlinear interactions. \lepto{} posses an additional electrosensory system, the passive or ampullary electrosensory system that responds to low-frequency exogeneous electric stimuli. The population of ampullary afferents is much less heterogeneous, and known for the much lower CVs of the ISI distributions (0.08--0.22)\citealp{Grewe2017}. Ampullary cells do not phase-lock to the EOD and the ISI are unimodally distributed (\subfigrefb{ampullary}{A}). According to the low irregularity of the baseline response, the power spectrum shows a very distinct peak at \fbase{} and harmonics of it. Since the cells do not fire phase locked to the EOD, there is no peak at \feod{} (\subfigrefb{ampullary}{B}). When driven by a noise stimulus with a low contrast (note: this is not an AM but an added stimulus to the self-generated EOD, \subfigref{ampullary}{C}), ampullary cells exhibit very pronounced bands when \fsum{} is equal to \fbase{} or its harmonic in the second-order susceptibility matrix implying that this cell is especially nonlinear at these frequency combinations (yellow diagonals, \subfigrefb{ampullary}{E, G}, light green). With higher stimulus contrasts these bands disappear (\subfigrefb{ampullary}{F}) and the projected diagonal is decreased and lacks the distinct peak at \fsum{} (\subfigrefb{ampullary}{G}, dark green). +\subsection*{Ampullary afferents exhibit strong nonlinear interactions} +Irrespective of the CV, neither cell shows the complete proposed structure of nonlinear interactions. \lepto{} posses an additional electrosensory system, the passive or ampullary electrosensory system that responds to low-frequency exogeneous electric stimuli. The population of ampullary afferents is much less heterogeneous, and known for the much lower CVs of the ISI distributions (0.08--0.22)\cite{Grewe2017}. Ampullary cells do not phase-lock to the EOD and the ISI are unimodally distributed (\subfigrefb{ampullary}{A}). According to the low irregularity of the baseline response, the power spectrum shows a very distinct peak at \fbase{} and harmonics of it. Since the cells do not fire phase locked to the EOD, there is no peak at \feod{} (\subfigrefb{ampullary}{B}). When driven by a noise stimulus with a low contrast (note: this is not an AM but an added stimulus to the self-generated EOD, \subfigref{ampullary}{C}), ampullary cells exhibit very pronounced bands when \fsum{} is equal to \fbase{} or its harmonic in the second-order susceptibility matrix implying that this cell is especially nonlinear at these frequency combinations (yellow diagonals, \subfigrefb{ampullary}{E, G}, light green). With higher stimulus contrasts these bands disappear (\subfigrefb{ampullary}{F}) and the projected diagonal is decreased and lacks the distinct peak at \fsum{} (\subfigrefb{ampullary}{G}, dark green). -\begin{figure*}[ht]%hp! -\includegraphics{ampullary} - \caption{\label{ampullary} Estimation of linear and nonlinear stimulus encoding in an ampullary afferent. Calculated based on the first frozen noise repeat. \textbf{A} Interspike interval (ISI) distribution of the cell's baseline activity. The CV of the ISIs is a dimensionless measure quantifying the response regularity. Zero CV would indicate perfect regularity. \textbf{B} Power-spectrum of the baseline response. \textbf{C} White noise stimulus (top) added to the fish's self-generated electric field and responses (spike raster in the lower traces). The stimulus contrast reflects the strength of the stimulus in relation to the own electric field amplitude. \textbf{D} Transfer function (first-order susceptibility, \eqnref{linearencoding_methods}) of the responses to 2\% (light green) and 20\% contrast (dark green) stimulation. \textbf{E, F} Absolute value of the second-order susceptibility, \eqnref{eq:susceptibility}, for the low, and high stimulus contrasts. Pink angles -- edges of the structure when \fone, \ftwo{} or \fsum{} are equal to \fbase{}. \textbf{G} Projected diagonals, calculated as the mean of the anti-diagonals of the matrices in \panel{E, F}. Gray dots mark \fbase{}. Horizontal dashed lines: medians of the projected diagonals. +\begin{figure*}[!ht] +\includegraphics[width=\columnwidth]{ampullary} + \caption{\label{ampullary} Estimation of linear and nonlinear stimulus encoding in an ampullary afferent. Calculated based on the first frozen noise repeat. \figitem{A} Interspike interval (ISI) distribution of the cell's baseline activity. The CV of the ISIs is a dimensionless measure quantifying the response regularity. Zero CV would indicate perfect regularity. \figitem{B} Power-spectrum of the baseline response. \figitem{C} White noise stimulus (top) added to the fish's self-generated electric field and responses (spike raster in the lower traces). The stimulus contrast reflects the strength of the stimulus in relation to the own electric field amplitude. \figitem{D} Transfer function (first-order susceptibility, \Eqnref{linearencoding_methods}) of the responses to 2\% (light green) and 20\% contrast (dark green) stimulation. \figitem{E, F} Absolute value of the second-order susceptibility, \Eqnref{eq:susceptibility}, for the low, and high stimulus contrasts. Pink angles -- edges of the structure when \fone, \ftwo{} or \fsum{} are equal to \fbase{}. \figitem{G} Projected diagonals, calculated as the mean of the anti-diagonals of the matrices in \panel{E, F}. Gray dots mark \fbase{}. Horizontal dashed lines: medians of the projected diagonals. } \end{figure*} -% \subsection{Full nonlinear structure visible only in P-unit models} -\subsection{Internal noise hides parts of the nonlinearity structure} -Traces of the proposed structure of second-order susceptibility are found in both ampullary and p-type electrosensory afferents. The nonlinerarity seems to depend on the CV, i.e. the level of intrinsic noise in the cells and in the next step we address whether the level of intrinsic masks the expected structures. One option to reduce the impact of this intrinsic noise is to average it out over many repetitions of the same stimulus. In the electrophysiological experiments we only have a limited recording duration and hence the number of stimulus repetitions is limited. To overcome these limitations we compared between an electrophysiologically recorded low-CV P-unit and its P-unit LIF model counterpart fitted to reproduce the P-unit behavior faithfully (see \figref{flowchart})\citealp{Barayeu2023}. In the recording depicted in \subfigrefb{model_and_data}{A} the cell was stimulated with a weak RAM stimulus with the contrast of 1\,$\%$ (red, \subfigrefb{model_and_data}\,\panel[i]{B}). The across trial average ($\n{}=11$) shows the diagonal band in second-order susceptibility at \fsumb{} (yellow diagonal in pink edges) and an additional nonlinear band appeared at \fsumehalf{} (\subfigrefb{model_and_data}{A}). The matched model reproduces the same diagonal at \fsumb{} (\subfigrefb{model_and_data}\,\panel[ii]{B}) but not the diagonal at \fsumehalf. By increasing the number of trials in the simulation, the signal-to-noise ratio and thus the estimation of the nonlinearity structures can be improved. Still, even with 1 million repetitions, no changes are observable in the nonlinearity structures (\subfigrefb{model_and_data}\,\panel[iii]{B}). The increased trial count, however goes along with generally smaller second-order susceptibility values (\subfigrefb{model_and_data}\,\panel[iii]{B}). +\subsection*{Internal noise hides parts of the nonlinearity structure} +Traces of the proposed structure of second-order susceptibility are found in both ampullary and p-type electrosensory afferents. The nonlinerarity seems to depend on the CV, i.e. the level of intrinsic noise in the cells and in the next step we address whether the level of intrinsic masks the expected structures. One option to reduce the impact of this intrinsic noise is to average it out over many repetitions of the same stimulus. In the electrophysiological experiments we only have a limited recording duration and hence the number of stimulus repetitions is limited. To overcome these limitations we compared between an electrophysiologically recorded low-CV P-unit and its P-unit LIF model counterpart fitted to reproduce the P-unit behavior faithfully (see \figref{flowchart})\cite{Barayeu2023}. In the recording depicted in \subfigrefb{model_and_data}{A} the cell was stimulated with a weak RAM stimulus with the contrast of 1\,$\%$ (red, \subfigrefb{model_and_data}\,\panel[i]{B}). The across trial average ($\n{}=11$) shows the diagonal band in second-order susceptibility at \fsumb{} (yellow diagonal in pink edges) and an additional nonlinear band appeared at \fsumehalf{} (\subfigrefb{model_and_data}{A}). The matched model reproduces the same diagonal at \fsumb{} (\subfigrefb{model_and_data}\,\panel[ii]{B}) but not the diagonal at \fsumehalf. By increasing the number of trials in the simulation, the signal-to-noise ratio and thus the estimation of the nonlinearity structures can be improved. Still, even with 1 million repetitions, no changes are observable in the nonlinearity structures (\subfigrefb{model_and_data}\,\panel[iii]{B}). The increased trial count, however goes along with generally smaller second-order susceptibility values (\subfigrefb{model_and_data}\,\panel[iii]{B}). %Note that this line doesn't appear in the susceptibility matrix of the model (\subfigrefb{model_and_data}{C})Each cell has an intrinsic noise level (\subfigrefb{model_and_data}{A}, bottom). %The signal component (purple) compensates for this total noise reduction, is not a weak signal anymore and -Based on the Novikov-Furutsu Theorem \citealp{Novikov1965, Furutsu1963} the intrinsic noise of a LIF model can be split up into several independent noise processes with the same correlation function. We make use of this and split the intrinsic noise $\xi$ into two parts: 90\% are then treated as signal ($\xi_{signal}$) while the remaining 10\% are treated as noise ($\xi_{noise}$, see methods for details). With this, the signal-to-noise ratio in the simulation can be arbitrarily varied and the combination of many repetitions and noise-split indeed reveals the triangular shape shown theoretically and for LIF models without carrier\citealp{Voronenko2017}(\subfigrefb{model_and_data}\,\panel[iii]{C}). +Based on the Novikov-Furutsu Theorem \cite{Novikov1965, Furutsu1963} the intrinsic noise of a LIF model can be split up into several independent noise processes with the same correlation function. We make use of this and split the intrinsic noise $\xi$ into two parts: 90\% are then treated as signal ($\xi_{signal}$) while the remaining 10\% are treated as noise ($\xi_{noise}$, see methods for details). With this, the signal-to-noise ratio in the simulation can be arbitrarily varied and the combination of many repetitions and noise-split indeed reveals the triangular shape shown theoretically and for LIF models without carrier\cite{Voronenko2017}(\subfigrefb{model_and_data}\,\panel[iii]{C}). %Adding an additional, independent, RAM stimulus to the simulation does not heavily influence the qualitative observations but the nonlinearity becomes weaker (compare \subfigrefb{model_and_data}\,\panel[iii]{C} and \panel[iv]{C}) -% Based on this a weak RAM signal as the input to the P-unit model (red in \subfigrefb{model_and_data}\,\panel[i]{A}) can be approximated by a model where no RAM stimulus is present (red) but instead the total noise of the model is split up into a reduced intrinsic noise component (gray) and a signal component (purple), maintaining the CV and \fbase{} as during baseline (\subfigrefb{model_and_data}\,\panel[i]{B}, see methods section \ref{intrinsicsplit_methods} for more details). This signal component (purple) can be used for the calculation of the second-order susceptibility. With the reduced noise component the signal-to-noise ratio increases and the number of stimulus realizations can be reduced. This noise split cannot be applied in experimentally measured cells. If this noise split is applied in the model with $\n{}=11$ stimulus realizations the nonlinearity at \fsumb{} is still present (\subfigrefb{model_and_data}\,\panel[iii]{B}). If instead, the RAM stimulus is drawn 1 million times the diagonal nonlinearity band is complemented by vertical and horizontal lines appearing at \foneb{} and \ftwob{} (\subfigrefb{model_and_data}\,\panel[iv]{B}). These nonlinear structures correspond to the ones observed in previous works \citealp{Voronenko2017}. If now a weak RAM stimulus is added (\subfigrefb{model_and_data}\,\panel[i]{C}, red), simultaneously the noise is split up into a noise and signal component (gray and purple) and the calculation is performed on the sum of the signal component and the RAM (red plus purple) only the diagonal band is present with $\n{}=11$ (\subfigrefb{model_and_data}\,\panel[iii]{C}) but also with 1\,million stimulus realizations (\subfigrefb{model_and_data}\,\panel[iv]{C}). +% Based on this a weak RAM signal as the input to the P-unit model (red in \subfigrefb{model_and_data}\,\panel[i]{A}) can be approximated by a model where no RAM stimulus is present (red) but instead the total noise of the model is split up into a reduced intrinsic noise component (gray) and a signal component (purple), maintaining the CV and \fbase{} as during baseline (\subfigrefb{model_and_data}\,\panel[i]{B}, see methods section \ref{intrinsicsplit_methods} for more details). This signal component (purple) can be used for the calculation of the second-order susceptibility. With the reduced noise component the signal-to-noise ratio increases and the number of stimulus realizations can be reduced. This noise split cannot be applied in experimentally measured cells. If this noise split is applied in the model with $\n{}=11$ stimulus realizations the nonlinearity at \fsumb{} is still present (\subfigrefb{model_and_data}\,\panel[iii]{B}). If instead, the RAM stimulus is drawn 1 million times the diagonal nonlinearity band is complemented by vertical and horizontal lines appearing at \foneb{} and \ftwob{} (\subfigrefb{model_and_data}\,\panel[iv]{B}). These nonlinear structures correspond to the ones observed in previous works \cite{Voronenko2017}. If now a weak RAM stimulus is added (\subfigrefb{model_and_data}\,\panel[i]{C}, red), simultaneously the noise is split up into a noise and signal component (gray and purple) and the calculation is performed on the sum of the signal component and the RAM (red plus purple) only the diagonal band is present with $\n{}=11$ (\subfigrefb{model_and_data}\,\panel[iii]{C}) but also with 1\,million stimulus realizations (\subfigrefb{model_and_data}\,\panel[iv]{C}). -In a high-CV P-unit we could not find such nonlinear structures, neither in the electrophysiologically recorded data nor in the respective model (not shown), corresponding to the theoretical predictions \citealp{Voronenko2017}. +In a high-CV P-unit we could not find such nonlinear structures, neither in the electrophysiologically recorded data nor in the respective model (not shown), corresponding to the theoretical predictions \cite{Voronenko2017}. -\begin{figure*}[hb!] +\begin{figure*}[!hb] \includegraphics[width=\columnwidth]{model_and_data} \caption{\label{model_and_data} High trial numbers and a reduced internal noise reveal the nonlinear structure in a LIF model with carrier. \figitem{A} \suscept{} surface of the electrophysiologically recorded low-CV P-unit (cell 2012-07-03-ak) driven with a weak stimulus with 1\% contrast. The plot shows that average \suscept{} surface ($N=11$). Pink edges mark the expected structure of enhanced nonlinearity. \figitem[i]{B} \textit{Standard condition} of the model simulation. The model is driven by a RAM stimulus (red trace). $\xi_{signal}$ and $\xi_{noise}$ components of the total intrinsic noise ($\xi$) are shown in the center and bottom. In this condition, $\xi_{noise} = \xi$. \figitem[ii]{B} \suscept{} surface of the cell's LIF model counterpart (cell 2012-07-03-ak, table~\ref{modelparams}) averaged over the same number of stimulus repetitions. \figitem[iii]{B} Same as \panel[ii]{B} but 1 million stimulus repetitions. \figitem[i - iii]{C} Same as in \panel[i - iii]{B} but in the \textit{noise split} condition: there is no external RAM signal driving the model. Parts (90\%) of the total intrinsic noise ($\xi$) are treated as signal (center trace, $\xi_{signal}$) while the remaining part is treated as $\xi_{noise}$ (see methods for details). } \end{figure*} -\subsection{Second-order susceptibility can explain nonlinear peaks in pure sinewave stimulation} -We calculated the second-order susceptibility surfaces at \fsum{} by extracting the respective spectral component of the P-unit responses to RAM stimuli. How does this relate to the pure sinewave situation that approximates the interference of real animals? The relevant signals are the the beat frequencies \bone{} and \btwo{} that arise from the interference of the receiving fish's own and either of the foreign EODs (\figref{motivation}). The response power spectrum showed peaks from nonlinear interaction at the sum of the two beat frequencies (orange marker, \subfigrefb{motivation}{D}) and at the difference between the two beat frequencies (rosa marker, \subfigrefb{motivation}{D}). In the example $\Delta f_{2}$ was similar to \fbase{}, corresponding to the horizontal line in the RAM-based second-order susceptibility matrix in \subfigrefb{model_and_data}\,\panel[iii]{C}. In the three-fish example, there was a second, less prominent, nonlinearity at the difference between the two beat frequencies (rosa circle, \subfigrefb{motivation}{D}), that cannot be explained with the second-order susceptibility matrix in \subfigrefb{model_and_data}\,\panel[iii]{C}, in which only the nonlinearity at \fsum{} in the response is addressed (\eqnref{eq:crosshigh}). -\begin{figure*}[h] +\subsection*{Second-order susceptibility can explain nonlinear peaks in pure sinewave stimulation} +We calculated the second-order susceptibility surfaces at \fsum{} by extracting the respective spectral component of the P-unit responses to RAM stimuli. How does this relate to the pure sinewave situation that approximates the interference of real animals? The relevant signals are the the beat frequencies \bone{} and \btwo{} that arise from the interference of the receiving fish's own and either of the foreign EODs (\figref{fig:motivation}). The response power spectrum showed peaks from nonlinear interaction at the sum of the two beat frequencies (orange marker, \subfigrefb{fig:motivation}{D}) and at the difference between the two beat frequencies (rosa marker, \subfigrefb{fig:motivation}{D}). In the example $\Delta f_{2}$ was similar to \fbase{}, corresponding to the horizontal line in the RAM-based second-order susceptibility matrix in \subfigrefb{model_and_data}\,\panel[iii]{C}. In the three-fish example, there was a second, less prominent, nonlinearity at the difference between the two beat frequencies (rosa circle, \subfigrefb{fig:motivation}{D}), that cannot be explained with the second-order susceptibility matrix in \subfigrefb{model_and_data}\,\panel[iii]{C}, in which only the nonlinearity at \fsum{} in the response is addressed (\Eqnref{eq:crosshigh}). +\begin{figure*}[!ht] \includegraphics[width=\columnwidth]{model_full} - \caption{\label{model_full} Second-order susceptibility of the model of an electrophysiologically recorded P-unit (see table~\ref{modelparams} for model parameters of 2012-07-03-ak). White lines -- coordinate axis. The nonlinearity at \fsum{} in the firing rate is quantified in the upper right and lower left quadrants (\eqnref{eq:susceptibility}). The nonlinearity at \fdiff{} in the firing rate is quantified in the upper left and lower right quadrants. Baseline firing rate $\fbase{}=120$\,Hz. Absolute value of the model second-order susceptibility, with $N=1$\,million stimulus realizations and the intrinsic noise split (see methods). The colored markers highlight the nonlinear effects found in \subfigrefb{motivation}{D}.} + \caption{\label{fig:model_full} Second-order susceptibility of the model of an electrophysiologically recorded P-unit (see table~\ref{modelparams} for model parameters of 2012-07-03-ak). White lines -- coordinate axis. The nonlinearity at \fsum{} in the firing rate is quantified in the upper right and lower left quadrants (\Eqnref{eq:susceptibility}). The nonlinearity at \fdiff{} in the firing rate is quantified in the upper left and lower right quadrants. Baseline firing rate $\fbase{}=120$\,Hz. Absolute value of the model second-order susceptibility, with $N=1$\,million stimulus realizations and the intrinsic noise split (see methods). The colored markers highlight the nonlinear effects found in \subfigrefb{fig:motivation}{D}.} \end{figure*} %section \ref{intrinsicsplit_methods}).\figitem{A} Absolute value of the second-order susceptibility of an electrophysiologically recorded P-unit. RAM stimulus realizations $N=11$. Diagonal bands appear when the sum of the frequencies \fsum{} or the difference \fdiff{} is equal to \fbase{}. \figitem{B} The diagonals, that were present in \panel{A}, are complemented by vertical and horizontal lines when \fone{} or \ftwo{} are equal to \fbase{}. Note that the different scale of the second-order susceptibility is associated with the higher signal-to-ratio in case of 1 million repeats in \panel{B}. -The second-order susceptibility can also be calculated for the full matrix including negative frequency components of the respective spectra in \eqnsref{eq:crosshigh} and (\ref{eq:susceptibility}). The resulting \suscept{} matrix is symmetric with respect to the origin and shows \suscept{} at \fsum{} in the upper-right and lower-left quadrants and \suscept{} for the differences \fdiff{} in the lower-right and upper-left quadrants \citealp{Voronenko2017} (\figref{model_full}). The vertical and horizontal lines at \foneb{} and \ftwob{} are very pronounced in the upper-right quadrant of \figref{model_full} for the nonlinearity at \fsum{} and extend into the upper-left quadrant (representing \fdiff) fading out towards more negative $f_1$ frequencies. \suscept{} values at \fsum{} match the \fsum{} peak seen in \figref{motivation}. +The second-order susceptibility can also be calculated for the full matrix including negative frequency components of the respective spectra in \Eqnsref{eq:crosshigh} and (\ref{eq:susceptibility}). The resulting \suscept{} matrix is symmetric with respect to the origin and shows \suscept{} at \fsum{} in the upper-right and lower-left quadrants and \suscept{} for the differences \fdiff{} in the lower-right and upper-left quadrants \cite{Voronenko2017} (\figref{fig:model_full}). The vertical and horizontal lines at \foneb{} and \ftwob{} are very pronounced in the upper-right quadrant of \figref{fig:model_full} for the nonlinearity at \fsum{} and extend into the upper-left quadrant (representing \fdiff) fading out towards more negative $f_1$ frequencies. \suscept{} values at \fsum{} match the \fsum{} peak seen in \figref{fig:motivation}. -The smaller \fdiff{} power spectral peak observed during pure sine-wave stimulation (\subfigrefb{motivation}{D}) can be explained by the fading out horizontal line in the upper-left quadrant (\figrefb{model_full}, \citealp{Schlungbaum2023}). Even though the second-order susceptibilities here were estimated form data and models with an modulated (EOD) carrier (\figrefb{model_full}) they are in good accordance with the second-order susceptibilities found in LIF models without a carrier\citealp{Voronenko2017, Schlungbaum2023}. +The smaller \fdiff{} power spectral peak observed during pure sine-wave stimulation (\subfigrefb{fig:motivation}{D}) can be explained by the fading out horizontal line in the upper-left quadrant (\figrefb{fig:model_full}, \cite{Schlungbaum2023}). Even though the second-order susceptibilities here were estimated form data and models with an modulated (EOD) carrier (\figrefb{fig:model_full}) they are in good accordance with the second-order susceptibilities found in LIF models without a carrier\cite{Voronenko2017, Schlungbaum2023}. -\begin{figure*}[ht!] +\begin{figure*}[!ht] \includegraphics[width=\columnwidth]{data_overview_mod} - \caption{\label{fig:data_overview_mod} Nonlinearity for a population of P-units (\panel{A, C}) and ampullary cells (\panel{B, D}). The nonlinearity index \nli{} is calculated as the maximal value in the range $\fbase{} \pm 5$\,Hz of the projected diagonal divided by its median (see \eqnref{eq:nli_equation}). Each recorded neuron contributes at maximum with two stimulus contrasts. \figitem{A, B} There is a negative correlation between the CV during baseline and \nli. \figitem{C, D} The \nli{} is plotted against the response modulation, (see methods), an indicator of the subjective stimulus strength for a cell. There is a negative correlation between response modulation and \nli. Restricting the analysis to the weakest stimulus that was presented to each unique neuron, does not change the results. The number of unique neurons is 222 for P-units and 45 for ampullary cells. + \caption{\label{fig:data_overview_mod} Nonlinearity for a population of P-units (\panel{A, C}) and ampullary cells (\panel{B, D}). The nonlinearity index \nli{} is calculated as the maximal value in the range $\fbase{} \pm 5$\,Hz of the projected diagonal divided by its median (see \Eqnref{eq:nli_equation}). Each recorded neuron contributes at maximum with two stimulus contrasts. \figitem{A, B} There is a negative correlation between the CV during baseline and \nli. \figitem{C, D} The \nli{} is plotted against the response modulation, (see methods), an indicator of the subjective stimulus strength for a cell. There is a negative correlation between response modulation and \nli. Restricting the analysis to the weakest stimulus that was presented to each unique neuron, does not change the results. The number of unique neurons is 222 for P-units and 45 for ampullary cells. } \end{figure*} -%\eqnref{response_modulation} +%\Eqnref{response_modulation} -\subsection{Low CVs and weak stimuli are associated with strong nonlinearity} -The nonlinear effects shown for single cell examples above are supported by the analysis of the pool of 222 P-units and 45 ampullary afferents recorded in 71 specimen. To compare across cells we expressed the second-order susceptibility as the nonlinearity index \nli{}, see \eqnref{eq:nli_equation}. The \nli{} characterizes the peakedness of the projections onto the diagonal of the \suscept{} matrices (e.g. \subfigref{fig:cells_suscept}{G} at \fbase{}). It assumes high values when the \fbase{} peak in the projected diagonal is pronounced relative to the median of the diagonal and is small when there is no distinct peak. The P-unit population \nli{} values depend weakly on the CV of the ISI distribution under baseline condition. Cells with lower baseline CVs have the tendency to exhibit a stronger nonlinearity than those that have high CVs (Pearson's $r=-0.16$, $p<0.01$). The two example P-units shown before (\figrefb{fig:cells_suscept} and \figrefb{fig:cells_suscept_high_CV}) are highlighted with dark circles in \subfigrefb{fig:data_overview_mod}{A, C}. The stimulus strength plays an important role. Several of the recorded neurons contribute with two dots to the data as their responses to the same RAM stimulus but with different contrasts were recorded. Higher stimulus contrasts lead to a stronger drive and thus stronger response modulations (see color code bar in \subfigref{fig:data_overview_mod}{A}, see methods). Since P-units are heterogeneous in their susceptibility to the stimulus\citealp{Grewe2017}, their responses to the same stimulus intensity vary a lot. When replotting the \nli{} against the response modulation (\subfigrefb{fig:data_overview_mod}{C}) a negative correlation is observed showing that cells that are strongly driven by the stimulus show less nonlinearity while those that are only weakly driven show higher nonlinearities. Whether or not a cell responds nonlinearly thus depends on both, the baseline CV (i.e. the internal noise) and the response strength. +\subsection*{Low CVs and weak stimuli are associated with strong nonlinearity} +The nonlinear effects shown for single cell examples above are supported by the analysis of the pool of 222 P-units and 45 ampullary afferents recorded in 71 specimen. To compare across cells we expressed the second-order susceptibility as the nonlinearity index \nli{}, see \Eqnref{eq:nli_equation}. The \nli{} characterizes the peakedness of the projections onto the diagonal of the \suscept{} matrices (e.g. \subfigref{fig:cells_suscept}{G} at \fbase{}). It assumes high values when the \fbase{} peak in the projected diagonal is pronounced relative to the median of the diagonal and is small when there is no distinct peak. The P-unit population \nli{} values depend weakly on the CV of the ISI distribution under baseline condition. Cells with lower baseline CVs have the tendency to exhibit a stronger nonlinearity than those that have high CVs (Pearson's $r=-0.16$, $p<0.01$). The two example P-units shown before (\figrefb{fig:cells_suscept} and \figrefb{fig:cells_suscept_high_CV}) are highlighted with dark circles in \subfigrefb{fig:data_overview_mod}{A, C}. The stimulus strength plays an important role. Several of the recorded neurons contribute with two dots to the data as their responses to the same RAM stimulus but with different contrasts were recorded. Higher stimulus contrasts lead to a stronger drive and thus stronger response modulations (see color code bar in \subfigref{fig:data_overview_mod}{A}, see methods). Since P-units are heterogeneous in their susceptibility to the stimulus\cite{Grewe2017}, their responses to the same stimulus intensity vary a lot. When replotting the \nli{} against the response modulation (\subfigrefb{fig:data_overview_mod}{C}) a negative correlation is observed showing that cells that are strongly driven by the stimulus show less nonlinearity while those that are only weakly driven show higher nonlinearities. Whether or not a cell responds nonlinearly thus depends on both, the baseline CV (i.e. the internal noise) and the response strength. %(Pearson's $r=-0.35$, $p<0.001$) %In a P-unit population where each cell is represented not by several contrasts but by the lowest recorded contrast, \nli{} significantly correlates with the CV during baseline ($r=-0.17$, $p=0.01$), the response modulation ($r=-0.35$, $p<0.001$) and \fbase{} ($r=-0.32$, $p<0.001$).%, $\n{}=221$*, $\n{}=221$******, $\n{}=221$ @@ -501,48 +584,46 @@ The population of ampullary cells is generally more homogeneous and have lower C %(Pearson's $r=-0.35$, $p < 0.01$) (Pearson's $r=-0.59$, $p < 0.0001$) -\section{Discussion} +\section*{Discussion} -Nonlinearities are ubiquitous in nervous systems, they are essential to extract certain features such as amplitude modulations of a carrier. Here, we analyzed the nonlinearity in primary electroreceptor afferents of the weakly electric fish \lepto{}. Under natural stimulus conditions, when the electric fields of two or more animals produce interference patterns that contain information of the context of the encounter, such nonlinear encoding is required to enable responding to AMs or also so-called envelopes. The work presented here is inspired by observations of electric fish interactions in the natural habitat. In these observed electrosensory cocktail-parties, a second male is intruding the conversations of a courting dyad. The courting male detects the intruder at a large distance where the foreign electric signal is strongly attenuated due to spatial distance and tiny compared to the signal emitted by the close-by female. Such a three-animal situation is exemplified in \figref{motivation} for a rather specific combination of interfering EODs. There, we observe that the electroreceptor response contains components that result from nonlinear interference (\subfigref{motivation}{D}) that might help to solve this detection task. Based on theoretical work we work out the circumstances under which electroreceptor afferents show such nonlinearities. +Nonlinearities are ubiquitous in nervous systems, they are essential to extract certain features such as amplitude modulations of a carrier. Here, we analyzed the nonlinearity in primary electroreceptor afferents of the weakly electric fish \lepto{}. Under natural stimulus conditions, when the electric fields of two or more animals produce interference patterns that contain information of the context of the encounter, such nonlinear encoding is required to enable responding to AMs or also so-called envelopes. The work presented here is inspired by observations of electric fish interactions in the natural habitat. In these observed electrosensory cocktail-parties, a second male is intruding the conversations of a courting dyad. The courting male detects the intruder at a large distance where the foreign electric signal is strongly attenuated due to spatial distance and tiny compared to the signal emitted by the close-by female. Such a three-animal situation is exemplified in \figref{fig:motivation} for a rather specific combination of interfering EODs. There, we observe that the electroreceptor response contains components that result from nonlinear interference (\subfigref{fig:motivation}{D}) that might help to solve this detection task. Based on theoretical work we work out the circumstances under which electroreceptor afferents show such nonlinearities. \subsection*{Theory applies to systems with and without carrier} -Theoretical work\citealp{Voronenko2017} explained analytically the occurrence of nonlinear products when a LIF model neuron is stimulated with pure sine-waves. To investigate whether the same mechanisms occur in electroreceptor afferents which are driven by AMs of a carrier and not by pure sine-waves, we followed the previous approach and quantified the second-order susceptibility from responses to white-noise stimuli\citealp{Voronenko2017,Neiman2011fish,Nikias1993}. We expected to see elevated second-order susceptibility where either of the foreign signals matches the baseline firing rate ($f_1=f_{base}$ or $f_2=f_{base}$) or when the sum equals the baseline firing rate of the neuron (\fsumb{}) creating a triangular pattern of elevated \suscept{} e.g.\,\subfigref{model_and_data}{C}{iii}. Indeed, we find traces of the same nonlinearities in the neuronal responses of p-type electroreceptor afferents. The nonlinear pattern observed in the experimental data, however, matches to the expectations only partially and only in a subset of neurons \figsref{fig:cells_suscept} and\,\ref{ampullary}. Nevertheless, the theory holds also for systems that are driven by AMs of a carrier and is thus more widely applicable. +Theoretical work\cite{Voronenko2017} explained analytically the occurrence of nonlinear products when a LIF model neuron is stimulated with pure sine-waves. To investigate whether the same mechanisms occur in electroreceptor afferents which are driven by AMs of a carrier and not by pure sine-waves, we followed the previous approach and quantified the second-order susceptibility from responses to white-noise stimuli\cite{Voronenko2017,Neiman2011fish,Nikias1993}. We expected to see elevated second-order susceptibility where either of the foreign signals matches the baseline firing rate ($f_1=f_{base}$ or $f_2=f_{base}$) or when the sum equals the baseline firing rate of the neuron (\fsumb{}) creating a triangular pattern of elevated \suscept{} e.g.\,\subfigref{model_and_data}{C}{iii}. Indeed, we find traces of the same nonlinearities in the neuronal responses of p-type electroreceptor afferents. The nonlinear pattern observed in the experimental data, however, matches to the expectations only partially and only in a subset of neurons \figsref{fig:cells_suscept} and\,\ref{ampullary}. Nevertheless, the theory holds also for systems that are driven by AMs of a carrier and is thus more widely applicable. \subsection*{Intrinsic noise limits nonlinear responses} Only those P-units that exhibit low coefficients of variation (CV) of the interspike-interval distribution (\figref{fig:cells_suscept}{A}) in their unperturbed baseline response show the expected nonlinerities (\subfigref{fig:data_overview_mod}{A}). Such low-CV cells are rare among the 221 P-units that we used in this study. The afferents of the passive electrosensory system, the ampullary cells, however have generally lower CVs and show a much clearer nonlinearity pattern than the low-CV P-unit exemplified here (compare \figsref{fig:cells_suscept} and \ref{ampullary}). The single ampullary cell featured in \figref{ampullary} is a representative of the majority of cells analyzed here. All ampullary cells have CVs below 0.4 with a median around 0.1 and the observed \nli{}s are 10-fold higher than in P-units. -The CV serves as a proxy for the intrinsic noise in the cells. In both cell types, we observe a negative correlation between \nli{} and the CV, indicating that it is the level of intrinsic noise that plays a role here. These findings are in line with previous studies that propose that noise linearizes the system\citealp{Roddey2000, Chialvo1997, Voronenko2017}. More intrinsic noise has been demonstrated to increase the CV and reduce nonlinear phase-locking in vestibular afferents\citealp{Schneider2011}. Reduced noise, on the other hand, has been associated with stronger nonlinearity in pyramidal cells of the ELL\citealp{Chacron2006}. Further support for the notion of noise limiting the nonlinearity comes from our P-unit LIF model that faithfully reproduces P-unit activity\citealp{Barayeu2023}. We can use this model and the noise-split according to the Novikov-Furutsu theorem\citealp{Novikov1965,Furutsu1963} to increase the signal-to-noise ratio in the cell while keeping the overall response variability constant (see methods). Treating 90\% of the total noise as signal and simulating large numbers of trial uncovers the full nonlinearity structure (\figref{model_and_data}) seen in LIF neurons and the analytical derivations when driven with sine-wave stimuli\citealp{Voronenko2017}. +The CV serves as a proxy for the intrinsic noise in the cells. In both cell types, we observe a negative correlation between \nli{} and the CV, indicating that it is the level of intrinsic noise that plays a role here. These findings are in line with previous studies that propose that noise linearizes the system\cite{Roddey2000, Chialvo1997, Voronenko2017}. More intrinsic noise has been demonstrated to increase the CV and reduce nonlinear phase-locking in vestibular afferents\cite{Schneider2011}. Reduced noise, on the other hand, has been associated with stronger nonlinearity in pyramidal cells of the ELL\cite{Chacron2006}. Further support for the notion of noise limiting the nonlinearity comes from our P-unit LIF model that faithfully reproduces P-unit activity\cite{Barayeu2023}. We can use this model and the noise-split according to the Novikov-Furutsu theorem\cite{Novikov1965,Furutsu1963} to increase the signal-to-noise ratio in the cell while keeping the overall response variability constant (see methods). Treating 90\% of the total noise as signal and simulating large numbers of trial uncovers the full nonlinearity structure (\figref{model_and_data}) seen in LIF neurons and the analytical derivations when driven with sine-wave stimuli\cite{Voronenko2017}. -\subsection{Noise stimulation approximates the real three-fish interaction} -Our analysis is based on the neuronal responses to white noise stimulus sequences. For the P-units, the stimulus was a random amplitude modulation (RAM) while it was a direct noise stimulus for the ampullary cells. In true white noise, each frequency component contributes with equal amplitude. This has the methodological advantage that all behaviorally relevant frequencies can be measured with a single stimulus presentation and is a widely used approach to characterize sensory coding\citealp{Borst1999}\notejg{better references?}. In contrast to the situation with individual frequencies (direct sine-waves or sinusoidal AMs) the total power of the stimulus is equally distributed on all frequencies leading to a weaker signal-to-noise ratio. This explains that the nonlinearity pattern in the electroreceptor recordings only partially matches the expectation \figsref{fig:cells_suscept},\,\ref{ampullary} while the single-frequency stimulation shows nonlinear interference when the individual stimulus frequencies ($f_1, f_2, \Delta f_1, \Delta f_2$) match the baseline firing rate (\figref{motivation}). With the noise-splitting trick, we could show that the full nonlinearity pattern is present but covered by the intrinsic noise. We thus conclude that the presence of the anti-diagonal pattern in the \suscept{} matrix is sufficient to conclude that the same nonlinear interactions happen here. This also validates the application of the white noise approach to characterize the full \suscept{} matrix instead of using all combinations of individual frequencies. +\subsection*{Noise stimulation approximates the real three-fish interaction} +Our analysis is based on the neuronal responses to white noise stimulus sequences. For the P-units, the stimulus was a random amplitude modulation (RAM) while it was a direct noise stimulus for the ampullary cells. In true white noise, each frequency component contributes with equal amplitude. This has the methodological advantage that all behaviorally relevant frequencies can be measured with a single stimulus presentation and is a widely used approach to characterize sensory coding\cite{Borst1999}\notejg{better references?}. In contrast to the situation with individual frequencies (direct sine-waves or sinusoidal AMs) the total power of the stimulus is equally distributed on all frequencies leading to a weaker signal-to-noise ratio. This explains that the nonlinearity pattern in the electroreceptor recordings only partially matches the expectation \figsref{fig:cells_suscept},\,\ref{ampullary} while the single-frequency stimulation shows nonlinear interference when the individual stimulus frequencies ($f_1, f_2, \Delta f_1, \Delta f_2$) match the baseline firing rate (\figref{fig:motivation}). With the noise-splitting trick, we could show that the full nonlinearity pattern is present but covered by the intrinsic noise. We thus conclude that the presence of the anti-diagonal pattern in the \suscept{} matrix is sufficient to conclude that the same nonlinear interactions happen here. This also validates the application of the white noise approach to characterize the full \suscept{} matrix instead of using all combinations of individual frequencies. -% The nonlinearity of ampullary cells in paddlefish \citealp{Neiman2011fish} has been previously accessed with bandpass limited white noise. +% The nonlinearity of ampullary cells in paddlefish \cite{Neiman2011fish} has been previously accessed with bandpass limited white noise. -% Here it was demonstrated that the second-order susceptibility for the two RAM noise input frequencies \fone{} and \ftwo{} can approximate a three-fish setting, where the driving force for the P-unit are two beats with frequencies \bone{} and \btwo{}. This was confirmed by a low-CV P-unit, where nonlinearities in the P-unit response occurred at the sum and difference of the beat frequencies for pure sine-wave stimulation (\figrefb{motivation}). In this P-unit the nonlinearity appeared in a three-wave setting with \bone{} being close to \fbase{}, corresponding to a frequency combination on the vertical line at \foneb{} in the second-order susceptibility (\subfigrefb{model_full}{B}). This implies that even if only the diagonal structure can be accessed with noise stimulation in the second-order susceptibility in an experiment (\subfigrefb{model_full}{A}), it can be taken as an indicator that the whole nonlinear structure should be present during pure sine-wave stimulation in the same cell. With this RAM stimulation was demonstrated to be an effective method to scan the three-fish or two-beat plane and estimate the whole theoretically predicted nonlinear structure in experimentally recorded cells. +% Here it was demonstrated that the second-order susceptibility for the two RAM noise input frequencies \fone{} and \ftwo{} can approximate a three-fish setting, where the driving force for the P-unit are two beats with frequencies \bone{} and \btwo{}. This was confirmed by a low-CV P-unit, where nonlinearities in the P-unit response occurred at the sum and difference of the beat frequencies for pure sine-wave stimulation (\figrefb{fig:motivation}). In this P-unit the nonlinearity appeared in a three-wave setting with \bone{} being close to \fbase{}, corresponding to a frequency combination on the vertical line at \foneb{} in the second-order susceptibility (\subfigrefb{fig:model_full}{B}). This implies that even if only the diagonal structure can be accessed with noise stimulation in the second-order susceptibility in an experiment (\subfigrefb{fig:model_full}{A}), it can be taken as an indicator that the whole nonlinear structure should be present during pure sine-wave stimulation in the same cell. With this RAM stimulation was demonstrated to be an effective method to scan the three-fish or two-beat plane and estimate the whole theoretically predicted nonlinear structure in experimentally recorded cells. -\subsection{Heterogeneity of P-units might influence nonlinearity}\notejg{Better title...} -The observed nonlinear effects might facilitate the detectability of faint signals during a three fish setting, the electrosensory cocktail party. They are, however, very specific with respect to the relation of stimulus frequencies and the P-unit baseline frequency. The EOD frequencies of the interacting fish would be drawn from the distributions of EOD frequencies found in male and female fish\citealp{Henninger2018, Henninger2020}\notejg{more refs on eod freq distributions?}. To be behaviorally relevant the faint signal detection would require reliable neuronal signalling irrespective of the individual EOD frequencies. -P-units, however are very heterogeneous in their baseline firing properties\citealp{Grewe2017, Hladnik2023}. The baseline firing rates vary in wide ranges (50--450\,Hz). This range covers substantial parts of the beat frequencies that may occur during animal interactions which is limited to frequencies below the Nyquist frequency (0 -- \feod/2)\citealp{Barayeu2023}. It is thus likely that there are P-units which match approximately to the specificities of the different encounters. +\subsection*{Heterogeneity of P-units might influence nonlinearity}\notejg{Better title...} +The observed nonlinear effects might facilitate the detectability of faint signals during a three fish setting, the electrosensory cocktail party. They are, however, very specific with respect to the relation of stimulus frequencies and the P-unit baseline frequency. The EOD frequencies of the interacting fish would be drawn from the distributions of EOD frequencies found in male and female fish\cite{Henninger2018, Henninger2020}\notejg{more refs on eod freq distributions?}. To be behaviorally relevant the faint signal detection would require reliable neuronal signalling irrespective of the individual EOD frequencies. +P-units, however are very heterogeneous in their baseline firing properties\cite{Grewe2017, Hladnik2023}. The baseline firing rates vary in wide ranges (50--450\,Hz). This range covers substantial parts of the beat frequencies that may occur during animal interactions which is limited to frequencies below the Nyquist frequency (0 -- \feod/2)\cite{Barayeu2023}. It is thus likely that there are P-units which match approximately to the specificities of the different encounters. -On the other hand, the nonlinearity was found only in low-CV P-units (with white noise stimulation). The CVs are also very heterogeneous (0.1--1.5\notejg{CVs aus Abbildung einfuegen}, \figref{fig:data_overview_mod}\panel{A}) in our sample. Only a small fraction of the P-units has a sufficiently low level of intrinsic noise and will exhibit nonlinear responses. The P-units project to the ELL\citealp{Krahe2014} and the integrating pyramidal cells in the different segments receive inputs in the range of 10s to 1000s of neurons\citealp{Maler2009a}. Since the receptive fields of the pyramidal neurons are patches of adjacent receptors on the fish surface \citealp{Bastian2002, Maler2009a, Haggard2023} and the input heterogeneity does not depend on the location of the receptor on the fish body \citealp{Hladnik2023} the pyramidal cell input population will be heterogeneous. Heterogeneity was shown to be generally advantageous for the encoding in this\citealp{Hladnik2023} and other systems\citealp{Padmanabhan2010, Beiran2018}. At the same time it contradicts the apparent need for a selective readout of low-CV cells to maintain information arising through nonlinear interactions. +On the other hand, the nonlinearity was found only in low-CV P-units (with white noise stimulation). The CVs are also very heterogeneous (0.1--1.5\notejg{CVs aus Abbildung einfuegen}, \figref{fig:data_overview_mod}\panel{A}) in our sample. Only a small fraction of the P-units has a sufficiently low level of intrinsic noise and will exhibit nonlinear responses. The P-units project to the ELL\cite{Krahe2014} and the integrating pyramidal cells in the different segments receive inputs in the range of 10s to 1000s of neurons\cite{Maler2009a}. Since the receptive fields of the pyramidal neurons are patches of adjacent receptors on the fish surface \cite{Bastian2002, Maler2009a, Haggard2023} and the input heterogeneity does not depend on the location of the receptor on the fish body \cite{Hladnik2023} the pyramidal cell input population will be heterogeneous. Heterogeneity was shown to be generally advantageous for the encoding in this\cite{Hladnik2023} and other systems\cite{Padmanabhan2010, Beiran2018}. At the same time it contradicts the apparent need for a selective readout of low-CV cells to maintain information arising through nonlinear interactions. A possible readout mechanism should be the topic of future studies that need to take into account that the nonlinearities are stronger in pure sine-wave stimulation and the fraction of cells that show it under naturalistic stimulation might be larger than expected from the distribution of CVs. -\subsection{Behavioral relevance of nonlinear interactions} -The behavioral relevance of the weak signal detection in P-units is evident from the courtship context observed in freely interacting animals\citealp{Henninger2018}. Outside courtship behavior, the encoding of secondary or social envelopes is a common need\citealp{Stamper2012Envelope}. In a previous study it was demonstrated that information about low-frequency secondary envelopes would not be present in P-units responses but would arise thorough nonlinear processing downstream in the ELL \citealp{Middleton2006,Middleton2007}. Based on our work we would predict that only a small subset of cells, with low CVs, should encode the social envelopes under white noise stimulation. An absence of low-CVs cells in the population analyzed in the previous studies could explain their conclusions. On the other hand, another study showed that P-units with strong nonlinearities, low firing rates and high CVs could indeed encode social envelopes\citealp{Savard2011}. These findings are in contrast to the previously mentioned work\citealp{Middleton2007} and, at first glance, also to our results. The missing link, that has not been considered in this work, might be the bursting of P-units, the repeated firing of spikes after one EOD period interleaved with longer intervals of quietness\citealp{Chacron2004}. Bursting was not explicitly addressed in the previous work, still the reported high CVs of the envelope encoding P-units indicate a higher rate of bursting\citealp{Savard2011}. How bursts influence the second-order susceptibility of will be addressed in following works (in preparation). - -Sinusoidal AMs are relevant in interactions with a few fish. We can understand the noise as the presence of many animals with individual EOD frequencies at the same time. Under noise stimulation, strong nonlinearities were demonstrated for weak stimuli but to decrease for stronger noise stimuli (\figrefb{fig:cells_suscept}). As long as the noise signal is weak, those fish are distant and the nonlinearity is maintained. An increasing stimulus amplitude would indicate that many fish are close to the receiver and a decrease of nonlinear effects can be observed. These findings imply that the nonlinear effects arising in the presence of three fish decline the more fish join. \lepto{} usually prefers small groups of fish\citealp{Stamper2010}. Thus, the described second-order susceptibility might still be behaviorally relevant under natural conditions. The decline of nonlinear effects when several fish are present might be an adaptive process reducing the number of frequencies represented in its primary sensory afferents to a minimum. Such representation would still leave room to create nonlinear effects at later processing steps in higher-order neurons. +\subsection*{Behavioral relevance of nonlinear interactions} +The behavioral relevance of the weak signal detection in P-units is evident from the courtship context observed in freely interacting animals\cite{Henninger2018}. Outside courtship behavior, the encoding of secondary or social envelopes is a common need\cite{Stamper2012Envelope}. In a previous study it was demonstrated that information about low-frequency secondary envelopes would not be present in P-units responses but would arise thorough nonlinear processing downstream in the ELL \cite{Middleton2006,Middleton2007}. Based on our work we would predict that only a small subset of cells, with low CVs, should encode the social envelopes under white noise stimulation. An absence of low-CVs cells in the population analyzed in the previous studies could explain their conclusions. On the other hand, another study showed that P-units with strong nonlinearities, low firing rates and high CVs could indeed encode social envelopes\cite{Savard2011}. These findings are in contrast to the previously mentioned work\cite{Middleton2007} and, at first glance, also to our results. The missing link, that has not been considered in this work, might be the bursting of P-units, the repeated firing of spikes after one EOD period interleaved with longer intervals of quietness\cite{Chacron2004}. Bursting was not explicitly addressed in the previous work, still the reported high CVs of the envelope encoding P-units indicate a higher rate of bursting\cite{Savard2011}. How bursts influence the second-order susceptibility of will be addressed in following works (in preparation). -The afferents of the passive electrosensory system, the ampullary cells, were found to exhibit much stronger nonlinearities than P-units (\figref{fig:data_overview_mod}). The adequate stimulus for this system is a direct stimulation not an amplitude modulation. In this sense the ampullary cells are closer to the LIF models used by Voroneko and colleagues\citealp{Voronenko2017} and we can thus expect that the same nonlinear mechanisms are at work here. For the ampullary system, the sinewave stimulation might is not as relevant as for the P-unit system. Ampullary cells encode low-frequency exogenous electric signals such as muscle potentials induced by prey movement\citealp{Kalmijn1974, Engelmann2010, Neiman2011fish}. The simultaneous muscle activity of a swarm of prey (such as \textit{Daphnia}) resembles Gaussian white noise\citealp{Neiman2011fish}, similar to the stimuli used here. Our results show some similarities with the analyses by Neiman and Russel\citealp{Neiman2011fish} who study the encoding in ampullary afferents in the paddlefish. There, the power spectrum of the spontaneous also shows a peak at the baseline frequency (internal oscillator) but also at the oscillation frequency measured at the epithelium and interactions of both. Most of the latter disappear in the case of external stimulation, though. Here we find only peaks at the baseline frequency of the neuron and its harmonics. There are interesting similarities and dissimilarities; stimulus encoding in the paddlefish as well as in the brown ghost is very linear for low frequencies and there are nonlinearities in both systems. Linear encoding in the paddlefish shows a gap in the spectrum at the frequency of the epitheliar oscillation, instead the nonlinear response is very pronounced there. In \lepto{}, the dominating frequency under baseline conditions is the baseline firing rate, and we do see increased nonlinearity in this frequency range. The baseline frequency, however, is outside the linear coding range\citealp{Grewe2017} while it is within the linear coding range in paddlefish\citealp{Neiman2011fish}. Interestingly, the nonlinear response in the paddlefish ampullaries increases with stimulus intensity while it disappears in our case (\figref{fig:data_overview_mod}\panel{D}). The population of ampullary cells is much more homogeneous with respect to the baseline rate (131$\pm$29\,Hz) and stimulus encoding properties than the P-units\citealp{Grewe2017}. This implies that, if the stimulus contains the appropriate frequency components that sum up to the baseline rate, there should be a nonlinear response at a frequency that is similar for the full population of ampullary cells (the baseline frequency) that is outside the linear coding range. Postsynaptic cells integrating ampullary input might be able to extract this nonlinear response from the input population. How such nonlinear effects in ampullary cells might influence prey detection should be addressed in further studies. +Sinusoidal AMs are relevant in interactions with a few fish. We can understand the noise as the presence of many animals with individual EOD frequencies at the same time. Under noise stimulation, strong nonlinearities were demonstrated for weak stimuli but to decrease for stronger noise stimuli (\figrefb{fig:cells_suscept}). As long as the noise signal is weak, those fish are distant and the nonlinearity is maintained. An increasing stimulus amplitude would indicate that many fish are close to the receiver and a decrease of nonlinear effects can be observed. These findings imply that the nonlinear effects arising in the presence of three fish decline the more fish join. \lepto{} usually prefers small groups of fish\cite{Stamper2010}. Thus, the described second-order susceptibility might still be behaviorally relevant under natural conditions. The decline of nonlinear effects when several fish are present might be an adaptive process reducing the number of frequencies represented in its primary sensory afferents to a minimum. Such representation would still leave room to create nonlinear effects at later processing steps in higher-order neurons. -\subsection{Conclusion} -We have demonstrated that there are pronounced nonlinear responses in the primary electrosensory afferences of the weakly electric fish \lepto{}, systems that are very often characterized using linear methods. The observed nonlinearities match the expectations from previous theoretical studies\citealp{Voronenko2017}. We can confirm that the theory applies also to systems that are encoding amplitude modulations of a carrier signal. Comparisons of P-units and ampullary cells showed that it is the level of intrinsic noise that determines how strongly nonlinear the system acts. Using the second-oder susceptibility estimated from the responses to white noise stimuli provides an easy way to determine the nonlinearity of the system under study. P-units share several features with mammalian -auditory nerve fibers and such nonlinear effects might also be expected in the auditory system during the encoding of amplitude modulations \citealp{Joris2004}. +The afferents of the passive electrosensory system, the ampullary cells, were found to exhibit much stronger nonlinearities than P-units (\figref{fig:data_overview_mod}). The adequate stimulus for this system is a direct stimulation not an amplitude modulation. In this sense the ampullary cells are closer to the LIF models used by Voroneko and colleagues\cite{Voronenko2017} and we can thus expect that the same nonlinear mechanisms are at work here. For the ampullary system, the sinewave stimulation might is not as relevant as for the P-unit system. Ampullary cells encode low-frequency exogenous electric signals such as muscle potentials induced by prey movement\cite{Kalmijn1974, Engelmann2010, Neiman2011fish}. The simultaneous muscle activity of a swarm of prey (such as \textit{Daphnia}) resembles Gaussian white noise\cite{Neiman2011fish}, similar to the stimuli used here. Our results show some similarities with the analyses by Neiman and Russel\cite{Neiman2011fish} who study the encoding in ampullary afferents in the paddlefish. There, the power spectrum of the spontaneous also shows a peak at the baseline frequency (internal oscillator) but also at the oscillation frequency measured at the epithelium and interactions of both. Most of the latter disappear in the case of external stimulation, though. Here we find only peaks at the baseline frequency of the neuron and its harmonics. There are interesting similarities and dissimilarities; stimulus encoding in the paddlefish as well as in the brown ghost is very linear for low frequencies and there are nonlinearities in both systems. Linear encoding in the paddlefish shows a gap in the spectrum at the frequency of the epitheliar oscillation, instead the nonlinear response is very pronounced there. In \lepto{}, the dominating frequency under baseline conditions is the baseline firing rate, and we do see increased nonlinearity in this frequency range. The baseline frequency, however, is outside the linear coding range\cite{Grewe2017} while it is within the linear coding range in paddlefish\cite{Neiman2011fish}. Interestingly, the nonlinear response in the paddlefish ampullaries increases with stimulus intensity while it disappears in our case (\figref{fig:data_overview_mod}\panel{D}). The population of ampullary cells is much more homogeneous with respect to the baseline rate (131$\pm$29\,Hz) and stimulus encoding properties than the P-units\cite{Grewe2017}. This implies that, if the stimulus contains the appropriate frequency components that sum up to the baseline rate, there should be a nonlinear response at a frequency that is similar for the full population of ampullary cells (the baseline frequency) that is outside the linear coding range. Postsynaptic cells integrating ampullary input might be able to extract this nonlinear response from the input population. How such nonlinear effects in ampullary cells might influence prey detection should be addressed in further studies. +\subsection*{Conclusion} +We have demonstrated that there are pronounced nonlinear responses in the primary electrosensory afferences of the weakly electric fish \lepto{}, systems that are very often characterized using linear methods. The observed nonlinearities match the expectations from previous theoretical studies\cite{Voronenko2017}. We can confirm that the theory applies also to systems that are encoding amplitude modulations of a carrier signal. Comparisons of P-units and ampullary cells showed that it is the level of intrinsic noise that determines how strongly nonlinear the system acts. Using the second-oder susceptibility estimated from the responses to white noise stimuli provides an easy way to determine the nonlinearity of the system under study. P-units share several features with mammalian +auditory nerve fibers and such nonlinear effects might also be expected in the auditory system during the encoding of amplitude modulations \cite{Joris2004}. +\section*{Methods} -\section{Methods} - -\subsection{Experimental subjects and procedures} +\subsection*{Experimental subjects and procedures} Within this project we re-evaluated datasets that were recorded between 2010 and 2023 at the Ludwig Maximilian University (LMU) M\"unchen and the Eberhard-Karls University T\"ubingen. All experimental protocols complied with national and European law and were approved by the respective Ethics Committees of the Ludwig-Maximilians Universität München (permit no. 55.2-1-54-2531-135-09) and the Eberhard-Karls Unversität Tübingen (permit no. ZP 1/13 and ZP 1/16). The final sample consisted of 221 P-units and 45 ampullary electroreceptor afferents recorded in 71 weakly electric fish of the species \lepto{}. The original electrophysiological recordings were performed on male and female weakly electric fish of the species \lepto{} that were obtained from a commercial supplier for tropical fish (Aquarium Glaser GmbH, Rodgau, @@ -551,29 +632,29 @@ Germany). The fish were kept in tanks with a water temperature of $25\pm1\,^\cir Before surgery, the animals were deeply anesthetized via bath application with a solution of MS222 (120\,mg/l, PharmaQ, Fordingbridge, UK) buffered with Sodium Bicarbonate (120\,mg/l). The posterior anterior lateral line nerve (pALLN) was exposed by making a small cut into the skin covering the nerve. The cut was placed dorsal of the operculum just before the nerve descends towards the anterior lateral line ganglion (ALLNG). Those parts of the skin that were to be cut were locally anesthetized by cutaneous application of liquid lidocaine hydrochloride (20\,mg/ml, bela-pharm GmbH). During the surgery water supply was ensured by a mouthpiece to maintain anesthesia with a solution of MS222 (100\,mg/l) buffered with Sodium Bicarbonate (100\,mg/l). After surgery fish were immobilized by intramuscular injection of from 25\,$\micro$l to 50\,$\micro$l of tubocurarine (5\,mg/ml dissolved in fish saline; Sigma-Aldrich). Respiration was then switched to normal tank water and the fish was transferred to the experimental tank. -\subsection{Experimental setup} -For the recordings fish were positioned centrally in the experimental tank, with the major parts of their body submerged into the water. Those body parts that were above the water surface were covered with paper tissue to avoid drying of the skin. Local analgesia was refreshed in intervals of two hours by cutaneous reapplication of Lidocaine (2\,\%; bela-pharm, Vechta, Germany) around the surgical wounds. Electrodes (borosilicate; 1.5\,mm outer diameter; GB150F-8P; Science Products, Hofheim, Germany) were pulled to a resistance of 50--100\,\mega\ohm{} (model P-97; Sutter Instrument, Novato, CA) and filled with 1\,M KCl solution. Electrodes were fixed in a microdrive (Luigs-Neumann, Ratingen, Germany) and lowered into the nerve (\figrefb{Setup}, blue triangle). Recordings of electroreceptor afferents were amplified and lowpass filtered at 10\,kHz (SEC-05, npi-electronics, Tamm, Germany, operated in bridge mode). All signals, neuronal recordings, recorded EOD and the generated stimulus, were digitized with sampling rates of 20 or 40\,kHz (PCI-6229, National Instruments, Austin, TX). RELACS (\url{www.relacs.net}) running on a Linux computer was used for online spike and EOD detection, stimulus generation, and calibration. Recorded data was then stored on the hard drive for offline analysis. +\subsection*{Experimental setup} +For the recordings fish were positioned centrally in the experimental tank, with the major parts of their body submerged into the water. Those body parts that were above the water surface were covered with paper tissue to avoid drying of the skin. Local analgesia was refreshed in intervals of two hours by cutaneous reapplication of Lidocaine (2\,\%; bela-pharm, Vechta, Germany) around the surgical wounds. Electrodes (borosilicate; 1.5\,mm outer diameter; GB150F-8P; Science Products, Hofheim, Germany) were pulled to a resistance of 50--100\,\mega\ohm{} (model P-97; Sutter Instrument, Novato, CA) and filled with 1\,M KCl solution. Electrodes were fixed in a microdrive (Luigs-Neumann, Ratingen, Germany) and lowered into the nerve (\figrefb{fig:setup}, blue triangle). Recordings of electroreceptor afferents were amplified and lowpass filtered at 10\,kHz (SEC-05, npi-electronics, Tamm, Germany, operated in bridge mode). All signals, neuronal recordings, recorded EOD and the generated stimulus, were digitized with sampling rates of 20 or 40\,kHz (PCI-6229, National Instruments, Austin, TX). RELACS (\url{www.relacs.net}) running on a Linux computer was used for online spike and EOD detection, stimulus generation, and calibration. Recorded data was then stored on the hard drive for offline analysis. -\subsection{Identification of P-units and ampullary cells} -The neurons were classified into cell types during the recording by the experimenter. P-units were classified based on baseline firing rates of 50--450\,Hz and a clear phase-locking to the EOD and their responses to amplitude modulations of their own EOD\citealp{Grewe2017, Hladnik2023}. Ampullary cells were classified based on firing rates of 80--200\,Hz absent phase-locking to the EOD and responses to low-frequency sinusoidal stimuli\citealp{Grewe2017}. We here selected only those cells of which the neuron's baseline activity as well as the responses to frozen noise stimuli were recorded. +\subsection*{Identification of P-units and ampullary cells} +The neurons were classified into cell types during the recording by the experimenter. P-units were classified based on baseline firing rates of 50--450\,Hz and a clear phase-locking to the EOD and their responses to amplitude modulations of their own EOD\cite{Grewe2017, Hladnik2023}. Ampullary cells were classified based on firing rates of 80--200\,Hz absent phase-locking to the EOD and responses to low-frequency sinusoidal stimuli\cite{Grewe2017}. We here selected only those cells of which the neuron's baseline activity as well as the responses to frozen noise stimuli were recorded. -\subsection{Electric field recordings} - The electric field of the fish was recorded in two ways: 1. we measured the so-called \textit{global EOD} with two vertical carbon rods ($11\,\centi\meter$ long, 8\,mm diameter) in a head-tail configuration (\figrefb{Setup}, green bars). The electrodes were placed isopotential to the stimulus. This signal was differentially amplified with a factor between 100 and 500 (depending on the recorded animal) and band-pass filtered (3 to 1500\,Hz pass-band, DPA2-FX; npi electronics, Tamm, Germany). 2. The so-called \textit{local EOD} was measured with 1\,cm-spaced silver wires located next to the left gill of the fish and orthogonal to the fish's longitudinal body axis (amplification 100 to 500 times, band-pass filtered with 3 to 1\,500\,Hz pass-band, DPA2-FX; npi-electronics, Tamm, Germany, \figrefb{Setup}, red markers). This local measurement recorded the combination of the fish's own field and the applied stimulus and thus serves as a proxy of the transdermal potential that drives the electroreceptors. +\subsection*{Electric field recordings} + The electric field of the fish was recorded in two ways: 1. we measured the so-called \textit{global EOD} with two vertical carbon rods ($11\,\centi\meter$ long, 8\,mm diameter) in a head-tail configuration (\figrefb{fig:setup}, green bars). The electrodes were placed isopotential to the stimulus. This signal was differentially amplified with a factor between 100 and 500 (depending on the recorded animal) and band-pass filtered (3 to 1500\,Hz pass-band, DPA2-FX; npi electronics, Tamm, Germany). 2. The so-called \textit{local EOD} was measured with 1\,cm-spaced silver wires located next to the left gill of the fish and orthogonal to the fish's longitudinal body axis (amplification 100 to 500 times, band-pass filtered with 3 to 1\,500\,Hz pass-band, DPA2-FX; npi-electronics, Tamm, Germany, \figrefb{fig:setup}, red markers). This local measurement recorded the combination of the fish's own field and the applied stimulus and thus serves as a proxy of the transdermal potential that drives the electroreceptors. -\subsection{Stimulation} -The stimulus was isolated from the ground (ISO-02V, npi-electronics, Tamm, Germany) and delivered via two horizontal carbon rods (30 cm length, 8 mm diameter) located $15\,\centi\meter$ laterally to the fish (\figrefb{Setup}, gray bars). The stimulus was calibrated with respect to the local EOD. +\subsection*{Stimulation} +The stimulus was isolated from the ground (ISO-02V, npi-electronics, Tamm, Germany) and delivered via two horizontal carbon rods (30 cm length, 8 mm diameter) located $15\,\centi\meter$ laterally to the fish (\figrefb{fig:setup}, gray bars). The stimulus was calibrated with respect to the local EOD. -\begin{figure*}[h!]%(\subfigrefb{beat_amplitudes}{B}). +\begin{figure*}[!ht]%(\subfigrefb{beat_amplitudes}{B}). \includegraphics[width=\columnwidth]{Settup} - \caption{\label{Setup} Electrophysiolocical recording setup. The fish, depicted as a black scheme and surrounded by isopotential lines, was positioned in the center of the tank. Blue triangle -- electrophysiological recordings were conducted in the posterior anterior lateral line nerve (pALLN). Gray horizontal bars -- electrodes for the stimulation. Green vertical bars -- electrodes to measure the \textit{global EOD} placed isopotential to the stimulus, i.e. recording fish's unperturbed EOD. Red dots -- electrodes to measure the \textit{local EOD} picking up the combination of fish's EOD and the stimulus. The local EOD was measured with a distance of 1 \,cm between the electrodes. All measured signals were amplified, filtered, and stored for offline analysis.} + \caption{\label{fig:setup} Electrophysiolocical recording setup. The fish, depicted as a black scheme and surrounded by isopotential lines, was positioned in the center of the tank. Blue triangle -- electrophysiological recordings were conducted in the posterior anterior lateral line nerve (pALLN). Gray horizontal bars -- electrodes for the stimulation. Green vertical bars -- electrodes to measure the \textit{global EOD} placed isopotential to the stimulus, i.e. recording fish's unperturbed EOD. Red dots -- electrodes to measure the \textit{local EOD} picking up the combination of fish's EOD and the stimulus. The local EOD was measured with a distance of 1 \,cm between the electrodes. All measured signals were amplified, filtered, and stored for offline analysis.} \end{figure*} -\subsection{White noise stimulation}\label{rammethods} +\subsection*{White noise stimulation}\label{rammethods} The fish were stimulated with band-limited white noise stimuli with a cut-off frequency of 150, 300 or 400\,Hz. The stimulus intensity is given as the contrast, i.e. the standard deviation of the white noise stimulus in relation to the fish's EOD amplitude. The contrast varied between 1 and 20\,$\%$. Only cell recordings with at least 10\,s of white noise stimulation were included for the analysis. When ampullary cells were recorded, the white noise was directly applied as the stimulus. To create random amplitude modulations (RAM) for P-unit recordings, the EOD of the fish was multiplied with the desired random amplitude modulation profile (MXS-01M; npi electronics). % and between 2.5 and 40\,$\%$ for \eigen -\subsection{Data analysis} Data analysis was performed with Python 3 using the packages matplotlib\cite{Hunter2007}, numpy\cite{Walt2011}, scipy\cite{scipy2020}, pandas\cite{Mckinney2010}, nixio\cite{Stoewer2014}, and thunderfish (\url{https://github.com/bendalab/thunderfish}). +\subsection*{Data analysis} Data analysis was performed with Python 3 using the packages matplotlib\cite{Hunter2007}, numpy\cite{Walt2011}, scipy\cite{scipy2020}, pandas\cite{Mckinney2010}, nixio\cite{Stoewer2014}, and thunderfish (\url{https://github.com/bendalab/thunderfish}). %sklearn\cite{scikitlearn2011}, @@ -643,7 +724,7 @@ The second-order susceptibility was calculated by dividing the higher-order cros % \end{split} % \end{equation} % \notejg{Wofuer genau brauchen wir equation 9?} -The absolute value of a second-order susceptibility matrix is visualized in \figrefb{model_full}. There the upper right and the lower left quadrants characterize the nonlinearity in the response $r(t)$ at the sum frequency of the two input frequencies. The lower right and upper left quadrants characterize the nonlinearity in the response $r(t)$ at the difference of the input frequencies. +The absolute value of a second-order susceptibility matrix is visualized in \figrefb{fig:model_full}. There the upper right and the lower left quadrants characterize the nonlinearity in the response $r(t)$ at the sum frequency of the two input frequencies. The lower right and upper left quadrants characterize the nonlinearity in the response $r(t)$ at the difference of the input frequencies. \paragraph{Nonlinearity index}\label{projected_method} We expect to see nonlinear susceptibility when $\omega_1 + \omega_2 = \fbase{}$. To characterize this we calculated the nonlinearity index (NLI) as @@ -655,9 +736,9 @@ For this index, the second-order susceptibility matrix was projected onto the di If the same frozen noise was recorded several times in a cell, each noise repetition resulted in a separate second-order susceptibility matrix. The mean of the corresponding \nli{} values was used for the population statistics in \figref{fig:data_overview_mod}. -\subsection{Leaky integrate-and-fire models}\label{lifmethods} +\subsection*{Leaky integrate-and-fire models}\label{lifmethods} -Leaky integrate-and-fire (LIF) models with a carrier were constructed to reproduce the specific firing properties of P-units \citealp{Chacron2001,Sinz2020}. The sole driving input into the P-unit model during baseline, i.e. when no external stimulus was given, is the fish's own EOD modeled as a cosine wave +Leaky integrate-and-fire (LIF) models with a carrier were constructed to reproduce the specific firing properties of P-units \cite{Chacron2001,Sinz2020}. The sole driving input into the P-unit model during baseline, i.e. when no external stimulus was given, is the fish's own EOD modeled as a cosine wave \begin{equation} \label{eq:eod} x(t) = x_{EOD}(t) = \cos(2\pi f_{EOD} t) @@ -703,27 +784,27 @@ Whenever the membrane voltage $V_m(t)$ crossed the spiking threshold $\theta=1$ % \label{eifnl} % f(V_m)= \Delta_V \text{e}^{\frac{V_m-1}{\Delta_V}} % \end{equation} -% \citealp{Fourcaud-Trocme2003}, where $\Delta_V$ was varied from 0.001 to 0.1. +% \cite{Fourcaud-Trocme2003}, where $\Delta_V$ was varied from 0.001 to 0.1. %, \figrefb{eif} -\begin{figure*}[hb!] +\begin{figure*}[!ht] \includegraphics[width=\columnwidth]{flowchart} \caption{\label{flowchart} - Flowchart of a LIF P-unit model with EOD carrier. The main steps of the model are depicted in the left model column (Model cell identifier 2012-07-03-ak, see table~\ref{modelparams} for model parameters). The three other columns show the relevant signals in three different settings. (i) the baseline situation, no external stimulus, only the animal's self-generated EOD (i.e. the carrier) is present (ii) RAM stimulation, the carrier is amplitude modulated with a weak (2\% contrast) stimulus, (iii) Noise split condition in which 90\% of the internal noise is used as a driving stimulus scaled with the correction factor $\rho$ (see text). Note: that the firing rate and the CV of the ISI distribution is the same in this and the baseline condition. \textbf{A} Thresholding: a simple linear threshold was applied to the EOD carrier, \eqnref{eq:eod}. The red line on top depicts the amplitude modulation (AM). \textbf{B} Dendritic low-pass filtering attenuates the carrier. \textbf{C} A Gaussian noise is added to the signal in B. Note the reduced internal noise amplitude in the noise split (iii) condition. \textbf{D} Spiking output of the LIF model in response to the addition of B and C. \textbf{E} Power spectra of the LIF neuron's spiking activity. Under the baseline condition (\panel[i]{D}) there are several peaks, from left to right, at the baseline firing rate $\fbase{}$, $f_{EOD} - \fbase{}$ $f_{EOD}$, and $f_{EOD} + \fbase{}$. In the stimulus driven regime (\panel[ii]{D}), there is only a peak at $f_{eod}$, while under the noise split condition (\panel[iii]{D}) again all peaks are present.} + Flowchart of a LIF P-unit model with EOD carrier. The main steps of the model are depicted in the left model column (Model cell identifier 2012-07-03-ak, see table~\ref{modelparams} for model parameters). The three other columns show the relevant signals in three different settings. (i) the baseline situation, no external stimulus, only the animal's self-generated EOD (i.e. the carrier) is present (ii) RAM stimulation, the carrier is amplitude modulated with a weak (2\% contrast) stimulus, (iii) Noise split condition in which 90\% of the internal noise is used as a driving stimulus scaled with the correction factor $\rho$ (see text). Note: that the firing rate and the CV of the ISI distribution is the same in this and the baseline condition. \figitem{A} Thresholding: a simple linear threshold was applied to the EOD carrier, \Eqnref{eq:eod}. The red line on top depicts the amplitude modulation (AM). \figitem{B} Dendritic low-pass filtering attenuates the carrier. \figitem{C} A Gaussian noise is added to the signal in B. Note the reduced internal noise amplitude in the noise split (iii) condition. \figitem{D} Spiking output of the LIF model in response to the addition of B and C. \figitem{E} Power spectra of the LIF neuron's spiking activity. Under the baseline condition (\panel[i]{D}) there are several peaks, from left to right, at the baseline firing rate $\fbase{}$, $f_{EOD} - \fbase{}$ $f_{EOD}$, and $f_{EOD} + \fbase{}$. In the stimulus driven regime (\panel[ii]{D}), there is only a peak at $f_{eod}$, while under the noise split condition (\panel[iii]{D}) again all peaks are present.} \end{figure*} -\subsection{Numerical implementation} -The model's ODEs were integrated by the Euler forward method with a time-step of $\Delta t = 0.05$\,ms. The intrinsic noise $\xi(t)$ (\eqnref{eq:LIF}, \subfigrefb{flowchart}{C}) was added by drawing a random number from a normal distribution $\mathcal{N}(0,\,1)$ with zero mean and standard deviation of one in each time step $i$. This number was multiplied with $\sqrt{2D}$ and divided by $\sqrt{\Delta t}$: +\subsection*{Numerical implementation} +The model's ODEs were integrated by the Euler forward method with a time-step of $\Delta t = 0.05$\,ms. The intrinsic noise $\xi(t)$ (\Eqnref{eq:LIF}, \subfigrefb{flowchart}{C}) was added by drawing a random number from a normal distribution $\mathcal{N}(0,\,1)$ with zero mean and standard deviation of one in each time step $i$. This number was multiplied with $\sqrt{2D}$ and divided by $\sqrt{\Delta t}$: \begin{equation} \label{eq:LIFintegration} V_{m_{i+1}} = V_{m_i} + \left(-V_{m_i} + \mu + \alpha V_{d_i} - A_i + \sqrt{\frac{2D}{\Delta t}}\mathcal{N}(0,\,1)_i\right) \frac{\Delta t}{\tau_m} \end{equation} -\subsection{Model parameters}\label{paramtext} +\subsection*{Model parameters}\label{paramtext} The eight free parameters of the P-unit model $\beta$, $\tau_m$, $\mu$, $D$, $\tau_A$, $\Delta_A$, $\tau_d$, and $t_{ref}$, were fitted to both the baseline activity (baseline firing rate, CV of ISIs, serial correlation of ISIs at lag one, and vector strength of spike coupling to EOD) and the responses to step-like increases and decreases in EOD amplitude (onset-state and steady-state responses, effective adaptation time constant). For each simulation, the start parameters $A$, $V_{d}$ and $V_{m}$ were drawn from a random starting value distribution, estimated from a 100\,s baseline simulation after an initial 100\,s of simulation that was discarded as a transient. -\subsection{Stimuli for the model} -The model neurons were driven with similar stimuli as the real neurons in the original experiments. To mimic the interaction with one or two foreign animals the receiving fish's EOD (\eqnref{eq:eod}) was normalized to an amplitude of one and the respective EODs of a second or third fish were added.%\eqnref{ eq.\,\ref{eq:eod} +\subsection*{Stimuli for the model} +The model neurons were driven with similar stimuli as the real neurons in the original experiments. To mimic the interaction with one or two foreign animals the receiving fish's EOD (\Eqnref{eq:eod}) was normalized to an amplitude of one and the respective EODs of a second or third fish were added.%\Eqnref{ eq.\,\ref{eq:eod} The random amplitude modulation (RAM) input to the model was created by drawing random amplitude and phases from Gaussian distributions for each frequency component in the range 0--300 Hz. An inverse Fourier transform was applied to get the final amplitude RAM time-course. The input to the model was then \begin{equation} @@ -735,10 +816,10 @@ From each simulation run, the first second was discarded and the analysis was ba % \subsection{Second-order susceptibility analysis of the model} % %\subsubsection{Model second-order nonlinearity} -% The second-order susceptibility in the model was calculated with \eqnref{eq:susceptibility}, resulting in matrices as in \figrefb{model_and_data} and \figrefb{model_full}. For this, the model neuron was presented the input $x(t)$ for 2\,s, with the first second being dismissed as the transient. The second-order susceptibility calculation was performed on the last second, resulting in a frequency resolution of 1\,Hz. +% The second-order susceptibility in the model was calculated with \Eqnref{eq:susceptibility}, resulting in matrices as in \figrefb{model_and_data} and \figrefb{fig:model_full}. For this, the model neuron was presented the input $x(t)$ for 2\,s, with the first second being dismissed as the transient. The second-order susceptibility calculation was performed on the last second, resulting in a frequency resolution of 1\,Hz. -\subsection{Model noise split into a noise and a stimulus component}\label{intrinsicsplit_methods} -According to the Novikov-Furutsu Theorem \citealp{Novikov1965, Furutsu1963} the total noise of a LIF model ($\xi$) can be split up into several independent noise processes with the same correlation function. Here we split the internal noise into two parts: (i) One part is treated as a driving input signal $\xi_{signal} = \sqrt{\rho \, 2D \,c_{signal}} \cdot \xi(t) $ and used to calculate the cross-spectra in \eqnref{eq:crosshigh} and (ii) the remaining noise $\xi_{noise} = \sqrt{2D \, c_{noise}} \cdot \xi(t)$ that is treated as pure noise. In this way the effective signal-to-noise ratio can be increased while maintaining the total noise in the system. $\rho$ a scaling factor that compensates (see below) for the signal transformations the amplitude modulation stimulus undergoes in the model, i.e. the threshold and the dendritic lowpass. In our case the model has a carrier (the fish's self-generated EOD) and we thus want to drive the model with an amplitude modulation stimulus +\subsection*{Model noise split into a noise and a stimulus component}\label{intrinsicsplit_methods} +According to the Novikov-Furutsu Theorem \cite{Novikov1965, Furutsu1963} the total noise of a LIF model ($\xi$) can be split up into several independent noise processes with the same correlation function. Here we split the internal noise into two parts: (i) One part is treated as a driving input signal $\xi_{signal} = \sqrt{\rho \, 2D \,c_{signal}} \cdot \xi(t) $ and used to calculate the cross-spectra in \Eqnref{eq:crosshigh} and (ii) the remaining noise $\xi_{noise} = \sqrt{2D \, c_{noise}} \cdot \xi(t)$ that is treated as pure noise. In this way the effective signal-to-noise ratio can be increased while maintaining the total noise in the system. $\rho$ a scaling factor that compensates (see below) for the signal transformations the amplitude modulation stimulus undergoes in the model, i.e. the threshold and the dendritic lowpass. In our case the model has a carrier (the fish's self-generated EOD) and we thus want to drive the model with an amplitude modulation stimulus %(1-c_{signal})\cdot\xi$c_{noise} = 1-c_{signal}$ %c_{signal} \cdot \xi @@ -766,12 +847,12 @@ According to the Novikov-Furutsu Theorem \citealp{Novikov1965, Furutsu1963} the -A big portion of the total noise was assigned to the signal component ($c_{signal} = 0.9$) and the remaining part to the noise component ($c_{noise} = 0.1$, \subfigrefb{flowchart}\,\panel[iii]{C}). For the application of the Novikov-Furutsu Theorem \citealp{Novikov1965, Furutsu1963} it is critical that both components add up to the initial 100\,$\%$ of the total noise. This is easily achieved in a model without a carrier if the condition $c_{signal}+c_{noise}=1$ is satisfied. The situation here is more complicated. To compensate for the transformations the signal undergoes before it enters the LIF core, $\xi_{signal}(t)$ was scaled up by the factor $\rho$ (\eqnref{eq:ram_split}, red in \subfigrefb{flowchart}\,\panel[iii]{A}). $\rho$ was found by bisecting the space of possible $\rho$ scaling factors by minimizing the error between the baseline CV (only carrier) and the CV during stimulation with noise split. +A big portion of the total noise was assigned to the signal component ($c_{signal} = 0.9$) and the remaining part to the noise component ($c_{noise} = 0.1$, \subfigrefb{flowchart}\,\panel[iii]{C}). For the application of the Novikov-Furutsu Theorem \cite{Novikov1965, Furutsu1963} it is critical that both components add up to the initial 100\,$\%$ of the total noise. This is easily achieved in a model without a carrier if the condition $c_{signal}+c_{noise}=1$ is satisfied. The situation here is more complicated. To compensate for the transformations the signal undergoes before it enters the LIF core, $\xi_{signal}(t)$ was scaled up by the factor $\rho$ (\Eqnref{eq:ram_split}, red in \subfigrefb{flowchart}\,\panel[iii]{A}). $\rho$ was found by bisecting the space of possible $\rho$ scaling factors by minimizing the error between the baseline CV (only carrier) and the CV during stimulation with noise split. % See section \ref{lifmethods} for model and parameter description. \begin{table*}[hp!] - \caption{\label{modelparams} Model parameters of LIF models, fitted to 2 electrophysiologically recorded P-units (\citealp{Ott2020}).} + \caption{\label{modelparams} Model parameters of LIF models, fitted to 2 electrophysiologically recorded P-units (\cite{Ott2020}).} \begin{center} \begin{tabular}{lrrrrrrrr} \hline @@ -784,84 +865,32 @@ A big portion of the total noise was assigned to the signal component ($c_{signa \end{table*}% 2013-01-08-aa % 2012-07-03-ak - - +% Either type in your references using +% \begin{thebibliography}{} +% \bibitem{} +% Text +% \end{thebibliography} +% +% or +% +% Compile your BiBTeX database using our plos2015.bst +% style file and paste the contents of your .bbl file +% here. See http://journals.plos.org/plosone/s/latex for +% step-by-step instructions. +% +\bibliography{references} \newpage - - - - - - -%Recording at the frequency combinations \bcsum{} and \bcdiff{} \fbasecorr{}at the burst-corrected firing rate - - - - - - - - - - - - - - - - - -\appendix -\setcounter{secnumdepth}{2} -\section{Appendix} - -\begin{figure*}[hp]%hp! -\includegraphics{cells_suscept_high_CV} - \caption{\label{fig:cells_suscept_high_CV} Response of experimentally measured P-units to RAM stimuli. Light purple -- low RAM contrast. Dark purple -- high RAM contrast. Noisy high-CV P-Unit. \figitem{A} Interspike intervals (ISI) distribution during baseline. \figitem{B} Baseline power spectrum of the firing rate. \figitem{C} Top: EOD carrier (gray) with RAM (red). Middle: Spike trains in response to a low RAM contrast. Bottom: Spike trains in response to a high RAM contrast. \figitem{D} First-order susceptibility (see \eqnref{linearencoding_methods}). \figitem{E} Absolute value of the second-order susceptibility, \eqnref{eq:susceptibility}, for the low RAM contrast. +\section*{Supporting information} +\paragraph*{S1 Second-order susceptibility of high-CV P-unit} +\label{S1:highcvpunit} +\begin{figure*}[!ht] +\includegraphics[width=\columnwidth]{cells_suscept_high_CV} + \caption{\label{fig:cells_suscept_high_CV} Response of experimentally measured P-units to RAM stimuli. Light purple -- low RAM contrast. Dark purple -- high RAM contrast. Noisy high-CV P-Unit. \figitem{A} Interspike intervals (ISI) distribution during baseline. \figitem{B} Baseline power spectrum of the firing rate. \figitem{C} Top: EOD carrier (gray) with RAM (red). Middle: Spike trains in response to a low RAM contrast. Bottom: Spike trains in response to a high RAM contrast. \figitem{D} First-order susceptibility (see \Eqnref{linearencoding_methods}). \figitem{E} Absolute value of the second-order susceptibility, \Eqnref{eq:susceptibility}, for the low RAM contrast. Pink lines -- edges of the structure when \fone, \ftwo{} or \fsum{} are equal to \fbase{}. \figitem{F} Absolute value of the second-order susceptibility for the higher RAM contrast. \figitem{G} Projected diagonals, calculated as the mean of the anti-diagonals of the matrices in \panel{E--F}. Gray dots: \fbase{}. Dashed lines: Medians of the projected diagonals. } \end{figure*} -\bibliographystyle{iscience} -%\bibliographystyle{apalike}%alpha}%}%alpha}%apalike} -%\bibliographystyle{elsarticle-num-names}%elsarticle-num-names} -%\ExecuteBibliographyOptions{sorting=nty} -%\bibliographystyle{authordate2} -\bibliography{journalsabbrv,references} -%\begin{thebibliography}{00} - -%\bibliographystyle{elsarticle-num-names} -%\bibliography{journalsabbrv,references} -%% \bibitem{label} -%% Text of bibliographic item - - - - - - - - - - - - -\clearpage -\renewcommand{\thesubsection}{\arabic{section}.\arabic{subsection}} - -\renewcommand{\thefigure}{S\arabic{figure}} -\setcounter{figure}{0} - -\renewcommand{\theequation}{S\arabic{equation}} -\setcounter{equation}{0} - -\newpage - - -\end{document} -\endinput -%% -%% End of file `elsarticle-template-num.tex'. +\end{document} \ No newline at end of file