High-resolution In Vivo Manual Segmentation Protocol for Human Hippocampal Subfields Using 3T Magnetic Resonance Imaging

Summary

The goal of this manuscript is to study the hippocampus and hippocampal subfields using MRI. The manuscript describes a protocol for segmenting the hippocampus and five hippocampal substructures: cornu ammonis (CA) 1, CA2/CA3, CA4/dentate gyrus, strata radiatum/lacunosum/moleculare, and subiculum.

Abstract

The human hippocampus has been broadly studied in the context of memory and normal brain function and its role in different neuropsychiatric disorders has been heavily studied. While many imaging studies treat the hippocampus as a single unitary neuroanatomical structure, it is, in fact, composed of several subfields that have a complex three-dimensional geometry. As such, it is known that these subfields perform specialized functions and are differentially affected through the course of different disease states. Magnetic resonance (MR) imaging can be used as a powerful tool to interrogate the morphology of the hippocampus and its subfields. Many groups use advanced imaging software and hardware (>3T) to image the subfields; however this type of technology may not be readily available in most research and clinical imaging centers. To address this need, this manuscript provides a detailed step-by-step protocol for segmenting the full anterior-posterior length of the hippocampus and its subfields: cornu ammonis (CA) 1, CA2/CA3, CA4/dentate gyrus (DG), strata radiatum/lacunosum/moleculare (SR/SL/SM), and subiculum. This protocol has been applied to five subjects (3F, 2M; age 29-57, avg. 37). Protocol reliability is assessed by resegmenting either the right or left hippocampus of each subject and computing the overlap using the Dice's kappa metric. Mean Dice's kappa (range) across the five subjects are: whole hippocampus, 0.91 (0.90-0.92); CA1, 0.78 (0.77-0.79); CA2/CA3, 0.64 (0.56-0.73); CA4/dentate gyrus, 0.83 (0.81-0.85); strata radiatum/lacunosum/moleculare, 0.71 (0.68-0.73); and subiculum 0.75 (0.72-0.78). The segmentation protocol presented here provides other laboratories with a reliable method to study the hippocampus and hippocampal subfields in vivo using commonly available MR tools.

Introduction

The hippocampus is a widely studied medial temporal lobe structure that is associated with episodic memory, spatial navigation, and other cognitive functions^10,31. Its role in neurodegenerative and neuropsychiatric disorders such as Alzheimer’s disease, schizophrenia, and bipolar disorder is well-documented^4,5,18,24,30. The goal of this manuscript is to provide additional detail to the manual segmentation protocol published previously³⁴ for human hippocampal subfields on high-resolution magnetic resonance (MR) images acquired at 3T. Additionally, the video component accompanying this manuscript will provide further assistance to researchers who wish to implement the protocol on their own datasets.

The hippocampus can be divided into subfields based on cytoarchitectonic differences observed in histologically-prepared post-mortem specimens^12,22. Such post-mortem specimens define the ground truth for the identification and study of hippocampal subfields; however preparations of this nature require specialized skills and equipment for staining, and are limited by the availability of fixed tissue, especially in diseased populations. In vivo imaging has the advantage of a much larger pool of subjects, and also presents the opportunity for follow-up studies and observing changes in populations. Although it has been shown that signal intensities in T2-weighted ex vivo MR images reflect cellular density¹³, it is still difficult to identify undisputed borders between subfields using solely MR signal intensities. As such, a number of different approaches for identifying histology-level detail on MR images have been developed.

