Histology-derived volumetric annotation of the human hippocampal subfields in postmortem MRI

Clinical phenotypes of Alzheimer's disease: investigating atrophy patterns and their pathological correlates

BACKGROUND

In Alzheimer's disease (AD), MRI atrophy patterns can distinguish between amnestic (typical) and non-amnestic (atypical) clinical phenotypes and are increasingly used for diagnosis and outcome measures in clinical trials. However, understanding how protein accumulation and other key features of neurodegeneration influence these imaging measurements, are lacking. The current study aimed to assess regional MRI patterns of cortical atrophy across clinical AD phenotypes, and their association with amyloid-beta (Aβ), phosphorylated tau (pTau), neuro-axonal degeneration and microvascular deterioration.

METHODS

Post-mortem in-situ 3DT1 3 T-MRI data was obtained from 33 AD (17 typical, 16 atypical) and 16 control brain donors. Additionally, ante-mortem 3DT1 3 T-MRI scans of brain donors were collected if available. Regional volumes were obtained from MRI scans using an atlas based parcellation software. Eight cortical brain regions were selected from formalin-fixed right hemispheres of brain donors and then immunostained for Aβ, pTau, neurofilament light, and collagen IV. Group comparisons and volume-pathology associations were analyzed using linear mixed models corrected for age, sex, post-mortem delay, and intracranial volume.

RESULTS

Compared to controls, both typical and atypical AD showed volume loss in the temporo-occipital cortex, while typical AD showed additional volume loss in the parietal cortex. Posterior cingulate volume was lower in typical AD compared to atypical AD (- 6.9%, p = 0.043). In AD, a global positive association between MRI cortical volume and Aβ load (βs = 0.21, p = 0.010), and a global negative association with NfL load (βs = - 0.18, p = 0.018) were observed. Regionally, higher superior parietal gyrus volume was associated with higher Aβ load in typical AD (βs = 0.47, p = 0.004), lower middle frontal gyrus volume associated with higher NfL load in atypical AD (βs = - 0.50, p < 0.001), and lower hippocampal volume associated with higher COLIV load in typical AD (βs = - 1.69, p < 0.001). Comparing post-mortem with ante-mortem scans showed minimal volume differences at scan-intervals within 2 years, highlighting the translational aspect of this study.

CONCLUSION

For both clinical phenotypes, cortical volume is affected by Aβ and neuro-axonal damage, but in opposing directions. Differences in volume-pathology relationships between clinical phenotypes are region-specific. The findings of this study could improve the interpretation of MRI datasets in heterogenous AD cohorts, both in research and clinical settings.

Mesoscale connectivity of the human hippocampus and fimbria revealed by ex vivo diffusion MRI

The human hippocampus is essential to cognition and emotional processing. Its function is defined by its connectivity. Although some pathways have been well-established, our knowledge about anterior-posterior connectivity and the distribution of fibers from major fiber bundles remains limited. Mesoscale (250 μm isotropic acquisition, upsampled to 125 μm) resolution MR images of the human temporal lobe afforded a detailed visualization of fiber tracts, including those that related anterior-posterior substructures defined as subregions (head, body, tail) and subfields (cornu ammonis 1-3, dentate gyrus) of the hippocampus. Fifty pathways were dissected between the head and body, highlighting an intricate mesh of connectivity between these two subregions. Along the body subregion, 12 lamellae were identified based on morphology and the presence of interlamellar fibers that appear to connect neighboring lamellae at the edge of the external limb of the granule cell layer (GCL). Translamellar fibers (i.e. longitudinal fibers crossing more than 2 lamellae) were also evident at the edge of the internal limb of the GCL. The dentate gyrus of the body was the main site of connectivity with the fimbria. Unique pathways were dissected within the fimbria that connected the body of the hippocampus with the amygdala and the temporal pole. A topographical segregation within the fimbria was determined by fibers' hippocampal origin, illustrating the importance of mapping the spatial distribution of fibers. Elucidating the detailed structural connectivity of the hippocampus is crucial to develop better diagnostic markers of neurological and psychiatric conditions, as well as to devise novel surgical interventions.

Curved multiplanar reformatting allows the accurate histological delineation of hippocampal subfields

BACKGROUND

Hippocampal subfields perform specific roles in normal cognitive functioning and have distinct vulnerabilities in neurological disorders. However, measurement of subfields with MRI is technically difficult in the head and tail of the hippocampus. Recent studies have utilized curved multiplanar reconstruction (CMPR) to improve subfield visualization in the head and tail, but this method has not yet been applied to histological data.

METHODS

We utilized BigBrain data, an open-source database of serially sectioned histological data for our analyses. The left hippocampus was segmented according to histological criteria by two raters in order to evaluate intra- and inter-rater reliability of histology-based segmentation throughout the long axis. Segmentation according to our previous protocol for the hippocampal body was then compared to these histological measurements to evaluate for histological validity. Agreement between segmentations was evaluated using Dice similarity coefficients (DSCs).

RESULTS

Intra-rater reliability (DSCs) of histological segmentation was excellent for all subfields: CA1 (0.8599), CA2 (0.7586), CA3/CA4/DG (0.8907), SLM (0.9123), subiculum (0.8149). Similarly, inter-rater reliability analysis demonstrated excellent agreement (DSCs) for all subfield locations: CA1 (0.8203), CA2 (0.7253), CA3/CA4/DG (0.8439), SLM (0.8700), subiculum (0.7794). Finally, histological accuracy (DSCs) for our previous protocol was excellent for all subfields: CA1 (0.8821), CA2 (0.8810), CA3/CA4/DG (0.9802), SLM (0.9879), subiculum (0.8774). When subfields in the hippocampus head, body, and tail were analyzed independently, DSCs also showed excellent agreement.

CONCLUSIONS

CMPR allows reliable subfield segmentation based on histological criteria throughout the hippocampal head, body, and tail. Our previous protocol for the hippocampal body can be applied to provide histologically valid subfield measurements throughout the entire hippocampal long axis.

Pathologic Substrates of Structural Brain Network Resilience and Topology in Parkinson Disease Decedents

BACKGROUND AND OBJECTIVES

In Parkinson disease (PD), α-synuclein spreading through connected brain regions leads to neuronal loss and brain network disruptions. With diffusion-weighted imaging (DWI), it is possible to capture conventional measures of brain network organization and more advanced measures of brain network resilience. We aimed to investigate which neuropathologic processes contribute to regional network topologic changes and brain network resilience in PD.

METHODS

Using a combined postmortem MRI and histopathology approach, PD and control brain donors with available postmortem in situ 3D T1-weighted MRI, DWI, and brain tissue were selected from the Netherlands Brain Bank and Normal Aging Brain Collection Amsterdam. Probabilistic tractography was performed, and conventional network topologic measures of regional eigenvector centrality and clustering coefficient, and brain network resilience (change in global efficiency upon regional node failure) were calculated. PSer129 α-synuclein, phosphorylated-tau, β-amyloid, neurofilament light-chain immunoreactivity, and synaptophysin density were quantified in 8 cortical regions. Group differences and correlations were assessed with rank-based nonparametric tests, with age, sex, and postmortem delay as covariates.

RESULTS

Nineteen clinically defined and pathology-confirmed PD (7 F/12 M, 81 ± 7 years) and 15 control (8 F/7 M, 73 ± 9 years) donors were included. With regional conventional measures, we found lower eigenvector centrality only in the parahippocampal gyrus in PD (d = -1.08, 95% CI 0.003-0.010, p = 0.021), which did not associate with underlying pathology. No differences were found in regional clustering coefficient. With the more advanced measure of brain network resilience, we found that the PD brain network was less resilient to node failure of the dorsal anterior insula compared with the control brain network (d = -1.00, 95% CI 0.0012-0.0015, p = 0.018). This change was not directly driven by neuropathologic processes within the dorsal anterior insula or in connected regions but was associated with higher Braak α-synuclein staging (rs = -0.40, p = 0.036).

DISCUSSION

Although our cohort might suffer from selection bias, our results highlight that regional network disturbances are more complex to interpret than previously believed. Regional neuropathologic processes did not drive regional topologic changes, but a global increase in α-synuclein pathology had a widespread effect on brain network reorganization in PD.

Automated deep learning segmentation of high-resolution 7 Tesla postmortem MRI for quantitative analysis of structure-pathology correlations in neurodegenerative diseases

