Functional connectivity abnormalities of the long-axis hippocampal subregions in schizophrenia during episodic memory

Past evidence suggests that hippocampal subregions, namely the anterior and posterior parts, may be engaged in distinct networks underlying the memory functions which may be altered in patients with schizophrenia. However, of the very few studies that have investigated the hippocampal longitudinal axis subdivisions functional connectivity in patients with schizophrenia, the majority was based on resting-state data, and yet, none aimed to examine these during an episodic memory task. A total of 41 patients with schizophrenia and 45 healthy controls were recruited for a magnetic resonance imaging protocol in which they performed an explicit memory task. Seed-based functional connectivity analysis was employed to assess connectivity abnormalities between hippocampal subregions and voxel-wise connectivity targets in patients with schizophrenia. We observed a significantly reduced connectivity between the posterior hippocampus and regions from the default mode network, but increased connectivity with the primary visual cortex, in patients with schizophrenia compared to healthy subjects. Increased connectivity between the anterior hippocampus and anterior temporal regions also characterized patients with schizophrenia. In the current study, we provided evidence and support for studying hippocampal subdivisions along the longitudinal axis in schizophrenia. Our results suggest that the abnormalities in hippocampal subregions functional connectivity reflect deficits in episodic memory that may be implicated in the pathophysiology of schizophrenia.


INTRODUCTION
Schizophrenia is a severe psychiatric disorder characterized by positive (e.g., delusions and hallucinations) and negative (e.g., amotivation, asociality) symptoms experienced by patients. This psychopathology is further associated with significant deficits in multiple cognitive domains including attention, executive functions, social cognition, processing speed, verbal and visual memory, as well as working memory [1][2][3] . Cognitive deficits have been observed early in the course of the illness, with cognitive scores of first-episode psychosis patients ranging from 1.0 to 1.5 standard deviations below the population average 3,4 . Importantly, there is vast literature establishing cognitive deficits to be stronger predictors of social and occupational dysfunction in schizophrenia than positive and negative symptoms 5,6 . Of these impairments, episodic memory deficits are among the most prominent and consequential for the functioning of individuals with schizophrenia 3,5 . This raises the need to better understand the pathophysiology of episodic memory deficits in schizophrenia.
The hippocampus is crucial for encoding and retrieving the details (e.g., context) of personal events 7 . In patients with schizophrenia, evidence from meta-analyses and literature reviews suggest significant hippocampal volume reduction when compared to healthy subjects [8][9][10][11] . The hippocampus is thought to be a crucial biomarker in the pathophysiology of schizophrenia 12 , given that structural deficits in this region (i) have been observed in the first episode of psychosis 13,14 , (ii) appear to predict progression to psychosis in individuals clinically at-risk 15 , and (iii) seem to further deteriorate after illness onset 16,17 . Noteworthy, several structural neuroimaging studies have shown volume reductions in the hippocampus to be associated with decreased episodic memory performance in schizophrenia 18,19 .
Hippocampal deficits in individuals with schizophrenia have also been observed in functional neuroimaging studies during resting-state activity 20 , emotional perception, and experience 21,22 , facial emotion processing 23 , as well as during episodic memory tasks [24][25][26][27][28][29] . Although less extensively studied, hippocampus-based functional connectivity studies have demonstrated abnormalities in schizophrenia patients when compared to healthy subjects. For instance, resting-state connectivity alterations were observed between the hippocampus and regions of the default mode network (DMN), including the medial prefrontal cortex, PCC, and the precuneus [30][31][32][33][34][35] . Moreover, during working memory tasks, connectivity alterations have been observed in patients with schizophrenia between the hippocampus and the PCC 36,37 , and executive regions such as the dorsolateral prefrontal cortex 38 and the inferior frontal gyrus 39 . To our knowledge, no study has employed an episodic memory task to examine the connectivity patterns of hippocampus subregions in patients with schizophrenia.
In the past decades, evidence from rodents, primates, and humans has supported the segmentation of the hippocampus along the longitudinal axis (i.e., anterior-posterior, ventral-dorsal, and temporal-septal). For instance, it has been shown that these two sub-regions differ substantially in their anatomical and functional connections. In fact, results suggest that the anterior subregion is preferentially connected to the perirhinal cortex, the lateral temporal cortex extending to the temporal pole, the amygdala, striatum, the pre-supplementary motor area, while the posterior portion is preferentially connected to the posterior parahippocampal cortex, the pregenual anterior and posterior cingulate cortices, precuneus, thalamus, inferior parietal lobe, and the cuneus/lingual gyrus [40][41][42][43][44][45] . However, differences in functional connectivity between both subregions with frontal regions (e.g., vmPFC and dmPFC) remain unclear [40][41][42][44][45][46] . These hippocampal subregions have long been suggested to be involved in memory encoding and retrieval respectively 47,48 . As the posterior hippocampus may also be implicated in encoding 43,[49][50][51] , the evidence in support of this dichotomy remains equivocal.
Only a few studies have examined the functional connectivity of hippocampal subregions in schizophrenia. In resting-state connectivity, the investigations performed thus far have shown deficits in functional connectivity mainly between the anterior hippocampus (aHippocampus) and DMN regions, such as the medial prefrontal cortex and the PCC 31,35,52,53 , as well as between the posterior hippocampus (pHippocampus) and the anterior cingulate cortex 31,35,52,53 . Alterations of connectivity between the aHippocampus and temporal regions 35,53,54 , as well as between the pHippocampus and the dorsolateral prefrontal cortex 54 , have also been observed, though less consistently. To our knowledge, only one task-based functional connectivity study has been published on hippocampal subregions in schizophrenia showing reduced connectivity from the posterior hippocampus to the inferior frontal gyrus during a complex working memory task 39 . It is striking that no research on functional connectivity of hippocampal subregions has been performed in schizophrenia using an episodic memory task, despite the hippocampus being well-known to play a key role in this cognitive domain 7 .
In view of the extant literature, the primary objective of the current study was to examine the functional connectivity of hippocampal subregions in schizophrenia during an episodic memory task. As a secondary objective, we sought to replicate the structural deficits that have been reported in these subregions in schizophrenia [55][56][57][58][59] . Finally, given that emotional experience may alter episodic memory performance in schizophrenia 60 and that the aHippocampus is more strongly connected to the limbic system than the posterior hippocampus 43 , we employed an episodic memory task with an emotional component in order to explore potential interactions between cognition and emotion.

