Homologous laminar organization of the mouse and human subiculum

The subiculum is the major output component of the hippocampal formation and one of the major brain structures most affected by Alzheimer’s disease. Our previous work revealed a hidden laminar architecture within the mouse subiculum. However, the rotation of the hippocampal longitudinal axis across species makes it unclear how the laminar organization is represented in human subiculum. Using in situ hybridization data from the Allen Human Brain Atlas, we demonstrate that the human subiculum also contains complementary laminar gene expression patterns similar to the mouse. In addition, we provide evidence that the molecular domain boundaries in human subiculum correspond to microstructural differences observed in high resolution MRI and fiber density imaging. Finally, we show both similarities and differences in the gene expression profile of subiculum pyramidal cells within homologous lamina. Overall, we present a new 3D model of the anatomical organization of human subiculum and its evolution from the mouse.

www.nature.com/scientificreports/ humans this axis is rotated into the posterior-anterior direction. Based on this rotation, the mouse dorsal SUB is generally believed to be homologous to the human posterior SUB and mouse ventral SUB is homologous to anterior human SUB. Functional evidence supports this view as the mouse dorsal SUB and human posterior SUB are involved in visuospatial navigation 10,[17][18][19] . In contrast, the mouse ventral SUB and human anterior SUB are related to limbic emotional processing and social behaviors 12,20,21 . Based on the anatomical and functional homology, we hypothesized that the laminar gene expression patterns that delineate mouse HGEA SUB subregions would demonstrate a similar relationship pattern in the corresponding parts of human SUB. Analyzing gene expression patterns from the Allen Human Brain Atlas in situ hybridization database, we identified genes with unique spatial distribution patterns that characterize relatively distinct and complementary lamina across the human posterior and anterior SUB. Particularly, the posterior human SUB contains three complementary gene expression layers whose distribution patterns strongly reflect the mouse SUBdd and ProSUB when accounting for the rotation of the longitudinal axis. Additionally, posterior SUB gene expression boundaries aligned with differences in fiber density microstructure demonstrated by ex vivo diffusion-weighted magnetic resonance imaging (DW-MRI). In the anterior SUB, we found these lamina continued within ventral parts of the SUB, but not in the dorsal parts of the anterior SUB, suggesting a difference between dorsal and ventral parts of the anterior SUB that was not previously known. Taking this new data together with our previous understanding of the mouse subiculum, we propose a new 3D model of the human subiculum and its homology to the rodent. Overall, our gene expression analysis provides a new understanding of the homology between mouse and human hippocampus that is critical to translational studies of hippocampal diseases such as Alzheimer's disease, hippocampal sclerosis, and epilepsy.