Postmortem MRI allows brain anatomy to be examined at high resolution and to link pathology measures with morphometric measurements. However, automated segmentation methods for brain mapping in postmortem MRI are not well developed, primarily due to limited availability of labeled datasets, and heterogeneity in scanner hardware and acquisition protocols. In this work, we present a high-resolution dataset of 135 postmortem human brain tissue specimens imaged at 0.3 mm3 isotropic using a T2w sequence on a 7T whole-body MRI scanner. We developed a deep learning pipeline to segment the cortical mantle by benchmarking the performance of nine deep neural architectures, followed by post-hoc topological correction. We evaluate the reliability of this pipeline via overlap metrics with manual segmentation in 6 specimens, and intra-class correlation between cortical thickness measures extracted from the automatic segmentation and expert-generated reference measures in 36 specimens. We also segment four subcortical structures (caudate, putamen, globus pallidus, and thalamus), white matter hyperintensities, and the normal appearing white matter, providing a limited evaluation of accuracy. We show generalizing capabilities across whole-brain hemispheres in different specimens, and also on unseen images acquired at 0.28 mm3 and 0.16 mm3 isotropic T2*w fast low angle shot (FLASH) sequence at 7T. We report associations between localized cortical thickness and volumetric measurements across key regions, and semi-quantitative neuropathological ratings in a subset of 82 individuals with Alzheimer's disease (AD) continuum diagnoses. Our code, Jupyter notebooks, and the containerized executables are publicly available at the project webpage (https://pulkit-khandelwal.github.io/exvivo-brain-upenn/).

Three-dimensional histology reveals dissociable human hippocampal long-axis gradients of Alzheimer's pathology

INTRODUCTION

Three-dimensional (3D) histology analyses are essential to overcome sampling variability and understand pathological differences beyond the dissection axis. We present Path2MR, the first pipeline allowing 3D reconstruction of sparse human histology without a magnetic resonance imaging (MRI) reference. We implemented Path2MR with post-mortem hippocampal sections to explore pathology gradients in Alzheimer's disease.

METHODS

Blockface photographs of brain hemisphere slices are used for 3D reconstruction, from which an MRI-like image is generated using machine learning. Histology sections are aligned to the reconstructed hemisphere and subsequently to an atlas in standard space.

RESULTS

Path2MR successfully registered histological sections to their anatomic position along the hippocampal longitudinal axis. Combined with histopathology quantification, we found an expected peak of tau pathology at the anterior end of the hippocampus, whereas amyloid-beta (Aβ) displayed a quadratic anterior-posterior distribution.

CONCLUSION

Path2MR, which enables 3D histology using any brain bank data set, revealed significant differences along the hippocampus between tau and Aβ.

HIGHLIGHTS

Path2MR enables three-dimensional (3D) brain reconstruction from blockface dissection photographs. This pipeline does not require dense specimen sampling or a subject-specific magnetic resonance (MR) image. Anatomically consistent mapping of hippocampal sections was obtained with Path2MR. Our analyses revealed an anterior-posterior gradient of hippocampal tau pathology. In contrast, the peak of amyloid-beta (Aβ) deposition was closer to the hippocampal body.

Immunohistochemical field parcellation of the human hippocampus along its antero-posterior axis

The primate hippocampus includes the dentate gyrus, cornu ammonis (CA), and subiculum. CA is subdivided into four fields (CA1-CA3, plus CA3h/hilus of the dentate gyrus) with specific pyramidal cell morphology and connections. Work in non-human mammals has shown that hippocampal connectivity is precisely patterned both in the laminar and longitudinal axes. One of the main handicaps in the study of neuropathological semiology in the human hippocampus is the lack of clear laminar and longitudinal borders. The aim of this study was to explore a histochemical segmentation of the adult human hippocampus, integrating field (medio-lateral), laminar, and anteroposterior longitudinal patterning. We provide criteria for head-body-tail field and subfield parcellation of the human hippocampus based on immunodetection of Rabphilin3a (Rph3a), Purkinje-cell protein 4 (PCP4), Chromogranin A and Regulation of G protein signaling-14 (RGS-14). Notably, Rph3a and PCP4 allow to identify the border between CA3 and CA2, while Chromogranin A and RGS-14 give specific staining of CA2. We also provide novel histological data about the composition of human-specific regions of the anterior and posterior hippocampus. The data are given with stereotaxic coordinates along the longitudinal axis. This study provides novel insights for a detailed region-specific parcellation of the human hippocampus useful for human brain imaging and neuropathology.

The prosubiculum in the human hippocampus: A rostrocaudal, feature-driven, and systematic approach

The hippocampal subfield prosubiculum (ProS), is a conserved neuroanatomic region in mouse, monkey, and human. This area lies between CA1 and subiculum (Sub) and particularly lacks consensus on its boundaries; reports have varied on the description of its features and location. In this report, we review, refine, and evaluate four cytoarchitectural features that differentiate ProS from its neighboring subfields: (1) small neurons, (2) lightly stained neurons, (3) superficial clustered neurons, and (4) a cell sparse zone. ProS was delineated in all cases (n = 10). ProS was examined for its cytoarchitectonic features and location rostrocaudally, from the anterior head through the body in the hippocampus. The most common feature was small pyramidal neurons, which were intermingled with larger pyramidal neurons in ProS. We quantitatively measured ProS pyramidal neurons, which showed (average, width at pyramidal base = 14.31 µm, n = 400 per subfield). CA1 neurons averaged 15.57 µm and Sub neurons averaged 15.63 µm, both were significantly different than ProS (Kruskal-Wallis test, p < .0001). The other three features observed were lightly stained neurons, clustered neurons, and a cell sparse zone. Taken together, these findings suggest that ProS is an independent subfield, likely with distinct functional contributions to the broader interconnected hippocampal network. Our results suggest that ProS is a cytoarchitecturally varied subfield, both for features and among individuals. This diverse architecture in features and individuals for ProS could explain the long-standing complexity regarding the identification of this subfield.

Microstructural integrity of the locus coeruleus and its tracts reflect noradrenergic degeneration in Alzheimer's disease and Parkinson's disease

BACKGROUND

Degeneration of the locus coeruleus (LC) noradrenergic system contributes to clinical symptoms in Alzheimer's disease (AD) and Parkinson's disease (PD). Diffusion magnetic resonance imaging (MRI) has the potential to evaluate the integrity of the LC noradrenergic system. The aim of the current study was to determine whether the diffusion MRI-measured integrity of the LC and its tracts are sensitive to noradrenergic degeneration in AD and PD.

METHODS

Post-mortem in situ T1-weighted and multi-shell diffusion MRI was performed for 9 AD, 14 PD, and 8 control brain donors. Fractional anisotropy (FA) and mean diffusivity were derived from the LC, and from tracts between the LC and the anterior cingulate cortex, the dorsolateral prefrontal cortex (DLPFC), the primary motor cortex (M1) or the hippocampus. Brain tissue sections of the LC and cortical regions were obtained and immunostained for dopamine-beta hydroxylase (DBH) to quantify noradrenergic cell density and fiber load. Group comparisons and correlations between outcome measures were performed using linear regression and partial correlations.

RESULTS

The AD and PD cases showed loss of LC noradrenergic cells and fibers. In the cortex, the AD cases showed increased DBH + immunoreactivity in the DLPFC compared to PD cases and controls, while PD cases showed reduced DBH + immunoreactivity in the M1 compared to controls. Higher FA within the LC was found for AD, which was correlated with loss of noradrenergic cells and fibers in the LC. Increased FA of the LC-DLPFC tract was correlated with LC noradrenergic fiber loss in the combined AD and control group, whereas the increased FA of the LC-M1 tract was correlated with LC noradrenergic neuronal loss in the combined PD and control group. The tract alterations were not correlated with cortical DBH + immunoreactivity.

CONCLUSIONS

In AD and PD, the diffusion MRI-detected alterations within the LC and its tracts to the DLPFC and the M1 were associated with local noradrenergic neuronal loss within the LC, rather than noradrenergic changes in the cortex.

Contributions of hippocampal subfields and subregions to episodic memory performance in healthy cognitive aging

In the present study we investigated whether hippocampal subfield (cornu ammonis 1-3, dentate gyrus, and subiculum) and anteroposterior hippocampal subregion (head,body, and tail) volumes can predict episodic memory function using high-field high resolution structural magnetic resonance imaging (MRI). We recruited 126 healthy participants (18-85 years). MRI datasets were collected on a 4.7 T system. Participants were administered the Wechsler Memory Scale (WMS-IV) to evaluate episodic memory function. Structural equation modeling was used to test the relationship between studied variables. We found that the volume of the dentate gyrus subfield and posterior hippocampus (body) showed a significant direct effect on visuospatial memory performance; additionally, an indirect effect of age on visuospatial memory mediated through these hippocampal subfield/subregion was significant. Logical and verbal memory were not significantly associated with hippocampal subfield or subregion volumes.

