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.

Patient-specific mutation of Dync1h1 in mice causes brain and behavioral deficits.

AIMS: Cytoplasmic dynein heavy chain (DYNC1H1) is a multi-subunit protein complex that provides motor force for movement of cargo on microtubules and traffics them back to the soma. In humans, mutations along the DYNC1H1 gene result in intellectual disabilities, cognitive delays, and neurologic and motor deficits. The aim of the study was to generate a mouse model to a newly identified de novo heterozygous DYNC1H1 mutation, within a functional ATPase domain (c9052C > T(P3018S)), identified in a child with motor deficits, and intellectual disabilities. RESULTS: P3018S heterozygous (HET) knockin mice are viable; homozygotes are lethal. Metabolic and EchoMRI™ testing show that HET mice have a higher metabolic rate, are more active, and have less body fat compared to wildtype mice. Neurobehavioral studies show that HET mice perform worse when traversing elevated balance beams, and on the negative geotaxis test. Immunofluorescent staining shows neuronal migration abnormalities in the dorsal and lateral neocortex with heterotopia in layer I. Neuron-subtype specific transcription factors CUX1 and CTGF identified neurons from layers II/III and VI respectively in cortical layer I, and abnormal pyramidal neurons with MAP2+ dendrites projecting downward from the pial surface. CONCLUSION: The HET mice are a good model for the motor deficits seen in the child, and highlights the importance of cytoplasmic dynein in the maintenance of cortical function and dendritic orientation relative to the pial surface. Our results are discussed in the context of other dynein mutant mice and in relation to clinical presentation in humans with DYNC1H1 mutations.

Synaptic integration of somatosensory and motor cortical inputs onto spiny projection neurons of mice caudoputamen.

The basal ganglia play pivotal roles in motor control and cognitive functioning. These nuclei are embedded in an anatomical loop: cortex to basal ganglia to thalamus back to cortex. We focus here on an essential synapse for descending control, from cortical layer 5 (L5) onto the GABAergic spiny projection neurons (SPNs) of the caudoputamen (CP). We employed genetic labeling to distinguish L5 neurons from somatosensory (S1) and motor (M1) cortices in large volume serial electron microscopy and electrophysiology datasets to better detail these inputs. First, M1 and S1 synapses showed a strong preference to innervate the spines of SPNs and rarely contacted aspiny cells, which are likely to be interneurons. Second, L5 inputs commonly converge from both areas onto single SPNs. Third, compared to unlabeled terminals in CP, those labeled from M1 and S1 show ultrastructural hallmarks of strong driver synapses: They innervate larger spines that were more likely to contain a spine apparatus, more often had embedded mitochondria, and more often contacted multiple targets. Finally, these inputs also demonstrated driver-like functional properties: SPNs responded to optogenetic activation from S1 and M1 with large EPSP/Cs that depressed and were dependent on ionotropic but not metabotropic receptors. Together, our findings suggest that individual SPNs integrate driver input from multiple cortical areas with implications for how the basal ganglia relay cortical input to provide inhibitory innervation of motor thalamus.

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.

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.

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.

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.

Modeling Reveals Human-Rodent Differences in H-Current Kinetics Influencing Resonance in Cortical Layer 5 Neurons.

While our understanding of human neurons is often inferred from rodent data, inter-species differences between neurons can be captured by building cellular models specifically from human data. This includes understanding differences at the level of ion channels and their implications for human brain function. Thus, we here present a full spiking, biophysically detailed multi-compartment model of a human layer 5 (L5) cortical pyramidal cell. Model development was primarily based on morphological and electrophysiological data from the same human L5 neuron, avoiding confounds of experimental variability. Focus was placed on describing the behavior of the hyperpolarization-activated cation (h-) channel, given increasing interest in this channel due to its role in pacemaking and differentiating cell types. We ensured that the model exhibited post-inhibitory rebound spiking considering its relationship with the h-current, along with other general spiking characteristics. The model was validated against data not used in its development, which highlighted distinctly slower kinetics of the human h-current relative to the rodent setting. We linked the lack of subthreshold resonance observed in human L5 neurons to these human-specific h-current kinetics. This work shows that it is possible and necessary to build human-specific biophysical neuron models in order to understand human brain dynamics.

Creative Destruction: A Basic Computational Model of Cortical Layer Formation.

One of the most characteristic properties of many vertebrate neural systems is the layered organization of different cell types. This cytoarchitecture exists in the cortex, the retina, the hippocampus, and many other parts of the central nervous system. The developmental mechanisms of neural layer formation have been subject to substantial experimental efforts. Here, we provide a general computational model for cortical layer formation in 3D physical space. We show that this multiscale, agent-based model, comprising two distinct stages of apoptosis, can account for the wide range of neuronal numbers encountered in different cortical areas and species. Our results demonstrate the phenotypic richness of a basic state diagram structure. Importantly, apoptosis allows for changing the thickness of one layer without automatically affecting other layers. Therefore, apoptosis increases the flexibility for evolutionary change in layer architecture. Notably, slightly changed gene regulatory dynamics recapitulate the characteristic properties observed in neurodevelopmental diseases. Overall, we propose a novel computational model using gene-type rules, exhibiting many characteristics of normal and pathological cortical development.

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.

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.