Some groups have made efforts to reconstruct and digitize histological datasets and then use these reconstructions along with image registration techniques to localize hippocampal subfield neuroanatomy on in vivo MR^{1,2,8,9,14,15,17,32}. Although this is an effective technique for mapping a version of the histological ground truth directly onto MR images, reconstructions of this nature are difficult to complete. Projects such as these are limited by the availability of intact medial temporal lobe specimens, histological techniques, data loss during histological processing, and the fundamental morphological inconsistencies between fixed and in vivo brains. Other groups have used high-field scanners (7T or 9.4T) in an effort to acquire in vivo or ex vivo images with a small enough (0.20-0.35 mm isotropic) voxel size to visualize spatially localized differences in image contrast that are used to infer boundaries between subfields^35,37. Even at 7T-9.4T and with such a small voxel size, the cytoarchitectonic characteristics of hippocampal subfields are not visible. As such, manual segmentation protocols have been developed that approximate the known histological boundaries on MR images. These protocols determine subfield boundaries by interpreting local image contrast differences and defining geometric rules (such as straight lines and angles) relative to visible structures. Although images taken at a high field strength are able to offer detailed insight into hippocampal subfields, high-field scanners are not yet common in clinical or research settings, so 7T and 9.4T protocols currently have limited applicability. Similar protocols have been developed for images collected on 3T and 4T scanners^{11,20,21,23,24,25,28,33}. Many of these protocols are based on images with sub-1mm voxels voxel dimensions in the coronal plane, but have large slice thicknesses (0.8-3 mm)^{11,20,21,23,25,28,33} or large inter-slice distances^20,28, both of which result in a significant measurement bias in the estimation of volumes of the individual subfields. Additionally, many of the existing 3T protocols exclude subfields in all or part of the hippocampal head or tail^20,23,25,33 or do not provide detailed segmentations of important substructures (i.e., combine the DG with CA2/CA3 or do not include the strata radiatum/lacunosum/moleculare of the CA)^{11,20,21,23,24,25,28,33}. There is therefore a need in the field for a detailed description of a protocol that can reliably identify relevant subfields throughout the head, body, and tail of the hippocampus that is based on a scanner commonly available in clinical and research settings. Efforts are currently underway by the Hippocampal Subfields Group (www.hippocampalsubfields.com) to harmonize the hippocampal subfield segmentation process between laboratories, similar to an existing harmonization effort for whole hippocampal segmentation⁶, and an initial paper comparing 21 existing protocols was recently published³⁸. The work from this group will further elucidate optimal segmentation procedures.

This manuscript provides detailed written and video instructions for reliably implementing the hippocampal subfield segmentation protocol described previously by Winterburn and colleagues³⁴ on high-resolution 3T MR images. The protocol has been implemented on five images of healthy controls for the whole hippocampus and five hippocampal subfields (CA1, CA2/CA3, CA4/dentate gyrus, strata radiatum/lacunosum/moleculare, and subiculum). These segmented images are available to the public online (cobralab.ca/atlases/Hippocampus). The protocol and the segmented images will be useful for groups who wish to study detailed hippocampal neuroanatomy in MR images.

Protocol

Study Participants

The protocol in this manuscript was developed for five representative high-resolution images collected from healthy volunteers (3F, 2M; age 29-57, avg. 37) who were free of neurological and neuropsychiatric disorders and cases of severe head trauma. All subjects were recruited at the Centre for Addiction and Mental Health (CAMH). The study was approved by the CAMH Research Ethics Board and was conducted in keeping with the Declaration of Helsinki. All subjects provided written, informed consent for data acquisition and sharing. For details about the acquisition sequence used to collect these images, please refer to Winterburn et al., 2013 and Park et al., 2014.^26,34 Images for all five subjects were checked for quality and retained. The hippocampus spanned an average of 118 coronal slices in these images.

1. Software Set-up

1. Open Display: From the terminal using the following command: Display image_name.mnc -label label_name.mnc. The program will open 3 windows: 3D visualization window, 3-orientation image viewing window, and a navigation window. The terminal will also be used to run the program. Enlarge the coronal view, as the segmentations will be performed coronally. Zoom in on the hippocampus. Select F (Segmenting) in the navigation window. Select F (XY Radius:0.1). The terminal window will prompt for the user to “Enter xy brush size: ”. Set to 0.1. This will set the size of your paintbrush. The user can now begin drawing the hippocampus onto the MR image.

2. Whole Hippocampus Manual Segmentation

1. Set-up: Using a T1-weighted image, scroll to the anterior-most coronal slice of the hippocampus. To advance slices in the anterior direction, use the '+' key; use the '-' key to move in the posterior direction.
2. Slice A: Anterior-Most Slice: Using the right-click on the mouse, draw the outer-most border of the hippocampal grey matter where it meets the surrounding temporal lobe white matter and use the high-intensity white matter of the alveus to assist with the superior border, where the hippocampus meets the amygdala^12,22. Use the E (Label Fill) key in the segmentation menu of the navigation window to fill in the label inside the border. Continue to apply these borders throughout the anterior hippocampal head.
3. Slice B: Hippocampal Head 1 (Figure 1B):
   1. Superior, inferior, lateral, medial borders: Continue to draw the borders as described in step 2.2, using the white matter of the temporal lobe and alveus as a guide.
   2. Supero-medial border: For this, using the axial view, draw a horizontal line from the anterior edge of the lateral hippocampus²⁹, and include anything below this line as hippocampus. NOTE: The supero-medial border becomes more ambiguous in these slices, where the grey matter of the hippocampus blends with the grey matter of the amygdala.