WarpDrive: Improving spatial normalization using manual refinements

Spatial normalization-the process of mapping subject brain images to an average template brain-has evolved over the last 20+ years into a reliable method that facilitates the comparison of brain imaging results across patients, centers & modalities. While overall successful, sometimes, this automatic process yields suboptimal results, especially when dealing with brains with extensive neurodegeneration and atrophy patterns, or when high accuracy in specific regions is needed. Here we introduce WarpDrive, a novel tool for manual refinements of image alignment after automated registration. We show that the tool applied in a cohort of patients with Alzheimer's disease who underwent deep brain stimulation surgery helps create more accurate representations of the data as well as meaningful models to explain patient outcomes. The tool is built to handle any type of 3D imaging data, also allowing refinements in high-resolution imaging, including histology and multiple modalities to precisely aggregate multiple data sources together.

Three-dimensional histology reveals dissociable human hippocampal long axis gradients of Alzheimer's pathology

INTRODUCTION

Three-dimensional (3D) histology analyses are essential to overcome sampling variability and understand pathological differences beyond the dissection axis. We present Path2MR, the first pipeline allowing 3D reconstruction of sparse human histology without an MRI reference. We implemented Path2MR with post-mortem hippocampal sections to explore pathology gradients in Alzheimer's Disease.

METHODS

Blockface photographs of brain hemisphere slices are used for 3D reconstruction, from which an MRI-like image is generated using machine learning. Histology sections are aligned to the reconstructed hemisphere and subsequently to an atlas in standard space.

RESULTS

Path2MR successfully registered histological sections to their anatomical position along the hippocampal longitudinal axis. Combined with histopathology quantification, we found an expected peak of tau pathology at the anterior end of the hippocampus, while amyloid-β displayed a quadratic anterior-posterior distribution.

CONCLUSION

Path2MR, which enables 3D histology using any brain bank dataset, revealed significant differences along the hippocampus between tau and amyloid-β.

Mapping mesoscale connectivity within the human hippocampus

The connectivity of the hippocampus is essential to its functions. To gain a whole system view of intrahippocampal connectivity, ex vivo mesoscale (100 μm isotropic resolution) multi-shell diffusion MRI (11.7T) and tractography were performed on entire post-mortem human right hippocampi. Volumetric measurements indicated that the head region was largest followed by the body and tail regions. A unique anatomical organization in the head region reflected a complex organization of the granule cell layer (GCL) of the dentate gyrus. Tractography revealed the volumetric distribution of the perforant path, including both the tri-synaptic and temporoammonic pathways, as well as other well-established canonical connections, such as Schaffer collaterals. Visualization of the perforant path provided a means to verify the borders between the pro-subiculum and CA1, as well as between CA1/CA2. A specific angularity of different layers of fibers in the alveus was evident across the whole sample and allowed a separation of afferent and efferent connections based on their origin (i.e. entorhinal cortex) or destination (i.e. fimbria) using a cluster analysis of streamlines. Non-canonical translamellar connections running along the anterior-posterior axis were also discerned in the hilus. In line with "dentations" of the GCL, mossy fibers were bunching together in the sagittal plane revealing a unique lamellar organization and connections between these. In the head region, mossy fibers projected to the origin of the fimbria, which was distinct from the body and tail region. Mesoscale tractography provides an unprecedented systems view of intrahippocampal connections that underpin cognitive and emotional processing.

The Brainbox-a tool to facilitate correlation of brain magnetic resonance imaging features to histopathology

Magnetic resonance imaging (MRI) has limitations in identifying underlying tissue pathology, which is relevant for neurological diseases such as multiple sclerosis, stroke or brain tumours. However, there are no standardized methods for correlating MRI features with histopathology. Thus, here we aimed to develop and validate a tool that can facilitate the correlation of brain MRI features to corresponding histopathology. For this, we designed the Brainbox, a waterproof and MRI-compatible 3D printed container with an integrated 3D coordinate system. We used the Brainbox to acquire post-mortem ex vivo MRI of eight human brains, fresh and formalin-fixed, and correlated focal imaging features to histopathology using the built-in 3D coordinate system. With its built-in 3D coordinate system, the Brainbox allowed correlation of MRI features to corresponding tissue substrates. The Brainbox was used to correlate different MR image features of interest to the respective tissue substrate, including normal anatomical structures such as the hippocampus or perivascular spaces, as well as a lacunar stroke. Brain volume decreased upon fixation by 7% (P = 0.01). The Brainbox enabled degassing of specimens before scanning, reducing susceptibility artefacts and minimizing bulk motion during scanning. In conclusion, our proof-of-principle experiments demonstrate the usability of the Brainbox, which can contribute to improving the specificity of MRI and the standardization of the correlation between post-mortem ex vivo human brain MRI and histopathology. Brainboxes are available upon request from our institution.

Mixed methodology in human brain research: integrating MRI and histology

Postmortem magnetic resonance imaging (MRI) can provide a bridge between histological observations and the in vivo anatomy of the human brain. Approaches aimed at the co-registration of data derived from the two techniques are gaining interest. Optimal integration of the two research fields requires detailed knowledge of the tissue property requirements for individual research techniques, as well as a detailed understanding of the consequences of tissue fixation steps on the imaging quality outcomes for both MRI and histology. Here, we provide an overview of existing studies that bridge between state-of-the-art imaging modalities, and discuss the background knowledge incorporated into the design, execution and interpretation of postmortem studies. A subset of the discussed challenges transfer to animal studies as well. This insight can contribute to furthering our understanding of the normal and diseased human brain, and to facilitate discussions between researchers from the individual disciplines.

Pentad: A reproducible cytoarchitectonic protocol and its application to parcellation of the human hippocampus

Introduction

The hippocampus is integral for learning and memory and is targeted by multiple diseases. Neuroimaging approaches frequently use hippocampal subfield volumes as a standard measure of neurodegeneration, thus making them an essential biomarker to study. Collectively, histologic parcellation studies contain various disagreements, discrepancies, and omissions. The present study aimed to advance the hippocampal subfield segmentation field by establishing the first histology based parcellation protocol, applied to n = 22 human hippocampal samples.

Methods

The protocol focuses on five cellular traits observed in the pyramidal layer of the human hippocampus. We coin this approach the pentad protocol. The traits were: chromophilia, neuron size, packing density, clustering, and collinearity. Subfields included were CA1, CA2, CA3, CA4, prosubiculum, subiculum, presubiculum, parasubiculum, as well as the medial (uncal) subfields Subu, CA1u, CA2u, CA3u, and CA4u. We also establish nine distinct anterior-posterior levels of the hippocampus in the coronal plane to document rostrocaudal differences.

Results

Applying the pentad protocol, we parcellated 13 subfields at nine levels in 22 samples. We found that CA1 had the smallest neurons, CA2 showed high neuronal clustering, and CA3 displayed the most collinear neurons of the CA fields. The border between presubiculum and subiculum was staircase shaped, and parasubiculum had larger neurons than presubiculum. We also demonstrate cytoarchitectural evidence that CA4 and prosubiculum exist as individual subfields.

Discussion

This protocol is comprehensive, regimented and supplies a high number of samples, hippocampal subfields, and anterior-posterior coronal levels. The pentad protocol utilizes the gold standard approach for the human hippocampus subfield parcellation.

Neurofilament light chain is increased in the parahippocampal cortex and associates with pathological hallmarks in Parkinson's disease dementia

BACKGROUND

Increased neurofilament levels in biofluids are commonly used as a proxy for neurodegeneration in several neurodegenerative disorders. In this study, we aimed to investigate the distribution of neurofilaments in the cerebral cortex of Parkinson's disease (PD), PD with dementia (PDD) and dementia with Lewy bodies (DLB) donors, and its association with pathology load and MRI measures of atrophy and diffusivity.

METHODS

Using a within-subject post-mortem MRI-pathology approach, we included 9 PD, 12 PDD/DLB and 18 age-matched control donors. Cortical thickness and mean diffusivity (MD) metrics were extracted respectively from 3DT1 and DTI at 3T in-situ MRI. After autopsy, pathological hallmarks (pSer129-αSyn, p-tau and amyloid-β load) together with neurofilament light-chain (NfL) and phosphorylated-neurofilament medium- and heavy-chain (p-NfM/H) immunoreactivity were quantified in seven cortical regions, and studied in detail with confocal-laser scanning microscopy. The correlations between MRI and pathological measures were studied using linear mixed models.