Results
All in situ hybridization data was downloaded from the online Allen Brain Atlas image database (www.brain -map. org). Compared to the mouse database which includes both coronal and sagittal tissue sections, the current Allen Human Brain Atlas database contains only a limited set of in situ hybridization data in coronally-sectioned tissue (www.human .brain -map.org). Two separate tissue blocks containing hippocampus/amygdala were found in four subjects as part of the Neurotransmitter Study: one block of tissue includes the anterior hippocampal pole and another block of tissue at a more intermediate/posterior hippocampal level. We observed relatively consistent gene expression patterns in all 4 cases although because of variability in tissue dissection and sectioning quality, it was difficult to relate corresponding sections across subjects. Therefore, we present data from the tissue series with the best histological quality (H0351.1010, 28 year old Hispanic male with no known cognitive impairment) to demonstrate their patterns across adjacent rostrocaudal sections (tissue index (TI) identifies adjacent series section number). We will first describe observed gene expression patterns in the human posterior SUB followed (yellow), CA1 (red), and SUB (green); see also Supplementary Movie 1 or use Schol-AR app). The longitudinal hippocampal axis (red in axes chart) in mouse is oriented dorsoventrally, whereas the human longitudinal axis is rotated into the anterior-posterior axis. In addition the anterior/posterior (septo-temporal; blue color) axis in mice is oriented in the superior-inferior direction in humans. (b) Sagittal view of human brain volume representing the spatial location of all in situ hybridization datasets from the Allen Human Brian Atlas (top). Two tissue blocks containing posterior (middle) and anterior (bottom) parts of the hippocampus (white arrows) are shown overlaid on sagittal 3 T structural MRI images. All images in (b) are from www.brain -map.org. Data viewable with Schol-AR augmented reality app, for details visit https ://www.ini.usc.edu/schol ar/downl oad.html.
by analysis of anterior SUB gene expression patterns. In addition, we provide evidence from ex vivo MRI imaging data that gene expression boundary delineations correspond to differences in human imaging microstructure. Finally, we provide a comparison of mouse and human gene expression within homologous SUB lamina.
Gene expression patterns in the posterior human SUB. Within posterior hippocampal tissue sections, in situ hybridization gene expression patterns reveal the laminar organization of SUB pyramidal neurons. Notably, neurons expressing Nts, Chrm2, and Htr2a are broadly distributed as three complementary gene expression patterns across the CA1, ProSUB, and SUB (Fig. 2). Nts is robustly expressed within a layer of superficial SUB pyramidal neurons (adjacent to the molecular layer), whereas Chrm2-expressing pyramidal neurons are distributed as a deep thin layer directly adjacent to the alveus white matter tract. The deep Chrm2 gene expression layer is thicker laterally where Chrm2 expression also continues into the presubiculum (PRE), but gradually becomes thinner as it extends medially into the ProSUB region. In contrast to the SUB, the superficial ProSUB neurons near the molecular layer robustly express Htr2a rather than Nts. Htr2a expression is continuous within CA1 and ProSUB, but ends near the SUB border. Together, Nts, Chrm2, and Htr2a expression patterns identify three separate gene expression domains across the human SUB and ProSUB with relatively distinct boundaries. The SUB and ProSUB can be distinguished by Nts vs. Htr2a expression although each region contains a common deep layer of Chrm2-expressing neurons. Along the longitudinal axis, the size and shape of the three SUB/ProSUB gene expression domains and their boundaries shift as the cytoarchitecture of the SUB and ProSUB changes (Fig. 3). However, the complementary arrangement of the gene expression domains to each other, as well as their relative position adjacent to the alveus and molecular layer, remains consistent. From the available consecutive tissue series covering 10 mm of the hippocampal longitudinal axis, the data suggests that the relatively segregated gene expression layers can be considered as continuous laminar sheets extending rostrocaudally across the whole SUB. Additionally, gene expression for Th and Pcp4 appear to mirror the distribution patterns of Nts and Chrm2, respectively (Fig. 3). Although it is unclear if these genes are expressed within the same individual cell types or shar a similar spatial distribution, Th and Pcp4 gene expression demonstrates that the SUB molecular domains represent differences in multiple combinatorial gene expression patterns.
Gene expression patterns in the anterior human SUB. Compared to the posterior hippocampus, the anterior hippocampus is structurally more complex including cytoarchitectural ridges and folding. At the anterior pole of the hippocampus, the transverse axis folds around the dentate gyrus and ultimately ends near the amygdala, an area commonly referred to as the hippocampal amygdala transition area (HATA). In coronal sections of the anterior hippocampus, the CA1/SUB is present on both dorsal and ventral sides of the hippocampal sulcus. In the ventral part of anterior SUB (Fig. 4), in situ hybridization reveals that Nts-, Htr2a-, and Chrm2expressing neurons are arranged in similar gene expression patterns to the posterior SUB levels, suggesting the ventral anterior SUB neurons are rostrocaudally continuous with neurons observed in the posterior SUB sections (Fig. 3). In contrast, the dorsal part of anterior SUB demonstrates a distinct combination of gene expression and cytoarchitecture. Examining the Nissl cytoarchitecture of this region reveals a distinct intermediate layer of densely packed and darkly stained pyramidal neurons that suggests a trilaminar organization (Fig. 4). Chrm2expressing neurons are distributed in the deep layer adjacent to the white matter whereas Htr2a-expressing neurons form a complementary gene expression pattern and are located within both the intermediate and superficial layer (Nts expression is mostly absent).
Examining tissue sections closer to the anterior pole reveals that the dorsal and ventral parts of the SUB are continuously adjoined. The layer of Chrm2-expressing neurons in the dorsal part of anterior SUB joins with the ventral part of SUB by continuing around the medial/anterior pole of the hippocampus (Fig. 5). In addition, the Nts-expressing neuronal layer also partially extends around the medial/anterior pole where these cells seem to be interposed between the layers of Chrm-and Htr2a-expressing neurons. At the anterior end of the SUB, Chrm2-expressing neurons are distributed as an outer ring layer with a more internal core of Nts-expressing neurons and a relative absence of Htr2a-expressing neurons. Overall, the distribution of gene expression patterns suggests that the SUB and ProSUB contain unique combinations of molecular domains composed of continuous sheets of distinct neuronal cell types spanning the entire hippocampal axis.