4. Slice C: Hippocampal Head 2 with Dentations: Depending on the subject, the dentations of the hippocampus may be visible for 3-4 slices (typically, they are more visible on T2-weighted versus T1-weighted images). In these slices, continue to use the white matter of the alveus and temporal lobe to guide border segmentation^12,22. For further details, follow steps 2.5.1-2.5.2.
5. Slice D: Hippocampal Head 3:
   1. Superior, inferior, lateral, medial borders: Draw the inferior border of the hippocampus at the white matter of the temporal lobe, the lateral border at the inferior horn of the lateral ventricle, the superior border, following the curve of the dentations, at the white matter of the alveus/fimbria, and the medial border at the hypointense region of the ambient cistern^12,22.
   2. Supero-medial and infero-medial borders: Continue to define the supero-medial border as described in step 2.3.2. Draw the inferior portion of the medial border where the hippocampus thins slightly and extends into the mildly hyperintense grey matter of the entorhinal cortex^12,22.
6. Slice E: Hippocampal Head 4 with Uncus: Continue to draw the inferior, lateral, and superior borders described in steps 2.5.1-2.5.2. Include the uncus (which is located medal to the main body of the hippocampus and is surrounded by low-intensity CSF) in the hippocampal segmentation^12,22.
7. Slice F: Hippocampal Body: Continue to draw the inferior, lateral, medial, and superior borders described in steps 2.5.1-2.5.2. Draw the infero-medial border at the point where the hippocampus thins as it transitions to entorhinal cortex/para-hippocampal gyrus^12,22. Do not include the low-intensity CSF of the vestigial hippocampal sulcus in the segmentation.
8. Slice G: Hippocampal Tail 1: Begin segmenting hippocampal tail-type slices when the crus of the fornix is first visible. Exclude the fascicular gyrus (a grey matter structure which blends with the hippocampus in parts of the hippocampal tail) from the segmentation by extrapolating the shape of the fascicular gyrus into the hippocampal tail from more anterior slices^12,22. This extrapolation is only possible for 2-3 slices, after which the two structures cannot be accurately distinguished; at this point, treat all visible grey matter in this area as hippocampus .
9. Slice H: Hippocampal Tail 2: Segment the low-intensity grey matter of the posterior hippocampal tail from the surrounding high-intensity white matter.
10. Slice I: Posterior-Most Slice: Segment the small remaining area of hippocampal grey matter from the surrounding white matter of the temporal lobe.

3. Hippocampal Subfield Manual Segmentation

1. Set-up: Using a T2-weighted image, scroll to the anterior-most coronal slice of the hippocampus (as in step 2.1). To change the color of the paintbrush, select D (Set Paint Lbl:) on the segmenting menu in the navigation window. The command terminal will prompt: “Enter current paint label:”. Enter a number between 1 and 255. Each number corresponds to a different label color.
2. Slice A: Anterior-Most Slice: Since subfield divisions are not yet visible in the anterior-most slice, draw a line dividing the visible hippocampal grey matter along its longest visible axis (which is not necessarily parallel to any of the cardinal axes) into two equal sections to approximate the true anatomy^12,22. Label the superior of these two sections as CA1 and the inferior section as subiculum by choosing a different colored label for each subfield^23,35.
3. Slice B: Hippocampal Head 1: Label the low-intensity area in the middle of the hippocampal formation as SR/SL/SM^13,37. When the bend along the inferior edge of the hippocampus becomes clear, use this landmark as the lateral border separating the subiculum from the CA1^12,22. Continue to follow the longest axis of the hippocampus to draw the CA1-subiculum border on the supero-medial tip³⁷.
4. Slice C: Hippocampal Head 2 with Dentations:
   1. SR/SL/SM, CA4/DG, and subiculum: Label the SR/SL/SM, CA4/DG, and subiculum as described for slice D (step 3.5.1).
   2. CA2/CA3 and CA1: Define the border between CA1 and CA2/CA3 as a 45º angle line extending in the supero-lateral direction from the most supero-lateral edge of the SR/SL/SM^12,22. Extend the CA2/CA3 medially along the superior edge to the trough between the dentations^12,22. Label the rest of the superior edge as CA1^12,22.
