Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

Differential Laminar Activation Dissociates Encoding and Retrieval in the Human Medial and Lateral Entorhinal Cortex

The hierarchically organized structures of the medial temporal lobe are critically important for episodic memory function. Accumulating evidence suggests dissociable information processing pathways are maintained throughout these structures including in the medial and lateral entorhinal cortex. Cortical layers provide an additional dimension of dissociation as the primary input to the hippocampus derives from layer 2 neurons in the entorhinal cortex, whereas the deeper layers primarily receive output from the hippocampus. Here, novel high-resolution T2-prepared functional MRI methods were successfully used to mitigate susceptibility artifacts typically affecting MRI signals in this region providing uniform sensitivity across the medial and lateral entorhinal cortex. During the performance of a memory task, healthy human subjects (age 25-33 years, mean age 28.2 ± 3.3 years, 4 female) showed differential functional activation in the superficial and deep layers of the entorhinal cortex associated with task-related encoding and retrieval conditions, respectively. The methods provided here offer an approach to probe layer-specific activation in normal cognition and conditions contributing to memory impairment.SIGNIFICANCE STATEMENT This study provides new evidence for differential neuronal activation in the superficial versus deep layers of the entorhinal cortex associated with encoding and retrieval memory processes, respectively, in cognitively normal adults. The study further shows that this dissociation can be observed in both the medial and the lateral entorhinal cortex. The study was achieved by using a novel functional MRI method allowing us to measure robust functional MRI signals in both the medial and lateral entorhinal cortex that was not possible in previous studies. The methodology established here in healthy human subjects lays a solid foundation for subsequent studies investigating layer-specific and region-specific changes in the entorhinal cortex associated with memory impairment in various conditions such as Alzheimer's disease.

Gradual changes in microarchitectural properties of cortex and juxtacortical white matter: Observed by anatomical and diffusion MRI

PURPOSE

Characterization of cerebral cortex is challenged by the complexity and heterogeneity of its cyto- and myeloarchitecture. This study evaluates quantitative MRI metrics, measured across two cortical depths and in subcortical white matter (WM) adjacent to cortex (juxtacortical WM), indicative of myelin content, neurite density, and diffusion microenvironment, for a comprehensive characterization of cortical microarchitecture.

METHODS

High-quality structural and diffusion MRI data (N = 30) from the Human Connectome Project were processed to compute myelin index, neurite density index, fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity from superficial cortex, deep cortex, and juxtacortical WM. The distributional patterns of these metrics were analyzed individually, correlated to one another, and were compared to established parcellations.

RESULTS

Our results supported that myeloarchitectonic and the coexisting cytoarchitectonic structures influence the diffusion properties of water molecules residing in cortex. Full cortical thickness showed myelination patterns similar to those previously observed in humans. Higher myelin indices with similar distributional patterns were observed in deep cortex whereas lower myelin indices were observed in superficial cortex. Neurite density index and other diffusion MRI derived parameters provided complementary information to myelination. Reliable and reproducible correlations were identified among the cortical microarchitectural properties and fiber distributional patterns in proximal WM structures.

CONCLUSION

We demonstrated gradual changes across the cortical sheath by assessing depth-specific cortical micro-architecture using anatomical and diffusion MRI. Mutually independent but coexisting features of cortical layers and juxtacortical WM provided new insights towards structural organizational units and variabilities across cortical regions and through depth.

Aberrant Cortical Layer Development of Brain Organoids Derived from Noonan Syndrome-iPSCs

Noonan syndrome (NS) is a genetic disorder mainly caused by gain-of-function mutations in Src homology region 2-containing protein tyrosine phosphatase 2 (SHP2). Although diverse neurological manifestations are commonly diagnosed in NS patients, the mechanisms as to how SHP2 mutations induce the neurodevelopmental defects associated with NS remain elusive. Here, we report that cortical organoids (NS-COs) derived from NS-induced pluripotent stem cells (iPSCs) exhibit developmental abnormalities, especially in excitatory neurons (ENs). Although NS-COs develop normally in their appearance, single-cell transcriptomic analysis revealed an increase in the EN population and overexpression of cortical layer markers in NS-COs. Surprisingly, the EN subpopulation co-expressing the upper layer marker SATB2 and the deep layer maker CTIP2 was enriched in NS-COs during cortical development. In parallel with the developmental disruptions, NS-COs also exhibited reduced synaptic connectivity. Collectively, our findings suggest that perturbed cortical layer identity and impeded neuronal connectivity contribute to the neurological manifestations of NS.