RESULTS
Sociodemographic, clinical, and task performance Healthy controls (HC) and schizophrenia patients (SZ) subjects did not differ in terms of age or gender. The SZ group attained significantly lower scores on the bloc design subset of the WAIS-III than the HC group, as well as SZ, reported lower education levels than HC (Table 1). Finally, a poorer performance was observed in SZ participants in memory accuracy during the retrieval condition of the functional magnetic resonance imaging (fMRI) task, across negative, positive, and neutral images (approximatively 25% difference between groups). Clinical characteristics of individuals with SZ can be found in Table 1. Furthermore, a significant between-group difference was observed on a number of removed volumes during the encoding task (p = 0.02) but not during retrieval conditions (p > 0.05).
Weighted seed-based connectivity analyses First, when examining the connectivity patterns of the anterior and posterior hippocampus between HC and SZ during the retrieval condition (e.g., across a block of Run-1) compared to the encoding condition (e.g., across a block of Run-2) (i.e., Diagnosis [HC versus SZ]*Condition [Retrieval versus Encoding]*Seeds [Hemisphere specific anterior versus posterior hippocampus] interaction], no significant difference was observed. We further examined the different connectivity patterns of the anterior and posterior hippocampus between HC and SZ (i.e., Diagnosis*Seeds) interaction within each condition separately.
In the Encoding condition, the Diagnosis by Left Hemisphere Seeds (aHippocampus versus pHippocampus) revealed significant connectivity differences with the intracalcarine cortex (ICC) and the dorsomedial prefrontal cortex (dmPFC). However, adding the number of removed volumes (motion scrubbing) as covariate only replicated the dmPFC (See Table 2 Table 2). Post hoc analyses showed reduced pHippocampus-dmPFC connectivity (F 1,84 = 30.00, p < 0.001), in SZ compared to HC. In within-group comparisons, we observed that pHippocampus-dmPFC connectivity was significantly reduced compared to aHippocampus-dmPFC in SZ (p < 0.001), but not in HC (p = 0.160). No significant target was observed for the diagnosis by right hemisphere seeds (anterior versus posterior subregions) during Encoding.
In the Retrieval condition, the diagnosis by left hemisphere seeds interaction resulted in significant connectivity differences in the precuneus/posterior cingulate cortex (PCC), the ventromedial prefrontal cortex (vmPFC), and a cluster that includes the supramarginal gyrus (SMG) and voxels spanning the inferior parietal lobule (see Table 2 and Fig. 1). Adding the subject's average of ART's composite motion measure on the remaining volumes as covariate replicated these results (Supplementary Table 2). Post hoc analyses revealed decreased connectivity between the left pHippocampus with the precuneus/PCC (F 1,84 = 10.06, p < 0.001), while weaker increased connectivity was observed between the left aHippocampus and the precuneus/ PCC (F 1,84 = 5.23, p = 0.025), in SZ compared to HC (see Fig. 2B). Moreover, we observed that the pHippocampus-vmPFC connectivity (F 1,84 = 20.69, p < 0.001) and the pHippocampus-SMG (F 1,84 = 16.48, p < 0.001) were reduced in SZ compared to HC (see Fig. 2).
Diagnosis by right hemisphere seeds during the retrieval condition suggested significant connectivity differences with the inferior temporal gyrus (ITG) and the ICC (see Table 2, Figs. 1 and 2C). These differences were also observed after adding the subject's average of ART's composite motion measure on the remaining volumes as covariate (Supplementary Table 2). Post hoc analyses revealed increased connectivity between the aHippocampus and the ITG (F 1,84 = 14.57, p < 0.001) in SZ compared to HC. Furthermore, increased connectivity between the pHippocampus and the ICC (F 1,84 = 12.94, p = 0.001) was observed in SZ compared to HC.  Within-group comparisons suggested that aHippocampus-ITG connectivity was significantly greater than pHippocampus-ITG in SZ (p < 0.001), but not in HC (p = 0.456) and that pHippocampus-ICC connectivity was greater than aHippocampus-ICC in SZ (p < 0.001), but not in HC (p = 0.412). Finally, the diagnosis by stimuli valence (i.e., POS, NEG, and NEU) interaction analysis revealed that none of the main connectivity findings were significantly altered by stimuli valence. Correlational analyses suggested that no clinical variables, including the level to antipsychotic medication (i.e., chlorpromazine equivalents) were significantly associated with our connectivity results (p > 0.341). Neurobiological correlates of memory accuracy in patients with schizophrenia Bivariate Pearson correlations revealed that in the SZ group, the aHippocampus-ITG and the pHippocampus-ICC connectivity, observed in the retrieval condition, were significantly associated with memory accuracy during the retrieval condition (r = 0.348, p = 0.03 and r = −0.470, p = 0.003, respectively). These relationships were not statistically significant in healthy subjects (r = −0.002 and r = −0.114, respectively). Moreover, the block design subset of the WAIS (r = 0.425, p = 0.014; available data for 35 patients) and education level (r = 0.345, p = 0.031) was also significantly associated with memory accuracy. No other significant relationship with memory accuracy was observed (e.g., GM volume of hippocampal subregions, medication, and clinical symptomatology).