RESULTS

Compared to controls, p-NfM/H immunoreactivity was increased in all cortical regions in PD and PDD/DLB, whereas NfL immunoreactivity was increased in the parahippocampal and entorhinal cortex in PDD/DLB. NfL-positive neurons showed degenerative morphological features and axonal fragmentation. The increased p-NfM/H correlated with p-tau load, and NfL correlated with pSer129-αSyn but more strongly with p-tau load in PDD/DLB. Lastly, neurofilament immunoreactivity correlated with cortical thinning in PD and with increased cortical MD in PDD/DLB.

CONCLUSIONS

Taken together, increased neurofilament immunoreactivity suggests underlying axonal injury and neurofilament accumulation in morphologically altered neurons with increased pathological burden. Importantly, we demonstrate that such neurofilament markers at least partly explain MRI measures that are associated with the neurodegenerative process.

Tensor image registration library: Deformable registration of stand-alone histology images to whole-brain post-mortem MRI data

BACKGROUND

Accurate registration between microscopy and MRI data is necessary for validating imaging biomarkers against neuropathology, and to disentangle complex signal dependencies in microstructural MRI. Existing registration methods often rely on serial histological sampling or significant manual input, providing limited scope to work with a large number of stand-alone histology sections. Here we present a customisable pipeline to assist the registration of stand-alone histology sections to whole-brain MRI data.

METHODS

Our pipeline registers stained histology sections to whole-brain post-mortem MRI in 4 stages, with the help of two photographic intermediaries: a block face image (to undistort histology sections) and coronal brain slab photographs (to insert them into MRI space). Each registration stage is implemented as a configurable stand-alone Python script using our novel platform, Tensor Image Registration Library (TIRL), which provides flexibility for wider adaptation. We report our experience of registering 87 PLP-stained histology sections from 14 subjects and perform various experiments to assess the accuracy and robustness of each stage of the pipeline.

RESULTS

All 87 histology sections were successfully registered to MRI. Histology-to-block registration (Stage 1) achieved 0.2-0.4 mm accuracy, better than commonly used existing methods. Block-to-slice matching (Stage 2) showed great robustness in automatically identifying and inserting small tissue blocks into whole brain slices with 0.2 mm accuracy. Simulations demonstrated sub-voxel level accuracy (0.13 mm) of the slice-to-volume registration (Stage 3) algorithm, which was observed in over 200 actual brain slice registrations, compensating 3D slice deformations up to 6.5 mm. Stage 4 combined the previous stages and generated refined pixelwise aligned multi-modal histology-MRI stacks.

CONCLUSIONS

Our open-source pipeline provides robust automation tools for registering stand-alone histology sections to MRI data with sub-voxel level precision, and the underlying framework makes it readily adaptable to a diverse range of microscopy-MRI studies.

Deep active learning for suggestive segmentation of biomedical image stacks via optimisation of Dice scores and traced boundary length

Manual segmentation of stacks of 2D biomedical images (e.g., histology) is a time-consuming task which can be sped up with semi-automated techniques. In this article, we present a suggestive deep active learning framework that seeks to minimise the annotation effort required to achieve a certain level of accuracy when labelling such a stack. The framework suggests, at every iteration, a specific region of interest (ROI) in one of the images for manual delineation. Using a deep segmentation neural network and a mixed cross-entropy loss function, we propose a principled strategy to estimate class probabilities for the whole stack, conditioned on heterogeneous partial segmentations of the 2D images, as well as on weak supervision in the form of image indices that bound each ROI. Using the estimated probabilities, we propose a novel active learning criterion based on predictions for the estimated segmentation performance and delineation effort, measured with average Dice scores and total delineated boundary length, respectively, rather than common surrogates such as entropy. The query strategy suggests the ROI that is expected to maximise the ratio between performance and effort, while considering the adjacency of structures that may have already been labelled - which decrease the length of the boundary to trace. We provide quantitative results on synthetically deformed MRI scans and real histological data, showing that our framework can reduce labelling effort by up to 60-70% without compromising accuracy.

Future research directions to identify risks and mitigation strategies for neurostructural, ocular, and behavioral changes induced by human spaceflight: A NASA-ESA expert group consensus report

A team of experts on the effects of the spaceflight environment on the brain and eye (SANS: Spaceflight-Associated Neuro-ocular Syndrome) was convened by NASA and ESA to (1) review spaceflight-associated structural and functional changes of the human brain and eye, and any interactions between the two; and (2) identify critical future research directions in this area to help characterize the risk and identify possible countermeasures and strategies to mitigate the spaceflight-induced brain and eye alterations. The experts identified 14 critical future research directions that would substantially advance our knowledge of the effects of spending prolonged periods of time in the spaceflight environment on SANS, as well as brain structure and function. They used a paired comparison approach to rank the relative importance of these 14 recommendations, which are discussed in detail in the main report and are summarized briefly below.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

Genetic Specificity of Hippocampal Subfield Volumes, Relative to Hippocampal Formation, Identified in 2148 Young Adult Twins and Siblings

The hippocampus is a complex brain structure with key roles in cognitive and emotional processing and with subregion abnormalities associated with a range of disorders and psychopathologies. Here we combine data from two large independent young adult twin/sibling cohorts to obtain the most accurate estimates to date of genetic covariation between hippocampal subfield volumes and the hippocampus as a single volume. The combined sample included 2148 individuals, comprising 1073 individuals from 627 families (mean age = 22.3 years) from the Queensland Twin IMaging (QTIM) Study, and 1075 individuals from 454 families (mean age = 28.8 years) from the Human Connectome Project (HCP). Hippocampal subfields were segmented using FreeSurfer version 6.0 (CA4 and dentate gyrus were phenotypically and genetically indistinguishable and were summed to a single volume). Multivariate twin modeling was conducted in OpenMx to decompose variance into genetic and environmental sources. Bivariate analyses of hippocampal formation and each subfield volume showed that 10%-72% of subfield genetic variance was independent of the hippocampal formation, with greatest specificity found for the smaller volumes; for example, CA2/3 with 42% of genetic variance being independent of the hippocampus; fissure (63%); fimbria (72%); hippocampus-amygdala transition area (41%); parasubiculum (62%). In terms of genetic influence, whole hippocampal volume is a good proxy for the largest hippocampal subfields, but a poor substitute for the smaller subfields. Additive genetic sources accounted for 49%-77% of total variance for each of the subfields in the combined sample multivariate analysis. In addition, the multivariate analyses were sufficiently powered to identify common environmental influences (replicated in QTIM and HCP for the molecular layer and CA4/dentate gyrus, and accounting for 7%-16% of total variance for 8 of 10 subfields in the combined sample). This provides the clearest indication yet from a twin study that factors such as home environment may influence hippocampal volumes (albeit, with caveats).

Quantitative and histologically validated measures of the entorhinal subfields in ex vivo MRI

Neuroimaging studies have routinely used hippocampal volume as a measure of Alzheimer's disease severity, but hippocampal changes occur too late in the disease process for potential therapies to be effective. The entorhinal cortex is one of the first cortical areas affected by Alzheimer's disease; its neurons are especially vulnerable to neurofibrillary tangles. Entorhinal atrophy also relates to the conversion from non-clinical to clinical Alzheimer's disease. In neuroimaging, the human entorhinal cortex has so far mostly been considered in its entirety or divided into a medial and a lateral region. Cytoarchitectonic differences provide the opportunity for subfield parcellation. We investigated the entorhinal cortex on a subfield-specific level-at a critical time point of Alzheimer's disease progression. While MRI allows multidimensional quantitative measurements, only histology provides enough accuracy to determine subfield boundaries-the pre-requisite for quantitative measurements within the entorhinal cortex. This study used histological data to validate ultra-high-resolution 7 Tesla ex vivo MRI and create entorhinal subfield parcellations in a total of 10 pre-clinical Alzheimer's disease and normal control cases. Using ex vivo MRI, eight entorhinal subfields (olfactory, rostral, medial intermediate, intermediate, lateral rostral, lateral caudal, caudal, and caudal limiting) were characterized for cortical thickness, volume, and pial surface area. Our data indicated no influence of sex, or Braak and Braak staging on volume, cortical thickness, or pial surface area. The volume and pial surface area for mean whole entorhinal cortex were 1131 ± 55.72 mm³ and 429 ± 22.6 mm² (mean ± SEM), respectively. The subfield volume percentages relative to the entire entorhinal cortex were olfactory: 18.73 ± 1.82%, rostral: 14.06 ± 0.63%, lateral rostral: 14.81 ± 1.22%, medial intermediate: 6.72 ± 0.72%, intermediate: 23.36 ± 1.85%, lateral caudal: 5.42 ± 0.33%, caudal: 10.99 ± 1.02%, and caudal limiting: 5.91 ± 0.40% (all mean ± SEM). Olfactory and intermediate subfield revealed the most extensive intra-individual variability (cross-subject variance) in volume and pial surface area. This study provides validated measures. It maps individuality and demonstrates human variability in the entorhinal cortex, providing a baseline for approaches in individualized medicine. Taken together, this study serves as a ground-truth validation study for future in vivo comparisons and treatments.