Comparison of gene expression-derived histological boundaries to human MRI imaging.
To determine if the regional boundaries identified by different gene expression patterns corresponded to other structural differences in connectivity microstructure, we compared our gene expression delineations to in vivo and ex vivo structural MRI and track density imaging (TDI) of hippocampal pathways (Fig. 6). At a rostrocaudal level similar to the posterior SUB described above, in vivo 7 T MRI images of human brain provide macroscopic resolution, sufficient to clearly distinguish the SUB molecular layer from the pyramidal layer, but difficult to distinguish finer details that are useful for delineating subregion boundaries (Fig. 6a,b). In comparison, ex vivo 16.4 T MRI of a dissected human hippocampal sample provides mesoscopic details. Hippocampal strata are clearly apparent, including the granule cell layer of the dentate gyrus, and the alveus can be distinguished from the angular bundle within the white matter (Fig. 6c). TDI analysis of the tissue microstructure in this sample reveals differences in the fiber pathway orientation and distribution (Fig. 6d,e) that notably align to our observed gene expression boundaries (Fig. 6f). The PRE contains dense bundles of dorsoventrally-oriented perforant path fibers that strongly contrast with more minimal fibers within the SUB. Within the boundaries of the Pro-SUB, thicker dorsoventally oriented fiber bundles are apparent that contrast it from the adjacent CA1 and SUB. Finally, a mediolateral fiber bundle can be observed running along the deep ProSUB and SUB just dorsal to the alveus that appears similar to the laminar distribution of Chrm2-expressing neurons. Using the Quantitative www.nature.com/scientificreports/ Imaging Toolkit (QIT) software 22 , we delineated the gene expression-based ROI across 600 μm anterior-posterior of the ex vivo MRI volume and quantified the fiber track density within each region (Fig. 6g). We found that each gene-expression defined regions of the SUB and ProSUB contained significantly different fiber density (Wilcoxon signed-ranked test (mean ± S.D.); SUB_3 = 29.2 ± 11.5, SUB_1 = 18.7 ± 6.5, SUB_4 = 38.1 ± 11.1; p < 0.0001 for all comparisons). Overall, the TDI fiber analysis provides evidence that microstructural differences in human hippocampal tissue correspond to differences in genetic cell type distribution when observed with high enough resolution.
Comparative analysis of the laminar organization in mouse and human SUB. In general, our observations of human SUB and ProSUB gene expression patterns are consistent with our previous analyses of gene expression patterns in the mouse SUB and we have maintained this color scheme in the human atlas drawings 6 (Fig. 7). The mouse HGEA defines 5 SUB subregions based on the representation and distribution of four distinct gene expression pyramidal sublayers. First, a dorsal and ventral part of dorsal subiculum (SUBdd and SUBdv) consists of a pyramidal layer with two gene expression sublayers (sublayers 1 and 4). The SUBdd most closely represents the area many studies refer to as the distal SUB, whereas ProSUB (sublayers 3 and 4) may closely correspond to proximal SUB in other studies. In contrast, the ventral subiculum (SUBv) has a trilaminar pyramidal layer (layers 2, 3, and 4). The thickness of these layers changes near the ventral tip of the subiculum (SUBvv) distinguishing this region from the SUBv. A major challenge for comparison analysis of mouse and human hippocampal gene expression data is the difference in orientation of the hippocampus between species. Because of the rotation of the longitudinal hippocampal axis across evolution, coronal mouse hippocampal sections are not in the same sectioning plane relative to the longitudinal axis as coronal human sections. Coronal However, several key aspects of SUB organization can be understood even comparing tissue that is sectioned at different angles relative to the longitudinal axis.
In both mouse dorsal SUB and human posterior SUB, Nts is strongly expressed in a superficial layer of SUB pyramidal neurons closest to the molecular layer, whereas Chrm2-expressing neurons are primarily located in the deepest part of the SUB, forming a continuous layer adjacent to the alveus white matter tract (Fig. 7a). Note that the position of the alveus and molecular layer are switched dorsoventrally in mouse vs. human tissue. This data suggests that the superficial SUB layer in humans (containing Nts-and Th-expressing neurons) is directly homologous to SUB sublayer 1 in the mouse HGEA, whereas the Chrm2-expressing deep layer in humans is homologous to mouse SUB sublayer 4. Therefore, Htr2a-expressing human SUB neurons likely correspond to the other two mouse SUB layers 2 and 3. However, Htr2a expression appears almost entirely absent in the mouse CA1 and SUB and other genes that are present in both mouse and human SUB (ex. Pcp4) are expressed in different combinatorial patterns (Fig. 7b). Overall, annotation of gene distribution within the SUB layers of both mouse and human subiculum suggests similarities and differences in combinatorial gene expression distribution (Supplementary Table 2). Together, this data suggests that while the anatomical laminar organization is evolutionarily conserved between mice and humans, there exists both conservation and divergence of the gene expression profile for each SUB layer and their unique cell types.
Although visualizing overall structure from thin-cut tissue sections can be difficult, the relationship of mouse and human SUB is easier to visualize in 3D and we have generated a model of the human SUB based on our understanding of the changes to the mouse SUB (Fig. 8a). In addition to the rotation of the longitudinal axis, one end of the hippocampus structure (the anterior pole in human corresponding to ventral pole in mouse) has folded back against the longitudinal axis. In comparison, the mouse and human SUB laminar organization www.nature.com/scientificreports/ demonstrate a similar subregional organization across the longitudinal axis despite its rotation (Fig. 8b). The area of the posterior human SUB containing gene expression layers 1, 3, and 4 correspond remarkably well to the SUBdd and ProSUB subregions in the dorsal part of the mouse SUB. In the anterior human SUB, the ventral area of SUB containing layers 1 and 4 is homologous to the mouse SUBdv whereas the dorsal part appears homologous to the trilaminar mouse SUBv/SUBvv. Additionally, the data from the human anterior pole is similar to the organization in the caudal parts of the mouse SUB where layers 1 and 4 partly curl around the ventral tip of the dentate gyrus. Overall, the data suggests the SUB has a conserved mammalian architecture that in humans has undergone expansion and structural folding.