5. Slice D: Hippocampal Head 3
   1. SR/SL/SM, CA4/DG, and subiculum: Label the dark SR/SL/SM band first, which will follow the curve of the CA1³⁷. Label any high-intensity grey matter inside of the SR/SL/SM as CA4/DG^{12,22,23,35,37}. This may not be a continuous region, as in Figure 2C. Continue to define the subiculum-CA1 border using the bend in the inferior hippocampus^12,22.
   2. CA2/CA3 and CA1: Continue to define the CA1 and CA2/CA3 border as in step 3.4.2. Extend the CA2/CA3 medially halfway along the superior edge of the hippocampus^12,22 and label the other half of the superior edge as CA1^12,22.
   3. Supero-medial hippocampal head: In this slice, divide the supero-medial hippocampal head vertically in half. Label the medial half as SR/SL/SM¹². Divide the lateral half in half again, this time horizontally. Label the superior portion as CA4/DG and the inferior portion as CA2/CA3¹².
6. Slice E: Hippocampal Head 4 with Uncus
   1. Lateral hippocampal head (subiculum): In the lateral portion of these slices, define the subiculum-CA1 border as a vertical line extending in the inferior direction from the most medial edge of the CA4/DG^12,22.
   2. Lateral hippocampal head (CA1, CA2/CA3, CA4/DG, SR/SL/SM.): Define the CA1-CA2/CA3 border in the same way as in step 3.4.2. Continue to label the SR/SL/SM as the low intensity region following the curve of the CA regions. Label the CA4/DG as the center cavity inside the SR/SL/SM, as in step 3.5.1.
   3. Uncal hippocampal head (SR/SL/SM): Label the uncus of the hippocampus for approximately 10 slices as the hippocampal head transitions into the hippocampal body. In the uncus, label the low intensity region in the center as SR/SL/SM (when this is difficult to see, approximate the anatomy by segmenting a line 2-3 voxels wide up the center of the uncus)¹².
   4. Uncal hippocampal head (CA2/CA3, CA4/DG): Draw a line at the superior edge of the SR/SL/SM section along infero-lateral/supero-medial axis of the uncus. Label all grey matter above this line as CA2/CA3¹². Label any unlabeled grey matter below this line (on either side of the SR/SL/SM) as CA4/DG¹².
7. Slice F: Hippocampal Body: Continue to apply the borders described in step 3.6.1-3.6.2.
8. Slice G: Hippocampal Tail 1: Continue to apply the rules described in step 3.6.1-3.6.2. The subiculum-CA1 border becomes a 45º angle line extending in the infero-medial direction from the medial edge of the CA4/DG^12,22.
9. Slice H: Hippocampal Tail 2: Once the fascicular gyrus can no longer be distinguished from the hippocampal formation, label the entire outer layer as CA1, the low-intensity area inside of this as SR/SL/SM (as in previous slices), and any remaining grey matter in the middle as CA4/DG^12,22.
10. Slice I: Posterior-Most Slice: Once the dark SR/SL/SM is no longer visible in the center of the hippocampal formation, label the entire structure as CA1^12,22.

4. Protocol Reliability

1. Resegment either the right or left hippocampus of each subject after waiting approximately one month from performing the original segmentation. Segment all of the subfields along the entire anterior-posterior length of the hippocampus, trying to follow the protocol rules as consistently as is possible.
2. Calculate the Dice’s kappa between the original and resegmented volumes:
   
   where k=Dice’s kappa and A and B are label volumes.

Results

Results from the protocol reliability test are summarized in Table 2. For the whole bilateral hippocampus, mean spatial overlap as measured by Dice’s kappa is 0.91 and ranges from 0.90 - 0.92. Subfield kappa values range from 0.64 (CA2/CA3) to 0.83 (CA4/dentate gyrus). Mean volumes for all subfields and the whole hippocampus are reported in Table 3. Volumes for the whole hippocampus range from 2456.72-3325.02 mm³. The CA2/CA3 is the smallest subfield at 208.33 mm³

Discussion

Hippocampal subfield segmentation in MR images is well-represented in the literature. However, existing protocols exclude portions of the hippocampus^20,23,33,35, apply only to fixed images³⁷, or require ultra-high field scanners for image acquisition^35,37. This manuscript offers a segmentation protocol that includes five major subdivisions (CA1, CA2/CA3, CA4/dentate gyrus, SR/SL/SM, and subiculum) of the hippocampus and spans the entire anterior-posterior length of the structure. The comp...

Materials

Name	Company	Catalog Number	Comments
Discovery MR750 3T	GE		Or equivalent 3T scanner
Minc Tool Kit	McConnell Brain Imaging Center, Montreal Neurological Institute		Open source: http://www.bic.mni.mcgill.ca/ServicesSoftware/ServicesSoftwareMincToolKit