Relationship between Neuronal Damage/Death and Astrogliosis in the Cerebral Motor Cortex of Gerbil Models of Mild and Severe Ischemia and Reperfusion Injury

Neuronal loss (death) occurs selectively in vulnerable brain regions after ischemic insults. Astrogliosis is accompanied by neuronal death. It can change the molecular expression and morphology of astrocytes following ischemic insults. However, little is known about cerebral ischemia and reperfusion injury that can variously lead to damage of astrocytes according to the degree of ischemic injury, which is related to neuronal damage/death. Thus, the purpose of this study was to examine the relationship between damage to cortical neurons and astrocytes using gerbil models of mild and severe transient forebrain ischemia induced by blocking the blood supply to the forebrain for five or 15 min. Significant ischemia tFI-induced neuronal death occurred in the deep layers (layers V and VI) of the motor cortex: neuronal death occurred earlier and more severely in gerbils with severe ischemia than in gerbils with mild ischemia. Distinct astrogliosis was detected in layers V and VI. It gradually increased with time after both ischemiae. The astrogliosis was significantly higher in severe ischemia than in mild ischemia. The ischemia-induced increase of glial fibrillary acidic protein (GFAP; a maker of astrocyte) expression in severe ischemia was significantly higher than that in mild ischemia. However, GFAP-immunoreactive astrocytes were apparently damaged two days after both ischemiae. At five days after ischemiae, astrocyte endfeet around capillary endothelial cells were severely ruptured. They were more severely ruptured by severe ischemia than by mild ischemia. However, the number of astrocytes stained with S100 was significantly higher in severe ischemia than in mild ischemia. These results indicate that the degree of astrogliosis, including the disruption (loss) of astrocyte endfeet following ischemia and reperfusion in the forebrain, might depend on the severity of ischemia and that the degree of ischemia-induced neuronal damage may be associated with the degree of astrogliosis.

Signals from posterior parietal area 5 to motor cortex during locomotion

Area 5 of the parietal cortex is part of the "dorsal stream" cortical pathway which processes visual information for action. The signals that area 5 ultimately conveys to motor cortex, the main area providing output to the spinal cord, are unknown. We analyzed area 5 neuronal activity during vision-independent locomotion on a flat surface and vision-dependent locomotion on a horizontal ladder in cats focusing on corticocortical neurons (CCs) projecting to motor cortex from the upper and deeper cortical layers and compared it to that of neighboring unidentified neurons (noIDs). We found that upon transition from vision-independent to vision-dependent locomotion, the low discharge of CCs in layer V doubled and the proportion of cells with 2 bursts per stride tended to increase. In layer V, the group of 2-bursters developed 2 activity peaks that coincided with peaks of gaze shifts along the surface away from the animal, described previously. One-bursters and either subpopulation in supragranular layers did not transmit any clear unified stride-related signal to the motor cortex. Most CC group activities did not mirror those of their noID counterparts. CCs with receptive fields on the shoulder, elbow, or wrist/paw discharged in opposite phases with the respective groups of pyramidal tract neurons of motor cortex, the cortico-spinal cells.

Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI

Despite extensive studies detecting laminar functional magnetic resonance imaging (fMRI) signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across cortical regions remain to be fully investigated in both evoked and resting conditions at different brain states. Here, we developed a multislice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 μm) and temporal resolution (100 ms) in anesthetized rats. Across different trials, we detected either laminar-specific positive or negative blood-oxygen-level-dependent (BOLD) responses in the surrounding cortical region adjacent to the most activated cortex under the evoked condition. Specifically, in contrast to typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for trials with negative BOLD in adjacent regions, which indicated brain state-dependent laminar-fMRI responses based on corticocortical interaction. Also, in resting-state (rs-) fMRI study, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from caudal to rostral slices. In summary, our study provided a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity in evoked and resting-state MS-LS.