Voxel-based morphometry
Forward stepwise multiple regression revealed a final model suggesting that the aHippocampus-ITG and pHippocampus-ICC connectivity remained statistically independent predictors of memory accuracy (β = 0.354, p = 0.021 and β = 0.437, p = 0.005, respectively), while the block design subset of the WAIS was no longer statistically significant (β = 0.295, p = 0.052). Education level was not included in the final selected model by the stepwise regression. These three variables factors showed low collinearity (VIF < 2). Furthermore, the model explained 37.2% of the variance (adjusted R-squared) associated with memory accuracy in patients with schizophrenia.

DISCUSSION
In the current study, we investigated the functional connectivity of distinct subregions of the hippocampus in patients with schizophrenia during an episodic memory task. We observed decreased pHippocampus-dmPFC connectivity in schizophrenia patients, relative to controls, during the encoding condition. Furthermore, we found reduced connectivity between the pHippocampus and DMN regions including the precuneus, the SMG, and the vmPFC, but increased pHippocampus-ICC connectivity in participants with schizophrenia during the retrieval condition. We also observed increased connectivity between the aHippocampus and anterior  temporal regions in patients with schizophrenia. The altered pHippocampus-ICC and aHippocampus-ITG connectivity was correlated with memory accuracy. The functional connectivity and behavioral results were complemented by the VBM analysis, which revealed reduced GM volume in the bilateral anterior and posterior hippocampus in patients versus healthy controls.
Although the number of studies investigating the functional connectivity of the hippocampal subregions in patients with schizophrenia is limited, similar results were observed by researchers using resting-state paradigms. In fact, some authors have found that patients showed disrupted functional connectivity between the aHippocampus and anterior temporal regions 35,52,53 . The former connectivity is particularly relevant since these anterior medial temporal lobe structures are known to play an important role in episodic memory [61][62][63] . In contrast to what we observed, the aforementioned studies found that only the anterior portion of the hippocampus was significantly disconnected from the precuneus/ PCC and the vmPFC 35,52,53 . A potential explanation for this discrepancy could be that we examined functional connectivity as related to episodic memory, while previous studies examined connectivity at rest. It is also worth mentioning that our results are congruent with structural and functional neuroimaging literature in healthy volunteers which consistently reported that the posterior subregion of the hippocampus is preferentially connected (posterior > anterior) to the posterior parahippocampal gyrus as well as DMN regions, while the anterior part is preferentially associated (anterior > posterior) with the perirhinal cortex, amygdala, and extended temporal pole [40][41][42][43][44][45] . This is further supported by studies using meta-analytical connectivity modeling and a reverse inference approach which suggest a clear Emotion (e.g., face monitor/discrimination, emotional valence, threat, anxiety, fear, and neutral stimuli and faces) to Cognition (e.g., world-centric behaviors including navigation, perceptual functioning, and processing information in its environmental context) gradient along the longitudinal axis of the hippocampus (i.e., anterior-posterior axis) 50,51 .
Although it has long been assumed that anterior-posterior subregions may play specific roles in encoding and retrieval 47,64,65 , this dichotomy is still a subject of debate 50,51,61,66 . Following recent results from the BrainMap and NeuroSynth databases that suggest equally distributed encoding/retrieval processes along the hippocampal axis 50 , our results support the importance of the pHippocampus in both the encoding and retrieval conditions as well as the aHippocampal subregions during retrieval condition. Indeed, in patients with schizophrenia, deficits in connectivity between pHippocampus and anterior DMN region has been observed during the encoding condition while the pHippocampus was altered mostly with posterior DMN during retrieval condition. Past evidence suggests that the posterior DMN regions are engaged in processes supporting successful episodic retrieval but are deactivated during encoding conditions [67][68][69][70] . This phenomenon, also called the encoding-retrieval flip pattern 67 , is however less clear for anterior DMN regions (i.e., vmPFC and dmPFC) 68,71-73 . Nonetheless, our results indicate that patients with schizophrenia exhibit deficits in connectivity between the pHippocampus and posterior DMN regions that may reflect failure in the encodingretrieval flip pattern.
Interestingly, we observed a significant increased pHippocampus-ICC connectivity during retrieval condition and this was found to be strongly associated with memory accuracy in schizophrenia. It is worth mentioning that connectivity alterations between hippocampal subregions and ICC have also been observed by other authors 31,54 . The ICC, including the primary visual cortex, is known to send projections to posterior parietal regions through the dorsal visual processing stream [74][75][76] . The posterior parahippocampal gyrus thus receives inputs from these posterior parietal regions (i.e., mainly caudal inferior parietal lobule) which, in turn, projects to the pHippocampus 77-79 , justifying the role of both structures in integrating contextual information of a stimuli as well as in retrieving the initial encoding context 80,81 . The alteration between the pHippocampus-ICC is coherent with a growing body of literature showing that the activity and connectivity of occipital regions involved in the early processing of visual stimuli are impaired in schizophrenia 82,83 .
Our study has a few limitations that need to be acknowledged. First, although the sample of the current study was sufficiently powered to detect differences in functional connectivity between patients and controls during episodic memory, it may have not been sufficiently powered to examine the interaction between emotional valence and memory conditions. We also acknowledge that some of our findings (or lack of findings) may be explained by the specificities of the episodic memory task that we used. Indeed, we used an incidental encoding task whereby participants were not explicitly asked to memorize the information that was presented to them in the first run. Furthermore, the retrieval condition consisted of a recognition memory condition, rather than a free recall condition. In both cases, this may explain why we did not observe impaired connectivity between hippocampal subregions and executive regions, such as the dlPFC. Another limitation has to do with the fact that the VBM was performed using ROIs and did not examine the whole brain. Finally, we acknowledge that antipsychotic treatment may have influenced results. Previous trials have shown that some second-generation antipsychotics produce small beneficial effects on episodic memory in schizophrenia 84 . At the neural level, the effects of antipsychotics on hippocampal volumes have, however, been inconsistent and mixed 85 . Although preliminary, a few fMRI studies have shown that second-generation antipsychotics may normalize hippocampal activity and connectivity during episodic memory in schizophrenia 86,87 . It remains unclear, however, how these latter observations may relate to our own findings considering that patients were treated with second-generation antipsychotics, and displayed hippocampal dysconnectivity nevertheless. Moreover, we found no significant association between antipsychotic dosage (e.g., chlorpromazine equivalents) and hippocampal dysconnectivity.
In the current study, we examined the functional connectivity of hippocampal subregions performed in schizophrenia using an episodic memory task. We observed a complex pattern of connectivity alterations from the anterior and posterior hippocampal subregions to the default-mode, temporal and occipital regions. Our results highlight the importance of the hippocampal subregions, which are known to play a crucial role in the rodent model of psychosis 88 and the pathophysiology of schizophrenia 12 . In the future, longitudinal fMRI studies will need to be performed in order to understand the role of hippocampal subregions in the psychosis spectrum, by pursuing research in individuals a high risk for psychosis, in patients with the first episode of psychosis, as well as schizophrenia patients who are treatment-resistant. Investigations in drug-free patients, as well as in unaffected firstdegree relatives are also warranted. Attention will also need to be paid to task design characteristics.