Discussion
Overall, this study provides many novel insights into the organization of the human SUB and its evolutionary relationship to the mouse. Previous studies have analyzed gene expression patterns to map anatomical boundaries between the human SUB and ProSUB, however these studies did not report expression patterns of Chrm2 and other genes which are present within a subset of neurons in both regions and demonstrate their laminar organization [7][8][9] . Based on the current data, we believe that the structure of mouse and human SUB are highly conserved with respect to two major changes: (1) the rotation of the longitudinal axis and (2) the folding of the CA1/SUB at the anterior pole. The mouse and human in situ hybridization data for Nts, Chrm2, and Htr2a suggests that the anatomical laminar organization of SUB pyramidal neuron cell types is conserved between mouse and human although many similarities and differences in their gene expression profiles can be observed. Generally, we find that the posterior SUB region in humans contains two distinct complementary gene expression layers (characterized by Nts and www.nature.com/scientificreports/ Chrm2 expression) that suggests a homology to HGEA sublayers 1 and 4 found in the mouse SUBdd region. This gene expression data supports previous evidence for structural and functional homology between the posterior human SUB and dorsal mouse SUB and extends this relationship to the cellular level. Additionally, human sublayers 1 and 4 run continuously across the longitudinal axis as previously observed in mice, suggesting that the anterior ventral part of the human SUB is homologous to the recently identified SUBdv mouse region 6 . While sublayer 4 is present in the deepest part of the SUB at all levels, gene expression demarcating sublayer 1 is absent from the ProSUB and the more posterior parts of the dorsal anterior SUB near the amygdala. Instead, Htr2a was strongly expressed in superficial ProSUB neurons, similar in anatomical position to mouse HGEA sublayer 3 (although Htr2a is not expressed in mouse sublayer 3). Despite this difference, the overall Htr2a expression pattern in humans suggests an anatomical correspondence to HGEA sublayers 2 and 3 in mouse. In the dorsal anterior SUB region, the Nissl and gene expression data suggest three layers: a deep Chrm2-expressing layer, an intermediate Htr2a-expressing layer with distinctive Nissl cytoarchitecture, and a superficial Htr2a-expressing layer. The trilaminar organization of the anterior dorsal SUB region (and its adjacent location to the amygdala) suggests a homology to the mouse HGEA SUBv/SUBvv regions which are composed of sublayers 2, 3, and 4. Although additional data is needed, the evidence presented here establishes a foundation for a translational HGEA map for describing mammalian hippocampal architecture across species. The observed gene expression patterns highlight an important distinction between the dorsal and ventral parts of the human anterior SUB. The ventral part of the anterior SUB appears to be homologous to the SUBdv in mice whereas the dorsal part of the anterior SUB is homologous to the SUBv/SUBvv in mice. Based on the connectivity data in the mouse HGEA, this difference between the dorsal and ventral parts of the anterior human SUB implies differences in connectivity and function. If homologous to the mouse SUBv/SUBvv, then the neurons in the dorsal part of the anterior SUB would send output to limbic brain regions such as the amygdala, hypothalamus and prefrontal cortex. In contrast, the ventral part of the anterior SUB would be connected with the retrosplenial cortex and anterior thalamus and more related to visuospatial navigation (similar to SUBdv). In the future, human imaging studies of the SUB should separately analyze the dorsal and ventral parts of the anterior SUB to dissociate their anatomical/functional relationships.
Gene expression pattern analysis provides a useful way to map the anatomical distribution of distinct classes of neuronal cell types. Human imaging segmentation of the hippocampus and SUB are limited by the relatively low resolution of MRI contrast and the identification and position of hippocampal boundaries is highly inconsistent across studies 23,24 . Due to the difficulty of drawing hippocampal boundaries on relatively homogenous MRI contrast, the current consensus among the community is to follow geometrical selection guidelines [25][26][27][28][29][30] . However, a gene-expression based atlas could provide insight as to mechanistic and cell type-specific architecture of the hippocampus. Gene expression mapping provides a cellular resolution approach that could define 'ground truth' delineations and then be applied to MRI datasets. In support of this approach, we have found evidence that the gene expression boundaries reflect anatomical differences in the microstructure of human tissue (Fig. 6). Although it is currently unclear if these differences in fiber orientation reflect differences in long range connectivity or cell type-specific morphologic differences, our previous study found a strong relationship between anatomical connectivity and gene expression within the mouse hippocampus 6 . Notably, our previous study also found that mouse sublayer 4 neurons had a high degree of interconnectivity along the transverse hippocampal axis and the fiber tract identified in the TDI fiber analysis suggests a similar connectivity pathway exists in humans. As the resolution of MRI imaging technology continues to improve and microstructural imaging becomes clinically feasible and quantitatively-specific 31-33 , we believe the data here provides an important foundation for an integrative approach to bridge the gap between MRI imaging and human tissue histology.

