- 1Institute of Medical Sciences, University of Toronto, Toronto, ON, Canada
- 2Centre for Addiction and Mental Health, Toronto, ON, Canada
- 3Campbell Family Mental Health Research Institute, Centre for Addiction and Mental Health (CAMH), Toronto, ON, Canada
- 4Department of Psychiatry, University of Toronto, Toronto, ON, Canada
- 5Departments of Pharmacology and Toxicology, University of Toronto, Toronto, ON, Canada
- 6Mouse Imaging Centre (MICe), Hospital for Sick Children, Toronto, ON, Canada
- 7Wellcome Centre for Integrative Neuroimaging, Oxford Centre for Functional MRI of the Brain, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
- 8Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada
Introduction: As the population skews toward older age, elucidating mechanisms underlying human brain aging becomes imperative. Structural MRI has facilitated non-invasive investigation of lifespan brain morphology changes, yet this domain remains uncharacterized in rodents despite increasing use as models of disordered human brain aging.
Methods: Young (2m, n = 10), middle-age (10m, n = 10) and old (22m, n = 9) mice were utilized for maturational (young vs. middle-age) and aging-related (middle-age vs. old mice) comparisons. Regional brain volume was averaged across hemispheres and reduced to 32 brain regions. Pairwise group differences in regional volume were tested using general linear models, with total brain volume as a covariate. Sample-wide associations between regional brain volume and Y-maze performance were assessed using logistic regression, residualized for total brain volume. Both analyses corrected for multiple comparisons. Structural covariance networks were generated using the R package “igraph.” Group differences in network centrality (degree), integration (mean distance), and segregation (transitivity, modularity) were tested across network densities (5–40%), using 5,000 (1,000 for degree) permutations with significance criteria of p < 0.05 at ≥5 consecutive density thresholds.
Results: Widespread significant maturational changes in volume occurred in 18 brain regions, including considerable loss in isocortex regions and increases in brainstem regions and white matter tracts. The aging-related comparison yielded 6 significant changes in brain volume, including further loss in isocortex regions and increases in white matter tracts. No significant volume changes were observed across either comparison for subcortical regions. Additionally, smaller volume of the anterior cingulate area (χ2 = 2.325, pBH = 0.044) and larger volume of the hippocampal formation (χ2 = −2.180, pBH = 0.044) were associated with poorer cognitive performance. Maturational network comparisons yielded significant degree changes in 9 regions, but no aging-related changes, aligning with network stabilization trends in humans. Maturational decline in modularity occurred (24–29% density), mirroring human trends of decreased segregation in young adulthood, while mean distance and transitivity remained stable.
Conclusion/Implications: These findings offer a foundational account of age effects on brain volume, structural brain networks, and working memory in mice, informing future work in facilitating translation between rodent models and human brain aging.
1. Introduction
As the world population skews sharply toward older age, with 21.1% projected to be 60 years or older by 2050, the risk of age-related physical and cognitive decline is expected to increase an already widespread societal burden (Chatterji et al., 2015). Understanding lifespan changes in brain maturation and aging is crucial for uncovering etiology driving the aging process and the emergence of age-related diseases, such as dementia (Hou et al., 2019).
Human structural magnetic resonance imaging (MRI) studies have begun to map trajectories of lifelong structural brain change. In humans, overall brain volume is observed to increase before a mid-childhood peak at 5.9 years, followed by a near-linear decrease across the lifespan (Bethlehem et al., 2022). On a regional level, maturational growth peaks earliest in primary sensory motor areas and latest in higher-order association areas (Giedd and Rapoport, 2010). Specifically, regional gray matter loss occurs most consistently in the cerebral cortex, with drastic loss in prefrontal regions accompanied by less extensive loss in parietal and temporal association cortices (Raz, 1997; Tisserand et al., 2002; Resnick et al., 2003). Moderately consistent findings of aging-related volume loss are reported for the insula, cerebellum, basal ganglia, and thalamus (Good et al., 2001; Alexander et al., 2012). However, findings of aging-related atrophy in limbic and prelimbic regions are more conflicting, particularly in the hippocampus, amygdala, and cingulate gyrus (Raz, 1997; Jernigan et al., 2001; Zheng et al., 2019). Compared to gray matter, white matter volume matures gradually, peaking at 50 years of age, before drastically declining after 60 years of age (Jernigan et al., 2001; Liu et al., 2016). In turn, age-related individual differences in regional gray and white matter properties have been mapped onto differences in cognition and may contribute to observed patterns of age-dependent cognitive decline (Raz and Rodrigue, 2006; Oschwald et al., 2019).
While structural MRI has facilitated progress in characterizing lifespan brain volume changes in humans, further leveraging these insights to study molecular mechanisms of aging is challenging due to the invasive nature of molecular brain research. The use of rodent models allows for more direct experimentation than human brain research, including transgenic and gene knock-out models to simulate pathology, drug trials, and more invasive biological interventions (Nadon, 2006). Moreover, the significantly shorter lifecycles in rodent models facilitate feasible lifespan and generational research (Folgueras et al., 2018). Ultimately, models are as useful as their ability to approximate healthy and pathological states in the human brain. However, while preclinical models are being increasingly implemented for Alzheimer’s disease, Parkinson’s disease, frontotemporal dementia, and Huntington’s disease (Ellenbroek and Youn, 2016; Dawson et al., 2018), normal lifespan changes in brain volume remain incompletely characterized in rodent models.
This incomplete characterization reflects, in part, evolving methodology such as specialized small-rodent scanners and non-standardized brain segmentation methods (Denic et al., 2011; Sawiak et al., 2013), that has contributed to a lag in the use of structural MRI in preclinical rodent models compared to its widespread use in human research (Feo and Giove, 2019). As a result, there are few lifespan brain volume studies in rodents. Whole-brain volumetric findings suggest that gray matter volume increases in rodents before plateauing at 2 months (Mengler et al., 2014). A recent study describes prominent regional volume loss in the isocortex (i.e., often referred to as neocortex in humans), most notably in the prefrontal cortex and temporal association areas, as well as the insula and perirhinal area. Further loss occurred in the cerebellum, while increases in volume occurred in somatosensory areas, thalamus, midbrain, septal nuclei, and specific hippocampal regions (Alexander et al., 2020). The sparse literature on lifespan brain volume changes in rodents poses a gap in the determination of their suitability as a model of healthy structural brain aging.
In this context, specific brain regions warrant particular attention due to their selective age-sensitivity and relevance to age-related cognitive phenotypes. These regions of interest (ROIs) include the anterior cingulate cortex (ACC), hippocampus, and orbitofrontal cortex (OFC). In humans, the ACC connects to both prefrontal and limbic regions, integrating emotional and cognitive functions, including attention, decision making, fear response, and memory (Bush et al., 2000). In rodents, the homologous anterior cingulate area subserves similar functions to those of the ACC (Frankland et al., 2004; Hosking et al., 2014; Zhang et al., 2014), comprising one of four subregions of the medial prefrontal cortex (Heidbreder and Groenewegen, 2003). Normal aging is largely associated with volume loss in the ACC (Good et al., 2001; Bergfield et al., 2010), although preservation has also been reported (Raz, 1997). Further evidence of aging-dependent changes in the ACC includes reductions in spine density and dendritic tree extent (Markham and Juraska, 2002), and decreases in glucose uptake and metabolism (Pardo et al., 2007).
The hippocampus is another region of particular relevance due to its sensitivity to aging-related changes. By integrating sensory, spatial, and temporal information, it contributes to memory formation, learning, and spatial navigation (Eichenbaum, 1999; Maguire et al., 1999). In humans, drastic hippocampal volume loss is observed in aging (Jernigan et al., 2001; Zheng et al., 2019), however, conflicting results have also been reported (Sullivan, 1995; Good et al., 2001). Furthermore, this region displays a unique pattern of continual neurogenesis and increase in brain volume into adulthood in both humans and rodents (Diamond et al., 1975; Sullivan, 1995; Sullivan et al., 2006), while reduced hippocampal neurogenesis in late life has been associated with decreased volume and memory performance (Kuhn et al., 1996; Driscoll et al., 2006; Bettio et al., 2017).
Lastly, the OFC in humans is a prefrontal cortical region that acts as a link between sensory integration, autonomic reactions, learning, and decision making for emotional and reward-related behaviors (Kringelbach, 2005). This age-sensitive region is of interest due to evidence of decreasing volume across aging (Convit et al., 2001; Tisserand et al., 2002). In rodents, the orbital area is homologous to a human agranular section of the OFC, serving similar function (Öngür and Price, 2000; Izquierdo, 2017). Better understanding the nature of age-dependent structural changes in these regions is likely to provide insight into age-dependent cognitive changes.
Network-based approaches to modeling brain structure can complement univariate analyses of regional brain volume to uncover mechanisms underlying normal aging. Structural covariance networks (SCNs) are comprised of nodes (brain regions) connected by edges which represent the statistical correlation between properties of adjacent nodes (Bullmore and Sporns, 2009). As such, they reflect synchronized changes in morphology (volume, cortical thickness) between brain regions (Alexander-Bloch et al., 2013), thought to be driven by mutually trophic influences (Ferrer et al., 1995; Pezawas et al., 2004) and experience-driven plasticity (Mechelli, 2005). SCNs are heritable (Mechelli, 2005), occur at the individual and population level (Tijms et al., 2012; Alexander-Bloch et al., 2013), and partially map on to structural and functional connectivity (Alexander-Bloch et al., 2013). SCNs also display consistent lifespan trajectories, reflecting coordinated neurodevelopmental changes across maturation and aging (Fan et al., 2011; Khundrakpam et al., 2013; Cao et al., 2016). Furthermore, perturbations to typical composition of SCNs have been observed in neurodegenerative diseases and psychiatric illnesses (Alexander-Bloch et al., 2013; Prasad et al., 2022).
While SCN analysis could provide additional insight into age-related changes in brain structure and organization, few studies have considered lifespan SCN trajectories in rodents. Existing literature has employed scaled sub-profile modeling SCN methods (Alexander et al., 2020), but none, to our knowledge, have leveraged whole-brain SCN and graph theoretical methods. The latter can be used to characterize different aspects of SCN topology, from which inferences regarding brain structure and function can be drawn (Bullmore and Sporns, 2009). Topological measures broadly index network (1) node centrality, identifying the “hubness” of highly connected and influential nodes; (2) integration, capturing the capacity for efficient network-wide communication and information processing; and (3) segregation, quantifying the presence of densely interconnected and clustered nodes facilitating local or specialized information processing (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010).
In the present study, we aim to use structural MRI and graph theoretical methods to provide a comprehensive characterization of lifespan brain volume changes in mice. By comparing young mice to middle-aged mice, and middle-aged mice to old mice, we distinguish between maturational and aging-related differences, respectively, in regional brain volume for 32 brain regions. Particular consideration is given to the three aforementioned a priori ROIs, for which there is evidence of age-sensitivity and relevance to age-related cognitive phenotypes (anterior cingulate area, hippocampal formation, and orbital area). We then map these regional brain volume changes to cognitive outcomes by establishing associations with working memory performance on the Y-maze task. Lastly, we further explore regional volumetric changes by analyzing maturational and aging-related differences in SCNs, examining network topology across the mouse lifespan.
2. Materials and methods
2.1. Animals
Two months-old (young), 10 months-old (middle-age) and 22 months-old (old) male C57BL/6 mice (n = 10/group, Charles River Laboratories, Senneville, Quebec, Canada) were housed under a 12/12-h light/dark cycle (7am–7pm light phase) at constant temperature (20–23°C) with free access to food and water. One mouse in the 22 months-old animal group was excluded from the study for developing cataract signs before behavioral testing. Testing was performed by an experimenter blinded to the animal group during the animal’s light cycle phase. Animal use and testing procedures were conducted in accordance with the Canadian Council for Animal Care with approval from the institutional animal care committee.
2.2. Behavioral testing
2.2.1. Y-maze
Testing was performed as in Prevot et al. (2019, 2021). The Y-maze was chosen in preference over other working memory tasks (e.g., Morris water maze or radial arm maze) as it is relatively less susceptible to confounding factors related to aging, such as weight, fatigue, and mobility differences (Matzel et al., 2008). The Y-maze apparatus consists of 3 identical black PVC arms (26 cm length × 8 cm width × 13 cm height) in the shape of a Y; each arm has a sliding door. Distal cues of different shape, form and color were placed on the walls around the room. Mice were allowed to freely explore the apparatus for 10 min on 2 consecutive days (habituation stage). The next day, mice performed a training session consisting of seven consecutive trials. Each trial began with placing the animal in the starting arm for 5 s, prior to the door of y-maze arm opening, allowing the animal to freely choose between the 2 goal arms, and choices were recorded. Upon entering the arm, the door was closed for a 30 s period (inter-trial interval or ITI) after which the animal was returned to the starting arm for the next trial (training stage). The following day, mice were subjected to a similar procedure as the training session, with the exception that the ITI being lengthened to 60 s ITI (testing stage). The mean percent alternation rate (number of alternation/number of trials × 100) was calculated as an index of working memory. ANOVA was performed to determine differences between groups, followed by the Wald test in a logistic regression model, to account for the non-normality of the data distribution and the non-continuous nature of the proportion alternation scores.
2.2.2. Locomotor activity
During the light cycle, locomotor activity of each mouse was measured for 1 h in an observational cage (30 cm × 30 cm) with bedding on the floor (Noldus Phenotyper®, Leesburg, VA, USA). Using Noldus EthoVision 10 tracking software, mice total distance traveled within the cage was analyzed.
2.3. MRI data collection and pre-processing
Twenty-four h after the last behavioral testing, animals were anesthetized and perfused with 4% paraformaldehyde containing 2 mM of ProHance (gadoteridol). Brains were then prepared and scanning was performed as described in Nikolova et al. (2018). In this case, however, we extracted brain volumes from a total of 280 regions (instead of 159, as in Wheeler et al., 2015), excluding ventricles. This parcellation method was adapted from previous rodent MRI studies (Dorr et al., 2008; Richards et al., 2011; Ullmann et al., 2013; Steadman et al., 2014). Deformation-based morphometry was used to derive absolute brain volumes (in mm3). For all MRI analyses, ROI volumes were averaged across hemispheres in order to minimize multiple testing burden.
2.4. Data reduction
MRI-derived volumes were extracted from 280 regions. To reduce the number of comparisons, these regions were combined according to their positions within a hierarchical tree, ultimately yielding 32 regions of interest (ROIs) (Refer to Supplementary Material 1). Of these ROIs, 5 are white matter tracts (cerebellar-related fiber tracts, cranial nerves, extrapyramidal fiber system, lateral forebrain bundle system, and medial forebrain bundle system) and 27 are gray matter divisions. The latter includes 5 brainstem areas (hypothalamus, medulla, midbrain, pons, and thalamus); 2 cerebral nuclei (pallidum, striatum); 18 cerebral cortical areas, including the cortical subplate, hippocampal formation, olfactory areas, and 15 isocortical areas (agranular insular area, anterior cingulate area, auditory areas, ectorhinal area, frontal pole, infralimbic area, orbital area, perirhinal area, posterior parietal association areas, prelimbic area, retrosplenial area, somatomotor areas, somatosensory areas, temporal association areas, visual areas); and 2 cerebellar divisions (cerebellar cortex, cerebellar nuclei).
2.5. Volumetric analyses
2.5.1. A priori ROI volumetric effects of age
Preliminary statistical analyses focused on the anterior cingulate area, hippocampal formation, and orbital area. For each ROI, we performed the following analyses. First, we tested pairwise group differences (young vs. middle-age, middle-age vs. old) in volume using separate general linear models including total brain volume as a covariate. Then, we tested the sample-wide associations of regional volume, residualized for total brain volume, with working memory performance. Working memory performance was evaluated using the Y maze task, and logistic regression was employed with volume measurement as a continuous independent variable and percent alteration expressed as a proportion score, bound between the values of 0 and 1, as the dependent variable (Long, 1997). To account for multiple comparisons, raw p-values from each analysis (n = 6, 3, respectively) were adjusted using the Bonferroni-Holm method (Holm, 1979), a correction commonly used in MRI analysis (Kerr et al., 2022; Straub et al., 2023). We opted for this relatively conservative correction method to ensure greater specificity of our findings. Notably, however, the Bonferroni-Holm method is less conservative than the Bonferroni correction, offering similar protection against Type I error coupled with reduced likelihood of Type II error (Eichstaedt et al., 2013).
2.5.2. Whole-brain volumetric effects of age
We repeated the aforementioned analyses brain-wide, correcting for multiple testing in all 32 ROIs using the Bonferroni-Holm method.
2.6. Structural covariance network analyses
2.6.1. A priori degree centrality analyses
We used the R package “igraph” to generate and compare age group-specific structural covariance matrices indexing correlations of volume between brain regions. Regional volumes were first adjusted for total brain volume, and all negative correlations were removed [to align with previous work (Nikolova et al., 2018; Misquitta et al., 2021)], resulting in group-specific structural covariance matrices that were unweighted and unsigned. For each group, a set of structural covariance networks was then defined by thresholding the covariance matrices using a range of density thresholds (0.05–0.40), indexing the top 5–40% strongest connections in sequential increments of 1%. We then tested pairwise group differences in degree centrality, a node’s level of interconnectedness within the network (Bullmore and Sporns, 2009), for each of the 3 a priori ROIs (anterior cingulate area, hippocampal formation, and orbital area) across each of the above-mentioned age group-specific structural covariance networks. For statistical analysis, permutation testing was performed at each density threshold (n = 5,000 permutations), yielding null distributions against which empirical two-tailed p-values were computed. Significance was defined as p < 0.05 at 5 or more consecutive density thresholds.
2.6.2. Brain-wide degree centrality analyses
We repeated the degree centrality analyses for all 32 brain regions. Permutation testing was performed at density thresholds using n = 1,000 permutations.
2.6.3. Global structural covariance network analyses
In addition, we conducted pairwise group comparisons (young vs. middle-age, middle-age vs. old) of modularity and transitivity, network segregation measures indexing densely internally connected sub-graphs within the network, and the tendency of nodes to form clusters, respectively; and mean distance, the average length of shortest paths between nodes, reflecting overall network efficiency (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010). Permutation testing was performed across a range of density thresholds (0.05–0.40), using n = 5,000 permutations. To further explore group differences in network modularity, a Walktrap community detection algorithm from the “igraph” package was used to identify densely connected sub-graphs within age-group specific networks. For each age group, the most prominent provincial hubs and connector hubs were identified in communities by calculating both the within-community degree z-score (the number of connections a given node has to other nodes within the same community) and the participation coefficient (the distribution of a node’s connections to other communities in the network) of nodes, as outlined in Guimerà et al. (2005) and adapted from Shizuka (2019). Provincial hubs drive exchange and integration of information within a single segregated community. As such, they are defined by a high within-community degree z-score and a low participation coefficient. In contrast, connector hubs facilitate sharing and integration of information between otherwise segregated communities and are therefore defined by both a high within-community degree z-score and a high participation coefficient (Sporns, 2016).
3. Results
3.1. Cognitive performance on Y-maze task and brain volume
Although old mice showed average percent alternations close to 50% chance compared to that of 70% in young animals, we found no significant difference in the percent of spontaneous alternation between young, middle-age, or old animal groups (ANOVA F2,26 = 1.96, p = 0.16, Figure 1A). Follow up analysis using the Wald test in a logistic regression model similarly found no significant difference in proportion of alternation between age groups (χ2 = 4.7, p = 0.095). We did not detect any statistical difference between groups in distance traveled in the locomotor activity test (F1,26 = 2.240, p = 0.117; Young: 6,311cm ± 133 cm; Middle-age: 8,437 cm ± 795cm; Old: 7,215 cm ± 1,141 cm). In the absence of between-group differences in cognitive performance, we then examined potential associations between working memory measured in the Y-maze test and brain volume of our 3 a priori ROIs and brain-wide (i.e., across all 32 regions) in the full sample. Independent of age group, better cognitive performance was associated with greater volume of the anterior cingulate area (χ2 = 2.325, pBH = 0.044, Figure 1B) and smaller volume of the hippocampal formation (χ2 = −2.180, pBH = 0.044, Figure 1C), while no association emerged with volume of the orbital area (χ2 = 0.855, pBH = 0.393). Significant brain-behavior associations are depicted in Figure 1, and corresponding test statistics are summarized in Table 1. At the whole-brain level, no association between Y-maze performance and regional volume survived correction for multiple comparisons. However, we identified several associations, wherein a “younger” brain phenotype was associated with better cognitive performance, regardless of chronological age (i.e., across age groups). Specifically, better cognitive performance was associated with larger volume of the agranular insular area (χ2 = 2.256, p = 0.024) and ectorhinal area (χ2 = 2.646, p = 0.008), both of which were larger in younger relative to middle-aged mice. Conversely, better cognitive performance was also associated with lower volume of the medial forebrain bundle system (χ2 = −1.968, p = 0.049), which was larger in older age.
Figure 1. Significant associations between brain volume and working memory performance. Percentage of correct alterations in the Y-maze task are depicted for the young, middle-age and old mice, with error bars representing standard error of the mean (SEM) (A). Greater volume of the anterior cingulate area was associated with better working memory performance on the Y-maze task (B). Smaller volume of the hippocampal formation was associated with better working memory performance on the Y-maze task (C). Associations survived correction for multiple comparisons across these a priori regions (p-value pictured) but not across brain-wide regions.
3.2. Volumetric changes across the lifespan
3.2.1. A priori ROI volumetric effects of age
We detected significant main effects of age on regional volume of the anterior cingulate area and orbital area, but not of the hippocampal formation. Specifically, we found a maturational decrease in volume in the anterior cingulate area, which was smaller in middle-aged mice than in young mice (t = −4.898, pBH = 0.001) (Figure 2A), but did not differ between middle-aged and old mice. We observed a different trajectory of volume loss in the orbital area, where volume was significantly lower in old mice than in middle-aged mice (t = −2.859, pBH = 0.038; Figure 2B), but did not differ between young and middle-aged mice, indicating aging-related volume loss. Results of all pairwise comparisons are summarized in Table 2.
Figure 2. Volumetric changes in maturational and aging-related comparisons. Negative and positive linear changes in regional brain volume are depicted between pairwise comparisons of young and middle-age mice (maturational) (A) and middle-age and old mice (aging-related) (B). Brain regions that underwent significant change in volume are denoted by asterisk (*).
3.2.2. Whole-brain volumetric effects of age
We found significant group differences in regional volume of 24 of 32 regions tested, including some of the a priori ROIs, as mentioned above (Figures 2A, B). In most of these regions (18/24), volume differed significantly between young mice and middle-aged mice, indicating a maturational effect. The maturational effects were characterized by increases in volume (i.e., middle-aged > young) and decreases in volume (i.e., middle-aged < young) in a roughly equivalent number of regions. Maturational increases in volume (t = 2.611–6.837, pBH ≤ 0.050) were observed in ten regions, including the cerebellar nuclei, pallidum, most of the brainstem areas (hypothalamus, medulla, midbrain, pons), and most of the white matter tracts (cerebellum-related fiber tracts, cranial nerves, lateral forebrain bundle system, medial forebrain bundle system) (Figure 2A). Maturational decreases in volume (t = 2.674–4.898, pBH ≤ 0.050) were observed in 10 regions, including the cerebellar cortex and 7 of 15 isocortical areas: anterior cingulate area, auditory areas, prelimbic area, somatosensory areas, somatomotor areas, temporal association areas, visual areas (Figure 2A).
Aging-related (i.e., middle-aged < old) increases in volume were observed in the lateral forebrain bundle system (t = 3.275, pBH = 0.020) and medial forebrain bundle system (t = 3.177, pBH = 0.022), white matter tracts in which we also found maturational increases in volume. Aging-related decreases (i.e., middle-aged > old) in volume were observed in two distinct isocortical areas, the ectorhinal cortex (t = −3.199, pBH = 0.022) and the orbital area (t = −2.859, pBH = 0.038) (Figure 2B). Results of all pairwise comparisons are summarized in Table 2.
3.3. Structural covariance network adaptations
3.3.1. A priori ROI degree centrality analyses
Among the 3 a priori ROIs, we found significant maturational group differences in regional degree centrality of the anterior cingulate area and hippocampal formation, but not of the orbital area (Figures 3A-C). Degree centrality in the anterior cingulate area was significantly lower in middle-aged mice than in young mice (23–40% density) indicating a maturational decrease in regional connectivity (Figures 3A, 4A-C). Conversely, degree centrality in the hippocampal formation was significantly higher in middle-aged mice than in young mice (8–16% density) indicating a maturational increase in regional connectivity (Figures 3B, 4D-F). In contrast, no significant degree centrality differences in a priori regions were detected in the aging-related comparisons (Figures 3D-F). Results of all pairwise comparisons are summarized in Supplementary Table 1.
Figure 3. Regional degree centrality differences in a priori ROIs for maturational and aging-related comparisons. Maturational comparisons in structural covariance networks for the a priori age-sensitive regions revealed significant decrease in degree centrality for the anterior cingulate area (A), a maturational increase in degree centrality in the hippocampal formation (B), and no significant differences observed in the orbital area (C). Aging-related comparisons in structural covariance networks for the a priori age-sensitive regions did not yield significant differences in degree (D-F). Differences in degree were determined across density thresholds of 0.05 and 0.40 at 5 consecutive significant density thresholds [statistical significance indicated by white circles (*p < 0.05)].
Figure 4. Circular degree plots of significant differences in regional connectivity in a priori ROIs. The anterior cingulate area networks in young (A), middle-age (B), and old (C) mice are depicted at network density of 24% for visualization. The hippocampal formation networks in young (D), middle-age (E), and old (F) mice are depicted at network density of 12% for visualization. Nodes are arranged clockwise by degree, the number of significant correlations a given brain region has with other regions within the network, with the highest-degree node at the top of the circle. Degree is further reflected in node size. Rank refers to the numerical value of a region’s degree, sorted largest to smallest in the network. The region of interest and its direct connections are highlighted in red.
3.3.2. Brain-wide degree centrality analyses
We identified group differences in degree centrality of 9 of 32 regions tested. In all of these regions, degree differed significantly between young mice and middle-aged mice, indicating maturational network connectivity changes. In addition to the above-mentioned maturational differences in degree centrality of the anterior cingulate area (decrease) and hippocampal formation (increase), we also detected significant maturational decreases in 3 additional isocortical regions (auditory areas, ectorhinal area, somatomotor areas) and significant maturational increases in 4 regions (cerebellar nuclei, cerebellum-related fiber tracts, cortical subplate, midbrain) across partially overlapping density thresholds (Figures 5A-G). No significant differences in regional degree were found between the middle-aged and old mice.
Figure 5. Brain-wide degree centrality differences in maturational and aging-related comparisons. Brain-wide maturational comparisons in structural covariance networks revealed significant decreases in degree in the auditory areas (A), ectorhinal area (B), and somatomotor areas (C), and maturational increases in degree in the cerebellar nuclei (D), cerebellum related fiber tracts (E), cortical subplate (F), and midbrain (G). No significant differences were found in any brain regions in the aging-related network comparisons. Changes in degree were determined across density thresholds of 0.05–0.40 at 5 consecutive significant density thresholds [statistical significance indicated by white circles (*p < 0.05)].
3.3.3. Global structural covariance network analyses
We detected a significant maturational decrease (i.e., middle-age < young) in global network modularity (density 24–29%) (Figure 6A), but no significant aging-related differences therein (Figure 6B). In addition, we did not detect any significant maturational or aging-related differences in global network transitivity or mean distance measures. Results of all pairwise comparisons are summarized in Supplementary Table 2.
Figure 6. Network modularity comparisons. A significant maturational decrease in modularity was observed between young and middle-age mice networks (A). Modularity did not differ significantly between middle-age and old mice networks (B). Changes in modularity were determined across density thresholds of 0.05 and 0.40. 5 consecutive significant density thresholds [statistical significance indicated by white circles (*p < 0.05)].
Network community modules generated by the Walktrap community detection algorithm revealed densely interconnected sub-graphs in the young, middle-age, and old networks (Figures 7A-C). Communities in the young network had densely connected within-community edges (black lines) and few inter-community connections (red lines), while the middle-age and old network communities had progressively fewer within-community connections and more inter-community connections. In the young mice network, the most prominent provincial hubs of their respective communities were the frontal pole, extrapyramidal fiber system, and olfactory areas, while the most prominent connector hubs were the ectorhinal area, temporal association areas. Moreover, in the middle-aged network, the frontal pole, posterior parietal association areas, hypothalamus, pallidum, and lateral forebrain bundle system were identified as provincial hubs, and the temporal association areas, cerebellar nuclei, and pons were identified as connector hubs. Lastly, in the old mice network, the midbrain was identified as the only provincial hub, while the pallidum, prelimbic area, agranular insular areas, and hippocampal formation were identified as connector hubs.
Figure 7. Network communities in mice. Communities of densely interconnected nodes were identified using a Walktrap community detection algorithm for the structural covariance networks of the young (A), middle-age (B), and old mice (C). Regions of interest are sorted by isocortical (green), subcortical (purple), white matter (blue), brainstem (orange), cerebellum (yellow), and cerebellar nuclei (red) regions. Provincial hub regions are denoted by “P,” and connector hub regions are denoted by “C” in the young “Y,” middle-age “M” and old “O” networks.
4. Discussion
Leveraging structural neuroimaging data in young, middle-aged, and old mice, we provide a comprehensive characterization of regional volume and SCN changes across the lifespan, and map them onto working memory performance in the Y-maze paradigm. We observed robust and widespread maturational changes in brain volume, including notable volume loss in isocortical regions accompanied by volume gain in white matter tracts and the brainstem. Aging-related volume changes were more discreet, with further volume loss in select regions of the isocortex, including the orbital area, ectorhinal area, auditory areas, and temporal association areas, coupled with volume gain in white matter tracts, and contrasted by relative preservation in subcortical and cerebellar regions. Additionally, poorer working memory performance on the Y-maze task was associated with smaller volume of the anterior cingulate area and larger volume of the hippocampal formation, regardless of age. Lastly, SCN analyses revealed widespread maturational changes in node-level degree centrality and a significant decrease in network modularity, but not in the aging-related comparison. Network transitivity and mean distance remained relatively stable across the mouse lifespan.
4.1. Brain volume changes across maturation and aging in mice
In the present study, the most prominent maturational changes in volume (young vs. middle-aged mice) occurred in the isocortex. This finding mirrors those in human literature, wherein protracted development in the cerebral cortex precedes widespread gray matter loss, most notably in the prefrontal cortex and to a lesser extent in the parietal and temporal association areas (Raz, 1997; Tisserand et al., 2002; Resnick et al., 2003). Similar to these trends, we detected significant volume loss in 7 of the 15 isocortical regions, including 2 of the 3 “prefrontal” regions in rodents (prelimbic area, anterior cingulate area) (Laubach et al., 2018) and the temporal association area. We did not observe significant age-related changes in the remaining prefrontal region (infralimbic area) or in the parietal association area. The aforementioned volume loss in the anterior cingulate area, an a priori ROI, aligns with previous rodent literature (Hamezah et al., 2017), along with most human findings (Raz, 1997; Good et al., 2001; Bergfield et al., 2010).
Aging-related volume changes, indexed by comparing middle-age vs. old mice, were also associated with considerable volume loss in the isocortex, with atrophy detected in 4 of the 15 associated regions. Notably, the temporal association area experienced both maturational and aging-related loss, suggesting lifelong volumetric decrease, as did the auditory areas. Moreover, the orbital area, another a priori ROI, experienced significant aging-related loss consistent with human literature (Convit et al., 2001; Tisserand et al., 2002), suggesting that this region may be particularly sensitive to late-life volume decrease in mice.
Beyond the isocortex, other areas demonstrated similar patterns to human trajectories of lifespan brain volume. Specifically, we found significant volume increase in 4 of 5 brainstem regions in the maturational comparison. Furthermore, these regions did not experience significant volume loss or gain in the aging-related comparison. These findings align with human trajectories of increasing brainstem volume in maturation (Walhovd et al., 2005), and accounts of relative stability in aging (Luft et al., 1999; Walhovd et al., 2011; Long et al., 2012). In addition, the maturational comparison of white matter tracts in the mice revealed 4 of 5 tracts significantly increasing in volume. This aligns with previous rodent studies (Calabrese and Johnson, 2013) as well as human trends of white matter gradually maturing into adulthood (Jernigan et al., 2001; Liu et al., 2016).
Although in the present study we described cortical maturational and aging-related volume changes in mice similar to those observed in humans, we also found deviations from human lifespan trajectories. Most notably, we did not detect significant maturational or aging-related volumetric changes in the 3 subcortical structures, including the hippocampal formation, an a priori ROI, or in cerebellar regions. This is contrary to evidence from human studies, which suggest that subcortical structures reach maximum volume early in development before undergoing a steep decline in the sixth decade (Dima et al., 2022). It is similarly contrary to evidence of continued hippocampal formation neurogenesis and growth into adulthood in both humans and rodents (Diamond et al., 1975; Sullivan, 1995) and of aging-related volumetric decline of the cerebellum in humans (Escalona et al., 1991; Jernigan et al., 2001; Raz et al., 2001). Although the absence of aging-related decrease in the hippocampal formation differs from accounts of drastic late-life volume loss in humans and rodents (Jernigan et al., 2001; Alexander et al., 2012; Hamezah et al., 2017; Zheng et al., 2019), it does align with contradictory findings of preservation in humans (Sullivan, 1995; Good et al., 2001). Moreover, it is possible that segmentation of the hippocampal formation may have increased sensitivity to changes in volume (Hackert et al., 2002). Finally, the observed significant aging-related volume increases in 2 of 5 white matter tracts differs from findings of drastic late-life white matter loss in human literature (Jernigan et al., 2001; Liu et al., 2016).
Overall, the volumetric changes found in the present study were more numerous and robust in the maturational comparisons than the aging-related comparisons. We speculate this may be attributable to the relatively shorter late-life period in mice compared to the protracted aging period humans experience. Furthermore, the controlled laboratory setting may have been protective against perturbations caused by adverse environmental variables contributing to aging-related decline. Lastly, it is also possible that our study was underpowered to detect aging-related changes due to small group size and greater inter-subject (i.e., significant within-group) variability.
4.2. Associations between working memory performance and regional brain volume
Poorer working memory performance in the Y-maze was associated with smaller volume of the anterior cingulate area. Notably, this region also underwent a significant maturational decrease in volume and degree centrality. The anterior cingulate area has previously been implicated in working memory in rodents (Teixeira et al., 2006), similar to the involvement of the ACC in working memory in humans (Kaneda and Osaka, 2008), including aging-related decline in functional activity of the ACC in working memory performance (Pardo et al., 2007). Furthermore, as one of three prefrontal regions in rodents, this association also supports the relationship between the prefrontal cortex and working memory decline (Head et al., 2002). In the a priori ROI analyses, poorer working memory performance was also associated with larger volume of the hippocampal formation, despite no significant volume differences identified over the lifespan. This association further supports the well-known involvement of the hippocampal formation in working memory in humans, and rodents (Yoon et al., 2008; Yonelinas, 2013). A meta-analysis of volumetric changes in the hippocampus and memory outcomes in humans demonstrated high variability in the structural associations in older adults, while youth and young adults trended toward worse performance with higher hippocampus volume (Van Petten, 2004). The association observed in our mice is difficult to interpret against the variable human literature, but notably does not align with rodent findings of smaller hippocampal formation volume and poorer working memory performance (Kadar et al., 1990; Hamezah et al., 2017).
We observed additional associations that did not survive multiple corrections. These included the agranular insular area and ectorhinal area, two additional isocortical regions where poorer working memory performance was associated with smaller brain volume. This corroborates the agranular insular area’s involvement in memory processes in rodents, including working memory (DeCoteau et al., 1997; Zhu et al., 2020). The ectorhinal area, which also experienced significant aging-related volume loss, borders the perirhinal cortex and comprises the posterior of the parahippocampal gyrus (Westerhaus and Loewy, 2001). The ectorhinal area is also implicated in memory (Tulving and Markowitsch, 1997), and has been observed to be sensitive to aging-related volume loss in humans (Raz and Rodrigue, 2006). Lastly, we found that poorer working memory performance was associated with larger volume of the medial forebrain bundle system, a white matter tract implicated in motivation and reward (Wise, 2005; Anthofer et al., 2015). While white matter volume generally shows less associations with cognition compared to those found in gray matter (Taki et al., 2011), it is notable that the significant aging-related volume increase in this region deviates from human patterns of white matter loss (Liu et al., 2016). Overall, these associations reflected a general trend of better working memory performance from structurally “younger” brain regions, and declining cognitive performance when regional brain structures were characteristic of older mice.
Although we observed numerous volumetric changes, few brain regions showed associations with working memory performance. This may be due to the relatively discreet significant regional volume changes between middle-age and old mice. It is also possible that older mice utilize cognitive reserve by employing alternative cognitive strategies to compensate for deficits through differential recruitment patterns of brain activity and structural reorganization (Stern, 2009; Eyler et al., 2011; Zatorre et al., 2012). Lastly, it is notable that aging-related decline in cognitive function is also driven by factors other than brain volume loss, including diminished synaptic, neurotransmitter, cell signaling, and mitochondrial function (Shankar, 2010).
4.3. Structural covariance network topology across maturation and aging in mice
Investigating SCN changes allowed us to consider structural synchronization suggestive of interactions among and between multiple brain regions. The existing SCN literature in rodents is limited (Alexander-Bloch et al., 2013), and to our knowledge no studies have explored lifespan SCN trajectories in mice using whole-brain SCN methodology. In our maturational comparison, 9 of the 32 regions experienced significant changes in degree centrality, with 5 regions in young mice experiencing initially highly connected and influential nodes in the network decrease in centrality into middle-age, including the anterior cingulate area. Four regions emerged as significantly connected nodes in the middle-age network, including the hippocampal formation. These differences in degree centrality largely coincided with significant changes in volume, whereby 3 of the 4 maturational degree centrality decreases co-occurred with significant maturational volume loss, and 3 of the 5 maturational increases in degree centrality accompanied significant maturational volume increase. This may reflect a general trend of higher or lower influence of nodes within the network being associated with directionally consistent volume changes. In contrast, the comparison between middle-age and old mice yielded no significant degree centrality differences across all 32 regions. These findings also offer insight into how mouse SCNs compare to established human lifespan SCN trajectories. In humans, the transition from childhood and adolescence to young adulthood brings a generalized homogenization of covariance patterns, as networks stabilize toward adult topology (Zielinski et al., 2010). Similarly, a recent large-scale study described pronounced stabilization of overall structural covariance by age 25 (Nadig et al., 2021). The widespread maturational differences in regional degree we observed may reflect developmental reorganization of the SCNs in the young mice, prior to stabilizing into adult SCN patterns in the middle-aged mice. The lack of degree centrality differences between middle-aged and old mice suggest that network stability persisted as mice got older.
In addition to regional connectivity, we investigated global measures of network integration and segregation to further characterize lifespan SCNs in mice. The lack of differences in mean distance across the lifespan of the mice may indicate that mice do not experience late-life reduction in network efficiency in normal aging that is typically observed in human aging (Khundrakpam et al., 2013; Betzel et al., 2014). Furthermore, the lack of maturational or aging-relating differences in network transitivity suggests clustering in mouse SCNs is relatively stable across the lifespan. However, we did observe a significant maturational decrease in modularity, whereby young mice had higher network modularity than middle-aged mice, but there were no differences in modularity in the aging-related comparison. Human maturational SCN trajectories demonstrate highly segregated networks favoring localized connections in youth and adolescence, before shifting toward more distributed and globally connected networks in young adulthood (Khundrakpam et al., 2013; Cao et al., 2016). The observed maturational decrease in modularity suggests mice may share a similar trajectory with humans, whereby network segregation decreases from youth to young adulthood. In older age, human SCN segregation has been posited to decrease, becoming further dispersed and disorganized (Chen et al., 2011; Jao et al., 2020), however, conflicting findings suggest that network segregation may actually increase in old age (Montembeault et al., 2012). In particular, modularity in structural and functional networks has drawn attention as a potential biomarker for aging (Aboud et al., 2018) as decreasing modularity may be associated with aging-related dedifferentiation of specialized cognitive processes (Chen et al., 2011; Onoda and Yamaguchi, 2013; Song et al., 2014). The observed lack of aging-related modularity loss suggests that normal mice aging may differ from SCN segregation trajectories present in human aging. Furthermore, this finding may also reflect the relatively discreet volumetric changes in our aging-related comparison.
To further probe our findings in network modularity, we explored concomitant changes in densely interconnected sub-graphs within the networks. The significant maturational decrease in modularity detected in our whole-brain analysis was reflected in the sub-graphs, as the high within-community connections and low inter-community connections in the young mice communities shifted to sparser within-community and increased inter-community connections in the middle-age communities. While the whole-brain aging-related comparison for modularity was not significant, the old-age communities more closely resembled an “old network phenotype,” with few within-community connections, high inter-community connections, and smaller, more numerous communities. Among the identified provincial hubs, the frontal pole emerged as the provincial hub for a community containing predominantly isocortical regions in both young and middle-age networks, suggesting its importance in facilitating specialized information processing among regions associated with higher-order sensory and cognitive processes. Only one provincial hub, the midbrain, emerged in the old-network communities, reflecting the overall reduced within-community connections. Lastly, the temporal association areas and pons were identified as connector hubs in both young and middle-aged networks, highlighting their importance in facilitating information exchange between segregated network communities.
Certain limitations in this study should be considered. Most notably, it is possible that that we were underpowered to detect volumetric changes, particularly between the middle-age and old groups, given our relatively small sample size. Future work to validate and further explore regional volume trajectories across the mouse lifespan will benefit from larger sample sizes in each age-group. Secondly, the study utilized a cross-sectional design, considering different groups of mice at each timepoint. Ideally, future studies could employ a longitudinal design, observing the same group of mice over time, thereby increasing the strength of inferences drawn from the maturational and aging-related comparisons. Furthermore, only male mice were available at the time of this study, therefore our results are not generalizable to female mice. Differences in gray matter volume have been observed between male and female mice, both in speed of gray matter development and variation in regional volume (Spring et al., 2007; Qiu et al., 2018). It is possible these inherent differences may influence maturational or aging-related brain volume outcomes. Furthermore, differences in brain structure and function may vary across mice strains. Therefore, the observed changes in brain volume in the Charles River C57BL/6 mice may not be fully generalizable to other strains of mice. Lastly, while working memory decline in aging is well-documented, an additional limitation of our work is that we used only one test, the Y-maze, to evaluate working memory. It is also important to acknowledge that multiple cognitive domains are affected by aging, and future work would benefit from exploring additional domains for a more nuanced picture of the association between volumetric changes and cognitive outcomes.
Despite these limitations, this study provides a detailed account of regional brain volume changes across the mouse lifespan and maps them to working memory performance, while also exploring complementary trends in structural covariance network topology. This comprehensive approach may facilitate future progress in translation between rodent models and human brain aging.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by the Centre for Addiction and Mental Health (CAMH), Animal Care Committee (ACC) Canadian institutional animal care committee.
Author contributions
KC: conceptualization, methodology, visualization, formal analysis (including behavioral and MRI data analysis), and writing – original draft. AM: conceptualization, methodology, formal analysis (including behavioral and MRI data analysis), and writing – review and editing. TP: investigation (behavior testing) and writing – review and editing. KM: investigation (behavioral testing). JE: investigation (MRI acquisition). JL: investigation (MRI acquisition). ES: conceptualization and writing – review and editing. YN: supervision, funding acquisition, conceptualization, methodology, and writing – review and editing. MB: funding acquisition, conceptualization, methodology, and writing – review and editing. All authors contributed to the article and approved the submitted version.
Funding
MB was supported by the CAMH discovery fund, the Canadian Institutes of Health Research (PGT165852, PI: MB), and the Campbell Family Mental Health Research Institute. YN was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and also Koerner New Scientist Award, as well as a Paul Garfinkel New Investigator Catalyst Award administered by the Centre for Addiction and Mental Health Foundation.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2023.1195748/full#supplementary-material
References
Aboud, K. S., Huo, Y., Kang, H., Ealey, A., Resnick, S. M., Landman, B. A., et al. (2018). Structural covariance across the lifespan: brain development and aging through the lens of inter-network relationships. Hum. Brain Mapp. 40, 125–136. doi: 10.1002/hbm.24359
Alexander, G. E., Lin, L., Yoshimaru, E. S., Bharadwaj, P. K., Bergfield, K. L., Hoang, L. T., et al. (2020). Age-related regional network covariance of magnetic resonance imaging gray matter in the rat. Front. Aging Neurosci. 12:267. doi: 10.3389/fnagi.2020.00267
Alexander, G., Ryan, L., Bowers, D., Foster, T., Bizon, J., Geldmacher, D., et al. (2012). Characterizing cognitive aging in humans with links to animal models. Front. Aging Neurosci. 4:21. doi: 10.3389/fnagi.2012.00021
Alexander-Bloch, A., Giedd, J. N., and Bullmore, E. (2013). Imaging structural co-variance between human brain regions. Nat. Rev. Neurosci. 14, 322–336. doi: 10.1038/nrn3465
Anthofer, J. M., Steib, K., Fellner, C., Lange, M., Brawanski, A., and Schlaier, J. (2015). DTI-based deterministic fibre tracking of the medial forebrain bundle. Acta Neurochir. 157, 469–477. doi: 10.1007/s00701-014-2335-y
Bergfield, K. L., Hanson, K. D., Chen, K., Teipel, S. J., Hampel, H., Rapoport, S. I., et al. (2010). Age-related networks of regional covariance in MRI gray matter: reproducible multivariate patterns in healthy aging. Neuroimage 49, 1750–1759. doi: 10.1016/j.neuroimage.2009.09.051
Bethlehem, R. A. I, Seidlitz, J., White, S. R., Vogel, J. W., Anderson, K. M., Adamson, C., et al. (2022). Brain charts for the human lifespan. Nature 604, 525–533. doi: 10.1038/s41586-022-04554-y
Bettio, L. E. B., Rajendran, L., and Gil-Mohapel, J. (2017). The effects of aging in the hippocampus and cognitive decline. Neurosci. Biobehav. Rev. 79, 66–86. doi: 10.1016/j.neubiorev.2017.04.030
Betzel, R. F., Byrge, L., He, Y., Goñi, J., Zuo, X.-N., and Sporns, O. (2014). Changes in structural and functional connectivity among resting-state networks across the human lifespan. Neuroimage 102, 345–357. doi: 10.1016/j.neuroimage.2014.07.067
Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Bush, G., Luu, P., and Posner, M. I. (2000). Cognitive and emotional influences in anterior cingulate cortex. Trends Cogn. Sci. 4, 215–222. doi: 10.1016/S1364-6613(00)01483-2
Calabrese, E., and Johnson, G. A. (2013). Diffusion tensor magnetic resonance histology reveals microstructural changes in the developing rat brain. Neuroimage 79, 329–339. doi: 10.1016/j.neuroimage.2013.04.101
Cao, M., Huang, H., Peng, Y., Dong, Q., and He, Y. (2016). Toward developmental connectomics of the human brain. Front. Neuroanat. 10:25. doi: 10.3389/fnana.2016.00025
Chatterji, S., Byles, J., Cutler, D., Seeman, T., and Verdes, E. (2015). Health, functioning, and disability in older adults—Present status and future implications. Lancet 385, 563–575. doi: 10.1016/S0140-6736(14)61462-8
Chen, Z. J., He, Y., Rosa-Neto, P., Gong, G., and Evans, A. C. (2011). Age-related alterations in the modular organization of structural cortical network by using cortical thickness from MRI. Neuroimage 56, 235–245. doi: 10.1016/j.neuroimage.2011.01.010
Convit, A., Wolf, O. T., de Leon, M. J., Patalinjug, M., Kandil, E., Caraos, C., et al. (2001). Volumetric analysis of the pre-frontal regions: findings in aging and schizophrenia. Psychiatry Res. Neuroimaging 107, 61–73. doi: 10.1016/S0925-4927(01)00097-X
Dawson, T. M., Golde, T. E., and Lagier-Tourenne, C. (2018). Animal models of neurodegenerative diseases. Nat. Neurosci. 21, 1370–1379. doi: 10.1038/s41593-018-0236-8
DeCoteau, W. E., Kesner, R. P., and Williams, J. M. (1997). Short-term memory for food reward magnitude: the role of the prefrontal cortex. Behav. Brain Res. 88, 239–249. doi: 10.1016/S0166-4328(97)00044-2
Denic, A., Macura, S. I., Mishra, P., Gamez, J. D., Rodriguez, M., and Pirko, I. (2011). MRI in rodent models of brain disorders. Neurotherapeutics 8, 3–18. doi: 10.1007/s13311-010-0002-4
Diamond, M. C., Johnson, R. E., and Ingham, C. A. (1975). Morphological changes in the young, adult and aging rat cerebral cortex, hippocampus, and diencephalon. Behav. Biol. 14, 163–174. doi: 10.1016/S0091-6773(75)90161-3
Dima, D., Modabbernia, A., Papachristou, E., Doucet, G. E., Agartz, I., Aghajani, M., et al. (2022). Subcortical volumes across the lifespan: data from 18,605 healthy individuals aged 3–90 years. Hum. Brain Mapp. 43, 452–469. doi: 10.1002/hbm.25320
Dorr, A. E., Lerch, J. P., Spring, S., Kabani, N., and Henkelman, R. M. (2008). High resolution three-dimensional brain atlas using an average magnetic resonance image of 40 adult C57Bl/6J mice. Neuroimage 42, 60–69. doi: 10.1016/j.neuroimage.2008.03.037
Driscoll, I., Howard, S. R., Stone, J. C., Monfils, M. H., Tomanek, B., Brooks, W. M., et al. (2006). The aging hippocampus: a multi-level analysis in the rat. Neuroscience 139, 1173–1185. doi: 10.1016/j.neuroscience.2006.01.040
Eichenbaum, H. (1999). The hippocampus and mechanisms of declarative memory. Behav. Brain Res. 103, 123–133. doi: 10.1016/S0166-4328(99)00044-3
Eichstaedt, K. E., Kovatch, K., and Maroof, D. A. (2013). A less conservative method to adjust for familywise error rate in neuropsychological research: the Holm’s sequential Bonferroni procedure. Neurorehabilitation 32, 693–696. doi: 10.3233/NRE-130893
Ellenbroek, B., and Youn, J. (2016). Rodent models in neuroscience research: is it a rat race? Dis. Models Mechan. 9, 1079–1087. doi: 10.1242/dmm.026120
Escalona, P. R., McDonald, W. M., Doraiswamy, P. M., Boyko, O. B., Husain, M. M., Figiel, G. S., et al. (1991). In vivo stereological assessment of human cerebellar volume: effects of gender and age. Am. J. Neuroradiol. 12, 927–929.
Eyler, L. T., Sherzai, A., Kaup, A. R., and Jeste, D. V. (2011). A review of functional brain imaging correlates of successful cognitive aging. Biol. Psychiatry 70, 115–122. doi: 10.1016/j.biopsych.2010.12.032
Fan, Y., Shi, F., Smith, J. K., Lin, W., Gilmore, J. H., and Shen, D. (2011). Brain anatomical networks in early human brain development. Neuroimage 54, 1862–1871. doi: 10.1016/j.neuroimage.2010.07.025
Feo, R., and Giove, F. (2019). Towards an efficient segmentation of small rodents brain: a short critical review. J. Neurosci. Methods 323, 82–89. doi: 10.1016/j.jneumeth.2019.05.003
Ferrer, I., Blanco, R., Carulla, M., Condom, M., Alcántara, S., Olivé, M., et al. (1995). Transforming growth factor-α immunoreactivity in the developing and adult brain. Neuroscience 66, 189–199. doi: 10.1016/0306-4522(94)00584-R
Folgueras, A. R., Freitas-Rodríguez, S., Velasco, G., and López-Otín, C. (2018). Mouse models to disentangle the hallmarks of human aging. Circ. Res. 123, 905–924. doi: 10.1161/CIRCRESAHA.118.312204
Frankland, P. W., Bontempi, B., Talton, L. E., Kaczmarek, L., and Silva, A. J. (2004). The involvement of the anterior cingulate cortex in remote contextual fear memory. Science 304, 881–883. doi: 10.1126/science.1094804
Giedd, J. N., and Rapoport, J. L. (2010). Structural MRI of pediatric brain development: what have we learned and where are we going? Neuron 67, 728–734. doi: 10.1016/j.neuron.2010.08.040
Good, C. D., Johnsrude, I. S., Ashburner, J., Henson, R. N. A., Friston, K. J., and Frackowiak, R. S. J. (2001). A voxel-based morphometric study of ageing in 465 normal adult human brains. Neuroimage 14, 21–36. doi: 10.1006/nimg.2001.0786
Guimerà, R., Mossa, S., Turtschi, A., and Amaral, L. A. N. (2005). The worldwide air transportation network: anomalous centrality, community structure, and cities’ global roles. Proc. Natl. Acad. Sci. U.S.A. 102, 7794–7799. doi: 10.1073/pnas.0407994102
Hackert, V. H., den Heijer, T., Oudkerk, M., Koudstaal, P. J., Hofman, A., and Breteler, M. M. B. (2002). Hippocampal head size associated with verbal memory performance in nondemented elderly. Neuroimage 17, 1365–1372. doi: 10.1006/nimg.2002.1248
Hamezah, H. S., Durani, L. W., Ibrahim, N. F., Yanagisawa, D., Kato, T., Shiino, A., et al. (2017). Volumetric changes in the aging rat brain and its impact on cognitive and locomotor functions. Exp. Gerontol. 99, 69–79. doi: 10.1016/j.exger.2017.09.008
Head, D., Raz, N., Gunning-Dixon, F., Williamson, A., and Acker, J. D. (2002). Age-related differences in the course of cognitive skill acquisition: The role of regional cortical shrinkage and cognitive resources. Psychol. Aging 17:72. doi: 10.1037/0882-7974.17.1.72
Heidbreder, C. A., and Groenewegen, H. J. (2003). The medial prefrontal cortex in the rat: evidence for a dorso-ventral distinction based upon functional and anatomical characteristics. Neurosci. Biobehav. Rev. 27, 555–579. doi: 10.1016/j.neubiorev.2003.09.003
Hosking, J. G., Cocker, P. J., and Winstanley, C. A. (2014). Dissociable contributions of anterior cingulate cortex and basolateral amygdala on a rodent cost/benefit decision-making task of cognitive effort. Neuropsychopharmacology 39, 1558–1567. doi: 10.1038/npp.2014.27
Hou, Y., Dan, X., Babbar, M., Wei, Y., Hasselbalch, S. G., Croteau, D. L., et al. (2019). Ageing as a risk factor for neurodegenerative disease. Nat. Rev. Neurol. 15, 565–581. doi: 10.1038/s41582-019-0244-7
Izquierdo, A. (2017). Functional heterogeneity within rat orbitofrontal cortex in reward learning and decision making. The J. Neurosci. 37, 10529–10540. doi: 10.1523/JNEUROSCI.1678-17.2017
Jao, C.-W., Yeh, J.-H., Wu, Y.-T., Lien, L.-M., Tsai, Y.-F., Chu, K.-E., et al. (2020). Alteration of the intra- and inter-lobe connectivity of the brain structural network in normal aging. Entropy 22:826. doi: 10.3390/e22080826
Jernigan, T. L., Archibald, S. L., Fennema-Notestine, C., Gamst, A. C., Stout, J. C., Bonner, J., et al. (2001). Effects of age on tissues and regions of the cerebrum and cerebellum. Neurobiol. Aging 22, 581–594. doi: 10.1016/S0197-4580(01)00217-2
Kadar, T., Silbermann, M., Brandeis, R., and Levy, A. (1990). Age-related structural changes in the rat hippocampus: correlation with working memory deficiency. Brain Res. 512, 113–120. doi: 10.1016/0006-8993(90)91178-J
Kaneda, M., and Osaka, N. (2008). Role of anterior cingulate cortex during semantic coding in verbal working memory. Neurosci. Lett. 436, 57–61. doi: 10.1016/j.neulet.2008.02.069
Kerr, W. T., Tatekawa, H., Lee, J. K., Karimi, A. H., Sreenivasan, S. S., O’Neill, J., et al. (2022). Clinical MRI morphological analysis of functional seizures compared to seizure-naïve and psychiatric controls. Epilepsy Behav. 134:108858. doi: 10.1016/j.yebeh.2022.108858
Khundrakpam, B. S., Reid, A., Brauer, J., Carbonell, F., Lewis, J., Ameis, S., et al. (2013). Developmental changes in organization of structural brain networks. Cereb. Cortex 23, 2072–2085. doi: 10.1093/cercor/bhs187
Kringelbach, M. L. (2005). The human orbitofrontal cortex: linking reward to hedonic experience. Nat. Rev. Neurosci. 6, 691–702. doi: 10.1038/nrn1747
Kuhn, H., Dickinson-Anson, H., and Gage, F. (1996). Neurogenesis in the dentate gyrus of the adult rat: age-related decrease of neuronal progenitor proliferation. J. Neurosci. 16, 2027–2033. doi: 10.1523/JNEUROSCI.16-06-02027.1996
Laubach, M., Amarante, L. M., Swanson, K., and White, S. R. (2018). What, if anything, is rodent prefrontal cortex? ENeuro 5, ENEURO.315–ENEURO.318. doi: 10.1523/ENEURO.0315-18.2018
Liu, H., Wang, L., Geng, Z., Zhu, Q., Song, Z., Chang, R., et al. (2016). A voxel-based morphometric study of age- and sex-related changes in white matter volume in the normal aging brain. Neuropsychiatr. Dis. Treat. 12, 453–465. doi: 10.2147/NDT.S90674
Long, J. S. (1997). Regression models for categorical and limited dependent variables. Thousand Oaks, CA: Sage Publications Inc, xxx + 297.
Long, X., Liao, W., Jiang, C., Liang, D., Qiu, B., and Zhang, L. (2012). Healthy aging: an automatic analysis of global and regional morphological alterations of human brain. Acad. Radiol. 19, 785–793. doi: 10.1016/j.acra.2012.03.006
Luft, A. R., Skalej, M., Schulz, J. B., Welte, D., Kolb, R., Bürk, K., et al. (1999). Patterns of age-related shrinkage in cerebellum and brainstem observed in vivo using three-dimensional MRI volumetry. Cereb. Cortex 9, 712–721. doi: 10.1093/cercor/9.7.712
Maguire, E. A., Burgess, N., and O’Keefe, J. (1999). Human spatial navigation: cognitive maps, sexual dimorphism, and neural substrates. Curr. Opin. Neurobiol. 9, 171–177. doi: 10.1016/S0959-4388(99)80023-3
Markham, J. A., and Juraska, J. M. (2002). Aging and sex influence the anatomy of the rat anterior cingulate cortex. Neurobiol. Aging 23, 579–588. doi: 10.1016/S0197-4580(02)00004-0
Matzel, L. D., Grossman, H., Light, K., Townsend, D., and Kolata, S. (2008). Age-related declines in general cognitive abilities of Balb/C mice are associated with disparities in working memory, body weight, and general activity. Learn. Mem. 15, 733–746. doi: 10.1101/lm.954808
Mechelli, A. (2005). Structural covariance in the human cortex. J. Neurosci. 25, 8303–8310. doi: 10.1523/JNEUROSCI.0357-05.2005
Mengler, L., Khmelinskii, A., Diedenhofen, M., Po, C., Staring, M., Lelieveldt, B. P. F., et al. (2014). Brain maturation of the adolescent rat cortex and striatum: changes in volume and myelination. Neuroimage 84, 35–44. doi: 10.1016/j.neuroimage.2013.08.034
Misquitta, K. A., Codeluppi, S. A., Knoch, J. K., Bansal, Y., Tomoda, T., Ellegood, J., et al. (2021). Behavioral and neurostructural changes associated with chronic amygdala hyperactivation. bioRxiv [Preprint]. doi: 10.1101/2021.09.11.459894
Montembeault, M., Joubert, S., Doyon, J., Carrier, J., Gagnon, J.-F., Monchi, O., et al. (2012). The impact of aging on gray matter structural covariance networks. Neuroimage 63, 754–759. doi: 10.1016/j.neuroimage.2012.06.052
Nadig, A., Seidlitz, J., McDermott, C. L., Liu, S., Bethlehem, R., Moore, T. M., et al. (2021). Morphological integration of the human brain across adolescence and adulthood. Proc. Natl. Acad. Sci. U.S.A. 118:e2023860118. doi: 10.1073/pnas.2023860118
Nadon, N. L. (2006). Exploiting the rodent model for studies on the pharmacology of lifespan extension. Aging Cell 5, 9–15. doi: 10.1111/j.1474-9726.2006.00185.x
Nikolova, Y. S., Misquitta, K. A., Rocco, B. R., Prevot, T. D., Knodt, A. R., Ellegood, J., et al. (2018). Shifting priorities: highly conserved behavioral and brain network adaptations to chronic stress across species. Transl. Psychiatry 8:26.
Öngür, D., and Price, J. L. (2000). The organization of networks within the orbital and medial prefrontal cortex of rats, monkeys and humans. Cereb. Cortex 10, 206–219. doi: 10.1093/cercor/10.3.206
Onoda, K., and Yamaguchi, S. (2013). Small-worldness and modularity of the resting-state functional brain network decrease with aging. Neurosci. Lett. 556, 104–108. doi: 10.1016/j.neulet.2013.10.023
Oschwald, J., Guye, S., Liem, F., Rast, P., Willis, S., Röcke, C., et al. (2019). Brain structure and cognitive ability in healthy aging: a review on longitudinal correlated change. Rev. Neurosci. 31, 1–57. doi: 10.1515/revneuro-2018-0096
Pardo, J. V., Lee, J. T., Sheikh, S. A., Surerus-Johnson, C., Shah, H., Munch, K. R., et al. (2007). Where the brain grows old: decline in anterior cingulate and medial prefrontal function with normal aging. Neuroimage 35, 1231–1237. doi: 10.1016/j.neuroimage.2006.12.044
Pezawas, L., Verchinski, B. A., Mattay, V. S., Callicott, J. H., Kolachana, B. S., Straub, R. E., et al. (2004). The brain-derived neurotrophic factor val66met polymorphism and variation in human cortical morphology. J. Neurosci. 24, 10099–10102. doi: 10.1523/JNEUROSCI.2680-04.2004
Prasad, K., Rubin, J., Mitra, A., Lewis, M., Theis, N., Muldoon, B., et al. (2022). Structural covariance networks in schizophrenia: a systematic review Part I. Schizophr. Res. 240, 1–21. doi: 10.1016/j.schres.2021.11.035
Prevot, T. D., Li, G., Vidojevic, A., Misquitta, K. A., Fee, C., Santrac, A., et al. (2019). Novel benzodiazepine-like ligands with various anxiolytic, antidepressant, or pro-cognitive profiles. Complex Psychiatry 5, 84–97. doi: 10.1159/000496086
Prevot, T. D., Sumitomo, A., Tomoda, T., Knutson, D. E., Li, G., Mondal, P., et al. (2021). Reversal of age-related neuronal atrophy by α5-GABAA receptor positive allosteric modulation. Cereb. Cortex 31, 1395–1408. doi: 10.1093/cercor/bhaa310
Qiu, L. R., Fernandes, D. J., Szulc-Lerch, K. U., Dazai, J., Nieman, B. J., Turnbull, D. H., et al. (2018). Mouse MRI shows brain areas relatively larger in males emerge before those larger in females. Nat. Commun. 9:2615. doi: 10.1038/s41467-018-04921-2
Raz, N. (1997). Selective aging of the human cerebral cortex observed in vivo: differential vulnerability of the prefrontal gray matter. Cereb. Cortex 7, 268–282. doi: 10.1093/cercor/7.3.268
Raz, N., and Rodrigue, K. M. (2006). Differential aging of the brain: patterns, cognitive correlates and modifiers. Neurosci. Biobehav. Rev. 30, 730–748. doi: 10.1016/j.neubiorev.2006.07.001
Raz, N., Gunning-Dixon, F., Head, D., Williamson, A., and Acker, J. D. (2001). Age and sex differences in the cerebellum and the ventral pons: a prospective mr study of healthy adults. Am. J. Neuroradiol. 22, 1161–1167.
Resnick, S. M., Pham, D. L., Kraut, M. A., Zonderman, A. B., and Davatzikos, C. (2003). Longitudinal magnetic resonance imaging studies of older adults: a shrinking brain. J. Neurosci. 23, 3295–3301. doi: 10.1523/JNEUROSCI.23-08-03295.2003
Richards, K., Watson, C., Buckley, R. F., Kurniawan, N. D., Yang, Z., Keller, M. D., et al. (2011). Segmentation of the mouse hippocampal formation in magnetic resonance images. Neuroimage 58, 732–740. doi: 10.1016/j.neuroimage.2011.06.025
Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003
Sawiak, S. J., Wood, N. I., Williams, G. B., Morton, A. J., and Carpenter, T. A. (2013). Voxel-based morphometry with templates and validation in a mouse model of Huntington’s disease. Magn. Reson. Imaging 31, 1522–1531. doi: 10.1016/j.mri.2013.06.001
Shizuka, D. (2019). Worked example: centrality and community structure of us air transportation network. Network analysis in R. Available online at: https://dshizuka.github.io/networkanalysis/example_usairports.html (accessed July 14, 2022).
Song, J., Birn, R. M., Boly, M., Meier, T. B., Nair, V. A., Meyerand, M. E., et al. (2014). Age-related reorganizational changes in modularity and functional connectivity of human brain networks. Brain Connect. 4, 662–676. doi: 10.1089/brain.2014.0286
Spring, S., Lerch, J. P., and Henkelman, R. M. (2007). Sexual dimorphism revealed in the structure of the mouse brain using three-dimensional magnetic resonance imaging. Neuroimage 35, 1424–1433. doi: 10.1016/j.neuroimage.2007.02.023
Steadman, P. E., Ellegood, J., Szulc, K. U., Turnbull, D. H., Joyner, A. L., Henkelman, R. M., et al. (2014). Genetic effects on cerebellar structure across mouse models of autism using a magnetic resonance imaging atlas. Autism Res. 7, 124–137. doi: 10.1002/aur.1344
Stern, Y. (2009). Cognitive reserve. Neuropsychologia 47, 2015–2028. doi: 10.1016/j.neuropsychologia.2009.03.004
Straub, S., El-Sanosy, E., Emmerich, J., Sandig, F. L., Ladd, M. E., and Schlemmer, H.-P. (2023). Quantitative magnetic resonance imaging biomarkers for cortical pathology in multiple sclerosis at 7 T. NMR Biomed. 36:e4847. doi: 10.1002/nbm.4847
Sullivan, E. (1995). Age-related decline in MRI volumes of temporal lobe gray matter but not hippocampus. Neurobiol. Aging 16, 591–606. doi: 10.1016/0197-4580(95)00074-O
Sullivan, E. V., Adalsteinsson, E., Sood, R., Mayer, D., Bell, R., McBride, W., et al. (2006). Longitudinal brain magnetic resonance imaging study of the alcohol-preferring rat. Part I: adult brain growth. Alcohol. Clin. Exp. Res. 30, 1234–1247. doi: 10.1111/j.1530-0277.2006.00145.x
Taki, Y., Kinomura, S., Sato, K., Goto, R., Wu, K., Kawashima, R., et al. (2011). Correlation between gray/white matter volume and cognition in healthy elderly people. Brain Cogn. 75, 170–176. doi: 10.1016/j.bandc.2010.11.008
Teixeira, C. M., Pomedli, S. R., Maei, H. R., Kee, N., and Frankland, P. W. (2006). Involvement of the anterior cingulate cortex in the expression of remote spatial memory. J. Neurosci. 26, 7555–7564. doi: 10.1523/JNEUROSCI.1068-06.2006
Tijms, B. M., Seriès, P., Willshaw, D. J., and Lawrie, S. M. (2012). Similarity-based extraction of individual networks from gray matter MRI scans. Cereb. Cortex 22, 1530–1541. doi: 10.1093/cercor/bhr221
Tisserand, D. J., Pruessner, J. C., Sanz Arigita, E. J., van Boxtel, M. P. J., Evans, A. C., Jolles, J., et al. (2002). Regional frontal cortical volumes decrease differentially in aging: an MRI study to compare volumetric approaches and voxel-based morphometry. Neuroimage 17, 657–669. doi: 10.1006/nimg.2002.1173
Tulving, E., and Markowitsch, H. J. (1997). Memory beyond the hippocampus. Curr. Opin. Neurobiol. 7, 209–216. doi: 10.1016/S0959-4388(97)80009-8
Ullmann, J. F. P., Watson, C., Janke, A. L., Kurniawan, N. D., and Reutens, D. C. (2013). A segmentation protocol and MRI atlas of the C57BL/6J mouse neocortex. Neuroimage 78, 196–203. doi: 10.1016/j.neuroimage.2013.04.008
Van Petten, C. (2004). Relationship between hippocampal volume and memory ability in healthy individuals across the lifespan: review and meta-analysis. Neuropsychologia 42, 1394–1413. doi: 10.1016/j.neuropsychologia.2004.04.006
Walhovd, K. B., Fjell, A. M., Reinvang, I., Lundervold, A., Dale, A. M., Eilertsen, D. E., et al. (2005). Effects of age on volumes of cortex, white matter and subcortical structures. Neurobiol. Aging 26, 1261–1270. doi: 10.1016/j.neurobiolaging.2005.05.020
Walhovd, K. B., Westlye, L. T., Amlien, I., Espeseth, T., Reinvang, I., Raz, N., et al. (2011). Consistent neuroanatomical age-related volume differences across multiple samples. Neurobiol. Aging 32, 916–932. doi: 10.1016/j.neurobiolaging.2009.05.013
Westerhaus, M. J., and Loewy, A. D. (2001). Central representation of the sympathetic nervous system in the cerebral cortex. Brain Re. 903, 117–127. doi: 10.1016/S0006-8993(01)02453-2
Wheeler, A. L., Wessa, M., Szeszko, P. R., Foussias, G., Chakravarty, M. M., Lerch, J. P., et al. (2015). Further neuroimaging evidence for the deficit subtype of schizophrenia: a cortical connectomics analysis. JAMA Psychiatry 72, 446–455. doi: 10.1001/jamapsychiatry.2014.3020
Wise, R. A. (2005). Forebrain substrates of reward and motivation. J. Comp. Neurol. 493, 115–121. doi: 10.1002/cne.20689
Yonelinas, A. P. (2013). The hippocampus supports high-resolution binding in the service of perception, working memory and long-term memory. Behav. Brain Res. 254, 34–44. doi: 10.1016/j.bbr.2013.05.030
Yoon, T., Okada, J., Jung, M. W., and Kim, J. J. (2008). Prefrontal cortex and hippocampus subserve different components of working memory in rats. Learn. Mem. 15, 97–105. doi: 10.1101/lm.850808
Zatorre, R. J., Fields, R. D., and Johansen-Berg, H. (2012). Plasticity in gray and white. Nat. Neurosci. 15, 528–536. doi: 10.1038/nn.3045
Zhang, S., Xu, M., Kamigaki, T., Hoang Do, J. P., Chang, W.-C., Jenvay, S., et al. (2014). Long-range and local circuits for top-down modulation of visual cortex processing. Science 345, 660–665. doi: 10.1126/science.1254126
Zheng, F., Liu, Y., Yuan, Z., Gao, X., He, Y., Liu, X., et al. (2019). Age-related changes in cortical and subcortical structures of healthy adult brains: a surface-based morphometry study: age-related study in healthy adult brain structure. J. Magn. Reson. Imaging 49, 152–163. doi: 10.1002/jmri.26037
Zhu, J., Cheng, Q., Chen, Y., Fan, H., Han, Z., Hou, R., et al. (2020). Transient delay-period activity of agranular insular cortex controls working memory maintenance in learning novel tasks. Neuron 105, 934–946.e5. doi: 10.1016/j.neuron.2019.12.008
Keywords: mouse model, aging, brain volume changes, MRI, Structural covariance network (SCN), cognition, working memory, lifespan
Citation: Clifford KP, Miles AE, Prevot TD, Misquitta KA, Ellegood J, Lerch JP, Sibille E, Nikolova YS and Banasr M (2023) Brain structure and working memory adaptations associated with maturation and aging in mice. Front. Aging Neurosci. 15:1195748. doi: 10.3389/fnagi.2023.1195748
Received: 28 March 2023; Accepted: 13 June 2023;
Published: 06 July 2023.
Edited by:
Ian M. McDonough, The University of Alabama, United StatesReviewed by:
Firat Kara, Mayo Clinic, United StatesAbbi R. Hernandez, University of Alabama at Birmingham, United States
Copyright © 2023 Clifford, Miles, Prevot, Misquitta, Ellegood, Lerch, Sibille, Nikolova and Banasr. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mounira Banasr, bW91bmlyYS5iYW5hc3JAY2FtaC5jYQ==; Yuliya S. Nikolova, eXVsaXlhLm5pa29sb3ZhQGNhbWguY2E=
†These authors share first authorship