Magnetization transfer weighted EPI facilitates cortical depth determination in native fMRI space.

The increased availability of ultra-high field scanners provides an opportunity to perform fMRI at sub-millimeter spatial scales and enables in vivo probing of laminar function in the human brain. In most previous studies, the definition of cortical layers, or depths, is based on an anatomical reference image that is collected by a different acquisition sequence and exhibits different geometric distortion compared to the functional images. Here, we propose to generate the anatomical image with the fMRI acquisition technique by incorporating magnetization transfer (MT) weighted imaging. Small flip angle binomial pulse trains are used as MT preparation, with a flexible duration (several to tens of milliseconds), which can be applied before each EPI segment without constraining the acquisition length (segment or slice number). The method's feasibility was demonstrated at 7T for coverage of either a small slab or the near-whole brain at 0.8 mm isotropic resolution. Tissue contrast was found to be similar to that obtained with a state-of-art anatomical reference based on MP2RAGE. This MT-weighted EPI image allows an automatic reconstruction of the cortical surface to support laminar analysis in native fMRI space, obviating the need for distortion correction and registration.

Integrated VASO and perfusion contrast: A new tool for laminar functional MRI.

Earlier research in cats has shown that both cerebral blood volume (CBV) and cerebral blood flow (CBF) can be used to identify layer-dependent fMRI activation with spatial specificity superior to gradient-echo blood-oxygen-level-dependent (BOLD) contrast (Jin and Kim, 2008a). CBF contrast of perfusion fMRI at ultra-high field has not been widely applied in humans to measure laminar activity due to its low sensitivity, while CBV contrast for fMRI using vascular space occupancy (VASO) has been successfully used. However, VASO can be compromised by interference of blood in-flow effects and a temporally limited acquisition window around the blood-nulling time point. Here, we proposed to use DANTE (Delay Alternating with Nutation for Tailored Excitation) pulse trains combined with 3D-EPI to acquire an integrated VASO and perfusion (VAPER) contrast. The signal origin of the VAPER contrast was theoretically evaluated with respect to its CBV and CBF contributions using a four-compartment simulation model. The feasibility of VAPER to measure layer-dependent activity was empirically investigated in human primary motor cortex at 7 â€‹T. We demonstrated this new tool, with its highly specified functional layer profile, robust reproducibility, and improved sensitivity, to allow investigation of layer-specific cortical functions.

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.

Differential gene regulatory plasticity between upper and lower layer cortical excitatory neurons.

Neocortical projection neurons consist of intracortical connected upper layer (UL, layer II-IV) neurons and subcortical connected lower layer (LL, layer V-VI) neurons. Afferent activity from the thalamus regulates layer-specific gene expression during postnatal development, which is critical for the formation of proper neocortical cytoarchitecture. Here, we show that activity-dependent gene regulation is confined to UL cortical neurons, but not LL neurons, and that this distinction is likely due to epigenetic modifications of chromatin. We found that the immediate early genes (IEGs), EGR1 and c-FOS, are downregulated in all cortical laminar layers in the absence of afferent activity in vivo. Transcriptional assays demonstrated that EGR1 and c-FOS are able to bind to the promoters of UL- and LL-specific genes to induce transcription. Furthermore, we discovered that LL neurons express higher levels of heterochromatin markers, such as H3K9m3 and H4K20m3, compared to UL neurons. Our results suggest that differential epigenetic modifications of chromatin is an intrinsic mechanism that underlies the different sensitivities of cortical neurons to activity-dependent gene regulation.

FIP4/Arfophilin-2 plays overlapping but distinct roles from FIP3/Arfophilin-1 in neuronal migration during cortical layer formation.

The class II Rab11 family-interacting proteins, FIP3 and FIP4, also termed Arfophilin-1 and Arfophilin-2, respectively, are endosomal proteins that function as dual effector proteins for Rab11 and ADP ribosylation factor (Arf) small GTPases. In the present study, we examined the expression and role of FIP4 in neuronal migration during cerebral layer formation. FIP4 mRNA was first weakly detected in post-mitotic migrating neurons in the upper intermediate zone, and expression was markedly increased in the cortical layer. Exogenously expressed FIP4 protein was localized to subpopulations of EEA1- and syntaxin 12-positive endosomes in migrating neurons, and was partially colocalized with FIP3. Knockdown of FIP4 by in utero electroporation significantly stalled transfected neurons in the lower cortical layer and decreased the speed of neuronal migration in the upper intermediate zone and in the cortical plate compared with control small hairpin RNA (shRNA)-transfected neurons. Furthermore, co-transfection of shRNA-resistant wild-type FIP4, but not wild type FIP3 or FIP4 mutants lacking the binding region for Rab11 or Arf, significantly improved the disturbed cortical layer formation caused by FIP4 knockdown. Collectively, our findings suggest that FIP4 and FIP3 play overlapping but distinct roles in neuronal migration downstream of Arf and Rab11 during cortical layer formation.

