Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI

Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI View ORCID Profile Marianna E. Schmidt , View ORCID Profile Iman Aganj , View ORCID Profile Jason Stockmann , View ORCID Profile Berkin Bilgic , Yulin Chang , W. Scott Hoge , Evgeniya Kirilina , View ORCID Profile Nikolaus Weiskopf , View ORCID Profile Shahin Nasr doi: https://doi.org/10.1101/2025.03.27.645795 Marianna E. Schmidt 1 Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences , Leipzig, Germany 2 Max Planck School of Cognition , Leipzig, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Marianna E. Schmidt Iman Aganj 3 Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital , Boston, MA, USA 4 Radiology Department, Harvard Medical School , Boston, MA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Iman Aganj Jason Stockmann 3 Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital , Boston, MA, USA 4 Radiology Department, Harvard Medical School , Boston, MA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jason Stockmann Berkin Bilgic 3 Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital , Boston, MA, USA 4 Radiology Department, Harvard Medical School , Boston, MA, USA 5 Harvard/MIT Health Sciences and Technology , Cambridge, MA, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Berkin Bilgic Yulin Chang 6 Siemens Medical Solutions USA Inc. , Malvern, PA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site W. Scott Hoge 7 Imaginostics, Inc. , Orlando, FL, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Evgeniya Kirilina 1 Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences , Leipzig, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Nikolaus Weiskopf 1 Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences , Leipzig, Germany 8 Felix Bloch Institute for Solid State Physics, Faculty of Physics and Earth System Sciences, Leipzig University , Leipzig, Germany 9 Wellcome Center for Human Neuroimaging, Institute of Neurology, University College London , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Nikolaus Weiskopf Shahin Nasr 3 Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital , Boston, MA, USA 4 Radiology Department, Harvard Medical School , Boston, MA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Shahin Nasr For correspondence: shahin.nasr{at}mgh.harvard.edu Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Despite their importance for visual perception, functional connectivity between ocular dominance columns (ODCs)in the human primary visual cortex (V1) remains largely unknown. Using high- resolution functional MRI (fMRI), we localized ODCs and assessed their resting-state functional connectivity (rs-FC) in 12 human adults. Consistent with anatomical studies in animals, we found stronger rs-FC in the middle compared to deep and superficial cortical depths and selectively stronger rs- FC between ODCs with alike compared to unalike ocular polarity. Beyond what was known from animal models, and consistent with human perceptual biases, intra- and interhemispheric rs-FC was stronger in peripheral (compared to central) and in dorsal (compared to ventral) V1 subregions. Lastly, rs-FC patterns correlated with ODC maps, suggesting that ODC maps can be predicted from rs-FC patterns within V1. These results highlight the heterogeneity in rs-FC between ODCs across cortical depths and V1 subfields, underscoring their likely association with human perceptual features. 1. Introduction In humans and many animals, the visual cortex is organized into neuronal columns that encode visual features. Despite their central role in shaping perception, the functional connectivity between cortical columns in the human visual cortex remains poorly understood. This gap in our knowledge stems primarily from the limited spatial resolution of conventional neuroimaging techniques relative to the size of cortical columns and layers. As a result, much of our understanding of the mesoscale organization of the visual cortex relies on anatomical studies in animals using invasive techniques that are not applicable to humans in vivo. To bridge this knowledge gap, we leveraged high-resolution functional MRI (fMRI) to investigate resting-state functional connectivity (rs-FC) between ocular dominance columns (ODCs), which are considered the building blocks of visual processing in the primary visual cortex (V1). Anatomical studies in non-human primates (NHPs) have revealed that ODCs as far as 8 mm apart are systematically interconnected via horizontal fibers 1 – 6 . These fibers are most dense in the middle cortical layers, particularly layer 4B 6 , preferentially connecting ODCs with alike rather than unalike ocular polarity 4 , 7 , 8 . The level of this bias in the number of horizontal fibers remains a matter of debate varying from 4% 4 to 15% 7 in favor of the ODC pairs with alike ocular polarity. However, the reason for this inconsistency could not be tested systematically due to limited sample sizes. Functional evidence for the biased structural connectivity between alike ODCs is provided by electrophysiological 9 – 11 and optical imaging studies 12 . According to these studies, neurons with alike rather than unalike ocular polarity exhibit stronger synchronized activity. Interestingly, this synchronous activity is detectable over distances up to 20 mm 12 indicating that the selective functional connectivity between ODCs extends far beyond the range of monosynaptic horizontal fibers and cannot be attributed to the overlapping receptive fields of the targeted neurons 11 . To our knowledge, evidence for homologous functional organization within human V1 has not yet been reported. This is mainly due to difficulties in studying ODCs based on conventional neuroimaging techniques that provide a low signal-to-noise ratio and poor spatial resolution (>2mm), inadequate for studying the mesoscale functional organization of the human visual cortex. Nonetheless, at coarser levels, fMRI studies have shown higher intra- and interhemispheric rs-FC between V1 subregions representing similar locations of the visual field 13 – 15 . Outside V1, high-resolution fMRI has successfully been used to study rs-FC between color-, stereo- as well as motion-selective cortical columns across extrastriate visual areas 16 – 19 . Findings from these studies have supported the hypothesis that cortical columns with alike stimulus preferences are more selectively connected to each other, and that this selectivity can be detected by high-resolution resting-state fMRI (see also 20 ). Using high-resolution fMRI in the past, multiple groups have successfully localized ODCs in humans 16 , 21 – 27 . In this study, we extended the application of high-resolution fMRI to study the rs-FC between ODCs, in the absence of any visual input (eyes closed), in 12 adult humans. We examined the influence of cortical depth, pairwise distance, relative ocular polarity (i.e., alike vs. unalike ocular preference), and strength of ocular preference on the mesoscale rs-FC within V1 and assessed the heterogeneity of the measurements across V1 subregions. These analyses were then extended to evaluate the mesoscale rs-FC between V1 subregions within and between the two hemispheres. Lastly, we tested the predictability of ODC maps based on the mesoscale rs-FC pattern within V1. In agreement with animal studies, we found selective mesoscale rs-FC between ODCs with stronger connectivity at middle compared to deep and superficial cortical depth levels. Beyond what was previously known from animal studies, we showed that these connectivity patterns are not uniform but vary across V1 subregions within and between hemispheres, consistent with the heterogeneity of visual processing across subfields. Finally, we demonstrated that rs-FC patterns correlate with ODC maps, providing an initial step toward predicting and segmenting the mesoscale functional organization of the human visual cortex. A preliminary version of the results has been already presented at a workshop 28 . 2. Results Studies in NHPs suggest that horizontal fibers in primate V1 are more numerous in the middle layer (layer 4B) compared to deeper (layers 5 and 6) and superficial cortical layers (layers 2 and 3) 6 . They also indicate a significant difference (4% to 15%) in the number of fibers connecting ODCs, favoring pairs with alike rather than unalike ocular polarity 4 , 7 , 8 . Here, we first examined whether these mesoscale features are also reflected in rs-FC between vertex pairs in human V1. Then, we extended our analysis to investigate the influence of ocular preference strength, the heterogeneity of these features across different V1 subregions and their presence in interhemispheric rs-FC measurements. To achieve these aims, we acquired high-resolution fMRI data at 1 mm isotropic resolution in the visual cortex of 12 adults ( Figure 1 , see Methods). ODCs were localized using dichoptic stimulation and resting- state fMRI data were collected in a separate session. Individual visual area borders and central, peripheral, dorsal, and ventral V1 subregions were defined using retinotopic mapping (eccentricity <10°), and rs-FC was quantified by correlating resting-state time series between cortical vertices within deep, middle and superficial cortical depth levels of V1. Connectivity was then assessed as a function of pairwise distance, ocular polarity (alike vs. unalike), cortical depth, pairwise ocular preference strength, visual field representation (central vs. peripheral, dorsal vs. ventral), and interhemispheric connections. Finally, rs-FC patterns were correlated with ODC maps to test their correspondence. Download figure Open in new tab Figure 1. Structural and functional MRI from one individual participant. A) Structural scan and reconstruction of cortical surfaces at different depth levels. For each participant, primary visual cortex (V1) was localized retinotopically along the calcarine sulcus (CS). Black lines indicate the borders of the stimulated portion of V1 (radius<10°). As a part of brain reconstruction, the cortical gray matter was divided into deep (blue), middle (orange), and superficial (green) depths with small gaps between them (see Methods). B) Differential ocular dominance column (ODC) map (beta values) overlaid on the flattened cortex. The ocular dominance response was derived by contrasting the response to stimulation of the dominant eye (DE; red to yellow) versus the non-dominant eye (NDE; blue to cyan). The map shows the stimulated portion of V1 (radius<10°) in the right hemisphere (RH) of a representative subject. Notably, a large activity patch along the dorsal part of the vertical meridian (white arrows), in which the ODCs could not be mapped due to technical limitations, was excluded from the analysis 23 . C) Mean spontaneous fMRI time course measured over the whole stimulated portion of V1. Consistent with previous studies, the signal intensity was weaker in deeper compared to more superficial depths. To exclude the potential signal instabilities at the beginning of the scan, the first 36 s were excluded from further analysis. D) Spontaneous fMRI time courses from three vertices sampled at deep, middle, and superficial layers, respectively, at a randomly selected single tangential location along the surface normal. Compared to the averaged signal, these individual laminar time courses exhibit lower temporal correlation, underscoring the likely presence of depth-specific fluctuations. 2.1. Rs-FC between alike and unalike ODCs Figure 2A illustrates the rs-FC levels, measured during the eyes-closed condition (see Methods), for vertex pairs with alike and unalike ocular polarity as a function of distance between the two vertices of the pair. A two-way RM-ANOVA (distance and relative ocular polarity), applied to the measured rs-FC levels, revealed that rs-FC was significantly stronger between alike than unalike vertex pairs (F(1, 11)=6.12, p =0.03). Here while rs-FC decreased significantly with increasing distance between vertex pair (F(9, 99)=29.94, p <10 - 3 ), the interaction between the effects of distance and relative ocular polarity remained non-significant (F(9, 99)=0.58, p =0.58). These results suggest that selective rs-FC between vertices with alike ocular polarity persists across relatively long distances (Supplementary Figure 1) Download figure Open in new tab Figure 2. The resting-state functional connectivity (rs-FC) between vertex pairs in the primary visual cortex and its selectivity in favor of the pair with alike ocular polarity as a function of distance. A) The distribution of rs-FC ( z -values) across subjects is displayed using split violin plots for vertex pairs with alike (orange) and unalike (violet) ocular polarity for 10 distance quantiles (with 1 being the closest but excluding vertices closer together than 3mm due to fMRI point spread function (see Methods) and 10 being the farthest distance (Supplementary Figure 1)). Rs-FC decreased with increasing distance. B) Selectivity of rs-FC, measured as the alike rs-FC minus the unalike rs-FC, is depicted similarly as violin plot for each distance quantile in percentage. To facilitate comparison with the previous findings in animal studies, values were normalized for overall rs-FC level. The red dashed line indicates 0. In both panels, the violins display the full distribution of the data, with thickness representing density and endpoints corresponding to the min and max outliers. The absence of interaction between the effects of distance and relative ocular polarity is also apparent in Figure 2B that illustrates the level of rs-FC selectivity, measured as rs-FC between alike – unalike vertex pairs, in each distance quantile. To make the rs-FC selectivity measurement similar to those in the previous anatomical studies 4 , 7 , 8 , we normalized it relative to the overall level of rs-FC in the same quantile. As demonstrated here, rs-FC selectivity remained on average positive even at the highest tested distance quantile, where vertex pairs were approximately 35 mm apart (Supplementary Figure 1), far exceeding the expected reach of a single horizontal fiber 1 , 4 , 5 . Notably, consistent with our previous reports in a smaller group of participants 23 , there was no significant difference between the number of vertices showing a preference for stimulation of the DE versus NDE ( p =0.10) or their ocular preference level as measure by an ocular dominance index (ODI; see Methods) ( p =0.53) (Supplementary Figure 2). Furthermore, a subsampling method assured us that the distance distribution was the same for alike and unalike vertex pairs (see Section 5.5.1 .). Thus, these findings cannot be attributed to an unbalanced distribution of ocular dominance activity in V1. 2.2. Reproducibility of rs-FC measurements across scan sessions We tested whether the selective rs-FC between alike ODCs could be reliably detected across different days. Four individuals were re-scanned on a separate day using a different ultra-high-field MRI scanner and a modified data acquisition protocol (see Section 5.3.2 .). Here again, a three-way RM-ANOVA (distance, relative ocular polarity, and scan session) applied to the measured rs-FC across different distance quantiles (Supplementary Figure 3) confirmed robust effects of relative ocular polarity (F(1, 3)=25.46, p =0.02), and distance (F(9, 27)=19.62, p =0.02), while revealing no significant main effect of scan session (F(1, 3)=0.69, p =0.47) and no interactions between the effects of scan session and the other two factors ( p >0.18). These results indicate that selective rs-FC between ODCs is reproducible over time and insensitive to scanner hardware and the exact acquisition parameters. Notably, the robustness of the ODC maps for our subjects has been tested and scrutinized in our previous studies 16 , 23 . 2.3. Selectivity of rs-FC as a function of ODI and cortical depth Studies in NHPs have suggested that the strength of ocular preference as measured by ODI may influence the selectivity of horizontal fibers 4 , 7 , 8 . Accordingly, here, we tested the hypothesis that rs-FC depends not only on the relative ocular polarity of the vertex pair but also on their ODI. Furthermore, we tested whether the level of rs-FC varied across cortical depth levels, as suggested by animal models 6 , and if rs- FC selectivity was detectable across all cortical depths. Considering the absence of significant interaction between the effects of relative ocular polarity and the distance between vertex pairs, data from all distance levels were averaged in all further analyses. Figure 3A illustrates the rs-FC between vertex pairs as a function of their ODI, relative ocular polarity, and cortical depth. A three-way RM-ANOVA revealed that rs-FC was significantly stronger between vertices with stronger ODIs than those with weaker ODIs (F(4, 44)=11.40, p <0.01) (Supplementary Table 1). As expected from our previous analysis, rs-FC was significantly stronger between vertex pairs with alike rather than unalike ocular polarity (F(1,11)=6.09, p =0.03). Here, a significant interaction between the effects of ODI and relative ocular polarity (F(4, 44)=6.08, p =0.02) demonstrated that rs-FC selectivity increases with the ODI of vertex pairs, reinforcing the link between functional connectivity and ocular dominance. This interaction is also apparent in Figure 3B that illustrates the level of rs-FC selectivity, measured as rs-FC between alike – unalike vertex pairs, without any normalization. Download figure Open in new tab Figure 3. The measured resting-state functional connectivity (rs-FC) between primary visual cortex vertex pairs and its selectivity as functions of cortical depth and ocular dominance index (ODI) averaged over all participants. A) The level of rs-FC between vertex pairs with alike (left) and unalike (right) ocular polarity measured across superficial (top), middle (middle), and deep (bottom) cortical depths. Each correlation matrix shows the mean rs-FC measured as z -values. In each matrix, x - and y -axes represent the ODIs of vertex pairs categorized in quantiles 1 to 10 representing low to high ODI beta values (Supplementary Figure 6). The measured rs-FC is generally higher at middle cortical depth and increases with ODI. B) Selectivity is defined as the difference between rs-FC of vertex pairs (alike minus unalike ocular polarity) without any normalization. The selectivity was generally higher, favoring pairs with alike ocular polarity and increased for vertex pairs with higher ODI. Moreover, consistent with the laminar organization of horizontal fibers in animals, rs-FC was highest in middle cortical depth compared to deeper and more superficial depths (F(2,22)=4.74, p =0.047). This effect was particularly pronounced in regions with stronger ODIs, resulting in a significant interaction between the effects of ODI and cortical depth (F(8, 88)=13.47, p <10 - 4 ). The lack of significant interaction between the effects of relative ocular polarity and cortical depth ( p =0.30) suggests that rs-FC selectivity remains equivalent across cortical depths, thus ruling out the possibility that rs-FC selectivity was arbitrarily stronger simply because the ODC map was always sampled from deep cortical depth (to minimize the blurring effect of pial veins). Similar results were also found after applying an LME model to the rs-FC data (see Supplementary Analysis, and Supplementary Table 2). Together, these results suggest that the influence of relative ocular polarity and cortical depth is contingent on higher ODIs. Considering the stronger fMRI signal in more superficial compared to deeper cortical depths across V1 (F(2,22)=74.23, p 49.98; p <10 - 4 ), these results also suggest that fMRI signal intensity does not dominate the strength of rs-FC measured at different cortical depths. Instead, other factors, such as the organization of horizontal fibers across cortical depths, may also influence the strength of rs-FC (see Discussion). 2.4. Heterogeneity in rs-FC Multiple previous studies suggested stronger global visual processing and greater sensitivity to low spatial frequency components —essential for global visual processing—in the lower compared to the upper visual field 29 – 32 . Given that sensitivity to lower spatial frequency components requires sampling from a larger visual field, it could be argued that rs-FC may be stronger in dorsal V1 subregions representing the lower visual field compared to ventral V1 subregions representing the upper visual field. To test this hypothesis, we divided V1 into two equivalent ROIs, representing the upper and lower portions of the visual field ( Figure 4A ), and compared the rs-FC pattern across them. Download figure Open in new tab Figure 4. Comparison of resting-state functional connectivity (rs-FC) and selectivity between ventral (V) and dorsal (D) primary visual cortex (V1) averaged across all participants. A) Location of dorsal and ventral V1 subregions, overlaid on the ocular dominance column (ODC) map in the left and right hemispheres (LH and RH, respectively) from one participant (other than the one shown in Figure 1 ). The stimulated part of V1 (dashed and dotted black lines) was divided into two equally sized areas along the calcarine sulcus (dashed white line). B) The measured rs-FC between (all) vertex pairs in dorsal (left) and ventral (right) V1 subregions irrespective of relative ocular polarity. The rs-FC level was on average higher in dorsal compared to ventral V1 subregion. C) The selectivity of rs-FC across the two subregions. In dorsal V1 subregion, the selectivity was more confined to sites with high ocular dominance index (ODI) whereas in the ventral V1 subregion, the selectivity was more uniformly distributed across ODI. The other details are the same as in Figure 3 . Figure 4B illustrates the level of rs-FC measured in dorsal V1 (representing the lower visual field) and ventral V1 (representing the upper visual field). Overall, consistent with our hypothesis, rs-FC was stronger in the dorsal compared to ventral V1. Additionally, rs-FC selectivity was more uniform — i.e., independent of ODIs — in ventral compared to dorsal subregions ( Figure 4C ). A four-way RM-ANOVA (ROI, relative ocular polarity, ODI, and cortical depth) (Supplementary Table 3) applied to the measured rs-FC revealed a significant main effect of ROI (F(1, 11)=29.26, p <10 - 3 ), as well as significant interactions between the effects of ROI, relative ocular polarity, and cortical depth (F(2, 22)=5.17, p =0.03) as well as between the effects of ROI, ODI, and cortical depth (F(8, 88)=4.38, p =0.02). These findings indicate significant heterogeneity in both the level of rs-FC and its selectivity across dorsal and ventral V1 subregions, as expected from the differential visual processing mechanism in these regions. A stronger preference for lower spatial frequency can also be found in V1 subregions representing peripheral compared to central visual fields 33 – 35 . Accordingly, we further tested if a similar heterogeneity can also be found between central and peripheral V1 subregions ( Figure 5A ). As demonstrated in Figure 5B , rs-FC was stronger in the peripheral compared to central V1 subregions across the superficial and middle (but not deep) cortical depths. The same pattern, but less pronounced, was found for the selectivity of rs-FC ( Figure 5C ). Here, the application of four-way RM-ANOVA (Supplementary Table 4) to the measured rs-FC showed a significant main effect of relative ocular polarity (F(1,11)=5.52; p =0.04) and a significant interaction between the effects of ROI and cortical depth (F(4, 44)=27.72, p =10 - 4 ). However, the two- and three-way interactions between the effects of ROI and relative ocular polarity remained non- significant (F0.53). These results suggested that the overall pattern of rs-FC, its dependency on ODI, relative ocular polarity, and cortical depth may vary between V1 subregions, and that this heterogeneity appears to be consistent with the known biases in human visual perception across visual fields. Download figure Open in new tab Figure 5. Comparison of resting-state functional connectivity (rs-FC) and selectivity between central (C) and peripheral (P) primary visual cortex (V1) averaged across all participants. A) Central (radius<4°) and peripheral V1 subregions are shown relative to the ocular dominance column (ODC) map in one participant (other than the one shown in Figures 1 and 4 ). B) The level of rs-FC was on average higher in peripheral V1 compared to central V1, which was more apparent in superficial and middle depths. C) Likewise, the selectivity of rs-FC was higher in the peripheral compared to the central V1 subregion in superficial and middle cortical depths. The other details are similar to Figure 4 . 2.5. Interhemispheric rs-FC between ODCs While monosynaptic callosal fibers directly linking V1 between hemispheres are limited 36 – 38 , numerous polysynaptic connections — primarily through feedback pathways from other callosally connected visual areas (e.g., V2) — may contribute to interhemispheric rs-FC in V1. Consistently, multiple studies have already shown interhemispheric rs-FC between V1 subregions that represent the same portion of visual field 13 – 15 . Here, we extended those findings by examining the dependence of interhemispheric rs-FC connecting V1 subregions (dorsal vs. ventral and central vs. peripheral) between the two hemispheres, on relative ocular polarity, ODI, and cortical depth. Because of the known retinotopic organization of interhemispheric connectivity, we limited our test to retinotopically corresponding areas, e.g., dorsal to dorsal and ventral to ventral V1. As illustrated in Figure 6A , interhemispheric rs-FC exhibited at least two regional differences between dorsal and ventral V1. First, similar to intrahemispheric rs-FC, dorsal V1 showed stronger interhemispheric rs-FC than ventral V1. Second, rs-FC selectivity ( Figure 6B ) followed distinct patterns in these subregions: in dorsal V1, vertex pairs with alike ocular polarity exhibited stronger rs-FC than those with unalike ocular polarity, whereas in ventral V1, this trend was reversed—rs-FC was stronger between vertex pairs with unalike ocular polarity than between those with alike ocular polarity. Download figure Open in new tab Figure 6. Interhemispheric resting-state functional connectivity (rs-FC) and its selectivity for vertex pairs located within dorsal and ventral primary visual cortex (V1) subregions averaged across all participants. A) Similar to the intrahemispheric rs-FC, the interhemispheric rs-FC was higher in the dorsal compared to ventral V1 subregion. The effects of cortical depth and ocular dominance index (ODI) that were observed for intrahemispheric rs-FC are preserved for interhemispheric rs-FC. B) In the dorsal V1 subregion, the selectivity of rs-FC was stronger and in favor of pairs with alike ocular polarity. Whereas in the ventral V1 subregions, the selectivity appeared to be weaker and reversed (i.e. in favor of vertex pairs with unalike ocular polarity). The other details are similar to Figure 4 . Consistently, a four-way RM-ANOVA (Supplementary Table 5) revealed a significant main effect of ROI (F(1, 11)=31.89, p <10 - 3 ), a significant interaction between ROI and relative ocular polarity (F(1,11)=6.67, p =0.03) and significant three-way interactions between the effects of ROI, relative ocular polarity, and ODI (F(4, 44)=5.35, p =0.02) as well as between ROI, ODI and cortical depth (F(8, 88)=3.16, p =0.04). Heterogeneous rs-FC was also observed when comparing interhemispheric rs-FC between central and peripheral V1 subregions. Specifically, again similar to intrahemispheric rs-FC, the relationship between rs-FC and cortical depth varied between central and peripheral regions ( Figure 7A ). Moreover, the overall pattern of rs-FC selectivity also differed— in the peripheral subregion, selectivity was more pronounced in superficial depths, favoring vertex pairs with alike ocular polarity ( Figure 7B ). In contrast, in the central subregion, we observed a weaker but opposite trend, with a slight preference for ODCs with unalike ocular polarity. A four-way RM-ANOVA on the measured rs-FC (Supplementary Table 6) revealed a significant interaction between the effect of ROI and cortical depth (F(2, 22)=24.77, p <10 - 3 ), as well as a significant three-way interaction between the effect of ROI, relative ocular polarity, and cortical depth (F(2, 22)=4.74, p =0.02) and four-way interaction between all factors (F(8,88)=3.03; p =0.048). Download figure Open in new tab Figure 7. Interhemispheric resting-state functional connectivity (rs-FC) and its selectivity for vertex pairs located within central and peripheral primary visual cortex (V1) subregions averaged across all participants. A) The level of interhemispheric rs-FC was, on average, higher in the peripheral compared to the central V1 subregion. Here again, this phenomenon resembles what we detected for the intrahemispheric rs-FC ( Figure 5 ). B) The selectivity of rs-FC was higher in peripheral compared to central V1 subregions mainly in superficial and middle cortical depths. The other details are similar to Figure 4 . These findings demonstrate that, like intrahemispheric rs-FC, interhemispheric V1 rs-FC is not uniform but instead varies systematically across V1 subregions, aligning with their specialized roles in visual processing. 2.6. Relating ODC maps to rs-FC Our findings suggest that for each vertex, the relative ocular polarity and the ODI of surrounding vertices may be related to (and even predictable based on) corresponding rs-FC measures. To rigorously test this hypothesis (H1), we pseudo-randomly selected 1,000 seed vertices (10.91% ± 1.76% of all vertices within our V1 ROI) in each hemisphere (see Section 5.5 .), generated their 2D rs-FC map with the other vertices, and quantified the correlation between the resultant rs-FC map and the ODC map centered around the seed point ( Figure 8A ). To mitigate the influence of the fMRI signal’s point spread function (PSF), we excluded the inner 3 mm radius from our analysis. We varied the outer radius from 4 to 10 mm in 0.5 mm increments and compared the measured correlation values against a stringent baseline, defined as the correlation between the same rs-FC pattern and a randomly selected, non-overlapping portion of the ODC map (H0 a ). These analyses were conducted separately for rs-FC measured within deep, middle, and superficial cortical depths. Download figure Open in new tab Figure 8. Correlation of resting-state functional connectivity (rs-FC) with the differential ocular dominance column (ODC) map. A) Schematic figure showing a map of rs-FC obtained by correlating the time-series of one seed vertex (x) with all other vertices (left) and the corresponding differential ODC map of that representative subject and hemisphere (right). Ring-shaped regions of the rs-FC map with varying outer radii and fixed inner radius of 3mm (white shaded inner circle, to exclude vertex pairs within the fMRI point spread function) were correlated with the ODC differential map at either the ring at the same location (H1) or a shifted location (H0 a ) or the same location but rotated by 180° (H0 b ). B) Resulting r -values ( y -axis) averaged across participants are plotted as a function of the ring radius ( x -axis) for the three different comparisons using rs-FC sampled at the three different cortical depths. The shaded areas correspond to the respective standard error of the mean. For all cortical depths and radii, r -values of the H1 were significantly higher than for the two null hypotheses (H0) indicating that the correlation is not only influenced by the overall topography of ODC patterns. As shown in Figure 8B , although the correlation between ODC and rs-FC maps decreased as the outer radius expanded, it remained higher than either baseline measure across all cortical depths. A three-way RM-ANOVA (measurement type (H1 vs. H0 a ), outer radius, cortical depth), applied to the measured correlation values yielded a significant effect of outer radius (F(12,132)=666.77, p <10 - 13 ) and measurement type (F(1, 11)=15.74, p <0.01), with no interaction between the two effects ( p= 0.13). There was no significant main effect of cortical depth ( p =0.31) or significant interaction between the effect of cortical depth and the other two independent factors ( p >0.29), indicating that the predictive relationship between rs-FC and ODC organization holds across different cortical depths. It has already been shown that the organization of ODCs may vary based on eccentricity in the visual field 39 . Given this possibility, it could be argued that the change in the location of the baseline disc could affect the measured correlation values, depending on the distance relative to the actual disc. To rule out this possibility, we conducted a control analysis using an alternative baseline in which the ODC map was only rotated 180° before computing its correlation with the rs-FC pattern (H0 b ). In this condition, the ODC and rs-FC rings remain overlapping relative to each other. Even under this control condition, an RM-ANOVA confirmed robust effects of outer radius (F(12,132)=856.01, p <10 - 14 ), measurement type (F(1, 11)=92.36, p <10 - 5 ), and the significant interaction between the effects of outer radius and measurement type (F(12, 132)=11.60), p 0.12), reinforcing the robustness of our findings across cortical depths. These results provide empirical support for the hypothesis that rs-FC patterns contain specific information about ODC organization, despite being derived from independent datasets collected in separate scan sessions. The persistence of this relationship across cortical depths and under rigorous control conditions underscores the fundamental link between rs-FC and the underlying columnar architecture of V1. 3. Discussion Rs-FC has been extensively used to investigate functional connections between predefined macroscale brain regions, primarily at the spatial scale of brain areas 40 , 41 , or to segment the brain into multiple co- fluctuating regions 42 – 44 . In this study, we provide compelling functional evidence for the selective rs-FC between ODCs — the building blocks of neuronal processing in V1. Our findings reveal key aspects of the mesoscale functional organization of V1, emphasizing heterogeneities across cortical depth and subregions, and offer a fresh perspective on its interhemispheric functional connectivity. Furthermore, we demonstrate new evidence for the predictability of ODC maps based on their rs-FC. 3.1. Consistency between rs-FC pattern and anatomical connections Our findings provide some of the first functional evidence of significantly stronger rs-FC between ODCs with alike compared to those with unalike ocular polarity, in humans, aligning with prior animal studies using postmortem histological methods 4 , 7 , 8 as well as optical imaging and electrophysiological recordings 9 – 12 . Furthermore, the observation that rs-FC is strongest in middle cortical layers—compared to deep and superficial layers—closely aligns with neuroanatomical studies that demonstrated the highest prevalence of horizontal fibers in layer 4B 6 . Despite the methodological differences between rs-FC and histological techniques based on using tracers, the overall level of rs-FC, and the level of rs-FC selectivity reported here ( Figure 2B ) is also comparable to the quantitative measure of this selectivity based on anatomical techniques. Specifically, three previous studies have reported 4% 4 , 8% 8 , and 15% 7 difference in the number of horizontal fibers connecting ODCs with alike rather than unalike relative ocular polarity. Given that at least two of these studies used the same tracer (biocytin) to detect the horizontal fibers 4 , 7 , the difference in the level of selectivity can at least be partly attributed to between-species (macaque nemestrina vs. macaque fascicularis) and between-individual variability. Beyond these functional parallels to structural measurements, our results highlight a strong influence of ODI on rs-FC strength. This effect was robustly detected despite rs-FC and ODI being measured during separate scans collected on different days. Notably, while the previous anatomical studies suggested that the selectivity of horizontal connections is higher in sites that show a stronger ODI 4 , 7 , 8 , to our knowledge, the link between ODI and the level of rs-FC has never been explicitly tested in animal studies. 3.2. The impact of common input on rs-FC To minimize the impact of visual input on rs-FC, we instructed the participants to keep their eyes closed during resting-state fMRI scans. This approach significantly reduces the impact of common sensory input, which is known to shape rs-FC patterns 45 , particularly in response to structured visual stimuli 11 , 46 – 49 . Nevertheless, one could still argue that relayed spontaneous activity into V1 via feedforward connections from the lateral geniculate nucleus (LGN) might still contribute to the selective rs-FC observed between ODCs with similar ocular polarity. Although input from the LGN may play a role in generating spontaneous activity within V1, this input follows a strict visuotopic organization 50 . However, given that rs-FC selectivity was robustly detected over large cortical distances—well beyond the receptive fields of individual V1 columns—the LGN is unlikely to be the driver of the observed phenomenon 9 – 11 . Moreover, to mitigate spatial blurring due to the fMRI PSF, we excluded vertex pairs closer than 3 mm, thereby further reducing the potential impact of LGN-driven signal correlations. These methodological choices strengthen the interpretation that the rs-FC structure we observed is not trivially inherited from early feedforward input. Instead, horizontal intracortical connections, are well-positioned to support long-range synchrony within V1 51 . In parallel, feedback projections from higher-order visual areas could reinforce or modulate these intrinsic patterns 52 , contributing additional structure to the resting-state signal based on learned or experience-dependent priors. However, considering the weak ocular preference of cortical column in extrastriate visual areas 53 , it is unlikely that the feedback from those regions is the main driver of the selective rs-FC between V1 ODCs. 3.3. Contribution of polysynaptic connectivity and synchronous neuronal activity to rs-FC Given these similarities, it is compelling to link the observed rs-FC patterns to the organization of horizontal fibers. While the spatial organization of rs-FC is thought to be at least partially shaped by underlying anatomical connections 41 , 54 – 56 , rs-FC is not merely a passive reflection of monosynaptic anatomical connectivity 57 . Instead, the existence of polysynaptic connections and synchronous neuronal activity across ODCs play crucial roles in shaping rs-FC. The role of polysynaptic connections in generating rs-FC between sites that are not directly connected has been well established 58 – 61 . This influence is also evident in our findings, which reveal rs-FC between vertex pairs separated by distances of up to 35 mm ( Figure 2 ), far exceeding the reach of a monosynaptic horizontal connection of up to 8mm 1 , 4 , 5 . Crucially, our results provide direct evidence that rs-FC selectivity is preserved even across polysynaptic pathways. Beyond structural connectivity, mounting evidence suggests that rs-FC strength depends on the amplitude of synchronous gamma-band neural activity between the seeded and targeted sites 62 – 67 . Notably, the first direct evidence of such synchronous activity came from simultaneous electrode recordings from V1 columns that share ocular 9 – 11 and orientation 10 , 11 , 48 , 68 preference—even when these columns are spatially distant and have non-overlapping receptive fields. Together, these results strongly support the hypothesis that polysynaptic connections within the visual cortex play a crucial role beyond local visual processing. These connections may be instrumental in integrating information across a large visual field, contributing to a more coherent visual perception. 3.4. Heterogeneity Previous studies have shown that global stimulus configuration is more prominent in the lower than in the upper visual field 29 , 32 , 69 . Neuroimaging research further reveals that dorsal cortical regions, representing the lower visual field, exhibit heightened sensitivity to visual features critical for global configuration encoding, such as lower spatial frequencies 30 . Extending these findings, we showed that rs-FC is significantly stronger in dorsal compared to ventral V1 subregions and in peripheral compared to central V1 subregions. Notably, these effects hold for both intrahemispheric and interhemispheric rs-FC (see also Section 3.5 .). This pattern aligns with the well- established higher sensitivity for lower spatial frequencies in dorsal versus ventral 31 , 70 , 71 and in peripheral vs. central V1 subregions 33 – 35 . Importantly, these findings cannot be explained by variations in gradient-echo BOLD signal intensity across these regions. To clarify, the BOLD signal is generally stronger in the superficial compared to middle and deep cortical depths ( Figure 1C ) and in central compared to peripheral V1 regions. Yet, rs-FC was strongest in the middle cortical depths and in peripheral V1 subregions, indicating that the observed connectivity patterns are not merely a byproduct of BOLD signal strength. Instead, these results highlight the fundamental role of V1 functional architecture in shaping rs-FC, independent of signal intensity biases. While our findings on the overall rs-FC pattern across V1 subregions align with the well-established heterogeneity in visual processing across visual fields, their direct link to visual perception remains uncertain without behavioral measurements. One promising direction is to investigate whether individual differences in mesoscale connectivity are predictive of behavioral performance. For instance, reduced connectivity between ODCs might be associated with poorer visual outcomes, such as diminished contrast sensitivity or stereoacuity. Consistent with this hypothesis, studies in animals have shown that strabismus may alter the overall patterns of functional connectivity 8 , 51 even though there is still a lack of analogous data in humans (but see 72 for preliminary findings). Such findings highlight the potential of rs- fMRI as a non-invasive marker of individual visual function and open avenues for future research into clinical populations with visual impairments. 3.5. Interhemispheric connections and their role Interhemispheric rs-FC is crucial for spatial integration around the vertical meridian. Despite its fundamental role in visual perception, our understanding of V1 interhemispheric FC remains limited. This gap in knowledge stems from earlier findings that, compared to downstream areas like V2, the callosal fibers linking the contralateral V1s are sparse and largely confined to a narrow zone along the V1-V2 border 36 – 38 . Given this and considering that V1 receives extensive feedback from V2—a binocular visual area with robust interhemispheric connections—we assumed that interhemispheric rs-FC in V1 should lack selectivity for ODCs. Our results overturn this assumption. We demonstrate that, particularly in more peripheral and dorsal regions, rs-FC remains significantly stronger between ODCs with alike rather than unalike ocular polarity. This raises two key possibilities: first, despite their limited numbers, V1 callosal connections may still facilitate selective rs-FC between ODCs in opposite hemispheres. Second, the pronounced rs-FC selectivity in superficial cortical layers is consistent with the hypothesis that indirect pathways, mediated by feedback connections from callosally-connected downstream areas, contribute to this phenomenon. Yet, it remains unclear how the eye-of-origin information persists within an indirect, polysynaptic network through binocular visual areas. Resolving this will be essential for understanding the mechanisms underlying interhemispheric functional organization in V1. 3.6. Mapping ODCs based on rs-FC Our findings demonstrate that the organization of ODCs and ODI have a direct influence on the pattern of rs-FC across V1. We also provide evidence suggesting a significant relationship between ODC maps and the rs-FC patterns, making a substantial contribution toward leveraging rs-FC for potential segmentation of the visual cortex at mesoscale levels 12 (see also 20 for a similar approach outside the visual cortex). To validate the specificity of this correspondence, we compared it to correlations with spatially shifted and rotated ODC maps, which preserve the stripy layout and offer a stricter and more interpretable baseline than a randomized ODC map. This said, ODCs alone are unlikely to drive the rs-FC pattern. Prior studies in NHPs suggest that orientation and wavelength (color) preference also modulate functional connectivity between cortical columns 10 – 12 , 68 , 73 , 74 . Consequently, accurate prediction of ODC maps from rs-FC requires disentangling the distinct contributions of these factors on the rs-FC pattern 12 . Future research must refine this approach to fully harness the potential of rs-FC in visual cortex mapping. 3.7. Limitations Although we leveraged an advanced fMRI technique with relatively high spatial resolution, our measurements only indirectly reflect neuronal activity mediated by the BOLD effect and are based on correlations (rather than causal inference). These measurements remain susceptible to partial volume effects and the spatial blurring introduced by the limited spatial resolution and BOLD PSF, which limit the spatial precision of rs-FC estimates 75 – 77 . Tangentially, we attempted to mitigate PSF-related artifacts by excluding vertex pairs closer than 3 mm in distance, thus reducing the impact of local signal spread. In the radial direction—i.e., across cortical depths—PSF effects are more difficult to constrain, particularly given the overall V1 thickness (∼2mm 78 ). This blurring can artificially enhance correlations between layers and obscure depth-specific connectivity patterns. Thus, the true depth-specific differences may in fact be larger than reported, as radial PSF and partial voluming likely attenuate these effects. Future studies employing higher spatial resolution or advanced deconvolution methods may help to better isolate radial connectivity patterns and provide more accurate estimates of depth-specific rs-FC 79 , 80 . Additionally, the inherent sensitivity of our fMRI method to pial veins results in stronger yet partially blurred signals in more superficial cortical depths 16 , 81 – 83 . In addition to the potential influence of pial veins, intracortical vasculature has been highlighted as a relevant factor in mesoscale fMRI analyses 84 . While vascular architecture can indeed introduce spatially structured signal correlations, intracortical vessels generally do not follow the spatial organization of mesoscale neuronal ensembles (e.g. ODCs in V1) 85 . Therefore, while we cannot fully exclude their influence, it is unlikely that they would selectively bias rs-FC between ocular dominance columns with alike eye preference. The observed selectivity in our results is more plausibly attributed to underlying neuronal connectivity rather than shared vascular dynamics. 4. Conclusion The introduction of higher-performance MRI scanners and technological advancements in hardware and software design enable us to now resolve neuronal responses at the fine spatial scale of cortical columns. This study demonstrates the power of high-resolution fMRI in capturing one of the most fine-grained features of cortical organization including the selective connectivity between ODCs. Given that this selectivity produces only a modest 4% to 15% difference in horizontal fiber density, our findings underscore the remarkable sensitivity of high-resolution fMRI in revealing the mesoscale functional organization of the human visual system. Thanks to the short duration of a resting-state scan session (<30 mins) compared to the time necessary for localizing ODCs based on dichoptic visual stimulation (∼3-4 hours 23 ) and the absence of cognitive load on the participant, the mesoscale rs-FC measurements present a potentially powerful tool for clinicians to assess visual system integrity at the mesoscale level. It may even enable such studies in patients where conventional methods of ODC mapping are impractical or impossible (e.g., in monocularly blind individuals). 5. Methods 5.1. Participants 12 subjects (4 females), aged 25-45 years, participated in this study. All participants had radiologically intact brains, intact vision, and no history of neuropsychological disorders. All experimental procedures conformed to the guidelines of the National Institutes of Health (NIH) and were approved by the Mass General Brigham Institutional Review Board. Written informed consent was obtained from all participants prior to all experiments. 5.2. Procedure Each participant was scanned in an ultra-high-field (7T) MRI scanner to collect resting-state data. To test the reproducibility of the findings, four of the participants—selected based on their availability—were rescanned on a different day using a different 7T scanner. For each participant, ODC and retinotopic maps were obtained from a separate set of scans 23 . Structural scans were acquired on a 3T scanner. 5.2.1. Resting-state measurements To collect spontaneous fMRI activity, each participant was scanned for 3 to 4 runs (512 s per run), except for one individual who was scanned for just one run due to time limitations. Participants were instructed to keep their eyes closed for the duration of the run but not to fall asleep. We communicated with participants between runs to ensure they remained awake during the scan. 5.2.2. ODC mapping ODC mapping was conducted on different days from the resting-state session. Details of the procedure are already described comprehensively in our recent publication 23 . Briefly, the ODC map was determined by contrasting the fMRI activity evoked in response to stimulation of the dominant eye (DE) vs. non- dominant eye (NDE) in a block design. The stimuli consisted of sparse (5%) horizontally moving random dots (-0.22° to 0.22°; 0.3 Hz), with red (50% of blocks) and green (remaining blocks) dots (0.09° × 0.09°; 56 cd/m²), presented against a black background. Participants viewed the stimuli through custom-made anaglyph spectacles (with red and green filters) mounted on the head coil. The stimuli spanned 20° × 26° of the visual field. Filter laterality (i.e., red-left vs. red-right) was counterbalanced between sessions and across participants. 5.2.3. Retinotopic mapping The borders of retinotopic visual areas were defined by mapping the horizontal vs. vertical meridians 86 , 87 . Stimuli consisted of flashing radial checkerboards within retinotopically defined wedge-shaped apertures, extending radially along either the horizontal or vertical meridians (polar angle = 30°) and presented against a gray background. Further details are provided in our previous publication 23 . 5.3. Imaging 5.3.1. Main experiments For the main resting-state fMRI measurements (but see also the reproducibility tests), and also for ODC and retinotopic mapping, we utilized a 7T whole-body Siemens scanner equipped with SC72 body gradients (70 mT/m gradient strength and max slew rate of 200 T/m/s), a custom-built 32-channel helmet radio-frequency (RF) receive coil array, and a birdcage volume RF transmit coil. Resting-state data were collected using single-shot 2D gradient-echo EPI, with 1.0 mm isotropic nominal voxel size, and the following protocol parameter values: TR=4000 ms, TE=28 ms, excitation flip angle=85°, field of view (FOV) matrix=192×192, BW=1184 Hz/pixel, echo-spacing=1 ms, 7/8 phase partial Fourier, 59 oblique- coronal slices covering occipital lobe, acceleration factor R=4 with GRAPPA reconstruction and FLEET- ACS data 88 with a 10° flip angle. For ODC and retinotopic mapping, we used a similar EPI protocol with TR=3000 ms, flip angle=78°, and 44 oblique-coronal slices. 5.3.2. Reproducibility tests To test the reproducibility of the resting-state results, 4 participants were scanned again (on a different day) in a Siemens MAGNETOM Terra whole-body MRI (80 mT/m gradient strength and max slew rate of 200 T/m/s). In contrast to the other 7T scanner (see above), this scanner provides passive and active shimming up to 3 rd order and is equipped with an in-house-built head RF coil system that use a quadrature birdcage transmit coil along with a close-fitting 64-ch RF receive array 89 . In these sessions, we used single-shot 2D gradient-echo EPI, with 1.0 mm isotropic voxel size, and the following protocol parameter values: TR=4000 ms, TE=25 ms, flip angle=80°, FOV matrix=180×180, BW=1158 Hz/Px, echo- spacing=1 ms, 7/8 phase partial Fourier, 66 oblique-coronal slices covering occipital cortex, and FLEET- ACS data 88 with 10° flip angle. We also used R=4 dual-polarity GRAPPA reconstruction for improved Nyquist ghost correction 90 . 5.3.3. Structural scans At ultra-high field (7T), T1-weighted anatomical brain images typically suffer from image contrast variations (bias field) over the FOV due to the combined effects of: (a) transmit RF field inhomogeneity caused by dielectric wavelength effects within the head 91 and (b) disrupted inversion pulses in areas where tissue susceptibility perturbs the static magnetic field 92 . To avoid this issue, we used anatomical images acquired at 3T on a Siemens TimTrio whole-body scanner, equipped with the standard vendor- supplied 32-channel RF head coil array, which provides relatively uniform image contrast over the full FOV. We used a 3D T1-weighted MPRAGE sequence with protocol parameter values: TR=2530 ms, TE=3.39 ms, TI=1100 ms, flip angle=7°, bandwidth=200 Hz/pix, echo spacing=8.2 ms, voxel size=1.0×1.0×1.0 mm³, and FOV=256×256×170 mm³. 5.4. Data analysis Functional and anatomical MRI data were pre-processed and analyzed using FreeSurfer (version 7.4.0; http://surfer.nmr.mgh.harvard.edu ) 93 and in-house MATLAB and Python code. As described in the following sections in detail, each subject was analyzed entirely in their native anatomical space to preserve the fine-grained structure of ODCs and cortical lamination. This native-space processing was maintained throughout all individual analyses. An overview of the processing steps applied to the ODC mapping and resting-state data is presented in Supplementary Figure 4. 5.4.1. Anatomical data 5.4.1.1. Defining cortical surfaces Based on each participant’s anatomical data, cortical surfaces were reconstructed using FreeSurfer’s recon-all pipeline, which includes intensity normalization, skull stripping, white and pial surface extraction, and topological correction using a deformable surface model constrained by neuroanatomical priors 94 . In this process, the standard pial surface was defined as the boundary between the gray matter (GM) and the surrounding cerebrospinal fluid (CSF), while the white matter (WM) surface was generated as the interface between the WM and GM. Additionally, nine intermediate surfaces were created using the equivolume approach, spaced at intervals of 10% of the cortical thickness. Finally, surfaces were flattened using FreeSurfer’s mris_flatten routine. 5.4.1.2. Defining cortical depths To examine the impact of cortical depth on rs-FC, the GM was divided into deep, middle, and superficial cortical depths ( Figure 1A ). The deep cortical depth was defined as the region between the WM-GM interface and the two adjacent surfaces above it. The superficial cortical depth was defined as the region between the GM-CSF interface and the two adjacent surfaces below it. The middle cortical depth included the three intermediate reconstructed surfaces. Notably, this approach left one intermediate surface gap between these depth levels to improve the separation of the three segments. To enhance the co- registration of functional and structural scans, all surfaces were upsampled to a 0.5 mm resolution using butterfly subdivision in mris_mesh_subdivide as implemented in FreeSurfer and outlined in our previously introduced pipeline 95 . 5.4.2. Functional data 5.4.2.1. Preprocessing All functional data were preprocessed with a pipeline for ultra-high-resolution fMRI which mostly uses FreeSurfer commands 95 . Briefly, the collected functional data were first upsampled to 0.5 mm isotropic resolution using trilinear interpolation using FreeSurfer’s mri_convert (see also 95 ). For each participant, functional data from each run were rigidly aligned (with 6 degrees of freedom) to their structural scan using rigid Boundary-Based Registration (Supplementary Figure 5) 96 . Head motion covariates were then derived for six directions. After motion correction, intensity normalization was applied to the fMRI data collected during ODC mapping, but not to the resting-state fMRI data (Supplementary Figure 4). Resting-state data were analyzed with an adjusted preprocessing pipeline designed to preserve as much of the original signal as possible while minimizing the impact of noise sources, which is particularly important when examining mesoscale functional connectivity. Specifically, before upsampling to 0.5 mm, the data were denoised using the NORDIC approach 97 , which was reported to effectively suppresses thermal noise in fMRI data 98 . Thermal noise poses a particular challenge for resting-state analyses, as it reduces temporal signal-to-noise ratio (tSNR) and can mask meaningful low- frequency BOLD fluctuations, thereby compromising the reliability of functional connectivity estimates. MP-PCA denoising, as implemented in NORDIC, has been shown to approximately double the tSNR and substantially increase correlation coefficients, with only a minimal increase in spatial smoothness (7%) 99 . After slice timing correction 93 , the data were detrended for quadratic trends and high-pass filtered using a cut-off frequency of 0.01 Hz 100 . To preserve spatial resolution on the surface, no tangential spatial smoothing was applied to the fMRI data. Instead, radial (intracortical) smoothing—perpendicular to the cortex and within cortical columns— was used 101 . Considering the cortical thickness of V1 (∼2mm 78 ), the 1 mm isotropic voxel size, and the fact that voxel boundaries do not precisely align with the gray–white matter interface, each cortical depth level likely receives signal contributions from approximately 1–2 voxels. For ODC and retinotopic mapping, to minimize the blurring impact of pial veins on functional map resolution 16 , 81 – 83 , radial smoothing was restricted to the deep cortical depth (as defined in 5.4.1.2). For rs-FC analysis, radial smoothing was applied either to the entire GM ( sections 2.1 - 2.2 ) or separately to the deep, middle, and superficial cortical depths ( sections 2.3 - 2.6 ). Notably, previous studies in NHPs 102 and humans 23 have shown that the ODC maps remain mostly unchanged across cortical depth. The stability and consistency of ODC maps derived from fMRI (across sessions and depth levels) are already demonstrated in our previous publications 16 , 23 . Notably, we avoided applying any method of spatial distortion correction to the fMRI data because EPI distortions are typically small in the occipital lobe 103 . This point is also demonstrated in Supplementary Figure 5, which highlights the alignment of EPI to anatomical scans in the occipital lobe, including V1. Moreover, distortion correction techniques generally cause spatial smoothing 95 . Considering that spatial resolution is critical for resolving fine-scale structures such as ODCs, we avoided such smoothing-inducing corrections. 5.4.2.2. Resting-state functional connectivity analysis The rs-FC between pairs of vertices within the flattened cortical surface mesh of V1 was measured by estimating Pearson’s partial correlation coefficients of their resting-state fMRI time series. For our main analysis, this was conducted separately within each hemisphere; however, we also evaluated rs-FC between hemispheres. We did not include correlations between vertex pairs located at different cortical depths considering the organization of horizontal fibers (i.e., primarily parallel to the cortical surface). To reduce the impact of physiological and motion-related confounds, we included the mean white matter signal and six head motion parameters as covariates in the partial correlation analyses. Moreover, to remove the impact of fMRI signal propagation, we excluded all vertex pairs with an inter-vertex distance less than 3 mm. This threshold is in line with prior estimates of the fMRI PSF at 7T 75 – 77 and is intended to reduce the influence of non-neuronal spatial autocorrelation. This approach allowed us to estimate functional connectivity while accounting for structured sources of physiological noise. Notably, we avoided using other more recently introduced methods that rely on specific assumptions about the correlated spontaneous activity evoked between voxels because they have not been tested at mesoscale levels104,105. To mitigate the influence of noise at the single-vertex level, we averaged correlation coefficients across runs, improving the signal-to-noise ratio. For our next analyses, which required linear statistical operations on correlation values, we used Fisher z -transformed correlation values that are expected to have a distribution closer to normal than Pearson’s correlation coefficients do and are thus more suitable for parametric calculations. 5.4.2.3. ODC and retinotopic mapping Details of the data analysis for mapping ODCs and retinotopic visual areas are described in our recent publication 23 . Briefly, after preprocessing (see above), a standard hemodynamic response model based on a gamma function was fitted to the signals from deep cortical depths to estimate the amplitude of the blood-oxygen-level-dependent (BOLD) response. To generate differential maps of DE vs. NDE activation, a vertex-wise general linear model was applied by computing the contrast of interest. The resulting beta and significance maps were projected onto each participant’s anatomical volumes and reconstructed cortical surfaces ( Figure 1B ). The reproducibility of the maps across scan sessions was tested and proved for all individual participants 16 , 23 . As noted above, only data from deep cortical depths were used for ODC and retinotopic mapping. 5.5. Statistical analysis As a second-order statistical analysis (see also Section 5.4.2 .2), we assessed the significance of the impact of the following independent parameters on rs-FC: (1) relative ocular polarity, (2) distance between vertices, (3) level of ocular preference (ocular dominance index (ODI)), (4) cortical depth, and (5) region of interest (ROI). These parameters were calculated as follows. 5.5.1. Distance between vertex pairs The distance between vertices in each pair was calculated as the Euclidean distance on the flattened cortical surface. To minimize the bias in the rs-FC resulting from the effect of the fMRI PSF, we excluded vertex pairs located closely together (<3 mm) 75 – 77 . Since vertex pairs with alike ocular polarity (see below) generally were closer together than unalike pairs, we used random subsampling to balance the distance distributions between alike and unalike vertex pairs. To that end, we divided the distance distribution of all vertex pairs into 100 quantiles, for each of which, we then determined the counts of alike and unalike vertex pairs and randomly discarded excess pairs from the larger group to equalize their numbers. Because V1 varied in size across participants, vertex pairs for each participant were divided into 10 distance quantiles, with 1 referring to the shortest distances and 10 referring to the longest distance. The quantile rank (rather than the absolute distance) was used for further statistical analysis. The median distance for each quantile is shown in Supplementary Figure 1. The distance between vertex pairs was then used as a numerical parameter in the statistical tests. 5.5.2. Relative ocular polarity of vertices To assess the effect of relative ocular polarity on rs-FC, a pair of vertices was classified as alike or unalike when the ocular preferences of the two vertices were the same (i.e. both vertices responded preferentially to the same eye) or opposite, respectively. The relative ocular polarity was then used as a categorical parameter in the statistical tests. Notably, the ocular preference of each vertex was determined based on the unthresholded ODC map ( Figure 1B ). Consistent with our previous findings in a smaller group of participants 23 , the number of vertices responding preferentially to DE and NDE stimulation, as well as the amplitude of ocular-preferring activity (i.e., the response to the preferring vs. non-preferring eye stimulation), did not differ significantly (Supplementary Figure 2). This result indicates that both eyes were represented equally within V1. 5.5.3. Ocular dominance index (ODI) To assess the impact of ocular preference level on rs-FC, we used a method similar to the one introduced by Hubel and Wiesel, 1962 and later adopted by others 106 – 108 . Specifically, for each vertex, the level of ocular preference — i.e., the estimated beta value for differential response to preferred vs. non-preferred eye stimulation —was measured during ODC mapping. For each subject, these beta values were divided into 10 quantiles, with 1 representing the weakest ocular preference and 10 representing the strongest. The quantile rank (rather than the absolute ocular preference) was then used as the ODI. The median beta value for each ODI is shown in Supplementary Figure 6. The effect of ODI on rs-FC between vertex pairs was tested by grouping vertices into 5 cumulative categories: (1) vertex pairs from quantiles 1 and 2, (2) quantiles 1 to 4, (3) quantiles 1 to 6, (4) quantiles 1 to 8, and (5) quantiles 1 to 10. 5.5.4. Cortical depths To test the impact of cortical depth on the level of rs-FC, the rs-FC data analysis was conducted separately for the deep, middle, and superficial depths of the GM. The cortical depth was then used as a categorical parameter in the statistical tests. 5.5.5. ROI In the main analysis (sections 2.1 to 2.3), the ROI included the stimulated portion of V1 (radius<10°) except for a small dorsal region along the V1-V2 border ( Figure 1B ) where ODCs could not be identified due to occlusion of the ventral visual field by the nose 23 . To evaluate the heterogeneity of rs-FC within V1, we manually divided the stimulated part of V1 into size-matched dorsal and ventral as well as central (radius<4°) and peripheral (4°<radius<10°) subregions and compared the rs-FC within these subdivisions (see sections 2.4 and 2.5 ). The ROI was then used as a categorical parameter in the statistical tests. 5.5.6. Statistical tests 5.5.6.1. Repeated measures (RM) ANOVA RM-ANOVA was primarily used to test the impact of independent parameters (see above) on the measured rs-FC between the vertex pairs. Specifically, we first tested the impact of relative ocular polarity (alike vs. unalike) and distance using a two-way RM-ANOVA ( Section 2.1 ). Second, for a subset of 4 participants, the reproducibility of the results was tested by adding the session number (categorical: 1 st vs. 2 nd ) to the independent parameters (i.e., distance and relative ocular polarity) and applying a three- way RM-ANOVA ( Section 2.2 ). The impact of ODI (5 cumulative categories, see Section 5.5.3 ) and cortical depth, along with relative ocular polarity, was tested by applying a three-way RM-ANOVA ( Section 2.3 ). Lastly, heterogeneity of intrahemispheric and interhemispheric rs-FC measurements was tested by including V1 subregions (ROIs) in the list of independent parameters and applying a four-way RM-ANOVA ( sections 2.3 and 2.4 ). Notably, in all tests except for interhemispheric rs-FC measurements, the rs-FC values from the two hemispheres were averaged to improve the signal-to-noise ratio. Considering that RM-ANOVA is particularly susceptible to violations of the sphericity assumption, caused by the correlation between measured values and unequal variance of differences between experimental conditions, results were corrected, when necessary (as determined using a Mauchly test), using the Greenhouse-Geisser method. A p-value of <0.05 was considered statistically significant. 5.5.6.2. Linear mixed effects (LME) model To ensure that overlap between the cumulative ODI categories did not skew the estimated significance of our results, we conducted a supplementary analysis based on using an LME model to assess the relationship between rs-FC and relative ocular polarity, ODI, and cortical depth accounting for subject- level variability (random intercepts by subject). In this analysis, rather than using the cumulative ODI categories, we used the product of BQ 1 and BQ 2 (i.e. BQ 12 ), where BQ 1 and BQ 2 represent the ODIs of the first and second vertices in a pair. Since this created a symmetric matrix of 10×10 categories, we excluded the redundant upper triangle, resulting in 55 joint ODI categories, which were then used in an LME model (see Section 5.5.6 .2). By using 55 non-overlapping joint ODI brackets, we ensured that the significance of the results was not overestimated due to datapoint (vertex pair) interdependencies. To apply the model the MATLAB command fitlme was used according to the following formula: rsFC ∼ BQ 12 + CD + ROP + BQ 12 * CD + BQ 12 * ROP + CD * ROP + BQ 12 * CD * ROP + (1|Subject), where BQ 12 represents the product of the ODIs of the first and second vertices in a pair, CD represents cortical depth, and ROP is relative ocular polarity. This approach enabled us to check for the consistency of the results obtained by ANOVA and LME. Notably, in contrast to ANOVA, LME takes the hierarchical structure of the data into account. This method further allowed us to quantify the independent contributions of the main factors and their interactions while controlling for variability across subjects (see Section 2.3 ). A p-value of <0.05 was considered statistically significant. 5.5.6. Correlation between rs-FC and ODC maps To evaluate how well rs-FC reflects the underlying ODC map, we conducted a correlation analysis. We selected 1000 seed vertices at random within V1 (far enough from the V1 border) and correlated each seed vertex’s resting-state fMRI time-series with the rest of the vertices in a 2D donut-shaped region (ring) centered at that seed vertex, generating 1000 2D rs-FC maps. To investigate the impact of distance, we conducted the analysis multiple times within concentric rings around each seed vertex with varying outer radius ranging from 4 to 10 mm. Similarly to our main analyses, we excluded nearby target vertices falling within an inner radius of 3 mm to minimize potential biases caused by the PSF of the BOLD signal 75 – 77 . Within each ring, we calculated the correlation between the rs-FC map and the underlying ODC map, which constitutes the correlation value under the alternative hypothesis (H1) that rs-FC and ODC maps are related. For each ring size, the resultant correlation values were averaged across the 1000 sampled sites. To determine the statistical significance of our results, we compared these correlations to null- hypothesis (H0) distributions (correlations between unrelated rs-FC and ODC maps) generated by correlating the rs-FC values in each ring with the ODC differential map in a spatially shifted (misaligned) ring (H0 a ) as well as a 180° rotated ring (H0 b ). With this approach, we tested whether the observed correlations were meaningful and not solely due to the inherent spatial structure of the ODC map. To ensure that interpolation and resampling artifacts affected both the null and alternative tests similarly, in all experiments, we resampled the rs-FC and ODC map values with nearest-neighbor interpolation on a set of points generated via 2D quasi-random Halton sampling 109 , 110 and distributed with a density of 10 points per mm². 5.6. Data availability The data and code will be shared upon reasonable request. Author’s contribution Conceptualization, formal analysis, visualization, writing – original draft: M.E.S., I.A. and S.N. Investigation, writing – review & editing: M.E.S., I.A., J.S., B.B., Y.C., W.S.H., E.K., N.W. and S.N. Methodology, software: M.E.S., I.A., J.S., B.B., Y.C., W.S.H. and S.N. Data curation, validation: M.E.S. and S.N. Supervision: E.K., N.W. and S.N. Funding acquisition: N.W. and S.N. Resources, project administration: S.N. Competing interests The authors declare no competing interests. Acknowledgment This work was supported by the National Eye Institute (R01EY030434 and R01EY029713), the National Institute on Aging (RF1AG068261), the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project no. 347592254 (WE 5046/4-2 and/or KI 1337/2-2), the Federal Ministry of Education and Research (BMBF) under support code 01ED2210, as well as by the MGH/HST Athinoula A. Martinos Center for Biomedical Imaging and the Max Planck School of Cognition. Crucial imaging resources were made available by the NIH Shared Instrumentation Grant S10RR019371. We thank Azma Mareyam for helping with hardware maintenance as well as John Kirsch and Daniel Park for their support in raw data export during this study. We thank Amanda Nabasaliza, Sarala N. Malladi, and Bryan Kennedy for their help with recruitment and data collection. We also thank Drs. Antony Morland and Jingyuan Chen for their comments. Funder Information Declared National Eye Institute , R01EY030434 , R01EY029713 National Institute on Aging , RF1AG068261 Deutsche Forschungsgemeinschaft, https://ror.org/018mejw64 , WE 5046/4-2 , KI 1337/2-2 Federal Ministry of Education and Research , 01ED2210 National Institutes of Health, https://ror.org/01cwqze88 , S10RR019371 Massachusetts General Hospital, https://ror.org/002pd6e78 Max Planck School of Cognition Footnotes ↵ + Co-first authors Updated paper structure: added a summary of key findings to the introduction and a methods overview paragraph to the results section. Abbreviations ODCs Ocular Dominance Columns FMRI Functional Magnetic Resonance Imaging RS-FC Resting-state Functional Connectivity DE Dominant Eye NDE Non-Dominant Eye ODI Ocular Dominance Index ROI Region Of Interest RM-ANOVA Repeated Measures Analysis of Variance References 1. ↵ Fisken , R. A. , Garey , L. J. & Powell , T. P . The intrinsic, association and commissural connections of area 17 on the visual cortex . Philos. Trans. R. Soc. Lond. B. Biol. Sci . 272 , 487 – 536 ( 1975 ). OpenUrl CrossRef PubMed Web of Science 2. Gilbert , C. D. & Wiesel , T. N . Morphology and intracortical projections of functionally characterised neurones in the cat visual cortex . Nature 280 , 120 – 125 ( 1979 ). OpenUrl CrossRef PubMed Web of Science 3. ↵ Rockland , K. S. & Lund , J. S . Intrinsic laminar lattice connections in primate visual cortex . J. Comp. Neurol . 216 , 303 – 318 ( 1983 ). OpenUrl CrossRef PubMed Web of Science 4. ↵ Yoshioka , T. , Blasdel , G. G. , Levitt , J. B. & Lund , J. S . Relation between patterns of intrinsic lateral connectivity, ocular dominance, and cytochrome oxidase-reactive regions in macaque monkey striate cortex . Cereb. Cortex 6 , 297 – 310 ( 1996 ). OpenUrl CrossRef PubMed Web of Science 5. ↵ Angelucci , A. et al. Circuits for local and global signal integration in primary visual cortex . J. Neurosci. Off. J. Soc. Neurosci . 22 , 8633 – 8646 ( 2002 ). OpenUrl Abstract / FREE Full Text 6. ↵ Lund , J. S. , Angelucci , A. & Bressloff , P. C . Anatomical substrates for functional columns in macaque monkey primary visual cortex . Cereb. Cortex 13 , 15 – 24 ( 2003 ). OpenUrl CrossRef PubMed Web of Science 7. ↵ Malach , R. , Amir , Y. , Harel , M. & Grinvald , A. Relationship between intrinsic connections and functional architecture revealed by optical imaging and in vivo targeted biocytin injections in primate striate cortex . Proc. Natl. Acad. Sci. U. S. A . 90 , 10469 – 10473 ( 1993 ). OpenUrl Abstract / FREE Full Text 8. ↵ Tychsen , L. , Wong , A. M.-F. & Burkhalter , A . Paucity of horizontal connections for binocular vision in V1 of naturally strabismic macaques: Cytochrome oxidase compartment specificity . J. Comp. Neurol . 474 , 261 – 275 ( 2004 ). OpenUrl CrossRef PubMed Web of Science 9. ↵ Chiu , C. & Weliky , M . Relationship of correlated spontaneous activity to functional ocular dominance columns in the developing visual cortex . Neuron 35 , 1123 – 1134 ( 2002 ). OpenUrl CrossRef PubMed Web of Science 10. ↵ Omer , D. B. , Fekete , T. , Ulchin , Y. , Hildesheim , R. & Grinvald , A . Dynamic patterns of spontaneous ongoing activity in the visual cortex of anesthetized and awake monkeys are different . Cereb. Cortex 29 , 1291 – 1304 ( 2019 ). OpenUrl CrossRef PubMed 11. ↵ Ts’o , D. Y. , Gilbert , C. D. & Wiesel , T. N . Relationships between horizontal interactions and functional architecture in cat striate cortex as revealed by cross-correlation analysis . J. Neurosci . 6 , 1160 – 1170 ( 1986 ). OpenUrl Abstract / FREE Full Text 12. ↵ Cai , X. et al. Mesoscale functional connectivity in macaque visual areas . NeuroImage 271 , 120019 ( 2023 ). 13. ↵ Dawson , D. A. et al. Partial correlation-based retinotopically organized resting-state functional connectivity within and between areas of the visual cortex reflects more than cortical distance . Brain Connect . 6 , 57 – 75 ( 2016 ). OpenUrl CrossRef PubMed 14. ↵ Mendola , J. D. , Lam , J. , Rosenstein , M. , Lewis , L. B. & Shmuel , A . Partial correlation analysis reveals abnormal retinotopically organized functional connectivity of visual areas in amblyopia . NeuroImage Clin . 18 , 192 – 201 ( 2018 ). OpenUrl PubMed 15. ↵ Genç , E. , Schölvinck , M. L. , Bergmann , J. , Singer , W. & Kohler , A . Functional connectivity patterns of visual cortex reflect its anatomical organization . Cereb. Cortex 26 , 3719 – 3731 ( 2016 ). OpenUrl CrossRef PubMed 16. ↵ Nasr , S. , Polimeni , J. R. & Tootell , R. B . Interdigitated color- and disparity-selective columns within human visual cortical areas V2 and V3 . J Neurosci 36 , 1841 – 57 ( 2016 ). OpenUrl Abstract / FREE Full Text 17. Kennedy , B. , Bex , P. , Hunter , D. G. & Nasr , S . Two fine-scale channels for encoding motion and stereopsis within the human magnocellular stream . Prog. Neurobiol . 220 , 102374 ( 2023 ). 18. Ai , H. , Lin , W. , Liu , C. , Chen , N. & Zhang , P . Mesoscale functional organization and connectivity of color, disparity, and naturalistic texture in human second visual area . eLife 13 , ( 2025 ). 19. ↵ Tootell , R. B. & Nasr , S . Scotopic vision is selectively processed in thick-type columns in human extrastriate cortex . Cereb. Cortex 31 , 1163 – 1181 ( 2021 ). OpenUrl CrossRef PubMed 20. ↵ Huber , L. et al. Layer-dependent functional connectivity methods . Prog. Neurobiol . 207 , 101835 ( 2021 ). 21. ↵ Cheng , K. , Waggoner , R. A. & Tanaka , K . Human ocular dominance columns as revealed by high-field functional magnetic resonance imaging . Neuron 32 , 359 – 374 ( 2001 ). OpenUrl CrossRef PubMed Web of Science 22. Menon , R. S. , Ogawa , S. , Strupp , J. P. & Uğurbil , K . Ocular dominance in human V1 demonstrated by functional magnetic resonance imaging . J. Neurophysiol . 77 , 2780 – 2787 ( 1997 ). OpenUrl CrossRef PubMed Web of Science 23. ↵ Nasr , S. et al. Differential impacts of strabismic and anisometropic amblyopia on the mesoscale functional organization of the human visual cortex . J. Neurosci . ( 2025 ) doi: 10.1523/JNEUROSCI.0745-24.2024 . OpenUrl Abstract / FREE Full Text 24. Yacoub , E. , Shmuel , A. , Logothetis , N. & Uğurbil , K . Robust detection of ocular dominance columns in humans using Hahn Spin Echo BOLD functional MRI at 7 Tesla . NeuroImage 37 , 1161 – 1177 ( 2007 ). OpenUrl CrossRef PubMed Web of Science 25. Zaretskaya , N. et al. Eye-selective fMRI activity in human primary visual cortex: Comparison between 3 T and 9.4 T, and effects across cortical depth . NeuroImage 220 , 117078 ( 2020 ). 26. Akbari , A. , Gati , J. S. , Zeman , P. , Liem , B. & Menon , R. S. Layer dependence of monocular and binocular responses in human ocular dominance columns at 7T using VASO and BOLD . bioRxiv ( 2023 ) doi: 10.1101/2023.04.06.535924 . OpenUrl Abstract / FREE Full Text 27. ↵ Haenelt , D. , et al. Decoding of columnar-level organization across cortical depth using BOLD- and CBV-fMRI at 7T . bioRxiv ( 2023 ) doi: 10.1101/2023.09.28.560016 . OpenUrl Abstract / FREE Full Text 28. ↵ Aganj , I. & Nasr , S . Selective functional connectivity between ocular dominance columns in the primary visual cortex . in Ophthalmic Medical Image Analysis (eds. Bhavna , A. , Chen , H. , Fang , H. , Fu , H. & Lee , C. S. ) 1 – 10 ( Springer Nature Switzerland, Cham , 2024 ). doi: 10.1007/978-3-031-73119-8_1 . OpenUrl CrossRef 29. ↵ Christman , S. D . Local-global processing in the upper versus lower visual fields . Bull. Psychon. Soc . 31 , 275 – 278 ( 1993 ). OpenUrl CrossRef 30. ↵ Nasr , S. & Tootell , R. B. H . Asymmetries in Global Perception Are Represented in Near- versus Far-Preferring Clusters in Human Visual Cortex . J. Neurosci . 40 , 355 – 368 ( 2020 ). OpenUrl Abstract / FREE Full Text 31. ↵ Niebauer , C. L. & Christman , S. D . Upper and lower visual field differences in categorical and coordinate judgments . Psychon. Bull. Rev . 5 , 147 – 151 ( 1998 ). OpenUrl CrossRef 32. ↵ Previc , F. H . Functional specialization in the lower and upper visual fields in humans: Its ecological origins and neurophysiological implications . Behav. Brain Sci . 13 , 519 – 542 ( 1990 ). OpenUrl CrossRef Web of Science 33. ↵ De Valois , R. L. , Albrecht , D. G. & Thorell , L. G. Spatial frequency selectivity of cells in macaque visual cortex . Vision Res . 22 , 545 – 559 ( 1982 ). OpenUrl CrossRef PubMed Web of Science 33. Schiller , P. H. , Finlay , B. L. & Volman , S. F . Quantitative studies of single-cell properties in monkey striate cortex . III. Spatial frequency. J. Neurophysiol . 39 , 1334 – 1351 ( 1976 ). OpenUrl PubMed 34. ↵ Tootell , R. B. , Silverman , M. S. , Hamilton , S. L. , Switkes , E. & De Valois , R. L . Functional anatomy of macaque striate cortex . V. Spatial frequency. J. Neurosci. Off. J. Soc. Neurosci . 8 , 1610 – 1624 ( 1988 ). OpenUrl CrossRef 35. ↵ Kennedy , H. , Dehay , C. & Bullier , J . Organization of the callosal connections of visual areas V1 and V2 in the macaque monkey . J. Comp. Neurol . 247 , 398 – 415 ( 1986 ). OpenUrl CrossRef PubMed Web of Science 36. Caspers , S. et al. Target sites for transcallosal fibers in human visual cortex – A combined diffusion and polarized light imaging study . Cortex 72 , 40 – 53 ( 2015 ). OpenUrl CrossRef PubMed 37. ↵ Clarke , S. & Miklossy , J . Occipital cortex in man: organization of callosal connections, related myelo- and cytoarchitecture, and putative boundaries of functional visual areas . J. Comp. Neurol . 298 , 188 – 214 ( 1990 ). OpenUrl CrossRef PubMed Web of Science 38. ↵ Adams , D. L. , Sincich , L. C. & Horton , J. C . Complete pattern of ocular dominance columns in human primary visual cortex . J. Neurosci . 27 , 10391 – 10403 ( 2007 ). OpenUrl Abstract / FREE Full Text 39. ↵ Biswal , B. , Zerrin Yetkin , F. , Haughton , V. M. & Hyde , J. S . Functional connectivity in the motor cortex of resting human brain using echo-planar mri . Magn. Reson. Med . 34 , 537 – 541 ( 1995 ). OpenUrl CrossRef PubMed Web of Science 40. ↵ Fox , M. D. & Raichle , M. E . Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging . Nat. Rev. Neurosci . 8 , 700 – 711 ( 2007 ). OpenUrl CrossRef PubMed Web of Science 41. ↵ Power , J. D. et al. Functional network organization of the human brain . Neuron 72 , 665 – 678 ( 2011 ). OpenUrl CrossRef PubMed Web of Science 42. Yeo , B. T. et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity . J. Neurophysiol . 106 , 1125 – 1165 ( 2011 ). OpenUrl CrossRef PubMed Web of Science 43. ↵ Gordon , E. M. et al. Generation and evaluation of a cortical area parcellation from resting-state correlations . Cereb. Cortex 26 , 288 – 303 ( 2016 ). OpenUrl CrossRef PubMed 44. ↵ Adams , P. et al. Corticothalamic interactions in the transfer of visual information . Philos. Trans. R. Soc. Lond. B. Biol. Sci . 357 , 1739 – 1752 ( 2002 ). OpenUrl CrossRef PubMed Web of Science 45. ↵ Nasr , S. , Kleinfeld , D. & Polimeni , J. R . The global configuration of visual stimuli alters co-fluctuations of cross-hemispheric human brain activity . J. Neurosci . 41 , 9756 – 9766 ( 2021 ). OpenUrl Abstract / FREE Full Text 46. Ontivero-Ortega , M. et al. Intra-V1 functional networks and classification of observed stimuli. Front . Neuroinformatics 18 , 1080173 ( 2024 ). 47. ↵ Engel , A. K. , König , P. , Gray , C. M. & Singer , W . Stimulus-dependent neuronal oscillations in cat visual cortex: Inter-columnar interaction as determined by cross-correlation analysis . Eur. J. Neurosci . 2 , 588 – 606 ( 1990 ). OpenUrl CrossRef PubMed Web of Science 48. ↵ Gray , C. M. , König , P. , Engel , A. K. & Singer , W . Oscillatory responses in cat visual cortex exhibit inter-columnar synchronization which reflects global stimulus properties . Nature 338 , 334 – 337 ( 1989 ). OpenUrl CrossRef PubMed Web of Science 49. ↵ Blasdel , G. G. & Lund , J. S . Termination of afferent axons in macaque striate cortex . J. Neurosci . 3 , 1389 – 1413 ( 1983 ). OpenUrl Abstract / FREE Full Text 50. ↵ Löwel , S. & Singer , W . Selection of intrinsic horizontal connections in the visual cortex by correlated neuronal activity . Science 255 , 209 – 212 ( 1992 ). OpenUrl Abstract / FREE Full Text 51. ↵ Angelucci , A. & Bressloff , P. C . Contribution of feedforward, lateral and feedback connections to the classical receptive field center and extra-classical receptive field surround of primate V1 neurons . Prog. Brain Res . 154 , 93 – 120 ( 2006 ). OpenUrl CrossRef PubMed Web of Science 52. ↵ Poggio , G. F. & Fischer , B . Binocular interaction and depth sensitivity in striate and prestriate cortex of behaving rhesus monkey . J. Neurophysiol . 40 , 1392 – 1405 ( 1977 ). OpenUrl CrossRef PubMed Web of Science 53. ↵ De Luca , M. , Beckmann , C. F. , De Stefano , N. , Matthews , P. M. & Smith , S. M . fMRI resting state networks define distinct modes of long-distance interactions in the human brain . NeuroImage 29 , 1359 – 1367 ( 2006 ). OpenUrl CrossRef PubMed Web of Science 54. van den Heuvel , M. P. , Mandl , R. C. W. , Kahn , R. S. & Hulshoff Pol , H. E. Functionally linked resting-state networks reflect the underlying structural connectivity architecture of the human brain . Hum. Brain Mapp . 30 , 3127 – 3141 ( 2009 ). OpenUrl CrossRef PubMed Web of Science 55. ↵ Wang , Z. et al. The relationship of anatomical and functional connectivity to resting-state connectivity in primate somatosensory cortex . Neuron 78 , 1116 – 1126 ( 2013 ). OpenUrl CrossRef PubMed Web of Science 56. ↵ Honey , C. J. et al. Predicting human resting-state functional connectivity from structural connectivity . Proc. Natl. Acad. Sci. U. S. A . 106 , 2035 – 2040 ( 2009 ). OpenUrl Abstract / FREE Full Text 57. ↵ Wong , A. M. F. , Lueder , G. T. , Burkhalter , A. & Tychsen , L . Anomalous retinal correspondence: Neuroanatomic mechanism in strabismic monkeys and clinical findings in strabismic children . J. Am. Assoc. Pediatr. Ophthalmol. Strabismus 4 , 168 – 174 ( 2000 ). OpenUrl 58. Reynaud , A. , Masson , G. S. & Chavane , F . Dynamics of local input normalization result from balanced short- and long-range intracortical interactions in area V1 . J. Neurosci . 32 , 12558 – 12569 ( 2012 ). OpenUrl Abstract / FREE Full Text 59. Roland , J. L. et al. On the role of the corpus callosum in interhemispheric functional connectivity in humans . Proc. Natl. Acad. Sci . 114 , 13278 – 13283 ( 2017 ). OpenUrl Abstract / FREE Full Text 60. ↵ Singer , W . Recurrent dynamics in the cerebral cortex: Integration of sensory evidence with stored knowledge . Proc. Natl. Acad. Sci . 118 , e2101043118 ( 2021 ). OpenUrl Abstract / FREE Full Text 61. ↵ Mateo , C. , Knutsen , P. M. , Tsai , P. S. , Shih , A. Y. & Kleinfeld , D . Entrainment of arteriole vasomotor fluctuations by neural activity is a basis of Blood-Oxygenation-Level-Dependent ‘resting-state’ connectivity . Neuron 96 , 936 – 948 .e3 ( 2017 ). OpenUrl CrossRef PubMed 62. Nir , Y. et al. Coupling between neuronal firing rate, gamma LFP, and BOLD fMRI is related to interneuronal correlations . Curr. Biol. CB 17 , 1275 – 1285 ( 2007 ). OpenUrl CrossRef PubMed 63. Scheeringa , R. et al. Neuronal dynamics underlying high- and low-frequency EEG oscillations contribute independently to the human BOLD signal . Neuron 69 , 572 – 583 ( 2011 ). OpenUrl CrossRef PubMed Web of Science 64. Schölvinck , M. L. , Maier , A. , Ye , F. Q. , Duyn , J. H. & Leopold , D. A . Neural basis of global resting-state fMRI activity . Proc. Natl. Acad. Sci. U. S. A . 107 , 10238 – 10243 ( 2010 ). OpenUrl Abstract / FREE Full Text 65. Shmuel , A. & Leopold , D. A . Neuronal correlates of spontaneous fluctuations in fMRI signals in monkey visual cortex: Implications for functional connectivity at rest . Hum. Brain Mapp . 29 , 751 – 761 ( 2008 ). OpenUrl CrossRef PubMed Web of Science 66. ↵ Nir , Y. et al. Interhemispheric correlations of slow spontaneous neuronal fluctuations revealed in human sensory cortex . Nat. Neurosci . 11 , 1100 – 1108 ( 2008 ). OpenUrl CrossRef PubMed Web of Science 67. ↵ Gilbert , C. D. & Wiesel , T. N . Columnar specificity of intrinsic horizontal and corticocortical connections in cat visual cortex . J. Neurosci . 9 , 2432 – 2442 ( 1989 ). OpenUrl Abstract / FREE Full Text 68. ↵ Levine , M. W. & McAnany , J. J . The relative capabilities of the upper and lower visual hemifields . Vision Res . 45 , 2820 – 2830 ( 2005 ). OpenUrl CrossRef PubMed Web of Science 69. ↵ Skrandies , W . The upper and lower visual field of man: Electrophysiological and functional differences. in Progress in Sensory Physiology (eds. Autrum , H. et al. ) 1 – 93 ( Springer, Berlin, Heidelberg , 1987 ). doi: 10.1007/978-3-642-71060-5_1 . OpenUrl CrossRef 70. ↵ Thomas , N. A. & Elias , L. J . Upper and lower visual field differences in perceptual asymmetries . Brain Res . 1387 , 108 – 115 ( 2011 ). OpenUrl CrossRef PubMed Web of Science 72. ↵ 71. Schmidt , M. E. , et al. Impact of amblyopia on mesoscale functional connectivity within the human primary visual cortex . in Investigative Ophthalmology & Visual Science vol. 66 2201 ( 2025 ). 72. ↵ Livingstone , M. S. & Hubel , D. H . Specificity of intrinsic connections in primate primary visual cortex . J. Neurosci . 4 , 2830 – 2835 ( 1984 ). OpenUrl Abstract / FREE Full Text 73. ↵ Sincich , L. C. & Blasdel , G. G . Oriented axon projections in primary visual cortex of the monkey . J. Neurosci . 21 , 4416 – 4426 ( 2001 ). OpenUrl Abstract / FREE Full Text 74. ↵ Olman , C. , Van de Moortele , P. & Ugurbil , K. Point spread function for gradient echo and spin echo BOLD fMRI at 7 Tesla . in Proceedings of the International Society for Magnetic Resonance in Medicine 1066 ( 2004 ). 75. Shmuel , A. , Yacoub , E. , Chaimow , D. , Logothetis , N. K. & Ugurbil , K . Spatio-temporal point-spread function of fMRI signal in human gray matter at 7 Tesla . Neuroimage 35 , 539 – 552 ( 2007 ). OpenUrl CrossRef PubMed Web of Science 76. ↵ Engel , S. A. , Glover , G. H. & Wandell , B. A . Retinotopic organization in human visual cortex and the spatial precision of functional MRI . Cereb Cortex 7 , 181 – 92 ( 1997 ). OpenUrl CrossRef PubMed Web of Science 77. ↵ Alvarez , I. , Parker , A. J. & Bridge , H . Normative cerebral cortical thickness for human visual areas . Neuroimage 201 , 116057 ( 2019 ). 78. ↵ Markuerkiaga , I. , Marques , J. P. , Gallagher , T. E. & Norris , D. G . Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function . J. Neurosci. Methods 353 , 109095 ( 2021 ). 79. ↵ Markuerkiaga , I. , Barth , M. & Norris , D. G . A cortical vascular model for examining the specificity of the laminar BOLD signal . NeuroImage 132 , 491 – 498 ( 2016 ). OpenUrl CrossRef PubMed 80. ↵ De Martino , F. et al. Cortical depth dependent functional responses in humans at 7T: improved specificity with 3D GRASE . PLoS One 8 , e60514 ( 2013 ). OpenUrl CrossRef PubMed 81. Koopmans , P. J. , Barth , M. & Norris , D. G . Layer-specific BOLD activation in human V1 . Hum. Brain Mapp . 31 , 1297 – 1304 ( 2010 ). OpenUrl CrossRef PubMed Web of Science 82. ↵ Polimeni , J. R. , Fischl , B. , Greve , D. N. & Wald , L. L . Laminar analysis of 7T BOLD using an imposed spatial activation pattern in human V1 . Neuroimage 52 , 1334 – 46 ( 2010 ). OpenUrl CrossRef PubMed Web of Science 83. ↵ Uludağ , K. & Blinder , P . Linking brain vascular physiology to hemodynamic response in ultra-high field MRI . NeuroImage 168 , 279 – 295 ( 2018 ). OpenUrl PubMed 84. ↵ Adams , D. L. , Piserchia , V. , Economides , J. R. & Horton , J. C . Vascular supply of the cerebral cortex is specialized for cell layers but not columns . Cereb. Cortex 25 , 3673 – 3681 ( 2015 ). OpenUrl CrossRef PubMed 85. ↵ Sereno , M. I. et al. Borders of multiple visual areas in humans revealed by functional magnetic resonance imaging . Science 268 , 889 – 93 ( 1995 ). OpenUrl Abstract / FREE Full Text 86. ↵ Tootell , R. B. et al. Functional analysis of V3A and related areas in human visual cortex . J Neurosci 17 , 7060 – 78 ( 1997 ). OpenUrl Abstract / FREE Full Text 87. ↵ Polimeni , J. R. et al. Reducing sensitivity losses due to respiration and motion in accelerated echo planar imaging by reordering the autocalibration data acquisition . Magn Reson Med ( 2015 ) doi: 10.1002/mrm.25628 . OpenUrl CrossRef PubMed 88. ↵ Mareyam , A. , Kirsch , J. E. , Chang , Y. , Madan , G. & Wald , L.L. . A 64-Channel 7T array coil for accelerated brain MRI. in Proceedings of the International Society for Magnetic Resonance in Medicine 0764 ( 2020 ). 89. ↵ Hoge , W. S. & Polimeni , J. R . ,article>Dual-polarity GRAPPA for simultaneous reconstruction and ghost correction of echo planar imaging data . Magn. Reson. Med . 76 , 32 – 44 ( 2016 ). OpenUrl CrossRef PubMed 90. ↵ Webb , A. G. & Collins , C. M . Parallel transmit and receive technology in high-field magnetic resonance neuroimaging . Int. J. Imaging Syst. Technol . 20 , 2 – 13 ( 2010 ). OpenUrl CrossRef 91. ↵ Wrede , K. H. et al. Caudal image contrast inversion in MPRAGE at 7 Tesla: problem and solution . Acad. Radiol . 19 , 172 – 178 ( 2012 ). OpenUrl CrossRef PubMed 93. ↵ 92. Fischl , B. FreeSurfer . Neuroimage 62 , 774 – 81 ( 2012 ). OpenUrl CrossRef PubMed Web of Science 93. ↵ Dale , A. M. , Fischl , B. & Sereno , M. I. Cortical surface-based analysis. I. Segmentation and surface reconstruction . Neuroimage 9 , 179 – 94 ( 1999 ). OpenUrl CrossRef PubMed Web of Science 94. ↵ Wang , J. , Nasr , S. , Roe , A. W. & Polimeni , J. R . Critical factors in achieving fine-scale functional MRI: Removing sources of inadvertent spatial smoothing . Hum. Brain Mapp . 43 , 3311 – 3331 ( 2022 ). OpenUrl CrossRef PubMed 95. ↵ Greve , D. N. & Fischl , B . Accurate and robust brain image alignment using boundary- based registration . Neuroimage 48 , 63 – 72 ( 2009 ). OpenUrl CrossRef PubMed Web of Science 96. ↵ Moeller , S. et al. NOise reduction with DIstribution Corrected (NORDIC) PCA in dMRI with complex-valued parameter-free locally low-rank processing . Neuroimage 226 , 117539 ( 2021 ). 97. ↵ Dowdle , L. T. et al. Evaluating increases in sensitivity from NORDIC for diverse fMRI acquisition strategies . NeuroImage 270 , 119949 ( 2023 ). 98. ↵ Kaptan , M. et al. Reliability of resting-state functional connectivity in the human spinal cord: Assessing the impact of distinct noise sources . NeuroImage 275 , 120152 ( 2023 ). 99. ↵ Chen , J. E. , Jahanian , H. & Glover , G. H . Nuisance regression of high-frequency functional magnetic resonance imaging data: denoising can be noisy . Brain Connect . 7 , 13 – 24 ( 2017 ). OpenUrl CrossRef PubMed 100. ↵ Blazejewska , A. I. , Fischl , B. , Wald , L. L. & Polimeni , J. R . Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data . Neuroimage 189 , 601 – 614 ( 2019 ). OpenUrl CrossRef PubMed 101. ↵ Tootell , R. , Hamilton , S. , Silverman , M. & Switkes , E . Functional anatomy of macaque striate cortex . I. Ocular dominance, binocular interactions, and baseline conditions. J. Neurosci . 8 , 1500 – 1530 ( 1988 ). OpenUrl PubMed 102. ↵ Renvall , V. , Witzel , T. , Wald , L. L. & Polimeni , J. R . Automatic cortical surface reconstruction of high-resolution T1 echo planar imaging data . NeuroImage 134 , 338 – 354 ( 2016 ). OpenUrl PubMed 103. Aquino , K. M. , Fulcher , B. D. , Parkes , L. , Sabaroedin , K. & Fornito , A . Identifying and removing widespread signal deflections from fMRI data: Rethinking the global signal regression problem . NeuroImage 212 , 116614 ( 2020 ). 104. Murphy , K. & Fox , M. D . Towards a consensus regarding global signal regression for resting state functional connectivity MRI . NeuroImage 154 , ( 2017 ). 105. ↵ Bi , H. et al. Neuronal responses in visual area V2 (V2) of macaque monkeys withstrabismic amblyopia . Cereb. Cortex 21 , 2033 – 2045 ( 2011 ). OpenUrl CrossRef PubMed Web of Science 106. Hubel , D. H. & Wiesel , T. N . Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex . J. Physiol . 160 , 106 – 154 ( 1962 ). OpenUrl CrossRef PubMed Web of Science 107. ↵ Kiorpes , L. , Kiper , D. C. , O’Keefe , L. P. , Cavanaugh , J. R. & Movshon , J. A . Neuronal correlates of amblyopia in the visual cortex of macaque monkeys with experimental strabismus and anisometropia . J. Neurosci . 18 , 6411 – 6424 ( 1998 ). OpenUrl Abstract / FREE Full Text 108. ↵ Aganj , I. , Yeo , B. T. T. , Sabuncu , M. R. & Fischl , B . On removing interpolation and resampling artifacts in rigid image registration . IEEE Trans. Image Process . 22 , 816 – 827 ( 2012 ). OpenUrl PubMed 109. ↵ Halton , J. H . On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals . Numer. Math . 2 , 84 – 90 ( 1960 ). OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted October 01, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI Marianna E. Schmidt , Iman Aganj , Jason Stockmann , Berkin Bilgic , Yulin Chang , W. Scott Hoge , Evgeniya Kirilina , Nikolaus Weiskopf , Shahin Nasr bioRxiv 2025.03.27.645795; doi: https://doi.org/10.1101/2025.03.27.645795 Share This Article: Copy Citation Tools Unraveling the mesoscale resting-state functional connectivity of ocular dominance columns in humans using high-resolution functional MRI Marianna E. Schmidt , Iman Aganj , Jason Stockmann , Berkin Bilgic , Yulin Chang , W. Scott Hoge , Evgeniya Kirilina , Nikolaus Weiskopf , Shahin Nasr bioRxiv 2025.03.27.645795; doi: https://doi.org/10.1101/2025.03.27.645795 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Neuroscience Subject Areas All Articles Animal Behavior and Cognition (7629) Biochemistry (17660) Bioengineering (13881) Bioinformatics (41911) Biophysics (21436) Cancer Biology (18578) Cell Biology (25482) Clinical Trials (138) Developmental Biology (13371) Ecology (19887) Epidemiology (2067) Evolutionary Biology (24302) Genetics (15599) Genomics (22483) Immunology (17728) Microbiology (40364) Molecular Biology (17163) Neuroscience (88537) Paleontology (666) Pathology (2830) Pharmacology and Toxicology (4821) Physiology (7637) Plant Biology (15129) Scientific Communication and Education (2045) Synthetic Biology (4290) Systems Biology (9817) Zoology (2269)