P Barttfeld

• S Abboud

H Lagercrantz

• U Adén

N Padilla

A D Edwards

L Cohen

M Sigman

• S Dehaene

G Dehaene-Lambertz

A lateral-to-mesial organization of human ventral visual cortex at birth

Keywords: Neonates, Brain, Functional connectivity, Language P. Barttfeld and S

In human adults, ventral extra-striate visual cortex contains a mosaic of functionally specialized areas, some responding preferentially to natural visual categories such as faces (fusiform face area) or places (parahippocampal place area) and others to cultural inventions such as written words and numbers (visual word form and number form areas). It has been hypothesized that this mosaic arises from innate biases in cortico-cortical connectivity. We tested this hypothesis by examining functional resting-state correlation at birth using fMRI data from full-term human newborns. The results revealed that ventral visual regions are functionally connected with their contra-lateral homologous regions and also exhibit distinct patterns of long-distance functional correlation with anterior associative regions. A mesial-to-lateral organization was observed, with the signal of the more lateral regions, including the sites of visual word and number form areas, exhibiting higher correlations with voxels of the prefrontal, inferior parietal and temporal cortices, including language areas. Finally, we observed hemispheric asymmetries in the functional correlation of key areas of the language network that may influence later adult hemispheric lateralization. We suggest that long-distance circuits present at birth constrain the subsequent functional differentiation of the ventral visual cortex.

Introduction

The human brain displays a remarkable consistency and reproducibility in its functional organization across individuals. This was first demonstrated in neuropsychological studies that mapped specific cognitive deficits to reproducible cortical lesion sites and more recently by brain-imaging research, which has revealed surprisingly similar functional subdivisions of the cortical surface across numerous individuals [START_REF] Glasser | A multi-modal parcellation of human cerebral cortex[END_REF][START_REF] Huth | Natural speech reveals the semantic maps that tile human cerebral cortex[END_REF]. Not only is this reproducible topological specialization evident in domains of evolutionary relevance, such as for the fusiform face area (FFA) and the parahippocampal place area (PPA) (Downing et al. 2006), but it also extends to cultural acquisitions such as reading. For example, the same region of the left fusiform gyrus, the visual form area (vWFA), is consistently activated whenever literate participants recognize written words, independently of case [START_REF] Dehaene | The unique role of the visual word form area in reading[END_REF], modality [START_REF] Striem-Amit | Reading with sounds: sensory substitution selectively activates the visual word form area in the blind[END_REF] or the age of reading acquisition [START_REF] Dehaene | How learning to read changes the cortical networks for vision and language[END_REF].

Such functional reproducibility suggests that cognitive processes are dependent on specific neural architectures [START_REF] Hannagan | Origins of the specialization for letters and numbers in ventral occipitotemporal cortex[END_REF]. Indeed, functional responses have been consistently related to specific micro-architectures (e.g. word-and face-specific activations are observed in the regions FG2 and FG4, whereas places are activating FG1 and FG3, [START_REF] Weiner | The cytoarchitecture of domain-specific regions in human high-level visual cortex[END_REF][START_REF] Caspers | Receptor architecture of visual areas in the face and word-form recognition region of the posterior fusiform gyrus[END_REF], although the exact relation between a given cortical layer organization and a computational process is still poorly specified. Anatomical connectivity, by orchestrating local and distant synchrony, might also drive the functional role of an area. This "biased connectivity" hypothesis was recently supported by two studies showing that the precise locations of the FFA [START_REF] Saygin | Anatomical connectivity patterns predict face selectivity in the fusiform gyrus[END_REF] and VWFA [START_REF] Saygin | Connectivity precedes function in the development of the visual word form area[END_REF]) can be predicted from their connections. The functional localization of the VWFA was predicted in 8-year-old readers by the anatomical fiber tracking performed at the age of 5 before they learned to read [START_REF] Saygin | Connectivity precedes function in the development of the visual word form area[END_REF].

Because all main fiber tracts are already present in fullterm human neonates [START_REF] Dubois | The early development of brain white matter: a review of imaging studies in fetuses, newborns and infants[END_REF][START_REF] Dubois | MRI and M/EEG studies of the white matter development in human fetuses and infants: review and opinion[END_REF]Kennedy et al. 2007;[START_REF] Kostović | The development of the subplate and thalamocortical connections in the human foetal brain[END_REF][START_REF] Marín-Padilla | Human motor cortex: development and cytoarchitecture[END_REF][START_REF] Takahashi | Linking oligodendrocyte and myelin dysfunction to neurocircuitry abnormalities in schizophrenia[END_REF], in the present paper, we asked whether the initial connectome, present at birth, may explain at least some aspects of the subsequent functional specialization. We focused here on the mesial-to-lateral organization of ventral visual cortex. The ventral visual pathway comprises a succession of functional areas extending from mesial to lateral regions including the mesial place area (PPA), the face fusiform area (FFA), the visual word form area (VWFA), the number form area (NFA) and a more lateral area responding to tools and objects. Furthermore, distinct face and place responses are already present in 6-month-old infants [START_REF] Deen | Organization of high-level visual cortex in human infants[END_REF]. Here, we examined whether a particular connectivity pattern of these regions is already in place at birth.