In Silico Assessment of Safety and Efficacy of Screw Placement for Pediatric Image-Guided Otologic Surgery

Introduction: Current high-accuracy image-guided systems for otologic surgery use fiducial screws for patient-to-image registration. Thus far, these systems have only been used in adults, and the safety and efficacy of the fiducial screw placement has not yet been investigated in the pediatric population.

Materials and Methods: In a retrospective study, CT image data of the temporal region from 11 subjects meeting inclusion criteria (8–48 months at the time of surgery) were selected, resulting in n = 20 sides. These datasets were investigated with respect to screw stability efficacy in terms of the cortical layer thickness, and safety in terms of the distance of potential fiducial screws to the dura mater or venous sinuses. All of these results are presented as distributions, thickness color maps, and with descriptive statistics. Seven regions within the temporal bone were analyzed individually. In addition, four fiducial screws per case with 4 mm thread-length were placed in an additively manufactured model according to the guidelines for robotic cochlear implantation surgery. For all these screws, the minimal distance to the dura mater or venous sinuses was measured, or if applicable how much they penetrated these structures.

Results: The cortical layer has been found to be mostly between 0.7–3.3 mm thick (from the 5^th to the 95^th percentile), while even thinner areas exist. The distance from the surface of the temporal bone to the dura mater or the venous sinuses varied considerably between the subjects and ranged mostly from 1.1–9.3 mm (from the 5^th to the 95^th percentile). From all 80 placed fiducial screws of 4 mm thread length in the pediatric subject younger than two years old, 22 touched or penetrated either the dura or the sigmoid sinus. The best regions for fiducial placement would be the mastoid area and along the petrous pyramid in terms of safety. In terms of efficacy, the parietal followed by the petrous pyramid, and retrosigmoid regions are most suited.

Conclusion: The current fiducial screws and the screw placement guidelines for adults are insufficiently safe or effective for pediatric patients.

Cross Laminar Traveling Components of Field Potentials due to Volume Conduction of Non-Traveling Neuronal Activity in Macaque Sensory Cortices

Field potentials (FPs) reflect neuronal activities in the brain, and often exhibit traveling peaks across recording sites. While traveling FPs are interpreted as propagation of neuronal activity, not all studies directly reveal such propagating patterns of neuronal activation. Neuronal activity is associated with transmembrane currents that form dipoles and produce negative and positive fields. Thereby, FP components reverse polarity between those fields and have minimal amplitudes at the center of dipoles. Although their amplitudes could be smaller, FPs are never flat even around these reversals. What occurs around the reversal has not been addressed explicitly, although those are rationally in the middle of active neurons. We show that sensory FPs around the reversal appeared with peaks traveling across cortical laminae in macaque sensory cortices. Interestingly, analyses of current source density did not depict traveling patterns but lamina-delimited current sinks and sources. We simulated FPs produced by volume conduction of a simplified 2 dipoles' model mimicking sensory cortical laminar current source density components. While FPs generated by single dipoles followed the temporal patterns of the dipole moments without traveling peaks, FPs generated by concurrently active dipole moments appeared with traveling components in the vicinity of dipoles by superimposition of individually non-traveling FPs generated by single dipoles. These results indicate that not all traveling FP are generated by traveling neuronal activity, and that recording positions need to be taken into account to describe FP peak components around active neuronal populations.SIGNIFICANCE STATEMENT Field potentials (FPs) generated by neuronal activity in the brain occur with fields of opposite polarity. Likewise, in the cerebral cortices, they have mirror-imaged waveforms in upper and lower layers. We show that FPs appear like traveling across the cortical layers. Interestingly, the traveling FPs occur without traveling components of current source density, which represents transmembrane currents associated with neuronal activity. These seemingly odd findings are explained using current source density models of multiple dipoles. Concurrently active, non-traveling dipoles produce FPs as mixtures of FPs produced by individual dipoles, and result in traveling FP waveforms as the mixing ratio depends on the distances from those dipoles. The results suggest that not all traveling FP components are associated with propagating neuronal activity.