A semi-automatic registration protocol to match ex-vivo high-field 7T MR images and histological slices in surgical samples from patients with drug-resistant epilepsy

BACKGROUND

MRI is a fundamental tool to detect brain structural anomalies and improvement in this technique has the potential to visualize subtle abnormalities currently undetected. Correlation between pre-operative MRI and histopathology is required to validate the neurobiological basis of MRI abnormalities. However, precise MRI-histology matching is very challenging with the surgical samples. We previously developed a coregistration protocol to match the in-vivo MRI with ex-vivo MRI obtained from surgical specimens. Now, we complete the process to successfully align ex-vivo MRI data with the proper digitalized histological sections in an automatic way.

NEW METHOD

The implemented pipeline is composed by the following steps: a) image pre-processing made of MRI and histology volumes conversion and masking; b) gross rigid body alignment between MRI volume and histology virtual slides; c) rigid alignment between each MRI section and histology slice and estimate of the correlation coefficient for each step to select the MRI slice that best matches histology; d) final linear registration of the selected slices.

RESULTS

This method is fully automatic, except for the first masking step, fast and reliable in comparison to the manual one, as assessed using a Bland-Altman plot.

COMPARISON WITH EXISTING METHODS

The visual assessment usually employed for choosing the best fitting ex-vivo MRI slice for each stained section takes hours and requires practice. Goubran et al. (2015) proposed an iterative registration protocol but its aim and methods were different from ours. No others similar methods are reported in the literature.

CONCLUSIONS

This protocol completes our previous pipeline. The ultimate goal will be to apply the entire process to finely investigate the relationship between clinical MRI data and histopathological features in patients with drug-resistant epilepsy.

Robust joint registration of multiple stains and MRI for multimodal 3D histology reconstruction: Application to the Allen human brain atlas

Joint registration of a stack of 2D histological sections to recover 3D structure ("3D histology reconstruction") finds application in areas such as atlas building and validation of in vivo imaging. Straightforward pairwise registration of neighbouring sections yields smooth reconstructions but has well-known problems such as "banana effect" (straightening of curved structures) and "z-shift" (drift). While these problems can be alleviated with an external, linearly aligned reference (e.g., Magnetic Resonance (MR) images), registration is often inaccurate due to contrast differences and the strong nonlinear distortion of the tissue, including artefacts such as folds and tears. In this paper, we present a probabilistic model of spatial deformation that yields reconstructions for multiple histological stains that that are jointly smooth, robust to outliers, and follow the reference shape. The model relies on a spanning tree of latent transforms connecting all the sections and slices of the reference volume, and assumes that the registration between any pair of images can be see as a noisy version of the composition of (possibly inverted) latent transforms connecting the two images. Bayesian inference is used to compute the most likely latent transforms given a set of pairwise registrations between image pairs within and across modalities. We consider two likelihood models: Gaussian (ℓ2 norm, which can be minimised in closed form) and Laplacian (ℓ1 norm, minimised with linear programming). Results on synthetic deformations on multiple MR modalities, show that our method can accurately and robustly register multiple contrasts even in the presence of outliers. The framework is used for accurate 3D reconstruction of two stains (Nissl and parvalbumin) from the Allen human brain atlas, showing its benefits on real data with severe distortions. Moreover, we also provide the registration of the reconstructed volume to MNI space, bridging the gaps between two of the most widely used atlases in histology and MRI. The 3D reconstructed volumes and atlas registration can be downloaded from https://openneuro.org/datasets/ds003590. The code is freely available at https://github.com/acasamitjana/3dhirest.

Amyloid-β, p-tau and reactive microglia are pathological correlates of MRI cortical atrophy in Alzheimer’s disease

Alzheimer’s disease is characterized by cortical atrophy on MRI and abnormal depositions of amyloid-beta, phosphorylated-tau and inflammation pathologically. However, the relative contribution of these pathological hallmarks to cortical atrophy, a widely used MRI biomarker in Alzheimer’s disease, is yet to be defined. Therefore, the aim of this study was to identify the histopathological correlates of MRI cortical atrophy in Alzheimer’s disease donors, and its typical amnestic and atypical non-amnestic phenotypes. Nineteen Alzheimer’s disease (of which 10 typical and 9 atypical) and 10 non-neurological control brain donors underwent post-mortem in situ 3T 3D-T1, from which cortical thickness was calculated with Freesurfer. Upon subsequent autopsy, 12 cortical brain regions from the right hemisphere and 9 from the left hemisphere were dissected and immunostained for amyloid-beta, phosphorylated-tau and reactive microglia, and percentage area load was calculated for each marker using ImageJ. In addition, post-mortem MRI was compared to ante-mortem MRI of the same Alzheimer’s disease donors when available. MRI-pathology associations were assessed using linear mixed models. Higher amyloid-beta load weakly correlated with higher cortical thickness globally (r = 0.22, P = 0.022). Phosphorylated-tau strongly correlated with cortical atrophy in temporal and frontal regions (−0.76 < r < −1.00, all P < 0.05). Reactive microglia load strongly correlated with cortical atrophy in the parietal region (r = −0.94, P < 0.001). Moreover, post-mortem MRI scans showed high concordance with ante-mortem scans acquired <1 year before death. In conclusion, distinct histopathological markers differently correlated with cortical atrophy, highlighting their different roles in the neurodegenerative process, and therefore contributing to the understanding of the pathological underpinnings of MRI atrophic patterns in Alzheimer’s disease. In our cohort, no or only subtle differences were found in MRI-pathology associations in Alzheimer’s disease phenotypes, indicating that the histopathological correlates of cortical atrophy in typical and atypical phenotypes might be similar. Moreover, we show that post-mortem in situ MRI can be used as proxy for ante-mortem in vivo MRI.

Tutorial: methods for three-dimensional visualization of archival tissue material

Analysis of three-dimensional patient specimens is gaining increasing relevance for understanding the principles of tissue structure as well as the biology and mechanisms underlying disease. New technologies are improving our ability to visualize large volume of tissues with subcellular resolution. One resource often overlooked is archival tissue maintained for decades in hospitals and research archives around the world. Accessing the wealth of information stored within these samples requires the use of appropriate methods. This tutorial introduces the range of sample preparation and microscopy approaches available for three-dimensional visualization of archival tissue. We summarize key aspects of the relevant techniques and common issues encountered when using archival tissue, including registration and antibody penetration. We also discuss analysis pipelines required to process, visualize and analyze the data and criteria to guide decision-making. The methods outlined in this tutorial provide an important and sustainable avenue for validating three-dimensional tissue organization and mechanisms of disease.

Constructing the rodent stereotaxic brain atlas: a survey

The stereotaxic brain atlas is a fundamental reference tool commonly used in the field of neuroscience. Here we provide a brief history of brain atlas development and clarify three key conceptual elements of stereotaxic brain atlasing: brain image, atlas, and stereotaxis. We also refine four technical indices for evaluating the construction of atlases: the quality of staining and labeling, the granularity of delineation, spatial resolution, and the precision of spatial location and orientation. Additionally, we discuss state-of-the-art technologies and their trends in the fields of image acquisition, stereotaxic coordinate construction, image processing, anatomical structure recognition, and publishing: the procedures of brain atlas illustration. We believe that the use of single-cell resolution and micron-level location precision will become a future trend in the study of the stereotaxic brain atlas, which will greatly benefit the development of neuroscience.

Accurate and Robust Alignment of Differently Stained Histologic Images Based on Greedy Diffeomorphic Registration

