Transversal functional connectivity and scene-specific processing in the human entorhinal-hippocampal circuitry
Abstract
Scene and object information reach the entorhinal-hippocampal circuitry in partly segregated cortical processing streams. Converging evidence suggests that such information-specific streams organize the cortical – entorhinal interaction and the circuitry’s inner communication along the transversal axis of hippocampal subiculum and CA1. Here, we leveraged ultra-high field functional imaging and advance Maass et al., 2015 who report two functional routes segregating the entorhinal cortex (EC) and the subiculum. We identify entorhinal subregions based on preferential functional connectivity with perirhinal Area 35 and 36, parahippocampal and retrosplenial cortical sources (referred to as ECArea35-based, ECArea36-based, ECPHC-based, ECRSC-based, respectively). Our data show specific scene processing in the functionally connected ECPHC-based and distal subiculum. Another route, that functionally connects the ECArea35-based and a newly identified ECRSC-based with the subiculum/CA1 border, however, shows no selectivity between object and scene conditions. Our results are consistent with transversal information-specific pathways in the human entorhinal-hippocampal circuitry, with anatomically organized convergence of cortical processing streams and a unique route for scene information. Our study thus further characterizes the functional organization of this circuitry and its information-specific role in memory function.
Editor's evaluation
Grande and colleagues provide important new insights into how different regions of the entorhinal cortex functionally interact with specific cortical brain areas and how, in turn, subregions of the entorhinal cortex interact with the hippocampus during 'scene' and 'object' processing. The study is well-motivated, well-designed, and provides convincing evidence using appropriate methodology. This paper is relevant to cognitive neuroscientists with an interest in the entorhinal cortex – hippocampal pathways and 'scene' and 'object' representation in the medial temporal lobe.
https://doi.org/10.7554/eLife.76479.sa0Introduction
Entorhinal and hippocampal subregions form a critical functional circuitry that binds cortical information into cohesive representations (Eichenbaum et al., 2007; Ritchey et al., 2015). The interaction of the entorhinal-hippocampal circuitry with large-scale cortical information streams and the circuitry’s inner communication are key to the formation of these cohesive representations. Here, we advance insight into how the human entorhinal cortex (EC) receives information from cortical streams and how information proceeds between the EC and the transversal axis of hippocampal subiculum and CA1 (here referred to as transversal sub/CA1 axis). These insights are relevant to our understanding of the circuitry’s fundamental role in cognitive functions such as episodic memory.
Large-scale cortical information streams, that originate in the visual ‘Where’ and ‘What’ pathways and process scene and object information (Berron et al., 2018; Haxby et al., 1991; Ranganath and Ritchey, 2012; Ritchey et al., 2015; Ungerleider and Haxby, 1994), map onto the EC in a complex manner and define functional EC subregions. [Note, in light of confusing nomenclature, here we adhere to scene and object information - elsewhere referred to as contextual, spatial or "Where" and content, non-spatial, item or "What" information, respectively.] Recent rodent research updates the former conception of a parallel mapping of scene and object information via parahippocampal and perirhinal cortices onto medial versus lateral EC subregions (cf. posterior-medial versus anterior-lateral EC subregions as the human homologues; Maass et al., 2015; Navarro Schröder et al., 2015). Instead of a strict parallel mapping, profound cross-projections exist from the parahippocampal cortex towards the perirhinal cortex and the lateral EC (Nilssen et al., 2019). In accordance, information seems to converge in the rodent lateral EC (Doan et al., 2019). The update, thus, implies a more complex functional organization than parallel scene and object information mapping. Moreover, this advance highlights the retrosplenial cortex as an additional source to convey information directly from the cortical scene processing stream onto the EC. The retrosplenial cortex projects to the medial EC and, like the parahippocampal cortex, is part of the scene processing stream (e.g. involved in scene translation; Vann et al., 2009; Nilssen et al., 2019; Witter et al., 2017). The update, furthermore, evokes the question how cortical sources of information uniquely map onto the EC and which kind of information is processed in the resulting functional EC subregions.
Within the entorhinal-hippocampal circuitry, an important direct way of communication exists between the EC and hippocampal subiculum and CA1. How functional EC subregions communicate towards the transversal sub/CA1 axis in humans is, however, unclear. Similarly, the extent to which specific scene and object information processing routes might emerge, despite information convergence in the EC, is unknown. On one hand, rodent research indicates a transversal organization where scene and object information is processed along two anatomically wired routes, the medial EC – distal subiculum – proximal CA1 route and the lateral EC – proximal subiculum – distal CA1 route, respectively (Witter et al., 2017; note sparse functional evidence in the subiculum: Ku et al., 2017; Cembrowski et al., 2018; but frequent reports in the rodent CA1 region: Henriksen et al., 2010; Nakamura et al., 2013; Igarashi et al., 2014; Nakazawa et al., 2016; Beer et al., 2018). Initial functional and structural connectivity data also indicate such a transversal connectivity profile in humans (Maass et al., 2015; Syversen et al., 2021). In accordance, scene information seems to be preferentially processed in the distal subiculum (Dalton et al., 2018; Dalton and Maguire, 2017; Zeidman et al., 2015) and hints exist for preferential object processing at the subiculum/CA1 border (Dalton et al., 2018). On the other hand, anatomical projections in the monkey show a longitudinal profile on top of the transversal profile with mainly the anterior-lateral and posterior-lateral entorhinal portions projecting to the distal subiculum – proximal CA1 and proximal subiculum – distal CA1, respectively (Witter and Amaral, 2020). According to information convergence in the EC, a recent report finds convergence along the rodent transversal CA1 axis (Vandrey et al., 2021). In humans, visual stream projections towards the entorhinal-hippocampal circuitry similarly suggest convergence of scene and object information in the subiculum/CA1 border region but preserved scene processing in the distal subiculum (Dalton and Maguire, 2017). A detailed examination of the latter hypothesis is, however, lacking. The diversity of findings emerging from the literature calls for a thorough investigation to elucidate whether multiple transversal processing routes exist within the human entorhinal-hippocampal circuitry.
To summarize, our conception of how information travels towards the entorhinal-hippocampal circuitry underwent key changes which warrant an extensive exploration of the circuitry’s functional organization. First, rodent research shows that there is no strict parallel mapping of cortical information from the perirhinal and parahippocampal cortex towards the EC. Second, information seems to converge already before the hippocampus. These changes add to several knowledge gaps. First, it is unclear in which subregions of the entorhinal-hippocampal circuitry scene and object information are processed. The general connectivity patterns in the human entorhinal-hippocampal circuitry have not yet been directly related to information processing. Moreover, it is unclear how scene information from the retrosplenial cortex maps onto the human EC as a critical source of the cortical scene processing stream. Hence, it is also unclear how retrosplenial information is communicated between the EC and the hippocampus. Finally, it remains elusive whether a transversal functional segregation can be extended towards the human CA1 region in analogy to the rodent literature.
Here, we leverage ultra-high field 7 Tesla functional imaging (fMRI) data and advance the earlier findings on human entorhinal subregions and a transversal intrinsic functional connectivity pattern in the subiculum (Maass et al., 2015). With a combination of functional connectivity and information processing analyses, we seek to answer two sets of questions. Regarding functional connectivity, we ask where the parahippocampal, perirhinal, and retrosplenial cortical sources uniquely map onto the human EC and how these functionally connected routes continue between EC subregions and the transversal sub/CA1 axis. Regarding information processing, we ask whether and where scene and object information are specifically processed in the EC and along the transversal sub/CA1 axis. We test the hypotheses of (1) a transversal functional connectivity pattern and (2) multiple information processing routes within the entorhinal-hippocampal circuitry. Thus, following the updated conception of a non-parallel cortical scene and object information mapping onto the EC in rodents, we will show how cortical information streams map onto the EC in humans. This mapping will then be our detailed starting point to investigate the functional connectivity and information processing within the entorhinal-hippocampal circuitry.
Results
We seek to comprehensively investigate functional connectivity within the entorhinal-hippocampal circuitry and the contribution of cortical scene and object information processing streams. In an initial step, we identified where cortical sources map uniquely onto the entorhinal cortex (building upon Maass et al., 2015). The identified entorhinal subregions are based on their voxel’s preferred intrinsic functional connectivity with the retrosplenial cortex, parahippocampal cortex, perirhinal Area 35 or Area 36 regions (‘sources’). Note, that we evaluate both perirhinal subregions, Area 35 and Area 36 as separate sources as accumulating evidence suggests their structural and functional distinction (e.g. Berron et al., 2021; Burwell, 2000; Suzuki and Naya, 2014; van Strien et al., 2009). Next, we evaluated the continuation of the functional connectivity streams within the entorhinal-hippocampal circuitry and examined the intrinsic functional connectivity pattern between the identified entorhinal subregions (‘seeds’) and hippocampal subiculum and CA1 in the hippocampal body. Therefore, temporal fluctuations of BOLD signal were correlated in a seed-to-voxel manner within each participant. The resulting statistical correlational maps were aligned between participants. Repeated measures ANOVAs were calculated on connectivity preferences with seeds and transversal segments as factors to determine statistical differences in connectivity topography. All functional connectivity analyses were performed on the dataset after task-related effects have been regressed out, creating a dataset that resembles resting-state data (Gavrilescu et al., 2008; Maass et al., 2015).
Finally, we identified the corresponding bias in scene (here operationalized with room stimuli) and object information processing within entorhinal subregions and along the transversal sub/CA1 axis using the same dataset. Therefore, we extracted parameter estimates from a mnemonic discrimination task with scene and object conditions from aligned statistical maps across participants. Repeated measures ANOVAs were calculated on parameter estimates in entorhinal subregions and transversal sub/CA1 segments to determine biases in information processing within the entorhinal-hippocampal circuitry.
In the following, we first describe the four obtained entorhinal seeds and display the intrinsic functional connectivity pattern with the entorhinal seed regions along the transversal sub/CA1 axis. Thereafter, we report the information processing characteristics of the entorhinal and hippocampal subregions. Note that all results have been obtained with independent analyses in the left and right hemispheres. The largely similar left hemisphere results can be found in appendix 1. Source data and statistical maps are provided under Grande, X., Berron, D. (2022). Open Science Framework. ID 9v3qp. Source Data from Functional Connectivity and Information Processing in the Entorhinal-Hippocampal Circuitry. https://osf.io/9v3qp.
Four cortical sources divide the EC in retrosplenial-, parahippocampal, Area 35- and Area 36-based seeds
The four entorhinal subregions that we later used as seeds to determine the topography of entorhinal-hippocampal connectivity are based on intrinsic functional connectivity preferences with either the parahippocampal cortex, the retrosplenial cortex, perirhinal Area 36 or Area 35. These cortical regions are in general concordance with Maass et al., 2015 but consider recent advances that put forward the retrosplenial cortex as a critical source from the cortical scene processing stream (Nilssen et al., 2019) and evaluate perirhinal Area 35 and 36 separately.
Based on functional connectivity preferences with the four sources - parahippocampal cortex (Source code 8), retrosplenial cortex (Source code 7, Area 36 (Source code 6), and Area 35 (Source code 5) - we obtained four entorhinal seeds. The seeds refer to different parts of the EC whose voxels expressed preferential functional connectivity to either cortical source. For the ECPHC-based seed, the majority of voxels can roughly be described as clustering in the posterior-medial entorhinal portion, for the ECRSC-based seed in the anterior-medial portion, for the ECArea35-based seed in the anterior-lateral portion and for the ECArea36-based seed in the posterior-lateral entorhinal portion (see appendix 2 for exact voxel counts). Note that both perirhinal-based entorhinal seeds extended along the anterior to posterior axis such that the ECArea35-based progressed more along the outer EC (i.e. laterally, with a main focus anteriorly) and the ECArea36-based along the inner EC (i.e. medially, with a main focus posteriorly, see Figure 1 and the medial reflection of the EC seeds). It is important to note that these are rough qualitative descriptions of the main clusters, without quantification or an established relationship to coherent cytoarchitectonic regions. We will therefore continue to refer to them as ECRSC-based, ECPHC-based, ECArea35-based and ECArea36-based seeds.
Distal subiculum is functionally connected with the ECPHC-based seed while the subiculum/CA1 border is connected with ECRSC-based and ECArea35-based seeds
Following the characterization of entorhinal seeds, we focused on the functional connectivity between these entorhinal subregions and hippocampal subiculum and CA1 to test the hypothesis of a transversal functional connectivity pattern. We predicted that while some EC subregions have a preference to functionally connect with the subiculum/CA1 border, others preferentially connect with the distal subiculum and proximal CA1. In the previous step, we identified EC subregions based on unique cortical source contributions. Therefore, our predictions remained in accordance with Maass et al., 2015: We expected that the EC subregion preferentially connected with the parahippocampal cortex (ECPHC-based seed) maps towards the distal subiculum and EC subregions connected with the perirhinal cortex (ECArea35-based seed, ECArea36-based seed) map towards the proximal subiculum, a mapping that we predicted to be extended towards the distal CA1.
When extracting estimates of connectivity preferences across individuals from proximal and distal hippocampal subfield segments for either entorhinal seed, repeated measures ANOVAs revealed significant seed X segments interaction effects along the transversal sub/CA1 axis (see Figure 2; subiculum: F(12,372) = 19.561; p<0.001; CA1: F(6,186) = 3.212; p=0.024).
In the subiculum, additional repeated measures ANOVAs showed that the ECArea35-based (F(4,124) = 8.913; pFDR <0.001), ECRSC-based (F(4,124) = 10.538; pFDR <0.001) and ECPHC-based (F(4,124) = 42.201; pFDR <0.001) seeds displayed a significant main effect across the transversal subiculum segments. These differential functional connectivity preferences across the transversal axis of the subiculum interacted significantly in a subsequent repeated measures ANOVA (ECPHC-based versus ECRSC-based seed preference interaction: F(4,124) = 46.452; pFDR <0.001; ECPHC-based versus ECArea35-based seed preference interaction: F(4,124) = 35.208; pFDR <0.001). This pattern provides statistical evidence for an increase in preferential functional connectivity with the ECPHC-based seed towards the distal portion of the subiculum while the preferential functional connectivity with the ECArea35-based as well as the ECRSC-based seeds rather increased towards the proximal portion of the subiculum.
In hippocampal CA1, additional repeated measures ANOVAs showed that the connectivity preference towards the ECRSC-based seed displays a significant main effect across the transversal axis of CA1 (F(2,62) = 10.489; pFDR <0.001). In distal CA1, the preferential functional connectivity with the ECRSC-based seed was higher than in the proximal portion of CA1. In right CA1, a similar but weaker transversal pattern was observed for connectivity preferences with the ECArea35-based (F(2,62) = 4.146; pFDR = 0.041; note in the left hemisphere a comparable transversal pattern was observed for the ECPHC-based and ECRSC-based portions, see appendix 1).
Thus, in the entorhinal-hippocampal circuitry, voxels in the distal subiculum were preferentially functionally connected with the ECPHC-based portion whereas voxels in the subiculum/CA1 border were preferentially connected with more anterior EC portions (ECRSC-based and ECArea35-based).
Distal subiculum and ECPHC-based exhibit higher functional activity in the scene condition while other subregions show no significant difference between conditions
Besides the intrinsic functional connectivity patterns within the entorhinal-hippocampal circuitry, we also examined the characteristics of scene and object information processing to test the hypothesis of multiple information processing routes within the entorhinal-hippocampal circuitry. We predicted a route of specific scene processing and another route of convergent information processing. Following the proposal by Dalton and Maguire, 2017 and the updated cross-projections from the scene to the object information processing stream (Nilssen et al., 2019), we expected scene processing in the distal subiculum. The updated parahippocampal cross-projections imply convergence wherever specific object processing had been expected previously. Thus, we explored whether any entorhinal-hippocampal subregions still process object information specifically. However, we largely expected to find evidence consistent with convergent processing of scene and object information within the entorhinal-hippocampal circuitry.
We first focused on the entorhinal seed regions. When extracting task-related parameter estimates from object and scene conditions, a repeated measures ANOVA showed a significant interaction between region and information type (object versus scene, F(3,93) = 20.9267; p<0.001). Post-hoc t-tests revealed that only in the ECPHC-based seed region functional activity in the scene condition was significantly higher than in the object condition (pFDR<0.001), while in the remaining three entorhinal seed regions no significant difference between scene and object conditions existed (ECArea35-based: pFDR = 0.9129; ECArea36-based: pFDR = 0.9129; ECRSC-based: pFDR = 0.5646; see Figure 3).
When extracting task-related parameter estimates for scene and object conditions from proximal and distal segments of hippocampal subregions within each participant, we found a significant interaction between transversal segments and information type only in the subiculum (F(4,124) = 15.994; p<0.001) and not in CA1 (F(2,62) = 2.553; p = 0.105) as revealed by repeated measures ANOVAs. Post-hoc T-tests showed significantly higher functional activity in the scene than object condition (both pFDR <0.001) only in the distal subiculum segments. In all other segments along the subiculum transversal axis, there was no significant difference in functional activity related to scene and object conditions (all pFDR = 0.1222; see Figure 4).
Discussion
This study aims to advance insight into the organizational principles of information processing within the entorhinal-hippocampal circuitry and the circuitry’s embedding in designated cortical processing streams. Leveraging ultra-high field 7 Tesla fMRI, we find a resemblance between the intrinsic functional connectivity pattern and subregional biases in scene information processing in the entorhinal-hippocampal circuitry. In the EC, we observe a topographical mapping of regions from the cortical scene and object information processing streams, including the retrosplenial, parahippocampal and perirhinal Area 35 and Area 36 cortices. This mapping continues to determine a transversal organization of information processing routes between the EC and the human hippocampal circuitry. Our results unify previous evidence and uncover novel features in the human brain that can be a window into the circuitry’s critical role in memory function and decline.
Scene information is processed within an ECPHC-based – distal subiculum route
We identified regions in the entorhinal-hippocampal circuitry that are dedicated to process scene information. These regions consisted of two functionally connected portions: the ECPHC-based and the distal subiculum. The subiculum showed a transversal difference in intrinsic functional connectivity with a preference to the ECPHC-based in its distal portions (of note: the ECPHC-based was defined by entorhinal voxels with preferential functional connectivity to the parahippocampal cortex). Importantly, the distal subiculum and the ECPHC-based were the only studied entorhinal-hippocampal subregions that exhibited functional activity specifically in the scene condition (see appendix 5 for information processing in cortical source regions).
These findings provide clear evidence for a hypothesized transversal difference in scene information processing within the human subiculum. Our data also replicate the earlier functional and structural connectivity reports in humans as well as anatomical findings of a route between posterior-medial EC (based on parahippocampal connectivity) and distal subiculum (Maass et al., 2015; Syversen et al., 2021; Witter et al., 2000). The scene information processing bias has mainly been previously reported for the EC (in rodents, operationalized by spatial processing conditions: Neunuebel et al., 2013; Keene et al., 2016; in humans, operationalized by scene stimulus conditions: Berron et al., 2018; Navarro Schröder et al., 2015; Reagh and Yassa, 2014; Schultz et al., 2015). In animal studies, the importance of the subiculum as a translator of hippocampal information towards the entorhinal and other cortical structures is increasingly acknowledged (O’Mara, 2006; Roy et al., 2017). We here contribute to the sparse investigations regarding the nature of information processed along the transversal axis of the subiculum (see Cembrowski et al., 2018; Ku et al., 2017). Our observation is in line with the hypothesis that the distal subiculum is more involved in processing scenes than objects based on previous findings in the human brain. While the subiculum in general was associated with scene discrimination (Hodgetts et al., 2017), a growing body of evidence relates particularly the medial hippocampus to scene processing. This entails two medial areas, the pre- and parasubiculum, that we attribute to the distal subiculum in our current segmentation. Especially the area that resembles the pre- (or here: distal) subiculum has been shown to be involved in scene construction (Dalton et al., 2018; Zeidman et al., 2015). Recently, a gradient with coarser voxel-wise autocorrelation signals in the medial hippocampus has been reported, a finding that implies larger representations in the distal subiculum (Bouffard et al., 2022). In the latter two studies, however, the authors did not specifically extract data from the transversal axis of hippocampal subfields. Our joint investigation of functional entorhinal-subiculum connectivity and type of information processing along the full transversal axis of the subiculum, is the first to show a clear preference of scene information towards the distal portion, in comparison to more proximal portions.
Information processing is consistent with convergence within the anterior entorhinal portions – subiculum/CA1 border route
Our data revealed a route that did not show differences in scene and object information processing. Both, the ECArea35-based and the ECRSC-based portion exhibited preferential functional connectivity with the subiculum/CA1 border. Comparable levels of functional activity in scene and object conditions along these entorhinal-hippocampal routes are consistent with information convergence.
While we again confirm earlier findings and previously stated hypotheses, several features in our data are fundamentally novel. First, we provide initial human evidence for a functional connection between the ECArea35-based and the subiculum/CA1 border. Non-primate and primate anatomical data as well as ex vivo and in vivo structural connectivity data in humans show the possibility of information flow along that route (Syversen et al., 2021; Witter et al., 2017; Witter and Amaral, 1991; Witter and Amaral, 2020). Our results now underpin a functional relevance of that connection beyond the subiculum (for the subiculum see Maass et al., 2015). Our findings are derived based on a voxel-wise analysis, unconstrained by a priori selection of regions-of-interest. We thereby confirm the long-held proposal of a transversal functional organization in human subiculum and CA1.
Convergence of scene and object information is compatible with recent rodent work that shows joint coding of scene and object information (notably operationalized as spatial and non-spatial information, respectively) along CA1 and within the lateral EC (Deshmukh, 2014; Doan et al., 2019; Vandrey et al., 2021; Wilson et al., 2013; Yeung et al., 2019). Note that a supplemental analysis of information processing in the cortical source regions showed indeed, specific object processing in perirhinal source regions (see appendix 5). The lack of increased object processing in the anterior EC subregions and subiculum/CA1 border is thus likely not a result of increased noise in the object condition. Instead, increased object processing in perirhinal cortical source regions indicates subsequent convergence in entorhinal-hippocampal subregions, as hypothesized based on the updated cortical mapping scheme onto the EC.
Our results cannot confirm previous reports about higher functional activity for object than scene processing within these areas in the human brain (Reagh and Yassa, 2014; Navarro Schröder et al., 2015; Berron et al., 2018; also indicated in Dalton et al., 2018 and Schultz et al., 2015). Neither did we observe proximodistal differences in CA1 for object versus scene information processing as suggested by several rodent studies (Beer et al., 2018; Henriksen et al., 2010; Nakamura et al., 2013; Nakazawa et al., 2016). Differences in experimental design and contrasts could have contributed to these discrepancies (i.e. specific object information processing versus convergence). Previous studies used a variety of different conditions to tackle scene and object information processing (e.g. objects in time versus objects in space in Beer et al., 2018 or imagined objects on a 2D versus 3D grid in Dalton et al., 2018). In contrast to the current data, previous human studies did not derive functional data from specific, functionally defined entorhinal portions in the same dataset. As most previous studies were conducted in the light of the ‘parallel mapping hypothesis’, the related assumptions influenced the examined subregions, which may have altered the extracted measures.
Regarding the human proximal CA1, a firm conclusion is limited with our data. First, the functional connectivity results varied between hemispheres. In both hemispheres, proximal CA1 showed a different connectivity profile compared to distal CA1. However, even though statistically not significant, the preferences at the group level indicated increased functional connectivity with the ECPHC-based portion in the right proximal CA1 but with the ECArea35-based portion in the left hemisphere. Second, we do not prove similar information processing along the transversal CA1 axis. Instead, we find no significant difference in information processing along the transversal CA1 axis. As indicated in the previous paragraph, we cannot rule out that our object versus scene processing conditions may not have been sensitive enough to tackle functional differences in CA1. Thus, future research will have to identify defining characteristics of information processing along the transversal CA1 axis in a less constraint manner to allow conclusions on distinct information processing in proximal CA1.
In addition, we observed an unreported resemblance in functional connectivity profiles of ECRSC-based and ECArea35-based portions in the anterior EC. The sources of these entorhinal portions are part of cortical scene and object processing streams, respectively (see also appendix 5 that shows increased scene processing in the retrosplenial and increased object processing in perirhinal cortical source regions). To our knowledge, the ECRSC-based portion has not yet been identified in earlier investigations. While anatomical projections from the retrosplenial to deep medial EC layers have been confirmed in rodents, they appear in the posterior EC (Czajkowski et al., 2013; Sugar et al., 2011). Very recently, Syversen et al., 2021 found structural connectivity between the human retrosplenial cortex and the medial EC, but again not in the anterior part of the EC. The EC segmentation of Syversen and colleagues, however, followed different rules which may have contributed to differences in the topographical evaluation of the region. Also, structural and functional connectivity methods may yield different results, in particular as we identified EC subregions with a different set of cortical source regions. Under the assumption that retrosplenial connectivity defines the medial EC (Witter et al., 2017), the mapping of the ECRSC-based to the subiculum/CA1 border opposes conventional views that the medial EC communicates with the distal subiculum and proximal CA1 (based on rodent anatomy – see e.g. Nilssen et al., 2019). Whether species differences exist in the retrosplenial cortex – EC – hippocampus connectivity pattern or whether functional and structural connectivity diverge needs further investigation in the future.
Relevance of the current findings for the functional organization of the entorhinal-hippocampal circuitry
The current findings (summarized in Figure 5) advance our insight into the organization of the entorhinal-hippocampal circuitry on multiple levels and contribute to cognitive and clinical research. Recent efforts to understand how the human entorhinal-hippocampal circuitry accomplishes conjunction and segregation of information largely focused on the longitudinal hippocampal axis (e.g. Brunec et al., 2018; Brunec et al., 2020; Robin and Moscovitch, 2017). The transversal axis of the hippocampus has been approached by studies in humans that did not directly relate connectivity findings to information processing and did not assess subfield-specific organization (Vos de Wael et al., 2018; Plachti et al., 2019; Kharabian Masouleh et al., 2020; Paquola et al., 2020; Bouffard et al., 2022; for an overview see Genon et al., 2021). Dalton and Maguire, 2017, however, made a relevant proposal based on visual processing pathways and information processing. In correspondence to our results, they proposed the subiculum/CA1 border as a point of convergence between scene and object information processing streams. While their conclusion was based on direct parahippocampal, retrosplenial and perirhinal connections to the hippocampus, we found that both, the ECArea35-based (that is connected with the cortical object processing stream) and the ECRSC-based (that is connected with the cortical scene processing stream) show connectivity with the subiculum/CA1 border (see also appendix V for information processing in cortical source regions). Convergence is potentially also achieved via recurrency within the entorhinal-hippocampal system and cortical regions (cf. Koster et al., 2018 for evidence on recurrency). These considerations are an exciting future research avenue and remain speculative based on the current data due to insufficient temporal resolution. We nevertheless hypothesize the existence of two processing routes: one that processes converged object and scene information and one that processes scene information specifically. Thus, scene and object information processing might converge before the hippocampus. This presumably occurs within the anterior EC, given object-specific and scene-specific processing take place in the cortical source regions of the ECArea35-based and ECRSC-based portions, respectively (see appendix 5). Here, objects may be bound together with their defining scene-like or contextual features (akin to the ‘object-in-location’ idea in Connor and Knierim, 2017; Knierim et al., 2014). In addition, the dedicated scene processing that we observe along the ECPHC-based – distal subiculum route, may functionally underpin ideas about an anatomically graded contextual scaffold that the hippocampus utilizes to incorporate detailed information from the object-in-scene route into meaningful chunks of cohesive memory representations (‘events’; Behrens et al., 2018; Clewett et al., 2019; Robin, 2018; Robin and Olsen, 2019).
For completeness, we noted differences in functional connectivity along the longitudinal axis of the subiculum. We observed, for instance, more widespread functional connectivity of the ECArea35-based in the posterior subiculum whereas functional connectivity with the ECPHC-based portion seems more prominent in the anterior subiculum. The latter is consistent with previous reports (Dalton et al., 2019). The former, however, needs to be explored further by taking different segmentation protocols and seed regions into account. Note, that Maass et al., 2015 did not report longitudinal differences in connectivity strength between the EC and the subiculum. Future work needs to investigate how these observations relate to the reported gradient in functional connectivity and information resolution along the hippocampal longitudinal axis (e.g. Brunec et al., 2018 but many more). Altogether, the functional organization indicates that when a memory is to be formed, some degree of information convergence happens already before the hippocampus, nevertheless keeping specific aspects of scene information separated. This conclusion is in accordance with the updated cortical mapping scheme onto the EC (Nilssen et al., 2019). The topographical specificity of our results supports the necessity of functionally assessing the entorhinal-hippocampal circuitry with high spatial resolution and investigate memory function at a subregional level (Lee et al., 2020). The features we identified can inform future hypotheses on how the hippocampus achieves the formation of cohesive representations that serve memory function.
From a clinical research perspective, it is remarkable that the current functional connectivity pattern resembles the topology of early cortical tau pathology in Alzheimer’s disease (Lace et al., 2009). An influential hypothesis suggests tau progression in Alzheimer’s disease along functionally connected pathways in the human brain (Franzmeier et al., 2020; Vogel et al., 2020). Earliest cortical tau pathology in Alzheimer’s disease accumulates in perirhinal Area 35 (also referred to as transentorhinal region) and the anterior-lateral EC before it can be found along the subiculum/CA1 border (Braak and Braak, 1995; Berron et al., 2021; Kaufman et al., 2018; Lace et al., 2009). The topology of early tau pathology in Alzheimer’s disease thus mirrors the regions that we find biased towards ECArea35-based connectivity (Braak and Braak, 1991; Lace et al., 2009; Roussarie et al., 2020). Tau pathology in Alzheimer’s disease is associated with memory impairment (Bejanin et al., 2017; Berron et al., 2021; Nelson et al., 2012) and information processing might be affected accordingly as reports have shown an association between Alzheimer’s related tau pathology and object memory in early disease stages (Berron et al., 2019; Maass et al., 2019). However, given our finding of activity patterns consistent with object – scene convergence in those subregions of the hippocampal-entorhinal circuitry that are affected by early tau pathology, object-in-scene memory tasks might have increased sensitivity to memory impairment. Moreover, both, the entorhinal portion based on retrosplenial connectivity (ECRSC-based) and the entorhinal portion based on Area 35 connectivity (ECArea35-based), are functionally connected to the subiculum/CA1 border. This overlapping functional connectivity pattern in the hippocampus might be a way along which tau and amyloid pathologies in Alzheimer’s disease could interact. This is consistent with early hypometabolism and cortical tau progression in the retrosplenial cortex and early amyloid in posterior parietal regions (Grothe et al., 2017; Palmqvist et al., 2017; Ziontz et al., 2021). The revealed functional connectivity and information processing profile may guide future hypotheses on the propagation of Alzheimer’s pathology and related functional and cognitive impairment.
Limitations
This study has a number of limitations. First, the biases in seed connectivity in the left hemisphere were generally weaker and proximal CA1 results were less consistent across hemispheres. We conducted all analyses independently for both hemispheres to allow internal replication of our findings, however, whether partially different effects indeed signal a lateralization of the entorhinal-hippocampal organization in humans or whether the task or another parameter influenced these observations, will require further research.
Second, while it is unlikely that our functional connectivity pattern is the result of spatial proximity, increased correlation between spatially adjacent regions is an inherent problem of functional connectivity analyses. Distances between seed and target regions differ and may determine patterns in the functional connectivity data. To diminish the influence of proximity, our smoothing kernel was smaller than two times the voxel size. It is important to stress moreover, that the pattern of our results is not easily explainable by spatial distance between seed and target regions. The ECArea35-based or ECRSC-based, for instance, are not adjacent to the subiculum/CA1 border. Furthermore, we observed roughly comparable results for neighboring seeds and targets (e.g. ECPHC-based and distal subiculum) when we performed the functional connectivity analyses with seed and source regions in the contralateral hemisphere.
Third, our perspective was entirely functional and we cannot conclude on the directionality of our results. Also, it was beyond the scope of this study to examine direct connectivity between the cortical sources and hippocampal subregions. To what extent there is correspondence to structural connectivity (Syversen et al., 2021) remains to be determined, considering differences in the experimental task constraints and contrasts. Note that as a first step towards an understanding of the system’s functional organization and to increase comparability with earlier studies, we assessed functional connectivity and information processing within the entorhinal-hippocampal circuitry with univariate methods. These allow relative comparisons between functional activity levels in different conditions. Consequently, we are neither able to assess what the EC is processing during the baseline condition, meaning the absolute level of functional activity, nor are we able to verify that information processing is similar across conditions in for example the ECArea35-based seed. Univariate methods, moreover, average the signal over regions of interest. To capture hidden voxel-wise patterns of activity that scale with the processing of certain representations, future studies could examine information pathways with multivariate methods that evaluate informational content in the activity pattern of voxels instead of in an averaged manner (Kragel et al., 2018; Kriegeskorte et al., 2008). Moreover, recent methodological advances can be employed in the future that study functional connectivity based on the underlying content representations between regions (Basti et al., 2020).
Fourth, our study was originally conducted within the assumption that (functional) connectivity profiles reveal functional subregions. Based on that approach, the medial EC is identified based on i.a. retrosplenial connectivity. We, therefore conclude a surprisingly anterior yet medial EC mapping of the retrosplenial cortex. This approach has been followed by Maass et al., 2015 and also in numerous anatomical connectivity studies in animals (see Witter et al., 2017). It is possible that species differences lead our ECRSC-based to be more anterior than one would expect based on animal studies. However, given that the medial subregion in the primate EC remains posterior (cf. posterior-medial EC homologue in Maass et al., 2015), another possibility is that our retrosplenial functional connectivity cluster maps onto the human anterior-lateral EC. Our data does not allow us to verify this latter option. It is unclear, however, why functional subregions in line with predictions from animal research can be identified for some cortical source-to-EC mappings (like the parahippocampal cortex) but not for others. In combination with closely matched histological or structural magnetic resonance imaging data, future work can further reveal the nature of retrosplenial mapping on the human EC.
In general, the quantification of the transversal connectivity pattern should be considered with some caution from the anatomist’s perspective. The segmentation of subregions on functional data is an approximation because the anatomical ground truth cannot be captured by any segmentation protocol (even histological data can lead to divergent opinions). This shortcoming is amplified by group comparisons that do not account for participant-specific anatomy. Future research is needed to evaluate how the functionally derived entorhinal seeds in this study relate to histologically derived entorhinal subregions (Oltmer et al., 2022) or entorhinal subregions based on structural connectivity (Syversen et al., 2021). For a dedicated comparison of subregions, it is essential to pay close attention to the segmentation of the EC itself. Note moreover, that we excluded the head and the tail of the hippocampus from our investigation. The head is highly complex in its subfield topography (Ding and Van Hoesen, 2015; Berron et al., 2017) and prevents clear hypotheses regarding a transversal pattern. For the tail we lack an established segmentation protocol (de Flores et al., 2020; DeKraker et al., 2018). In the future, advanced segmentation methods and evaluations in the participant-space will improve this issue and reveal the organization in more detail.
Conclusion
In sum, leveraging ultra-high field functional imaging, we provide a comprehensive in vivo exploration of the functional organization within the human entorhinal and hippocampal subregions and the circuitry`s embedding within cortical information processing streams. Within the entorhinal and hippocampal subiculum, our data partially support a continuation of cortical object and scene information processing with convergence in anterior and lateral entorhinal portions (ECArea35-based, ECRSC-based, ECArea36-based), proximal subiculum and CA1, while the posterior-medial entorhinal portion (ECPHC-based) and distal subiculum process scene information specifically. Topographically, this organization of information processing overlaps with our identified pattern of functional connectivity. The data yield spatially organized information processing along functionally connected subregions in the human EC and transversal sub/CA1 axis. Our high-resolution approach revealed unknown characteristics of functional connectivity and scene processing within the human entorhinal-hippocampal circuitry. These aid our understanding of how cortical information comes together and is further communicated within the entorhinal-hippocampal circuitry, underpinning the formation of cohesive memory representations. We provide essential insights for basic and clinical research that we believe to be crucial for the development of future hypotheses on memory function and decline.
Methods
The current data is part of a larger study that examines exercise effects on cognition. The data that is subject to the current study have been acquired during the baseline measurement before any intervention took place. In the following, we focus on the study setup and methodological aspects of direct relevance for the current questions and data analyses.
Participants
In total, 32 healthy participants (15 female) with a mean age of 25.5 years (range 19–35 years, standard deviation 4.3 years) were included in the current data analyses. All participants were right-handed, finished education on A-level (German Abitur or comparable) and reported absence of any neurological or psychiatric diseases. General exclusion criteria determined by the 7 Tesla MR scanning procedure were applied (e.g. metallic implants, tinnitus, known metabolic disorders). All participants gave informed consent prior to participation and received a monetary compensation. The study received approval by the ethics committee of Otto-von-Guericke University, Magdeburg (Germany) under reference number 128/14.
Task
While functional images were acquired, participants engaged in a mnemonic discrimination task (see Berron et al., 2018). The object-scene task consisted of 64 objects and 64 rooms. In two runs, participants encoded always two stimuli, two 3D rendered objects in the object condition or two 3D rendered rooms in the scene condition and subsequently identified the following two same or similar stimuli as novel or old. Ten scrambled images were presented in blocks at the beginning and end of each run and served as baseline condition. All stimuli were presented for three seconds. In the recognition phase, participants had to respond during that time. Each stimulus was followed by a noise stimulus to prevent after-image and pop-out effects. The short alternating encoding/recognition sequences were embedded in an event-related design.
Data acquisition
All MRI data was acquired with a 7 Tesla Siemens MR machine (Erlangen, Germany) using a 32-channel head coil. First structural images were obtained. A whole-brain MPRAGE volume was acquired with isotropic voxel size of 0.6 mm, TR 2500ms; TE 2.8ms, 288 slices in an interleaved manner (FOV 384x384x288). Thereafter, a partial structural T2*- weighted volume (TR 8000ms; TE 76ms, interleaved, 55 slices, FOV 512x512x55), orientated orthogonal to the main longitudinal hippocampal axis was obtained with a resolution of 0.4x0.4 mm in-plane and a slice thickness of 1 mm.
The subsequent acquisition of functional data took place in two runs à 14 min (332 volumes each) employing echo-planar imaging (EPI). The volumes were partial (40 slices, TR 2400ms, TE 22ms, FOV 216x216x40, interleaved slice acquisition), oriented along the longitudinal axis of the hippocampus and acquired with an isotropic voxel size of 1 mm.
All EPIs were distortion corrected with a point-spread function method and motion corrected during online reconstruction (Zaitsev et al., 2004).
Data analyses
Preprocessing
Preprocessing and statistical modeling of fMRI data was performed with SPM12 (Wellcome Department of Cognitive Neuroscience, University College, London UK; Penny et al., 2011). The individual functional images were slice time corrected and smoothed with a full-width half-maximum Gaussian kernel of 1.5 mm. To preserve a high level of anatomical specificity, smoothing was performed with a kernel smaller than two times the voxel size. The artifact detection toolbox ARTrepair (Mozes & Whitfield-Gabrieli, 2011) was subsequently used to identify outliers regarding mean image intensity and motion between scans (threshold in global intensity: 1.3%; movement threshold: 0.3 mm). Identified outliers are included as spike regressors in subsequent statistical modeling.
Task effects in the functional data were removed by fitting general linear models (with regressors for all task conditions, outliers and movement parameters) to the data. The obtained residual images were saved for the intrinsic functional connectivity analyses. Note that task-related parameter estimates were extracted for the final information processing analysis, as described later.
Structural data processing and segmentation
Structural template calculation (T1-weighted) and segmentation
To examine and illustrate group-level results later on, a group specific T1-weighted template was calculated using ANTS buildtemplateparallel.sh (Avants et al., 2010). For illustration purposes and to aid group analyses, in addition, the T1 template was manually segmented into subregions subiculum and CA1 within the hippocampal body with ITK-SNAP (Yushkevich et al., 2006) based on the segmentation rules described in Berron et al., 2017. The first slice in each hemisphere that did not contain the uncus anymore, served as start of the hippocampal body in all hippocampal subregions. The last segmented slice was the one at which both, the inferior and superior colliculi had completely disappeared, applied for each hemisphere separately. Moreover, to evaluate results across the transversal sub/CA1 axis, the subiculum masks in each hemisphere were cut in five equally wide segments from medial to lateral within each coronal image. As the CA1 region gets more and more tilted towards the hippocampal tail, the three transversal CA1 segments were determined based on manual segmentation following a geometrical rule. Therefore, the two outer borders along the transversal axis of CA1 were connected with a line. From the middle point of that line, two straight lines were drawn in a 60° angle to determine roughly equally sized transversal CA1 segments within each coronal slice and hemisphere (a figure displaying the cuts, the procedure for the CA1 segments and the numbers of voxels within each segment can be found in appendix 9). Related to the overall size of the subregions, we opted to build five subiculum and three CA1 segments along the transversal axis from proximal to distal ends.
Segmentation of individual regions of interest
We manually segmented regions of interest (ROI) in the medial temporal lobe according to the segmentation protocol by Berron et al., 2017. Based on individual T2-weighted images, the parahippocampal cortex, Area 35, Area 36 and the EC are delineated (see appendix 8 for quality assurance measures). Moreover, we ran a Freesurfer 6.0 segmentation on the group T1 template to segment the isthmus cingulate cortex as retrosplenial mask (Desikan et al., 2006; Fischl, 2012). Note here that Syversen et al., 2021 used a similar region, however excluded the most superior part. For individual retrosplenial masks, the obtained mask was co-registered from the group T1 template space to the individual T1 space by making use of the alignment matrices obtained during above described T1 group template calculation (see appendix 7 for co-registration procedure and alignment assessment). For this alignment process we used ANTS WarpImageMultiTransform.sh (Avants et al., 2011). The retrosplenial, parahippocampal and perirhinal Area 36 and Area 35 regions served as cortical source regions for an initial functional connectivity analysis that we conducted to obtain functional subregions within the entorhinal cortex (see upcoming paragraphs and appendix 2).
Co-registration of individual structural data to functional data space
For later functional data extraction, the individual T1-weighted and T2-weighted structural images were co-registered and resliced to the echo-planar images. Therefore, ANTS was used to transfer the T2-weighted structural image to the participant’s T1 space (Avants et al., 2011). For the co-registration between individual T1-weighted and EPIs, FSL epi_reg was applied (Jenkinson and Smith, 2001). All subsequently segmented individual masks were co-registered to the participant’s functional (echo-planar) images using the obtained warping matrices (see appendix 7 for co-registration procedure and alignment assessment). ANTS WarpImageMultiTransform.sh was applied for T2 to T1 co-registration and FSL flirt was used for T1 to echo-planar image co-registration (Avants et al., 2011; Jenkinson and Smith, 2001).
ROI preparation for seed regions in functional connectivity analyses
All masks that served as source and seed regions throughout the functional connectivity analyses (retrosplenial, parahippocampal, perirhinal Area 36 and Area 35 and the later defined entorhinal subregions) were thresholded according to mean intensity to prevent signal dropout and thus a distortion of the average functional signal extracted from seed regions for the connectivity analysis. Therefore, we followed Libby et al., 2012 and Maass et al., 2015, to remove all voxels from each ROI that showed a mean intensity over time of less than two standard deviations from the mean intensity across all voxels. The thresholding was performed before each seed-to-voxel functional connectivity analysis.
Functional connectivity analyses at the participant level
Two different functional connectivity analyses were performed that build upon the approach by Maass et al., 2015. The first analysis served to identify functional subregions (‘seeds’) within the entorhinal cortex that uniquely connect with functionally and clinically relevant cortical sources. The second, core analysis, then evaluated the intrinsic functional connectivity pattern between these entorhinal seeds and hippocampal subiculum and CA1. Both functional connectivity analyses were performed on residuals of task-related functional data, creating a dataset that resembles resting-state data (Gavrilescu et al., 2008; Maass et al., 2015). In the following we describe the analysis procedure in detail. Note, that all analyses were conducted independently in both hemispheres.
To determine functional entorhinal seed regions we first performed a seed-to-voxel semipartial correlation analysis (Whitfield-Gabrieli and Nieto-Castanon, 2012) between the individually extracted residuals from retrosplenial, parahippocampal and perirhinal Area 36 and Area 35 sources as well as entorhinal voxels. The regions we call cortical sources served as seeds in that analysis. Note that the semipartial correlations calculate the variance in a voxel that is uniquely explained by the source, excluding contributions from other sources. Please refer to appendix 6 for more details on this functional semipartial correlation analysis. To obtain entorhinal seeds for the core functional connectivity analysis, first, we calculated one-sample T-tests across participants on the individually obtained and aligned, standardized beta maps for each source, respectively. Second, the four resulting statistical maps (one for each source) have been thresholded at T>3.1. Each entorhinal voxel now was attributed to be preferentially connected with one of the four source regions, based on the voxel’s maximum T value across the thresholded one-sample T-test maps. Those voxels that did not reach the threshold of T>3.1 in any of the four statistical maps have not been attributed to be preferentially connected with any of the four cortical sources. Finally, across hemispheres we selected for each source preference an equal number of these highest preference voxels across all T-tests (the number is determined by the hemisphere with the lowest relevant number of voxels). This procedure yielded four entorhinal subregions, one containing the entorhinal voxels that preferentially functionally connect with the retrosplenial (1530 voxels), one containing the entorhinal voxels that preferentially functionally connect with the parahippocampal cortex (145 voxels) and one each that contained the preferentially functionally connected voxels with perirhinal Area 35 (298 voxels) and Area 35 (751 voxels), respectively. All four entorhinal seed masks were determined on group level and co-registered to each participant. They then served as seed regions for the core functional connectivity analysis between entorhinal cortex seeds and hippocampal subregions.
For the core functional connectivity analysis (entorhinal seeds-to-hippocampal subregion voxels), an analogous seed-to-voxel semipartial correlation analysis was performed on the individual residual functional imaging data using the CONN toolbox (Whitfield-Gabrieli and Nieto-Castanon, 2012). Note again that the semipartial correlations calculate the variance in a voxel that is uniquely explained by the seed, excluding contributions from other seeds. Now the four entorhinal subregions served as seeds and functional connectivity was examined with the whole brain (later masked by the hippocampal subregion masks). For each functional connectivity analysis, mean time series were extracted from the respective seed region and entered as regressor of interest. White matter and CSF time series, realignment parameters and outliers served as regressor of no interest. The functional data from the residuals was band-pass filtered (0.01–0.1 Hz) and semipartial correlations were obtained between the seed timeseries and all other brain voxel’s timeseries. The obtained beta maps contained Fisher-transformed correlation coefficients and were used for subsequent group analyses.
Alignment between participants
To be able to perform group statistics on the resulting topography in the beta maps, the individual data was aligned to group space. Here, the T1 template image served as reference space. Using the inverse of the previously obtained individual warping matrices from individual T1 to EPI, first the standardized beta maps were co-registered from EPI to individual T1 space. In a further step, the statistical maps were then aligned between the individual T1 space and the group T1 template space, by making use of the alignment matrices obtained during above described T1 group template calculation. For this alignment process we used ANTS WarpImageMultiTransform.sh (Avants et al., 2011).
Functional connectivity analysis at group level
To investigate the functional connectivity profile between the four entorhinal seeds and the subiculum and CA1 subregion across individuals, we evaluated connectivity preferences to either seed within all transversal segments of the subiculum and CA1 target regions. Therefore, mean values for connectivity estimates to either entorhinal cortex seed were extracted from the group aligned but participant-specific beta maps out of each transversal segment, averaged along all coronal slices. Note, that segment-based extraction is necessary due to the varying number of sagittal slices that cover the respective regions along the longitudinal axis of the hippocampal body. Based on these participant-level connectivity results, connectivity preference plots for all four entorhinal seeds have been created to depict tendencies along the transversal sub/CA1 axis.
A hierarchical repeated-measures ANOVA testing procedure was employed to reveal significant differences in the transversal hippocampal connectivity patterns between entorhinal seed regions. Therefore, in a first step, an overall repeated measures ANOVA (4 seed X transversal segments) was performed per target region (subiculum and CA1 in both hemispheres) to reveal whether significant differences in seed connectivity estimates exist across the transversal axis of the respective target region (subiculum or CA1). If the overall seed X transversal segment interaction effect was significant (false-discovery-rate corrected according to Benjamini and Hochberg, 1995), in a second step, one-way repeated measures ANOVAs have been performed for each seed to identify those entorhinal seeds that indeed show a differential connectivity pattern across the transversal axis of the target region (all false-discovery-rate corrected according to Benjamini and Hochberg, 1995). If more than one seed main effect was significant, finally we determined whether these seeds exhibit an opposing connectivity pattern across the transversal axis of subiculum or CA1, respectively, by evaluating the pair-wise seed X transversal segment interaction effects on the extracted connectivity estimates.
For a more detailed topographical display of the entorhinal-hippocampal connectivity results, we calculated one-sample t-tests on the aligned, standardized beta maps that we obtained in the first-level analyses for each seed respectively. Crucially, the resulting group-level one-sample T-test statistical maps were only used to display results but not for any further statistical inference. To depict the topography of the respective voxel-wise seed preferences, the resulting group-level T-maps were thresholded with T > |0.001| and masked with the respective subregion of interest. To depict general tendencies in the connectivity profile, for each voxel in the region of interest the preferred seed connectivity was determined by attributing it to the seed with the highest T value across the one-sample T-test maps. The resulting maps were depicted in 3D plots, generated with ITK-SNAP (Yushkevich et al., 2006) that provide an overview of each voxel’s preference for the respective seed functional connectivity at a glance.
Functional analysis of content-related activity at participant level
To investigate whether scene and object information is differentially processed within entorhinal seed regions and along the transversal sub/CA1 axis, the results from the initially fitted general linear models (used to remove task effects) were examined. Contrast estimates were calculated between the beta estimates obtained from task conditions in which individuals saw indoor rooms (scene) versus objects on the screen and conditions in which individuals saw the scrambled stimuli (baseline). The resulting contrast value maps for object > baseline and scene > baseline were then co-registered to the T1 group template space. Subsequently, individual mean contrast estimates have been extracted from the four entorhinal seed regions and from those transversal segments that had previously been used for the evaluation of the intrinsic functional connectivity results (three or five segments in CA1 and subiculum, respectively).
With repeated measures ANOVAs (content condition X entorhinal region or content condition X transversal hippocampal segment), we investigated whether contrast estimates differed under scene and object conditions in the respective regions. Effect of interest thus, was the interaction between the content condition and the subregion or segment, respectively. Post-hoc paired-samples T-test were performed if the respective interaction effect was significant, to reveal in which subregion or segment functional activity between scene and object processing conditions differed significantly from each other (all false-discovery-rate corrected according to Benjamini and Hochberg, 1995).
Appendix 1
Left hemisphere results
Four cortical sources divide the left entorhinal cortex (EC) in retrosplenial-, parahippocampal-, Area 35- and Area 36-based seeds
Based on functional connectivity preferences to the sources parahippocampal cortex (Source code 4), retrosplenial cortex (Source code 3), Area 36 (Source code 2) and Area 35 (Source code 1), we obtained four entorhinal seeds. The majority of voxels can roughly be described as clustering in the posterior-medial entorhinal portion for the ECPHC-based, the anterior-medial (and posterior-medial) portion for the ECRSC-based seed, the anterior-lateral portion for the ECArea35-based and the posterior-lateral portion for the ECArea36-based seed (see appendix 2 for exact voxel counts). Note that both perirhinal-based entorhinal seeds extend along the anterior to posterior axis such that the ECArea35-based progresses more along deep entorhinal portions (with a main focus anteriorly) and the ECArea36-based along superficial entorhinal portions (with a main focus posteriorly, see Appendix 1—figure 1 and the medial reflection of the EC seeds).
Left distal subiculum is functionally connected with the ECPHC-based seed while the subiculum/CA1 border is connected with ECRSC-based and ECArea35-based seeds
When extracting estimates of connectivity preferences across individuals from proximal and distal hippocampal subfield segments for either entorhinal seed, repeated measures ANOVAs revealed significant seed X segments interaction effects along the transversal axis of the left subiculum and CA1 (subiculum: F(12,372) = 4.609; p<0.001; CA1: F(6,186) = 2.458; p=0.047; see Appendix 1—figure 2).
-
Appendix 1—figure 2—source data 1
- https://cdn.elifesciences.org/articles/76479/elife-76479-app1-fig2-data1-v3.zip
In the left subiculum, additional repeated measures ANOVAs showed that the ECArea35-based (F(4,124) = 4.489; pFDR = 0.025), and ECPHC-based (F(4,124) = 8.701; pFDR <0.001) seeds displayed a significant main effect across the transversal subiculum axis. Here, the transversal preference to the ECRSC-based entorhinal seed does not survive FDR correction (F(4,124) = 4.489; Huynh-Field uncorrected p=.05), shows however the same tendency as in the right hemisphere. The differential functional connectivity preferences for the ECArea35-based and ECPHC-based seed interacted significantly across the transversal axis, as shown in a subsequent repeated measures ANOVA (F(4,124) = 10.795; pFDR <0.001).
In the left CA1, additional repeated measures ANOVAs showed that the connectivity preference towards the ECRSC-based seed displayed a significant main effect across the transversal CA1 axis (F(2,62) = 6.753; p=0.024). In the distal CA1, the preferential functional connectivity with the ECPHC-based seed was higher than in the proximal portion of CA1. In the left CA1 a similar but weaker transversal pattern was observed for connectivity preferences with the ECArea36-based (F(2,62) = 3.841; pFDR = 0.051) and ECPHC-based seed regions (F(2,62) = 3.468; pFDR = 0.051).
Left distal subiculum and ECPHC-based exhibit higher functional activity in the scene condition while other subregions show no significant difference between conditions
For the characteristics of information processing, we first focus on the left entorhinal seed regions. When extracting task-related parameter estimates for object and scene information conditions, a repeated measures ANOVA showed a significant interaction between region and information type (object versus scene; F(3,93) = 9.772; p<0.001). Post-hoc t-tests revealed that only in the ECPHC-based seed region functional activity for scene information was significantly higher than for object information (pFDR = 0.003), while in the remaining three left entorhinal seed regions no significant difference between object and scene conditions existed (all pFDR = 0.5776; see Appendix 1—figure 3).
-
Appendix 1—figure 3—source data 1
- https://cdn.elifesciences.org/articles/76479/elife-76479-app1-fig3-data1-v3.zip
In the left hippocampal subregions, extracting the task-related parameter estimates for object and scene conditions from proximal and distal segments within each participant showed a significant interaction between transversal segments and information type in the subiculum (F(4,124) = 7.697; p<0.001), not however in CA1 (F(2,62) = 1.1925; p = 0.3042) as revealed by repeated measures ANOVAs. Post-hoc T-tests showed that only in the distal subiculum segments and in the middle segment significantly more scene than object information was processed (from most distal to middle segment pFDR <0.001; pFDR = 0.0015; pFDR = 0.0274). In all other segments along the transversal axis, no significant difference in functional activity related to object and scene conditions existed (from medial to proximal: pFDR = 0.1009; pFDR = 0.2435; see Appendix 1—figure 4).
-
Appendix 1—figure 4—source data 1
- https://cdn.elifesciences.org/articles/76479/elife-76479-app1-fig4-data1-v3.zip
Appendix 2
Quantitative assessment of entorhinal seeds
To assess the main location of each cortical source preferences within the EC, we cut the left and right EC in four quadrants. This was performed in T1 template space. First, the middle slice of all coronal slices which capture the EC was determined separately for each hemisphere. This slice was used to cut the EC in quadrants I, III and II, IV. Second, the middle slice of all axial slices which capture the EC was determined. This slice served to cut the EC in quadrants I, II and III, IV (see Appendix 2—figure 1). Note, to determine the most superior axial slice, the most posterior coronal level of the EC was used. Subsequently, we counted the number of voxels that have been assigned to each of the four cortical source regions after the initial functional connectivity analyses (that served to determined EC seeds). Averaged across hemispheres, most voxels assigned to the retrosplenial source are in EC quadrant I, most voxels assigned to the Area 35 source in EC quadrant II, most voxels assigned to the parahippocampal cortex in EC quadrant III and most voxels assigned to Area 36 in EC quadrant IV (see Appendix 2—table 1 for detailed voxel counts). Note that these quadrants do not refer to anatomically defined EC subregions.
Appendix 3
Functional connectivity gradients by source and seed region
Appendix 4
Superior and inferior view on voxel-wise functional connectivity preferences to entorhinal seeds
Appendix 5
Object and scene processing in cortical source regions
To examine whether lower parameter estimates for object processing could be due to increased noise in this condition, we evaluated object and scene processing in the four cortical source regions. Therefore, we extracted parameter estimates for the object versus baseline and the scene versus baseline contrast from the retrosplenial and parahippocampal cortex and from perirhinal Area 36 and Area 35, respectively. All parameter estimates were extracted from the previously segmented regions of interests, coregistered to the individual EPI space.
Repeated-measures ANOVAs in both hemispheres showed a significant interaction effect between condition and region (right: F(3,93) = 60.4229; p<0.001; left: F(3,93) = 47.3421; p<0.001). Subsequent paired-samples T-tests show significantly more functional activity in the object than scene condition in Area 36 (bilateral: pFDR < 0.001) and the left Area 35 (pFDR= 0.0011). No significant difference between object and scene conditions is observed in the right Area 35 (pFDR = 0.9821). There is a significant effect of more functional activity in the scene than object condition in the parahippocampal (bilateral: pFDR<0.001) and retrosplenial cortex (bilateral: pFDR<0.001, see Appendix 5—figure 1).
-
Appendix 5—figure 1—source data 1
- https://cdn.elifesciences.org/articles/76479/elife-76479-app5-fig1-data1-v3.zip
The increased object processing in adjacent cortical source regions indicates that noise differences across conditions are not likely to cause the lack of increased object processing within entorhinal seed regions and hippocampal subregions.
Appendix 6
Functional connectivity analysis to determine entorhinal seeds
Before performing the core functional connectivity analysis between entorhinal seeds and hippocampal voxels, we had to determine the entorhinal seeds, that is, the functional subregions of the EC. We largely followed Maass et al., 2015 approach to assure comparability of results. The seeds were determined based on their functional connectivity with functionally relevant sources from the cortical object and scene information processing streams, that are the perirhinal Area 35 and Area 36, the parahippocampal cortex and the retrosplenial cortex (see Nilssen et al., 2019).
The CONN toolbox (Whitfield-Gabrieli and Nieto-Castanon, 2012) was applied to perform a seed-to-voxel semipartial correlation analysis on the residual fMRI data between the retrosplenial, parahippocampal, Area 35 and Area 36 sources and the voxels within the segmented EC mask of each individual (see the description of the core functional connectivity analysis for the precise parameters). The resulting z-transformed correlation maps were then aligned for each participant to the group template T1 space and subjected to four one-sample T-tests (one for each source preference map) to reveal significant clusters of entorhinal connectivity preferences per source across all other entorhinal seeds, respectively. The functional subregions in the EC that we identify on group level generally overlap for the preferences towards the perirhinal cortex (Area 35 and Area 36) and towards the parahippocampal cortex with the findings by Maass et al., 2015. The exact procedure to determine the entorhinal seeds for further analysis is described in the main article.
Appendix 7
Co-registration procedure and alignment assessment
Appendix 8
Quality assurance measures of manually segmented regions-of-interest
The individual regions of interest were segmented by the same two experienced raters that also segmented a subsample of our data (24 hemispheres of 22 participants) for a previous publication (Berron et al., 2017). Quality assurance measures were calculated for that subsample. Regarding intra-rater reliability, the dice similarity coefficients are above 0.88 for all segmented regions (region-specific means (SD) are as follows: PHC 0.93 (0.03); Area 36 0.91 (0.02); Area 35 0.88 (0.02); EC 0.91 (0.01)). The intraclass-correlation coefficients for intra-rater reliability are all above 0.95 (PHC 0.99; Area 36 0.96; Area 35 0.97; EC 0.98). For the inter-rater reliability, dice similarity coefficients are above 0.84 for all segmented regions (region-specific means (SD) are as follows: PHC 0.86 (0.12); Area 36 0.91 (0.02); Area 35 0.84 (0.05); EC 0.87 (0.02)). The intraclass-correlation coefficients for inter-rater reliability are all above 0.78 (PHC 0.94; Area 36 0.88; Area 35 0.87; EC 0.94; see Berron et al., 2017).
Appendix 9
Metrics for transversal subiculum and CA1 segments
Transversal subiculum and CA1 segments were cut on the group template T1 images. The average number of voxels contained in each subiculum segment was 460.8 voxels for the left subiculum (standard deviation 104.36) and 458 voxels for the right subiculum (standard deviation 75.09). For the left CA1 the average equals 360 voxels (standard deviation 27.58) and 335 voxels for the right CA1 segments (standard deviation 3.56, see Appendix 9—table 1 for segment-specific values and Appendix 9—figure 1 for an illustration).
Data availability
Source data that contain numerical data used to generate Figure 2, Figure 3, Figure 4, Appendix 1 Figure 2, Appendix 1 Figure 3, Appendix 1 Figure 4, Appendix 5 Figure 1 as well as group-level statistical maps (referred to as Source Code 1-16) that underlie Figure 1, Figure 2, Appendix 1 Figure 1, Appendix 1 Figure 2, Appendix 3 Figure 1, Appendix 3 Figure 2 and Appendix 4 Figure 1 have been provided under: Open Science Framework. ID 9v3qp. https://osf.io/9v3qp.
-
Open Science FrameworkID 9v3qp. Source Data from Functional Connectivity and Information Processing in the Entorhinal-Hippocampal Circuitry.
References
-
Controlling the false discovery rate: a practical and powerful approach to multiple testingJournal of the Royal Statistical Society 57:289–300.https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
-
Higher CSF tau levels are related to hippocampal hyperactivity and object mnemonic discrimination in older adultsThe Journal of Neuroscience 39:8788–8797.https://doi.org/10.1523/JNEUROSCI.1279-19.2019
-
Neuropathological stageing of alzheimer-related changesActa Neuropathologica 82:239–259.https://doi.org/10.1007/BF00308809
-
Staging of alzheimer’s disease-related neurofibrillary changesNeurobiology of Aging 16:271–278.https://doi.org/10.1016/0197-4580(95)00021-6
-
Integration and differentiation of hippocampal memory tracesNeuroscience and Biobehavioral Reviews 118:196–208.https://doi.org/10.1016/j.neubiorev.2020.07.024
-
The parahippocampal region: corticocortical connectivityAnnals of the New York Academy of Sciences 911:25–42.https://doi.org/10.1111/j.1749-6632.2000.tb06717.x
-
Integration of objects and space in perception and memoryNature Neuroscience 20:1493–1503.https://doi.org/10.1038/nn.4657
-
The pre/parasubiculum: a hippocampal hub for scene-based cognition?Current Opinion in Behavioral Sciences 17:34–40.https://doi.org/10.1016/j.cobeha.2017.06.001
-
Differentiable processing of objects, associations, and scenes within the hippocampusThe Journal of Neuroscience 38:8146–8159.https://doi.org/10.1523/JNEUROSCI.0263-18.2018
-
BookSpatial and nonspatial representations in the lateral entorhinal cortexIn: Derdikman D, Knierim JJ, editors. Space,Time and Memory in the Hippocampal Formation. Springer Vienna. pp. 127–152.https://doi.org/10.1007/978-3-7091-1292-2_6
-
Organization and detailed parcellation of human hippocampal head and body regions based on a combined analysis of cyto- and chemoarchitectureThe Journal of Comparative Neurology 523:2233–2253.https://doi.org/10.1002/cne.23786
-
The medial temporal lobe and recognition memoryAnnual Review of Neuroscience 30:123–152.https://doi.org/10.1146/annurev.neuro.30.051606.094328
-
The many dimensions of human hippocampal organization and (dys) functionTrends in Neurosciences 44:977–989.https://doi.org/10.1016/j.tins.2021.10.003
-
Ultra-high-field fmri reveals a role for the subiculum in scene perceptual discriminationThe Journal of Neuroscience 37:3150–3159.https://doi.org/10.1523/JNEUROSCI.3225-16.2017
-
A global optimisation method for robust affine registration of brain imagesMedical Image Analysis 5:143–156.https://doi.org/10.1016/s1361-8415(01)00036-6
-
Functional correlates of the lateral and medial entorhinal cortex: objects, path integration and local-global reference framesPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 369:20130369.https://doi.org/10.1098/rstb.2013.0369
-
Representational similarity analysis-connecting the branches of systems neuroscienceFrontiers in Systems Neuroscience 2:4.https://doi.org/10.3389/neuro.06.004.2008
-
Parallel processing streams in the hippocampusCurrent Opinion in Neurobiology 64:127–134.https://doi.org/10.1016/j.conb.2020.03.004
-
Memory retrieval along the proximodistal axis of CA1Hippocampus 26:1140–1148.https://doi.org/10.1002/hipo.22596
-
Correlation of alzheimer disease neuropathologic changes with cognitive status: a review of the literatureJournal of Neuropathology and Experimental Neurology 71:362–381.https://doi.org/10.1097/NEN.0b013e31825018f7
-
Controlling hippocampal output: the central role of subiculum in hippocampal information processingBehavioural Brain Research 174:304–312.https://doi.org/10.1016/j.bbr.2006.08.018
-
BookStatistical Parametric Mapping: The Analysis of Functional Brain ImagesElsevier.
-
Two cortical systems for memory-guided behaviourNature Reviews. Neuroscience 13:713–726.https://doi.org/10.1038/nrn3338
-
Cortico-hippocampal systems involved in memory and cognition: the PMAT frameworkProgress in Brain Research 219:45–64.https://doi.org/10.1016/bs.pbr.2015.04.001
-
Details, GIST and schema: hippocampal–neocortical interactions underlying recent and remote episodic and spatial memoryCurrent Opinion in Behavioral Sciences 17:114–123.https://doi.org/10.1016/j.cobeha.2017.07.016
-
Spatial scaffold effects in event memory and imaginationWiley Interdisciplinary Reviews. Cognitive Science 9:e1462.https://doi.org/10.1002/wcs.1462
-
Scenes facilitate associative memory and integrationLearning & Memory 26:252–261.https://doi.org/10.1101/lm.049486.119
-
The role of the human entorhinal cortex in a representational account of memoryFrontiers in Human Neuroscience 9:628.https://doi.org/10.3389/fnhum.2015.00628
-
The perirhinal cortexAnnual Review of Neuroscience 37:39–53.https://doi.org/10.1146/annurev-neuro-071013-014207
-
What’ and ‘where’ in the human brainCurrent Opinion in Neurobiology 4:157–165.https://doi.org/10.1016/0959-4388(94)90066-3
-
The anatomy of memory: an interactive overview of the parahippocampal-hippocampal networkNature Reviews. Neuroscience 10:272–282.https://doi.org/10.1038/nrn2614
-
What does the retrosplenial cortex do?Nature Reviews. Neuroscience 10:792–802.https://doi.org/10.1038/nrn2733
-
Entorhinal cortex of the monkey: V-projections to the dentate gyrus, hippocampus, and subicular complexThe Journal of Comparative Neurology 307:437–459.https://doi.org/10.1002/cne.903070308
-
Architecture of the entorhinal cortex a review of entorhinal anatomy in rodents with some comparative notesFrontiers in Systems Neuroscience 11:46.https://doi.org/10.3389/fnsys.2017.00046
-
The entorhinal cortex of the monkey: VI. organization of projections from the hippocampus, subiculum, presubiculum, and parasubiculumThe Journal of Comparative Neurology 529:828–852.https://doi.org/10.1002/cne.24983
-
Object-in-place memory predicted by anterolateral entorhinal cortex and parahippocampal cortex volume in older adultsJournal of Cognitive Neuroscience 31:711–729.https://doi.org/10.1162/jocn_a_01385
-
Investigating the functions of subregions within anterior hippocampusCortex; a Journal Devoted to the Study of the Nervous System and Behavior 73:240–256.https://doi.org/10.1016/j.cortex.2015.09.002
-
Hippocampal connectivity with retrosplenial cortex is linked to neocortical tau accumulation and memory functionThe Journal of Neuroscience 41:8839–8847.https://doi.org/10.1523/JNEUROSCI.0990-21.2021
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (Project-ID 42589994)
- Magdalena M Sauvage
- Emrah Düzel
HORIZON EUROPE Marie Sklodowska-Curie Actions (843074)
- David Berron
European Union’s Horizon 2020 Framework Programme for Research and Innovation (785907 (HBP SGA2) and 945539 (HBP SGA3))
- Emrah Düzel
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the Leibniz Institute for Neurobiology in Magdeburg for providing access to the 7 Tesla MRI Scanner. We are grateful for the support of Anne Hochkeppler and Regina Schwarzer with manual segmentations of the medial temporal lobe subregions and for insightful discussions regarding functional connectivity with Yi Chen. This work was supported by the European Union’s Horizon 2020 Research and Innovation Programme under grant agreements 785907 (HBP SGA2) and 945539 (HBP SGA3) as well as the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 42589994. DB has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 843074.
Ethics
Human subjects: Informed consent and consent to publish was obtained from human participants. The study received approval by the ethics committee of Otto-von-Guericke University, Magdeburg (Germany) under reference number 128/14.
Copyright
© 2022, Grande et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,153
- views
-
- 316
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.
-
- Cell Biology
- Neuroscience
The assembly and maintenance of neural circuits is crucial for proper brain function. Although the assembly of brain circuits has been extensively studied, much less is understood about the mechanisms controlling their maintenance as animals mature. In the olfactory system, the axons of olfactory sensory neurons (OSNs) expressing the same odor receptor converge into discrete synaptic structures of the olfactory bulb (OB) called glomeruli, forming a stereotypic odor map. The OB projection neurons, called mitral and tufted cells (M/Ts), have a single dendrite that branches into a single glomerulus, where they make synapses with OSNs. We used a genetic method to progressively eliminate the vast majority of M/T cells in early postnatal mice, and observed that the assembly of the OB bulb circuits proceeded normally. However, as the animals became adults the apical dendrite of remaining M/Ts grew multiple branches that innervated several glomeruli, and OSNs expressing single odor receptors projected their axons into multiple glomeruli, disrupting the olfactory sensory map. Moreover, ablating the M/Ts in adult animals also resulted in similar structural changes in the projections of remaining M/Ts and axons from OSNs. Interestingly, the ability of these mice to detect odors was relatively preserved despite only having 1–5% of projection neurons transmitting odorant information to the brain, and having highly disrupted circuits in the OB. These results indicate that a reduced number of projection neurons does not affect the normal assembly of the olfactory circuit, but induces structural instability of the olfactory circuitry of adult animals.