Considerations for investigations of translational mammalian neuronal cell types.
Recently, a translational gene expression study comparing mouse and human suggested a conservation of cortical cell types with divergent features across mammals 34 . In the hippocampus, CALB1 has previously been shown to be selectively expressed in the human dentate gyrus (DG), whereas CALB1 in the mouse is expressed within the DG, CA1, and CA2 35 . These variations raise important considerations as to how to classify cell types across species: Do the CALB1-positive CA1 and CA2 neurons represent a distinctly murine neuronal cell type not found in humans or are the homologous cell types still present in the human, but simply no longer express CALB1? Although the current study focused on the SUB, we also observed several instances where the combinatorial expression of specific genes is different between mice and human across a variety of other hippocampal structures. Cross-species differences at the single-cell transcriptome level may underlie the failure to successfully translate mouse research to the clinic. Drugs that are effective in mice because they bind a receptor in specific target neurons may not work in humans because (1) either the homologous human neuron doesn't express that receptor or (2) the receptor is expressed in additional neurons causing unexpected side effects that were not observed in the mouse experiments. Because of these potential pitfalls in cross species drug development, it is critical that we understand the homologous single-cell gene expression profiles between mouse and human neuronal cell types. Ultimately, a convergence of data modalities describing anatomical, molecular, and physiological characteristics is imperative to accurately defining cell types across different mammalian nervous systems. The evidence presented here provides a foundation for understanding the organization of mouse and human subiculum, but further characterization is necessary for a comprehensive map of spatial gene expression patterns similar to the mouse HGEA 6 .
Conclusions. This study establishes a new understanding of the organization of the human SUB and provides a translational perspective that corroborates and extends the anatomical and functional homology across mouse and human SUB. We find that both mouse and human SUB pyramidal neurons contain a hidden laminar