Histopathologic assessment routinely provides rich microscopic information about tissue structure and disease process. However, the sections used are very thin, and essentially capture only 2D representations of a certain tissue sample. Accurate and robust alignment of sequentially cut 2D slices should contribute to more comprehensive assessment accounting for surrounding 3D information. Towards this end, we here propose a two-step diffeomorphic registration approach that aligns differently stained histology slides to each other, starting with an initial affine step followed by estimating a deformation field. It was quantitatively evaluated on ample (n = 481) and diverse data from the automatic non-rigid histological image registration challenge, where it was awarded the second rank. The obtained results demonstrate the ability of the proposed approach to robustly (average robustness = 0.9898) and accurately (average relative target registration error = 0.2%) align differently stained histology slices of various anatomical sites while maintaining reasonable computational efficiency (<1 min per registration). The method was developed by adapting a general-purpose registration algorithm designed for 3D radiographic scans and achieved consistently accurate results for aligning high-resolution 2D histologic images. Accurate alignment of histologic images can contribute to a better understanding of the spatial arrangement and growth patterns of cells, vessels, matrix, nerves, and immune cell interactions.

Imaging of the human subthalamic nucleus

The human subthalamic nucleus (STN) is a small lens shaped iron rich nucleus, which has gained substantial interest as a target for deep brain stimulation surgery for a variety of movement disorders. The internal anatomy of the human STN has not been fully elucidated, and an intensive debate, discussing the level of overlap between putative limbic, associative, and motor zones within the STN is still ongoing. In this chapter, we have summarized anatomical information obtained using different neuroimaging modalities focusing on the anatomy of the STN. Additionally, we have highlighted a number of major challenges faced when using magnetic resonance imaging (MRI) approaches for the visualization of small iron rich deep brain structures such as the STN. In vivo MRI and postmortem microscopy efforts provide valuable complementary information on the internal structure of the STN, although the results are not always fully aligned. Finally, we provide an outlook on future efforts that could contribute to the development of an integrative research approach that will help with the reconciliation of seemingly divergent results across research approaches.

Mesoscale diffusion magnetic resonance imaging of the ex vivo human hippocampus

Mesoscale diffusion magnetic resonance imaging (MRI) endeavors to bridge the gap between macroscopic white matter tractography and microscopic studies investigating the cytoarchitecture of human brain tissue. To ensure a robust measurement of diffusion at the mesoscale, acquisition parameters were arrayed to investigate their effects on scalar indices (mean, radial, axial diffusivity, and fractional anisotropy) and streamlines (i.e., graphical representation of axonal tracts) in hippocampal layers. A mesoscale resolution afforded segementation of the pyramidal cell layer (CA1-4), the dentate gyrus, as well as stratum moleculare, radiatum, and oriens. Using ex vivo samples, surgically excised from patients with intractable epilepsy (n = 3), we found that shorter diffusion times (23.7 ms) with a b-value of 4,000 s/mm2 were advantageous at the mesoscale, providing a compromise between mean diffusivity and fractional anisotropy measurements. Spatial resolution and sample orientation exerted a major effect on tractography, whereas the number of diffusion gradient encoding directions minimally affected scalar indices and streamline density. A sample temperature of 15°C provided a compromise between increasing signal-to-noise ratio and increasing the diffusion properties of the tissue. Optimization of the acquisition afforded a system's view of intra- and extra-hippocampal connections. Tractography reflected histological boundaries of hippocampal layers. Individual layer connectivity was visualized, as well as streamlines emanating from individual sub-fields. The perforant path, subiculum and angular bundle demonstrated extra-hippocampal connections. Histology of the samples confirmed individual cell layers corresponding to ROIs defined on MR images. We anticipate that this ex vivo mesoscale imaging will yield novel insights into human hippocampal connectivity.

A multimodal computational pipeline for 3D histology of the human brain

Ex vivo imaging enables analysis of the human brain at a level of detail that is not possible in vivo with MRI. In particular, histology can be used to study brain tissue at the microscopic level, using a wide array of different stains that highlight different microanatomical features. Complementing MRI with histology has important applications in ex vivo atlas building and in modeling the link between microstructure and macroscopic MR signal. However, histology requires sectioning tissue, hence distorting its 3D structure, particularly in larger human samples. Here, we present an open-source computational pipeline to produce 3D consistent histology reconstructions of the human brain. The pipeline relies on a volumetric MRI scan that serves as undistorted reference, and on an intermediate imaging modality (blockface photography) that bridges the gap between MRI and histology. We present results on 3D histology reconstruction of whole human hemispheres from two donors.

Ex vivo mesoscopic diffusion MRI correlates with seizure frequency in patients with uncontrolled mesial temporal lobe epilepsy

The role of hippocampal connectivity in mesial temporal lobe epilepsy (mTLE) remains poorly understood. The use of ex vivo hippocampal samples excised from patients with mTLE affords mesoscale diffusion magnetic resonance imaging (MRI) to identify individual cell layers, such as the pyramidal (PCL) and granule cell layers (GCL), which are thought to be impacted by seizure activity. Diffusion tensor imaging (DTI) of control (n = 3) and mTLE (n = 7) hippocampi on an 11.7 T MRI scanner allowed us to reveal intra‐hippocampal connectivity and evaluate how epilepsy affected mean (MD), axial (AD), and radial diffusivity (RD), as well as fractional anisotropy (FA). Regional measurements indicated a volume loss in the PCL of the cornu ammonis (CA) 1 subfield in mTLE patients compared to controls, which provided anatomical context. Diffusion measurements, as well as streamline density, were generally higher in mTLE patients compared to controls, potentially reflecting differences due to tissue fixation. mTLE measurements were more variable than controls. This variability was associated with disease severity, as indicated by a strong correlation (r = 0.87) between FA in the stratum radiatum and the frequency of seizures in patients. MD and RD of the PCL in subfields CA3 and CA4 also correlated strongly with disease severity. No correlation of MR measures with disease duration was evident. These results reveal the potential of mesoscale diffusion MRI to examine layer‐specific diffusion changes and connectivity to determine how these relate to clinical measures. Improving the visualization of intra‐hippocampal connectivity will advance the development of novel hypotheses about seizure networks.

Multimodal cross-registration and quantification of metric distortions in marmoset whole brain histology using diffeomorphic mappings

Whole brain neuroanatomy using tera-voxel light-microscopic data sets is of much current interest. A fundamental problem in this field is the mapping of individual brain data sets to a reference space. Previous work has not rigorously quantified in-vivo to ex-vivo distortions in brain geometry from tissue processing. Further, existing approaches focus on registering unimodal volumetric data; however, given the increasing interest in the marmoset model for neuroscience research and the importance of addressing individual brain architecture variations, new algorithms are necessary to cross-register multimodal data sets including MRIs and multiple histological series. Here we present a computational approach for same-subject multimodal MRI-guided reconstruction of a series of consecutive histological sections, jointly with diffeomorphic mapping to a reference atlas. We quantify the scale change during different stages of brain histological processing using the Jacobian determinant of the diffeomorphic transformations involved. By mapping the final image stacks to the ex-vivo post-fixation MRI, we show that (a) tape-transfer assisted histological sections can be reassembled accurately into 3D volumes with a local scale change of 2.0 ± 0.4% per axis dimension; in contrast, (b) tissue perfusion/fixation as assessed by mapping the in-vivo MRIs to the ex-vivo post fixation MRIs shows a larger median absolute scale change of 6.9 ± 2.1% per axis dimension. This is the first systematic quantification of local metric distortions associated with whole-brain histological processing, and we expect that the results will generalize to other species. These local scale changes will be important for computing local properties to create reference brain maps.

Dentate gyrus volume deficit in schizophrenia

BACKGROUND

Schizophrenia is associated with robust hippocampal volume deficits but subregion volume deficits, their associations with cognition, and contributing genes remain to be determined.

METHODS

Hippocampal formation (HF) subregion volumes were obtained using FreeSurfer 6.0 from individuals with schizophrenia (n = 176, mean age ± s.d. = 39.0 ± 11.5, 132 males) and healthy volunteers (n = 173, mean age ± s.d. = 37.6 ± 11.3, 123 males) with similar mean age, gender, handedness, and race distributions. Relationships between the HF subregion volume with the largest between group difference, neuropsychological performance, and single-nucleotide polymorphisms were assessed.

RESULTS