A Protocol for Cortical Type Analysis of the Human Neocortex Applied on Histological Samples, the Atlas of Von Economo and Koskinas, and Magnetic Resonance Imaging

The human cerebral cortex is parcellated in hundreds of areas using neuroanatomy and imaging methods. Alternatively, cortical areas can be classified into few cortical types according to their degree of laminar differentiation. Cortical type analysis is based on the gradual and systematic variation of laminar features observed across the entire cerebral cortex in Nissl stained sections and has profound implications for understanding fundamental aspects of evolution, development, connections, function, and pathology of the cerebral cortex. In this protocol paper, we explain the general principles of cortical type analysis and provide tables with the fundamental features of laminar structure that are studied for this analysis. We apply cortical type analysis to the micrographs of the Atlas of the human cerebral cortex of von Economo and Koskinas and provide tables and maps with the areas of this Atlas and their corresponding cortical type. Finally, we correlate the cortical type maps with the T1w/T2w ratio from widely used reference magnetic resonance imaging scans. The analysis, tables and maps of the human cerebral cortex shown in this protocol paper can be used to predict patterns of connections between areas according to the principles of the Structural Model and determine their level in cortical hierarchies. Cortical types can also predict the spreading of abnormal proteins in neurodegenerative diseases to the level of cortical layers. In summary, cortical type analysis provides a theoretical and practical framework for directed studies of connectivity, synaptic plasticity, and selective vulnerability to neurologic and psychiatric diseases in the human neocortex.

Parallel RNA and DNA analysis after deep sequencing (PRDD-seq) reveals cell type-specific lineage patterns in human brain

Elucidating the lineage relationships among different cell types is key to understanding human brain development. Here we developed parallel RNA and DNA analysis after deep sequencing (PRDD-seq), which combines RNA analysis of neuronal cell types with analysis of nested spontaneous DNA somatic mutations as cell lineage markers, identified from joint analysis of single-cell and bulk DNA sequencing by single-cell MosaicHunter (scMH). PRDD-seq enables simultaneous reconstruction of neuronal cell type, cell lineage, and sequential neuronal formation ("birthdate") in postmortem human cerebral cortex. Analysis of two human brains showed remarkable quantitative details that relate mutation mosaic frequency to clonal patterns, confirming an early divergence of precursors for excitatory and inhibitory neurons, and an "inside-out" layer formation of excitatory neurons as seen in other species. In addition our analysis allows an estimate of excitatory neuron-restricted precursors (about 10) that generate the excitatory neurons within a cortical column. Inhibitory neurons showed complex, subtype-specific patterns of neurogenesis, including some patterns of development conserved relative to mouse, but also some aspects of primate cortical interneuron development not seen in mouse. PRDD-seq can be broadly applied to characterize cell identity and lineage from diverse archival samples with single-cell resolution and in potentially any developmental or disease condition.

Lamina-specific AMPA receptor dynamics following visual deprivation in vivo

Regulation of AMPA receptor (AMPAR) expression is central to synaptic plasticity and brain function, but how these changes occur in vivo remains elusive. Here, we developed a method to longitudinally monitor the expression of synaptic AMPARs across multiple cortical layers in awake mice using two-photon imaging. We observed that baseline AMPAR expression in individual spines is highly dynamic with more dynamics in primary visual cortex (V1) layer 2/3 (L2/3) neurons than V1 L5 neurons. Visual deprivation through binocular enucleation induces a synapse-specific and depth-dependent change of synaptic AMPARs in V1 L2/3 neurons, wherein deep synapses are potentiated more than superficial synapses. The increase is specific to L2/3 neurons and absent on apical dendrites of L5 neurons, and is dependent on expression of the AMPAR-binding protein GRIP1. Our study demonstrates that specific neuronal connections, across cortical layers and even within individual neurons, respond uniquely to changes in sensory experience.

Temporal refinement of sensory-evoked activity across layers in developing mouse barrel cortex