METHODS Sample
Forty-one individuals with schizophrenia (SZ) in a stable state (i.e., no psychotic relapse within the last 2 months and no change in their antipsychotic medication within the month preceding the study) and 45 healthy controls (HC) were recruited for the purposes of this study. Schizophrenia was diagnosed according to DSM-IV-TR 88 using the Structured Clinical Interview for DSM-IV (SCID), and participants were all treated with second-generation antipsychotics. Urine drug screenings were performed. Psychiatric symptoms were evaluated with the Positive and Negative Syndrome Scale (PANSS 89 ), while IQ was estimated based on the vocabulary, similarities, and block design subtests of the Wechsler J.R. Dugré et al.
abbreviated scale of intelligence (WASI 90 ). Antipsychotic dosage was calculated using chlorpromazine equivalents 91,92 . Control participants were screened using the non-patient edition of the SCID in order to rule out any Axis-I psychiatric disorders. General exclusion criteria were age younger than 18 years or older than 45 years, past or present neurological disorder, substance use disorder (during the year preceding the study for participants; lifetime for controls), or contraindications for magnetic resonance imaging (MRI) (e.g., cardiac pacemaker, aneurysm clip).
The ethics committees of the Centre de recherche de l'Institut en Santé Mentale de Montréal and the Regroupement de Neuroimagerie du Québec approved the study protocol. In agreement with the Declaration of Helsinki, we obtained written informed consent from participants before the experiment. The ability of participants with SZ to give informed consent was established using the guidelines of the Canadian Psychiatric Association. The data are not publicly available as they contain information that could compromise research participant privacy/consent.