This study found a significant group by region interaction on hippocampal subregion volumes. Compared to healthy volunteers, individuals with schizophrenia had significantly smaller dentate gyrus (DG) (Cohen's d = -0.57), Cornu Ammonis (CA) 4, molecular layer of the hippocampus, hippocampal tail, and CA 1 volumes, when statistically controlling for intracranial volume; DG (d = -0.43) and CA 4 volumes remained significantly smaller when statistically controlling for mean hippocampal volume. DG volume showed the largest between group difference and significant positive associations with visual memory and speed of processing in the overall sample. Genome-wide association analysis with DG volume as the quantitative phenotype identified rs56055643 (β = 10.8, p < 5 × 10-8, 95% CI 7.0-14.5) on chromosome 3 in high linkage disequilibrium with MOBP. Gene-based analyses identified associations between SLC25A38 and RPSA and DG volume.

CONCLUSIONS

This study suggests that DG dysfunction is fundamentally involved in schizophrenia pathophysiology, that it may contribute to cognitive abnormalities in schizophrenia, and that underlying biological mechanisms may involve contributions from MOBP, SLC25A38, and RPSA.

MRI-based classification of the anatomical variants of the hippocampal head

PURPOSE

The head of the hippocampus (H) is classically described as having two to four digitations both in ex vivo specimens and in vivo MR coronal images. The aim of this study was to develop and evaluate a new MR-based classification of the anatomical variants of the hippocampal head in a large sample population of healthy subjects.

METHODS

MR images of the brain of 238 young healthy subjects (138 men and 100 women; age range 18-39) were analyzed. The head of the H was identified on coronal reformatted 3D T1 weighted MR images. The frequencies were reported for hemisphere and sex. Inter-rater reliability was assessed.

RESULTS

Eight variants of the hippocampal head were described. Class 0 (11.4%) indicated a total absence of sulci. This class was further subdivided as follows: 0A (one digitation, 10.1%) and 0B (no digitations or "null variant", 1.3%). Class 1 (25.6%) presented a single sulcus and was further subdivided into four types according to the location and the width of the sulcus [1A (8.8%), 1B (12.8%), 1C (1.3%), and 1D (2.7%)]. Class 2 (63.0%, the most frequent and the classical variant) had two symmetrical sulci and three digitations. Statistically significant differences between the two hemispheres were observed only in women and overall. Differences in prevalence between sexes were not observed.

CONCLUSIONS

The large study population allowed the description of a novel morphological classification of the different anatomical variants of normal H in the coronal plane. This classification could reduce the risk of misinterpreting normal anatomical variants as pathological.

Diffeomorphic Registration With Intensity Transformation and Missing Data: Application to 3D Digital Pathology of Alzheimer's Disease

This paper examines the problem of diffeomorphic image registration in the presence of differing image intensity profiles and sparsely sampled, missing, or damaged tissue. Our motivation comes from the problem of aligning 3D brain MRI with 100-micron isotropic resolution to histology sections at 1 × 1 × 1,000-micron resolution with multiple varying stains. We pose registration as a penalized Bayesian estimation, exploiting statistical models of image formation where the target images are modeled as sparse and noisy observations of the atlas. In this injective setting, there is no assumption of symmetry between atlas and target. Cross-modality image matching is achieved by jointly estimating polynomial transformations of the atlas intensity. Missing data is accommodated via a multiple atlas selection procedure where several atlas images may be of homogeneous intensity and correspond to "background" or "artifact." The two concepts are combined within an Expectation-Maximization algorithm, where atlas selection posteriors and deformation parameters are updated iteratively and polynomial coefficients are computed in closed form. We validate our method with simulated images, examples from neuropathology, and a standard benchmarking dataset. Finally, we apply it to reconstructing digital pathology and MRI in standard atlas coordinates. By using a standard convolutional neural network to detect tau tangles in histology slices, this registration method enabled us to quantify the 3D density distribution of tauopathy throughout the medial temporal lobe of an Alzheimer's disease postmortem specimen.

Multimodal Parcellations and Extensive Behavioral Profiling Tackling the Hippocampus Gradient

The hippocampus displays a complex organization and function that is perturbed in many neuropathologies. Histological work revealed a complex arrangement of subfields along the medial-lateral and the ventral-dorsal dimension, which contrasts with the anterior-posterior functional differentiation. The variety of maps has raised the need for an integrative multimodal view. We applied connectivity-based parcellation to 1) intrinsic connectivity 2) task-based connectivity, and 3) structural covariance, as complementary windows into structural and functional differentiation of the hippocampus. Strikingly, while functional properties (i.e., intrinsic and task-based) revealed similar partitions dominated by an anterior-posterior organization, structural covariance exhibited a hybrid pattern reflecting both functional and cytoarchitectonic subdivision. Capitalizing on the consistency of functional parcellations, we defined robust functional maps at different levels of partitions, which are openly available for the scientific community. Our functional maps demonstrated a head-body and tail partition, subdivided along the anterior-posterior and medial-lateral axis. Behavioral profiling of these fine partitions based on activation data indicated an emotion-cognition gradient along the anterior-posterior axis and additionally suggested a self-world-centric gradient supporting the role of the hippocampus in the construction of abstract representations for spatial navigation and episodic memory.

Characterization of hippocampal subfields using ex vivo MRI and histology data: Lessons for in vivo segmentation

Hippocampal subfield segmentation on in vivo MRI is of great interest for cognition, aging, and disease research. Extant subfield segmentation protocols have been based on neuroanatomical references, but these references often give limited information on anatomical variability. Moreover, there is generally a mismatch between the orientation of the histological sections and the often anisotropic coronal sections on in vivo MRI. To address these issues, we provide a detailed description of hippocampal anatomy using a postmortem dataset containing nine specimens of subjects with and without dementia, which underwent a 9.4 T MRI and histological processing. Postmortem MRI matched the typical orientation of in vivo images and segmentations were generated in MRI space, based on the registered annotated histological sections. We focus on the following topics: the order of appearance of subfields, the location of subfields relative to macroanatomical features, the location of subfields in the uncus and tail and the composition of the dark band, a hypointense layer visible in T2-weighted MRI. Our main findings are that: (a) there is a consistent order of appearance of subfields in the hippocampal head, (b) the composition of subfields is not consistent in the anterior uncus, but more consistent in the posterior uncus, (c) the dark band consists only of the CA-stratum lacunosum moleculare, not the strata moleculare of the dentate gyrus, (d) the subiculum/CA1 border is located at the middle of the width of the hippocampus in the body in coronal plane, but moves in a medial direction from anterior to posterior, and (e) the variable location and composition of subfields in the hippocampal tail can be brought back to a body-like appearance when reslicing the MRI scan following the curvature of the tail. Our findings and this publicly available dataset will hopefully improve anatomical accuracy of future hippocampal subfield segmentation protocols.

Dissociable age and memory relationships with hippocampal subfield volumes in vivo:Data from the Irish Longitudinal Study on Ageing (TILDA)

The heterogeneous specialisation of hippocampal subfields across memory functions has been widely shown in animal models. Yet, few in vivo studies in humans have explored correspondence between hippocampal subfield anatomy and memory performance in ageing. Here, we used a well-validated automated MR segmentation protocol to measure hippocampal subfield volumes in 436 non-demented adults aged 50+. We explored relationships between hippocampal subfield volume and verbal episodic memory, as indexed by word list recall at immediate presentation and following delay. In separate multilevel models for each task, we tested linearity and non-linearity of associations between recall performance and subfield volume. Fully-adjusted models revealed that immediate and delayed recall were both associated with cubic fits with respect to volume of subfields CA1, CA2/3, CA4, molecular layer, and granule cell layer of dentate gyrus; moreover, these effects were partly dissociable from quadratic age trends, observed for subiculum, molecular layer, hippocampal tail, and CA1. Furthermore, analyses of semantic fluency data revealed little evidence of robust associations with hippocampal subfield volumes. Our results show that specific hippocampal subfields manifest associations with memory encoding and retrieval performance in non-demented older adults; these effects are partly dissociable from age-related atrophy, and from retrieval of well-consolidated semantic categories.

Matching ex vivo MRI With Iron Histology: Pearls and Pitfalls

Iron levels in the brain can be estimated using newly developed specific magnetic resonance imaging (MRI) sequences. This technique has several applications, especially in neurodegenerative disorders like Alzheimer's disease or Parkinson's disease. Coupling ex vivo MRI with histology allows neuroscientists to better understand what they see in the images. Iron is one of the most extensively studied elements, both by MRI and using histological or physical techniques. Researchers were initially only able to make visual comparisons between MRI images and different types of iron staining, but the emergence of specific MRI sequences like R2* or quantitative susceptibility mapping meant that quantification became possible, requiring correlations with physical techniques. Today, with advances in MRI and image post-processing, it is possible to look for MRI/histology correlations by matching the two sorts of images. For the result to be acceptable, the choice of methodology is crucial, as there are hidden pitfalls every step of the way. In order to review the advantages and limitations of ex vivo MRI correlation with iron-based histology, we reviewed all the relevant articles dealing with the topic in humans. We provide separate assessments of qualitative and quantitative studies, and after summarizing the significant results, we emphasize all the pitfalls that may be encountered.