Rhythmic whisking behavior in rodents fully develops during a critical period about 2 weeks after birth, in parallel with the maturation of other sensory modalities and the onset of exploratory locomotion. How whisker-related sensory processing develops during this period in the primary somatosensory cortex (S1) remains poorly understood. Here, we characterized neuronal activity evoked by single- or dual-whisker stimulation patterns in developing S1, before, during and after the occurrence of active whisking. Employing multi-electrode recordings in all layers of barrel cortex in urethane-anesthetized mice, we find layer-specific changes in multi-unit activity for principal and neighboring barrel columns. While whisker stimulation evoked similar early responses (0-50 ms post-stimulus) across development, the late response (50-150 ms post-stimulus) decreased in all layers with age. Furthermore, peak onset times and the duration of the late response decreased in all layers across age groups. Responses to paired-pulse stimulation showed increases in spiking precision and in paired-pulse ratios in all cortical layers during development. Sequential activation of two neighboring whiskers with varying stimulus intervals evoked distinct response profiles in the activated barrel columns, depending on the direction and temporal separation of the stimuli. In conclusion, our findings indicate that the temporal sharpening of sensory-evoked activity coincides with the onset of active whisking.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

Meningeal Bmps Regulate Cortical Layer Formation

Neuronal connectivity in the cortex is determined by the laminar positioning of neurons. An important determinant of laminar positioning is likely to be the control of leading process behavior during migration, maintaining their tips directed toward the pia. In this study, we provide evidence that pial bone morphogenetic protein (Bmp) signaling regulates cortical neuronal migration during cortical layer formation. Specific disruption of pial Bmp ligands impaired the positioning of early-born neurons in the deep layer; further, cell-autonomous inhibition of Smad4, a core nuclear factor mediating Bmp signaling, in the cortical radial glial cells or postmitotic cortical neurons also produced neuronal migration defects that blurred the cortical layers. We found that leading processes were abnormal and that this was accompanied by excess dephosphorylated cofilin-1, an actin-severing protein, in Smad4 mutant neurons. This suggested that regulation of cofilin-1 might transduce Bmp signaling in the migrating neurons. Ectopic expression of a phosphorylation-defective form of cofilin-1 in the late-born wild-type neurons led them to stall in the deep layer, similar to the Smad4 mutant neurons. Expression of a phosphomimetic variant of cofilin-1 in the Smad4 mutant neurons rescued the migration defects. This suggests that cofilin-1 activity underlies Bmp-mediated cortical neuronal migration. This study shows that cofilin-1 mediates pial Bmp signaling during the positioning of cortical neurons and the formation of cortical layers.

Laminar Module Cascade from Layer 5 to 6 Implementing Cue-to-Target Conversion for Object Memory Retrieval in the Primate Temporal Cortex

The cerebral cortex computes through the canonical microcircuit that connects six stacked layers; however, how cortical processing streams operate in vivo, particularly in the higher association cortex, remains elusive. By developing a novel MRI-assisted procedure that reliably localizes recorded single neurons at resolution of six individual layers in monkey temporal cortex, we show that transformation of representations from a cued object to a to-be-recalled object occurs at the infragranular layer in a visual cued-recall task. This cue-to-target conversion started in layer 5 and was followed by layer 6. Finally, a subset of layer 6 neurons exclusively encoding the sought target became phase-locked to surrounding field potentials at theta frequency, suggesting that this coordinated cell assembly implements cortical long-distance outputs of the recalled target. Thus, this study proposes a link from local computation spanning laminar modules of the temporal cortex to the brain-wide network for memory retrieval in primates.

Identity of neocortical layer 4 neurons is specified through correct positioning into the cortex

Many cell-intrinsic mechanisms have been shown to regulate neuronal subtype specification in the mammalian neocortex. However, how much cell environment is crucial for subtype determination still remained unclear. Here, we show that knockdown of Protocadherin20 (Pcdh20), which is expressed in post-migratory neurons of layer 4 (L4) lineage, caused the cells to localize in L2/3. The ectopically positioned "future L4 neurons" lost their L4 characteristics but acquired L2/3 characteristics. Knockdown of a cytoskeletal protein in the future L4 neurons, which caused random disruption of positioning, also showed that those accidentally located in L4 acquired the L4 characteristics. Moreover, restoration of positioning of the Pcdh20-knockdown neurons into L4 rescued the specification failure. We further suggest that the thalamocortical axons provide a positional cue to specify L4 identity. These results suggest that the L4 identity is not completely determined at the time of birth but ensured by the surrounding environment after appropriate positioning.