To this aim, we relied on resting-state functional correlation, i.e. the synchronization of the BOLD responses across distant regions at rest due to the spontaneous co-activity of different nodes of the same network [START_REF] Smith | Functional connectomics from resting-state fMRI[END_REF][START_REF] Glasser | A multi-modal parcellation of human cerebral cortex[END_REF][START_REF] Raichle | The restless brain: how intrinsic activity organizes brain function[END_REF]. Using this approach in full-term and even pre-term neonates, several studies have reported the presence of the main resting-state networks described in adults [START_REF] Fransson | Resting-state networks in the infant brain[END_REF][START_REF] Fransson | Spontaneous brain activity in the newborn brain during natural sleep-an fMRI study in infants born at full term[END_REF][START_REF] Fransson | The functional architecture of the infant brain as revealed by resting-state fMRI[END_REF][START_REF] Doria | Emergence of resting state networks in the preterm human brain[END_REF][START_REF] Smyser | Longitudinal analysis of neural network development in preterm infants[END_REF][START_REF] Smyser | Functional connectivity MRI in infants: exploration of the functional organization of the developing brain[END_REF] as also more specific connectivity, such as thalamocortical connectivity [START_REF] Toulmin | Specialization and integration of functional thalamocortical connectivity in the human infant[END_REF] or from multisensory areas [START_REF] Sours | Cortical multisensory connectivity is present near birth in humans[END_REF]. Our goal was to go beyond the whole-brain description and determine whether long-distance connections specifically arising from ventral visual areas already target distinct regions. We aimed to evaluate two contrasting theories for the origins of adult functional specialization. The first of which, the "biased connectivity" hypothesis [START_REF] Hannagan | Origins of the specialization for letters and numbers in ventral occipitotemporal cortex[END_REF][START_REF] Saygin | Anatomical connectivity patterns predict face selectivity in the fusiform gyrus[END_REF] as explained above, predicts that mesial-to-lateral differences in functional correlation are already present at birth, in approximate correspondence to the functional mosaic later found in older children and adults. We then contrasted it to the latter hypothesis according to which brain connectivity is initially diffuse, overabundant and largely isotropic and only becomes refined at a later stage during infancy and childhood through exposure to a given environment [START_REF] Johnson | Interactive specialization: a domain-general framework for human functional brain development?[END_REF][START_REF] Elman | Rethinking innateness: a connectionist perspective on development[END_REF].

The "biased connectivity" hypothesis makes specific predictions for which we searched in our data. Given the reproducibility of the location of the visual word and number form areas, we wondered whether these regions were already more connected with the temporo-frontal language network and the inferior parietal "number sense" regions than other mesial regions, which would explain their propensity to acquire the role of processing arbitrary shapes to represent phonemes or numbers. Similarly, we wondered whether hemispheric biases in connectivity at birth already predated the later lateralized responses to faces and words in the right and left fusiform gyrus [START_REF] Weiner | The cytoarchitecture of domain-specific regions in human high-level visual cortex[END_REF]. We thus analyzed whether key regions of the linguistic network (i.e. pars triangularis of the inferior frontal gyrus, superior and middle temporal gyrus) are differently related to the left and right ventral temporal regions.

Methods

Subjects

Twenty-four healthy full-term neonates scanned for research purposes were included in this study. The fMRI data were provided by the Karolinska University Hospital, Stockholm (14 neonates) and by the Centre for the Developing Brain, King's College London (10 neonates) and have already been published [START_REF] Fransson | Spontaneous brain activity in the newborn brain during natural sleep-an fMRI study in infants born at full term[END_REF][START_REF] Doria | Emergence of resting state networks in the preterm human brain[END_REF]. All MR examinations were carried out according to the ethical guidelines and declarations of the Declaration of Helsinki (1975). Details can be found in the seminal papers. Six subjects were subsequently rejected due to excessive movement artefacts (see "Preprocessing" section). We thus analyzed 16 girls and 9 boys aged 41.1 weeks (std: 1.08) in average at scan.

MR image acquisition

For the Swedish cohort, infants were scanned during natural sleep. Functional MRI scans were acquired on a Philips Intera 1.5 T scanner equipped with a 6-channel receiveonly head coil. Echo planar imaging [time repetition (TR)/ time echo (TE)/flip angle = 2000 ms/50 ms/80°, voxel size = 2.8 × 2.8 × 4.5 mm, 20 slices] of the infant brain was performed for 10 min (300 EPI volumes). For the London cohort, infants were sedated with chloral hydrate (25-50 mg/ kg) and images were acquired at 3 T (Achieva, Philips, Best, Netherlands) with an 8-channel phased array head coil. Echo planar imaging (TR/TE/flip angle = 1500 ms/45 ms/90º, 2.5 × 2.5 × 3.25 mm, 22 slices) was performed for 6.4 min. In both groups, T1-weigthed images were acquired in each neonate. Although it is unconventional to merge data obtained from two different MRI scanners for combined group-level inference, we were confronted with the scarcity of the publically available data at this age and the small size of the groups. Note that our first-level comparisons were computed within-subject searching for differences in functional connectivity between voxels. By merging the two groups, we probably increased the variance between subjects but we benefited from more degrees of freedom. To capture the group differences, a group-factor was declared in our statistical analyses.

Preprocessing