Experimental procedure
In the first run, participants passively viewed blocks of pictures while in the MRI machine. The stimuli were selected from the International Affective Picture System (IAPS) 93 based on normative valence and were matched for content (e.g., people, animals, and landscapes). The images presented differed in their valence, with each image category being shown in separate blocks lasting 48.5 s, resulting in three experimental conditions: positive (POS), negative (NEG), and Neutral (NEU) content. Although subjects were aware that a memory task would follow the encodingprocessing run, they were not explicitly instructed to remember the images. Instead, to ensure that participants were attentive to the images presented during encoding, they were asked to indicate with the press of a button whether they saw a person or part of a person in the picture. This task served as an incidental learning procedure, which allowed for secondary analyses of potential interactions between learning (i.e., encoding) and emotion. This encoding task was then followed by an unrelated mental rotation task lasting 15 min as a means of separating both incidental encoding and subsequent recognition memory.
The retrieval portion consisted of viewing 48.5 s blocks of emotionally positive, negative, and neutral pictures similar to the incidental encoding task. During this second run, however, 50% of the stimuli in each block originated from the encoding task (previously viewed), while 50% were novel (never viewed before). The order of presentation of stimuli was randomized. There were 16-s periods of rest separating the blocks from one another. Each block contained ten images and each block type was repeated four times. Each picture appeared for 3000 ms followed by a blank screen with a fixation point for an average of 1.75 s (ranging from 1 to 2.5 s and giving an average interstimulus interval (ISI) of 4.75 s). During this recognition memory task, participants were to determine, by pressing the correct button, which of the stimuli were old and which were new. This task allowed for primary analyses on recognition memory regardless of emotional valence, and secondary analyses on potential interactions between retrieval and emotion.