Laminar Differences in Dendritic Structure of Pyramidal Neurons in the Juvenile Rat Somatosensory Cortex

Pyramidal cell structure varies between different cortical areas and species, indicating that the cortical circuits that these cells participate in are likely to be characterized by different functional capabilities. Structural differences between cortical layers have been traditionally reported using either the Golgi method or intracellular labeling, but the structure of pyramidal cells has not previously been systematically analyzed across all cortical layers at a particular age. In the present study, we investigated the dendritic architecture of complete basal arbors of pyramidal neurons in layers II, III, IV, Va, Vb, and VI of the hindlimb somatosensory cortical region of postnatal day 14 rats. We found that the characteristics of basal dendritic morphologies are statistically different in each cortical layer. The variations in size and branching pattern that exist between pyramidal cells of different cortical layers probably reflect the particular functional properties that are characteristic of the cortical circuit in which they participate. This new set of complete basal dendritic arbors of 3D-reconstructed pyramidal cell morphologies across each cortical layer will provide new insights into interlaminar information processing in the cerebral cortex.

Linking contemporary high resolution magnetic resonance imaging to the von Economo legacy: A study on the comparison of MRI cortical thickness and histological measurements of cortical structure

The cerebral cortex is a distinctive part of the mammalian nervous system, displaying a spatial variety in cyto-, chemico-, and myelinoarchitecture. As part of a rich history of histological findings, pioneering anatomists von Economo and Koskinas provided detailed mappings on the cellular structure of the human cortex, reporting on quantitative aspects of cytoarchitecture of cortical areas. Current day investigations into the structure of human cortex have embraced technological advances in Magnetic Resonance Imaging (MRI) to assess macroscale thickness and organization of the cortical mantle in vivo. However, direct comparisons between current day MRI estimates and the quantitative measurements of early anatomists have been limited. Here, we report on a simple, but nevertheless important cross-analysis between the histological reports of von Economo and Koskinas on variation in thickness of the cortical mantle and MRI derived measurements of cortical thickness. We translated the von Economo cortical atlas to a subdivision of the commonly used Desikan-Killiany atlas (as part of the FreeSurfer Software package and a commonly used parcellation atlas in studies examining MRI cortical thickness). Next, values of "width of the cortical mantle" as provided by the measurements of von Economo and Koskinas were correlated to cortical thickness measurements derived from high-resolution anatomical MRI T1 data of 200+ subjects of the Human Connectome Project (HCP). Cross-correlation revealed a significant association between group-averaged MRI measurements of cortical thickness and histological recordings (r = 0.54, P < 0.001). Further validating such a correlation, we manually segmented the von Economo parcellation atlas on the standardized Colin27 brain dataset and applied the obtained three-dimensional von Economo segmentation atlas to the T1 data of each of the HCP subjects. Highly consistent with our findings for the mapping to the Desikan-Killiany regions, cross-correlation between in vivo MRI cortical thickness and von Economo histology-derived values of cortical mantle width revealed a strong positive association (r = 0.62, P < 0.001). Linking today's state-of-the-art T1-weighted imaging to early histological examinations our findings indicate that MRI technology is a valid method for in vivo assessment of thickness of human cortex.

Neural dynamics and information representation in microcircuits of motor cortex