FMRI images were all pre-processed using SPM8 (Wellcome Institute of Imaging Neuroscience, London, UK) implemented in Matlab (Mathworks, Natick, MA, USA). Pre-processing involved within-subject realignment to account for head motion, spatial normalization to a template and smoothing: all images were realigned to the mean image of the resting state session using linear transformations. The mean EPI image was spatially normalized to a 7-week-old infant template image for whom a sulci-based parcellation is available [START_REF] Kabdebon | Anatomical correlations of the international 10-20 sensor placement system in infants[END_REF], thus providing anatomically defined region of interests to study functional connectivity. Brain growth and sulcation massively develop during the last trimester of gestation but decelerate after term-age. Although tertiary sulcation begins around full-term birth and continues to develop until adulthood, all main structures important for normalization are visible (see [START_REF] Kabdebon | Anatomical correlations of the international 10-20 sensor placement system in infants[END_REF], for quantitative measures of inter-age variability of the location of the main sulci during the first semester of life). The transformation parameters obtained from this spatial normalization were subsequently applied to every functional image acquired during the resting state scanning session with a target voxel size of 2 × 2 × 2 mm 3 . Then, spatially normalized images were smoothed with an isotropic 5-mm 3 full-width half-maximum Gaussian kernel. Next, we regressed out the effect of the following nuisance variables from the voxel time series using the REST toolbox for Matlab (http://restf mri.net): the six movement parameters from the head motion correction stage along with the global signal of the entire brain and four nuisance signals obtained from a 3-mm sphere located in the ventricles and white matter. Finally, all functional MRI time series were band-pass filtered (0.01-0.08 Hz).

Spontaneous movement during sleep in infants often consists of sudden hyperextensions along the body axis which produced spiking artefacts in the time series and might have, therefore, disrupted anterior-posterior correlations as a result of these pitch movements. Because subject motion is particularly deleterious on long-distance correlations [START_REF] Power | Recent progress and outstanding issues in motion correction in resting state fMRI[END_REF], we removed the artifacted time points through a semi-automatic method. First, we deleted brain volumes detected by the ART toolbox for SPM8 as contaminated by high levels of movement. Second, we visually inspected whether all traces of movement were removed from the voxel time series. When needed, more brain volumes were discarded until no trace of movement was visually detected in the time series. Concerns about such scrubbing have been raised in the literature as it decreases the degrees of freedom and destroys the temporal time-series when too many successive brain volumes are removed (for a discussion see [START_REF] Power | Recent progress and outstanding issues in motion correction in resting state fMRI[END_REF]. Furthermore, given the non-stationarity of the functional correlations, some stages during which subjects might naturally move may be under-represented. These concerns are particularly critical when different populations with different levels of agitation are compared, such as infants with adults. However, this concern does not apply to our study, since we are comparing the connectivity between different voxels that have been similarly processed within subjects. In the end, six subjects whose remaining time-series were shorter than 130 scans were discarded to avoid noisy correlation estimations. All included subjects had at least 138 brain volumes (mean number of volumes: 218, SD: 47, min: 132, max: 300). We also estimated the temporal signal-to-noise ratio (tSNR) for every voxel of every infant, as the inverse of the coefficient of variation in time (tSNR = mean/stdev of the times series). There was no difference between the two groups [t(22) = 1.25, p = 0.22].

Clustering analysis

To parcel the ventral areas in clusters of voxels with similar patterns of functional correlation, we applied the following data-driven clustering method. First, we performed a functional correlation analysis (Brain Voyager, Neuroelf v1.0, http://neuro elf.net/) seeded from each voxel belonging to a mask of the ventral stream of the visual cortex, i.e., ventrally to the Calcarine sulcus expanding towards the anterior parahippocampal gyrus and limited to y = -14 in the template (1484 voxels, corresponding to the union of the clusters presented in Fig. 1). In other words, in each individual neonate, time-series were extracted from each functional voxel within the mask z-normalized and used as general linear model (GLM) predictors. For each voxel in the mask, random-effects inference was then performed on the group level. Second, each resulting group-level map was masked by a gray-matter mask, thresholded (voxelwise p < 0.01 uncorrected), binarized and reshaped into a numerical vector. We then represented the functional correlation from all seeds as a matrix of 1484 columns (the number of seeds) and 14 808 rows (the number of gray-matter brain voxels). Third, we used an agglomerative hierarchical clustering algorithm to cluster the seed voxels with a similar spatial pattern of functional correlation (i.e. a similar matrix columns). A hierarchical cluster tree was constructed using the Ward linkage measure on Euclidean distances (minimum variance algorithm merging at each step the pair of clusters that leads to the minimal increase in total within-cluster variance). The final number of clusters was optimized using the silhouette technique which quantifies the ratio between the largest inter-cluster distance over the smallest distance to a point outside the cluster. An assessment of all solutions in the range of 2-15 clusters showed a local maximum for 4 clusters. Fourth, in order to study the brain networks correlated with each of these four clusters, we performed a group-level analysis similar to the one described above, but using the average time-series from all the voxels of each cluster computed in each individual neonate. We also computed partial correlation coefficients for each cluster in order to regress out the time-series of all other clusters and thus obtain its unique correlation pattern (Fig. 1). Significance thresholds were set to q = 0.01 (FDR-corrected for multiple comparisons) when using partial correlation and q = 0.001 (FDR-corrected for multiple comparisons) when using Pearson's correlation.

Seed analysis

Clustering analyses provide a coarse grouping of voxels with similar connectivity tendencies but blind to anatomical boundaries. To reveal a finer dissection of the ventral connectivity, we used the anatomical parcellation proposed by [START_REF] Kabdebon | Anatomical correlations of the international 10-20 sensor placement system in infants[END_REF] on a two-month-old infant template to define defined four ROIs, based on the gray matter of 4 sulci bordering the lingual, fusiform (internal part), fusiform (external part) and the inferior temporal gyrus. Those ROIs were used as seeds for the correlation analysis (Fig. 2a). For each subject and each ROI, we defined the time-series T as the average time-series across all the voxels within the ROI. Then, we calculated the Pearson correlation between T and the time-series of all the voxels in the brain, V. Through this procedure, we obtained four correlation brain volumes or maps per subject s and ROI r (C s,r ), which were Fishertransformed (Matlab function atanh) to produce Z s,r . To test for significant correlations across subjects between T and all voxels in the brain, we performed a one-way t test on each ROI r, setting the significance value at p < 0.01, FDR-corrected voxelwise. We also performed a within-subjects ANOVA between maps so as to compare pairs of Z s,r for each individual subject, setting the significance value at p < 0.01 for visualization purposes (a more strict threshold of p < 0.001 did not qualitatively change the results).

Regression analysis

Given that the boundaries between ROIs in the previous analysis were based on sulci that might not follow cytoarchitectonic boundaries [START_REF] Fischl | Cortical folding patterns and predicting cytoarchitecture[END_REF], we also performed an analysis free of this constraint through voxel-wise regression to explore the mesial-to-lateral gradient in a ROI-independent manner. We constructed a mask comprising both left and right ventral regions (Fig. 4a). We reduced this mask to a x,y grid along the left-right and antero-posterior axes, merging together all voxels along the z axis (dorsal-ventral axis), as the mosaic of activations to visual categories that we wanted to explore is mainly evident along these two axes [START_REF] Weiner | The cytoarchitecture of domain-specific regions in human high-level visual cortex[END_REF]. We computed the correlation map for each x,y pair and the time series of all voxels in the brain in each subject and obtained a total of 276 (for left hemisphere) and 258 (for right hemisphere) correlation maps per subject. To identify regions whose correlation with voxels in the ROI would depend on their x and y coordinates, we performed a first-level regression analysis (using SPM) for each subject. We then submitted the x-related contrast to a second-level analysis, in order to quantify the influence of the x coordinate on functional correlation. We modeled the cohort origin (i.e. England or Sweden) explicitly as blocks in the analysis. The resulting second-level beta map and its associated t-map revealed the relation between the x axis and functional correlation: the higher the beta and the t value, the more connected each voxel was with increasing values of x (i.e. from medial to lateral localization). Conversely, the smaller the effect the more connected each voxel was with decreasing values of x. We performed this procedure for the left and the right hemispheres, inverting the x values when necessary to keep the convention that increasing x values correspond to lateral locations. We focused here on the medial-lateral axis (x coordinate) along which specific activation for the different visual categories are ordered whereas the anterior-posterior axis (y coordinate) is more sensitive to the hierarchical analyses within a category, more difficult to segregate in particular across subjects [START_REF] Tanaka | Inferotemporal cortex and object vision[END_REF]. We wondered whether the differences we observed between mesial and lateral regions were related to their level of maturation which might have created different delays in neurovascular coupling. As a control procedure aimed at accounting for this possibility, we repeated the analysis including the average T2 value of the template brain at the same locations as a regressor and thus a proxy for maturation index [START_REF] Leroy | Atlas-free surface reconstruction of the cortical grey-white interface in infants[END_REF]).

Analysis of hemispheric differences

Because the VWFA and the FFA are, respectively, left and right lateralized in adults, we studied whether differences in functional correlation were observed early in development, between the ventral areas and key regions of the linguistic network. Therefore, we recovered the following anatomically pre-defined ROIs from the infant atlas provided by [START_REF] Kabdebon | Anatomical correlations of the international 10-20 sensor placement system in infants[END_REF]: the pars triangularis of the inferior frontal gyrus (region number 4 in the atlas), the middle temporal gyrus (region 27) and the superior temporal gyrus (region 28). We truncated regions 27 and 28 along the anterior-posterior axis at the level of the deepest angle of the planum temporale so as to keep the ROIs' most anterior parts. This was done to avoid any overlap or vicinity of the seeds with the visual stream that could have biased the analysis with spurious local correlations.

For every subject and every pair of left-right regions, (e.g. left and right pars triangularis), we performed two seed analyses, the first between the left ROI time-series and the rest of the brain and the second using the right ROI timeseries. We thus obtained two functional correlation maps per subject and ROI pair. To test for hemispheric laterality biases, we performed a series of t test (for each value of x coordinate) on the left-right difference collapsing all y coordinate values between values of y coordinate of 9 and 13. The significance threshold was set at p < 0.05, Bonferroni corrected for multiple comparisons (42 comparisons, corresponding to the 42 voxels along the x coordinate).

Results

Parcellation of visual ventral regions

We first performed a data-driven clustering of the voxels in the bilateral ventral regions according to their functional correlation pattern. Voxels were clustered together or separately as a function of the similarity of the pattern of correlation with the entire brain. This clustering procedure was combined with a partial correlation analysis in order to display the unique functional fingerprint of each cluster (i.e. not shared with other clusters). We identified four main clusters: C1, or Medial Occipital, posterior, an exclusively posterior occipital group connected with the frontal pole, the right inferior temporal gyrus and paracentral regions consisting of the right post-central gyrus, paracentral lobule and supplementary motor area; C2, or Medial Occipital, anterior, connected with neighboring mesial areas (lingual, fusiform, parahipocampal gyri, and the cuneus) and more distant primary regions, such as the central sulcus, Heschl's gyrus, planum temporale, bilateral posterior insula and the center of the right STAP region (the STAP is a region of the superior temporal sulcus under Heschl's gyrus, deeper in the right than in the left hemisphere in almost all humans, Leroy et al. 2015); C3, or Anterior Lateral, a clusterwhich had the most extended unique functional correlation, connected to the inferior and middle temporal lobes, the inferior parietal regions comprising the inter-parietal sulcus, and the precentral gyri; C4, or Posterior Lateral, a cluster connected to the dorsal visual steam, cuneus and left dorso-lateral prefrontal region (Fig. 1 for functional correlation using Pearson's correlation).

Functional correlation based on anatomical ROIs

The parcels defined above correspond to an optimal statistical solution given our data but may differ from the anatomical reality. We thus sought to confirm the above conclusions using four anatomical ROIs progressively aligned along the medial-external axis of the left ventral temporal cortex in order to determine how they connected to the rest of the brain. As expected from the previous analysis, the observed functional correlation was not limited to neighboring ipsilateral areas or to homologous regions in the right hemisphere (Fig. 2). Rather, the connections crossed hemispheres to join distant associative areas in the right hemisphere. We observed a different correlation pattern depending on ROI position (Fig. 3a,b): the most mesial ROI (ROI1, or Lingual gyrus) functionally related to the superior temporal region as well as central motor and sensory regions, whereas the most lateral ROI (ROI4, inferior temporal gyrus) was principally correlated with the whole temporal lobe, the inferior parietal region and the prefrontal region. ROIs 1 and 4 thus shared similarities with the clusters C2 and C3 of the previous analysis.

Pairwise comparisons between the patterns of functional correlation of the different ROIs (Fig. 2) revealed that ROI1 (lingual gyrus) was most strongly correlated with the lower part of the central sulcus and Heschl's gyrus bilaterally. ROI2, the medial part of the fusiform gyrus presented only local correlations which extended more anteriorly in the left fusiform and inferior temporal gyri relative to ROI1. A sharp change occurred with ROI3, the lateral part of the fusiform gyrus, which displayed correlations with the ventral regions in both hemispheres and moreover with the supra-marginal and angular gyri, also in both hemispheres. This change corresponds to the mid-fusiform sulcus, separating the internal and external part of the fusiform gyrus [START_REF] Weiner | The mid-fusiform sulcus: a landmark identifying both cytoarchitectonic and functional divisions of human ventral temporal cortex[END_REF]. Finally, ROI4 (inferior temporal gyrus) displayed the most extended correlation pattern, exhibiting stronger correlations with the whole temporal lobe and prefrontal regions relatively to the other ROIs and notably to ROI3. It is noticeable that the significant changes of functional correlation patterns between the lateral part of the fusiform (ROI3) and the inferior temporal gyrus (ROI4) spared the parietal and frontal regions.

A possible confound is a similar neurovascular coupling in regions at the same maturational stage. While there is no indication that equally maturating regions will show a correlated functional activity in time, to account for this possibility we repeated the analysis including the average T2 value of the template brain at the same location as a regressor and thus as a proxy for maturation index [START_REF] Leroy | Atlas-free surface reconstruction of the cortical grey-white interface in infants[END_REF]. Because dendritic development impairs water motion in the gray matter and myelination in the white matter, T2 signal becomes progressively darker during the first year of life and can hence be used as a proxy for neural maturation [START_REF] Leroy | Atlas-free surface reconstruction of the cortical grey-white interface in infants[END_REF]. We observed virtually the same functional correlation pattern, suggesting that maps presented in Fig. 2 are not directly related to a maturational neuro-vascular gradient.

A third, voxel-wise, analysis aimed to verify that the mesial-to-lateral organization of the fusiform gyrus was not biased by our arbitrary, although anatomical, ROI definition. To this aim, we examined the functional correlation between each x,y location in a ventral region encompassing the calcarine sulci, lingual gyrus, fusiform gyrus and inferior temporal gyrus and examined the dependency of this functional correlation pattern on the x-coordinate of the voxel of origin (see "Methods"). In agreement with our previous ROI analysis, functional correlation was sensitive to the x-coordinates (Fig. 4): Heschl's gyrus and posterior insula signal showed a greater functional correlation with mesial voxels in the ventral visual cortex (thus appearing as a negative correlation in Fig. 4b), whereas associative areas of the temporal lobe (sparing the superior temporal region) in both hemispheres, the left angular and supramarginal gyri and both prefrontal regions showed a greater x coordinate (MNI)

x coordinate (MNI)

x coordinate (MNI)

x coordinate (MNI)

x coordinate (MNI) functional correlation with lateral voxels (thus appearing as a positive correlation in Fig. 4b; see Table S2).

Hemispheric asymmetries

Finally, to study putative asymmetries in the manner in which right and left ventral streams relate to the language network, we performed a within-subject comparison of the functional correlation maps arising from three pairs of symmetrical anatomical ROIs located at key language-related areas: pars triangularis in the inferior frontal gyrus (IFG), superior and middle temporal gyri (STG and MTG). show to what extent each site in the left and right ventral visual cortex was more strongly correlated to the left or right IFG, STG and MTG. Unsurprisingly, STG and MTG exhibited a stronger correlation with the ipsilateral than with the contralateral hemisphere, but this relation was stronger in the left hemisphere at a peak close to the location of the VWFA in adults (x = -36 and y = -30, MNI coordinates). After correction for multiple comparisons, this peak remained significant only for the MTG (Fig. 4 g,j). More surprisingly, the left IFG signal was more correlated than the right IFG with a right lateral ventral spot at x = 24 and y = -30 (MNI coordinates).

Discussion

We explored the pattern of functional correlations of the ventral visual cortex with the rest of the brain in sleeping full-term neonates. Both of our analyses, whether datadriven (Fig. 1) or based on anatomically defined regions (Fig. 2), revealed a clear partition of the ventral regions supporting the hypothesis of an early functional organization. Mesial ventral areas were notably correlated with other primary regions whereas the most lateral regions were more correlated with the whole temporal lobe, the inferior parietal region and the prefrontal region with a transition occurring between ROI 2 and 3 (Fig. 2). This boundary roughly corresponds to the mid-fusiform sulcus which splits the fusiform gyrus [START_REF] Weiner | The mid-fusiform sulcus: a landmark identifying both cytoarchitectonic and functional divisions of human ventral temporal cortex[END_REF]. This sulcus has been proposed as a robust cyto-architectonic landmark between FG1 and FG3 on one side and FG2 and FG4 on the other side, and functional landmark between the place area and the FFA/ VWFA in adults [START_REF] Weiner | The mid-fusiform sulcus: a landmark identifying both cytoarchitectonic and functional divisions of human ventral temporal cortex[END_REF][START_REF] Weiner | The cytoarchitecture of domain-specific regions in human high-level visual cortex[END_REF]. Finally, studying hemispheric differences in functional correlations, two results emerged: (1) when comparing intra-hemispheric correlations, we observed a stronger correlation between the middle temporal gyrus and the fusiform region in the left than in the right hemisphere, and (2) for interhemispheric correlations, we observed a stronger correlation between a right lateral ventral area and the contralateral left IFG than with the ipsilateral right IFG. This organization of the associative visual regions is all the more remarkable given that vision is very limited in utero, contrary to the other senses which are regularly stimulated through the fetus' spontaneous movements, internal and external sounds crossing the maternal abdomen and a tasty and odorant amniotic liquid reflecting the mother's meals. This early organization of the visual areas at birth, beyond primary areas, dismisses the old but persisting idea of a blank slate newborn [START_REF] Johnson | Interactive specialization: a domain-general framework for human functional brain development?[END_REF][START_REF] Elman | Rethinking innateness: a connectionist perspective on development[END_REF]. Although the human infant appears immature at birth due to his poor motor repertoire, early learning is thus supported by a pre-existing cortical architecture now extensively documented by a decade of functional brain imaging studies in infants [START_REF] Bartocci | Pain activates cortical areas in the preterm newborn brain[END_REF][START_REF] Arichi | Somatosensory cortical activation identified by functional MRI in preterm and term infants[END_REF][START_REF] Mahmoudzadeh | Syllabic discrimination in premature human infants prior to complete formation of cortical layers[END_REF][START_REF] Mahmoudzadeh | Functional maps at the onset of auditory inputs in very early preterm human neonates[END_REF][START_REF] Sours | Cortical multisensory connectivity is present near birth in humans[END_REF][START_REF] Toulmin | Specialization and integration of functional thalamocortical connectivity in the human infant[END_REF].

Our working hypothesis went beyond a mere organization at birth, and we studied whether the functional specificity described at later ages in ventral temporal areas for visual categories such as words, faces, tools and places might partly be based on a pre-existing pattern of functional correlation of these regions with the rest of the brain ("biased connectivity" hypothesis) contrasting with the hypothesis of an initially diffuse connectivity only secondarily refined through exposure to a given environment [START_REF] Johnson | Interactive specialization: a domain-general framework for human functional brain development?[END_REF][START_REF] Elman | Rethinking innateness: a connectionist perspective on development[END_REF]. Anatomically, long-range pathways, in particular those connecting frontal regions with temporal and occipital regions, can be observed in post-mortem brains from the second trimester of gestation on [START_REF] Takahashi | Linking oligodendrocyte and myelin dysfunction to neurocircuitry abnormalities in schizophrenia[END_REF]. Even if the tracts are initially not well myelinated, the structural architecture of the human brain set up during the last month of gestation offers the possibility of longdistance synchrony. Previous studies of functional correlations in full-term and preterm neonates have indeed reported a set of maps of functionally synchronized regions similar to those observed in adults [START_REF] Fransson | Resting-state networks in the infant brain[END_REF][START_REF] Fransson | Spontaneous brain activity in the newborn brain during natural sleep-an fMRI study in infants born at full term[END_REF][START_REF] Fransson | The functional architecture of the infant brain as revealed by resting-state fMRI[END_REF][START_REF] Doria | Emergence of resting state networks in the preterm human brain[END_REF][START_REF] Smyser | Longitudinal analysis of neural network development in preterm infants[END_REF][START_REF] Smyser | Functional connectivity MRI in infants: exploration of the functional organization of the developing brain[END_REF] Here, we explored whether a finer-grained pattern of functional correlation was present within ventral visual cortex, correlating with the functional specificity for visual categories observed in adults. We focused here on the medial-lateral axis along which category-selective domains are ordered rather than the posterior-anterior axis which is dependent on the hierarchy of visual analyses within a domain [START_REF] Tanaka | Inferotemporal cortex and object vision[END_REF][START_REF] Arcaro | A hierarchical, retinotopic protoorganization of the primate visual system at birth[END_REF][START_REF] Vinckier | Hierarchical coding of letter strings in the ventral stream: dissecting the inner organization of the visual word-form system[END_REF]. Although this last axis is probably more difficult to explore because its organization might depend on local computations along the visual hierarchy, Figs. 1 and 4b are indicative of another axis of organization that should be explored in further studies.

Our results indicate that the most lateral regions are more connected to distal associative regions, notably in the frontal lobe, than mesial regions are (compare cluster 3 and ROI 4 (inferior temporal gyrus) to cluster 2 and ROI 1 (lingual gyrus) in Figs. 1, 2, respectively). The privileged connectivity of such lateral regions residing at the top of the hierarchy of early visual regions may thus explain why the cultural symbols of letters and digits later develop at this location. This aspect of our results corroborate prior findings in older children, which have indicated that anatomical connectivity prior to reading acquisition (at 5 years of age) of ventral visual voxels to temporal and frontal regions was predictive of those voxels which would later exhibit a preference for written words at 8 years of age, once children had learned to read [START_REF] Saygin | Connectivity precedes function in the development of the visual word form area[END_REF]. In both cases, connectivity (anatomical or functional) precedes functionality, thus supporting the "biased connectivity" hypothesis for the origins of ventral visual specialization.

By contrast, we found that mesial areas are functionally more connected to sensorimotor areas, Heschl's gyrus and the posterior insula. This observation could be mediated by direct connections between visual primary areas and auditory and parietal cortices which have been reported in newborn kittens [START_REF] Innocenti | Development of projections from auditory to visual areas in the cat[END_REF]) and adult macaques [START_REF] Rockland | Multisensory convergence in calcarine visual areas in macaque monkey[END_REF]. Alternatively, the observed functional correlation between mesial visual regions and other primary cortices (auditory areas and motor sensory areas at the level of the mouth and tongue) could also be driven by thalamic synchrony. [START_REF] Alcauter | Development of thalamocortical connectivity during infancy and its cognitive correlations[END_REF], studying thalamo-cortical connectivity during infancy, reported an overlap of the connectivity between primary visual, auditory and sensory-cortices within the same thalamus cluster in neonates, whereas distinct clusters were observed at one year of age. Contrarily, [START_REF] Toulmin | Specialization and integration of functional thalamocortical connectivity in the human infant[END_REF] described a more segregated thalamo-cortical connectivity from birth on, especially concerning the primary cortices. We may speculate that an initial stage of multisensory integration might facilitate learning by enhancing responses to cooccurring cues that need to be integrated into higher-level representations. For example, it might explain the observations that neonates visually recognize their mother only if they have been simultaneously exposed to her face and voice [START_REF] Sai | The role of the mother's voice in developing mother's face preference: Evidence for intermodal perception at birth[END_REF] and that 3-month-olds more robustly identify their mother vs. a stranger's face if her voice is simultaneously presented [START_REF] Burnham | Visual recognition of mother by young infants: facilitation by speech[END_REF].

This mesial-lateral division of functional correlation was similarly observed in the left and right hemispheres (see Fig. 4b) and may correspond to the cyto-architectonic division highlighted in adults by [START_REF] Weiner | The mid-fusiform sulcus: a landmark identifying both cytoarchitectonic and functional divisions of human ventral temporal cortex[END_REF][START_REF] Weiner | The cytoarchitecture of domain-specific regions in human high-level visual cortex[END_REF], and to the functional parcellation of activations revealed in older infants looking to faces/objects vs. places images [START_REF] Deen | Organization of high-level visual cortex in human infants[END_REF]. In both hemispheres, long-distance functional correlation was not limited to within-hemispheric areas, but also recruited contra-lateral heterotopic areas encompassing the inferior parietal and dorsal frontal regions. Although roughly similar, interesting differences were observed between the two hemispheres. Because the VWFA is systematically lateralized to the hemisphere involved in spoken language processing [START_REF] Pinel | Beyond hemispheric dominance: brain regions underlying the joint lateralization of language and arithmetic to the left hemisphere[END_REF][START_REF] Pinel | Genetic and environmental influences on the visual word form and fusiform face areas[END_REF][START_REF] Cai | The left ventral occipito-temporal response to words depends on language lateralization but not on visual familiarity[END_REF], which is usually in the left hemisphere, we examined the functional correlation at birth between three key regions of the spoken language network and ventral visual areas (Fig. 4c-i). The STG and MTG were strongly connected with both the ipsi-and contra-lateral visual regions, but the left MTG was significantly more related with the left lateral extra-striate visual region than its right counterpart. Thus, a stronger path already exists at birth between the spoken language network and the future location of the VWFA. It may explain why reading acquisition always recruits this particular location within the left ventral visual cortex [START_REF] Hannagan | Origins of the specialization for letters and numbers in ventral occipitotemporal cortex[END_REF].

The significant trans-hemispheric correlation between the left inferior frontal region and the right lateral ventral visual cortex was more surprising. This cluster is more mesial than the cluster discussed above (~ 24 vs -38 mm in our mask x-coordinates). Contrary to temporal regions, left and right inferior frontal regions were weakly correlated with ventral voxels in our mask except for this particular relation between this right cluster and the left frontal region. Note that in Fig. 4b, in which we studied the functional correlation of the left and right ventral temporal regions, we also observed a crossed relation yet in the reverse direction between the left temporal region and the right inferior frontal region. Does this functional correlation pattern correspond to a genuine, anatomically validated pathway? In rare patients with a developmental agenesis of the corpus callosum, sigmoid bundles connecting the right frontal region with the left contra-lateral occipital region (and, even more rarely, the converse pattern) have been detected by diffusion-tensor imaging [START_REF] Bénézit | Organising white matter in a brain without corpus callosum fibres[END_REF][START_REF] Tovar-Moll | Neuroplasticity in human callosal dysgenesis: a diffusion tensor imaging study[END_REF]). On one hand, these cases may correspond to aberrant bundles. On the other hand, the absence of a corpus callosum may reveal a smaller tract, otherwise masked by the large corpus callosum. Indeed, in a post-mortem tracing study in a single adult human brain, Di [START_REF] Virgilio | Direct interhemispheric visual input to human speech areas[END_REF] were able to identify the degenerating axons in a patient who died a few weeks after a stroke in the right inferior temporal region and reported inter-hemispheric mono-synaptic connections between the right inferior temporal region and the lefthemispheric Wernicke's and Broca's areas. Thus, interhemispheric connections can be heterotopic and relate widely different regions, in agreement with the present results. Our functional correlation pattern may also be related to a rare study of face processing in infancy by Tzourio-Mazoyer and collaborators (2002). Using positron emission tomography, these authors reported activations in left inferior frontal and left superior temporal cortex, together with right fusiform cortex, when 2-month-olds were looking at faces. Efficient connections between right-hemispheric fusiform face representations and left-hemispheric language areas might explain why vocal imitations can be so easily elicited when infants view articulating faces [START_REF] Bristow | Hearing faces: how the infant brain matches the face it sees with the speech it hears[END_REF][START_REF] Chen | Auditory-oral matching behavior in newborns[END_REF]. Future studies in larger groups should reveal whether these observations are anecdotal or reflect robust structural connections present in most of the neonates.

Although we studied a relatively small number of infants here, we observe a clear parcellation of the visual areas that speak to major theories on the origins of functional specialization in the adult human brain. The present results show that the neonate connectivity landscape is very far from a "blooming buzzing confusion" [START_REF] Embree | The principles of psychology[END_REF]. As such, our results refute theories of an extreme empiricist nature, that assume that the infant brain is largely devoid of initial architecture and that connectivity is initially diffuse, over-abundant and largely isotropic, only becoming refined at a later stage during infancy and childhood through exposure to a given environment [START_REF] Johnson | Interactive specialization: a domain-general framework for human functional brain development?[END_REF][START_REF] Elman | Rethinking innateness: a connectionist perspective on development[END_REF][START_REF] Quartz | The neural basis of cognitive development: a constructivist manifesto[END_REF]. Ventral visual cortex is not just "trained" by exposure to objects, faces, houses, letters or digits; on the contrary, it contains functionally specialized areas already present at the dawn of infancy [START_REF] Deen | Organization of high-level visual cortex in human infants[END_REF][START_REF] Biagi | BOLD response selective to flow-motion in very young infants[END_REF], supported by well-defined fiber tracts. Several previous studies have already shown how at a global scale, an organized functional and anatomical connectivity is present in the infant brain [START_REF] Fransson | Resting-state networks in the infant brain[END_REF][START_REF] Fransson | Spontaneous brain activity in the newborn brain during natural sleep-an fMRI study in infants born at full term[END_REF][START_REF] Fransson | The functional architecture of the infant brain as revealed by resting-state fMRI[END_REF][START_REF] Doria | Emergence of resting state networks in the preterm human brain[END_REF][START_REF] Smyser | Longitudinal analysis of neural network development in preterm infants[END_REF][START_REF] Smyser | Functional connectivity MRI in infants: exploration of the functional organization of the developing brain[END_REF][START_REF] Van Den Heuvel | The neonatal connectome during preterm brain development[END_REF]. Here we show a refined organization within a single region with a gradient of preferred relations with local and distant regions, in accordance with the hypothesis of an early organization ("protomaps") determined during gestation [START_REF] Rakic | Specification of cerebral cortical areas[END_REF][START_REF] Mahmoudzadeh | Functional maps at the onset of auditory inputs in very early preterm human neonates[END_REF][START_REF] Toulmin | Specialization and integration of functional thalamocortical connectivity in the human infant[END_REF]. These maps provide neonates with a blueprint on which learning capitalizes. This conclusion is in agreement with a recently published study in macaque newborns. Also using functional connectivity, Arcaro and Livingstone (2017) revealed a proto-organization in the macaque visual areas prior to the emergence of selective responses.

While indisputable differences in functional correlation were found between mesial and lateral regions, we did not recover a full adult-like organization. For instance, given the systematic location of word-and number-form areas in adults (VWFA and NFA), with the NFA systematically occupying a location lateral to the VWFA [START_REF] Shum | A brain area for visual numerals[END_REF][START_REF] Amalric | Origins of the brain networks for advanced mathematics in expert mathematicians[END_REF][START_REF] Hannagan | Origins of the specialization for letters and numbers in ventral occipitotemporal cortex[END_REF], we could have expected a differential pattern of connectivity, with the most lateral ventral visual voxels projecting preferentially to the inferior parietal lobe which has been implicated in number sense; however, such was not the case. That being said, our group was small merging two different groups of infants and our method was crude in the sense that we did not know precisely which cortical sectors, in our infants, would later become specialized for letters vs. digits, unlike what had been done in older children [START_REF] Saygin | Connectivity precedes function in the development of the visual word form area[END_REF]. Future studies with an improved methodology might reveal that the infant brain contains even finer-grained biases that predate the fine differences between VWFA and NFA seen in adults. Alternatively, it is possible that such patterns of functional specialization emerge through exposure and learning and that initially only a rough division between mesial and lateral regions is visible. In line with this idea, [START_REF] Deen | Organization of high-level visual cortex in human infants[END_REF] reported clear differences between fMRI activations for faces vs. places, but not for objects vs. faces in 4-6 month-old infants. Similarly, [START_REF] Gomez | Microstructural proliferation in human cortex is coupled with the development of face processing[END_REF] underlined the protracted development of the fusiform region sensitive to face relative to more mesial areas. The acquisition of longitudinal data in the same children from infancy to childhood, although extremely difficult, may be necessary to answer this question.

Fig. 1

 1 Fig. 1 Clustering of ventral visual cortex according to long-distance whole-brain functional correlations in newborns. a Four clusters were identified and are presented by different colors on axial slices of the infant template. b Partial correlation maps showing the distinct pat-

Fig. 2

 2 Fig. 2 Functional connectivity of the four ROIS located along the medial-lateral axis. a The four ROIs used in this analysis as seeds, ordered from mesial to lateral, are presented on a horizontal and coro-

Fig. 3 A

 3 Fig. 3 A mesial-lateral gradient of functional correlation in the infant ventral visual cortex. The connectivity patterns of each pair of the ROIs presented in Fig. 2 are compared two by two. Each column presents the t-contrast row ROI > column ROI on a lateral views and b ventral views. For illustration purposes, p value was set at p < 0.01,

Fig. 4

 4 Fig.4Hemispheric differences in long-distance functional correlations. a The mask containing all voxels used as seeds in the regression analysis is presented on a horizontal slice. To quantify to what extent functional correlations vary with the x-coordinate, a connectivity map was computed for each voxel in this mask and was subsequently entered into a regression analysis with the x coordinate as a regressor. b Functional correlation pattern as a function of the x-axis location of the seed in the left and right hemispheres. The more positive the value, the stronger the correlation with lateral voxels (e.g. the right inferior frontal region with a lateral left seed), whereas the more negative the value the stronger the correlation with mesial voxels (e.g. Heschl's gyrus on both sides). c-k Analyses of hemispheric biases.

Figure S1 Clustering

 S1 Figure S1

Acknowledgements Authors would like to thank Elvis Dohmatob, Isabelle Denghien, François Leroy, Jessica Dubois and Parvaneh Adibpour for their constructive remarks during this work.