Functional connectivity along the anterior-posterior axis of hippocampal subfields in the ageing human brain

While age‐related volumetric changes in human hippocampal subfields have been reported, little is known about patterns of subfield functional connectivity (FC) in the context of healthy ageing. Here we investigated age‐related changes in patterns of FC down the anterior–posterior axis of each subfield. Using high resolution structural MRI we delineated the dentate gyrus (DG), CA fields (including separating DG from CA3), the subiculum, pre/parasubiculum, and the uncus in healthy young and older adults. We then used high resolution resting state functional MRI to measure FC in each group and to directly compare them. We first examined the FC of each subfield in its entirety, in terms of FC with other subfields and with neighboring cortical regions, namely, entorhinal, perirhinal, posterior parahippocampal, and retrosplenial cortices. Next, we analyzed subfield to subfield FC within different portions along the hippocampal anterior–posterior axis, and FC of each subfield portion with the neighboring cortical regions of interest. In general, the FC of the older adults was similar to that observed in the younger adults. We found that, as in the young group, the older group displayed intrinsic FC between the subfields that aligned with the tri‐synaptic circuit but also extended beyond it, and that FC between the subfields and neighboring cortical areas differed markedly along the anterior–posterior axis of each subfield. We observed only one significant difference between the young and older groups. Compared to the young group, the older participants had significantly reduced FC between the anterior CA1‐subiculum transition region and the transentorhinal cortex, two brain regions known to be disproportionately affected during the early stages of age‐related tau accumulation. Overall, these results contribute to ongoing efforts to characterize human hippocampal subfield connectivity, with implications for understanding hippocampal function and its modulation in the ageing brain.

Cytoarchitectonically-driven MRI atlas of nonhuman primate hippocampus: Preservation of subfield volumes in aging

Identification of primate hippocampal subfields in vivo using structural MRI imaging relies on variable anatomical guidelines, signal intensity differences, and heuristics to differentiate between regions (Yushkevich et al., 2015a). Thus, a clear anatomically-driven basis for subfield demarcation is lacking. Recent work, however, has begun to develop methods to use ex vivo histology or ex vivo MRI (Adler et al., 2014; Iglesias et al., 2015) that have the potential to inform subfield demarcations of in vivo images. For optimal results, however, ex vivo and in vivo images should ideally be matched within the same healthy brains, with the goal to develop a neuroanatomically-driven basis for in vivo structural MRI images. Here, we address this issue in young and aging rhesus macaques (young n = 5 and old n = 5) using ex vivo Nissl-stained sections in which we identified the dentate gyrus, CA3, CA2, CA1, subiculum, presubiculum, and parasubiculum guided by morphological cell properties (30 μm thick sections spaced at 240 μm intervals and imaged at 161 nm/pixel). The histologically identified boundaries were merged with in vivo structural MRIs (0.625 × 0.625 × 1 mm) from the same subjects via iterative rigid and diffeomorphic registration resulting in probabilistic atlases of young and old rhesus macaques. Our results indicate stability in hippocampal subfield volumes over an age range of 13 to 32 years, consistent with previous results showing preserved whole hippocampal volume in aged macaques (Shamy et al., 2006). Together, our methods provide a novel approach for identifying hippocampal subfields in non-human primates and a potential 'ground truth' for more accurate identification of hippocampal subfield boundaries on in vivo MRIs. This could, in turn, have applications in humans where accurately identifying hippocampal subfields in vivo is a critical research goal.

Volumetric comparison of hippocampal subfields extracted from 4-minute accelerated vs. 8-minute high-resolution T2-weighted 3T MRI scans

The hippocampus has been widely studied using neuroimaging, as it plays an important role in memory and learning. However, hippocampal subfield information is difficult to capture by standard magnetic resonance imaging (MRI) techniques. To facilitate morphometric study of hippocampal subfields, ADNI introduced a high resolution (0.4 mm in plane) T2-weighted turbo spin-echo sequence that requires 8 min. With acceleration, the protocol can be acquired in 4 min. We performed a comparative study of hippocampal subfield volumes using standard and accelerated protocols on a Siemens Prisma 3T MRI in an independent sample of older adults that included 10 cognitively normal controls, 9 individuals with subjective cognitive decline, 10 with mild cognitive impairment, and 6 with a clinical diagnosis of Alzheimer's disease (AD). The Automatic Segmentation of Hippocampal Subfields (ASHS) software was used to segment 9 primary labeled regions including hippocampal subfields and neighboring cortical regions. Intraclass correlation coefficients were computed for reliability tests between 4 and 8 min scans within and across the four groups. Pairwise group analyses were performed, covaried for age, sex and total intracranial volume, to determine whether the patterns of group differences were similar using 4 vs. 8 min scans. The 4 and 8 min protocols, analyzed by ASHS segmentation, yielded similar volumetric estimates for hippocampal subfields as well as comparable patterns of differences between study groups. The accelerated protocol can provide reliable imaging data for investigation of hippocampal subfields in AD-related MRI studies and the decreased scan time may result in less vulnerability to motion.

On variational solutions for whole brain serial-section histology using a Sobolev prior in the computational anatomy random orbit model

This paper presents a variational framework for dense diffeomorphic atlas-mapping onto high-throughput histology stacks at the 20 μm meso-scale. The observed sections are modelled as Gaussian random fields conditioned on a sequence of unknown section by section rigid motions and unknown diffeomorphic transformation of a three-dimensional atlas. To regularize over the high-dimensionality of our parameter space (which is a product space of the rigid motion dimensions and the diffeomorphism dimensions), the histology stacks are modelled as arising from a first order Sobolev space smoothness prior. We show that the joint maximum a-posteriori, penalized-likelihood estimator of our high dimensional parameter space emerges as a joint optimization interleaving rigid motion estimation for histology restacking and large deformation diffeomorphic metric mapping to atlas coordinates. We show that joint optimization in this parameter space solves the classical curvature non-identifiability of the histology stacking problem. The algorithms are demonstrated on a collection of whole-brain histological image stacks from the Mouse Brain Architecture Project.

New developments in neural tracing techniques have motivated the widespread use of histology as a modality for exploring the circuitry of the brain. Automated mapping of pre-labeled atlases onto modern large datasets of histological imagery is a critical step for elucidating the brain’s neural circuitry and shape. This task is challenging as histological sections are imaged independently and the reconstruction of the unsectioned volume is nontrivial. Typically, neuroanatomists use reference volumes of the same subject (e.g. MRI) to guide reconstruction. However, obtaining reference imagery is often non-standard, as in high-throughput animal models like mouse histology. Others have proposed using anatomical atlases as guides, but have not accounted for the intrinsic nonlinear shape difference from atlas to subject. Our method addresses these limitations by jointly optimizing reconstruction informed by an atlas simultaneously with the nonlinear change of coordinates that encapsulates anatomical variation. This accounts for intrinsic shape differences and enables rigorous, direct comparisons of atlas and subject coordinates. Using simulations, we demonstrate that our method recovers the reconstruction parameters more accurately than atlas-free models and innately produces accurate segmentations from simultaneous atlas mapping. We also demonstrate our method on the Mouse Brain Architecture dataset, successfully mapping and reconstructing over 1000 brains.

Mapping Histological Slice Sequences to the Allen Mouse Brain Atlas Without 3D Reconstruction

Histological brain slices are widely used in neuroscience to study the anatomical organization of neural circuits. Systematic and accurate comparisons of anatomical data from multiple brains, especially from different studies, can benefit tremendously from registering histological slices onto a common reference atlas. Most existing methods rely on an initial reconstruction of the volume before registering it to a reference atlas. Because these slices are prone to distortions during the sectioning process and often sectioned with non-standard angles, reconstruction is challenging and often inaccurate. Here we describe a framework that maps each slice to its corresponding plane in the Allen Mouse Brain Atlas (2015) to build a plane-wise mapping and then perform 2D nonrigid registration to build a pixel-wise mapping. We use the L2 norm of the histogram of oriented gradients difference of two patches as the similarity metric for both steps and a Markov random field formulation that incorporates tissue coherency to compute the nonrigid registration. To fix significantly distorted regions that are misshaped or much smaller than the control grids, we train a context-aggregation network to segment and warp them to their corresponding regions with thin plate spline. We have shown that our method generates results comparable to an expert neuroscientist and is significantly better than reconstruction-first approaches. Code and sample dataset are available at sites.google.com/view/brain-mapping.