The brain has to analyze and respond to external events that can change rapidly from time to time, suggesting that information processing by the brain may be essentially dynamic rather than static. The dynamical features of neural computation are of significant importance in motor cortex that governs the process of movement generation and learning. In this paper, we discuss these features based primarily on our recent findings on neural dynamics and information coding in the microcircuit of rat motor cortex. In fact, cortical neurons show a variety of dynamical behavior from rhythmic activity in various frequency bands to highly irregular spike firing. Of particular interest are the similarity and dissimilarity of the neuronal response properties in different layers of motor cortex. By conducting electrophysiological recordings in slice preparation, we report the phase response curves (PRCs) of neurons in different cortical layers to demonstrate their layer-dependent synchronization properties. We then study how motor cortex recruits task-related neurons in different layers for voluntary arm movements by simultaneous juxtacellular and multiunit recordings from behaving rats. The results suggest an interesting difference in the spectrum of functional activity between the superficial and deep layers. Furthermore, the task-related activities recorded from various layers exhibited power law distributions of inter-spike intervals (ISIs), in contrast to a general belief that ISIs obey Poisson or Gamma distributions in cortical neurons. We present a theoretical argument that this power law of in vivo neurons may represent the maximization of the entropy of firing rate with limited energy consumption of spike generation. Though further studies are required to fully clarify the functional implications of this coding principle, it may shed new light on information representations by neurons and circuits in motor cortex.

Monocular inhibition reveals temporal and spatial changes in gene expression in the primary visual cortex of marmoset

We investigated the time course of the expression of several activity-dependent genes evoked by visual inputs in the primary visual cortex (V1) in adult marmosets. In order to examine the rapid time course of activity-dependent gene expression, marmosets were first monocularly inactivated by tetrodotoxin (TTX), kept in darkness for two days, and then exposed to various length of light stimulation. Activity-dependent genes including HTR1B, HTR2A, whose activity-dependency were previously reported by us, and well-known immediate early genes (IEGs), c-FOS, ZIF268, and ARC, were examined by in situ hybridization. Using this system, first, we demonstrated the ocular dominance type of gene expression pattern in V1 under this condition. IEGs were expressed in columnar patterns throughout layers II–VI of all the tested monocular marmosets. Second, we showed the regulation of HTR1B and HTR2A expressions by retinal spontaneous activity, because HTR1B and HTR2A mRNA expressions sustained a certain level regardless of visual stimulation and were inhibited by a blockade of the retinal activity with TTX. Third, IEGs dynamically changed its laminar distribution from half an hour to several hours upon a stimulus onset with the unique time course for each gene. The expression patterns of these genes were different in neurons of each layer as well. These results suggest that the regulation of each neuron in the primary visual cortex of marmosets is subjected to different regulation upon the change of activities from retina. It should be related to a highly differentiated laminar structure of marmoset visual systems, reflecting the functions of the activity-dependent gene expression in marmoset V1.

Single axon branching analysis in rat thalamocortical projection from the anteroventral thalamus to the granular retrosplenial cortex

The granular retrosplenial cortex (GRS) in the rat has a distinct microcolumn-type structure. The apical tufts of dendritic bundles at layer I, which are formed by layer II neurons, co-localize with patches of thalamic terminations from anteroventral (AV) thalamic nucleus. To further understand this microcolumn-type structure in the GRS, one of remaining questions is whether this structure extends into other layers, such as layers III/IV. Other than layer I, previous tracer injection study showed that AV thalamic nucleus also projects to layer III/IV in the GRS. In this study, we examined the morphology of branches in the GRS from the AV thalamus in single axon branch resolution in order to determine whether AV axon branches in layer III/IV are branches of axons with extensive branch in layer I, and, if so, whether the extent of these arborizations in layer III/IV vertically matches with that in layer I. For this purpose, we used a small volume injection of biotinylated dextran-amine into the AV thalamus and reconstructing labeled single axon branches in the GRS. We found that the AV axons consisted of heterogeneous branching types. Type 1 had extensive arborization occurring only in layer Ia. Type 2 had additional branches in III/IV. Types 1 and 2 had extensive ramifications in layer Ia, with lateral extensions within the previously reported extensions of tufts from single dendritic bundles (i.e., 30–200 μm; mean 78 μm). In type 2 branches, axon arborizations in layer III/IV were just below to layer Ia ramifications, but much wider (148–533 μm: mean, 341 μm) than that in layer Ia axon branches and dendritic bundles, suggesting that layer-specific information transmission spacing existed even from the same single axons from the AV to the GRS. Thus, microcolumn-type structure in the upper layer of the GRS was not strictly continuous from layer I to layer IV. How each layer and its components interact each other in different spatial scale should be solved future.