Unbiased Quantification of Subplate Neuron Loss following Neonatal Hypoxia-Ischemia in a Rat Model.

BACKGROUND: Cellular targets of neonatal hypoxia-ischemia (HI) include both oligodendrocyte and neuronal lineages with differences in the patterns of vulnerable cells depending upon the developmental stage at which the injury occurs. Injury to the developing white matter is a characteristic feature of human preterm brain injury. Data are accumulating, however, for neuronal injury in the developing cerebral cortex. In the most widely used rodent model of preterm HI brain injury, conflicting data have been reported regarding the sensitivity of subplate neurons to early neonatal HI, with some reports of selective vulnerability and others that find no increased loss of subplate neurons in comparison with other cortical layers. Methods used to identify subplate neurons and quantify their numbers vary across studies. OBJECTIVE: To use recently developed cortical layer-specific markers quantified with definitive stereologic methods to determine the magnitude and specificity of subplate neuron cell loss following neonatal HI in a rodent model. METHODS: Postnatal day 2 (P2) rats underwent right common carotid artery coagulation followed by 2-3 h of hypoxia (5.6% oxygen). Categorically moderately injured brains were stained with subplate and cortical layer III-V markers (Complexin3 and Foxp1, respectively) at P8 and P21 (Foxp1 only). An Optical Fractionator was used to quantify subplate and middle/lower cortical neuronal numbers and these were compared across groups (naive control, hypoxia hemisphere, and HI hemisphere). RESULTS: Following HI at P2 in rats, the total Complexin3-expressing subplate neuron number decreases significantly in the HI hemisphere compared with naive controls or hypoxia alone (HI vs. control 26,747 ± 7,952 vs. 35,468 ± 8,029, p = 0.04; HI vs. hypoxia, 26,747 ± 7,952 vs. 40,439 ± 7,363, p = 0.003). In contrast, the total Foxp1-expressing layer III-V cell number did not differ across the 3 conditions at P8 (HI vs. control 1,195,085 ± 436,609 vs. 1,234,640 ± 178,540, p = 0.19; HI vs. hypoxia, 1,195,085 ± 436,609 vs. 1,289,195 ± 468,941, p = 0.35) and at P21 (HI vs. control 1,265,190 ± 48,089 vs. 1,195,632 ± 26,912, p = 0.19; HI vs. hypoxia, 1,265,190 ± 48,089 vs. 1,309,563 ± 41,669, p = 0.49). CONCLUSIONS: There is significant biological variability inherent in both the subplate neuron cell number and the pattern and severity of cortical injury following HI at P2 in rats. Despite this variability, the subplate neuron cell number is lower following P2 HI in animals with mild or moderate cortical injury, whereas the middle-to-lower-layer cortical neuronal number is unchanged. In more severe cases, neurons are lost from the lower cortical layers, suggesting a relative vulnerability of subplate neurons.

The Regulation of Corticofugal Fiber Targeting by Retinal Inputs.

Corticothalamic projection systems arise from 2 main cortical layers. Layer V neurons project exclusively to higher-order thalamic nuclei, while layer VIa fibers project to both first-order and higher-order thalamic nuclei. During early postnatal development, layer VIa and VIb fibers accumulate at the borders of the dorsal lateral geniculate nucleus (dLGN) before they innervate it. After neonatal monocular enucleation or silencing of the early retinal activity, there is premature entry of layer VIa and VIb fibers into the dLGN contralateral to the manipulation. Layer V fibers do not innervate the superficial gray layer of the superior colliculus during the first postnatal week, but also demonstrate premature entry to the contralateral superficial gray layer following neonatal enucleation. Normally, layer V driver projections to the thalamus only innervate higher-order nuclei. Our results demonstrate that removal of retinal input from the dLGN induces cortical layer V projections to aberrantly enter, arborize, and synapse within the first-order dLGN. These results suggest that there is cross-hierarchical corticothalamic plasticity after monocular enucleation. Cross-hierarchical rewiring has been previously demonstrated in the thalamocortical system (Pouchelon et al. 2014), and now we provide evidence for cross-hierarchical corticothalamic rewiring after loss of the peripheral sensory input.