We performed constrained spherical deconvolution (CSD) and particle-filtering tractography (PFT) with anatomical priors27,28 on DWI data of 39 healthy participants of the BIL&GIN database (Fig. 1A)29. All tractograms were registered to MNI space30, and the SLS was extracted for each individual leveraging the regions of interest (ROIs) of the JHU cortical atlas (Fig. 1B)31. Individual SLS were then concatenated. The N39 SLS was then parceled into sub-SLS using a cortex-to-cortex pairing approach, in which each frontal gyrus was systematically paired with every other gyrus of the same hemisphere as defined by the JHU template (Fig. 1C). sub-SLS from the right hemisphere were flipped and concatenated with their left hemisphere homologous (Fig. 1D). Out of the 119 possible ROI combinations, the cortex-to-cortex pairing approach led to the extraction of 96 non-empty sub-SLS tractograms. Of these, 45 contained more than 200 streamlines and were subsequently processed with QuickBundlesX clustering for the initial exclusion of implausible streamlines based on their geometrical features (Fig. 1D)32. These sub-SLS revealed a defined wiring of connections linking the 4 frontal gyri of the lateral convexity (i.e., IFG, MFG, superior frontal gyrus (SFG), and PrCG) to several regions covering the lateral surface of the parietal (i.e., angular gyrus (AG), postcentral gyrus (PoCG), superior parietal gyrus (SPG) and supramarginal gyrus (SMG)), temporal (i.e., ITG, MTG, STG and temporal pole (T-pole)) and occipital cortices (i.e., inferior (IOG) and middle occipital gyri (MOG)), and a restricted section of the basal temporal cortex (i.e., fusiform gyrus (FuG)) (Fig. 2). Connections involving the FuG, IOG, and T-pole were null when paired with the SFG. On the medial surface, connectivity was confined to the precuneus (PrCu). Conversely, 74 ROI combinations, consisting of 33 empty tractograms (Fig. 2, black circles) and 41 with less than 200 streamlines in both hemispheres (Fig. 2, gray circles), were deemed non-viable. They included the entire connectivity of the fronto-orbital region (lateral fronto-orbital gyrus (LFOG), middle fronto-orbital gyrus (MFOG), and rectus gyrus (RG)) as well as connections to the cuneus (Cu), lingual gyrus (LG), parahippocampal gyrus (PHG), entorhinal cortex (ENT), and superior occipital gyrus (SOG). Additionally, as noted earlier, connections between the SFG and the IOG, FuG, and T-pole were null.
Sub-SLS tractography templates description
This section is dedicated to the anatomical description of the 45 sub-SLS tractography templates obtained with the cortex-to-cortex pairing approach and after the removal of streamlines classified as implausible based on their geometrical features, as represented in Fig. 2. For this description, templates are grouped based on the frontal gyrus from which they originate.
Inferior frontal gyrus
Connections between the IFG and the inferior parietal lobule (IPL), namely AG and SMG, primarily originate from pars opercularis (pOp) and pars triangularis (pTri), with only a few streamlines reaching the dorsal and posterior pars orbitalis (pOrb). These connections course longitudinally within the deep WM of the central region of the hemisphere, running parallel the Sylvian fissure before resurfacing laterally to reach their parietal target. Connections to the AG mostly terminate in its anterior part, while those to the SMG are more evenly distributed. The IFG-PoCG template follows a similar but shorter course, parallelling the Sylvian fissure. For this connection, terminations in the IFG are more densely distributed in pOp and become sparser in pTri, with no streamline reaching pOrb. Posterior terminations cover only the most ventral part of the PoCG. The two templates representing connections to the dorso-medial parietal cortex, namely IFG-SPG and IFG-PrCu, follow a similar anatomical course to the connections described above up to the level of the PoCG, where they shift dorsally to reach their respective targets. These connections arise primarily from pOp, with only a few streamlines reaching pTri, particularly for the IFG-PrCu template, and reach anterior SPG and anterior PrCu. Connections between the IFG and the temporal lobe follow the same longitudinal course as the IFG-parietal connections up to the deep WM at the level of the AG. At this point, streamlines arch around the posterior ramus of the Sylvian fissure, then proceed ventrally with a transverse orientation. After entering the deep WM of the temporal lobe, they proceed towards their target temporal gyrus and start fanning to resurface laterally. While their fanning enlarges both anteriorly and posteriorly for connections to the ITG, MTG, STG, and FuG, connections to the T-pole follow a postero-anterior direction as they traverse the deep WM of the dorsal temporal lobe, parallelling the Sylvian fissure ventrally. Connections reaching the STG and the MTG arise uniformly from both pOp and pTri, with fewer streamlines originating from pOrb, while connections with the FuG, ITG, and T-pole almost exclusively originate from pOp. Posteriorly, connections to the FuG cover only the mid-segment of the ROI. While connections to the ITG, MTG, and STG distribute equally along the posterior ⅔ of the whole gyri, connections to the T-pole correspond to the anterior continuation of these connections. For what concerns connections with the occipital cortex, the IFG-MOG template follows a course similar to IFG-parietal connections, while the IFG-IOG template is more akin to IFG-temporal connections. Indeed, the IFG-MOG template has an overall longitudinal orientation, whereas IFG-IOG fibers arch around the Sylvian fissure, spanning the boundary between the temporal and occipital cortices before fanning out posteriorly to reach the anterior portion of IOG.
Middle frontal gyrus
Connections between the MFG and the IPL originate from the posterior ⅔ of the MFG ROI. After leaving the gyrus, streamlines travel posteriorly within the deep WM of the central region, maintaining a longitudinal orientation along the z-axis (streamlines course at the same height from the posterior end of the MFG to the AG and SMG). Streamlines resurface laterally at the level of the gyrus they terminate into. Posterior terminations are homogeneously distributed across both the AG and SMG ROIs. Templates of connections from the MFG to the PoCG, SPG, and PrCu exhibit a more compact course compared to connections to the IPL, and course slightly dorsal to them. Streamlines connecting to the SPG display more anterior terminations in MFG, covering the posterior ⅔ of the ROI similarly to the IPL connections. In contrast, streamlines to the PoCG and PrCu tend to terminate more posteriorly in the MFG, within the posterior third of the region. Streamlines of the MFG-PoCG template resurface laterally after passing the level of the central sulcus, with posterior terminations covering only the middle third and the inferior part of the dorsal third of the PoCG ROI. Connections from the MFG to the SPG and PrCu shift slightly vertically after passing the level of the postcentral sulcus, reaching the anterior part of the SPG and the entire PrCu ROI. Connections between the MFG and the ITG, MTG, and STG are more widespread, covering the posterior ⅔ of the MFG ROI. In contrast, connections to the FuG and the T-pole tend to arise more posteriorly in MFG, with the MFG-T-pole template covering only the very posterior aspect of the frontal ROI. After coursing longitudinally in the deep WM of the frontal and parietal lobes, MFG-temporal streamlines arch around the Sylvian fissure at the level of the IPL and enter the temporal cortex along a transverse course. Connections to the FuG terminate in the mid-segment of the FuG ROI, while connections to the ITG, MTG, and STG span the posterior ⅔ of the lateral temporal gyri. Connections to the T-pole traverse the deep WM of the superior temporal lobe and extend anteriorly. For what concerns connections with the occipital cortex, the MFG-MOG template resembles in its course MFG-parietal connections, while MFG-IOG conforms with MFG-temporal connections, arching around the Sylvian fissure and coursing ventrally with a transverse orientation before fanning posteriorly to reach the anterior end of the IOG ROI. The MFG-MOG template terminates in the anterior MOG.
Superior frontal gyrus
Connections between the SFG and the parietal lobe originate solely from the posterior third of the SFG. Templates reaching the PoCG, SPG, PrCu, and AG follow a longitudinal course through the deep WM of the dorsal central region. Contrarily, connections reaching the SMG leave the SFG ventrally, follow a main transverse orientation up to the level of the deep WM below the ventral PrCG and abruptly turn longitudinally before traveling posteriorly to the anterior SMG. Terminations in the PoCG cover only the most dorsal portion of the ROI, while connections to the SPG and PrCu cover the entire gyri. The AG is mainly covered in its antero-dorsal part. Connectivity between the SFG and the occipital cortex is represented by the sole SFG-MOG template. It resembles the SFG-AG template in its course, but extends further posteriorly to the anterior MOG. SFG-temporal connections involve only the posterior halves of the ITG, MTG, and STG. These streamlines leave the posterior half of the SFG ventrally and traverse the frontal lobe with a transverse orientation, up to the level of the dorsal IFG. Here, they abruptly switch course to a longitudinal orientation, reaching the deep WM of the IPL before arching into the temporal lobe. Streamlines fan laterally, both posteriorly and anteriorly, at the level of their respective temporal gyrus.
Precentral gyrus
The PrCG-PoCG template represents the only pattern of connectivity between adjacent gyri of the SLS. It consists of U-shaped fibers that span the entire extent of both gyri. Connections between the PrCG and the IPL arise from the ventral half of the PrCG ROI. After coursing briefly through the deep WM of the central region, streamlines resurface laterally, terminating homogeneously along the entire cortical surface of both the AG and SMG. The fanning of these streamlines is larger anteriorly to accommodate for the involvement of the PrCG ROI, which develops primarily along the z-axis. Connections between the PrCG and the SPG or the PrCu originate more dorsally in the PrCG than those to the IPL. Streamlines connecting to the SPG originate from the dorsal half of the PrCG, while those reaching the PrCu originate from the most dorsal third of this region. Streamlines of the PrCG-SPG template group in the deep WM below the PoCG and then turn slightly vertically to reach the anterior half of the SPG, while PrCG-PrCu connections remain dorsal and cover the anterior ⅔ of PrCu. Connections between the PrCG and the temporal lobe arise from the ventral half of the PrCG. Streamlines connecting to the ITG and the MTG are uniformly distributed along this portion of the PrCG ROI, those to the T-pole and the FuG arise more dorsally compared to the others, and those to STG seed only in the most ventral end of the region. Streamlines of the PrCG-temporal templates travel longitudinally up to the deep WM below the IPL, arch around the Sylvian fissure, and follow a transverse course until they reach the deep WM of their target gyrus. Connections to the FuG terminate in the mid-segment of the ROI, while those to the ITG, MTG, and STG cover the posterior ⅔ of the respective gyri. Connections to the T-pole extend more anteriorly, paralleling the Sylvian fissure ventrally, and cover the anterior third of the temporal lobe. As previously described for IFG and MFG, the PrCG-MOG template follows a similar course to PrCG-parietal connections, while streamlines of the PrCG-IOG template arch into the temporal lobe, crossing the boundary between temporal and occipital cortices, and fan posteriorly.
Sub-SLS ex-vivo dissection
Layer-by-layer Klingler microdissection of 6 left and 4 right hemispheres resulted in the identification of 24 dissection layers (henceforth referred to as epochs – abbreviated as “epo” when referring to a specific dissection layer) representing connections of the SLS. In some specimens, (abbreviated as “spc” when referring to a specific specimen) both fronto-parietal and fronto-temporal connectivity were exposed while proceeding with a latero-medial dissection approach (e.g., Fig. 3, spc-01, spc-06, spc-07, spc-10, spc-12, spc-15, and spc-19). Occasionally, fronto-parietal connectivity could be exposed on two subsequent epochs of the same specimen: a lateral one revealing more superficial and shorter fibers, and a medial one revealing deeper and longer ones (e.g., Fig. 3, spc-02 and spc-07). 3D digital models of each epoch were reconstructed with photogrammetry and aligned to the radiological space of each specimen25,33. A detailed anatomical analysis of the trajectory of the dorsal fibers leaving the frontal cortex and connecting to any other cortex of the same hemisphere was carried out on these models. Each sub-SLS connection visible on the texture of the models was manually annotated with the CloudCompare “Segment Tool” following a cortex-to-cortex pairing approach, compatible with the one adopted for the parcellation of tractography data. The identified sub-SLS annotated on the different epochs are listed and represented in Fig. 3.
Cortex-to-cortex annotation of the SLS connections on the specimens resulted in the identification of 22 sub-SLS, all of which were also previously reconstructed with tractography and represented by sub-SLS templates. These annotations depict connections between the IFG with the IPL (i.e., AG and SMG) as well as with the ITG and MTG. Connections linking the MFG and the PrCG to the lateral aspect of the parietal (i.e., AG, SMG, PoCG, and SPG) and temporal cortices (i.e., STG, MTG, and ITG) were also traced, alongside those to the IOG. The SFG-AG connection was the only sub-SLS identified from the SFG. Fibers reaching the FuG were visible just in one specimen, where they connected to the PrCG.
While the manual annotation of sub-SLS connections on the dissection models confirmed the existence of these connectivity patterns, the concurrent analysis of all the sub-SLS annotations across epochs revealed several organizational principles of the WM fibers of the SLS. Shorter fibers were found to course more superficially, while longer fibers ran more deeply. Specifically, fronto-parietal connections, being shorter by nature, were found on more superficial epochs compared to fronto-temporal fibers, which required the removal of overlaying tissue (e.g., spc-01 and spc-02 epo-06 vs. epo-07 in Fig. 3). While this observation aligns with the classical description of the SLF/AF complex34, we also noted that this organizational principle holds within individual connections defined by specific termination territories. For instance, MFG-AG connection in spc-02 was identified on two consecutive epochs (see Fig. 3, spc-02, epo-05 and epo-06 in orange). On epo-05, MFG-AG appears as a short connection involving only posterior MFG, while on epo-06, deeper fibers extend to reach anterior MFG. This latero-medial organizational pattern contributes to a layered architecture of WM fibers of the SLS, which also respects a ventro-dorsal organization. Ventral regions connect to other ventral regions through fibers coursing ventrally, while dorsal regions are connected to other dorsal regions through fibers coursing dorsally. This is particularly evident in fronto-parietal epochs. For instance, in epo-06 of spc-01 connections between the IFG and the AG course more ventrally compared to those between the MFG and the AG. Conversely, connections between the PrCG and the AG remain dorsal if originating from dorsal PrCG and ventral if arising from ventral PrCG. As seen in epo-07 of spc-02, this organizational layout also applies to fronto-temporal connections. In the case of PrCG-ITG fibers, those arising from the dorsal part of the PrCG after the splitting of the gyrus remain more external as they arch into the temporal lobe (i.e., farther from the posterior end of the Sylvian fissure) compared to fibers originating from the ventral part of the PrCG, which run closer to the Sylvian fissure. Thus, the ventro-dorsal organization in fronto-temporal connections translates into a mirrored pattern relative to the Sylvian fissure. Specifically, connections more ventral in the frontal cortex tend to connect to more dorsal regions in the temporal cortex, following the shape of the Sylvian fissure closely. Interestingly, as fronto-temporal fibers of the SLS arise more dorsally in the frontal cortex and thus course farther from the Sylvian fissure, the shift from a longitudinal course in the frontal and parietal cortices to a transverse course in the temporal cortex becomes less steep. Eventually, this results in an inversion in the fanning of temporal fibers – from more anteriorly directed for internal connections to more posteriorly directed for external connections. This trend culminates in fibers bending posteriorly to reach the occipital cortex, as shown in epo-05 of spc-16.
Sub-SLS templates anatomical evaluation with BraDiPho
All the 45 sub-SLS tractography templates were non-linearly registered to the radiological space of the 10 3D models of ex-vivo dissection for anatomical evaluation (Fig. 1E). Of these, 22 sub-SLS templates obtained with tractography were matched with their corresponding anatomical annotations from ex-vivo dissection, as described in the previous section. Quantitative overlap scores between tractography reconstructions and the corresponding manual annotations on the photogrammetric models of ex-vivo dissection are reported in Supplementary Data 1 and visualized in Supplementary Fig. 1. Across all matched sub-SLS templates, the average overlap between modalities was 88.15% ± 11.60% (mean ± SD; min 67.8%; max 100%). It is important to note that a degree of spatial mismatch may be expected, as tractography is registered to the T1 of the specimens before starting the dissection and minor shifts in tissue conformation occur during the dissection process, particularly as tissues relax following removal of overlying structures.
Based on the representation of the sub-SLS in the ex-vivo dissection models, the tractography templates were further refined. Initially, the selection of plausible streamlines was made without reference to the anatomical features of the sub-SLS, relying instead on plausibility judgment based on the geometrical characteristics of the streamlines. Moreover, this process involved broad clusters encompassing a wide array of streamlines, sometimes intermingling both plausible and implausible ones, possibly leading to the rejection of some plausible streamlines along with implausible ones. While rejecting some plausible streamlines would typically not be problematic in most cases due to the redundancy of tractography, with the correction for density bias in creation of the templates it could result in the loss of relevant reconstructions. Therefore, once anatomical information was obtained through the one-to-one comparison of in-vivo tractography and ex-vivo dissection data, the selection of plausible streamlines was re-performed in Tractome. This allowed for a more selective, connection-specific filtering of implausible streamlines35,36. The re-evaluation of the sub-SLS templates resulted in a more comprehensive and anatomically informed representation of WM connections, with improved coverage of their anatomical extent as suggested by the dissection data. These anatomically plausible and validated sub-SLS templates are shown in green in Fig. 1. A graphical summary of the anatomical evaluation process with BraDiPho is provided for each template in Supplementary Fig. 1.
The remaining 23 sub-SLS tractography templates were not found in our dissection data. These templates were classified as anatomically plausible but not validated or anatomically implausible based on the information gathered through dissection. The anatomically plausible but not validated templates (Fig. 2, yellow templates) represent connections that were reconstructed with tractography but could not be confirmed in our dissections, although their anatomical course remains consistent with anatomically plausible and validated connections. For instance, the existence of very medial connections, such as those from the MFG/PrCG/SFG to the PrCu, between SFG and dorsal PoCG, and between SFG and SPG, could not be proven, but they could not be excluded either. Indeed, these connections could not be exposed by our latero-medial dissection approach, which was not optimal for exposing such WM pathways. Similarly, connections to the MOG were not found in our dissections, as the dissection of purely longitudinal fibers stopped at the level of the AG. Dissection of the temporal cortex also did not extend to the T-pole, and rarely revealed connections to the basal surface of the hemisphere. The only specimen showing connections to the FuG (i.e., PrCG-FuG, see Fig. 3) was spc-06. Indeed, its unique anatomical configuration showing a “bulge” of the FuG interposing between the ITG and the IOG on the ventral lateral aspect of the convexity allowed the exposure of fibers connecting to this region through a latero-medial dissection approach. Additionally, connections with the IFG were only observed in 3 specimens, leaving anatomically plausible connections such as IFG-PoCG, IFG-STG, and IFG-IOG unproven.
On the other hand, we classified as anatomically implausible those sub-SLS templates showing abrupt changes of direction (Fig. 2, red templates). For instance, sub-SLS templates between the SFG and the temporal cortex (ITG, MTG, and STG) leave the posterior half of the SFG ventrally and traverse the posterior frontal cortex transversally up to the IFG before joining the other fronto-temporal connections of the SLS. Similarly, connections of the SFG-SMG template leave the SFG ventrally, travel vertically up to the level of the IFG, and then abruptly bend longitudinally toward the SMG. A similar but opposed shift in direction is seen in the IFG-SPG and IFG-PrCu templates, which travel longitudinally from the IFG up to the postcentral sulcus, before shifting vertically to reach the dorsal and medial parietal cortices. No such abrupt changes in direction were observed in the ex-vivo dissection of the SLS, therefore, these templates were classified as anatomically implausible. These directional shifts are likely caused by the intersection of longitudinal fibers with the anterior transverse system in the frontal lobe or with projection fibers of the pyramidal tract in the parietal lobe, leading the tractography algorithm to follow undesired directions.
sub-SLS templates replication study
The replication analysis of the sub-SLS tractography templates carried out on the three HCP test-retest tractography datasets distributed on brainlife.io37 with the Fast Streamline Search approach38 revealed a good agreement across tractography reconstructions derived from two different DWI acquisition protocols (i.e., BIL&GIN vs. HCP). Overall, we achieved 0.77 ± 0.27 average wDSC (i.e., dice coefficient on voxels weighted for streamline density) for test and test_run_2 and 0.78 ± 0.26 for retest. When considering only anatomically plausible and validated templates, mean wDSC spiked to 0.92 ± 0.08 for all HCP datasets, while it attested to 0.72 ± 0.29, 0.71 ± 0.3, and 0.71 ± 0.31 for test, test_run_2, and retest respectively, when considering anatomically plausible but not validated templates. On the other hand, wDSC dropped to 0.39 ± 0.22, 0.37 ± 0.19, and 0.43 ± 0.18 when considering sub-SLS templates deemed anatomically implausible. All measurements, including further streamline- and voxel-based similarity measurements on top of wDSC (i.e., bundle adjacency, dice coefficient, density correlation, overlap and overreach) are available in Supplementary Data 2, also reporting mean and standard deviations calculated considering all the templates together as well as filtered according to the classification of the templates (anatomically plausible and validated in green, anatomically plausible but not validated in yellow, and anatomically implausible in red). Single tractography reconstructions and voxel masks showing the overlap between the sub-SLS provided in the present work and sub-SLS extracted from each HCP dataset are displayed in Supplementary Fig. 2 with the respective wDSC. Of note, while on average wDSC dropped for anatomically implausible templates, the visualization provided in Supplementary Fig. 2 shows that implausible streamlines following the exact same course of the ones extracted on the N39 BIL&GIN dataset could also be found in the three HCP datasets. This reinforces the importance of relying on direct anatomical explorations to make tractography-based inferences on the architecture of WM pathways, independently from the characteristic of DWI and tractography.
sub-SLS templates characterization
A numerical description of each anatomically plausible sub-SLS tractography template—both validated and not validated (Fig. 2, green and yellow templates, respectively)—is provided through shape analysis (Supplementary Data 3)39. Additionally, we provide the mean coordinates of both anterior and posterior terminations for each of these templates in MNI space, with “anterior” indicating frontal terminations (Supplementary Data 4).
sub-SLS templates topological organization
The study of the SLS wiring through ex-vivo dissection and the annotation of different sub-SLS on photogrammetric models revealed key organizational principles of the WM fibers of this system. We analyzed the sub-SLS tractography templates to test for the occurrence of the same medio-lateral and dorso-ventral patterns of fiber distribution identified in the dissection study. The combined visualization of centroids of the sub-SLS templates originating in the IFG, MFG, and SFG shows that IFG connections start, course, and terminate more ventrally compared to MFG connections, which in turn are more ventral compared to SFG connections (Fig. 4, middle column, top row, in yellow, orange, and bordeaux, respectively). As previously described in the dissection data, this ordered layout is respected when streamlines bend into the temporal cortex, and defines their distance from the Sylvian fissure when arching around it (i.e., the more dorsal in the frontal cortex, the farther they course from the Sylvian fissure). Connections originating in the PrCG, which parallels the central sulcus and spans the most posterior part of the frontal lobe with a dorso-vental configuration, overlay to those arising from the other frontal gyri (Fig. 4, third column).
The dorso-ventral organization of fibers observed by considering all the frontal gyri together is respected also at the level of the single gyrus (Fig. 4, icons dedicated to IFG, MFG, SFG, and PrCG in the first and second columns). Connections arising in the PrCG provide an illustrative example. As it can be observed in Fig. 4, (second column, PrCG icon) fronto-parietal connections, depicted in green, are organized according to a ventro-dorsal layout, with a sequential organization orderly displaying, from ventral to dorsal, PrCG-SMG, PrCG-AG, PrCG-PoCG (whose centroid lies in the middle of the gyrus, as the connection spans the whole extent of both PrCG and PoCG), PoCG-SPG, and PrCG-PrCu. This organization respects the arrangement of cortical gyri in the parietal cortex, and drives the distribution of the fibers. The same principle applies to connections between the PrCG and both the lateral aspect of the temporal cortex and the occipital lobe. Connections from the PrCG to the STG terminate more ventrally in the PrCG compared to those reaching the MTG, which, in turn, are more ventral than those reaching the ITG. Accordingly, as we proceed from ventral to dorsal in the PrCG – and inversely from dorsal to ventral in the temporal cortex -, streamlines shift from coursing close to the Sylvian fissure to following a more distant trajectory. Similarly, connections to the IOG arise more ventrally in the PrCG compared to those reaching the MOG, and this dorso-ventral organization persists along their entire course. Connections to the T-pole and to the FuG demonstrate a similar organization, forming an additional deeper layer. Indeed, as seen in all the icons of Fig. 4 dedicated to each single gyrus, fronto-parietal connections, depicted in green, are the most superficial, followed by connections to the lateral temporal cortex and to the occipital lobe, which are themselves more superficial compared to connections reaching the anterior temporal lobe (i.e., T-pole) and the basal temporal cortex (i.e., FuG).
This medio-lateral layering reflects the length of streamlines, with shorter connections coursing more superficially and longer connections traveling deeper. This principle is shown in Fig. 5A, B, where streamlines are color-coded based on their distance from the cortex. More superficial fibers (yellow) cover shorter distances, while longer fibers (dark red), travel more deeply. Supplementary Fig. 3 provides the same color-coded representation for each sub-SLS template, further demonstrating that this organizational principle applies broadly across the SLS, even when considering only part of the system. Statistical analysis revealed a significant positive correlation between the length of each streamline and its distance from the cortex (Pearson’s r = 0.689, p < 0.001) (Fig. 5C).