Study population and design
This study was conducted in accordance with the ethical guidelines of the Declaration of Helsinki, and use of data from the Iki-Iki Health Promotion Project (Iki-Iki study) was approved by the Ethics Committee of Hirosaki University School of Medicine (authorization number 2019-064-1). Written informed consent was obtained from all participants.
The Iki-Iki Health Promotion Project was established in 2016 as a population-based prospective study of cerebro- and cardiovascular diseases and dementia in an older Japanese population from the Iwaki area of Hirosaki City, western Aomori Prefecture, Japan. In 2016 and 2017, 2390 residents aged >64 years participated in the screening survey. Of these 2390 residents, 2226 (93.1%) underwent brain MRI. We excluded 50 participants with image distortions (7 with metal artifacts, 13 with excessive motion artifacts, 30 for whom brain volume or structural covariance intra-networks could not be measured accurately for various reasons), and 43 participants without available MRI data (43 without T1WI).
Of the 2133 participants, 218 (9.78 %) were diagnosed with MCI28 and 11 participants were diagnosed with AD by the National Institute of Neurologic and Communicative Diseases and Stroke/Alzheimer Disease and Related Disorders Association29. The remaining 1904 (89.3%) patients were CNOA. Eleven participants with AD were also excluded. Thus, 2122 participants (218 with MCI and 1904 CNOA) were enrolled in the present study (Table 1).
MRI acquisition
All brain MRI data were obtained using the same protocol on a single 3T MRI scanner (Signa EXCITE 3T; GE Healthcare, Wankesha, WI, USA) with an 8-channel brain phased-array coil. Original T1WI were acquired in the steady state using a 3D fast-spoiled gradient recalled sequence with the following parameters: repetition time, 10 ms; echo time, 4.1 ms; inversion time, 700 ms; flip angle, 10; field-of view, 24 cm; section thickness, 1.2 mm; and resolution, 1.0 × 1.0 × 1.2 mm. All images were corrected for distortion due to gradient non-linearity using Grad Warp software 30 and for intensity inhomogeneity using the “N3” function 31.
Image processing
The preprocessing of images was identical to the procedure adopted for classical voxel-based morphometry (VBM) analyses using SPM12 software (Statistical Parametric Mapping 12; Institute of Neurology, London, UK) 32,33. The structural images in native space were spatially normalized, segmented into grey matter (GM), white matter, and cerebrospinal fluid images, and modulated using the Diffeomorphic Anatomical Registration Through Exponential Lie Algebra (DARTEL) toolbox in SPM12 34. Ashburner proposed DARTEL as an alternative method of normalization in the SPM package 32. To preserve the gray and white matter volumes within each voxel, we modulated the images using Jacobian determinants derived from spatial normalization using DARTEL. In contrast to conventional VBM processing, to maintain a high spatial resolution, the voxel size was set at 1.2 mm isotropic voxel size, which is normally converted to a 1.5 mm isotropic voxel size.
Furthermore, the resulting modulated GM images were smoothed using a 3-mm full width at half maximum Gaussian kernel. After the smoothing process, we extracted the hippocampal image, defined by automated anatomical labeling35 using the WFU PickAtlas version 3.0.436,37.
To identify structural networks among hippocampal voxels, SBM analysis was performed using the GIFT toolbox (http://icatb.sourceforge.net) 14. The minimum description length principle was used to estimate the number of independent components. The minimum description length yielded 16 reliable incident components. We performed ICA using a neural network algorithm (Infomax) to minimize the mutual information of the network outputs and identify naturally grouping and maximally independent sources 38. Independent component analysis was repeated 20 times in ICASSO to ensure the stability of the estimated components.
As a result, we obtained a matrix in which the 2122 rows represented 2122 subjects (1904 CNOA and 218 MCI), and each column indicated a voxel. This matrix was decomposed into two matrices using ICA. The first matrix, called the “mixing matrix,” comprises one subject per row and IC per column. The mixing matrix involved loading coefficients that demonstrated how each structural component contributed to the 2122 participants and thus contained information about the relationship between each participant and each component. In other words, the loading coefficients reflected the contribution of each participant to the specific brain network components. The second matrix was named the “source matrix” and specified the relationship between the ICs and the voxels.
To visualize the independent components, the source matrix was reshaped back to a 3D image, scaled to unit standard deviations (Z maps), and the threshold was set at Z > 2.0.
ICV, total GMV, total WMV, and bilateral HV were calculated using the 40-brain LPBA40 atlas 39. We then calculated the GMV/ICV, WMV/ICV, and HV/ICV ratios as indicators of hippocampal atrophy. The analysis was conducted using the CAT12 toolbox (C. Gaser, Structural Brain Mapping Group, Jena University Hospital, Jena, Germany; http://dbm.neuro.uni-jena.de/cat/) implemented within SPM12 software 32,33.
Statistical analyses
All statistical analyses were performed using EZR (Saitama Medical Center, Jichi Medical University, Saitama, Japan) 40. To compare the demographic characteristics between patients with MCI and CNOA, a Mann-Whitney U test was performed to assess differences in age. The chi-squared test was used for sex comparisons. Nominal variables were expressed as percentages and continuous variables as medians (interquartile ranges) or ranges based on distribution.
Logistic regression analysis was used to investigate the association of cognitive status (MCI and CNOA) with hippocampal intra-network connectivity and the HV/ICV ratio. Logistic regression analysis was also used to assess whether hippocampal intra-networks could predict MCI.
Analyses were adjusted for age, sex, education level (less than high school, high school or equivalent, college or graduate, or professional school), and self-reported medical history (hypertension, hyperlipidemia, or diabetes). The diagnosis group (MCI and CNOA) was entered as an independent variable, and all loading coefficients were calculated to indicate hippocampal intra-network connectivity. To confirm whether the loading coefficients of networks can be independent of hippocampal volume for the prediction of MCI, we used not only the loading coefficients of networks but also the HV/ICV ratio as predictive variables (Model 2 in Table 2).