Temporal dynamics and spatial specificity of arterial and venous blood volume changes during visual stimulation: implication for BOLD quantification

Determination of compartment-specific cerebral blood volume (CBV) changes is important for understanding neurovascular physiology and quantifying blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI). In isoflurane-anesthetized cats, we measured the spatiotemporal responses of arterial CBV (CBV(a)) and total CBV (CBV(t)) induced by a 40-second visual stimulation, using magnetization transfer (MT)-varied BOLD and contrast-agent fMRI techniques at 9.4 T. To determine the venous CBV (CBV(v)) change, we calculated the difference between CBV(t) and CBV(a) changes. The dynamic response of CBV(a) was an order of magnitude faster than that of CBV(v), while the magnitude of change under steady-state conditions was similar between the two. Following stimulation offset, ΔCBV(a) showed small poststimulus undershoots, while ΔCBV(v) slowly returned to baseline. The largest CBV(a) and CBV(t) response occurred after 10 seconds of simulation in cortical layer 4, which we identified as the stripe of Gennari by T(1)-weighted MRI. The CBV(v) response, however, was not specific across the cortical layers during the entire stimulation period. Our data indicate that rapid, more-specific arterial vasodilation is followed by slow, less-specific venous dilation. Our finding implies that the contribution of CBV(v) changes to BOLD signals is significant for long, but not short, stimulation periods.

Axonal netrin-Gs transneuronally determine lamina-specific subdendritic segments

Axons from a distinct group of neurons make contact with dendritic trees of target neurons in clearly segregated and laminated patterns, thereby forming functional units for processing multiple inputs of information in the vertebrate central nervous system. Whether and how dendrites acquire lamina-specific properties corresponding to each pathway is not known. We show here that vertebrate-specific membrane-anchored members of the UNC-6/netrin family, netrin-G1 and -G2, organize the lamina/pathway-specific differentiation of dendrites. Netrin-G1 and -G2 distribute on axons of different pathways and specifically interact with receptors NGL-1 and -2, respectively. In the hippocampus, parietal cortex, and piriform cortex, NGL-1 is concentrated in the dendritic segments corresponding to the lamina-specific termination of netrin-G1-positive axons, and NGL-2 is concentrated in distinct dendritic segments corresponding to the termination of netrin-G2-positive axons. In netrin-G1- and -G2-deficient mice, in which axonal path-finding is normal, the segmental distribution of NGL-1 and -2 is selectively disrupted, and the individual receptors are diffused along the dendrites. These findings indicate that transneuronal interactions of netrin-Gs and their specific receptors provide a molecular basis for the axonal innervation-dependent mechanism of postsynaptic membrane organization, and provide insight into the formation of the laminar structure within the dendrites.

Spatial and temporal structure of receptive fields in primate somatosensory area 3b: effects of stimulus scanning direction and orientation

This is the third in a series of studies of the neural representation of tactile spatial form in somatosensory cortical area 3b of the alert monkey. We previously studied the spatial structure of >350 fingerpad receptive fields (RFs) with random-dot patterns scanned in one direction () and at varying velocities (). Those studies showed that area 3b RFs have a wide range of spatial structures that are virtually unaffected by changes in scanning velocity. In this study, 62 area 3b neurons were studied with three to eight scanning directions (58 with four or more directions). The data from all three studies are described accurately by an RF model with three components: (1) a single, central excitatory region of short duration, (2) one or more inhibitory regions, also of short duration, that are adjacent to and nearly synchronous with the excitation, and (3) a region of inhibition that overlaps the excitation partially or totally and is temporally delayed with respect to the first two components. The mean correlation between the observed RFs and the RFs predicted by this three-component model was 0.81. The three-component RFs also predicted orientation sensitivity and preferred orientation to a scanned bar accurately. The orientation sensitivity was determined most strongly by the intensity of the coincident RF inhibition in relation to the excitation. Both orientation sensitivity and this ratio were stronger in the supragranular and infragranular layers than in layer IV.