Methods
Allen Brain Atlas in situ hybridization data. All human and mouse in situ hybridization images used in the analysis of this study were downloaded from the Allen Brain Atlas website (www.brain -map.org). Complete detailed information about histological processing and hybridization can be found in the Allen Institute white paper (Mouse: "In situ hybridization data production", Nov. 2011; http://help.brain -map.org/displ ay/mouse brain /Docum entat ion ; Human: "In Situ Hybridization in the Human Brian Atlas", Oct. 2013 v.7; http://help.brain -map.org/displ ay/human brain /Docum entat ion). Human hippocampal in situ hybridization image datasets are part of the Neurotransmitter Study (four postmortem subjects (H0351.1009, a 57 year old Caucasian male with hypertension; H0351.1010, 28 year old Hispanic male control; H0351.1012, 31 year Caucasian male control; H0351.1016, 55 year old Caucasian male control). All 88 probes (see Supplementary Table 1 gene list) and additional histological stains are performed in sequential tissue sections in order from anterior to posterior (Nissl stained sections occur every 25 sections). Each tissue section is hybridized to one unique probe and the procedure is repeated every 100 sections such that the rostrocaudal sequence of each gene expression tissue series is spaced 2 mm apart (100 sections x 20 μm section thickness). All image data used for analysis was openly published online and downloaded from www. brain -map.org.
Post-mortem tissue and ex vivo imaging parameters. The hippocampal specimen used for ex vivo imaging (a 55-year-old male with normal cognitive function) was acquired post-mortem from the Alzheimer's Disease Research Center (ADRC) Neuropathology Core at the Keck Medical Center of USC (Los Angeles, CA; NIA P50 AG05142). All experiments were reviewed and approved by the Institutional Review Board (IRB) of the University of Southern California (USC) and data was collected in accordance with the relevant guidelines of the USC IRB. Patient consent was obtained by the USC ADRC and de-identified prior to distribution. The specimens were immediately fixed in 10% phosphate-buffered formalin, and later dissected at the level of the lateral geniculate nucleus, including the hippocampal and parahippocampal gyri. A 1 cm extent of the hippocampal tissue was submitted for ex vivo MRI analysis.
DW-MRI was acquired using 3D Stejskal-Tanner 38 spin-echo sequence to achieve high signal-to-noise ratio (SNR) at (150 μm) 3 isotropic resolution. Data was acquired with three diffusion weightings: b-values of 1000, 3000 and 5000 s/mm 2 , with 20, 30 and 45 diffusion encoding gradient directions with distinct optimized spherically even distribution 39 . In addition, a total of 6 unweighted images (B0) were acquired. 3D DW-MRI spin-echo was acquired using FOV = 25.5 × 20 × 23.4 cm and matrix size = 170 × 130 × 156. TE, pulse duration and separation times were fixed across shells to avoid time-dependent effect on diffusion signal. Total DW-MRI scan time was 48 h and 22 min (also at 22 °C).
Ex vivo MRI data analysis. Image volumes were corrected for Gibbs ringing (using Trapezoid windowing) and N4 field bias (using ANTs software: http://stnav a.githu b.io/ANTs/). DW-MRI data quality and signal values were visually inspected for quality control. Track-density imaging (TDI) 40 was performed using single shell DW-MRI on the shell with b-value of 5000 s/mm 2 , using MRtrix software (version 0.2.12; http://jdtou rnier .githu b.io/mrtri x-0.2/index .html). Voxels with FA > 0.7 were segmented and the spherical harmonic decompositions of all the resulting profiles were then averaged to estimate the response function. We then applied constrained spherical deconvolution 41 to estimate the fiber orientation distribution in each voxel using a maximum spherical harmonic of order 6. Then, 500,000 streamlines were generated using probabilistic tractography tool 42 with Figure 7. Homologous SUB laminar organization with similarities and differences in gene expression. (a) Comparison of Nts and Chrm2 expression in the mouse dorsal subiculum (specifically SUBdd in the HGEA nomenclature, top row) and human posterior SUB (bottom row). Despite the difference in cell packing density, the laminar organization of the SUB pyramidal cells appears conserved across species. In the mouse HGEA, the Nts-expressing cells are located in SUB layer 1 (SUB_1), whereas Chrm2 expression is located in SUB layer 4 (SUB_4). In both mouse and human, the Nts-expressing cells are located superficially near the molecular layer (m) whereas the Chrm2-expressing cells are located deep near the alveus (alv). The position of the alveus and molecular layer (as well as the laminar organization of the pyramidal layer) is dorsoventrally inverted due to the rotation of the hippocampal axis between the two species. (b) In contrast to the highly similar Nts and Chrm2 expression patterns, Htr2a and Pcp4 are differentially expressed between the mouse (top, sagittal-cut images, HGEA atlas to left)) and human SUB (bottom, coronal sections). Htr2a expression is not strongly expressed in the mouse SUB (some minor expression in ventral CA1), but strongly expressed in human CA1 and SUB (bottom left). Pcp4′s combinatorial expression pattern in SUB is bilaminar in the mouse (SUB_2 and SUB_4, top right) but only expressed in the deep layer of human SUB cells (corresponding to SUB_4 only, bottom right). www.nature.com/scientificreports/ the following parameters: curvature = 0.075, cutoff = 0.1, minlength = 1, length = 15, step = 0.015. TDI with voxel size of (100 mm) 3 was then derived from generated streamlines. Gene expression-based domain volumes within ProSUB and SUB were manually delineated on MRI dataset from images spanning 600 μm anterior-posterior using the Quantitative Imaging Toolkit (QIT) software (SUB_1 = 379 voxels, SUB_3 = 499 voxels, SUB_4 = 276 voxels) 22 . The track density index is computed as the number of reconstructed streamlines in each voxel generated from the TDI technique. The average track density index of these regions was averaged across all voxels and compared using Wilcoxon signed-rank test (mean ± S.D.).
In vivo high resolution hippocampal imaging. A 32-year-old female was scanned using T2-weighted turbo spin echo sequences with 2 mm slice thicknesses and 340 µm (interpolated to 170 µm) in-plane resolution, 4 averages, resulting to total scan time of 13 min 43,44 . We used a single-channel quadrature transmit radiofrequency (RF) coil and a 32-channel receive array coil (Nova Medical Inc., MA). All methods were carried out in accordance with relevant guidelines and regulations and approved by the institutional review board of the University of Southern California. Informed consent was obtained from the volunteer, and the image was anonymized.  www.nature.com/scientificreports/ Human SUB structural modeling and 3D rendering. The human SUB 3D model was created by morphing the anatomical positions of the previously generated HGEA mouse SUB 3D model 6 to align with the distribution of gene expression lamina observed in the coronal in situ hybridization tissue section dataset. Using Autodesk 3ds Max animation software, vertex positions were assigned to the mouse SUB model and the longitudinal axis of the model was rotated to the orientation to correspond to the position of a reference hippocampal MRI volume and the anterior hippocampal pole was folded back to the approximate rostrocaudal level where dorsal sections of the SUB are present in the human hippocampus sections. Finally, the Morpher modifier, which is able to morph a model between several different shapes with corresponding vertices and orientation, was used to create an animated morph of the model transition state from the mouse shape and orientation to the human shape and orientation (Supplementary Movie 2).