Establishing a baseline for human cortical folding morphological variables: a multicenter study =============================================================================================== * Fernanda Hansen Pacheco de Moraes * Victor B. B. Mello * Fernanda Tovar-Moll * Bruno Mota ## Abstract Differences in the way human cerebral cortices fold have been correlated to health, disease, development, and aging. But to obtain a deeper understating of the mechanisms that generate such differences it is useful to derive one’s morphometric variables from first principles. This work explores one such set of variables that arise naturally from a model for universal self-similar cortical folding that was validated on comparative neuroanatomical data. We aim to establish a baseline for these variables across the human lifespan using a heterogeneous compilation of cross-sectional datasets, as the first step to extend the model to incorporate the time evolution of brain morphology. We extracted the morphological features from structural MRI of 3650 subjects: 3095 healthy controls (CTL) and 555 Alzheimer’s Disease (AD) patients from 9 datasets, which were harmonized with a straightforward procedure to reduce the uncertainty due to heterogeneous acquisition and processing. The unprecedented possibility of analyzing such a large number of subjects in this framework allowed us to compare CTL and AD subjects’ lifespan trajectories, testing if AD is a form of accelerated aging at the brain structural level. After validating this baseline from development to aging, we estimate the variables’ uncertainties and show that Alzheimer’s Disease is similar to premature aging when measuring global and local degeneration. This new methodology may allow future studies to explore the structural transition between healthy and pathological aging and may be essential to generate data for the cortical folding process simulations. **Significance Statement** Understating Cortical folding is of increasing interest in neurosciences as it has been used to discriminate disease in humans while integrating pieces of knowledge from compared neuroanatomy and neuroproliferations programs. Here we propose estimating the baseline of cortical folding variables from multi-site MRI human images, evaluating the changing rate of its independent variables through the human lifespan, and proposing a simple harmonization procedure to combine multicentric datasets. Finally, we present a practical application of these techniques comparing Alzheimer’s Disease and Cognitive Unimpaired Controls based on the estimated changing rates. **Highlights** * Baseline of independent cortical folding variables from 3650 multi-site human MRI * Propose a simple harmonization procedure to combine multicentric datasets * Evaluate the changing rate of independent variables through the human lifespan * Practical application comparing Alzheimer’s Disease and Controls rates Keywords * Cortical Folding * Aging * Alzheimer’s Disease * Harmonization ## Introduction Mapping the human brain development and aging longitudinally from a unique dataset is barely impossible due to our extensive lifespan. Suppose one wants to understand the time evolution of the brain morphology through the whole lifespan. In this case, it is mandatory to combine multiple datasets and devise methods that allow a fair comparison among them. The recent advent of large heterogeneous datasets (to ensure the legit results across several populations) and data curation innovations allowed neuroscience to explore the brain’s changes on an enormous scale, from functional data to structural studies. Here, we propose a novel methodology of combining structural MRI from different acquisition sites and equipment, acknowledging their heterogeneity and providing insights about cortical folding during development and aging, with a practical application of these results in Alzheimer’s Disease, the most common dementia worldwide. This effort is essential to provide experimental information about the evolution span of the cortex morphology and build a cortical folding theory. In the past years, the study of cortical folding in humans and other mammals was extensively promoted by the application of translational science, open-access databases, and data science tools which allowed the field to grow in multiple and sometimes convergent directions [1–4]. Due to its biological background, cortical folding measurements have been included in human brain structure analysis, discriminating disease from healthy controls [5, 6], describing its correlation with cognition [7], and relation to aging [8, 9]. Despite the lack of consensus defining a unique theory that explains cortical folding in every scale, Mota & Herculano-Houzel [10] proposed a cortical folding model respected by more than 50 mammals brain hemispheres. It predicts a power-law relationship between cortical thickness (T), exposed (A*E*), and total area (A*T*) (Equation 1). ![Formula][1] Where *α* is the fractal index and universal constant, with a theoretical value of 1.25 and calculated as 1.305 [10] when considering a heterogeneous dataset of different species, thus, the usual geometrical variables used to describe the cortex are not independent in this model. Wang et al. continued this work, validating the theory for different groups of humans [11] in specific regions of interest [12], obtaining new evidence that the cortex is indeed a fractal structure. In 2021, Wang et al. showed that from this model one could derive a more natural set of nearly independent morphological variables, K, S, and I, that could be used to improve disease discrimination from regular and typical structural changes of the brain, like aging [13]. In this framework, the morphology of each cortex is expressed as a point in a three-dimensional abstract space, with each coordinate component corresponding to the log-value of an independent morphometric variable. The use of log-values guarantees that linear combinations of the basis vectors correspond to power-law relations. Conversely, as long as any new variables are expressible as power laws, different but equivalent sets of variables can thus be related to each other by a change of Cartesian coordinate bases. As a starting point, we use the log-values of the commonly morphometric variables, the total area log10 AT, exposed area log10 AE and average thickness log10 T2. We then derive our new variables thusly: *K* = *log*10k (Equation 2) is a near-invariant quantity obtained by isolating *k* in Equation 1, and is associated with conserved viscoelastic [REF] properties of cortical matter. S (Equation 3), also dimensionless, encapsulates the aspect of brain shape that changes more significantly across cortices, and I (Equation 4), representing brain isometric volume and carries the information about overall cortical size.1. Those three variables combined may help distinguish pathological events similar to age effects, as AD. ![Formula][2] ![Formula][3] ![Formula][4] Baseline values of those cortical morphological variables and their inherent uncertainties over the human lifespan were never defined due to methodological limitations of combining multicenter studies of brain structural images [14, 15]. However, both are crucial for applied studies. Uncertainties determine measuring limitations and allow a proper comparison between different cohorts. Baseline values allow clinical applications and, when necessary, simulations of healthy control data. Multi-center studies have become more common in the last few decades due to improved data sharing and the open science trend. However, combining MRI data from multiple experiments is not a trivial task. Part of these uncertainties is explained by type A errors, i.e., random, which appear as a consequence of the natural fluctuation within the species and from the data acquisition itself. Theoretically, these random errors are the same if one considers multiple experiments. Confounding components could be added when gathering MRI images with differences in acquisition parameters, acquisition equipment, and versions of the post-processing software [15, 16]. Together, these effects can add type B errors to the data, ultimately reflecting a systematic difference in morphological variables calculated from different databases and increasing the data spread. In practice, the natural fluctuations are convoluted with the type B errors, enhancing the data spread. Recently, a vast literature in primary morphological variables from MR images, [17–19], and more complex ones [14], explored the limitations of compiling images acquired with different scanners and protocols and its implication in statistical analysis results. The suggested solution for these limitations is a procedure called data harmonization, in which there is an “explicit removal of site-related effects in multi-site data” [18]. It is essential to add that one cannot, in general, untangle the contributions of natural fluctuation and data acquisition to random errors. Also, it is not possible to distinguish between type B errors of the acquisition and processing. For this, one would need a diligent procedure to track all differences between acquisition and processing as well as multiple images of the same individual for each sample. This manuscript deals only with the first case, estimating the type A errors due to acquisition and processing using three subsequent structural images of the subjects2. The uncertainties associated with the variables (K, S, and I) are summarized in Equation 5. ![Formula][5] In this manuscript, we developed a simplistic harmonization method that allows a multi-site analysis, determining K, S, and I baseline values, their uncertainties, and the ratio of changing rates through years for cross-sectional images of healthy controls. Although this work focused is these novel morphological variables, we also provided the baseline of usual morphological variables such as thickness *T* and Gyrification Index GI. A robust estimation framework of a baseline time function is important to provide future clinical studies enough information to compare brain structure trajectories in case of a pathological investigation. We then, as a clinical application and extension of de Moraes et al. [20], and Wang et al. [12], verify if Alzheimer’s Disease (AD) is similar to a premature accelerated aging brain in terms of the independent morphological variables, K, S, and I, by comparing their rates. Finally, we have determined the time evolution of the fractal dimension *α* proposed in Mota & Herculano-Houzel’s model [10], interpreting all results within the framework of this theory. ## Materials and Methods Participants and data included in this analysis were either acquired by the D’Or Institute for Research and Education (IDOR), Rio de Janeiro, RJ, Brazil, used by [20] (approved by the Hospital Copa D’Or Research Ethics Committee under protocol number CAAE 47163715.0.0000.5249), published by Wang et al. [11, 12] or from open-access databases, as AHEAD [21] and AOMIC (PIOP01 and PIOP02) [22]. The total number of subjects is 3095 healthy controls from 4 to 96 years old and 555 Alzheimer’s Disease subjects from 56 to 92 years old. To investigate the type B errors during image acquisition and variances in repeated measures performed in the same conditions, we included the first 50 subjects from the AOMIC ID1000 [22], described in Supplementary Material. AHEAD and HCPr900 only include age range instead of the actual age due to local Ethics Committee rules. To overcome the lack of a exact number, we assumed the interval’s mean age as a reference. Datasets’ demographics are summarized in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T1) Table 1: Summary for each dataset. Mean value ± standard deviation. Diagnostic: CTL – Healthy control; AD – Alzheimer’s Disease. Details in Supplementary Material, Table 4. The structural images from IDOR, AHEAD, AOMIC ID1000, PIOP1, and PIOP2 were processed in FreeSurfer v6.0.0 [27] with the longitudinal pipeline [28] without manual intervention at the surfaces [29]. The FreeSurfer localGI pipeline generates the external surface and calculates each vertex’s local Gyrification Index (localGI) [30]. Values of Average Thickness, Total Area, Exposed Area, and Local Gyrification Index were extracted with Cortical Folding Analysis Tool [31]. We defined as ROI the whole hemisphere, frontal, temporal, occipital, and lateral lobes (based on FreeSurfer definition of lobes). The lobes’ area measurements were corrected by their integrated Gaussian Curvature, removing the partition size effect and directly comparing lobes and hemisphere cortical folding. Datasets’ demographics, acquisition, and processing information are summarized in Supplementary Material, Table 4. ### Statistical data analysis The data harmonization procedure is based on the time evolution of the basic morphological variables *T*, A*T* or A*E*, modeled as an exponential decay. Using Linear Mixed Models (LMM), we subtracted the sample shift, estimated as the Standard Deviation of Sample Random Effects (log10(*Variable*) ∼ *Age* × *ROI* + (1|*Sample* : *ROI*)) in Healthy Controls. From this procedure we directly calculate the variables K, S and I and determine for all morphological variables of interest (also *T*, A*T*, A*E*,GI = *A**t**/A**e*) their baseline and uncertainties. This harmonization procedure was validated as the most probable description of the data using the Bayesian comparison framework described in the appendix B. More than that, the Bayesian approach cross-validated all results presented here. Multiple comparisons of means were made with ANOVA, correlations between cortical folding independent variables, and age estimated with Pearson r and post hoc evaluations with Tukey multiple comparisons of means. The statistical significance threshold was set at *α* = 0.05, and multiple corrections (Bonferroni) were applied when needed. Healthy Control subjects are used as a reference in the necessary normalization and harmonization procedures. A standard linear regression was used to obtain the fractal index *α*. All analyses were done considering the whole cortical hemisphere and the lobes (see appendix C) as most diseases imply local or nonuniform global structural damage. LMM, standard linear regression and statistics derived from their results were analyzed with RStudio (v1.4.1717 and R v4.1.1). The Bayesian model comparison was developed in Python3 with the package PyMC3. ### Data and code availability Data from previous publications are available from their repositories: Wang et al., 2016 [11] and Wang et al., 2019 [12]. Data from AHEAD [21] and AOMIC [22] datasets with the cortical folding variables extracted by these authors are available [32]. Data acquired at IDOR does not have clearance for public sharing patients’ information that could lead to identification due to local Ethics Committee approval restrictions but will be shared upon request. Informed consent was obtained from all subjects for the published datasets as part of the original data acquisition. ## Results ### Multisite harmonization The analysis starts with data harmonization, a crucial step to combine information from multiple sites allowing a fair comparison among datasets. Figure 1 shows the result of the harmonization procedure and its reflection on the primary morphological variables. The LMMs are used to fit the data assuming a single angular coefficient for all samples and different linear coefficients due to type B errors. We considered the variables trend with age linear based on the correlation between *T*, A*T*, and A*E* with Age for the Healthy Control group calculated from our data. There was a significant correlation for T (Pearson r = −0.49; p < 0.0001), A*T* (Pearson r = −0.11; p < 0.0001), and A*E* (Pearson r = −0.05; p < 0.0001). As shown in appendix B, the harmonization was validated as the most probable description of this data. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F1) Figure 1: Basic morphological variables through age for Healthy Controls (A) raw data, (B) after harmonizing, removing the estimated residual from the Linear Mixed Model. This procedure results in the harmonized values of the novel morphological variables K, S and I. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F2) Figure 2: K, S and I through age for Healthy Controls after harmonizing in the primary morphological variables T, A*T*, and A*E* (removing the estimated residual from the Linear Mixed Model). ### Baseline, rate estimates, and Alzheimer’s Disease diagnostic discrimination An immediate consequence of the harmonization is defining a baseline value and its uncertainty for any morphological variable considered in this work. The LMM provide quantitative information about the trends of these variables through age. The estimated uncertainties values are summarized in Table 2 and summaries of linear mixed models are available in Supplementary Material A, Table 5. The so-called Natural Fluctuation represents human diversity and was considered the Residual Standard Deviation of the models. Moreover, we estimated the random component of the experimental error derived from the multiple-image acquisition and processing, described in Supplementary Material C.1. We considered the Random Effect of the Sample category (intercept) Standard Deviation as the type B component of the experimental error. The *σ**X* is the resultant uncertainty, composed of the previously described errors, and is the square root sum of each one’s squared values. Thereby, *σ**K* = 0.026 (4.9%), *σ**S* = 0.14 (1.6%), and *σ**I* = 0.091 (0.9%) with the corresponding fraction of their healthy controls intercept. View this table: [Table 2:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T2) Table 2: Summary of the estimated uncertainty components for each morphological variable of interest. The baseline determined by our method is a standard that can be used to compare Health Controls with any atypical group (in terms of brain morphology). It is directly calculated using LMM from the time evolution of the independent morphological variables as shown in 2. Considering Alzheimer’s Disease, we used the baseline to verify if the age trends of K, S, and I are different between the diseased and Healthy Control groups. We extracted the rates of change for each variable and diagnostic, summarized in Table 3 and displayed in Figure 3. The slopes were compared with post hoc pairwise mean comparisons. View this table: [Table 3:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T3) Table 3: Changing rate per year for each variable (after harmonization) and diagnostic. Mean value ± standard deviation. Rate [%] was estimated as the fraction from the mean value. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F3) Figure 3: Fitted values extracted from the linear regression model after data harmonization. Bars represent the 95% confidence interval. Complete summary of linear models in Supplementary Material A, Table 5. (A) Concerning K, Alzheimer’s Disease (AD) has a shallow slope, meaning small changes with Age. Compared to healthy controls, AD values of K are almost constant and similar to older subjects (AD slope p-value < 0.0001 and CTL slope p-value < 0.0001; pairwise comparison estimate −0.00090, p < 0.0001). (B) For S, the AD and CTL patterns have similar intercepts, but statistically different slopes (AD slope p-value = 0.004 and CTL slope p-value < 0.0001; pairwise comparison estimate 0.0013, p = 0.004). (C) For I, that reflects brain volume, CTL has a decreasing volume with aging, while AD has a smaller slope (AD slope p-value < 0.0001 and CTL slope p-value < 0.0001; pairwise comparison estimate −0.0017, p < 0.0001). The curves for the AD group were extended to early ages to ease the comparison. There is a significant difference between Control and AD slopes in K (pairwise comparison estimate −0.00090, p < 0.0001), S (pairwise comparison estimate 0.0013, p = 0.004) and I (pairwise comparison estimate −0.0017, p < 0.0001) as expected from Figure 3. For K, the AD curve is virtually flat as if it had reached a plateau, while for I, we see a constant reduction from brain isometric volume, with a reduced intercept. Comparing the mean values within AD and CTL, there is a significant difference for K (pairwise comparison estimate 0.042, p < 0.0001, representing 7.9% of the CTL mean) value, S (pairwise comparison estimate −0.14, p < 0.0001, 1.5%) and I (pairwise comparison estimate 0.11, p < 0.0001, 1.0%). In order to complete the full description of the cortical morphology within the framework established by the Mota & Herculano-Houzel’s model, we separated the data from both Health controls and AD patients per decade and fitted the fractal dimension *α* theorized to be 1.25. The result is shown in Figure 4 ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F4) Figure 4: Slope *α* derived from the Cortical Folding model from Mota & Herculano-Houzel through age for Healthy Controls and AD subjects. Data points from (0, 5] and [95, 100] years old were omitted due to the reduced data with 95% Confidence Interval. The numbers on top of each point indicates the number of subjects at each interval. As most diseases imply local or nonuniform global structural damage, we expanded the results of trajectories to the lobes. We included ROI as a fixed effect in the linear mixed model equation with harmonized variables: *Variable* ∼ *Age* × *Diagnostic* × *ROI* + (1 *Sample* : *Diagnostic* : *ROI*). Rates are summarized in Supplementary Material C. We discriminated AD and healthy aging (CTL) trends with K at the Frontal (p adj < 0.0001), Occipital (p adj < 0.0001), Parietal (p adj < 0.0001), and Temporal lobes (p adj < 0.0001), as expected from [20]. Despite the rates being smaller for AD than CTL, the results mimic the whole brain with Alzheimer’s Disease trajectories appearing as a plateau during aging. The Frontal and Parietal lobe shows the least folded pattern (K) in AD (within lobe comparison are described in Supplementary Material C). Also, there are significant differences with S CTL and AD at the Frontal (p adj = 0.00035) and Parietal lobes (p adj < 0.0001); and among I, between CTL and AD at the Frontal (p adj < 0.0001), Occipital (p adj = 0.0031), Parietal (p adj < 0.0001), and Temporal lobes (p adj < 0.0001)). ## Discussion We analyzed a combination of datasets aiming to establish K, S, and I baseline through human lifespan. The suggested model included complexity at two levels: accounting for heterogeneous acquisition and processing and heterogeneous brain structures across healthy and disease subjects. The hypothesis is that non-homogeneous methodology (inclusion criteria, acquisition, and processing) implies more significant uncertainty due to type B errors. Thereby, we suggested a harmonization procedure that allows the comparison of data from multisite reducing this kind of experimental errors. This result is innovative by its simplicity and focuses on expanding our independent morphological variables knowledge. The main caveat of our procedure is the impossibility of accounting for all differences between the samples and using this to define a universal harmonization according to each acquisition protocol. Our method does not depend on one reference sample or a gold standard definition. The removal of the shift for each sample is carried out based on the available healthy control data, which can always be added to the analysis. In other words, our method does not account for a possible global type B error that could shift all samples equally. By looking at the uncertainty in cortical thickness estimation, our results are compatible with previous work. Our findings of 0.02 mm (Table 2) of uncertainty in cortical thickness are in agreement with a variability less than 0.03 mm found in a test-retest analysis by Han et al. [33]. Moreover, they suggest that uncertainties of 0.15 mm across platforms of manufacturers and 0.17 mm across field strength, at this work, coupled with processing uncertainties and estimated in 0.14 mm. A most recent study from Frangou et al. [34] suggests a Mean Inter-individual variation of 0.07 ± 0.06 mm for hemisphere’s cortical thickness at all age ranges in the study compared to 0.10 mm from our findings. This compatibility suggests that the possible global type B error is close to 0; thus the harmonization could be regarded as being close to the universal one. The typical values for the K, S, and I independent variables through the whole life span and their related uncertainties were obtained by this work. These curves not only can be used to discriminate pathological from healthy aging in cross-sectional harmonized data but also are a powerful tool to morphologically describe non-typical cortex. Estimating the uncertainties is crucial to constructing the baseline curves and was carefully developed by this work. We estimated three types of uncertainties with linear mixed models and repeated measures: type B errors (combination of multiple datasets), repeated measures, and the natural variation in the human species. The error from acquisition and processing were coupled at the individual and sample levels. We expect future works to disassociate the uncertainties in detailed factors (test-retest analysis, across platforms, field strengths, and acquisition parameters with the same sample of subjects) to better measure morphological variables of interest. We have demonstrated the usage of the baseline curves by analyzing AD data. In terms of K, a measure of axonal tension, the AD curve is virtually flat as if it had reached a plateau of reduced brain tension. Comparing with Healthy Controls, this plateau seems to be the same experienced by the oldest subjects. Interpreting in terms of the Mota & Herculano-Houzel cortex model, AD can be regarded as a form of accelerated aging from the axonal tension perspective. On the other hand, the evolution of S, I, and even the fractal dimension *α* suggest that the AD cortex is not morphologically similar to the health one, with a very different shape, isometric volume, and time evolution. However, it is interesting to note that both S and I converge to a value compatible with the oldest health subjects. In the future, correlating these results with the elasto-mechanic properties of the AD cortex will help to understand better the physical limits of Mota & Herculano-Houzel model and expand it theoretically. A significant result that could only be achieved by combining multiple data is the value of the fractal dimension *α* through human lifespan and its evolution through atypical conditions. We extended the analysis in de Moraes et al. [35] and verified how the slope *α* behaves with healthy aging. We found that the slope is compatible with the theoretical prediction for subjects between 20 and 60 years old. Defining this range of applicability of the theory is essential to understand the period of life in which the model’s basic assumptions are still valid. After, we confirmed the findings that the slope escapes the model after 60 years old, diverging from the theoretical value of 1.25. Before 15 years old, we have found a hint that the slope is also inferior to 1.25, in disagreement with the model. Along, a non-linear pattern is suggested by the K baseline function before 20 years. One of the leading hypotheses to explain the deviation of the measured *α* from the expected value in both cases (before 15 and after 60 years old) is the breakdown of the homogeneity on the cortical surface. The suggestion of a homogeneity breakdown makes it particularly interesting to further study AD in the age range compatible with 1.25. Understanding the regimes outside the model is a fundamental step towards the theoretical understanding of cortical folding. Coupled with analysis of neuropathologies such as AD, these empirical results improve the cortical folding theory. More than discriminating disease and healthy subjects, understanding the pathology trajectory is vital to defining landmarks, indicating future injury, and collaborating to develop efficient treatment. Notably, in Alzheimer’s Disease and its prodromal form, Mild Cognitive Impairment, subsequent phases of damage may start at least 20 years before the diagnostic [36]. The structural damage is an initial state than cognitive alterations, as memory and clinical functions [37]. As pointed out by Fjell et al. [38], healthy aging also contributes to reductions in cortical thickness, and understanding the overlap between healthy aging and pathological aging triggers is crucial to segregate both events. Nevertheless, the possible similarity of structural damage in Alzheimer’s Disease and aging is explored by multiple previous studies [13, 39–41]. We hypothesize that the structural changes due to AD are reflected in measures of the brain’s global structure and shape. Therefore, cortical folding derived variables K, S, and I are sensitive to this perturbation. We verified this hypothesis by analyzing the difference in aging trajectories between healthy and pathological brains. At a more refined grain scale, the lobes, the results are compelling to confirm [13] suggestion that K, S, and I together would be a powerful tool to discriminate similarly but distinct events, as AD and aging. In this case, the difference between events is in the cortical structure’s changing velocity or rate. Then with the rates for K, S, and I at the lobes, we could explore whether a lobe would unfold faster. For AD, the temporal lobe has the biggest unfolding (K) and shrinking rate (I), while healthy aging is more aggressive at the Parietal lobe, followed by the temporal and frontal lobes. The natural expansion of this work is to include Mild Cognitive Impairment (MCI) subjects. As seen before, the MCI is an intermediary step between healthy aging and AD in the disease progression, cognitive impairment, and structural changes [20, 36, 37]. Also, we intend to continue this work by applying this approach to a longitudinal dataset composed of MRI structural images with MCI subjects that eventually convert to AD and compare to the cortical folding components of non-converters. Including this data would allow having a more reasonable aging trajectory by having the transition period from healthy aging to pathological aging. Further, we expect to explore K, S, and I function with age for the age range comprising newborns to 25 years old in future studies. This study will help better understand the value of the slope *α* at this age, revealing if there is an inflection point that can be connected to the development/aging transition in human brains [42], such as suggested by the data. This work successfully achieved its goal of combining multiple samples with a simplistic harmonization procedure that adjusts K, S, and I values for the 3650 subjects. Hereafter we estimated K, S, and I typical values and aging rates to discriminate against Alzheimer’s Disease subjects based on their changing rates. To validate the proposed methodology, we analyzed results for Cortical Thickness, which are in concordance with previously published findings. At a glance, the significant differences found in aging rates for the cross-sectional data are suggestive that a brain with Alzheimer’s Disease has premature and accelerated aging in terms of morphology. ## Data Availability Data from previous publications are available from their repositories: Wang et al., 2016 and Wang et al., 2019. Data from AHEAD and AOMIC datasets with the cortical folding variables extracted by these authors are available. Data acquired at IDOR does not have clearance for public sharing patients' information that could lead to identification due to local Ethics Committee approval restrictions but will be shared upon request. [https://zenodo.org/record/5750619](https://zenodo.org/record/5750619) [https://www.pnas.org/doi/pdf/10.1073/pnas.1610175113](https://www.pnas.org/doi/pdf/10.1073/pnas.1610175113) [https://zenodo.org/record/2626779](https://zenodo.org/record/2626779) ## Author Contributions F.H.P.M., V.B.B.M, and B.M. designed research; F.H.P.M., V.B.B.M, F.T.M., and B.M. performed research; F.H.P.M., V.B.B.M and B.M. analyzed data; and F.H.P.M., V.B.B.M and B.M. wrote the paper. All authors reviewed the manuscript. ## Supplementary Material View this table: [Table 4:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T4) Table 4: Summary for each dataset. Mean value ± standard deviation. Cross-sectional analyses were performed compilating all samples into one heterogeneous dataset. View this table: [Table 5:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T5) Table 5: Linear mixed models summary. Raw data, with no harmonization. 4364 observations; 11 Diagnostic:Sample; 8 Sample. Equation: *Variable* ∼ *Age × Diagnostic* + (1|*Sample*: *Diagnostic*). ### A Supplementary Tables, Figures and Results ### B Bayesian model selection In order to show that the harmonization procedure used were the most suitable, one can compare it with other models to describe the data. For example, the angular coefficient could also be a random effect for the LMM, varying across the samples. Using Bayesian model selection [43], one can show that the LMM used in this work is the most probable description to the data. Considering our dataset comprised by *J* samples where ![Graphic][6] denotes the data of a morphological variable in the *j*-th sample. The degree of belief in a model *M* that describes a normally distributed data for each age is given by the posterior probability: ![Formula][7] ![Formula][8] The term *µ*(*x**i*(*j*), *A*(*j*)*β*) is given by the model *M* proposed to describe the data, with the set {*A*(*j*)*β*} being the free parameters of the model defined within a range Δ*A*(*j*)*β*. A prior probability *P* (*A*(*j*)*β*|*I*) is assigned for all free parameters reflecting the prior knowledge about their values. The basic model of this work assumes that the log of morphological variables follow *µ* = *ax* + *b*(*j*), with the dispersion *σ*(*j*) for each samples. Uninformative priors were used for the free parameters. This integral can be easily done numerically using Python3 with PyMC3 package or similar. The advantage of this approach is that it naturally penalizes the insertion of unnecessary new parameters avoiding over-fitting. With multiple models *M**i*, one defines the odds as ![Graphic][9], the probability of the model *M*1 is given by: ![Formula][10] The alternative models considered by this work were: i) *µ* = *ax*+*b*(*j*), with the same dispersion *σ*(*j*); ii) *µ* = *a*(*j*)*x*+*b*(*j*), with the dispersion *σ*(*j*) for each sample; iii) *µ* = *a*(*j*)*x* + *b*(*j*), with the same dispersion *σ*(*j*) for each sample; iv) We also tested the simplest non-linear trend *µ* = *cx*2 + *ax* + *b*(*j*), with the dispersion *σ*(*j*) for each sample; We have found that the LMM used in the work is the most probable. The code implementing this approach is given by the authors. It is important to note the all parameters calculated using the LMM approach can be obtained in this framework by a slight modification of equation 7. Instead of marginalizing over all the parameters, by left one out, we obtain *P* (*A**j*1|*D*) which gives directly the most probable value for the parameter *A**j*1 and its credible region. ### C Expanded results to lobes We expanded the analyses of trajectories to the lobes by including ROI as a fixed effect in the linear mixed model equation: *Variable* ∼ *Age* × *Diagnostic* × *ROI* + (1|*Sample* : *Diagnostic* : *ROI*). Rates are summarized in Table 6 and Figure 5. Differences within lobes and rate differences within diagnostics are described in Figure 6. View this table: [Table 6:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T6) Table 6: Changing rate per year for each variable, diagnostic and lobe. Mean value *±* standard deviation. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F5) Figure 5: Fitted values for (A) K, (B) S and (I) and age for all for the Frontal lobe (“F”), Parietal lobe (“P”), Occipital lobe (“O”) and Temporal lobe (“T”) after data harmonization. Bars represents 95% confidence interval. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F6) Figure 6: Difference in slopes for each lobe within diagnostics after data harmonization. Bars represent a 95% confidence interval. (A) K, (B) S, and (C) I. #### C.1 Repeated measurements The random error is an essential fraction of the global uncertainties, as it is the variance that occurred after multiples measures in the same condition. To estimate the reproducibility within subjects, we processed (as described in Materials and Methods) the first 50 subjects with all three T1w images from the AOMIC ID1000 [22] dataset (described in Supplementary Material Table 7). The AOMIC ID1000 comprises multiple Magnetic Resonance Imaging protocols, especially three structural T1w images acquired during the same scan with the same protocol in a Philips Intera 3T. The select sample has the advantage of being composed of only healthy human subjects within a narrow age range, from 20 to 26 years old. We hypothesized that should be no difference of means within the three runs, despite the variation for each subject (Figure 7). View this table: [Table 7:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T7) Table 7: Summary for the 50 subjects of AOMIC ID1000 dataset. Mean value *±* standard deviation. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F7) Figure 7: Cortical folding model with the traced path for each hemisphere of each subject (gray line). The sample respects the model with slope *α* = 1.23 *±* 0.06, 95% confidence interval = (1.12,1.34). To estimate the uncertainty/variation of repeated measures, one could estimate the distribution of standard deviations of the data. We calculated the standard deviation for each subject and hemisphere across the three images, for the Cortical Thickness, K, S, and I. We compared the means for each run with a Repeat Measure ANOVA. There is no significant difference of means through runs (1, 2 and 3) (Figure 8): Cortical Thickness [mm], F(2, 198) = 0.481, p = 0.62, eta2[g] = 0.000192; GI, F(2, 186) = 0.196, p = 0.81, eta2[g] = 0.00001; K, F(2, 198) = 0.138, p = 0.87, eta2[g] = 0.0000144; S, F(2, 185) = 0.352, p = 0.69, eta2[g] = 0.0000798; I, F(2, 198) = 0.775, p = 0.46, eta2[g] = 0.0000523. The estimate standard deviations are: Cortical Thickness [mm], 0.019 *±* 0.013 mm; GI, 0.0069 *±* 0.0047; K, 0.0017 *±* 0.0012; S, 0.014 *±* 0.01; I, 0.0064 *±* 0.0043 (Figure 9). ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F8) Figure 8: Trajectories for each hemisphere of each subject through the acquisition runs. There is no significant difference in means. (A) Cortical Thickness [mm]: F(2, 198) = 0.481, p = 0.62, eta2[g] = 0.000192; (B) K: F(2, 198) = 0.138, p = 0.87, eta2[g] = 0.0000144; (C) S: F(2, 185) = 0.352, p = 0.69, eta2[g] = 0.0000798; (D) I: F(2, 198) = 0.775, p = 0.46, eta2[g] = 0.0000523. ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/15/2022.03.10.22272228/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/F9) Figure 9: Standard Deviation (SD) distribution for each variable. The solid line represents the mean value, and the dashed line the standard deviation of the distribution. (A) Cortical Thickness: 0.019*±*0.013 mm; (B) K: 0.0017*±*0.0012; (C) S: 0.014*±*0.01; (D) I: 0.0064*±*0.0043. These results suggest high reliability in subsequently acquired images processed with FreeSurfer v6.0, in concordance with [44] finds for reliability in FreeSurfer processing comparing volumes measurements. We suggest that the estimated standard errors accurately estimate the uncertainty due to a repeated acquisition and image processing within the same site and parameters. ### D multicenter influence on cortical folding model and its independent morphological components We estimated the influence of multicenter datasets in estimating K, S, and I typical values and rates per year across the human lifespan, including as random effects at the linear mixed model. Using these models allows one to calculate the residual variance arising from the methodological heterogeneity and the variance from each category, the different included data samples in this manuscript. The systematic shifts, estimated from the Linear Mixed Models (LMM, *Variable* ∼ *Age* × *Diagnostic* × *ROI* + (1|*Sample* : *Diagnostic* : *ROI*)), can be used as inputs to multisite harmonization procedures, reducing the influence of heterogeneous methodology (Figure 1). The harmonization consists of subtracting the estimated shift from the raw data of each sample (*X**harmonized* = *X**raw* − *shift**X*). It reduces the systematic error due to the samples, and consecutively, the total uncertainty of the variable (Table 8). A step-by-step code is available [35]. View this table: [Table 8:](http://medrxiv.org/content/early/2022/03/15/2022.03.10.22272228/T8) Table 8: Comparison of estimated error after the multisite harmonization. Future studies can progress on the limitation by detailing influences and increasing specificity in methodology, such as the scan field strength, model and manufacturer, and decoupling acquisition imaging processing components. ## Acknowledgements This work was funded by the Research Support Foundation of the State of Rio de Janeiro (FAPERJ), National Council for Scientific and Technological Development (CNPq), and intramural grants from Instituto D’Or de Pesquisa e Ensino (IDOR). Bruno Mota is supported by Instituto Serrapilheira (grant Serra-1709-16981) and CNPq (PQ 2017 312837/2017-8). We are thankful to the research team and volunteers for participating in this research project. The authors have no competing interests. ## Footnotes * Correct 3rd author's name in the system, not the pdf. * 1 One may understand this as similar to “Principal Component Analysis”, a change of basis on the data representation, but using scaling considerations rather than covariance eigenvectors to define the principal directions. The new set of variables are combinations of the usual geometrical variables used to describe the cortex. The principal component K is given by the power-law relation, which contains the correlation between the variables. The variables S and I are orthonormal to K. * 2 Structural T1w images from AOMIC ID1000 dataset, described in Material and Methods. * Received March 10, 2022. * Revision received March 14, 2022. * Accepted March 15, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1]. David C. Van Essen. A 2020 view of tension-based cortical morphogenesis. Proceedings of the National Academy of Sciences, 117(52):32868–32879, December 2020. Publisher: National Academy of Sciences Section: Biological Sciences. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzUyLzMyODY4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMTUvMjAyMi4wMy4xMC4yMjI3MjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 2. [2]. Cristina Llinares-Benadero and Víctor Borrell. Deconstructing cortical folding: genetic, cellular and mechanical determinants. Nature Reviews. Neuroscience, 20(3):161–176, March 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41583-018-0112-2&link_type=DOI) 3. [3]. Christopher R. Madan and Elizabeth A. Kensinger. Cortical complexity as a measure of age-related brain atrophy. NeuroImage, 134:617–629, July 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2016.04.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27103141&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) 4. [4]. Bruno Mota, Sandra E. Dos Santos, Lissa Ventura-Antunes, Débora Jardim-Messeder, Kleber Neves, Rodrigo S. Kazu, Stephen Noctor, Kelly Lambert, Mads F. Bertelsen, Paul R. Manger, Chet C. Sherwood, Jon H. Kaas, and Suzana Herculano-Houzel. White matter volume and white/gray matter ratio in mammalian species as a consequence of the universal scaling of cortical folding. Proceedings of the National Academy of Sciences, 116(30):15253–15261, July 2019. Publisher: National Academy of Sciences Section: PNAS Plus. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzMwLzE1MjUzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMTUvMjAyMi4wMy4xMC4yMjI3MjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. [5]. Bo Cao, Benson Mwangi, Ives Cavalcante Passos, Mon-Ju Wu, Zafer Keser, Giovana B. Zunta-Soares, Dianping Xu, Khader M. Hasan, and Jair C. Soares. Lifespan Gyrification Trajectories of Human Brain in Healthy Individuals and Patients with Major Psychiatric Disorders. Scientific Reports, 7(1):1–8, March 2017. Number: 1 Publisher: Nature Publishing Group. 6. [6]. Sander Lamballais, Elisabeth J. Vinke, Meike W. Vernooij, M. Arfan Ikram, and Ryan L. Muetzel. Cortical gyrification in relation to age and cognition in older adults. NeuroImage, 212:116637, May 2020. 7. [7]. Christian Núñez, Antonio Callén, Federica Lombardini, Yaroslau Compta, Christian Stephan-Otto, and Alzheimer’s Disease Neuroimaging Initiative. Different Cortical Gyrification Patterns in Alzheimer’s Disease and Impact on Memory Performance. Annals of Neurology, 88(1):67–80, July 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25741&link_type=DOI) 8. [8]. Larson J. Hogstrom, Lars T. Westlye, Kristine B. Walhovd, and Anders M. Fjell. The structure of the cerebral cortex across adult life: age-related patterns of surface area, thickness, and gyrification. Cerebral Cortex (New York, N.Y.: 1991), 23(11):2521–2530, November 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhs231&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22892423&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325761100001&link_type=ISI) 9. [9]. Christopher R. Madan. Age-related decrements in cortical gyrification: Evidence from an accelerated longitudinal dataset. The European Journal of Neuroscience, 53(5):1661–1671, March 2021. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ejn.15039&link_type=DOI) 10. [10]. Bruno Mota and Suzana Herculano-Houzel. Cortical folding scales universally with surface area and thickness, not number of neurons. Science (New York, N.Y.), 349(6243):74–77, July 2015. 11. [11]. Yujiang Wang, Joe Necus, Marcus Kaiser, and Bruno Mota. Universality in human cortical folding in health and disease. Proceedings of the National Academy of Sciences of the United States of America, 113(45):12820–12825, November 2016. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEzLzQ1LzEyODIwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMTUvMjAyMi4wMy4xMC4yMjI3MjIyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. [12]. Yujiang Wang, Joe Necus, Luis Peraza Rodriguez, Peter Neal Taylor, and Bruno Mota. Human cortical folding across regions within individual brains follows universal scaling law. Communications Biology, 2(1):191, May 2019. 13. [13]. Yujiang Wang, Karoline Leiberg, Tobias Ludwig, Bethany Little, Joe H Necus, Gavin Winston, Sjoerd B Vos, Jane de Tisi, John S Duncan, Peter N Taylor, and Bruno Mota. Independent components of human brain morphology. NeuroImage, 226:117546, February 2021. 14. [14]. Jean-Philippe Fortin, Nicholas Cullen, Yvette I. Sheline, Warren D. Taylor, Irem Aselcioglu, Philip A. Cook, Phil Adams, Crystal Cooper, Maurizio Fava, Patrick J. McGrath, Melvin McInnis, Mary L. Phillips, Madhukar H. Trivedi, Myrna M. Weissman, and Russell T. Shinohara. Harmonization of cortical thickness measurements across scanners and sites. NeuroImage, 167:104–120, February 2018. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) 15. [15]. Ed H. B. M. Gronenschild, Petra Habets, Heidi I. L. Jacobs, Ron Mengelers, Nico Rozendaal, Jim van Os, and Machteld Marcelis. The Effects of FreeSurfer Version, Workstation Type, and Macintosh Operating System Version on Anatomical Volume and Cortical Thickness Measurements. PLOS ONE, 7(6):e38234, June 2012. Publisher: Public Library of Science. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0038234&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22675527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) 16. [16]. B. C. Dickerson, E. Fenstermacher, D. H. Salat, D. A. Wolk, R. P. Maguire, R. Desikan, J. Pacheco, B. T. Quinn, A. Van der Kouwe, D. N. Greve, D. Blacker, M. S. Albert, R. J. Killiany, and B. Fischl. Detection of cortical thickness correlates of cognitive performance: Reliability across MRI scan sessions, scanners, and field strengths. NeuroImage, 39(1):10–18, January 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2007.08.042&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17942325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251406000002&link_type=ISI) 17. [17]. Hugo G. Schnack, Neeltje E. M. van Haren, Rachel M. Brouwer, G. Caroline M. van Baal, Marco Picchioni, Matthias Weisbrod, Heinrich Sauer, Tyrone D. Cannon, Matti Huttunen, Claude Lepage, D. Louis Collins, Alan Evans, Robin M. Murray, René S. Kahn, and Hilleke E. Hulshoff Pol. Mapping reliability in multicenter MRI: Voxel-based morphometry and cortical thickness. Human Brain Mapping, 31(12):1967–1982, 2010. _eprint: [https://onlinelibrary.wiley.com/doi/pdf/10.1002/hbm.20991](https://onlinelibrary.wiley.com/doi/pdf/10.1002/hbm.20991). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.20991&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21086550&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285259500012&link_type=ISI) 18. [18]. Raymond Pomponio, Guray Erus, Mohamad Habes, Jimit Doshi, Dhivya Srinivasan, Elizabeth Mamourian, Vishnu Bashyam, Ilya M. Nasrallah, Theodore D. Satterthwaite, Yong Fan, Lenore J. Launer, Colin L. Masters, Paul Maruff, Chuanjun Zhuo, Henry Völzke, Sterling C. Johnson, Jurgen Fripp, Nikolaos Koutsouleris, Daniel H. Wolf, Raquel Gur, Ruben Gur, John Morris, Marilyn S. Albert, Hans J. Grabe, Susan M. Resnick, R. Nick Bryan, David A. Wolk, Russell T. Shinohara, Haochang Shou, and Christos Davatzikos. Harmonization of large MRI datasets for the analysis of brain imaging patterns throughout the lifespan. NeuroImage, 208:116450, March 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016j.neuroimage.2019.116450&link_type=DOI) 19. [19]. R.A.I. Bethlehem, J. Seidlitz, S.R. White, J.W. Vogel, K.M. Anderson, C. Adamson, S. Adler, G.S. Alexopoulos, E. Anagnostou, A. Areces-Gonzalez, D.E. Astle, B. Auyeung, M. Ayub, G. Ball, S. Baron-Cohen, R. Beare, S.A. Bedford, V. Benegal, F. Beyer, J. Bin Bae, J. Blangero, M. Blesa Cábez, J.P. Boardman, M. Borzage, J.F. Bosch-Bayard, N. Bourke, V.D. Calhoun, M.M. Chakravarty, C. Chen, C. Chertavian, G. Chetelat, Y.S. Chong, J.H. Cole, A. Corvin, E. Courchesne, F. Crivello, V.L. Cropley, J. Crosbie, N. Crossley, M. Delarue, S. Desrivieres, G. Devenyi, M.A. Di Biase, R. Dolan, K.A. Donald, G. Donohoe, K. Dunlop, A.D. Edwards, J.T. Elison, C.T. Ellis, J.A. Elman, L. Eyler, D.A. Fair, P.C. Fletcher, P. Fonagy, C.E. Franz, L. Galan-Garcia, A. Gholipour, J. Giedd, J.H. Gilmore, D.C. Glahn, I. Goodyer, P.E. Grant, N.A. Groenewold, F.M. Gunning, R.E. Gur, R.C. Gur, C.F. Hammill, O. Hansson, T. Hedden, A. Heinz, R. Henson, K. Heuer, J. Hoare, B. Holla, A.J. Holmes, R. Holt, H. Huang, K. Im, J. Ipser, C.R. Jack, A.P. Jackowski, T. Jia, K.A. Johnson, P.B. Jones, D.T. Jones, R. Kahn, H. Karlsson, L. Karlsson, R. Kawashima, E.A. Kelley, S. Kern, K. Kim, M.G. Kitzbichler, W.S. Kremen, F. Lalonde, B. Landeau, S. Lee, J. Lerch, J.D. Lewis, J. Li, W. Liao, D.P. Linares, C. Liston, M.V. Lombardo, J. Lv, C. Lynch, T.T. Mallard, M. Marcelis, R.D. Markello, B. Mazoyer, P. McGuire, M.J. Meaney, A. Mechelli, N. Medic, B. Misic, S.E. Morgan, D. Mothersill, J. Nigg, M.Q.W. Ong, C. Ortinau, R. Ossenkoppele, M. Ouyang, L. Palaniyappan, L. Paly, P.M. Pan, C. Pantelis, M.M. Park, T. Paus, Z. Pausova, A. Pichet Binette, K. Pierce, X. Qian, J. Qiu, A. Qiu, A. Raznahan, T. Rittman, C.K. Rollins, R. Romero-Garcia, L. Ronan, M.D. Rosenberg, D.H. Rowitch, G.A. Salum, T.D. Satterthwaite, H.L. Schaare, R.J. Schachar, A.P. Schultz, G. Schumann, M. Schöll, D. Sharp, R.T. Shinohara, I. Skoog, C.D. Smyser, R.A. Sperling, D.J. Stein, A. Stolicyn, J. Suckling, G. Sullivan, Y. Taki, B. Thyreau, R. Toro, K.A. Tsvetanov, N.B. Turk-Browne, J.J. Tuulari, C. Tzourio, É. Vachon-Presseau, M.J. Valdes-Sosa, P.A. Valdes-Sosa, S.L. Valk, T. van Amelsvoort, S.N. Vandekar, L. Vasung, L.W. Victoria, S. Villeneuve, A. Villringer, P. E. Vértes, K. Wagstyl, Y.S. Wang, S.K. Warfield, V. Warrier, E. Westman, M.L. Westwater, H.C. Whalley, A.V. Witte, N. Yang, B.T.T. Yeo, H.J. Yun, A. Zalesky, H.J. Zar, A. Zettergren, J.H. Zhou, H. Ziauddeen, A. Zugman, X.N. Zuo, AIBL, Alzheimer’s Disease Neuroimaging Initiative, Alzheimer’s Disease Repository Without Borders Investigators, ASRB, CALM Team, Cam-CAN, 3R-BRAIN CCNP, COBRE, ENIGMA Developmental Brain Age working group, FinnBrain, Harvard Aging Brain Study, KNE96 IMAGEN, NSPN, OASIS-3, OpenPain Project, POND, VETSA, Bullmore, E.T. The PREVENT-AD Research Group, and A.F. Alexander-Bloch. Brain charts for the human lifespan. bioRxiv, page 2021.06.08.447489, January 2021. 20. [20]. Fernanda Hansen Pacheco de Moraes, Felipe Sudo, Marina Carneiro Monteiro, Bruno R. P. de Melo, Paulo Mattos, Bruno Mota, and Fernanda Tovar-Moll. Independent morphological variables correlate with Aging, Mild Cognitive Impairment, and Alzheimer’s Disease. medRxiv, page 2022.01.10.22268812, January 2022. 21. [21]. Anneke Alkemade, Martijn J Mulder, Josephine M Groot, Bethany R Isaacs, Nikita van Berendonk, Nicky Lute, Scott JS Isherwood, Pierre-Louis Bazin, and Birte U Forstmann. The Amsterdam Ultra-high field adult lifespan database (AHEAD): A freely available multimodal 7 Tesla submillimeter magnetic resonance imaging database. NeuroImage, 221:117200, November 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2020.117200&link_type=DOI) 22. [22]. Lukas Snoek, Maite M. van der Miesen, Tinka Beemsterboer, Andries van der Leij, Annemarie Eigenhuis, and H. Steven Scholte. The Amsterdam Open MRI Collection, a set of multimodal MRI datasets for individual difference analyses. Scientific Data, 8(1):85, March 2021. 23. [23]. Clifford R Jack, Matt A Bernstein, Bret J Borowski, Jeffrey L Gunter, Nick C Fox, Paul M Thompson, Norbert Schuff, Gunnar Krueger, Ronald J Killiany, Charles S DeCarli, Anders M Dale, and Michael W Weiner. Update on the MRI Core of the Alzheimer’s Disease Neuroimaging Initiative. Alzheimer’s & dementia: the journal of the Alzheimer’s Association, 6(3):212–220, May 2010. 24. [24]. Matthew F. Glasser, Stamatios N. Sotiropoulos, J. Anthony Wilson, Timothy S. Coalson, Bruce Fischl, Jesper L. Andersson, Junqian Xu, Saad Jbabdi, Matthew Webster, Jonathan R. Polimeni, David C. Van Essen, and Mark Jenkinson. The minimal preprocessing pipelines for the Human Connectome Project. NeuroImage, 80:105–124, October 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2013.04.127&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23668970&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322416000011&link_type=ISI) 25. [25]. Kate Brody Nooner, Stanley J. Colcombe, Russell H. Tobe, Maarten Mennes, Melissa M. Benedict, Alexis L. Moreno, Laura J. Panek, Shaquanna Brown, Stephen T. Zavitz, Qingyang Li, Sharad Sikka, David Gutman, Saroja Bangaru, Rochelle Tziona Schlachter, Stephanie M. Kamiel, Ayesha R. Anwar, Caitlin M. Hinz, Michelle S. Kaplan, Anna B. Rachlin, Samantha Adelsberg, Brian Cheung, Ranjit Khanuja, Chaogan Yan, Cameron C. Craddock, Vincent Calhoun, William Courtney, Margaret King, Dylan Wood, Christine L. Cox, A. M. Clare Kelly, Adriana Di Martino, Eva Petkova, Philip T. Reiss, Nancy Duan, Dawn Thomsen, Bharat Biswal, Barbara Coffey, Matthew J. Hoptman, Daniel C. Javitt, Nunzio Pomara, John J. Sidtis, Harold S. Koplewicz, Francisco Xavier Castellanos, Bennett L. Leventhal, and Michael P. Milham. The NKI-Rockland Sample: A Model for Accelerating the Pace of Discovery Science in Psychiatry. Frontiers in Neuroscience, 6, October 2012. 26. [26]. Daniel S. Marcus, Anthony F. Fotenos, John G. Csernansky, John C. Morris, and Randy L. Buckner. Open Access Series of Imaging Studies (OASIS): Longitudinal MRI Data in Nondemented and Demented Older Adults. Journal of cognitive neuroscience, 22(12):2677–2684, December 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1162/jocn.2009.21407&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19929323&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281129400001&link_type=ISI) 27. [27]. Bruce Fischl. FreeSurfer. NeuroImage, 62(2):774–781, August 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.01.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22248573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306390600031&link_type=ISI) 28. [28]. Martin Reuter and Bruce Fischl. Avoiding Asymmetry-Induced Bias in Longitudinal Image Processing. NeuroIm-age, 57(1):19–21, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2011.02.076&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21376812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291624100004&link_type=ISI) 29. [29]. Christopher S. McCarthy, Avinash Ramprashad, Carlie Thompson, Jo-Anna Botti, Ioana L. Coman, and Wendy R. Kates. A comparison of FreeSurfer-generated data with and without manual intervention. Frontiers in Neuro-science, 9, 2015. 30. [30]. M. Schaer, M. B. Cuadra, L. Tamarit, F. Lazeyras, S. Eliez, and J. Thiran. A Surface-Based Approach to Quantify Local Cortical Gyrification. IEEE Transactions on Medical Imaging, 27(2):161–170, February 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TMI.2007.903576&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18334438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252815200002&link_type=ISI) 31. [31]. Yujiang Wang, Tobias Ludwig, and Bruno Mota. Analysis pipeline to extract cortical morphology measures from Freesurfer for cortical folding analysis. [https://zenodo.org/record/3608675](https://zenodo.org/record/3608675), October 2019. 32. [32]. Fernanda Hansen Pacheco de Moraes. OpenData for the Independent morphological variables correlate to aging, Mild Cognitive Impairment, and Alzheimer’s Disease manuscript, December 2021. Type: dataset. 33. [33]. Xiao Han, Jorge Jovicich, David Salat, Andre van der Kouwe, Brian Quinn, Silvester Czanner, Evelina Busa, Jenni Pacheco, Marilyn Albert, Ronald Killiany, Paul Maguire, Diana Rosas, Nikos Makris, Anders Dale, Bradford Dickerson, and Bruce Fischl. Reliability of MRI-derived measurements of human cerebral cortical thickness: The effects of field strength, scanner upgrade and manufacturer. NeuroImage, 32(1):180–194, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2006.02.051&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16651008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239530700020&link_type=ISI) 34. [34]. Sophia Frangou, Amirhossein Modabbernia, Gaelle E. Doucet, Efstathios Papachristou, Steven CR Williams, Ingrid Agartz, Moji Aghajani, Theophilus N. Akudjedu, Anton Albajes-Eizagirre, Dag Alnæs, Kathryn I. Alpert, Micael Andersson, Nancy Andreasen, Ole A. Andreassen, Philip Asherson, Tobias Banaschewski, Nuria Bargallo, Sarah Baumeister, Ramona Baur-Streubel, Alessandro Bertolino, Aurora Bonvino, Dorret I. Boomsma, Stefan Borgwardt, Josiane Bourque, Daniel Brandeis, Alan Breier, Henry Brodaty, Rachel M. Brouwer, Jan K. Buitelaar, Geraldo F. Busatto, Randy L. Buckner, Vincent Calhoun, Erick J. Canales-Rodríguez, Dara M. Cannon, Xavier Caseras, Francisco X. Castellanos, Simon Cervenka, Tiffany M. Chaim-Avancini, Christopher RK Ching, Vincent P. Clark, Patricia Conrod, Annette Conzelmann, Benedicto Crespo-Facorro, Fabrice Crivello, Eveline AM Crone, Anders M. Dale, Cristopher Davey, Eco JC de Geus, Lieuwe de Haan, Greig I. de Zubicaray, Anouk den Braber, Erin W. Dickie, Annabella Di Giorgio, Nhat Trung Doan, Erlend S. Dørum, Stefan Ehrlich, Susanne Erk, Thomas Espeseth, Helena Fatouros-Bergman, Simon E. Fisher, Jean-Paul Fouche, Barbara Franke, Thomas Frodl, Paola Fuentes-Claramonte, David C. Glahn, Ian H. Gotlib, Hans- Jörgen Grabe, Oliver Grimm, Nynke A. Groenewold, Dominik Grotegerd, Oliver Gruber, Patricia Gruner, Rachel E. Gur, Ruben C. Gur, Ben J. Harrison, Catharine A. Hartman, Sean N. Hatton, Andreas Heinz, Dirk J. Heslenfeld, Derrek P. Hibar, Ian B. Hickie, Beng-Choon Ho, Pieter J. Hoekstra, Sarah Hohmann, Avram J. Holmes, Martine Hoogman, Norbert Hosten, Fleur M. Howells, Hilleke E. Hulshoff Pol, Chaim Huyser, Neda Jahanshad, Anthony James, Jiyang Jiang, Erik G. Jönsson, John A. Joska, Rene Kahn, Andrew Kalnin, Ryota Kanai, Sim Kang, Marieke Klein, Tatyana P. Klushnik, Laura Koenders, Sanne Koops, Bernd Krämer, Jonna Kuntsi, Jim Lagopoulos, Luisa Lázaro, Irina Lebedeva, Won Hee Lee, Klaus-Peter Lesch, Christine Lochner, Marise WJ Machielsen, Sophie Maingault, Nicholas G. Martin, Ignacio Martínez-Zalacaín, David Mataix-Cols, Bernard Mazoyer, Colm McDonald, Brenna C. McDonald, Andrew M. McIntosh, Katie L. McMahon, Genevieve McPhilemy José, M. Menchón, Sarah E. Medland, Andreas Meyer-Lindenberg, Jilly Naaijen, Pablo Najt, Tomohiro Nakao, Jan E. Nordvik, Lars Nyberg, Jaap Oosterlaan, Víctor Ortiz-García de la Foz, Yannis Paloyelis, Paul Pauli, Giulio Pergola, Edith Pomarol-Clotet, Maria J. Portella, Steven G. Potkin, Joaquim Radua, Andreas Reif, Joshua L. Roffman, Pedro GP Rosa, Matthew D. Sacchet, Perminder S. Sachdev, Raymond Salvador, Pascual Sánchez-Juan, Salvador Sarró, Theodore D. Satterthwaite, Andrew J. Saykin, Mauricio H. Serpa, Lianne Schmaal, Knut Schnell, Gunter Schumann, Jordan W. Smoller, Iris Sommer, Carles Soriano-Mas, Dan J. Stein, Lachlan T. Strike, Suzanne C. Swagerman, Christian K. Tamnes, Henk S. Temmingh, Sophia I. Thomopoulos, Alexander S. Tomyshev, Diana Tordesillas-Gutiérrez, Julian N. Trollor, Jessica A. Turner, Anne Uhlmann, Odille A. van den Heuvel, Dennis van den Meer, Nic JA van der Wee, Neeltje EM van Haren, Dennis van ‘t Ent, Theo GM van Erp, Ilya M. Veer, Dick J. Veltman, Henry Völzke, Henrik Walter, Esther Walton, Lei Wang, Yang Wang, Thomas H. Wassink, Bernd Weber, Wei Wen, John D. West, Lars T. Westlye, Heather Whalley, Lara M. Wierenga, Katharina Wittfeld, Daniel H. Wolf, Margaret J. Wright, Kun Yang, Yulyia Yoncheva, Marcus V. Zanetti, Georg C. Ziegler, Paul M. Thompson, and Danai Dima. Cortical Thickness Trajectories across the Lifespan: Data from 17,075 healthy individuals aged 3-90 years. bioRxiv, page 2020.05.05.077834, May 2020. Publisher: Cold Spring Harbor Laboratory Section: New Results. 35. [35]. Fernanda Hansen P. de Moraes, Victor B. B. Mello, and Bruno Mota. Deaging and harmonization -Cortical Folding in Humans v1.1, September 2021. 36. [36]. Rudimar L. Frozza, Mychael V. Lourenco, and Fernanda G. De Felice. Challenges for Alzheimer’s Disease Therapy: Insights from Novel Mechanisms Beyond Memory Defects. Frontiers in Neuroscience, 12, 2018. Publisher: Frontiers. 37. [37]. Clifford R. Jack, David S. Knopman, William J. Jagust, Leslie M. Shaw, Paul S. Aisen, Michael W. Weiner, Ronald C. Petersen, and John Q. Trojanowski. Hypothetical model of dynamic biomarkers of the Alzheimer’s pathological cascade. The Lancet Neurology, 9(1):119–128, January 2010. 38. [38]. Anders M. Fjell, Linda McEvoy, Dominic Holland, Anders M. Dale, and Kristine B. Walhovd. What is normal in normal aging? Effects of aging, amyloid and Alzheimer’s disease on the cerebral cortex and the hippocampus. Progress in Neurobiology, 117:20–40, June 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pneurobio.2014.02.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24548606&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) 39. [39]. Jennifer Pacheco, Joshua O. Goh, Michael A. Kraut, Luigi Ferrucci, and Susan M. Resnick. Greater cortical thinning in normal older adults predicts later cognitive impairment. Neurobiology of Aging, 36(2):903–908, February 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neurobiolaging.2014.08.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25311277&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) 40. [40]. Benjamin Gutierrez Becker, Tassilo Klein, and Christian Wachinger. Gaussian process uncertainty in age estimation as a measure of brain abnormality. NeuroImage, 175:246–258, July 2018. 41. [41]. W. Huizinga, D. H. J. Poot, M. W. Vernooij, G. V. Roshchupkin, E. E. Bron, M. A. Ikram, D. Rueckert, W. J. Niessen, and S. Klein. A spatio-temporal reference model of the aging brain. NeuroImage, 169:11–22, April 2018. 42. [42]. Nick J. Davis. Quantifying the trajectory of gyrification changes in the aging brain (Commen-tary on Madan, 2021). European Journal of Neuroscience, n/a(n/a), April 2021. _eprint: [https://onlinelibrary.wiley.com/doi/pdf/10.1111/ejn.15220](https://onlinelibrary.wiley.com/doi/pdf/10.1111/ejn.15220). 43. [43]. Phil Gregory. Bayesian Logical Data Analysis for the Physical Sciences: A Comparative Approach with Mathematica® Support. Cambridge University Press, 2005. 44. [44]. Lucas D. Eggert, Jens Sommer, Andreas Jansen, Tilo Kircher, and Carsten Konrad. Accuracy and Reliability of Automated Gray Matter Segmentation Pathways on Real and Simulated Structural Magnetic Resonance Images of the Human Brain. PLOS ONE, 7(9):e45081, September 2012. Publisher: Public Library of Science. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0045081&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23028771&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F15%2F2022.03.10.22272228.atom) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/inline-graphic-1.gif [7]: /embed/graphic-15.gif [8]: /embed/graphic-16.gif [9]: /embed/inline-graphic-2.gif [10]: /embed/graphic-17.gif