MRI acquisition parameters
Blood oxygen level-dependent (BOLD) data were acquired using a T2weighted gradient echo-planar imaging (EPI) sequence [repetition time (TR) = 3000 ms, echo time (TE) = 30 ms, flip angle = 90°, matrix 64 × 64; voxel size = 3.5 mm 3 ; 41 axial slices] on a 3.0 Tesla TRIO-TIM MRI system. The functional slices were angled parallel to the AC-PC line. An inline retrospective motion correction algorithm was employed while the EPI images were acquired. Individual high-resolution co-planar anatomical images were also acquired using a three-dimensional, spoiled gradientecho sequence (TR = 19 ms; TE = 4.92 ms; FA = 25 o ; matrix size: 256 × 256; voxel size; 1 mm 3 ; 176 sagittal slices). fMRI data preprocessing Functional images were realigned, corrected for motion artifacts with the Artifact Detection Tool 94 (ART, setting a threshold of 0.9 mm subject ART's composite motion and a global signal threshold of Z = 5) with the implemented in CONN Toolbox 95 , high-passes filtered (>0.008 Hz) and coregistered to the corresponding anatomical image. The anatomical images were segmented (into GM, white matter, and cerebrospinal fluid) and normalized to the Montreal Neurological Institute (MNI) stereotaxic space. Functional images were then normalized based on structural data, spatially smoothed with a 6 mm full-width-at-half-maximum (FWHM) 3D isotropic Gaussian kernel, and resampled to 2 mm 3 voxels. For the preprocessing, the anatomical component-based noise correction method (aCompCor strategy 96 ), was employed to remove confounding effects from the BOLD time series, such as the physiological noise originating from the white matter and cerebrospinal fluid. This method was found to increase the validity and sensitivity of analyses 97 .

Weighted seed-based connectivity analyses
Hippocampus subregions in MNI space were constructed using maximum probability maps through the Anatomy Toolbox 98 . Voxels were included in the hippocampus (also comprising CA1-4 and the dentate gyrus 99 ) if their likelihood of being located in the hippocampus was at least 60%. Voxels were excluded if they had any probability of belonging to the amygdala. These steps were taken in order to prevent false positives due to the proximity of the aHippocampus and amygdala in functional connectivity/ GM volume differences. The hippocampus maps were then binarized and sectioned into anterior (y = −11 to −21) and posterior (y = −32 to −43) subregions along the y-axis, with a gap (y = −22 to −31), to avoid any overlap 40 .
Physiological noise, realignment parameters, and movement artifacts were regressed out as confounding effects from the BOLD time-series at each voxel. The main activation effects of the conditions were also regressed out to avoid spurious connectivity due to task co-activation. The residual time-series were weighted by the appropriate hemodynamic response function-convolved regressor to derive task condition-specific time-series for weighted functional connectivity analyses 95 . In the first-level analysis, weighted seed-based connectivity maps were calculated with a weighted least-squares linear model (bivariate regression) between the time-course of each region-of-interest (ROI) seed-to-voxel, for each subject.
In the second-level analysis, between-group differences were investigated (SZ versus HC), thus allowing to search for connections that were altered in SZ. Since we were interested in examining differences in the distinct connectivity profiles of the anterior and posterior hippocampus 40 between the two groups and task conditions, we performed a group (HC versus SZ) by task (Encoding versus Retrieval) by seed (aHippocampus versus pHippocampus, each hemisphere separately) interaction analysis. A two-sample t test was performed between groups on the number of removed volumes (motion scrubbing) during the task. Healthy subjects (M = 6.74, SD = 13.5) and SZ (M = 15.4, SD = 16.76) did differ significantly from each other (t = 2.63, p = 0.01). Therefore, the number of removed volumes was added as a covariate in the second-level analysis if the twogroups significantly differed (p < 0.05). Second-level analyses were also performed using the average of the ART composite motion measure on the remaining volumes (Supplementary Table 1). Since more than 88% of total volumes remained after motion scrubbing (>481/545 volumes), for every participant, no participant was excluded. The peak threshold was set at p < 0.005 (two-sided) with a cluster threshold corrected for the false discovery rate of p < 0.05.
Finally, regression coefficients of significant targets were extracted using MARSeille Boîte À Région d'Intérêt (MarsBaR) toolbox (sourceforge.net/ projects/marsbar/). Within-group post hoc analyses were performed to assess the preferential connectivity targets between aHippocampus versus pHippocampus. In order to assess whether specific emotions may have driven functional connectivity results, we performed post hoc analyses using IBM SPSS Statistics 25 for Windows. Moreover, relationships between functional connectivity regression coefficients and clinical correlates were tested with Pearson's correlation. Imaging results were visualized with BrainNet 100 , data were displayed using ggplot2 101 , and figures were built with Inkscape (inkscape.org).

Voxel-based morphometry
Anatomical images were preprocessed using the Computational Anatomy Toolbox for SPM12 (CAT12) 102 . The CAT12 preprocessing pipeline includes image realignment and tissue segmentation into gray matter (GM), white matter, and cerebrospinal fluid. The segmented scans were normalized to a 1.5 mm isotropic predefined adult template provided by the CAT12 toolbox (in MNI space), using the Diffeomorphic Anatomic Registration Through Exponentiated Lie algebra algorithm (DARTEL). The deformation parameters were estimated using the nonlinear spatial registration of DARTEL and applied to the tissue segmentations. Modulation with the Jacobian determinant of the deformation was executed to preserve the total amount of GM signal following spatial normalization. The segmented and modulated GM images were then smoothed using a 6-mm FWHM Gaussian kernel.
Retrospective quality assurance was performed using a weighted image quality rating (IQR) based on noise, inhomogeneities, and image resolution, provided by the CAT12 toolbox. Individuals were excluded from further analysis if the IQR was at >2 standard deviations below the mean percent rating of the sample (M = 77.8%, SD = 5.94%). This resulted in the exclusion of three HC and four SZ participants.
In order to perform analysis on hippocampal subregions and their respective connectivity targets, the forward deformation DARTEL field of an ICBM152 T1 image (MNI space) was saved and applied to each binary masked ROI. Two-sample t tests were then performed on resulting images with total intracranial volume and IQR to remove their potential confounding effects on ROIs. Moreover, control over the family wise error rate through nonparametric permutation TFCE approach 103 was performed with 5000 permutations within each ROIs separately (SVC). Standard parameter values E = ½, H = 2 were used for the TFCE method. Significant differences were detected at a p < 0.05 (TFCE-SVC). Relationships between GM volume of hippocampal subregions and clinical correlates were performed with Pearson's correlations.

Neurobiological correlates of memory accuracy in patients with schizophrenia
In order to better understand the correlates of memory accuracy, as defined as percentage images well identified, in patients with schizophrenia, bivariate Pearson's correlations were first performed between neurobiological features (connectivity and VBM analyses) and memory accuracy during the retrieval condition. The statistically significant neurobiological factors were then entered in a multiple linear regression model. Clinical variables (e.g., medication, age of onset, clinical symptomatology, and IQ) associated with memory accuracy (p < 0.05) were entered as covariates in the multiple linear regression. A stepwise (forward selection) regression model was used to prevent overfitting.

Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.

DATA AVAILABILITY
The data are not publicly available as they contain information that could compromise research participant privacy/consent. The data that support the findings of this study are available upon reasonable request from the corresponding author. S.P. but are only redistributable to researchers engaged in IRB approved research collaborations.

CODE AVAILABILITY
Embedded pipelines in toolboxes were used to generate results with no custom code necessary for conclusions or replication: CONN Toolbox version 18b (https://web. conn-toolbox.org/) and CAT12 (http://www.neuro.uni-jena.de/cat/).