<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD with MathML3 v1.3 20210610//EN"  "JATS-archivearticle1-3-mathml3.dtd"><article xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.3"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic" pub-type="epub">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">91722</article-id><article-id pub-id-type="doi">10.7554/eLife.91722</article-id><article-id pub-id-type="doi" specific-use="version">10.7554/eLife.91722.4</article-id><article-version article-version-type="publication-state">version of record</article-version><article-categories><subj-group subj-group-type="display-channel"><subject>Research Article</subject></subj-group><subj-group subj-group-type="heading"><subject>Neuroscience</subject></subj-group></article-categories><title-group><article-title>Brain state and cortical layer-specific mechanisms underlying perception at threshold</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Morton</surname><given-names>Mitchell P</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0001-9185-1633</contrib-id><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="other" rid="fund5"/><xref ref-type="other" rid="fund6"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author"><name><surname>Denagamage</surname><given-names>Sachira</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author"><name><surname>Blume</surname><given-names>Isabel J</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author"><name><surname>Reynolds</surname><given-names>John H</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0001-6988-4607</contrib-id><xref ref-type="aff" rid="aff3">3</xref><xref ref-type="other" rid="fund1"/><xref ref-type="other" rid="fund7"/><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author"><name><surname>Jadi</surname><given-names>Monika P</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0003-1092-5026</contrib-id><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="aff" rid="aff4">4</xref><xref ref-type="aff" rid="aff5">5</xref><xref ref-type="other" rid="fund3"/><xref ref-type="other" rid="fund4"/><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes"><name><surname>Nandy</surname><given-names>Anirvan S</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0002-4225-5349</contrib-id><email>anirvan.nandy@yale.edu</email><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="aff" rid="aff5">5</xref><xref ref-type="aff" rid="aff6">6</xref><xref ref-type="aff" rid="aff7">7</xref><xref ref-type="other" rid="fund1"/><xref ref-type="other" rid="fund2"/><xref ref-type="other" rid="fund8"/><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><label>1</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Department of Neuroscience, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff><aff id="aff2"><label>2</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Interdepartmental Neuroscience Program, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff><aff id="aff3"><label>3</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03xez1567</institution-id><institution>Systems Neurobiology Laboratories, The Salk Institute for Biological Studies</institution></institution-wrap><addr-line><named-content content-type="city">La Jolla</named-content></addr-line><country>United States</country></aff><aff id="aff4"><label>4</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Department of Psychiatry, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff><aff id="aff5"><label>5</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Wu Tsai Institute, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff><aff id="aff6"><label>6</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Department of Psychology, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff><aff id="aff7"><label>7</label><institution-wrap><institution-id institution-id-type="ror">https://ror.org/03v76x132</institution-id><institution>Kavli Institute for Neuroscience, Yale University</institution></institution-wrap><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Gold</surname><given-names>Joshua I</given-names></name><role>Reviewing Editor</role><aff><institution-wrap><institution-id institution-id-type="ror">https://ror.org/00b30xv10</institution-id><institution>University of Pennsylvania</institution></institution-wrap><country>United States</country></aff></contrib><contrib contrib-type="senior_editor"><name><surname>Gold</surname><given-names>Joshua I</given-names></name><role>Senior Editor</role><aff><institution-wrap><institution-id institution-id-type="ror">https://ror.org/00b30xv10</institution-id><institution>University of Pennsylvania</institution></institution-wrap><country>United States</country></aff></contrib></contrib-group><pub-date publication-format="electronic" date-type="publication"><day>18</day><month>11</month><year>2024</year></pub-date><volume>12</volume><elocation-id>RP91722</elocation-id><history><date date-type="sent-for-review" iso-8601-date="2023-10-10"><day>10</day><month>10</month><year>2023</year></date></history><pub-history><event><event-desc>This manuscript was published as a preprint.</event-desc><date date-type="preprint" iso-8601-date="2023-10-10"><day>10</day><month>10</month><year>2023</year></date><self-uri content-type="preprint" xlink:href="https://doi.org/10.1101/2021.11.16.468866"/></event><event><event-desc>This manuscript was published as a reviewed preprint.</event-desc><date date-type="reviewed-preprint" iso-8601-date="2023-12-11"><day>11</day><month>12</month><year>2023</year></date><self-uri content-type="reviewed-preprint" xlink:href="https://doi.org/10.7554/eLife.91722.1"/></event><event><event-desc>The reviewed preprint was revised.</event-desc><date date-type="reviewed-preprint" iso-8601-date="2024-08-14"><day>14</day><month>08</month><year>2024</year></date><self-uri content-type="reviewed-preprint" xlink:href="https://doi.org/10.7554/eLife.91722.2"/></event><event><event-desc>The reviewed preprint was revised.</event-desc><date date-type="reviewed-preprint" iso-8601-date="2024-10-21"><day>21</day><month>10</month><year>2024</year></date><self-uri content-type="reviewed-preprint" xlink:href="https://doi.org/10.7554/eLife.91722.3"/></event></pub-history><permissions><copyright-statement>© 2023, Morton et al</copyright-statement><copyright-year>2023</copyright-year><copyright-holder>Morton et al</copyright-holder><ali:free_to_read/><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><ali:license_ref>http://creativecommons.org/licenses/by/4.0/</ali:license_ref><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife-91722-v1.pdf"/><self-uri content-type="figures-pdf" xlink:href="elife-91722-figures-v1.pdf"/><abstract><p>Identical stimuli can be perceived or go unnoticed across successive presentations, producing divergent behavioral outcomes despite similarities in sensory input. We sought to understand how fluctuations in behavioral state and cortical layer and cell class-specific neural activity underlie this perceptual variability. We analyzed physiological measurements of state and laminar electrophysiological activity in visual area V4 while monkeys were rewarded for correctly reporting a stimulus change at perceptual threshold. Hit trials were characterized by a behavioral state with heightened arousal, greater eye position stability, and enhanced decoding performance of stimulus identity from neural activity. Target stimuli evoked stronger responses in V4 in hit trials, and excitatory neurons in the superficial layers, the primary feed-forward output of the cortical column, exhibited lower variability. Feed-forward interlaminar population correlations were stronger on hits. Hit trials were further characterized by greater synchrony between the output layers of the cortex during spontaneous activity, while the stimulus-evoked period showed elevated synchrony in the feed-forward pathway. Taken together, these results suggest that a state of elevated arousal and stable retinal images allow enhanced processing of sensory stimuli, which contributes to hits at perceptual threshold.</p></abstract><kwd-group kwd-group-type="author-keywords"><kwd>visual attention</kwd><kwd>visual cortex</kwd><kwd>area V4</kwd><kwd>laminar</kwd><kwd>threshold perception</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>Rhesus macaque</kwd></kwd-group><funding-group><award-group id="fund1"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>R01 EY021827</award-id><principal-award-recipient><name><surname>Reynolds</surname><given-names>John H</given-names></name><name><surname>Nandy</surname><given-names>Anirvan S</given-names></name></principal-award-recipient></award-group><award-group id="fund2"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>R01 EY032555</award-id><principal-award-recipient><name><surname>Nandy</surname><given-names>Anirvan S</given-names></name></principal-award-recipient></award-group><award-group id="fund3"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>R01 EY034605</award-id><principal-award-recipient><name><surname>Jadi</surname><given-names>Monika P</given-names></name></principal-award-recipient></award-group><award-group id="fund4"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>R00 EY025026</award-id><principal-award-recipient><name><surname>Jadi</surname><given-names>Monika P</given-names></name></principal-award-recipient></award-group><award-group id="fund5"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000065</institution-id><institution>National Institute of Neurological Disorders and Stroke</institution></institution-wrap></funding-source><award-id>T32 NS007224</award-id><principal-award-recipient><name><surname>Morton</surname><given-names>Mitchell P</given-names></name></principal-award-recipient></award-group><award-group id="fund6"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000065</institution-id><institution>National Institute of Neurological Disorders and Stroke</institution></institution-wrap></funding-source><award-id>T32 NS041228</award-id><principal-award-recipient><name><surname>Morton</surname><given-names>Mitchell P</given-names></name></principal-award-recipient></award-group><award-group id="fund7"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>P30 EY019005</award-id><principal-award-recipient><name><surname>Reynolds</surname><given-names>John H</given-names></name></principal-award-recipient></award-group><award-group id="fund8"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000053</institution-id><institution>National Eye Institute</institution></institution-wrap></funding-source><award-id>P30 EY026878</award-id><principal-award-recipient><name><surname>Nandy</surname><given-names>Anirvan S</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>States of elevated arousal and stable retinal images allow enhanced processing of threshold stimuli in the laminar microcircuits of area V4.</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>publishing-route</meta-name><meta-value>prc</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Physical properties of stimuli strongly influence perception such that low-intensity stimuli are detected infrequently. As intensity increases, detection probability remains low until some perceptual threshold is crossed, after which stimuli are perceived robustly. A psychometric function (<xref ref-type="bibr" rid="bib76">Prins and Kingdom, 2018</xref>; <xref ref-type="bibr" rid="bib97">Watson, 1979</xref>; <xref ref-type="bibr" rid="bib100">Wichmann and Hill, 2001</xref>) mathematically describes this property of perception. Only within a narrow range around the perceptual threshold do stimuli lead to significant trial-to-trial perceptual variance. While many studies present stimuli at threshold (<xref ref-type="bibr" rid="bib34">Herman et al., 2019</xref>; <xref ref-type="bibr" rid="bib43">Levitt, 1971</xref>; <xref ref-type="bibr" rid="bib74">Pins and Ffytche, 2003</xref>; <xref ref-type="bibr" rid="bib79">Ress and Heeger, 2003</xref>), few have probed the laminar cortical microcircuit mechanisms that underlie successful or unsuccessful perception under these conditions (<xref ref-type="bibr" rid="bib52">McCormick et al., 2020</xref>; <xref ref-type="bibr" rid="bib93">van Vugt et al., 2018</xref>).</p><p>Prior studies have characterized how perceived stimuli trigger stronger information propagation from earlier visual areas to higher-order visual and frontal regions (<xref ref-type="bibr" rid="bib34">Herman et al., 2019</xref>; <xref ref-type="bibr" rid="bib93">van Vugt et al., 2018</xref>). This information propagation and sensory processing are strongly influenced by brain states such as arousal and attention (<xref ref-type="bibr" rid="bib52">McCormick et al., 2020</xref>; <xref ref-type="bibr" rid="bib31">Harris and Thiele, 2011</xref>). Arousal has long been recognized for its role in modulating cortical activity (<xref ref-type="bibr" rid="bib44">Livingstone and Hubel, 1981</xref>; <xref ref-type="bibr" rid="bib51">McCormick and Bal, 1997</xref>; <xref ref-type="bibr" rid="bib96">Vinck et al., 2015</xref>) and affecting performance in various sensory tasks (<xref ref-type="bibr" rid="bib2">Aston-Jones and Cohen, 2005</xref>; <xref ref-type="bibr" rid="bib54">McGinley et al., 2015</xref>; <xref ref-type="bibr" rid="bib103">Yerkes and Dodson, 1908</xref>). In visual area V4, a key intermediate region in the ventral visual processing stream (<xref ref-type="bibr" rid="bib27">Goodale and Milner, 1992</xref>; <xref ref-type="bibr" rid="bib67">Mountcastle, 1997</xref>; <xref ref-type="bibr" rid="bib81">Roe et al., 2012</xref>), attention strongly modulates neural activity (<xref ref-type="bibr" rid="bib13">Desimone and Duncan, 1995</xref>; <xref ref-type="bibr" rid="bib49">McAdams and Maunsell, 1999</xref>; <xref ref-type="bibr" rid="bib65">Moran and Desimone, 1985</xref>; <xref ref-type="bibr" rid="bib80">Reynolds et al., 2000</xref>). Attention increases the firing rates of V4 neurons, enhances the reliability of individual neuron firing, and reduces correlated fluctuations among pairs of neurons (<xref ref-type="bibr" rid="bib49">McAdams and Maunsell, 1999</xref>; <xref ref-type="bibr" rid="bib7">Cohen and Maunsell, 2009</xref>; <xref ref-type="bibr" rid="bib58">Mitchell et al., 2007</xref>; <xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>). Brain state dynamics impact both cortical and subcortical structures, contributing to behavior (<xref ref-type="bibr" rid="bib26">Ghosh and Maunsell, 2021</xref>; <xref ref-type="bibr" rid="bib106">Zénon and Krauzlis, 2012</xref>). Fluctuations in attention are reflected in the on- and off-state dynamics of V4 ensembles, which have been shown to correlate with behavioral performance (<xref ref-type="bibr" rid="bib19">Engel et al., 2016</xref>; <xref ref-type="bibr" rid="bib92">van Kempen et al., 2021</xref>).</p><p>The visual cortex has a columnar architecture in which multiple cell classes (<xref ref-type="bibr" rid="bib8">Connors and Gutnick, 1990</xref>; <xref ref-type="bibr" rid="bib46">Markram et al., 2004</xref>; <xref ref-type="bibr" rid="bib57">Migliore and Shepherd, 2005</xref>; <xref ref-type="bibr" rid="bib101">Wonders and Anderson, 2006</xref>; <xref ref-type="bibr" rid="bib105">Zeng and Sanes, 2017</xref>) across the cortical layers (<xref ref-type="bibr" rid="bib67">Mountcastle, 1997</xref>; <xref ref-type="bibr" rid="bib15">Douglas and Martin, 2004</xref>) form distinct sub-populations. These sub-populations have unique and stereotyped patterns of connectivity, thus forming a canonical microcircuit that orchestrates the encoding and flow of information (<xref ref-type="bibr" rid="bib16">Douglas and Martin, 2007</xref>; <xref ref-type="bibr" rid="bib36">Hirsch and Martinez, 2006</xref>). Moreover, these sub-populations contribute uniquely to sensory processing and are differentially modulated by brain states (<xref ref-type="bibr" rid="bib58">Mitchell et al., 2007</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib50">McCormick et al., 1985</xref>; <xref ref-type="bibr" rid="bib73">Pettine et al., 2019</xref>). While it has been shown that attentional modulation varies across cortical layers (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib73">Pettine et al., 2019</xref>; <xref ref-type="bibr" rid="bib55">Mehta et al., 2000a</xref>; <xref ref-type="bibr" rid="bib56">Mehta et al., 2000b</xref>; <xref ref-type="bibr" rid="bib5">Buffalo et al., 2011</xref>; <xref ref-type="bibr" rid="bib20">Ferro et al., 2021</xref>; <xref ref-type="bibr" rid="bib98">Westerberg et al., 2021</xref>; <xref ref-type="bibr" rid="bib99">Westerberg et al., 2022</xref>), the role of these sub-populations in attentive perception at threshold remains poorly understood. Moreover, the influence of physiological states, which may be responsible for different outcomes at threshold, on these sub-populations has not been studied in detail.</p><p>Here, we examine the neural mechanisms that regulate perception at threshold. We specifically focus on the columnar microcircuit mechanisms within area V4. We hypothesized that minor fluctuations in behavioral state, such as arousal and visual sensitivity, and in the activity of neural sub-populations across the layers of the visual cortex, result in different perceptual outcomes at threshold. Specifically, we hypothesized that output layer (II–III, V–VI) sub-populations, ones projecting to higher cortical areas and subcortical structures, would show evidence of improved capacity for stimulus representation during successful perception. We also hypothesized that such successful events would be accompanied by improved information propagation throughout a cortical column. We find that differences in behavioral states and lamina-specific neural states characterize correct and incorrect trials at threshold and explain perceptual variability.</p></sec><sec id="s2" sec-type="results"><title>Results</title><p>To study the neural dynamics responsible for determining whether a stimulus presented at perceptual threshold is perceived, we analyzed behavioral and cortical layer-specific neural data from area V4, collected while monkeys performed a cued attention task (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>). Monkeys were trained to detect an orientation change in one of two Gabor stimuli that were presented concurrently at two spatial locations, and to report having seen the change by making an eye movement to the changed stimulus. Prior to a block of trials, monkeys were cued as to which of the two spatial locations was likely to undergo the orientation change (95% valid cue; presented at the start of each block). During a trial, ‘non-target’ stimuli at a fixed reference orientation were repeatedly presented. Non-targets were turned on for 200 ms at the two spatial locations, and then turned off for a variable interval (200–400 ms). At a random time (1–5 s, mean 3 s) a ‘target’ stimulus, differing in orientation from the non-targets, was presented at one of the locations. If the monkey reported having detected the orientation change by making an eye movement to the location of the target stimulus, it received a juice reward (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, ‘hit’ trial). If the monkey failed to detect the orientation change and instead continued to maintain fixation on the center of the monitor, it was not rewarded (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, ‘miss’ trial). In this study, we focused exclusively on trials in which the target stimulus was presented at the cued location (95% of trials). All figures relate exclusively to trials in which the change occurred at the cued location.</p><fig-group><fig id="fig1" position="float"><label>Figure 1.</label><caption><title>Orientation change detection task at perceptual threshold.</title><p>(<bold>A</bold>) Schematic of task structure. The monkey initiated a trial by fixating on the center of the screen. Two Gabor stimuli (represented by oriented lines) were presented for 200 ms and then turned off for 200–400 ms. This was repeated until, at an unpredictable time, one of the stimuli changed orientation. The monkey could report having seen the change by making an eye movement to the location of the target stimulus to receive a reward (hit trials). If the monkey failed to report the orientation change and maintained fixation on the center of the screen it was not rewarded (miss trials). Before a block of trials, the monkey was cued as to which stimulus was likely to undergo the change (95% valid cue). In 5% of trials the orientation change occurred at the other location (foil trials). Circles indicating the cued location and receptive field are drawn for figure reference only and were not presented during the task. (<bold>B</bold>) Example behavioral psychometric function from one recording session and attention condition. Behavioral performance (hit rate, circles) is presented as a function of orientation change. Data was fitted with a logistic function. The threshold condition, trials with performance halfway between the upper and lower asymptotes of the logistic function, is indicated by the orange box. Error bars represent standard deviation calculated with a jackknife procedure (20 jackknives). The square symbol indicates foil trial performance.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig1-v1.tif"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 1.</label><caption><title>Laminar recordings in V4 (modified from <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>).</title><p>(<bold>A</bold>) An artificial dura (AD) chamber is shown over dorsal V4 in the right hemisphere of Monkey A. The native dura mater was resected and replaced with a silicone-based AD, thereby providing an optically clear window into the cortex. Scale bar = 5 mm. (<bold>B</bold>) An enlarged view of the boxed region in <bold>A</bold> clearly shows the sulci and the microvasculature. sts = superior temporal sulcus, lu = lunate sulcus, io = inferior occipital sulcus. Area V4 lies on the pre-lunate gyrus between the superior temporal and lunate sulci. Scale bar = 2 mm. (<bold>C</bold>) Electrophysiology setup: a plastic stabilizer with a circular aperture is secured in place inside the chamber such that the aperture is centered over the pre-lunate gryus. A 16-channel linear array electrode (electrode spacing 150 μm) is positioned over the center of the gyrus and lowered into the cortex under microscopic guidance. The microvasculature pattern was used as a reference to target different cortical sites across recording sessions. (<bold>D</bold>) Example recording session in Monkey C depicting 12 single-unit waveforms (mean ± s.e.m.) isolated along the cortical column. Teal waveforms correspond to narrow-spiking putative interneurons and gold waveforms correspond to broad-spiking putative excitatory units. (<bold>E</bold>) Stimulus triggered local field potentials (LFPs) obtained by flashing 30 ms high-contrast ring stimuli in the receptive field (RF) of a V4 cortical column. LFP traces averaged across all stimulus repeats are shown color-coded as being part of either the superficial (green), input (gray), or deep (pink) layers. Layer assignment was done after current source density (CSD) analysis. (<bold>F</bold>) CSD calculated as the second spatial derivative of the stimulus-triggered LFPs and displayed as a colored map. The x-axis represents time from stimulus onset; the y-axis represents cortical depth oriented such that the pial surface is at the top and the white matter is at the bottom. Red hues represent current sink, blue hues represent current source. The input layer is identified as the first current sink followed by a reversal to current source. The superficial and deep layers have the opposite sink-source pattern. The CSD map has been spatially smoothed for visualization. (<bold>G</bold>) Stacked contour plots show spatial RFs mapped along each contact point in the laminar probe. The spatial RFs were obtained by applying reverse correlation to the LFP power evoked by sparse pseudo-random sequences of Gabor stimuli. The RFs are well aligned, indicating perpendicular penetration down a cortical column. Zero depth represents the center of the input layer as estimated from the CSD. (<bold>H</bold>) Contrast response functions – spikes rate as a function of stimulus contrast – are shown for two example units identified in a single recording session in Monkey A. Red and blue traces correspond to the attend-in to RF and attend-away from RF conditions, respectively. The dotted lines represent the corresponding background firing rates. The dashed lines are four standard deviations above baseline. A unit was considered as visually responsive, if the contrast response functions exceeded this threshold in both the attention conditions. Mean ± s.e.m. Panels are reproduced from <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig1-figsupp1-v1.tif"/><permissions><copyright-statement>© 2017, Elsevier Inc</copyright-statement><copyright-year>2017</copyright-year><copyright-holder>Elsevier Inc</copyright-holder><license><license-p>Panels are reprinted from <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>, with permission from Elsevier. It is not covered by the CC-BY 4.0 licence and further reproduction of this panel would need permission from the copyright holder.</license-p></license></permissions></fig><fig id="fig1s2" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 2.</label><caption><title>Additional psychometric function examples.</title><p>Example behavioral psychometric functions from nine recording sessions in the attend-in (red) and attend-away (blue) conditions. Behavioral performance (circles) is presented as a function of orientation change. Data was fitted with a logistic function. The threshold of each fitted logistic function is indicated by the vertical dashed line. The square symbol indicates foil trial performance, and the star symbol indicates catch trial performance.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig1-figsupp2-v1.tif"/></fig><fig id="fig1s3" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 3.</label><caption><title>Psychometric function parameters.</title><p>(<bold>A</bold>) Distribution of the orientation threshold of the fitted psychometric functions in the attend-in (red solid) and attend-away (blue hatch) conditions for each recording session. The threshold of the sigmoid may not correspond exactly to an orientation change tested in the experiment. (<bold>B</bold>) Distribution of threshold orientation changes used as the threshold condition in the rest of this study (red = attend-in, blue = attend-away, circle = Monkey A, plus = Monkey C). (<bold>C</bold>) Distribution of slopes for the fitted psychometric functions (same conventions as in <bold>A</bold>). (<bold>D</bold>) Distribution of guess rate for the fitted psychometric functions (same conventions as in <bold>A</bold>). (<bold>E</bold>) Distribution of lapse rate for the fitted psychometric functions (same conventions as in <bold>A</bold>). (<bold>F</bold>) Distribution of catch trial performance in each session (same conventions as in <bold>A</bold>). (<bold>G</bold>) Performance on orientation change-matched regular trials (target presented at cued location) and foil trials (target presented at uncued location) for each session (same conventions as in <bold>B</bold>). (<bold>H</bold>) Distribution of target presentation time for hit (orange) and miss (gray) trials. (<bold>I</bold>) Performance in threshold trials in each recording session split by quartiles (see Materials and methods). Error bars represent standard deviation.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig1-figsupp3-v1.tif"/></fig></fig-group><p>On each trial, the magnitude of the orientation change was drawn from a distribution that spanned multiple levels of difficulty. We fit the behavioral data with a logistic function (<xref ref-type="bibr" rid="bib76">Prins and Kingdom, 2018</xref>) and defined the threshold condition as the orientation change that was closest to the 50% threshold of the fitted psychometric function for that session (<xref ref-type="fig" rid="fig1">Figure 1B</xref>, Materials and methods; see <xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref> for additional examples and <xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3A–G</xref> for logistic fit parameters). We selected this subset of trials for further analysis, since the constant target stimuli in these trials were equally likely to be perceived or not perceived. Target presentation times were not different between hit and miss trials (<xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3H</xref>; p=0.15, Wilcoxon rank sum test). There was a slight difference in threshold trial performance based on time in the session (<xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3I</xref>, p&lt;0.01, permutation test). Performance in the middle of the recording session (second and third quartiles) was higher than in the beginning and end of the session (first and fourth quartiles). Monkeys initiated a median of 905 trials per session (range: 651–1086).</p><p>While monkeys performed this task, single- and multi-unit activity and local field potentials (LFPs) were recorded in area V4 using 16-channel linear array electrodes (Plexon Inc, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A–E</xref>). The array was inserted perpendicular to the cortical surface and spanned the cortical layers. We used current source density (CSD) analysis (<xref ref-type="bibr" rid="bib62">Mitzdorf, 1985</xref>) to estimate the boundaries between the superficial (I–III), input (IV), and deep (V–VI) cortical layers (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1E and F</xref>), and assign individual neurons their layer identity (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>). Single units were classified as either broad-spiking (putative excitatory neurons) or narrow-spiking (putative inhibitory neurons) on the basis of their waveform width using previously published techniques (peak-to-trough duration; <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1D</xref>; see Materials and methods; <xref ref-type="bibr" rid="bib8">Connors and Gutnick, 1990</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib50">McCormick et al., 1985</xref>; <xref ref-type="bibr" rid="bib40">Kawaguchi, 1993</xref>; <xref ref-type="bibr" rid="bib71">Nowak et al., 2003</xref>). Eye position and pupil diameter were also recorded (ISCAN ETL-200). When analyzing pupil diameter and eye position data, we considered all trials in the threshold condition in which the change occurred at the cued location, regardless of whether the cued location was within the receptive field (RF) of the recorded neurons. For all electrophysiological analyses, we only considered trials in which the cued stimulus was within the RF of the recorded neurons, and the stimulus change occurred at the cued location.</p><p>To assess the behavioral impact of variations in arousal and retinal image stability across trials at the threshold condition, we compared pupil diameter and microsaccade incidence across trial outcomes. Larger pupil diameter is thought to be a proxy for elevated alertness and arousal (<xref ref-type="bibr" rid="bib52">McCormick et al., 2020</xref>; <xref ref-type="bibr" rid="bib2">Aston-Jones and Cohen, 2005</xref>; <xref ref-type="bibr" rid="bib54">McGinley et al., 2015</xref>; <xref ref-type="bibr" rid="bib3">Beatty and Lucero-Wagoner, 2000</xref>; <xref ref-type="bibr" rid="bib35">Hess and Polt, 1964</xref>; <xref ref-type="bibr" rid="bib78">Reimer et al., 2014</xref>; <xref ref-type="bibr" rid="bib91">Tang and Higley, 2020</xref>). We found that hit trials were associated with larger pupil diameters compared to miss trials, both before and during non-target and target stimulus presentations (<xref ref-type="fig" rid="fig2">Figure 2A</xref>). We quantified this difference in the estimation statistics framework (<xref ref-type="bibr" rid="bib6">Calin-Jageman and Cumming, 2019</xref>; <xref ref-type="bibr" rid="bib37">Ho et al., 2019</xref>) by comparing effect sizes and using bootstrapping to estimate uncertainty in the differences. We found that the mean of the distribution of pupil diameters associated with hit trials is greater than that associated with misses (<xref ref-type="fig" rid="fig2">Figure 2B</xref>; complementary null hypothesis testing results in <xref ref-type="supplementary-material" rid="supp1">Supplementary file 1a</xref>). Prior work has shown that the optimal state for sensory performance occurs at intermediate levels of arousal, with states of low and hyper arousal associated with decreased performance (<xref ref-type="bibr" rid="bib2">Aston-Jones and Cohen, 2005</xref>; <xref ref-type="bibr" rid="bib54">McGinley et al., 2015</xref>; <xref ref-type="bibr" rid="bib103">Yerkes and Dodson, 1908</xref>; <xref ref-type="bibr" rid="bib10">Cools and D’Esposito, 2011</xref>; <xref ref-type="bibr" rid="bib68">Murphy et al., 2011</xref>; <xref ref-type="bibr" rid="bib77">Rajagovindan and Ding, 2011</xref>; <xref ref-type="bibr" rid="bib94">Vijayraghavan et al., 2007</xref>). In both hit and miss trials, the mean pupil diameter was close to the optimal arousal state for perceptual performance (<xref ref-type="fig" rid="fig2">Figure 2C</xref>; <xref ref-type="bibr" rid="bib54">McGinley et al., 2015</xref>). The average differences in pupil diameter across hit and miss trials reflect differences within the optimal state of intermediate arousal. All results were held for individual animals (<xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2A–C</xref>). Our results thus suggest that hits are more likely to occur during periods of greater arousal.</p><fig-group><fig id="fig2" position="float"><label>Figure 2.</label><caption><title>Hit trials have larger pupil diameter whereas microsaccades more often precede misses.</title><p>(<bold>A</bold>) Normalized pupil diameter for hit and miss trials in the threshold condition. 0 ms corresponds to non-target and target stimulus onset. Mean ± s.e.m. (<bold>B</bold>) Distribution of pupil diameter values associated with hit and miss trials. Pupil diameter was averaged from 100 ms before to 100 ms after non-target and target stimulus onset. Violin plots were generated using kernel smoothing (see Materials and methods). Error bars represent 95% confidence intervals for the mean of each distribution, and the mean difference (blue, right axes). <italic>Inset</italic>: zoomed-in view of the mean difference between hit and miss trials. Black bar represents a 95% confidence interval of the mean difference. Shaded region reflects the distribution of the bootstrapped estimation of the mean difference. (<bold>C</bold>) Histogram of mean pupil diameter around the time of non-target and target stimulus onset (calculated as in <bold>B</bold>). Orange and gray lines represent the mean pupil diameter for hit and miss trials, respectively. (<bold>D</bold>) <italic>Left:</italic> Hit rate for trials with (387 trials) and without (<italic>right</italic>, 1336 trials) a microsaccade detected in the time window 0–400 ms before target onset. <italic>Right:</italic> Bootstrapped estimation of the mean difference in hit proportion in trials with vs without a pre-target microsaccade. Same conventions as in B.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig2-v1.tif"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><label>Figure 2—figure supplement 1.</label><caption><title>Microsaccades are preferentially directed toward the target in correct trials and have a slight correlation with pupil diameter.</title><p>Data is presented for all trials, regardless of orientation change (not just the threshold condition). (<bold>A</bold>) The histograms represent the direction of microsaccades relative to the attended stimulus (left column) or unattended stimulus (right column) in correct (top row) and incorrect (bottom row) trials. Black lines represent the mean (solid) and 99.5% confidence interval (dashed) of the bootstrapped null distribution estimated by pooling correct and incorrect microsaccades. *<italic>p</italic> &lt; 0.005. <italic>Inset:</italic> Schematic for calculation of relative microsaccade direction. Microsaccade is represented by the gray arrow. (<bold>B</bold>) Scatterplot of microsaccade rate vs mean normalized pupil diameter shows a small but statistically significant relationship between the two quantities (<italic>r</italic><sup>2</sup> = 0.006, <italic>p</italic> &lt; 0.001). Each dot is color-coded by trial length. <italic>Inset:</italic> Histogram of mean pupil diameter for all trials, including those less than –2 or greater than 2 (red lines). Only four of the trials with pupil diameters less than –2 or greater than 2 had detected microsaccades. (<bold>C</bold>) Boxplot of microsaccade rates for all trials. (<bold>D</bold>) Boxplot of the time between microsaccades recorded in the same trial. 4% of microsaccades separated by more than 538 ms were excluded as outliers (see Materials and methods).</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig2-figsupp1-v1.tif"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><label>Figure 2—figure supplement 2.</label><caption><title>Single monkey pupil diameter and microsaccade data.</title><p>(<bold>A</bold>) Normalized pupil diameter for hit and miss trails in the threshold condition separated by monkey. 0 ms corresponds to non-target and target stimulus onset. Mean ± s.e.m. (<bold>B</bold>) Distribution of pupil diameter values associated with hit and miss trials separated by monkey. Pupil diameter was averaged from 100 ms before to 100 ms after non-target and target stimulus onset. Violin plots were generated using kernel smoothing (see Materials and methods). Error bars represent 95% confidence intervals for the mean of each distribution, and the mean difference (blue, right axes). (<bold>C</bold>) Histogram of mean pupil diameter around the time of non-target and target stimulus onset for each monkey (calculated as in <bold>B</bold>). Orange and gray lines represent the mean pupil diameter for hit and miss trials, respectively. (<bold>D</bold>) Hit rate for trials with and without a microsaccade detected in the time window 0–400 ms before target onset for each monkey.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig2-figsupp2-v1.tif"/></fig></fig-group><p>Microsaccades, small fixational eye movements of &lt;1° in amplitude that occur during normal fixation, are associated with periods of decreased visual sensitivity due to unstable retinal images (<xref ref-type="bibr" rid="bib14">Dicke et al., 2008</xref>; <xref ref-type="bibr" rid="bib108">Zuber and Stark, 1966</xref>). Microsaccades have been linked to suppressed neural responses in visual areas during perceptual tasks, impairing fine visual discrimination and behavioral performance (<xref ref-type="bibr" rid="bib4">Beeler, 1967</xref>; <xref ref-type="bibr" rid="bib30">Hafed and Krauzlis, 2010</xref>). We grouped trials in the threshold condition based on whether a microsaccade occurred in a 400 ms window preceding the onset of the target stimulus. Most trials with a pre-target microsaccade were misses, whereas the majority of trials without a microsaccade in this window were hits (<xref ref-type="fig" rid="fig2">Figure 2D</xref>; see <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2D</xref> for individual animal plots). There is a strong link between microsaccade direction and attention deployment (<xref ref-type="bibr" rid="bib45">Lowet et al., 2018</xref>; <xref ref-type="bibr" rid="bib18">Engbert and Kliegl, 2003</xref>; <xref ref-type="bibr" rid="bib29">Hafed and Clark, 2002</xref>; <xref ref-type="bibr" rid="bib28">Gowen et al., 2007</xref>; <xref ref-type="bibr" rid="bib25">Galfano et al., 2004</xref>; <xref ref-type="bibr" rid="bib72">Pastukhov and Braun, 2010</xref>; <xref ref-type="bibr" rid="bib104">Yu et al., 2022</xref>). Consistent with previous reports we also find that microsaccades toward the attended stimulus were overrepresented in correct trials (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A</xref>, upper left). Conversely, microsaccades toward the attended stimulus were underrepresented in incorrect trials (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A</xref>, lower left). There was a very low but statistically significant negative correlation between pupil diameter and microsaccade rate (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref>, <italic>r</italic><sup>2</sup> = 0.006, <italic>p</italic> &lt; 0.001). Microsaccade rates and inter-microsaccade times are reported in <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1C and D</xref>. Overall, these results suggest that successful trials at threshold are significantly more likely to occur during a state of greater arousal and improved visual sensitivity.</p><p>Having established that hit trials are more likely to occur in states of elevated arousal and visual stability, we investigated whether hits are characterized by differential information processing in V4. We first examined the ability to discriminate target stimuli from non-target stimuli using the firing rates of single- and multi-unit V4 neurons in each of the three identified cortical layers (superficial, input, and deep). A linear decoder could better discriminate targets from non-targets in hits compared to misses (<xref ref-type="fig" rid="fig3">Figure 3A</xref>; see <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref> for individual animal plots), suggesting differences in firing rates across these trial types. This improved stimulus discriminability was consistent across all three layers (<xref ref-type="fig" rid="fig3">Figure 3A</xref>).</p><fig-group><fig id="fig3" position="float"><label>Figure 3.</label><caption><title>Target stimuli evoke higher firing rates in hit trials.</title><p>Rows correspond to different layers (top = superficial, middle = input, bottom = deep). (<bold>A</bold>) Performance for decoding targets from non-targets from single units and multi-units in each layer. Points in the left section of each plot show the decoding performance for each of the 20 different cross-validations. The right section for each layer shows the bootstrapped estimation of the difference between decoding performance between hits and misses. Half-violin plots show the bootstrapped distribution of the difference, and black dots and bars represent the mean and 95% confidence intervals of the difference in decoding performance. Chance levels, determined by shuffling target and non-target identity, were subtracted from the raw decoding performance values. (<bold>B</bold>) Non-target population (single- and multi-unit) PSTH of visually responsive neurons for the hit (orange) and miss (dark-gray) trials in the threshold condition (mean ± s.e.m.). The horizontal black bar indicates the time and duration of stimulus presentation. (<bold>C</bold>) As in B but for target stimuli. The star indicates the time at which firing rates in the input layer first differ significantly between hit and miss trials. Vertical lines represent the mean time at which firing rates for each neuron rise above the 95% confidence interval of their baseline activity (see also <xref ref-type="fig" rid="fig3s3">Figure 3—figure supplement 3C</xref>). (<bold>D</bold>) Bootstrapped estimation of the paired mean difference in target stimulus-evoked firing rate between hit and miss trials in the time window 60–260 ms (red dotted box in C) after target stimulus onset. Shaded regions represent the bootstrapped estimation of the paired mean difference in firing rate (hit-miss), and black lines are 95% confidence intervals. Plots include data from both single and multi-units, separated by layer (top = superficial, middle = input, bottom = deep). (<bold>E</bold>) As in D, bootstrapped estimation of the paired mean difference in firing rate for hit trials compared to miss trials in the target stimulus-evoked period, but only for single units broken up by cell class (gold = broad, teal = narrow).</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig3-v1.tif"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><label>Figure 3—figure supplement 1.</label><caption><title>Single monkey decoding performance.</title><p>(<bold>A</bold>) Performance for decoding targets from non-targets from single units and multi-units in the superficial layer for individual monkeys. Points in the left section of each plot show the decoding performance for each of the 20 different cross-validations. The right section for each layer shows the bootstrapped estimation of the difference between decoding performance in hits and misses. Half-violin plots show the bootstrapped distribution of the difference, and black dots and bars represent the mean and 95% confidence intervals of the difference in decoding performance. Chance levels, determined by shuffling target and non-target identity, were subtracted from the raw decoding performance values. (<bold>B</bold>) Same as <bold>A</bold> but for the input layer. (<bold>C</bold>) Same as <bold>A</bold> but for the deep layers.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig3-figsupp1-v1.tif"/></fig><fig id="fig3s2" position="float" specific-use="child-fig"><label>Figure 3—figure supplement 2.</label><caption><title>Single monkey firing rate data.</title><p>(<bold>A</bold>) Single monkey non-target population (single- and multi-unit) PSTH of visually responsive neurons for the hit (orange) and miss (dark-gray) trials in the threshold condition (mean ± s.e.m.). (<bold>B</bold>) As in <bold>A</bold> but for target stimuli. Red box indicates the time window used for analysis in <xref ref-type="fig" rid="fig3">Figure 3C</xref>.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig3-figsupp2-v1.tif"/></fig><fig id="fig3s3" position="float" specific-use="child-fig"><label>Figure 3—figure supplement 3.</label><caption><title>Firing rates for individual neurons and reaction time in threshold condition.</title><p>(<bold>A</bold>) Target stimulus-evoked normalized firing rates in hit and miss trials for each recorded single- and multi-unit cluster in hit and miss trials. Clusters are divided by layer: left = superficial, middle = input, right = deep. Related to <xref ref-type="fig" rid="fig3">Figure 3B</xref>. Each line represents the mean firing rate in response to target stimuli in hit and miss trials for a given unit. Data is color-coded by unit type (gold = broad, teal = narrow, gray = multi-unit). See Materials and methods for normalization method. (<bold>B</bold>) Mean reaction time (time from target stimulus onset to saccade onset) for hit trials in the threshold condition in both attend-in (red) and attend-away (blue) trials for each recording session. (<bold>C</bold>) Mean time at which single- and multi-unit firing rates rise significantly above pre-stimulus baseline in response to target stimuli in hit and miss trials. Mean ± s.e.m.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig3-figsupp3-v1.tif"/></fig></fig-group><p>Elevated stimulus-evoked firing rates would indicate a stronger representation of the stimulus that could cause this improved discriminability in hits. We compared the firing rates of all neurons (single- and multi-units) recorded in each cortical layer across hit and miss trials. For non-target stimuli, firing rates were equivalent for hits and misses in both the pre-stimulus (0–200 ms before stimulus onset) and stimulus-evoked (60–260 ms following stimulus onset) periods (<xref ref-type="fig" rid="fig3">Figure 3B</xref>; see <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2</xref> for individual animal plots). For the target stimulus, firing rates were once again equivalent in the pre-stimulus period, but hit trials were characterized by elevated firing across cortical layers in the stimulus-evoked period (<xref ref-type="fig" rid="fig3">Figure 3C and D</xref>). Broad- and narrow-spiking neurons in both the input and deep layers respond more to target stimuli in hit trials, and trend toward elevated firing rates in the superficial layers during hits (<xref ref-type="fig" rid="fig3">Figure 3E</xref>). The average firing rate in response to target stimuli for each neuron is shown in <xref ref-type="fig" rid="fig3s3">Figure 3—figure supplement 3A</xref> for both hit and miss trials. It is important to note that the stimuli presented to the animals were identical for both hits and misses. Moreover, the responses to the target stimuli occur early, and elevated firing in hits emerges at the time of expected V4 response latencies (70–100 ms; <xref ref-type="fig" rid="fig3s3">Figure 3—figure supplement 3C</xref>), and thus cannot be attributed to elevated levels of firing due to subsequent saccade planning in these trials (<xref ref-type="fig" rid="fig3s3">Figure 3—figure supplement 3B</xref>; expected &gt;200 ms; <xref ref-type="bibr" rid="bib89">Steinmetz and Moore, 2014</xref>).</p><p>Variability in response reflects how reliably information is encoded by a neural population. Lower baseline variability can enhance the ability of neurons to encode stimulus differences. We calculated the Fano factor, a mean-normalized measure of trial-to-trial variability in firing, for single units in our population (<xref ref-type="fig" rid="fig4">Figure 4A</xref>; see <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref> for individual animal plots). We find that broad-spiking units in the superficial layer exhibited lower Fano factor during the pre-stimulus period in hit trials (0–60 ms before non-target stimulus onset, <xref ref-type="fig" rid="fig4">Figure 4B</xref>), indicating this population of neurons fired more reliably when the animal correctly detected the orientation change. This was not the case for broad-spiking neurons in other layers (<xref ref-type="fig" rid="fig4">Figure 4B</xref>) or narrow-spiking neurons (<xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2</xref>).</p><fig-group><fig id="fig4" position="float"><label>Figure 4.</label><caption><title>Broad-spiking neurons in the superficial layer have decreased variability in hit trials.</title><p>(<bold>A</bold>) Rows correspond to different layers (top = superficial, middle = input, bottom = deep). The Fano factor of broad-spiking putative excitatory neurons for the hit and miss trials in the threshold condition (mean ± s.e.m.). There is a significant decrease in variability for the hit trials prior to stimulus onset only in the superficial layer. 0 ms corresponds to non-target stimulus onset. The average Fano factor within a 60 ms time window (red dashed box) prior to non-target stimulus onset is plotted in <bold>B</bold>. (<bold>B</bold>) <italic>Top:</italic> Fano factor modulation index for each broad-spiking neurons recorded in each layer, averaged in the 60 ms preceding non-target stimulus onset. <italic>Bottom:</italic> Bootstrapped estimation of the mean difference of the Fano factor modulation index from zero in each of the three layers. Colored curves represent the estimated bootstrapped distribution. Black dots and lines reflect the mean and 95% confidence intervals of the distributions.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig4-v1.tif"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><label>Figure 4—figure supplement 1.</label><caption><title>Single monkey Fano factor data.</title><p>Rows correspond to different layers (top = superficial, middle = input, bottom = deep) and columns correspond to monkeys. The Fano factor of broad-spiking putative excitatory neurons for the hit and miss trials in the threshold condition in each monkey (mean ± s.e.m.).</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig4-figsupp1-v1.tif"/></fig><fig id="fig4s2" position="float" specific-use="child-fig"><label>Figure 4—figure supplement 2.</label><caption><title>Narrow-spiking neurons do not have decreased variability in hit trials.</title><p><italic>Top:</italic> Fano factor modulation index for each narrow-spiking neurons recorded in each layer, averaged in the 60 ms preceding non-target stimulus onset. <italic>Bottom:</italic> Bootstrapped estimation of the mean difference of the Fano factor modulation index from zero in each of the three layers. Colored curves represent the estimated bootstrapped distribution. Black dots and lines reflect the mean and 95% confidence intervals of the distributions.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig4-figsupp2-v1.tif"/></fig></fig-group><p>We next wanted to test how the relationship between spiking activity and LFPs may differ across hits and misses. Spike-LFP synchrony can reflect cortical processing and both within- and inter-areal coordination (<xref ref-type="bibr" rid="bib24">Fries, 2009</xref>; <xref ref-type="bibr" rid="bib23">Fries et al., 2008</xref>; <xref ref-type="bibr" rid="bib87">Siapas et al., 2005</xref>). We calculated the PPC (<xref ref-type="bibr" rid="bib95">Vinck et al., 2010</xref>), a frequency-resolved measure of spike-LFP phase-locking, for single and multi-units relative to their local LFP signal during the pre-stimulus period (0–200 ms before non-target stimulus onset, <xref ref-type="fig" rid="fig5">Figure 5A</xref>; see <xref ref-type="fig" rid="fig5s2">Figure 5—figure supplement 2</xref> for individual animal plots). We averaged PPC values at low (3–12 Hz), medium (15–25 Hz), and high (30–80 Hz) frequency bands (superficial and input: <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1A and B</xref>; deep: <xref ref-type="fig" rid="fig5">Figure 5B</xref>, <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1</xref>). Deep-layer neurons exhibit reduced low-frequency phase-locking in hit trials than in misses (<xref ref-type="fig" rid="fig5">Figure 5B</xref>). This desynchronization during hits is similar to prior reports of desynchronization due to the deployment of attention (<xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>).</p><fig-group><fig id="fig5" position="float"><label>Figure 5.</label><caption><title>Deep-layer neurons are phase-locked to low-frequency rhythms in miss trials.</title><p>(<bold>A</bold>) Pairwise phase consistency (PPC) of single and multi-units in each layer to the local field potential (LFP) signal recorded from the same channel in hit and miss trials at threshold. PPC was calculated in the pre-stimulus period (0–200 ms before stimulus onset). Dashed red line indicates a PPC of 0, below which there is no consistent relationship between spikes and LFP phase. (<bold>B</bold>) Bootstrapped estimation plot for the paired mean difference in PPC for deep-layer neurons over three frequency bands: 3–12 Hz, 15–25 Hz, 30–80 Hz. Curves represent the bootstrapped distribution for the paired difference, and black dots and vertical lines represent the mean and 95% confidence intervals for the paired mean difference.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig5-v1.tif"/></fig><fig id="fig5s1" position="float" specific-use="child-fig"><label>Figure 5—figure supplement 1.</label><caption><title>Additional pairwise phase consistency (PPC) data.</title><p>(<bold>A–B</bold>) <italic>Top:</italic> Raw PPC values calculated for clusters recorded in the superficial (<bold>A</bold>) and (<bold>B</bold>) input layers in hit and miss trials, averaged into three frequency bands, 3–12 Hz, 15–25 Hz, and 30–80 Hz. PPC was calculated using the local field potential (LFP) recorded on the same channel as the spikes. <italic>Bottom:</italic> Bootstrapped estimation of the paired mean difference in PPC across hit and miss trials for each frequency band. Note that although there appears to be a difference in high-frequency PPC in the superficial layer, this population does not have significantly positive PPC in either condition, indicating that there is no phase-locking in either hits or misses. (<bold>C</bold>) Raw PPC values for neurons recorded in the deep layer in hit and miss trials, averaged into the same three frequency bands. Related to <xref ref-type="fig" rid="fig5">Figure 5B</xref>.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig5-figsupp1-v1.tif"/></fig><fig id="fig5s2" position="float" specific-use="child-fig"><label>Figure 5—figure supplement 2.</label><caption><title>Single monkey pairwise phase consistency (PPC) data.</title><p>Rows correspond to different layers (top = superficial, middle = input, bottom = deep) and columns correspond to monkeys. PPC of single and multi-units in each layer to the local field potential (LFP) signal recorded from the same channel in hit (orange) and miss (gray) trials at threshold. PPC was calculated in the pre-stimulus period (0–200 ms before stimulus onset). Dashed red line indicates a PPC of 0, below which there is no consistent relationship between spikes and LFP phase.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig5-figsupp2-v1.tif"/></fig></fig-group><p>Our results at the individual neuron or neural-sub-population levels suggest enhanced processing of perceived stimuli. However, it is the concerted activity among neural sub-populations that ultimately determine information flow through the laminar cortical circuit. We turned to canonical correlation analysis (CCA) to investigate the strength of feed-forward communication across layers (<xref ref-type="bibr" rid="bib61">Mitra, 2007</xref>). CCA has previously been used to describe interactions among multiple cortical areas (<xref ref-type="bibr" rid="bib86">Semedo et al., 2022</xref>). We performed CCA on each pair of layers: input to superficial, input to deep, and superficial to deep, where the two elements in each pair correspond to the upstream and downstream layers respectively (<xref ref-type="fig" rid="fig6">Figure 6A</xref>). We refer to the results of CCA as population correlations. Interlaminar feed-forward population correlations were higher in hits than in misses in both the pre-stimulus and stimulus-evoked periods (<xref ref-type="fig" rid="fig6">Figure 6B and C</xref>). This suggests that feed-forward information flow through the column is more effective in hits than in misses.</p><fig id="fig6" position="float"><label>Figure 6.</label><caption><title>Hit trials are characterized by stronger feed-forward interlaminar population correlations.</title><p>(<bold>A</bold>) Canonical correlation analysis (CCA)-based population correlation as a function of time and interlaminar delay during the pre-stimulus and stimulus-evoked periods in hit and miss trials in an example session. (<bold>B</bold>) Mean feed-forward population correlation in each session. Color indicates the monkey (blue = Monkey A, yellow = Monkey C). (<bold>C</bold>) Bootstrapped estimation plot for the paired mean difference in population correlation for each pair of layers and time window (pre-stimulus or stimulus-evoked). Curves represent the bootstrapped distribution for the paired difference, and black dots and vertical lines represent the mean and 95% confidence intervals for the paired mean difference.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig6-v1.tif"/></fig><p>To further investigate interlaminar communication, we analyzed interlaminar synchrony as signatures of differential information flow between hit and miss trials. Spike-spike coherence (SSC) is a frequency-resolved measure of the degree to which two spike trains fluctuate together (<xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>; <xref ref-type="bibr" rid="bib60">Mitra and Pesaran, 1999</xref>). We measured interlaminar SSC for spike trains from pairs of cortical layers, each spike train being comprised of all recorded action potentials in a given layer (see Materials and methods). We computed interlaminar SSC separately for hit and mis trials in both the pre-stimulus (0–200 ms before non-target stimulus onset, <xref ref-type="fig" rid="fig7">Figure 7A</xref>) and non-target stimulus-evoked (60–260 ms after non-target stimulus onset, <xref ref-type="fig" rid="fig7">Figure 7C</xref>) periods, matching the firing rates across hit and miss trials separately for the pre-stimulus and non-target stimulus-evoked conditions (see <xref ref-type="fig" rid="fig7s1">Figure 7—figure supplement 1</xref> for individual animal plots). We averaged SSC for each pair of layers across three frequency bands, 3–12 Hz, 15–25 Hz, and 30–80 Hz (<xref ref-type="fig" rid="fig7">Figure 7B and D</xref>).</p><fig-group><fig id="fig7" position="float"><label>Figure 7.</label><caption><title>Greater interlaminar coherence in hit trials in the pre-stimulus and non-target stimulus-evoked periods.</title><p>Rows correspond to different pairs of layers (top = superficial-input, middle = superficial-deep, bottom = input-deep). (<bold>A</bold>) Multi-unit interlaminar spike-spike coherence (SSC) calculated in the 200 ms before non-target stimulus onset in hit and miss trials (solid lines, mean ± s.e.m.). Firing rates were matched across hit and miss trials. Dashed lines represent coherence calculated with shuffled trial identities (mean ± s.e.m.). (<bold>B</bold>) Bootstrapped estimation plot for the paired mean difference in SSC for each pair of layers averaged over three frequency bands: 3–12 Hz, 15–25 Hz, 30–80 Hz. Curves represent the bootstrapped distribution for the paired difference, and black dots and vertical lines represent the mean and 95% confidence intervals for the paired mean difference. (<bold>C</bold>) Interlaminar SSC in the non-target stimulus-evoked period (60–260 ms after stimulus onset). Same conventions as in <bold>A</bold>. (<bold>D</bold>) Bootstrapped estimation plot for the paired mean difference in SSC for each pair of layers averaged over three frequency bands. Same conventions as in <bold>B</bold>.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig7-v1.tif"/></fig><fig id="fig7s1" position="float" specific-use="child-fig"><label>Figure 7—figure supplement 1.</label><caption><title>Single animal interlaminar coherence.</title><p>Rows correspond to different pairs of layers (top = superficial-input, middle = superficial-deep, bottom = input-deep). (<bold>A</bold>) Single monkey multi-unit interlaminar spike-spike coherence (SSC) calculated in the 200 ms before non-target stimulus onset in hit and miss trials (solid lines, mean ± s.e.m.). Firing rates were matched across hit and miss trials. Dashed lines represent coherence calculated with shuffled trial identities (mean ± s.e.m.). (<bold>B</bold>) Interlaminar SSC in the non-target stimulus-evoked period in each monkey (60–260 ms after stimulus onset). Same conventions as in A.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig7-figsupp1-v1.tif"/></fig></fig-group><p>Overall, hit trials have greater interlaminar SSC compared to misses at almost all frequencies (<xref ref-type="fig" rid="fig7">Figure 7B and D</xref>). In the pre-stimulus period, the strongest SSC difference between hits and misses was observed between the superficial and deep layers across all frequencies (<xref ref-type="fig" rid="fig7">Figure 7B</xref>, middle panel). This implies greater synchrony of the output layers of the cortex during hit trials. In contrast, this pattern was directionally the same during the non-target stimulus-evoked period, but stronger in the other layer pairs, with greater SSC differences being found in pairs that involve the input layer (<xref ref-type="fig" rid="fig7">Figure 7D</xref>, top and bottom). This may reflect a higher degree of stimulus-driven feed-forward information propagation during hit trials. When comparing across time (pre-stimulus vs non-target stimulus-evoked), layers, and frequency band, there was a significant interaction effect of layer pair and time window (three-way ANOVA, <italic>p</italic> = 0.0075).</p><p>Finally, we sought to compare the predictive power of our results on the monkey’s perceptual performance. We created a generalized linear model (GLM) to regress behavioral outcome from the pupil diameter, number of microsaccades in the pre-target window, and average target-evoked multi-unit firing rate in each of the three layers (see Materials and methods; <xref ref-type="bibr" rid="bib11">Davis et al., 2020</xref>). Other reported measures (Fano factor, PPC, interlaminar population correlations, and SSC) that we could not estimate reliably on a single-trial basis were not considered in the GLM analysis. Pre-target microsaccades were by far the strongest predictor of performance (<inline-formula><mml:math id="inf1"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>w</mml:mi><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mn>1.3116</mml:mn><mml:mo>;</mml:mo><mml:mspace width="thinmathspace"/><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>6.0757</mml:mn><mml:mi>e</mml:mi><mml:mo>−</mml:mo><mml:mn>08</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>). Input layer firing rate also significantly predicted perception (<inline-formula><mml:math id="inf2"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>w</mml:mi><mml:mi>e</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0.3276</mml:mn><mml:mo>;</mml:mo><mml:mspace width="thinmathspace"/><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>0.020068</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>). Superficial firing rate, deep firing rate, and pupil diameter were not significant predictors (<xref ref-type="supplementary-material" rid="supp1">Supplementary file 1b</xref>, all <italic>p</italic> &gt; 0.5). This indicates that, among the variables that we could estimate reliably on a single-trial basis, stable retinal images in the pre-target window are critical for behavioral performance, and elevated firing in the input layer is the most reliable physiological signature of a perceived stimulus. GLM fit parameters can be found in <xref ref-type="supplementary-material" rid="supp1">Supplementary file 1c</xref>.</p></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>We investigated the physiological processes responsible for variable behavioral outcomes at perceptual threshold. Controlling for <italic>both</italic> the attentive instruction (thus minimizing large-scale attentional effects) and the stimulus condition that elicited performance at a threshold level allowed us to examine the physiological and neural correlates that underlie correct vs incorrect behavioral outcomes. While this study cannot disentangle the independent roles of behavioral state fluctuations and neural fluctuations in determining behavioral outcomes, evidence suggests that differences in both are associated with hits. We found multiple lines of evidence which suggest that a state of higher arousal and eye position stability and the accompanying enhanced processing of visual stimuli contributes to accurate perception in hit trials (<xref ref-type="fig" rid="fig8">Figure 8</xref>).</p><fig id="fig8" position="float"><label>Figure 8.</label><caption><title>Conceptual model for stimulus processing at perceptual threshold.</title><p>(<bold>A</bold>) Hit trials have a larger pupil diameter and fewer pre-target microsaccades, reflecting a state of increased arousal and greater eye position stability. Conversely, miss trials show decreased arousal and eye position stability. (<bold>B</bold>) In the spontaneous pre-stimulus period, hits are characterized by decreased variability in superficial layer broad-spiking neurons, which we hypothesize is reflective of lower membrane potential (V<sub>m</sub>) variability (inset). Hit trials are also characterized by greater synchrony between the superficial and deep layers (indicated by thicker arrows), which could be reflecting a stronger top-down influence on the cortical column. (<bold>C</bold>) In the stimulus-evoked period there is greater interlaminar synchrony between pairs that include the input layer (represented by thicker arrows), which we propose reflects improved feed-forward propagation of information. We propose these state differences in hits contribute to elevated firing rates in response to target stimuli, particularly in the superficial layers (inset), resulting in a higher-fidelity output to downstream areas. E=excitatory; I=inhibitory; s=superficial; i=input; d=deep.</p></caption><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-fig8-v1.tif"/></fig><p>Pupil diameter is elevated in hit trials (<xref ref-type="fig" rid="fig2">Figure 2A–C</xref>; <xref ref-type="fig" rid="fig8">Figure 8A</xref>), and prior studies have shown that pupil diameter is strongly linked to arousal and alertness (<xref ref-type="bibr" rid="bib3">Beatty and Lucero-Wagoner, 2000</xref>; <xref ref-type="bibr" rid="bib35">Hess and Polt, 1964</xref>; <xref ref-type="bibr" rid="bib91">Tang and Higley, 2020</xref>). This provides evidence that a state of higher arousal may contribute to improved sensory processing. The much lower hit rate in trials with a microsaccade preceding the target (<xref ref-type="fig" rid="fig2">Figure 2D</xref>; <xref ref-type="fig" rid="fig8">Figure 8A</xref>) and our GLM analysis show that stability of retinal images is critical for accurate discrimination at threshold. It is unlikely that these two measures are reflecting the same phenomenon, as there is a very weak correlation between them over the course of a trial (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref>).</p><p>There is a strong link between oculomotor control and attentional deployment (<xref ref-type="bibr" rid="bib64">Moore and Zirnsak, 2017</xref>; <xref ref-type="bibr" rid="bib84">Schafer and Moore, 2011</xref>; <xref ref-type="bibr" rid="bib63">Moore and Fallah, 2001</xref>). In this study, hits and misses differ in their behavioral responses, with hit trials being characterized by a saccade to the target stimulus. Almost all of our neural results reflect differences around the time of non-target stimulus presentations during which the monkeys maintained fixation at the center of the screen and, therefore, were hundreds of milliseconds prior to saccade planning and execution in the case of hit trials. Trials in which saccades were made to non-target stimuli were excluded from analysis, as were trials in which the monkey made a saccade to the target too soon after its presentation to have been a behavioral response to stimulus perception (see Materials and methods). The analysis of microsaccade occurrence focused on the window just before target stimulus presentation and before monkeys could begin oculomotor planning. Only the analysis of neural responses to target stimuli appears in conjunction with divergent oculomotor behavior between the hit (saccade) and miss (no saccade) trials. However, here too firing rates diverge much earlier, particularly in the input layer, than would be consistent with the effects of saccade planning (<xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1B</xref>; <xref ref-type="bibr" rid="bib89">Steinmetz and Moore, 2014</xref>).</p><p>Non-target stimulus contrasts were slightly different between hits and misses (mean: 33.1% in hits, 34.0% in misses, permutation test, <inline-formula><mml:math id="inf3"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>0.02</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>), but the contrast of the target was higher in hits compared to misses (mean: 38.7% in hits, 27.7% in misses, permutation test, <inline-formula><mml:math id="inf4"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1.6</mml:mn><mml:mspace width="thinmathspace"/><mml:mi>e</mml:mi><mml:mo>−</mml:mo><mml:mn>31</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>). To control for potential effects of stimulus contrast, firing rates were first normalized by contrast before performing the analyses reported in <xref ref-type="fig" rid="fig3">Figure 3</xref>. For all other results, we considered only non-target stimuli, which had very minor differences in contrast (&lt;1%) across hits and misses. In fact, this minor difference was in the opposite direction of our results with mean contrast being slightly higher for misses. While we cannot completely rule out any other effects of stimulus contrast, the normalization in <xref ref-type="fig" rid="fig3">Figure 3</xref> and minor differences for non-target stimuli should minimize them.</p><p>A body of evidence (see <xref ref-type="bibr" rid="bib48">Martinez-Conde et al., 2013</xref>, for review) suggests that microsaccades directed toward a target stimulus reflect attention-related processing and performance (<xref ref-type="bibr" rid="bib45">Lowet et al., 2018</xref>; <xref ref-type="bibr" rid="bib18">Engbert and Kliegl, 2003</xref>; <xref ref-type="bibr" rid="bib29">Hafed and Clark, 2002</xref>; <xref ref-type="bibr" rid="bib28">Gowen et al., 2007</xref>; <xref ref-type="bibr" rid="bib25">Galfano et al., 2004</xref>; <xref ref-type="bibr" rid="bib72">Pastukhov and Braun, 2010</xref>; <xref ref-type="bibr" rid="bib104">Yu et al., 2022</xref>). In our dataset, during the pre-target period, microsaccades toward the attended stimulus were overrepresented in correct trials (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A</xref>, upper left). Conversely, microsaccades toward the attended stimulus were underrepresented in incorrect trials (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A</xref>, lower left). Microsaccades directed toward the location of the eventual target may reflect elevated attentional deployment that can compensate for reduced sensitivity due to a higher incidence of microsaccades.</p><p>Our electrophysiological findings and their laminar patterns associated with hit trials <italic>within</italic> a cued attention state mirror several previous findings that are associated with the deployment of covert spatial attention. Attention has long been known to increase firing rates in V4 (<xref ref-type="bibr" rid="bib49">McAdams and Maunsell, 1999</xref>; <xref ref-type="bibr" rid="bib58">Mitchell et al., 2007</xref>; <xref ref-type="bibr" rid="bib88">Spitzer et al., 1988</xref>), and there is evidence that this increase occurs in all cortical layers in V4 (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>). We find improved target vs non-target discriminability in hits (<xref ref-type="fig" rid="fig3">Figure 3A</xref>) across all cortical layers. Additionally, elevated target-evoked firing rates in hits occur across all layers in conjunction with elevated arousal (<xref ref-type="fig" rid="fig3">Figure 3B–D</xref>; <xref ref-type="fig" rid="fig8">Figure 8C</xref>). Attention reduces the variability in the firing of V4 neurons, and this reduction is thought to contribute to the improved information coding capacity of a population of neurons (<xref ref-type="bibr" rid="bib7">Cohen and Maunsell, 2009</xref>; <xref ref-type="bibr" rid="bib58">Mitchell et al., 2007</xref>; <xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib66">Moreno-Bote et al., 2014</xref>). The reduction in Fano factor among broad-spiking superficial-layer neurons in hit trials mirrors the effects of attention (<xref ref-type="fig" rid="fig4">Figure 4</xref>). Multiple lines of evidence suggest broad- and narrow-spiking correspond to putative excitatory and inhibitory neurons respectively. Narrow-spiking neurons exhibit higher firing rates, which corresponds well with inhibitory interneuron (<xref ref-type="bibr" rid="bib8">Connors and Gutnick, 1990</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib50">McCormick et al., 1985</xref>; <xref ref-type="bibr" rid="bib9">Contreras and Palmer, 2003</xref>; <xref ref-type="bibr" rid="bib21">Foehring et al., 1991</xref>; <xref ref-type="bibr" rid="bib75">Povysheva et al., 2006</xref>). Repolarization times in broad-spiking neurons are also longer, as they are in excitatory pyramidal neurons (<xref ref-type="bibr" rid="bib50">McCormick et al., 1985</xref>; <xref ref-type="bibr" rid="bib71">Nowak et al., 2003</xref>; <xref ref-type="bibr" rid="bib33">Hasenstaub et al., 2005</xref>). Since these neurons are putative projection neurons to downstream cortical areas, this reduction in Fano factor may indicate increased reliability in stimulus encoding that could contribute to hits. Our finding is also in agreement with previous reports of higher variability in representations of unperceived stimuli in humans (<xref ref-type="bibr" rid="bib85">Schurger et al., 2010</xref>). Synchronous neural activity appears to modulate perceptual and cognitive ability in a variety of contexts (<xref ref-type="bibr" rid="bib1">Abbas et al., 2018</xref>; <xref ref-type="bibr" rid="bib22">Fries et al., 2001</xref>; <xref ref-type="bibr" rid="bib82">Rohenkohl et al., 2018</xref>; <xref ref-type="bibr" rid="bib102">Worden et al., 2000</xref>). We found that deep-layer neurons exhibit less low-frequency phase-locking in hit trials (<xref ref-type="fig" rid="fig5">Figure 5</xref>). This is consistent with prior studies that find an attention-mediated reduction in the power spectrum of the spike-triggered-averaged LFP (<xref ref-type="bibr" rid="bib22">Fries et al., 2001</xref>).</p><p>In examining interlaminar population synchrony, we found that hit trials were characterized by stronger feed-forward interactions across the cortical column (<xref ref-type="fig" rid="fig6">Figure 6</xref>). This state of improved interlaminar information flow could be a result of neuromodulatory or top-down processes that maintain the cortex in a state of sustained depolarization corresponding to a state of higher arousal during hits (<xref ref-type="bibr" rid="bib52">McCormick et al., 2020</xref>; <xref ref-type="bibr" rid="bib54">McGinley et al., 2015</xref>). Our examination of interlaminar synchrony revealed two interesting and complementary patterns: hits were associated with greater coherence between the superficial and deep layers during spontaneous activity in the pre-stimulus period (<xref ref-type="fig" rid="fig7">Figure 7A and B</xref>; <xref ref-type="fig" rid="fig8">Figure 8B</xref>); in contrast, we found enhanced coherence between the input layer and both the output layers (superficial and deep) in the stimulus-evoked period during hits (<xref ref-type="fig" rid="fig7">Figure 7C and D</xref>; <xref ref-type="fig" rid="fig8">Figure 8C</xref>). Increased superficial-deep coherence in the pre-stimulus period could be the result of the same neuromodulatory or top-down processes. Increased synchrony between the input layer and the output layers during the stimulus-evoked period provides further evidence of stronger information propagation through the cortical circuit, and hence with improved stimulus detection (<xref ref-type="bibr" rid="bib47">Marshel et al., 2019</xref>). In contrast to broad global synchrony or local correlated fluctuations, which may signal a default state of minimal processing or decreased information coding capacity (<xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>; <xref ref-type="bibr" rid="bib90">Steriade et al., 1993</xref>; <xref ref-type="bibr" rid="bib42">Krosigk von et al., 1993</xref>; <xref ref-type="bibr" rid="bib107">Zohary et al., 1994</xref>), these patterns of interlaminar coherence that we found suggest that successful perception at threshold is mediated by pathway-specific modulation of information flow through the laminar cortical circuit.</p><p>Prior studies showing decreased correlations under attention typically do not contain laminar information (<xref ref-type="bibr" rid="bib7">Cohen and Maunsell, 2009</xref>; <xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>) or only consider decreased correlations within a layer (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>). In contrast, the correlation and synchrony analyses presented here are interlaminar, which we expect could reflect improved information processing in a column, similar to principles of communication across areas (<xref ref-type="bibr" rid="bib86">Semedo et al., 2022</xref>).</p><p>Taken together, our results provide insight about how information about a threshold stimulus may successfully propagate through a cortical column and influence sensory perception. Lower baseline variability among broad-spiking superficial layer neurons and decreased low-frequency synchronous activity in the deep layers could be indicative of improved capacity to encode sensory information. Higher target-evoked firing rates and elevated interlaminar synchrony could enhance the propagation of this encoded signal. These results associate pre-stimulus baseline state differences with enhanced cortical processing in the stimulus-evoked period.</p><p>Several studies have examined how information flow differs for perceived and unperceived stimuli at a more macroscopic scale (<xref ref-type="bibr" rid="bib34">Herman et al., 2019</xref>; <xref ref-type="bibr" rid="bib93">van Vugt et al., 2018</xref>). <xref ref-type="bibr" rid="bib93">van Vugt et al., 2018</xref>, recorded from three brain regions, V1, V4, and dorsolateral prefrontal cortex, while a monkey performed a stimulus detection task at threshold. Their work supports the model that feed-forward propagation of sensory information from the visual cortex to the PFC causes a non-linear ‘ignition’ of association areas resulting in conscious perception (<xref ref-type="bibr" rid="bib12">Dehaene and Changeux, 2011</xref>). <xref ref-type="bibr" rid="bib34">Herman et al., 2019</xref>, found that conscious human perception triggers a wave of activity propagation from occipital to frontal cortex while switching off default mode and other networks. Our study provides insight into the functions of the cortical microcircuit at the columnar level that could reflect these large-scale sweeping activity changes in perception.</p><p>Overall, we identified substantial layer-specific differences in cortical activity between hits and misses at perceptual threshold, leading to the following conceptual model (<xref ref-type="fig" rid="fig8">Figure 8</xref>). During spontaneous activity, the state of elevated arousal and eye position stability during hits (<xref ref-type="fig" rid="fig8">Figure 8A</xref>) is manifested by increased interlaminar synchrony between the superficial and deep layers (<xref ref-type="fig" rid="fig8">Figure 8B</xref>, thicker orange arrows), which we propose is due to top-down influences. We predict that decreased firing variability in broad-spiking neurons in the superficial layer is caused by a lower variability in membrane potential closer to the action potential threshold among these neurons (<xref ref-type="fig" rid="fig8">Figure 8B</xref>, inset). Elevated feed-forward propagation in the stimulus-evoked period (<xref ref-type="fig" rid="fig8">Figure 8C</xref>) and a membrane potential closer to action potential threshold could both contribute to higher firing rates in the output layers of the cortex (<xref ref-type="fig" rid="fig8">Figure 8C</xref>, inset), and are indicative of greater fidelity of stimulus processing in hits. These physiological differences in the laminar microcircuit likely contribute to successful perceptual discrimination at threshold.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><table-wrap id="keyresource" position="anchor"><label>Key resources table</label><table frame="hsides" rules="groups"><thead><tr><th align="left" valign="top">Reagent type (species) or resource</th><th align="left" valign="top">Designation</th><th align="left" valign="top">Source or reference</th><th align="left" valign="top">Identifiers</th><th align="left" valign="top">Additional information</th></tr></thead><tbody><tr><td align="left" valign="top">Other</td><td align="left" valign="top">Monkey</td><td align="left" valign="top">This paper</td><td align="left" valign="top"/><td align="left" valign="top">Species (<italic>Macaca mulatta</italic>)</td></tr><tr><td align="left" valign="top">Software, algorithm</td><td align="left" valign="top">MATLAB</td><td align="left" valign="top">MathWorks</td><td align="left" valign="top"/><td align="left" valign="top">R2019a</td></tr><tr><td align="left" valign="top">Software, algorithm</td><td align="left" valign="top">Cortex</td><td align="left" valign="top">NIMH</td><td align="left" valign="top"/><td align="left" valign="top"><ext-link ext-link-type="uri" xlink:href="http://www.cortex.salk.edu/">http://www.cortex.salk.edu/</ext-link></td></tr></tbody></table></table-wrap><sec id="s4-1"><title>Surgical procedures</title><p>Surgical procedures have been described in detail previously (<xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>; <xref ref-type="bibr" rid="bib70">Nassi et al., 2015</xref>; <xref ref-type="bibr" rid="bib83">Ruiz et al., 2013</xref>). In brief, an MRI-compatible low-profile titanium chamber was placed over the pre-lunate gyrus, on the basis of preoperative MRI in two rhesus macaques (right hemisphere in Monkey A, left hemisphere in Monkey C). The native dura mater was then removed, and a silicone-based optically clear artificial dura (AD) was inserted, resulting in an optical window over dorsal V4 (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A and B</xref>). Antibiotic (amikacin or gentamicin) soaked gauze was placed in the chamber between recording sessions to prevent bacterial growth. All procedures were approved by the Institutional Animal Care and Use Committee and conformed to NIH guidelines (Salk Institute protocol number 14-00014).</p></sec><sec id="s4-2"><title>Electrophysiology</title><p>At the beginning of each recording session a plastic insert, with an opening for targeting electrodes, was lowered into the chamber and secured. This served to stabilize the recording site against cardiac pulsations. Neurons were recorded from cortical columns in dorsal V4 using 16-channel linear array electrodes (‘laminar probes’, Plexon Inc, Plexon V-probe). The laminar probes were mounted on adjustable X-Y stages attached to the recording chamber and positioned over the center of the pre-lunate gyrus under visual guidance through a microscope (Zeiss Inc, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1C</xref>). This ensured that the probes were maximally perpendicular to the surface of the cortex and thus had the best possible trajectory to make a perpendicular penetration down a cortical column. Across recording sessions, the probes were positioned over different sites along the center of the gyrus in the parafoveal region of V4 with RF eccentricities between 2° and 7° of visual angle. Care was taken to target cortical sites with no surface microvasculature, with surface microvasculature used as reference so that the same cortical site was not targeted across recording sessions. The probes were advanced using a hydraulic microdrive (Narishige Inc) to first penetrate the AD and then through the cortex under microscopic visual guidance. Probes were advanced until the point that the top-most electrode (toward the pial surface) registered LFP signals. At this point, the probe was retracted by about 100–200 μm to ease the dimpling of the cortex due to the penetration. This procedure greatly increased the stability of the recordings and increased the neuronal yield in the superficial electrodes.</p><p>The distance from the tip of the probes to the first electrode contact was either 300 or 700 μm. The inter-electrode distance was 150 μm, thus minimizing the possibility of recording the same neural spikes in adjacent recording channels. Electrical signals were recorded extracellularly from each channel. These were then amplified, digitized, and filtered either between 0.5 Hz and 2.2 kHz (LFPs) or between 250 Hz and 8 kHz (spikes) and stored using the Multichannel Acquisition Processor system (MAP system, Plexon Inc). Spikes and LFPs were sampled at 40 and 10 kHz, respectively. LFP signals were further low-pass filtered with a sixth-order Butterworth filter with 300 Hz cut-off and down-sampled to 1 kHz for further analysis. Spikes were classified as either multi-unit clusters or isolated single units using the Plexon Offline Sorter software program. Single units were identified based on two criteria: (1) if they formed an identifiable cluster, separated from noise and other units, when projected into the principal components of waveforms recorded on that electrode and (2) if the inter-spike interval distribution had a well-defined refractory period. Single units were classified as either narrow-spiking (putative interneurons) or broad-spiking (putative pyramidal cells) based on methods described in detail previously (<xref ref-type="bibr" rid="bib58">Mitchell et al., 2007</xref>; <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>). Specifically, only units with waveforms having a clearly defined peak <italic>preceded</italic> by a trough were potential candidates. The distribution of trough-to-peak duration was clearly bimodal (Hartigan’s Dip Test, <inline-formula><mml:math id="inf5"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>0.012</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>) (<xref ref-type="bibr" rid="bib32">Hartigan and Hartigan, 1985</xref>). Units with trough-to-peak duration less than 225 μs were classified as narrow-spiking units; units with trough-to-peak duration greater than 225 μs were classified as broad-spiking units (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1D</xref>; teal = narrow, gold = broad).</p><p>Data was collected over 32 sessions (23 sessions in Monkey A, 9 in Monkey C), yielding a total of 413 single units (146 narrow-spiking, 267 broad-spiking) and 296 multi-unit clusters. Per session unit yield was considerably higher in Monkey C compared to Monkey A, resulting in a roughly equal contribution of both monkeys toward the population data.</p></sec><sec id="s4-3"><title>Task, stimuli, and inclusion criteria</title><p>Stimuli were presented on a computer monitor placed 57 cm from the eye. Eye position was continuously monitored with an infrared eye tracking system (ISCAN ETL-200). Trials were aborted if eye position deviated more than 1° (degree of visual angle [‘dva’]) from fixation. Experimental control was handled by NIMH Cortex software (<ext-link ext-link-type="uri" xlink:href="http://www.cortex.salk.edu/">http://www.cortex.salk.edu/</ext-link>). Eye position (all sessions) and pupil diameter (18/32 sessions) data were concurrently recorded and stored using the MAP system.</p><sec id="s4-3-1"><title>Receptive field mapping</title><p>At the beginning of each recording session, neuronal RFs were mapped using subspace reverse correlation in which Gabor (eight orientations, 80% luminance contrast, spatial frequency 1.2 cycles/°, Gaussian half-width 2°) or ring stimuli (80% luminance contrast) appeared at 60 Hz while the monkeys maintained fixation. Each stimulus appeared at a random location selected from an 11×11 grid with 1° spacing in the appropriate visual quadrant. Spatial receptive maps were obtained by applying reverse correlation to the evoked LFP signal at each recording site. For each spatial location in the 11×11 grid, we calculated the time-averaged power in the stimulus-evoked LFP (0–200 ms after each stimulus flash) at each recording site. The resulting spatial map of LFP power was taken as the spatial RF at the recording site. For the purpose of visualization, the spatial RF maps were smoothed using spline interpolation and displayed as stacked contours plots of the smoothed maps (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1G</xref>). All RFs were in the lower visual quadrant (lower-left in Monkey A, lower-right in Monkey C) and with eccentricities between 2 and 7 dva.</p></sec><sec id="s4-3-2"><title>CSD mapping</title><p>In order to estimate the laminar identity of each recording channel, we used a CSD mapping procedure (<xref ref-type="bibr" rid="bib62">Mitzdorf, 1985</xref>). Monkeys maintained fixation while 100% luminance contrast ring stimuli were flashed (30 ms) centered at the estimated RF overlap region across all channels. The size of the ring was scaled to about three-quarters of the estimated diameter of the RF. CSD was calculated as the second spatial derivative of the flash-triggered LFPs (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1E</xref>). The resulting time-varying traces of current across the cortical layers can be visualized as CSD maps (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1F</xref>; maps have been spatially smoothed with a Gaussian kernel for aid in visualization). Red regions depict current sinks in the corresponding region of the cortical laminae; blue regions depict current sources. The input layer (Layer 4) was identified as the first current sink followed by a reversal to current source. The superficial (Layers 1–3) and deep (Layers 5–6) layers had opposite sink-source patterns. LFPs and spikes from the corresponding recording channels were then assigned to one of three layers: superficial, input, or deep.</p></sec><sec id="s4-3-3"><title>Attention task</title><p>In the main experiment, monkeys had to perform an attention-demanding orientation change detection task (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). While the monkey maintained fixation, two achromatic Gabor stimuli (orientation optimized per recording session, spatial frequency 1.2 cycles/°, Gaussian half-width 2°, 6 contrasts randomly chosen from an uniform distribution of luminance contrasts, Monkey A: <inline-formula><mml:math id="inf6"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>10</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>18</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>26</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>34</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>42</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>50</mml:mn><mml:mi mathvariant="normal">%</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:math></inline-formula>, Monkey C: <inline-formula><mml:math id="inf7"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>8</mml:mn><mml:mspace width="thinmathspace"/><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>20</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>28</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>36</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>44</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>52</mml:mn><mml:mo>,</mml:mo><mml:mspace width="thinmathspace"/><mml:mn>60</mml:mn><mml:mi mathvariant="normal">%</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:math></inline-formula> or <inline-formula><mml:math id="inf8"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mspace width="thinmathspace"/><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>30</mml:mn><mml:mo>,</mml:mo><mml:mn>40</mml:mn><mml:mo>,</mml:mo><mml:mn>50</mml:mn><mml:mo>,</mml:mo><mml:mn>60</mml:mn><mml:mo>,</mml:mo><mml:mn>70</mml:mn><mml:mo>,</mml:mo><mml:mn>80</mml:mn><mml:mi mathvariant="normal">%</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:math></inline-formula>) were flashed on for 200 ms and off for a variable period chosen from a uniform distribution between 200 and 400 ms. One of the Gabor’s was flashed at the RF overlap region, the other at a location of equal eccentricity across the vertical meridian. The range of stimulus contrasts was the same at both locations. At the beginning of a block of trials, the monkey was spatially cued (‘instruction trials’) to covertly attend to one of these two spatial locations. During these instruction trials, the stimuli were only flashed at the spatially cued location. No further spatial cue was presented during the rest of the trials in a block. At an unpredictable time drawn from a truncated exponential distribution (minimum 1 s, maximum 5 s, mean 3 s), one of the two stimuli changed in orientation. The monkey was rewarded for making a saccade to the location of orientation change. The monkey was rewarded for only those saccades where the saccade onset time was within a window of 100–400 ms after the onset of the orientation change. The orientation change occurred at the cued location with 95% probability and at the uncued location with 5% probability (‘foil trials’). We controlled task difficulty by varying the degree of orientation change (<inline-formula><mml:math id="inf9"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:mrow><mml:mtext>ori</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:math></inline-formula>), which was randomly chosen from one of the following: 1°, 2°, 3°, 4°, 6°, 8°, 10°, and 12°. The orientation change in the foil trials was fixed at 4°. These foil trials allowed us to assess the extent to which the monkey was using the spatial cue, with the expectation that there would be an impairment in performance and slower reaction times compared to the case in which the change occurred at the cued location. If no change occurred before 5 s, the monkey was rewarded for maintaining fixation (‘catch trials’, 13% of trials). We refer to all stimuli at the baseline orientation as ‘non-targets’ and the stimulus flash with the orientation change as the ‘target’. Monkeys initiated a median of 905 trials (range of 651–1086).</p></sec><sec id="s4-3-4"><title>Inclusion criteria</title><p>Of the 413 single units, we included only a subset of neurons that were visually responsive for further analysis. For each neuron we calculated its baseline firing-rate for each attention condition (attend into RF [‘attend-in’ or ‘IN’], attend away from RF [‘attend-away’ or ‘AWAY’]) from a 200 ms window before a stimulus flash. We also calculated the neuron’s contrast response function for both attention conditions (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1H</xref>). This was calculated as the firing rate over a window between 60 and 200 ms after stimulus onset and averaged across all stimulus flashes (restricted to non-targets) of a particular contrast separately for each attention condition. A neuron was considered visually responsive if any part of the contrast response curves exceeded the baseline rate by four standard deviations for both attention conditions. This left us with 274 single units (84 narrow-spiking, 190 broad-spiking) and 217 multi-unit clusters for further analysis.</p></sec></sec><sec id="s4-4"><title>Data analysis</title><sec id="s4-4-1"><title>Behavioral analysis</title><p>For each orientation change condition <inline-formula><mml:math id="inf10"><mml:msub><mml:mrow><mml:mo>∆</mml:mo></mml:mrow><mml:mrow><mml:mtext>ori</mml:mtext></mml:mrow></mml:msub></mml:math></inline-formula>, we calculated the hit rate as the ratio of the number of trials in which the monkey correctly identified the target by making a saccadic eye movement to the location of the target over the number of trials in which the target was presented. The hit rate as a function of <inline-formula><mml:math id="inf11"><mml:msub><mml:mrow><mml:mo>∆</mml:mo></mml:mrow><mml:mrow><mml:mtext>ori</mml:mtext></mml:mrow></mml:msub></mml:math></inline-formula>,yields a behavioral psychometric function (<xref ref-type="fig" rid="fig1">Figure 1B</xref>). We performed this analysis independently for each recording day for each monkey, yielding a similar but distinct psychometric function for every session. Psychometric functions were fitted with a smooth logistic function (<xref ref-type="bibr" rid="bib76">Prins and Kingdom, 2018</xref>). Error bars were obtained by a jackknife procedure (20 jackknives, 5% of trials left out for each jackknife). Performance for the foil trials were calculated similarly as the hit rate for trials in which the orientation change occurred at the uncued location (<xref ref-type="fig" rid="fig1">Figure 1B</xref>, square symbol). When fitting the psychometric function, we did not include the contrast of the target stimulus as a variable, and fits were calculated by including trials with target stimuli of all tested contrasts. For each fitted psychometric function in both the attend-in and attend-away conditions, we calculated the threshold of the fitted logistic function (i.e. the <inline-formula><mml:math id="inf12"><mml:msub><mml:mrow><mml:mo>∆</mml:mo></mml:mrow><mml:mrow><mml:mtext>ori</mml:mtext></mml:mrow></mml:msub></mml:math></inline-formula> at which performance was mid-way between the lower and upper asymptotes). Because the threshold of the fitted function always lies somewhere on the axis of <inline-formula><mml:math id="inf13"><mml:msub><mml:mrow><mml:mo>∆</mml:mo></mml:mrow><mml:mrow><mml:mtext>ori</mml:mtext></mml:mrow></mml:msub></mml:math></inline-formula>, but not exactly at an orientation change presented to the subject, we then defined the threshold condition as the subset of trials in which the orientation change of the target stimulus was closest to the threshold of the fitted function (<xref ref-type="fig" rid="fig1">Figure 1B</xref>). We restricted further analysis to this threshold condition. For this threshold condition we identified the trials in which the monkey correctly identified the target as ‘hit’ trials and those in which the monkey failed to identify the target as ‘miss’ trials. Analysis of behavior, pupil diameter, and microsaccades was conducted on both the attend-in and attend-away conditions; all electrophysiological analysis was applied only to the attend-in condition.</p><p>To compare the effect of target timing across hits and misses, we determined the time between trial initiation and target presentation in all trials in the threshold condition. For comparison purposes, the contrast values of all target and non-target stimuli presented in the threshold condition were compared in hit and miss trials using a permutation test.</p></sec><sec id="s4-4-2"><title>Pupil diameter</title><p>The raw pupil diameter measurements from the infrared eye-tracking system could differ across days due to external factors such as display monitor illumination. To control for this, we normalized the raw data by a z-score procedure separately for each session (using the mean and standard deviation of all measurements during the session). We analyzed normalized pupil diameter traces for hit and miss trials in the threshold condition, over a time window from 100 ms before to 100 ms after all stimulus presentations (non-target and target), excluding the first stimulus presentation in a trial. The first stimulus was excluded to avoid pupil diameter changes due to the pupillary near response caused by acquiring fixation (<xref ref-type="bibr" rid="bib53">McDougal and Gamlin, 2015</xref>). The pupil diameter was averaged over this time period and compared across conditions using bootstrap estimation and t-test. Distribution violin plots were generated using kernel density estimation (<xref ref-type="bibr" rid="bib38">Hoffmann, 2015</xref>) (bandwidth(hit)=0.0801, bandwidth(miss)=0.0648). When analyzing pupil diameter, we did not include the contrast of the target and non-target stimuli as a variable, and all analyses were performed by including stimuli of all tested contrasts.</p></sec><sec id="s4-4-3"><title>Microsaccade analysis</title><p>Saccadic eye movements were detected using ClusterFix (<xref ref-type="bibr" rid="bib41">König and Buffalo, 2014</xref>). We identified microsaccades by filtering for eye movements with amplitudes between 0.1 and 1 dva. We then split all trials in the threshold condition into two groups: those in which a microsaccade was detected in the 400 ms preceding the target stimulus presentation, and those without a detected microsaccade. We calculated the hit rate for trials within those two groups. For all trials in which a target stimulus was presented at the attended location, we determined the direction of all microsaccades in the 400 ms period preceding target presentation, relative to both the attended and unattended stimuli. The relative microsaccade direction was defined as the angle between two vectors: the one defined by the eye positions at the beginning and end of the microsaccade, and the vector from the initial eye position to the center of the stimulus (calculated separately for attended and unattended stimuli). Relative microsaccade directions were grouped into 12 bins from 0° to 360°. The distribution of relative microsaccade directions were calculated separately for correct and incorrect trials, relative to both the attended and unattended stimuli (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A</xref>).</p><p>We next created a null distribution of relative microsaccade direction. This was done by pooling together microsaccades from correct and incorrect trials and then sampling with replacement from this pooled data (bootstrap procedure [<xref ref-type="bibr" rid="bib17">Efron and Tibshirani, 1993</xref>]; 1000 samples). The number of microsaccades chosen for each sample was the same as the number in correct or incorrect trials, respectively. These bootstrapped samples were used to create 99.5% confidence intervals for the count of microsaccades expected in each of the 12 bins. A bin was considered significantly different from chance if it’s true count fell outside this confidence interval.</p><p>We calculated microsaccade rate for an entire trial by dividing the total number of detected microsaccades in the whole trial by the trial length (4592 total trials). The Pearson correlation between microsaccade rate and mean normalized pupil diameter (see above) for the trial was calculated for all trials with pupil diameter data, regardless of trial type or outcome (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref>). Not pictured in <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref> but included in correlation analysis were trials with a mean normalized pupil diameter greater than 2 or less than –2 (~4% of trials). Only four of these trials were longer than 1 s, out of which two trials contained detected microsaccades. The mean pupil diameter in these trials is shown in <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref>, inset. Inter-microsaccade time was calculated as the time between microsaccade onset of microsaccades detected in the same trial. 4% of microsaccades separated by &gt;538 ms are excluded from <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1D</xref> as they were more than 1.5× the interquartile range above the third quartile. When analyzing microsaccades, we did not include the contrast of the target and non-target stimuli as a variable, and all analyses were performed by including stimuli of all tested contrasts.</p></sec><sec id="s4-4-4"><title>Decoding analysis</title><p>For each single- or multi-unit neuron, we extracted spike counts from 60 to 260 ms following all non-target or target stimulus onsets in the threshold condition. Using these spike counts, we fit a Poisson distribution to estimate the mean firing rate for each neuron in each of four stimulus conditions (non-targets and targets in hit and miss trials). Responses to stimuli of all contrasts were pooled together to create the Poisson fit. We then created a pseudo-population of neurons in each layer. We generated spike counts drawn from the fitted Poisson distributions to create synthetic spike counts for target and non-target stimuli in hits and misses separately (1000 repeats). We used linear discriminant analysis to decode target from non-target stimuli. Decoders were trained separately for hit and miss trials. The procedure was repeated for a 20-fold cross-validation. We calculated the chance performance by training the decoder on data generated from all trials in the threshold condition (both hits and misses) and with shuffled labels (target or non-target).</p></sec><sec id="s4-4-5"><title>Firing rate</title><p>Firing rates were normalized per neuron to that neuron’s maximum stimulus-evoked response to each contrast before being combined across contrasts and trial types. PSTH of firing rates were calculated in 30 ms bins shifted in 5 ms increments. For each tested contrast level, we calculated each neuron’s maximum average firing rate in any bin from 0 to 260 ms after stimulus onset in hit and miss trials together. We divided all bins in the entire PSTH for that contrast by this maximum firing rate to obtain the normalized firing rate for each contrast. We then combined PSTHs across contrasts and split the data between hit and miss trials. We averaged stimulus-evoked firing rates from 60 to 260 ms following non-target or target stimulus presentations. We used bootstrapped estimation to compare firing rates in hit and miss trials in a paired comparison. This was done for all single- and multi-unit clusters, as well as broad- and narrow-spiking single units in each layer. Firing rates were also compared across hit and miss trials by paired t-test for each group. To calculate the time of firing rate divergence, we compared the difference in each single neuron or multi-unit firing rate in the two conditions over time. At each time point, we performed a Wilcoxon rank sum test comparing the firing rates across hits and misses, and defined the divergence point as the first time the firing rates were significantly different. Divergence was calculated separately for each layer. To determine the time at which firing rates rose above baseline levels, we used bootstrapping to estimate 95% confidence intervals for each neuron’s pre-stimulus firing rate (0–100 ms before target or non-target stimulus onset). We then calculated the target PSTH for each neuron in 30 ms bins shifted in 5 ms increments. We defined the response latency as the first time bin in which the neuron’s firing rate in the PSTH exceeded the upper limit of the 95% confidence interval of baseline firing. We calculated the response latency independently for hit and miss trials for each neuron.</p></sec><sec id="s4-4-6"><title>Fano factor</title><p>Trial-to-trial variability was estimated by the Fano factor, which is the ratio of the variance of the spike counts across trials over the mean of the spike counts for each broad- and narrow-spiking single unit. The Fano factor was calculated over non-overlapping 20 ms time bins in a window from 200 ms prior to each non-target flash onset to 200 ms after each non-target flash onset for hit and miss trials in the threshold condition. To compare across conditions, we calculated the Fano factor modulation index (MI), defined as<disp-formula id="equ1"><mml:math id="m1"><mml:mrow><mml:mrow><mml:mi mathvariant="normal">M</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">I</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>F</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:mi>F</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>F</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>F</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula></p><p>where <inline-formula><mml:math id="inf14"><mml:msub><mml:mrow><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math id="inf15"><mml:msub><mml:mrow><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> represent the Fano factor for a given unit in hit and miss trials respectively at each point in time with respect to non-target stimulus onset. The Fano factor MI was averaged from 0 to 60 ms prior to non-target stimulus onset (disregarding stimulus contrast) and compared across trial types in the threshold condition for each sub-population.</p></sec><sec id="s4-4-7"><title>Pairwise phase consistency</title><p>We calculated PPC (<xref ref-type="bibr" rid="bib95">Vinck et al., 2010</xref>) for single and multi-units in the non-target pre-stimulus period (0–200 ms preceding onset regardless of stimulus contrast) in trials in the threshold condition. Although PPC is unbiased by spike count, we set a threshold of 50 spikes for analysis so that only units with enough spikes for a reliable estimate of PPC were included (superficial: <inline-formula><mml:math id="inf16"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>26</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>, input: <inline-formula><mml:math id="inf17"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>41</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>, deep: <inline-formula><mml:math id="inf18"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>64</mml:mn></mml:mrow></mml:mstyle></mml:math></inline-formula>). LFP phase was calculated using Morlet wavelets. PPC for each unit was calculated for the phase of the LFP recorded on the same channel and averaged in three frequency bands (3–12 Hz, 15–25 Hz, and 30–80 Hz). PPC was calculated separately for hit and miss trials and compared across trial outcomes by t-test, corrected for multiple comparisons.</p></sec><sec id="s4-4-8"><title>Canonical correlation analysis</title><p>We used CCA (<xref ref-type="bibr" rid="bib86">Semedo et al., 2022</xref>; <xref ref-type="bibr" rid="bib39">Hotelling, 1992</xref>) to capture the correlation between layers at different time periods and with different amounts of temporal delays. We considered all possible combinations of pairs of layers. We took two windows of activity, one in each layer, in either the pre-stimulus period (0–200 ms before non-target stimulus onset) or stimulus-evoked period (60–260 ms after non-target stimulus onset), including responses to stimuli of all contrasts. Window length was 50 ms and the window was advanced in 10 ms steps. The activity within each window was then binned using 10 ms bins. We reported correlation associated with the first two canonical pairs, and calculated it separately for hit and miss trials. For each pair of layers, we limited CCA to sessions in which there were at least two neurons recorded in both layers. Feed-forward (FF) layer pairs were defined as follows: input to superficial, input to deep, and superficial to deep. The correlations along the FF signaling pathways (<inline-formula><mml:math id="inf19"><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>) were calculated as the mean correlation at positive delays:<disp-formula id="equ2"><mml:math id="m2"><mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:munder><mml:mo>∑</mml:mo><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:munder><mml:mi>C</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mfrac></mml:mstyle></mml:mrow></mml:math></disp-formula></p><p>where t is the time following response onset, <inline-formula><mml:math id="inf20"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mstyle></mml:math></inline-formula> is the interlaminar delay involved between windows from two layers, <inline-formula><mml:math id="inf21"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mstyle></mml:math></inline-formula> is the corresponding correlation. <inline-formula><mml:math id="inf22"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:mrow><mml:mo>&gt;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:math></inline-formula> is the number of positive delays investigated. The values in <xref ref-type="fig" rid="fig6">Figure 6B</xref> correspond to <inline-formula><mml:math id="inf23"><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> in each session.</p></sec><sec id="s4-4-9"><title>Spike-spike coherence</title><p>For each recording session, all spikes recorded from visually responsive single and multi-units in each layer were combined into a single spike train for that layer (layer multi-unit). Separately for both the pre-stimulus and non-target stimulus-evoked periods, we randomly deleted spikes from the layer multi-unit with a higher firing rate so that the firing rates were matched across hit and miss trials. SSC was calculated for each of the three possible pairs of layer multi-units in each session for both the pre-stimulus (0–200 ms preceding stimulus onset) and non-target stimulus evoked period (60–260 ms following non-target stimulus onset, including all stimuli regardless of contrast) separately for hit and miss trials using Chronux (NW = 1; K=1; <ext-link ext-link-type="uri" xlink:href="http://chronux.org">http://chronux.org</ext-link>; <xref ref-type="bibr" rid="bib61">Mitra, 2007</xref>; <xref ref-type="bibr" rid="bib60">Mitra and Pesaran, 1999</xref>). To control for differences in firing rates across hit and miss trials we used a rate matching procedure (<xref ref-type="bibr" rid="bib59">Mitchell et al., 2009</xref>). For estimation statistics, interlaminar SSC values were calculated for each frequency and subsequently averaged across three frequency bands: 3–12 Hz, 5–15 Hz, and 30–80 Hz and compared across hit and miss trials for each pair of layers in each recording session. For null hypothesis testing, we calculated the SSC modulation index, defined as<disp-formula id="equ3"><mml:math id="m3"><mml:mrow><mml:mrow><mml:mi mathvariant="normal">M</mml:mi><mml:mi mathvariant="normal">I</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>h</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:math></disp-formula></p><p>The SSC MI was calculated for each frequency and subsequently averaged across three frequency bands: 3–12 Hz, 5–15 Hz, and 30–80 Hz. MI values for each frequency band were compared to zero by t-test, Bonferroni corrected for multiple comparisons. We tested for interaction effects with a three-factor ANOVA, with frequency, pair of layers, and time window (pre or post stimulus) as factors. We calculated a shuffled distribution of SSC by shuffling the trial identities of the spikes in one of the layers in the pair. We then calculated SSC with the shuffled trial identities. This procedure was repeated 10 times to create the shuffled distribution.</p></sec><sec id="s4-4-10"><title>GLM quantification</title><p>To compare how well our results can predict behavioral performance, we fit a GLM to the response of the monkeys in trials in the threshold condition (<xref ref-type="bibr" rid="bib11">Davis et al., 2020</xref>). We included five regressors in our analysis: (1) average pupil diameter during the trial, (2) number of microsaccades in the pre-target window (0–400 ms before target stimulus onset), and average target-evoked multi-unit firing rate in the (3) superficial, (4) input, and (5) deep layers. We calculated the average target-evoked firing rate by averaging the firing rate of all single- and multi-units in a given layer 60–260 ms after target stimulus onset in each trial. In order to be able to compare weights across regressors, each regressor was transformed into a z-score before being included in the model. We fit the GLM using a logit link function, using the predictors to regress the categorical binary trial outcome (hit or miss). A total of 309 trials were included in the GLM. Stimulus contrast was not a variable in the GLM.</p></sec></sec></sec></body><back><sec sec-type="additional-information" id="s5"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="COI-statement" id="conf1"><p>No competing interests declared</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>Formal analysis, Writing - original draft</p></fn><fn fn-type="con" id="con2"><p>Software</p></fn><fn fn-type="con" id="con3"><p>Formal analysis</p></fn><fn fn-type="con" id="con4"><p>Funding acquisition, Writing - review and editing</p></fn><fn fn-type="con" id="con5"><p>Conceptualization</p></fn><fn fn-type="con" id="con6"><p>Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing - review and editing</p></fn></fn-group><fn-group content-type="ethics-information"><title>Ethics</title><fn fn-type="other"><p>All procedures were approved by the Institutional Animal Care and Use Committee and conformed to NIH guidelines.(Salk Institute protocol number 14-414 00014).</p></fn></fn-group></sec><sec sec-type="supplementary-material" id="s6"><title>Additional files</title><supplementary-material id="supp1"><label>Supplementary file 1.</label><caption><title>Supplementary tables.</title><p>(a) Corresponding null hypothesis testing results. Null hypothesis testing results corresponding to the estimation statistics-based analyses. (b) Generalized linear model (GLM) coefficient values. Coefficients and significance values of the variables used in the GLM analysis. (c) GLM summary. Additional summary statistics of the GLM analysis.</p></caption><media xlink:href="elife-91722-supp1-v1.docx" mimetype="application" mime-subtype="docx"/></supplementary-material><supplementary-material id="mdar"><label>MDAR checklist</label><media xlink:href="elife-91722-mdarchecklist1-v1.docx" mimetype="application" mime-subtype="docx"/></supplementary-material></sec><sec sec-type="data-availability" id="s7"><title>Data availability</title><p>This study used data from a previously published study: <xref ref-type="bibr" rid="bib69">Nandy et al., 2017</xref>. Data used in this study has been deposited in Zenodo: <ext-link ext-link-type="uri" xlink:href="https://zenodo.org/records/13968327">https://zenodo.org/records/13968327</ext-link>.</p><p>The following dataset was generated:</p><p><element-citation publication-type="data" specific-use="isSupplementedBy" id="dataset1"><person-group person-group-type="author"><name><surname>Morton</surname><given-names>M</given-names></name><name><surname>Denagamage</surname><given-names>S</given-names></name><name><surname>Blume</surname><given-names>I</given-names></name><name><surname>Reynolds</surname><given-names>J</given-names></name><name><surname>Jadi</surname><given-names>M</given-names></name><name><surname>Nandy</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2024">2024</year><data-title>Brain state and cortical layer-specific mechanisms underlying perception at threshold</data-title><source>Zenodo</source><pub-id pub-id-type="doi">10.5281/zenodo.13968327</pub-id></element-citation></p></sec><ack id="ack"><title>Acknowledgements</title><p>This research was supported by NIH R01 EY021827 to JHR and ASN, NIH R01 EY032555, NARSAD Young Investigator Grant, Ziegler Foundation Grant, Yale Orthwein Scholar Funds &amp; Lawrence Family Young Investigator Funds to ASN, NIH R01 EY034605 and NIH R00 EY025026 to MPJ, and by NEI core grants for vision research P30 EY019005 to the Salk Institute and P30 EY026878 to Yale University. MM was supported by training grants T32-NS007224 and T32-NS041228 to Yale University. We would like to thank Catherine Williams and Mat LeBlanc for their excellent animal care.</p></ack><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abbas</surname><given-names>AI</given-names></name><name><surname>Sundiang</surname><given-names>MJM</given-names></name><name><surname>Henoch</surname><given-names>B</given-names></name><name><surname>Morton</surname><given-names>MP</given-names></name><name><surname>Bolkan</surname><given-names>SS</given-names></name><name><surname>Park</surname><given-names>AJ</given-names></name><name><surname>Harris</surname><given-names>AZ</given-names></name><name><surname>Kellendonk</surname><given-names>C</given-names></name><name><surname>Gordon</surname><given-names>JA</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Somatostatin interneurons facilitate hippocampal-prefrontal synchrony and prefrontal spatial encoding</article-title><source>Neuron</source><volume>100</volume><fpage>926</fpage><lpage>939</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2018.09.029</pub-id><pub-id pub-id-type="pmid">30318409</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aston-Jones</surname><given-names>G</given-names></name><name><surname>Cohen</surname><given-names>JD</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>An integrative theory of locus coeruleus-norepinephrine function: adaptive gain and optimal performance</article-title><source>Annual Review of Neuroscience</source><volume>28</volume><fpage>403</fpage><lpage>450</lpage><pub-id pub-id-type="doi">10.1146/annurev.neuro.28.061604.135709</pub-id><pub-id pub-id-type="pmid">16022602</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Beatty</surname><given-names>J</given-names></name><name><surname>Lucero-Wagoner</surname><given-names>B</given-names></name></person-group><year iso-8601-date="2000">2000</year><chapter-title>The pupillary system</chapter-title><person-group person-group-type="editor"><name><surname>Beatty</surname><given-names>J</given-names></name></person-group><source>Handbook of Psychophysiology</source><publisher-name>Cambridge University Press</publisher-name><fpage>142</fpage><lpage>162</lpage></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beeler</surname><given-names>GW</given-names></name></person-group><year iso-8601-date="1967">1967</year><article-title>Visual threshold changes resulting from spontaneous saccadic eye movements</article-title><source>Vision Research</source><volume>7</volume><fpage>769</fpage><lpage>775</lpage><pub-id pub-id-type="doi">10.1016/0042-6989(67)90039-9</pub-id><pub-id pub-id-type="pmid">5608665</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buffalo</surname><given-names>EA</given-names></name><name><surname>Fries</surname><given-names>P</given-names></name><name><surname>Landman</surname><given-names>R</given-names></name><name><surname>Buschman</surname><given-names>TJ</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Laminar differences in gamma and alpha coherence in the ventral stream</article-title><source>PNAS</source><volume>108</volume><fpage>11262</fpage><lpage>11267</lpage><pub-id pub-id-type="doi">10.1073/pnas.1011284108</pub-id><pub-id pub-id-type="pmid">21690410</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Calin-Jageman</surname><given-names>RJ</given-names></name><name><surname>Cumming</surname><given-names>G</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Estimation for better inference in neuroscience</article-title><source>eNeuro</source><volume>6</volume><elocation-id>ENEURO.0205-19.2019</elocation-id><pub-id pub-id-type="doi">10.1523/ENEURO.0205-19.2019</pub-id><pub-id pub-id-type="pmid">31453316</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen</surname><given-names>MR</given-names></name><name><surname>Maunsell</surname><given-names>JHR</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Attention improves performance primarily by reducing interneuronal correlations</article-title><source>Nature Neuroscience</source><volume>12</volume><fpage>1594</fpage><lpage>1600</lpage><pub-id pub-id-type="doi">10.1038/nn.2439</pub-id><pub-id pub-id-type="pmid">19915566</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Connors</surname><given-names>BW</given-names></name><name><surname>Gutnick</surname><given-names>MJ</given-names></name></person-group><year iso-8601-date="1990">1990</year><article-title>Intrinsic firing patterns of diverse neocortical neurons</article-title><source>Trends in Neurosciences</source><volume>13</volume><fpage>99</fpage><lpage>104</lpage><pub-id pub-id-type="doi">10.1016/0166-2236(90)90185-d</pub-id><pub-id pub-id-type="pmid">1691879</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Contreras</surname><given-names>D</given-names></name><name><surname>Palmer</surname><given-names>L</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>Response to contrast of electrophysiologically defined cell classes in primary visual cortex</article-title><source>The Journal of Neuroscience</source><volume>23</volume><fpage>6936</fpage><lpage>6945</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.23-17-06936.2003</pub-id><pub-id pub-id-type="pmid">12890788</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cools</surname><given-names>R</given-names></name><name><surname>D’Esposito</surname><given-names>M</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Inverted-U-shaped dopamine actions on human working memory and cognitive control</article-title><source>Biological Psychiatry</source><volume>69</volume><fpage>e113</fpage><lpage>e125</lpage><pub-id pub-id-type="doi">10.1016/j.biopsych.2011.03.028</pub-id><pub-id pub-id-type="pmid">21531388</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Davis</surname><given-names>ZW</given-names></name><name><surname>Muller</surname><given-names>L</given-names></name><name><surname>Martinez-Trujillo</surname><given-names>J</given-names></name><name><surname>Sejnowski</surname><given-names>T</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="2020">2020</year><article-title>Spontaneous travelling cortical waves gate perception in behaving primates</article-title><source>Nature</source><volume>587</volume><fpage>432</fpage><lpage>436</lpage><pub-id pub-id-type="doi">10.1038/s41586-020-2802-y</pub-id><pub-id pub-id-type="pmid">33029013</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dehaene</surname><given-names>S</given-names></name><name><surname>Changeux</surname><given-names>J-P</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Experimental and theoretical approaches to conscious processing</article-title><source>Neuron</source><volume>70</volume><fpage>200</fpage><lpage>227</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2011.03.018</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Desimone</surname><given-names>R</given-names></name><name><surname>Duncan</surname><given-names>J</given-names></name></person-group><year iso-8601-date="1995">1995</year><article-title>Neural mechanisms of selective visual attention</article-title><source>Annual Review of Neuroscience</source><volume>18</volume><fpage>193</fpage><lpage>222</lpage><pub-id pub-id-type="doi">10.1146/annurev.ne.18.030195.001205</pub-id><pub-id pub-id-type="pmid">7605061</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dicke</surname><given-names>PW</given-names></name><name><surname>Chakraborty</surname><given-names>S</given-names></name><name><surname>Thier</surname><given-names>P</given-names></name></person-group><year iso-8601-date="2008">2008</year><article-title>Neuronal correlates of perceptual stability during eye movements</article-title><source>The European Journal of Neuroscience</source><volume>27</volume><fpage>991</fpage><lpage>1002</lpage><pub-id pub-id-type="doi">10.1111/j.1460-9568.2008.06054.x</pub-id><pub-id pub-id-type="pmid">18333969</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Douglas</surname><given-names>RJ</given-names></name><name><surname>Martin</surname><given-names>KAC</given-names></name></person-group><year iso-8601-date="2004">2004</year><article-title>Neuronal circuits of the neocortex</article-title><source>Annual Review of Neuroscience</source><volume>27</volume><fpage>419</fpage><lpage>451</lpage><pub-id pub-id-type="doi">10.1146/annurev.neuro.27.070203.144152</pub-id><pub-id pub-id-type="pmid">15217339</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Douglas</surname><given-names>RJ</given-names></name><name><surname>Martin</surname><given-names>KAC</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Mapping the matrix: the ways of neocortex</article-title><source>Neuron</source><volume>56</volume><fpage>226</fpage><lpage>238</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2007.10.017</pub-id><pub-id pub-id-type="pmid">17964242</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Efron</surname><given-names>B</given-names></name><name><surname>Tibshirani</surname><given-names>RJ</given-names></name></person-group><year iso-8601-date="1993">1993</year><chapter-title>An introduction to the bootstrap</chapter-title><person-group person-group-type="editor"><name><surname>Efron</surname><given-names>B</given-names></name></person-group><source>Champman and Hall/CRC</source><publisher-name>Springer</publisher-name><fpage>1</fpage><lpage>10</lpage><pub-id pub-id-type="doi">10.1007/978-1-4899-4541-9</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Engbert</surname><given-names>R</given-names></name><name><surname>Kliegl</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>Microsaccades uncover the orientation of covert attention</article-title><source>Vision Research</source><volume>43</volume><fpage>1035</fpage><lpage>1045</lpage><pub-id pub-id-type="doi">10.1016/s0042-6989(03)00084-1</pub-id><pub-id pub-id-type="pmid">12676246</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Engel</surname><given-names>TA</given-names></name><name><surname>Steinmetz</surname><given-names>NA</given-names></name><name><surname>Gieselmann</surname><given-names>MA</given-names></name><name><surname>Thiele</surname><given-names>A</given-names></name><name><surname>Moore</surname><given-names>T</given-names></name><name><surname>Boahen</surname><given-names>K</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>Selective modulation of cortical state during spatial attention</article-title><source>Science</source><volume>354</volume><fpage>1140</fpage><lpage>1144</lpage><pub-id pub-id-type="doi">10.1126/science.aag1420</pub-id><pub-id pub-id-type="pmid">27934763</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ferro</surname><given-names>D</given-names></name><name><surname>van Kempen</surname><given-names>J</given-names></name><name><surname>Boyd</surname><given-names>M</given-names></name><name><surname>Panzeri</surname><given-names>S</given-names></name><name><surname>Thiele</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2021">2021</year><article-title>Directed information exchange between cortical layers in macaque V1 and V4 and its modulation by selective attention</article-title><source>PNAS</source><volume>118</volume><elocation-id>e2022097118</elocation-id><pub-id pub-id-type="doi">10.1073/pnas.2022097118</pub-id><pub-id pub-id-type="pmid">33723059</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Foehring</surname><given-names>RC</given-names></name><name><surname>Lorenzon</surname><given-names>NM</given-names></name><name><surname>Herron</surname><given-names>P</given-names></name><name><surname>Wilson</surname><given-names>CJ</given-names></name></person-group><year iso-8601-date="1991">1991</year><article-title>Correlation of physiologically and morphologically identified neuronal types in human association cortex in vitro</article-title><source>Journal of Neurophysiology</source><volume>66</volume><fpage>1825</fpage><lpage>1837</lpage><pub-id pub-id-type="doi">10.1152/jn.1991.66.6.1825</pub-id><pub-id pub-id-type="pmid">1812219</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fries</surname><given-names>P</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name><name><surname>Rorie</surname><given-names>AE</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2001">2001</year><article-title>Modulation of oscillatory neuronal synchronization by selective visual attention</article-title><source>Science</source><volume>291</volume><fpage>1560</fpage><lpage>1563</lpage><pub-id pub-id-type="doi">10.1126/science.1055465</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fries</surname><given-names>P</given-names></name><name><surname>Womelsdorf</surname><given-names>T</given-names></name><name><surname>Oostenveld</surname><given-names>R</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2008">2008</year><article-title>The effects of visual stimulation and selective visual attention on rhythmic neuronal synchronization in macaque area V4</article-title><source>The Journal of Neuroscience</source><volume>28</volume><fpage>4823</fpage><lpage>4835</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.4499-07.2008</pub-id><pub-id pub-id-type="pmid">18448659</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fries</surname><given-names>P</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Neuronal gamma-band synchronization as a fundamental process in cortical computation</article-title><source>Annual Review of Neuroscience</source><volume>32</volume><fpage>209</fpage><lpage>224</lpage><pub-id pub-id-type="doi">10.1146/annurev.neuro.051508.135603</pub-id><pub-id pub-id-type="pmid">19400723</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Galfano</surname><given-names>G</given-names></name><name><surname>Betta</surname><given-names>E</given-names></name><name><surname>Turatto</surname><given-names>M</given-names></name></person-group><year iso-8601-date="2004">2004</year><article-title>Inhibition of return in microsaccades</article-title><source>Experimental Brain Research</source><volume>159</volume><fpage>400</fpage><lpage>404</lpage><pub-id pub-id-type="doi">10.1007/s00221-004-2111-y</pub-id><pub-id pub-id-type="pmid">15480591</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghosh</surname><given-names>S</given-names></name><name><surname>Maunsell</surname><given-names>JHR</given-names></name></person-group><year iso-8601-date="2021">2021</year><article-title>Single trial neuronal activity dynamics of attentional intensity in monkey visual area V4</article-title><source>Nature Communications</source><volume>12</volume><elocation-id>2003</elocation-id><pub-id pub-id-type="doi">10.1038/s41467-021-22281-2</pub-id><pub-id pub-id-type="pmid">33790282</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodale</surname><given-names>MA</given-names></name><name><surname>Milner</surname><given-names>AD</given-names></name></person-group><year iso-8601-date="1992">1992</year><article-title>Separate visual pathways for perception and action</article-title><source>Trends in Neurosciences</source><volume>15</volume><fpage>20</fpage><lpage>25</lpage><pub-id pub-id-type="doi">10.1016/0166-2236(92)90344-8</pub-id><pub-id pub-id-type="pmid">1374953</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gowen</surname><given-names>E</given-names></name><name><surname>Abadi</surname><given-names>RV</given-names></name><name><surname>Poliakoff</surname><given-names>E</given-names></name><name><surname>Hansen</surname><given-names>PC</given-names></name><name><surname>Miall</surname><given-names>RC</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Modulation of saccadic intrusions by exogenous and endogenous attention</article-title><source>Brain Research</source><volume>1141</volume><fpage>154</fpage><lpage>167</lpage><pub-id pub-id-type="doi">10.1016/j.brainres.2007.01.047</pub-id><pub-id pub-id-type="pmid">17313940</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hafed</surname><given-names>ZM</given-names></name><name><surname>Clark</surname><given-names>JJ</given-names></name></person-group><year iso-8601-date="2002">2002</year><article-title>Microsaccades as an overt measure of covert attention shifts</article-title><source>Vision Research</source><volume>42</volume><fpage>2533</fpage><lpage>2545</lpage><pub-id pub-id-type="doi">10.1016/s0042-6989(02)00263-8</pub-id><pub-id pub-id-type="pmid">12445847</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hafed</surname><given-names>ZM</given-names></name><name><surname>Krauzlis</surname><given-names>RJ</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>Microsaccadic suppression of visual bursts in the primate superior colliculus</article-title><source>The Journal of Neuroscience</source><volume>30</volume><fpage>9542</fpage><lpage>9547</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.1137-10.2010</pub-id><pub-id pub-id-type="pmid">20631182</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Harris</surname><given-names>KD</given-names></name><name><surname>Thiele</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Cortical state and attention</article-title><source>Nature Reviews. Neuroscience</source><volume>12</volume><fpage>509</fpage><lpage>523</lpage><pub-id pub-id-type="doi">10.1038/nrn3084</pub-id><pub-id pub-id-type="pmid">21829219</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hartigan</surname><given-names>JA</given-names></name><name><surname>Hartigan</surname><given-names>PM</given-names></name></person-group><year iso-8601-date="1985">1985</year><article-title>The dip test of unimodality</article-title><source>The Annals of Statistics</source><volume>13</volume><elocation-id>1176346577</elocation-id><pub-id pub-id-type="doi">10.1214/aos/1176346577</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hasenstaub</surname><given-names>A</given-names></name><name><surname>Shu</surname><given-names>Y</given-names></name><name><surname>Haider</surname><given-names>B</given-names></name><name><surname>Kraushaar</surname><given-names>U</given-names></name><name><surname>Duque</surname><given-names>A</given-names></name><name><surname>McCormick</surname><given-names>DA</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>Inhibitory postsynaptic potentials carry synchronized frequency information in active cortical networks</article-title><source>Neuron</source><volume>47</volume><fpage>423</fpage><lpage>435</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2005.06.016</pub-id><pub-id pub-id-type="pmid">16055065</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Herman</surname><given-names>WX</given-names></name><name><surname>Smith</surname><given-names>RE</given-names></name><name><surname>Kronemer</surname><given-names>SI</given-names></name><name><surname>Watsky</surname><given-names>RE</given-names></name><name><surname>Chen</surname><given-names>WC</given-names></name><name><surname>Gober</surname><given-names>LM</given-names></name><name><surname>Touloumes</surname><given-names>GJ</given-names></name><name><surname>Khosla</surname><given-names>M</given-names></name><name><surname>Raja</surname><given-names>A</given-names></name><name><surname>Horien</surname><given-names>CL</given-names></name><name><surname>Morse</surname><given-names>EC</given-names></name><name><surname>Botta</surname><given-names>KL</given-names></name><name><surname>Hirsch</surname><given-names>LJ</given-names></name><name><surname>Alkawadri</surname><given-names>R</given-names></name><name><surname>Gerrard</surname><given-names>JL</given-names></name><name><surname>Spencer</surname><given-names>DD</given-names></name><name><surname>Blumenfeld</surname><given-names>H</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>A switch and wave of neuronal activity in the cerebral cortex during the first second of conscious perception</article-title><source>Cerebral Cortex</source><volume>29</volume><fpage>461</fpage><lpage>474</lpage><pub-id pub-id-type="doi">10.1093/cercor/bhx327</pub-id><pub-id pub-id-type="pmid">29194517</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hess</surname><given-names>EH</given-names></name><name><surname>Polt</surname><given-names>JM</given-names></name></person-group><year iso-8601-date="1964">1964</year><article-title>Pupil size in relation to mental activity during simple problem-solving</article-title><source>Science</source><volume>143</volume><fpage>1190</fpage><lpage>1192</lpage><pub-id pub-id-type="doi">10.1126/science.143.3611.1190</pub-id><pub-id pub-id-type="pmid">17833905</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hirsch</surname><given-names>JA</given-names></name><name><surname>Martinez</surname><given-names>LM</given-names></name></person-group><year iso-8601-date="2006">2006</year><article-title>Laminar processing in the visual cortical column</article-title><source>Current Opinion in Neurobiology</source><volume>16</volume><fpage>377</fpage><lpage>384</lpage><pub-id pub-id-type="doi">10.1016/j.conb.2006.06.014</pub-id><pub-id pub-id-type="pmid">16842989</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ho</surname><given-names>J</given-names></name><name><surname>Tumkaya</surname><given-names>T</given-names></name><name><surname>Aryal</surname><given-names>S</given-names></name><name><surname>Choi</surname><given-names>H</given-names></name><name><surname>Claridge-Chang</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Moving beyond P values: data analysis with estimation graphics</article-title><source>Nature Methods</source><volume>16</volume><fpage>565</fpage><lpage>566</lpage><pub-id pub-id-type="doi">10.1038/s41592-019-0470-3</pub-id><pub-id pub-id-type="pmid">31217592</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hoffmann</surname><given-names>H</given-names></name></person-group><year iso-8601-date="2015">2015</year><source>Violin.m - Simple Violin Plot Using Matlab Default Kernel Density Estimation</source><publisher-name>University of Bonn</publisher-name></element-citation></ref><ref id="bib39"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hotelling</surname><given-names>H</given-names></name></person-group><year iso-8601-date="1992">1992</year><chapter-title>Relations between two sets of variates</chapter-title><person-group person-group-type="editor"><name><surname>Kotz</surname><given-names>S</given-names></name><name><surname>Johnson</surname><given-names>NL</given-names></name></person-group><source>Breakthroughs in Statistics: Methodology and Distribution</source><publisher-name>Springer</publisher-name><fpage>162</fpage><lpage>190</lpage><pub-id pub-id-type="doi">10.1007/978-1-4612-4380-9_14</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kawaguchi</surname><given-names>Y</given-names></name></person-group><year iso-8601-date="1993">1993</year><article-title>Groupings of nonpyramidal and pyramidal cells with specific physiological and morphological characteristics in rat frontal cortex</article-title><source>Journal of Neurophysiology</source><volume>69</volume><fpage>416</fpage><lpage>431</lpage><pub-id pub-id-type="doi">10.1152/jn.1993.69.2.416</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>König</surname><given-names>SD</given-names></name><name><surname>Buffalo</surname><given-names>EA</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>A nonparametric method for detecting fixations and saccades using cluster analysis: removing the need for arbitrary thresholds</article-title><source>Journal of Neuroscience Methods</source><volume>227</volume><fpage>121</fpage><lpage>131</lpage><pub-id pub-id-type="doi">10.1016/j.jneumeth.2014.01.032</pub-id><pub-id pub-id-type="pmid">24509130</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Krosigk von</surname><given-names>M</given-names></name><name><surname>Bal</surname><given-names>T</given-names></name><name><surname>McCormick</surname><given-names>DA</given-names></name></person-group><year iso-8601-date="1993">1993</year><article-title>Cellular mechanisms of a synchronized oscillation in the thalamus</article-title><source>Science</source><volume>261</volume><fpage>361</fpage><lpage>364</lpage><pub-id pub-id-type="doi">10.1126/science.8392750</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Levitt</surname><given-names>H</given-names></name></person-group><year iso-8601-date="1971">1971</year><article-title>Transformed up-down methods in psychoacoustics</article-title><source>The Journal of the Acoustical Society of America</source><volume>49</volume><fpage>467</fpage><lpage>477</lpage><pub-id pub-id-type="pmid">5541744</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Livingstone</surname><given-names>MS</given-names></name><name><surname>Hubel</surname><given-names>DH</given-names></name></person-group><year iso-8601-date="1981">1981</year><article-title>Effects of sleep and arousal on the processing of visual information in the cat</article-title><source>Nature</source><volume>291</volume><fpage>554</fpage><lpage>561</lpage><pub-id pub-id-type="doi">10.1038/291554a0</pub-id><pub-id pub-id-type="pmid">6165893</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowet</surname><given-names>E</given-names></name><name><surname>Gomes</surname><given-names>B</given-names></name><name><surname>Srinivasan</surname><given-names>K</given-names></name><name><surname>Zhou</surname><given-names>H</given-names></name><name><surname>Schafer</surname><given-names>RJ</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Enhanced neural processing by covert attention only during microsaccades directed toward the attended stimulus</article-title><source>Neuron</source><volume>99</volume><fpage>207</fpage><lpage>214</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2018.05.041</pub-id><pub-id pub-id-type="pmid">29937279</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Markram</surname><given-names>H</given-names></name><name><surname>Toledo-Rodriguez</surname><given-names>M</given-names></name><name><surname>Wang</surname><given-names>Y</given-names></name><name><surname>Gupta</surname><given-names>A</given-names></name><name><surname>Silberberg</surname><given-names>G</given-names></name><name><surname>Wu</surname><given-names>C</given-names></name></person-group><year iso-8601-date="2004">2004</year><article-title>Interneurons of the neocortical inhibitory system</article-title><source>Nature Reviews. Neuroscience</source><volume>5</volume><fpage>793</fpage><lpage>807</lpage><pub-id pub-id-type="doi">10.1038/nrn1519</pub-id><pub-id pub-id-type="pmid">15378039</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marshel</surname><given-names>JH</given-names></name><name><surname>Kim</surname><given-names>YS</given-names></name><name><surname>Machado</surname><given-names>TA</given-names></name><name><surname>Quirin</surname><given-names>S</given-names></name><name><surname>Benson</surname><given-names>B</given-names></name><name><surname>Kadmon</surname><given-names>J</given-names></name><name><surname>Raja</surname><given-names>C</given-names></name><name><surname>Chibukhchyan</surname><given-names>A</given-names></name><name><surname>Ramakrishnan</surname><given-names>C</given-names></name><name><surname>Inoue</surname><given-names>M</given-names></name><name><surname>Shane</surname><given-names>JC</given-names></name><name><surname>McKnight</surname><given-names>DJ</given-names></name><name><surname>Yoshizawa</surname><given-names>S</given-names></name><name><surname>Kato</surname><given-names>HE</given-names></name><name><surname>Ganguli</surname><given-names>S</given-names></name><name><surname>Deisseroth</surname><given-names>K</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Cortical layer-specific critical dynamics triggering perception</article-title><source>Science</source><volume>365</volume><elocation-id>eaaw5202</elocation-id><pub-id pub-id-type="doi">10.1126/science.aaw5202</pub-id><pub-id pub-id-type="pmid">31320556</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martinez-Conde</surname><given-names>S</given-names></name><name><surname>Otero-Millan</surname><given-names>J</given-names></name><name><surname>Macknik</surname><given-names>SL</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>The impact of microsaccades on vision: towards a unified theory of saccadic function</article-title><source>Nature Reviews. Neuroscience</source><volume>14</volume><fpage>83</fpage><lpage>96</lpage><pub-id pub-id-type="doi">10.1038/nrn3405</pub-id><pub-id pub-id-type="pmid">23329159</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McAdams</surname><given-names>CJ</given-names></name><name><surname>Maunsell</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="1999">1999</year><article-title>Effects of attention on orientation-tuning functions of single neurons in macaque cortical area V4</article-title><source>The Journal of Neuroscience</source><volume>19</volume><fpage>431</fpage><lpage>441</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.19-01-00431.1999</pub-id><pub-id pub-id-type="pmid">9870971</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McCormick</surname><given-names>DA</given-names></name><name><surname>Connors</surname><given-names>BW</given-names></name><name><surname>Lighthall</surname><given-names>JW</given-names></name><name><surname>Prince</surname><given-names>DA</given-names></name></person-group><year iso-8601-date="1985">1985</year><article-title>Comparative electrophysiology of pyramidal and sparsely spiny stellate neurons of the neocortex</article-title><source>Journal of Neurophysiology</source><volume>54</volume><fpage>782</fpage><lpage>806</lpage><pub-id pub-id-type="doi">10.1152/jn.1985.54.4.782</pub-id><pub-id pub-id-type="pmid">2999347</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McCormick</surname><given-names>DA</given-names></name><name><surname>Bal</surname><given-names>T</given-names></name></person-group><year iso-8601-date="1997">1997</year><article-title>Sleep and arousal: thalamocortical mechanisms</article-title><source>Annual Review of Neuroscience</source><volume>20</volume><fpage>185</fpage><lpage>215</lpage><pub-id pub-id-type="doi">10.1146/annurev.neuro.20.1.185</pub-id><pub-id pub-id-type="pmid">9056712</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McCormick</surname><given-names>DA</given-names></name><name><surname>Nestvogel</surname><given-names>DB</given-names></name><name><surname>He</surname><given-names>BJ</given-names></name></person-group><year iso-8601-date="2020">2020</year><article-title>Neuromodulation of brain state and behavior</article-title><source>Annual Review of Neuroscience</source><volume>43</volume><fpage>391</fpage><lpage>415</lpage><pub-id pub-id-type="doi">10.1146/annurev-neuro-100219-105424</pub-id><pub-id pub-id-type="pmid">32250724</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McDougal</surname><given-names>DH</given-names></name><name><surname>Gamlin</surname><given-names>PD</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Autonomic control of the eye</article-title><source>Comprehensive Physiology</source><volume>5</volume><fpage>439</fpage><lpage>473</lpage><pub-id pub-id-type="doi">10.1002/cphy.c140014</pub-id><pub-id pub-id-type="pmid">25589275</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McGinley</surname><given-names>MJ</given-names></name><name><surname>David</surname><given-names>SV</given-names></name><name><surname>McCormick</surname><given-names>DA</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Cortical membrane potential signature of optimal states for sensory signal detection</article-title><source>Neuron</source><volume>87</volume><fpage>179</fpage><lpage>192</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2015.05.038</pub-id><pub-id pub-id-type="pmid">26074005</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mehta</surname><given-names>AD</given-names></name><name><surname>Ulbert</surname><given-names>I</given-names></name><name><surname>Schroeder</surname><given-names>CE</given-names></name></person-group><year iso-8601-date="2000">2000a</year><article-title>Intermodal selective attention in monkeys. I: distribution and timing of effects across visual areas</article-title><source>Cerebral Cortex</source><volume>10</volume><fpage>343</fpage><lpage>358</lpage><pub-id pub-id-type="doi">10.1093/cercor/10.4.343</pub-id><pub-id pub-id-type="pmid">10769247</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mehta</surname><given-names>AD</given-names></name><name><surname>Ulbert</surname><given-names>I</given-names></name><name><surname>Schroeder</surname><given-names>CE</given-names></name></person-group><year iso-8601-date="2000">2000b</year><article-title>Intermodal selective attention in monkeys. II: physiological mechanisms of modulation</article-title><source>Cerebral Cortex</source><volume>10</volume><fpage>359</fpage><lpage>370</lpage><pub-id pub-id-type="doi">10.1093/cercor/10.4.359</pub-id><pub-id pub-id-type="pmid">10769248</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Migliore</surname><given-names>M</given-names></name><name><surname>Shepherd</surname><given-names>GM</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>Opinion: an integrated approach to classifying neuronal phenotypes</article-title><source>Nature Reviews. Neuroscience</source><volume>6</volume><fpage>810</fpage><lpage>818</lpage><pub-id pub-id-type="doi">10.1038/nrn1769</pub-id><pub-id pub-id-type="pmid">16276357</pub-id></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mitchell</surname><given-names>JF</given-names></name><name><surname>Sundberg</surname><given-names>KA</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Differential attention-dependent response modulation across cell classes in macaque visual area V4</article-title><source>Neuron</source><volume>55</volume><fpage>131</fpage><lpage>141</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2007.06.018</pub-id><pub-id pub-id-type="pmid">17610822</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mitchell</surname><given-names>JF</given-names></name><name><surname>Sundberg</surname><given-names>KA</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Spatial attention decorrelates intrinsic activity fluctuations in macaque area V4</article-title><source>Neuron</source><volume>63</volume><fpage>879</fpage><lpage>888</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2009.09.013</pub-id><pub-id pub-id-type="pmid">19778515</pub-id></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mitra</surname><given-names>PP</given-names></name><name><surname>Pesaran</surname><given-names>B</given-names></name></person-group><year iso-8601-date="1999">1999</year><article-title>Analysis of dynamic brain imaging data</article-title><source>Biophysical Journal</source><volume>76</volume><fpage>691</fpage><lpage>708</lpage><pub-id pub-id-type="doi">10.1016/S0006-3495(99)77236-X</pub-id><pub-id pub-id-type="pmid">9929474</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Mitra</surname><given-names>P</given-names></name></person-group><year iso-8601-date="2007">2007</year><source>Observed Brain Dynamics</source><publisher-name>Oxford University Press</publisher-name></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mitzdorf</surname><given-names>U</given-names></name></person-group><year iso-8601-date="1985">1985</year><article-title>Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena</article-title><source>Physiological Reviews</source><volume>65</volume><fpage>37</fpage><lpage>100</lpage><pub-id pub-id-type="doi">10.1152/physrev.1985.65.1.37</pub-id><pub-id pub-id-type="pmid">3880898</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moore</surname><given-names>T</given-names></name><name><surname>Fallah</surname><given-names>M</given-names></name></person-group><year iso-8601-date="2001">2001</year><article-title>Control of eye movements and spatial attention</article-title><source>PNAS</source><volume>98</volume><fpage>1273</fpage><lpage>1276</lpage><pub-id pub-id-type="doi">10.1073/pnas.021549498</pub-id></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moore</surname><given-names>T</given-names></name><name><surname>Zirnsak</surname><given-names>M</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Neural mechanisms of selective visual attention</article-title><source>Annual Review of Psychology</source><volume>68</volume><fpage>47</fpage><lpage>72</lpage><pub-id pub-id-type="doi">10.1146/annurev-psych-122414-033400</pub-id><pub-id pub-id-type="pmid">28051934</pub-id></element-citation></ref><ref id="bib65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moran</surname><given-names>J</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="1985">1985</year><article-title>Selective attention gates visual processing in the extrastriate cortex</article-title><source>Science</source><volume>229</volume><fpage>782</fpage><lpage>784</lpage><pub-id pub-id-type="doi">10.1126/science.4023713</pub-id><pub-id pub-id-type="pmid">4023713</pub-id></element-citation></ref><ref id="bib66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moreno-Bote</surname><given-names>R</given-names></name><name><surname>Beck</surname><given-names>J</given-names></name><name><surname>Kanitscheider</surname><given-names>I</given-names></name><name><surname>Pitkow</surname><given-names>X</given-names></name><name><surname>Latham</surname><given-names>P</given-names></name><name><surname>Pouget</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Information-limiting correlations</article-title><source>Nature Neuroscience</source><volume>17</volume><fpage>1410</fpage><lpage>1417</lpage><pub-id pub-id-type="doi">10.1038/nn.3807</pub-id><pub-id pub-id-type="pmid">25195105</pub-id></element-citation></ref><ref id="bib67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mountcastle</surname><given-names>VB</given-names></name></person-group><year iso-8601-date="1997">1997</year><article-title>The columnar organization of the neocortex</article-title><source>Brain</source><volume>120 ( Pt 4)</volume><fpage>701</fpage><lpage>722</lpage><pub-id pub-id-type="doi">10.1093/brain/120.4.701</pub-id><pub-id pub-id-type="pmid">9153131</pub-id></element-citation></ref><ref id="bib68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Murphy</surname><given-names>PR</given-names></name><name><surname>Robertson</surname><given-names>IH</given-names></name><name><surname>Balsters</surname><given-names>JH</given-names></name><name><surname>O’connell</surname><given-names>RG</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Pupillometry and P3 index the locus coeruleus-noradrenergic arousal function in humans</article-title><source>Psychophysiology</source><volume>48</volume><fpage>1532</fpage><lpage>1543</lpage><pub-id pub-id-type="doi">10.1111/j.1469-8986.2011.01226.x</pub-id><pub-id pub-id-type="pmid">21762458</pub-id></element-citation></ref><ref id="bib69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nandy</surname><given-names>AS</given-names></name><name><surname>Nassi</surname><given-names>JJ</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Laminar organization of attentional modulation in macaque visual area V4</article-title><source>Neuron</source><volume>93</volume><fpage>235</fpage><lpage>246</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2016.11.029</pub-id><pub-id pub-id-type="pmid">27989456</pub-id></element-citation></ref><ref id="bib70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nassi</surname><given-names>JJ</given-names></name><name><surname>Avery</surname><given-names>MC</given-names></name><name><surname>Cetin</surname><given-names>AH</given-names></name><name><surname>Roe</surname><given-names>AW</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Optogenetic activation of normalization in alert macaque visual cortex</article-title><source>Neuron</source><volume>86</volume><fpage>1504</fpage><lpage>1517</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2015.05.040</pub-id><pub-id pub-id-type="pmid">26087167</pub-id></element-citation></ref><ref id="bib71"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nowak</surname><given-names>LG</given-names></name><name><surname>Azouz</surname><given-names>R</given-names></name><name><surname>Sanchez-Vives</surname><given-names>MV</given-names></name><name><surname>Gray</surname><given-names>CM</given-names></name><name><surname>McCormick</surname><given-names>DA</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>Electrophysiological classes of cat primary visual cortical neurons in vivo as revealed by quantitative analyses</article-title><source>Journal of Neurophysiology</source><volume>89</volume><fpage>1541</fpage><lpage>1566</lpage><pub-id pub-id-type="doi">10.1152/jn.00580.2002</pub-id><pub-id pub-id-type="pmid">12626627</pub-id></element-citation></ref><ref id="bib72"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pastukhov</surname><given-names>A</given-names></name><name><surname>Braun</surname><given-names>J</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>Rare but precious: microsaccades are highly informative about attentional allocation</article-title><source>Vision Research</source><volume>50</volume><fpage>1173</fpage><lpage>1184</lpage><pub-id pub-id-type="doi">10.1016/j.visres.2010.04.007</pub-id><pub-id pub-id-type="pmid">20382176</pub-id></element-citation></ref><ref id="bib73"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pettine</surname><given-names>WW</given-names></name><name><surname>Steinmetz</surname><given-names>NA</given-names></name><name><surname>Moore</surname><given-names>T</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Laminar segregation of sensory coding and behavioral readout in macaque V4</article-title><source>PNAS</source><volume>116</volume><fpage>14749</fpage><lpage>14754</lpage><pub-id pub-id-type="doi">10.1073/pnas.1819398116</pub-id><pub-id pub-id-type="pmid">31249141</pub-id></element-citation></ref><ref id="bib74"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pins</surname><given-names>D</given-names></name><name><surname>Ffytche</surname><given-names>D</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>The neural correlates of conscious vision</article-title><source>Cerebral Cortex</source><volume>13</volume><fpage>461</fpage><lpage>474</lpage><pub-id pub-id-type="doi">10.1093/cercor/13.5.461</pub-id><pub-id pub-id-type="pmid">12679293</pub-id></element-citation></ref><ref id="bib75"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Povysheva</surname><given-names>NV</given-names></name><name><surname>Gonzalez-Burgos</surname><given-names>G</given-names></name><name><surname>Zaitsev</surname><given-names>AV</given-names></name><name><surname>Kröner</surname><given-names>S</given-names></name><name><surname>Barrionuevo</surname><given-names>G</given-names></name><name><surname>Lewis</surname><given-names>DA</given-names></name><name><surname>Krimer</surname><given-names>LS</given-names></name></person-group><year iso-8601-date="2006">2006</year><article-title>Properties of excitatory synaptic responses in fast-spiking interneurons and pyramidal cells from monkey and rat prefrontal cortex</article-title><source>Cerebral Cortex</source><volume>16</volume><fpage>541</fpage><lpage>552</lpage><pub-id pub-id-type="doi">10.1093/cercor/bhj002</pub-id><pub-id pub-id-type="pmid">16033926</pub-id></element-citation></ref><ref id="bib76"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Prins</surname><given-names>N</given-names></name><name><surname>Kingdom</surname><given-names>FAA</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Applying the model-comparison approach to test specific research hypotheses in psychophysical research using the palamedes toolbox</article-title><source>Frontiers in Psychology</source><volume>9</volume><elocation-id>1250</elocation-id><pub-id pub-id-type="doi">10.3389/fpsyg.2018.01250</pub-id><pub-id pub-id-type="pmid">30083122</pub-id></element-citation></ref><ref id="bib77"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rajagovindan</surname><given-names>R</given-names></name><name><surname>Ding</surname><given-names>M</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>From prestimulus alpha oscillation to visual-evoked response: an inverted-U function and its attentional modulation</article-title><source>Journal of Cognitive Neuroscience</source><volume>23</volume><fpage>1379</fpage><lpage>1394</lpage><pub-id pub-id-type="doi">10.1162/jocn.2010.21478</pub-id><pub-id pub-id-type="pmid">20459310</pub-id></element-citation></ref><ref id="bib78"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reimer</surname><given-names>J</given-names></name><name><surname>Froudarakis</surname><given-names>E</given-names></name><name><surname>Cadwell</surname><given-names>CR</given-names></name><name><surname>Yatsenko</surname><given-names>D</given-names></name><name><surname>Denfield</surname><given-names>GH</given-names></name><name><surname>Tolias</surname><given-names>AS</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Pupil fluctuations track fast switching of cortical states during quiet wakefulness</article-title><source>Neuron</source><volume>84</volume><fpage>355</fpage><lpage>362</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2014.09.033</pub-id><pub-id pub-id-type="pmid">25374359</pub-id></element-citation></ref><ref id="bib79"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ress</surname><given-names>D</given-names></name><name><surname>Heeger</surname><given-names>DJ</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>Neuronal correlates of perception in early visual cortex</article-title><source>Nature Neuroscience</source><volume>6</volume><fpage>414</fpage><lpage>420</lpage><pub-id pub-id-type="doi">10.1038/nn1024</pub-id><pub-id pub-id-type="pmid">12627164</pub-id></element-citation></ref><ref id="bib80"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reynolds</surname><given-names>JH</given-names></name><name><surname>Pasternak</surname><given-names>T</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name></person-group><year iso-8601-date="2000">2000</year><article-title>Attention increases sensitivity of V4 neurons</article-title><source>Neuron</source><volume>26</volume><fpage>703</fpage><lpage>714</lpage><pub-id pub-id-type="doi">10.1016/s0896-6273(00)81206-4</pub-id><pub-id pub-id-type="pmid">10896165</pub-id></element-citation></ref><ref id="bib81"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Roe</surname><given-names>AW</given-names></name><name><surname>Chelazzi</surname><given-names>L</given-names></name><name><surname>Connor</surname><given-names>CE</given-names></name><name><surname>Conway</surname><given-names>BR</given-names></name><name><surname>Fujita</surname><given-names>I</given-names></name><name><surname>Gallant</surname><given-names>JL</given-names></name><name><surname>Lu</surname><given-names>H</given-names></name><name><surname>Vanduffel</surname><given-names>W</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Toward a unified theory of visual area V4</article-title><source>Neuron</source><volume>74</volume><fpage>12</fpage><lpage>29</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2012.03.011</pub-id><pub-id pub-id-type="pmid">22500626</pub-id></element-citation></ref><ref id="bib82"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rohenkohl</surname><given-names>G</given-names></name><name><surname>Bosman</surname><given-names>CA</given-names></name><name><surname>Fries</surname><given-names>P</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Gamma Synchronization between V1 and V4 Improves Behavioral Performance</article-title><source>Neuron</source><volume>100</volume><fpage>953</fpage><lpage>963</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2018.09.019</pub-id><pub-id pub-id-type="pmid">30318415</pub-id></element-citation></ref><ref id="bib83"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ruiz</surname><given-names>O</given-names></name><name><surname>Lustig</surname><given-names>BR</given-names></name><name><surname>Nassi</surname><given-names>JJ</given-names></name><name><surname>Cetin</surname><given-names>A</given-names></name><name><surname>Reynolds</surname><given-names>JH</given-names></name><name><surname>Albright</surname><given-names>TD</given-names></name><name><surname>Callaway</surname><given-names>EM</given-names></name><name><surname>Stoner</surname><given-names>GR</given-names></name><name><surname>Roe</surname><given-names>AW</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Optogenetics through windows on the brain in the nonhuman primate</article-title><source>Journal of Neurophysiology</source><volume>110</volume><fpage>1455</fpage><lpage>1467</lpage><pub-id pub-id-type="doi">10.1152/jn.00153.2013</pub-id></element-citation></ref><ref id="bib84"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schafer</surname><given-names>RJ</given-names></name><name><surname>Moore</surname><given-names>T</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Selective attention from voluntary control of neurons in prefrontal cortex</article-title><source>Science</source><volume>332</volume><fpage>1568</fpage><lpage>1571</lpage><pub-id pub-id-type="doi">10.1126/science.1199892</pub-id><pub-id pub-id-type="pmid">21617042</pub-id></element-citation></ref><ref id="bib85"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schurger</surname><given-names>A</given-names></name><name><surname>Pereira</surname><given-names>F</given-names></name><name><surname>Treisman</surname><given-names>A</given-names></name><name><surname>Cohen</surname><given-names>JD</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>Reproducibility distinguishes conscious from nonconscious neural representations</article-title><source>Science</source><volume>327</volume><fpage>97</fpage><lpage>99</lpage><pub-id pub-id-type="doi">10.1126/science.1180029</pub-id><pub-id pub-id-type="pmid">19965385</pub-id></element-citation></ref><ref id="bib86"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Semedo</surname><given-names>JD</given-names></name><name><surname>Jasper</surname><given-names>AI</given-names></name><name><surname>Zandvakili</surname><given-names>A</given-names></name><name><surname>Krishna</surname><given-names>A</given-names></name><name><surname>Aschner</surname><given-names>A</given-names></name><name><surname>Machens</surname><given-names>CK</given-names></name><name><surname>Kohn</surname><given-names>A</given-names></name><name><surname>Yu</surname><given-names>BM</given-names></name></person-group><year iso-8601-date="2022">2022</year><article-title>Feedforward and feedback interactions between visual cortical areas use different population activity patterns</article-title><source>Nature Communications</source><volume>13</volume><elocation-id>1099</elocation-id><pub-id pub-id-type="doi">10.1038/s41467-022-28552-w</pub-id><pub-id pub-id-type="pmid">35232956</pub-id></element-citation></ref><ref id="bib87"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Siapas</surname><given-names>AG</given-names></name><name><surname>Lubenov</surname><given-names>EV</given-names></name><name><surname>Wilson</surname><given-names>MA</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>Prefrontal phase locking to hippocampal theta oscillations</article-title><source>Neuron</source><volume>46</volume><fpage>141</fpage><lpage>151</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2005.02.028</pub-id><pub-id pub-id-type="pmid">15820700</pub-id></element-citation></ref><ref id="bib88"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Spitzer</surname><given-names>H</given-names></name><name><surname>Desimone</surname><given-names>R</given-names></name><name><surname>Moran</surname><given-names>J</given-names></name></person-group><year iso-8601-date="1988">1988</year><article-title>Increased attention enhances both behavioral and neuronal performance</article-title><source>Science</source><volume>240</volume><fpage>338</fpage><lpage>340</lpage><pub-id pub-id-type="doi">10.1126/science.3353728</pub-id><pub-id pub-id-type="pmid">3353728</pub-id></element-citation></ref><ref id="bib89"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Steinmetz</surname><given-names>NA</given-names></name><name><surname>Moore</surname><given-names>T</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Eye movement preparation modulates neuronal responses in area V4 when dissociated from attentional demands</article-title><source>Neuron</source><volume>83</volume><fpage>496</fpage><lpage>506</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2014.06.014</pub-id><pub-id pub-id-type="pmid">25033188</pub-id></element-citation></ref><ref id="bib90"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Steriade</surname><given-names>M</given-names></name><name><surname>Nuñez</surname><given-names>A</given-names></name><name><surname>Amzica</surname><given-names>F</given-names></name></person-group><year iso-8601-date="1993">1993</year><article-title>A novel slow (&lt; 1 Hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components</article-title><source>The Journal of Neuroscience</source><volume>13</volume><fpage>3252</fpage><lpage>3265</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.13-08-03252.1993</pub-id><pub-id pub-id-type="pmid">8340806</pub-id></element-citation></ref><ref id="bib91"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tang</surname><given-names>L</given-names></name><name><surname>Higley</surname><given-names>MJ</given-names></name></person-group><year iso-8601-date="2020">2020</year><article-title>Layer 5 circuits in V1 differentially control visuomotor behavior</article-title><source>Neuron</source><volume>105</volume><fpage>346</fpage><lpage>354</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2019.10.014</pub-id><pub-id pub-id-type="pmid">31757603</pub-id></element-citation></ref><ref id="bib92"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>van Kempen</surname><given-names>J</given-names></name><name><surname>Gieselmann</surname><given-names>MA</given-names></name><name><surname>Boyd</surname><given-names>M</given-names></name><name><surname>Steinmetz</surname><given-names>NA</given-names></name><name><surname>Moore</surname><given-names>T</given-names></name><name><surname>Engel</surname><given-names>TA</given-names></name><name><surname>Thiele</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2021">2021</year><article-title>Top-down coordination of local cortical state during selective attention</article-title><source>Neuron</source><volume>109</volume><fpage>894</fpage><lpage>904</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2020.12.013</pub-id><pub-id pub-id-type="pmid">33406410</pub-id></element-citation></ref><ref id="bib93"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>van Vugt</surname><given-names>B</given-names></name><name><surname>Dagnino</surname><given-names>B</given-names></name><name><surname>Vartak</surname><given-names>D</given-names></name><name><surname>Safaai</surname><given-names>H</given-names></name><name><surname>Panzeri</surname><given-names>S</given-names></name><name><surname>Dehaene</surname><given-names>S</given-names></name><name><surname>Roelfsema</surname><given-names>PR</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>The threshold for conscious report: Signal loss and response bias in visual and frontal cortex</article-title><source>Science</source><volume>360</volume><fpage>537</fpage><lpage>542</lpage><pub-id pub-id-type="doi">10.1126/science.aar7186</pub-id></element-citation></ref><ref id="bib94"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vijayraghavan</surname><given-names>S</given-names></name><name><surname>Wang</surname><given-names>M</given-names></name><name><surname>Birnbaum</surname><given-names>SG</given-names></name><name><surname>Williams</surname><given-names>GV</given-names></name><name><surname>Arnsten</surname><given-names>AFT</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Inverted-U dopamine D1 receptor actions on prefrontal neurons engaged in working memory</article-title><source>Nature Neuroscience</source><volume>10</volume><fpage>376</fpage><lpage>384</lpage><pub-id pub-id-type="doi">10.1038/nn1846</pub-id><pub-id pub-id-type="pmid">17277774</pub-id></element-citation></ref><ref id="bib95"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vinck</surname><given-names>M</given-names></name><name><surname>van Wingerden</surname><given-names>M</given-names></name><name><surname>Womelsdorf</surname><given-names>T</given-names></name><name><surname>Fries</surname><given-names>P</given-names></name><name><surname>Pennartz</surname><given-names>CMA</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>The pairwise phase consistency: a bias-free measure of rhythmic neuronal synchronization</article-title><source>NeuroImage</source><volume>51</volume><fpage>112</fpage><lpage>122</lpage><pub-id pub-id-type="doi">10.1016/j.neuroimage.2010.01.073</pub-id><pub-id pub-id-type="pmid">20114076</pub-id></element-citation></ref><ref id="bib96"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vinck</surname><given-names>M</given-names></name><name><surname>Batista-Brito</surname><given-names>R</given-names></name><name><surname>Knoblich</surname><given-names>U</given-names></name><name><surname>Cardin</surname><given-names>JA</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Arousal and locomotion make distinct contributions to cortical activity patterns and visual encoding</article-title><source>Neuron</source><volume>86</volume><fpage>740</fpage><lpage>754</lpage><pub-id pub-id-type="doi">10.1016/j.neuron.2015.03.028</pub-id><pub-id pub-id-type="pmid">25892300</pub-id></element-citation></ref><ref id="bib97"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Watson</surname><given-names>AB</given-names></name></person-group><year iso-8601-date="1979">1979</year><article-title>Probability summation over time</article-title><source>Vision Research</source><volume>19</volume><fpage>515</fpage><lpage>522</lpage><pub-id pub-id-type="doi">10.1016/0042-6989(79)90136-6</pub-id><pub-id pub-id-type="pmid">483579</pub-id></element-citation></ref><ref id="bib98"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Westerberg</surname><given-names>JA</given-names></name><name><surname>Sigworth</surname><given-names>EA</given-names></name><name><surname>Schall</surname><given-names>JD</given-names></name><name><surname>Maier</surname><given-names>A</given-names></name></person-group><year iso-8601-date="2021">2021</year><article-title>Pop-out search instigates beta-gated feature selectivity enhancement across V4 layers</article-title><source>PNAS</source><volume>118</volume><elocation-id>e2103702118</elocation-id><pub-id pub-id-type="doi">10.1073/pnas.2103702118</pub-id><pub-id pub-id-type="pmid">34893538</pub-id></element-citation></ref><ref id="bib99"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Westerberg</surname><given-names>JA</given-names></name><name><surname>Schall</surname><given-names>MS</given-names></name><name><surname>Maier</surname><given-names>A</given-names></name><name><surname>Woodman</surname><given-names>GF</given-names></name><name><surname>Schall</surname><given-names>JD</given-names></name></person-group><year iso-8601-date="2022">2022</year><article-title>Laminar microcircuitry of visual cortex producing attention-associated electric fields</article-title><source>eLife</source><volume>11</volume><elocation-id>e72139</elocation-id><pub-id pub-id-type="doi">10.7554/eLife.72139</pub-id><pub-id pub-id-type="pmid">35089128</pub-id></element-citation></ref><ref id="bib100"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wichmann</surname><given-names>FA</given-names></name><name><surname>Hill</surname><given-names>NJ</given-names></name></person-group><year iso-8601-date="2001">2001</year><article-title>The psychometric function: I. Fitting, sampling, and goodness of fit</article-title><source>Perception &amp; Psychophysics</source><volume>63</volume><fpage>1293</fpage><lpage>1313</lpage><pub-id pub-id-type="doi">10.3758/BF03194544</pub-id></element-citation></ref><ref id="bib101"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wonders</surname><given-names>CP</given-names></name><name><surname>Anderson</surname><given-names>SA</given-names></name></person-group><year iso-8601-date="2006">2006</year><article-title>The origin and specification of cortical interneurons</article-title><source>Nature Reviews. Neuroscience</source><volume>7</volume><fpage>687</fpage><lpage>696</lpage><pub-id pub-id-type="doi">10.1038/nrn1954</pub-id><pub-id pub-id-type="pmid">16883309</pub-id></element-citation></ref><ref id="bib102"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Worden</surname><given-names>MS</given-names></name><name><surname>Foxe</surname><given-names>JJ</given-names></name><name><surname>Wang</surname><given-names>N</given-names></name><name><surname>Simpson</surname><given-names>GV</given-names></name></person-group><year iso-8601-date="2000">2000</year><article-title>Anticipatory biasing of visuospatial attention indexed by retinotopically specific alpha-band electroencephalography increases over occipital cortex</article-title><source>The Journal of Neuroscience</source><volume>20</volume><elocation-id>RC63</elocation-id><pub-id pub-id-type="doi">10.1523/JNEUROSCI.20-06-j0002.2000</pub-id><pub-id pub-id-type="pmid">10704517</pub-id></element-citation></ref><ref id="bib103"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yerkes</surname><given-names>RM</given-names></name><name><surname>Dodson</surname><given-names>JD</given-names></name></person-group><year iso-8601-date="1908">1908</year><article-title>The relation of strength of stimulus to rapidity of habit‐formation</article-title><source>Journal of Comparative Neurology and Psychology</source><volume>18</volume><fpage>459</fpage><lpage>482</lpage><pub-id pub-id-type="doi">10.1002/cne.920180503</pub-id></element-citation></ref><ref id="bib104"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname><given-names>G</given-names></name><name><surname>Herman</surname><given-names>JP</given-names></name><name><surname>Katz</surname><given-names>LN</given-names></name><name><surname>Krauzlis</surname><given-names>RJ</given-names></name></person-group><year iso-8601-date="2022">2022</year><article-title>Microsaccades as a marker not a cause for attention-related modulation</article-title><source>eLife</source><volume>11</volume><elocation-id>e74168</elocation-id><pub-id pub-id-type="doi">10.7554/eLife.74168</pub-id><pub-id pub-id-type="pmid">35289268</pub-id></element-citation></ref><ref id="bib105"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zeng</surname><given-names>H</given-names></name><name><surname>Sanes</surname><given-names>JR</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Neuronal cell-type classification: challenges, opportunities and the path forward</article-title><source>Nature Reviews. Neuroscience</source><volume>18</volume><fpage>530</fpage><lpage>546</lpage><pub-id pub-id-type="doi">10.1038/nrn.2017.85</pub-id><pub-id pub-id-type="pmid">28775344</pub-id></element-citation></ref><ref id="bib106"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zénon</surname><given-names>A</given-names></name><name><surname>Krauzlis</surname><given-names>RJ</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Attention deficits without cortical neuronal deficits</article-title><source>Nature</source><volume>489</volume><fpage>434</fpage><lpage>437</lpage><pub-id pub-id-type="doi">10.1038/nature11497</pub-id><pub-id pub-id-type="pmid">22972195</pub-id></element-citation></ref><ref id="bib107"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zohary</surname><given-names>E</given-names></name><name><surname>Shadlen</surname><given-names>MN</given-names></name><name><surname>Newsome</surname><given-names>WT</given-names></name></person-group><year iso-8601-date="1994">1994</year><article-title>Correlated neuronal discharge rate and its implications for psychophysical performance</article-title><source>Nature</source><volume>370</volume><fpage>140</fpage><lpage>143</lpage><pub-id pub-id-type="doi">10.1038/370140a0</pub-id></element-citation></ref><ref id="bib108"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zuber</surname><given-names>BL</given-names></name><name><surname>Stark</surname><given-names>L</given-names></name></person-group><year iso-8601-date="1966">1966</year><article-title>Saccadic suppression: Elevation of visual threshold associated with saccadic eye movements</article-title><source>Experimental Neurology</source><volume>16</volume><fpage>65</fpage><lpage>79</lpage><pub-id pub-id-type="doi">10.1016/0014-4886(66)90087-2</pub-id></element-citation></ref></ref-list></back><sub-article article-type="editor-report" id="sa0"><front-stub><article-id pub-id-type="doi">10.7554/eLife.91722.4.sa0</article-id><title-group><article-title>eLife Assessment</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Gold</surname><given-names>Joshua I</given-names></name><role specific-use="editor">Reviewing Editor</role><aff><institution>University of Pennsylvania</institution><country>United States</country></aff></contrib></contrib-group><kwd-group kwd-group-type="evidence-strength"><kwd>Solid</kwd></kwd-group><kwd-group kwd-group-type="claim-importance"><kwd>Useful</kwd></kwd-group></front-stub><body><p>This <bold>useful</bold> study by Nandy and colleagues examined relationships between behavioral state, neural activity in cortical area V4, and trial-by-trial variability in the ability to detect weak visual stimuli. They present <bold>solid</bold> evidence indicating that certain changes in arousal and eye-position stability, along with patterns of synchrony in the activity of neurons in different layers of V4, can show modest correspondences to changes in the ability to correctly detect a stimulus. These findings are likely to be of interest to those who seek a deeper understanding of circuit mechanisms that underlie perception.</p></body></sub-article><sub-article article-type="referee-report" id="sa1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.91722.4.sa1</article-id><title-group><article-title>Reviewer #1 (Public review):</article-title></title-group><contrib-group><contrib contrib-type="author"><anonymous/><role specific-use="referee">Reviewer</role></contrib></contrib-group></front-stub><body><p>Summary:</p><p>In this study, Nandy and colleagues examine neural, physiological and behavioral correlates of perceptual variability in monkeys performing a visual change detection task. They used a laminar probe to record from area V4 while two macaque monkeys detected a small change in stimulus orientation that occurred at a random time in one of two locations, focusing their analysis on stimulus conditions where the animal was equally likely to detect (hit) or not-detect (miss) a briefly presented orientation change (target). They discovered two behavioral and physiological measures that are significantly different between hit and miss trials - pupil size tends to be slightly larger on hits vs. misses, and monkeys are more likely to miss the target on trials in which they made a microsaccade shortly before target onset. They also examined multiple measures of neural activity across the cortical layers and found some measures that are significantly different between hits and misses.</p><p>Strengths:</p><p>Overall the study is well executed and the analyses are appropriate (with some possible caveats discussed below).</p><p>Weaknesses:</p><p>I have two remaining concerns. First, with the exception of the pre-target microsaccades, the correlates of perceptual variability (differences between hits and misses) appear to be weak and disconnected. The GLM analysis of the predictive power of trial outcome based on the behavioral and neural measures is only discussed at the end of the paper. This analysis shows that some of the measures have no significant predictive power, while others cannot be examined using the GLM analysis because these measures cannot be estimated in single trials. Given these weak and disconnected effects, my overall sense is that the current results provide a limited advance to our understanding of the neural basis of perceptual variability.</p><p>In addition, because the authors combine data across stimulus contrasts, I am somewhat uneasy about the possible confounding effect of contrast. As expected, stimulus contrast affected the probability of hits vs. misses. Independently, contrast may have affected some of the physiological measurements. Therefore, showing that contrast is not the source of the covariations between the physiological/behavioral measurements and perception can be challenging, and I am not convinced that the authors have ruled this out as a possible confound. It is unclear why the authors had to vary contrast in the first place, and why the analyses had to be done by combining the data across contrasts or by ignoring contrast as a variable (e.g., in the GLM analysis).</p></body></sub-article><sub-article article-type="author-comment" id="sa2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.91722.4.sa2</article-id><title-group><article-title>Author response</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Morton</surname><given-names>Mitchell P</given-names></name><role specific-use="author">Author</role><aff><institution>Yale University</institution><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib><contrib contrib-type="author"><name><surname>Denagamage</surname><given-names>Sachira</given-names></name><role specific-use="author">Author</role><aff><institution>Yale University</institution><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib><contrib contrib-type="author"><name><surname>Blume</surname><given-names>Isabel J</given-names></name><role specific-use="author">Author</role><aff><institution>Yale University</institution><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib><contrib contrib-type="author"><name><surname>Reynolds</surname><given-names>John H</given-names></name><role specific-use="author">Author</role><aff><institution>Salk Institute for Biological Studies</institution><addr-line><named-content content-type="city">La Jolla</named-content></addr-line><country>United States</country></aff></contrib><contrib contrib-type="author"><name><surname>Jadi</surname><given-names>Monika P</given-names></name><role specific-use="author">Author</role><aff><institution>Yale University</institution><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib><contrib contrib-type="author"><name><surname>Nandy</surname><given-names>Anirvan S</given-names></name><role specific-use="author">Author</role><aff><institution>Yale University</institution><addr-line><named-content content-type="city">New Haven</named-content></addr-line><country>United States</country></aff></contrib></contrib-group></front-stub><body><p>The following is the authors’ response to the previous reviews.</p><disp-quote content-type="editor-comment"><p>As you can see from the assessment (which is unchanged from before) and the reviews included below, the reviewers felt that the revisions did not yet address all of the major concerns. There was agreement that the strength of evidence would be upgraded to &quot;solid&quot; by addressing, at minimum, the following:</p><p>(1) Which of the results are significant for individual monkeys; and</p><p>(2) How trials from different target contrasts were analyzed</p></disp-quote><p>In this revision, we have addressed the two primary editorial recommendations:</p><p>(1) We apologize if this information was not clear in the previous version. We have updated Table 1 to highlight clearly the significant results for individual monkeys. Six of our key results – pupil diameter (Fig 2B), microsaccades (Fig 2D), decoding performance for narrow-spiking units (Fig 3A), decoding performance for broad-spiking units (Fig 3B), target-evoked firing rate for all units (Fig 3E) and target-evoked firing rate for broad-spiking units (Fig 3F) – are significant for individual animals and therefore gives us high confidence regarding our results. Please also note that we present all results for individual animals in the Supplementary figures accompanying each main figure.</p><p>(2) We have updated the manuscript and methods to explain how trials of each contrast were included in each analysis, and how contrast normalization was performed for the analysis in Figure 3. In addition, we discuss this point in the Discussion section, which we quote below:</p><p>“Non-target stimulus contrasts were slightly different between hits and misses (mean: 33.1% in hits, 34.0% in misses, permutation test, 𝑝 = 0.02), but the contrast of the target was higher in hits compared to misses (mean: 38.7% in hits, 27.7% in misses, permutation test, 𝑝 = 1.6 𝑒 − 31). To control for potential effects of stimulus contrast, firing rates were first normalized by contrast before performing the analyses reported in Figure 3. For all other results, we considered only non-target stimuli, which had very minor differences in contrast (&lt;1%) across hits and misses. In fact, this minor difference was in the opposite direction of our results with mean contrast being slightly higher for misses. While we cannot completely rule out any other effects of stimulus contrast, the normalization in Figure 3 and minor differences for non-target stimuli should minimize them.”</p><disp-quote content-type="editor-comment"><p><bold>Reviewer #1 (Public Review):</bold></p><p>Summary:</p><p>In this study, Nandy and colleagues examine neural, physiological and behavioral correlates of perceptual variability in monkeys performing a visual change detection task. They used a laminar probe to record from area V4 while two macaque monkeys detected a small change in stimulus orientation that occurred at a random time in one of two locations, focusing their analysis on stimulus conditions where the animal was equally likely to detect (hit) or not-detect (miss) a briefly presented orientation change (target). They discovered two behavioral and physiological measures that are significantly different between hit and miss trials - pupil size tends to be slightly larger on hits vs. misses, and monkeys are more likely to miss the target on trials in which they made a microsaccade shortly before target onset. They also examined multiple measures of neural activity across the cortical layers and found some measures that are significantly different between hits and misses.</p><p>Strengths:</p><p>Overall the study is well executed and the analyses are appropriate (though several issues still need to be addressed as discussed in Specific Comments).</p></disp-quote><p>Thank you.</p><disp-quote content-type="editor-comment"><p>Weaknesses:</p><p>My main concern with this study is that, with the exception of the pre-target microsaccades, the correlates of perceptual variability (differences between hits and misses) appear to be weak, potentially unreliable and disconnected. The GLM analysis of predictive power of trial outcome based on the behavioral and neural measures is only discussed at the end of the paper. This analysis shows that some of the measures have no significant predictive power, while others cannot be examined using the GLM analysis because these measures cannot be estimated in single trials. Given these weak and disconnected effects, my overall sense is that the current results provide limited advance to our understanding of the neural basis of perceptual variability.</p></disp-quote><p>Please see our response above to item #1 of the editorial recommendation. Six of our key results are individually significant in both animals giving us high confidence about the reliability and strength of our results.</p><p>Regarding the reviewer’s comment about the GLM, we note (also stated in the manuscript) that among the measures that we could estimate reliably on a single trial basis, two of these – pre-target microsaccades and input-layer firing rates – were reliable signatures of stimulus perception at threshold. This analysis does not imply that the other measures – Fano Factor, PPC, inter-laminar population correlations, SSC (which are all standard tools in modern systems neuroscience, and which cannot be estimated on a single-trial basis) – are irrelevant. Our intent in including the GLM analyses was to <italic>complement</italic> the results reported from these across-trial measures (Figs 4-7) with the predictive power of single-trial measures.</p><p>While no study is entirely complete in itself, we have attempted to synthesize our results into a conceptual model as depicted in Fig 8.</p><disp-quote content-type="editor-comment"><p><bold>Reviewer #2 (Public Review):</bold></p><p>Strengths:</p><p>The experiments were well-designed and executed with meticulous control. The analyses of both behavioural and electrophysiological data align with the standards in the field.</p></disp-quote><p>Thank you.</p><disp-quote content-type="editor-comment"><p>Weaknesses:</p><p>Many of the findings appear to be subtle differences and incremental compared to previous literature, including the authors' own work. While incremental findings are not necessarily a problem, the manuscript lacks clear statements about the extent to which the dataset, analysis, and findings overlap with the authors' prior research. For example, one of the main findings, which suggests that V4 neurons exhibit larger visual responses in hit trials (as shown in Fig. 3), appears to have been previously reported in their 2017 paper.</p></disp-quote><p>We respectfully disagree with the assessment that the findings reported here are incremental over the results reported in our prior study (Nandy et al,. 2017). In the previous study, we compared the laminar profile of neural modulation due to the deployment of attention i.e. the main comparison points were the attend-in and the attend-away conditions while controlling for visual stimulation. In this study, we go one step further and home in on the attend-in condition and investigate the differences in the laminar profile of neural activity (and two additional physiological measures: pupil and microsaccades) when the animal either correctly reports or fails to report a stimulus with equal probability. We thus control for <italic>both</italic> the visual stimulation and the cued attention state of the animal. While there are parallels to our previous results (as the reviewer correctly noted), the results reported here cannot be trivially predicted from our previous results. Please also note that we discuss our new results in the context of prior results, from both our group and others, in the manuscript (lines 310-332).</p><disp-quote content-type="editor-comment"><p>Furthermore, the manuscript does not explore potentially interesting aspects of the dataset. For instance, the authors could have investigated instances where monkeys made 'false' reports, such as executing saccades towards visual stimuli when no orientation change occurred, which allows for a broader analysis that considers the perceptual component of neural activity over pure sensory responses. Overall, lacking broad interest with the current form.</p></disp-quote><p>We appreciate the reviewer’s feedback on analyzing false alarm trials. Our focus for this study was to investigate the behavioral and neural correlates accompanying a correct or incorrect perception of a target stimulus presented at perceptual threshold. False alarm trials, by definition, do not include a target presentation. Moreover, false alarm rates rapidly decline with duration into a trial, with high rates during the first non-target presentation and rates close to zero by the time of the eighth presentation (see figure). Investigating false alarms will thus involve a completely different form of analysis than we have undertaken here. We therefore feel that while analyzing false alarm trials will be an interesting avenue to pursue in the future, it is outside the scope of the present study.</p><fig id="sa2fig1" position="float"><label>Author response image 1.</label><graphic mimetype="image" mime-subtype="tiff" xlink:href="elife-91722-sa2-fig1-v1.tif"/></fig></body></sub-article></article>