Impact of social isolation on grey matter structure and cognitive functions: A population-based longitudinal neuroimaging study
Abstract
Background:
Social isolation has been suggested to increase the risk to develop cognitive decline. However, our knowledge on causality and neurobiological underpinnings is still limited.
Methods:
In this preregistered analysis, we tested the impact of social isolation on central features of brain and cognitive ageing using a longitudinal population-based magnetic resonance imaging (MRI) study. We assayed 1992 cognitively healthy participants (50–82years old, 921women) at baseline and 1409 participants after~6y follow-up.
Results:
We found baseline social isolation and change in social isolation to be associated with smaller volumes of the hippocampus and clusters of reduced cortical thickness. Furthermore, poorer cognitive functions (memory, processing speed, executive functions) were linked to greater social isolation, too.
Conclusions:
Combining advanced neuroimaging outcomes with prevalent lifestyle characteristics from a well-characterized population of middle- to older aged adults, we provide evidence that social isolation contributes to human brain atrophy and cognitive decline. Within-subject effects of social isolation were similar to between-subject effects, indicating an opportunity to reduce dementia risk by promoting social networks.
Funding:
European Union, European Regional Development Fund, Free State of Saxony, LIFE-Leipzig Research Center for Civilization Diseases, University of Leipzig, German Research Foundation.
Editor's evaluation
This is a large-scale study on how social isolation influences brain structure and cognition. Baseline social isolation and change in social isolation were associated with smaller hippocampus volumes, reduced cortical thickness, and poorer cognitive function. The level of evidence is strong and converging cross-sectionally and longitudinally.
https://doi.org/10.7554/eLife.83660.sa0Introduction
Over 50 million humans suffer from dementia today. In just 20 years, this number will likely double. Already now, dementia’s global annual costs exceed one trillion US dollars (Prince et al., 2015), and its detrimental effects on the lives of the afflicted make it a major contributor to the world’s burden of disease (Vos et al., 2020).
Research on pharmacological interventions targeting dementia pathogenesis has not yielded any result with a clear clinical benefit yet (Knopman et al., 2021), and available drugs targeting cognitive symptoms offer at most a minor alleviation (Knight et al., 2018). Henceforth, prevention is of cardinal importance and potentially modifiable risk factors are our most promising target (Livingston et al., 2020).
Systematic reviews and meta-analyses have concluded that social isolation, the objective lack of social contact, is such a risk factor for dementia (Kuiper et al., 2015; Penninkilampi et al., 2018) and its main feature cognitive decline (Evans et al., 2019; Kuiper et al., 2016). Assuming causal relationships, Livingston et al. calculated population-attributable fractions for risk factors for dementia and concluded that 3.5% of cases could be attributed to social isolation. This is almost as many as to obesity, hypertension, and diabetes combined (Livingston et al., 2020).
Risk factors of later dementia development often affect the structural brain changes dementia is characterized by: vascular degeneration, amyloid plaques, tau fibrillary tangles, neural degeneration, and grey matter loss. Neuroimaging correlates of these brain changes have been observed multiple years prior to symptom onset in autosomal-dominant dementia (Gordon et al., 2018) and can already be detected in cognitively healthy persons using neuroimaging (Ewers et al., 2011; Jack et al., 2013). Thus, brain magnetic resonance imaging (MRI) can be a potent dementia-risk indicator (Wang et al., 2019), might offer pivotal guidance to identify patients for intensive dementia prevention (Ten Kate et al., 2018), and serve as secondary outcome for intervention trials (Stephen et al., 2019). Still, the link between brain structure and social connection, the umbrella term encompassing social isolation, social support, and loneliness, has not received much attention (Wassenaar et al., 2019). Some studies have linked low social connection to an elevated ‘brain age’ gap estimate (de Lange et al., 2021), changes in microstructural (Molesworth et al., 2015; Spreng et al., 2020; Tian et al., 2014), and volumetric measures in brain regions including the hippocampus and the prefrontal cortex (Blumen and Verghese, 2019; Cotton et al., 2020; Düzel et al., 2019; James et al., 2012; Schurz et al., 2021; Shen et al., 2022; Spreng et al., 2020; Taebi et al., 2020); however, these cross-sectional designs render conclusions about causality difficult. In a longitudinal study using a small sample of 70 participants (37 at follow-up) > 80 years old, microstructural deteriorations and a larger total white matter hyperintensity volume correlated with decreases in predominantly social activities (Köhncke et al., 2016). Furthermore, it suggested that white matter changes mediated the positive association between social activities and perceptual speed (Köhncke et al., 2016). Mortimer et al., 2012 conducted a small randomized controlled trial (RCT) with older adults and found increased total brain volumes and cognitive function in participants after a social interaction intervention compared to a non-intervention control group.
Taken together, the current evidence suggests social isolation to have an adverse effect on brain health. Still, data from longitudinal studies are required to distinguish between from within-participant effects on brain structure and cognitive function and to gain insights into temporal dynamics and causal relationships. Furthermore, to pointedly leverage the power of such datasets for an improved understanding of the effect of social isolation, conceptual clarity regarding the dimensions of social connection is pivotal but still lacking.
Moreover, no solid evidence on the mechanistic underpinnings of the relationship between social isolation and accelerated brain ageing exists. Several mutually non-exclusive, partly overlapping theories are used to explain the beneficial effects of social interaction (Hultsch et al., 1999; Kawachi and Berkman, 2001). Amongst them, the stress-buffering hypothesis puts forward the beneficial effects of social support in strenuous times on mental, cognitive, and immunological health (Kawachi and Berkman, 2001), yet this mediating effect has not been explored regarding brain measures.
Longitudinal population-based neuroimaging studies now offer reliable sample sizes to gain knowledge on effect sizes and disentangle correlation from causation to better understand the impact of social isolation on brain and cognitive ageing. In this pre-registered analysis, we aimed to determine the relationship between social isolation, measured using the Lubben Social Network Scale (LSNS-6, Lubben et al., 2006), and brain structure and cognitive functions, measured using FreeSurfer segmentations on advanced high-resolution MRI at 3 Tesla and neuropsychological testings, in a large well-characterized longitudinal sample of mid- to late-life individuals (n > 1900) from the Health Study of the Leipzig Research Centre for Civilization Diseases (LIFE) (Engel et al., 2023).
To this end, we applied linear mixed effects modelling and structural equation modelling to predict volume of the hippocampus, a focal point of age-related atrophy and Alzheimer’s disease pathology (Rodriguez et al., 2020), by baseline social isolation and change in social isolation over time. Analogously, we modelled memory performance, processing speed, and executive function, as well as whole-brain vertex-wise cortical thickness. Significance was evaluated based on frequentist p-values and Bayes factors, and we adjusted for control variables including age in all models. Details on MRI preprocessing and predefined statistical analyses were preregistered at https://osf.io/8h5v3/.
We hypothesized that both baseline and change in social isolation would correlate with smaller hippocampal volume, cognitive functions (memory, processing speed, executive functions), and cortical thickness. Additionally, we hypothesized interaction effects of baseline social isolation with change in age in the same direction. Moreover, we aimed to test a mediating role of chronic stress as well as hippocampal volume on cognition in these models and explored possible gender differences in stratified analyses.
Results
We included all individuals equal to or over the age of 50 with available neuroimaging of LIFE (Engel et al., 2023) due to the accelerated volume shrinkage starting at about 50 y of age in the hippocampus (Fjell et al., 2013). To avoid reverse causation, we further excluded cognitive impairment or prior brain pathology such as history of stroke, neurodegenerative disease, or brain tumours. In total, we analysed 1335 participants at baseline and 912 participants at follow-up with a mean age of 67 and 73 y, respectively, thereof 51% women and an ~6 y mean change in age at follow-up. For various sensitivity analyses, we reincluded participants that did not meet our preregistered inclusion criteria from the entire sample of 1992 participants at baseline and 1409 at follow-up. The sample displayed a high prevalence of cardiovascular risk factors, with 60% hypertension and <20% diabetes, and 11–13% had no tertiary education (Table 1).
Individuals exhibited LSNS scores ranging across the whole spectrum, with an average score of 16 and 19.7% scoring below the accepted threshold of 12, indicating elevated risk of social isolation, similar to other populations (Lubben et al., 2006). Individual trends in social isolation are depicted in Appendix 1—figure 1. Note that for further analyses, LSNS values were calculated as 30 – LSNS to make larger values indicate greater social isolation and coefficients should thus be interpreted accordingly. Hippocampus volumes derived from T1-weighted high-resolution anatomical MRI scans at 3T (Reuter et al., 2012) showed shrinkage with higher age of about –0.75% per year (Figure 1, left panel), similar to previous estimates (Fjell et al., 2013). To test the effects of social isolation on hippocampal volume, we conducted hierarchical linear mixed effects models adjusting for confounding effects of age, gender, and random effects of the individual in a first model (model 1), and additionally for cardiovascular risk factors in a second model (model 2). We differentiated within- and between-subject effects (van de Pol and Wright, 2009) of social isolation and investigated the interaction effect of baseline LSNS and change in age to test whether participants that are socially more isolated at baseline experienced more pronounced age-related changes. Please see osf.io/8h5v3/and ‘Methods’ for details.
In our sample, social isolation was positively correlated with not living alone, being married, the number of persons living in the participants’ dwelling, being gainfully employed, younger baseline age, and less change in age but no to gender or having a migration background. See Appendix 1—tables 1 and 2 for descriptive statistics and details of the associations. To contextualize the observed link to SES, a comparison of SES category frequencies in LIFE-Adult and a fully representative sample (Lampert et al., 2013b) is provided in Appendix 1—table 3.
Social isolation and hippocampal volume
We found that both stronger baseline social isolation (values for models 1/2: β = −5.6/–5.7 mm3/point on the LSNS (pt), FDR-corrected q-value (q) = 0.0034/0.0078) and increases in social isolation (β = −4.7/–4.5 mm3/pt, q = 0.0035/0.0066) significantly predict smaller hippocampal volumes independent of confounders (Table 2, Figures 1—3). Significance of these findings is further underlined by Bayes factors of 15–19 for baseline social isolation and of 2–3 for change in social isolation. The effect size of one point on the LSNS is equivalent to a 2.5-month difference in baseline age.
Social isolation and cognitive functions
In analogous linear mixed effects models, we tested the effects of social isolation on cognition, measured using domain-specific composite scores based on z-scored results of the trail-making test (TMT A and B) and the CERAD-plus test battery (CERAD – Consortium to Establish a Registry for Alzheimer’s Disease, RRID:SCR_003016) assessed under standardized conditions (Beyer et al., 2017). Overall, stronger baseline social isolation and to a lesser extent increases in social isolation were linked to worse cognitive performance (Table 3, Figure 1). Specifically, stronger social isolation at baseline significantly predicted lower executive functions (β = −0.028/–0.017 SD/pt, q = 9.6e-09/0.0014) and lower processing speed (β = −0.018/–0.017 SD/pt, q = 1.2e-05/3e-04). The link to lower memory (β = −0.014/–0.008 SD/pt, q = 0.0016/0.0914) was strong in model 1 but did not survive FDR-correction when controlling for additional covariates. Increases in social isolation over time significantly predicted lower memory in models 1 and 2 (β = −0.019/–0.0018 SD/pt, q = 0.0034/0.0142) but not processing speed (β = −0.007/–0.008 SD/pt, q = 0.238/0.198) and executive functions (β = –0.003/0.001 SD/pt, q = 0.41/0.69). Very high Bayes factors corroborate and substantiate the evidence for the negative effect of baseline social isolation on cognitive functions. Figures 2 and 3 allow comparisons of these effects with other predictors for the different dependent variables.
We did not observe interaction effects of social isolation on hippocampal volume or cognitive performance with age. Appendix 1—tables 4–6 provide a comprehensive summary of all LMEs and predictors including covariates.
Social isolation and cortical thickness
To explore whether social isolation affects regional cortical thickness, we conducted whole-brain vertex-wise linear mixed effects analyses on FreeSurfer-derived 3D cortical maps (Reuter et al., 2012). In model 1, we found a total of eight clusters of significantly decreased cortical thickness associated with stronger baseline social isolation after FDR-correction with an α-level of 5% (Figure 4). The clusters were located in the left precuneus, cuneus, precentral gyrus and posterior cingulate gyrus, and right supramarginal gyrus and cuneus. Increases in social isolation over time were linked to decreased cortical thickness in one cluster in the right superior frontal gyrus (Figure 5). When additionally controlling for cardiovascular covariates (model 2), no significant clusters were detected. Table 4 lists these clusters, their locations, and sizes.
Mediation analyses
Turning to the stress-buffering hypothesis, we investigated whether perceived stress, measured using the Trierer Inventar zum chronischen Stress (TICS) (Schulz and Schlotz, 1999), mediated the relationship of social isolation and hippocampal volume. Moreover, we investigated whether hippocampal volume mediated the association between social isolation and cognitive functions. Specifically, we investigated the indirect path resulting from the regressions of follow-up mediator on baseline LSNS and follow-up dependent variable on baseline mediator.
Neither the mediation analyses with chronic stress as a mediator (n = 51 complete observations) nor the mediation analyses with hippocampal volume as a mediator (n = 341–360) yielded significant results. Due to the requirements of the model design and over 50% missingness in the stress questionnaire, the sample sizes of the mediation analyses were gravely diminished. Details on the mediation analyses are provided in Appendix 1—table 7.
Sensitivity analyses
In addition to these pre-registered analyses, we conducted sensitivity analyses to test the robustness of our results on hippocampal volume and cognitive functions. These included possible effects of the Covid-19 pandemic, effects related to the definition of exclusion criteria or confounder specificities. Analyses accounting for (a) potential effects of measurements before compared to during the Covid-19 pandemic, (b) reducing the exclusion criteria (i.e. not excluding cognitively impaired participants, participants taking centrally active medication, and participants with recent cancer treatment), (c) only including participants with two timepoints and using mean and within scores, (d) using a hypertension cut-off of 140 mmHg, (e) using an MMSE cut-off of <27, (f) additionally controlling for physical activity, (g) additionally controlling for sleep quality, and (h) standardizing cognitive functions using the baseline mean rather than the grand mean confirmed the regression coefficients of our models in terms of direction and size (Appendix 1—tables 8–15).
Moreover, we found that treating social isolation as a dichotomous variable, using the standard LSNS cut-off of 12 points, led to results very similar to those of our analyses with continuous LSNS scores (Appendix 1—tables 16–18). However, we found no evidence for an interaction effect of continuous and categorical LSNS variables (Appendix 1—tables 19 and 20).
Of note, neuroscience has historically neglected sex and gender differences, predominantly resulting in increased misdiagnoses of and relatively worse treatments for women (Shansky and Murphy, 2021). Therefore, we recalculated analyses in gender-stratified samples (n women = 1110 observations, n male = 1137 observations) to test for differences in the effects of social isolation (Appendix 1—table 21). No clear pattern of difference emerged between women and men. A minor observable difference was that the effect of social isolation on hippocampal volume was mostly driven by baseline social isolation amongst women and by change in social isolation amongst men. This pattern was reversed for other outcomes, though.
In order to further investigate the nature of the correlations, we calculated bivariate latent change score (BLCS) models. In these models we simultaneously tested for an effect of baseline social isolation on change in cognitive functions or hippocampal volume and vice versa (see Appendix 1—figure 2 for a visualization). The BLCS models did not produce solid evidence regarding directionality (Appendix 1—table 22). As in the mediation analyses, the design requirements of the BLCS resulted in smaller sample sizes (n = 362–585 complete observations).
Discussion
In this pre-registered study, we investigated the associations of social isolation with brain structure and cognition in a large cognitively healthy mid- to late-life longitudinal sample. In line with our pre-specified hypotheses, we showed a significant link between stronger baseline social isolation and increases in social isolation over the course of ~6 y and smaller hippocampal volumes. Both predictors had an effect size per point on the LSNS comparable to a 2.5 months difference in baseline age in this age range. Simply put, assuming that if everything else remained stable, the difference between having 1 or 3–4 close and supportive friends is comparable to a 1-year difference in hippocampal ageing. Furthermore, we found significant associations of stronger baseline social isolation with lower executive functions, memory, and processing speed. The link to executive functions was particularly strong with an effect size larger than a 1-year difference in baseline age. For increases in social isolation, confidence intervals were wider but effect sizes, except for executive functions, were similar in magnitude to that of baseline social isolation. In multiple sensitivity analyses, we showed the robustness of these findings. Neither applying less exclusion criteria, only including participants with two timepoints, nor controlling for the impact of the ongoing pandemic changed our results substantially. Moreover, we found clusters of decreased cortical thickness in the cuneus, precuneus, precentral, posterior cingulate, supramarginal, and middle frontal gyrus associated with social isolation cross-sectionally or longitudinally. Mediation analysis in smaller sample sizes testing potential effects of social isolation through lowering adverse effects of stress revealed no significant effects.
Hippocampal volume
Our findings indicate that social isolation contributes to grey matter loss in the hippocampus, a focal point of atrophy in mild cognitive impairment (Devanand et al., 2007), and Alzheimer’s dementia (Fox et al., 1996).
Notably, not only baseline social isolation (a between-subject effect) but also change in social isolation (a within-subject effect) significantly predicted hippocampal volume. Through the employment of statistical LMEs, we were able to distinguish and study effects at these different levels (van de Pol and Wright, 2009) and the design helped us to avoid fallacious inferences from single-level data (Robinson, 1950) to which simple linear regressions would have been susceptible. Specifically for the study of social isolation as a risk factor for dementia, it is crucial to disentangle between- and within-subject effects. Social isolation has both been described as a trait (Noonan et al., 2021), implying it to be an invariant between-subject characteristic and as a potential target for interventions (Hussenoeder and Riedel-Heller, 2018), implying it to be a modifiable within-subject effect. The finding of a significant within-subject effect of change in social isolation therefore offers hope for modifiability as it implies that the observed associations are not (exclusively) the effect of an invariant trait. Thus, our data point towards that reducing social isolation could help to maintain hippocampus integrity in ageing.
However, this assumes a causal effect of social isolation. As associations with social isolation could also have resulted from reverse causation through health selection, that is, that participants with accelerated brain ageing are more likely to become socially isolated, this assumption needs careful consideration. Bayes factors imply the absence of an interaction effect of baseline social isolation with change in age and the bivariate latent change score models did not provide evidence in favour of causality in the hypothesized direction either. However, neither did they provide evidence for reverse causality. This inconclusiveness might result from our reduced follow-up sample size and thus related lower power, especially in the latent score models. For example, data from the English Longitudinal Study of Aging from >6000 older adults measured at up to 6 two-year intervals supports the assumed causality of social isolation with regards to memory performance (Read et al., 2020). In our study, presence of considerable effect sizes and the high statistical confidence in these estimates on multiple outcomes in our healthy sample without cognitive impairment speaks against the competing hypothesis of reverse causality through health selection and in favour of a causal role of social isolation. Furthermore, the lack of any strong increase in effect size when including health-impaired participants or decrease when applying more stringent exclusion criteria for cognitive health corroborate this interpretation. Still, overall these results only add a modicum of corroboration to the case for a causal role of social isolation.
Cognitive functions
Baseline social isolation, and to a lesser extent, change in social isolation, were significantly associated with cognitive performance, that is, executive functions, processing speed, and memory, all of which undergo decline in (pathological) ageing (Blazer et al., 2015). Again, our results thus imply a detrimental role of social isolation on cognitive functions. We could however not observe that social isolation lowered memory performance through reductions in hippocampal volume, a hypothesis raised by considerations of the central role of the hippocampus in memory (Buzsáki and Moser, 2013). Similarly, we could not find evidence that social isolation affected hippocampal volume through higher chronic stress measured with questionnaires, a hypothesis put forward by the stress buffering theory (Kawachi and Berkman, 2001). These latter analyses suffered from small sample sizes and a limited number of timepoints. Nonetheless, the lack of any significant link between chronic stress and social isolation (see Appendix 1—table 2) is hard to align with the stress-buffering hypothesis in spite of the missingness in the TICS.
Cortical thickness
Overall, comparing our brain morphometric results with those of existing cross-sectional studies on social isolation, detected brain regions coincide. A rather small-sampled study did not find a link between social isolation and grey matter volumes (Lin et al., 2020) but James et al., 2012 (occipital lobe), Blumen and Verghese, 2019 (hippocampus, precuneus, medial frontal gyrus) and Shen et al., 2022 (hippocampus, right supramarginal gyrus) found decreased volumes in regions we detected, too.
Several of the cortical regions identified in our study (precuneus) belong to the pattern of exacerbated regional atrophy found in Alzheimer’s disease. Furthermore, we detected regions known for increased cortical thinning in the healthy process of ageing (cuneus) and both in healthy and pathological ageing (supramarginal gyrus) (Bakkour et al., 2013; Pini et al., 2016). This indicates an aggravating role of social isolation in cortical thinning that may contribute to normal and accelerated brain ageing processes. However, the findings of lower cortical thickness must be interpreted cautiously due to the limited consistency between cross-sectional and longitudinal effects and the exploratory approach of whole-brain analyses.
Limitations
A limitation of this study is its uncertain generalizability to the general population because the sample was probably affected by selection and attrition bias common to longitudinal studies (Chatfield et al., 2005). Attrition bias might have mostly affected the mediation and BLCS models that thus offered reduced interpretability, despite the comparatively large neuroimaging cohort. However, the LMEs were mostly unscathed by this problem due to their ability to make use of datapoints of participants with only one full observation. In addition, our population represents a WEIRD sample (i.e. western, educated, industrialized, rich, democratic) which might skew our understanding of how social isolation affects brain health (Laird, 2021). As we found higher SES to be associated with lower LSNS scores, this relatively high SES sample might have led to underestimation of the detrimental effects of social isolation and increases in social isolation in the ageing process. Considering hippocampus segmentations, it has been argued that FreeSurfer systematically overestimates volumes compared to manual volumetry; however, this difference did barely emerge in participants over the age of 50 (Wenger et al., 2014). A further limitation are ceiling effects in the CERAD word list memory task in healthy adults, potentially limiting the sensitivity to detect subtle differences. In addition, time of day during testing might have affected cognitive performance (Schmidt et al., 2007), yet we did not control for this. Almost all cognitive tests were performed between 9 a.m. and 1 p.m., though. Covariance of social isolation with other variables such as hypertension or diabetes could have influenced the results. However, note that all variance inflation factors (VIFs) were acceptable, indicating low reason for concern regarding multicollinearity. Lastly, inferences from our results regarding dementia aetiology must be made with caution as we did not investigate clinically diagnosed dementia patients.
In quantitative studies, despite its importance in shaping the research process and conclusions, for example, in functional MRI analysis (Botvinik-Nezer et al., 2020), researchers’ influence is often disregarded. In the supplementary text, we offer a brief reflexivity section to make relevant influences on this study transparent and to shortly discuss the value of reflexivity for quantitative science.
Implications for public health and future work
This pre-registered large-scale population neuroimaging analysis adds robust support to the view that social isolation is associated with accelerated brain ageing and cognitive decline in non-demented adults in mid- to late-life. Our findings further imply that social contact protects from detrimental processes and thereby preserves brain structure and function. Henceforth, targeting social isolation through tailored strategies might contribute to maintaining brain health into old age.
We showed that the established LSNS cut-off can be employed by clinicians to identify subjects likely to suffer adverse effects due to social isolation. However, the absence of evidence for more pronounced negative effects of less social contact amongst those that are deemed socially isolated by the cut-off renders a public health strategy focused on high-risk individuals questionable.
While we could not observe significant contributions of physical activity or sleep quality measured using questionnaires in a smaller subsample on brain and cognitive outcomes, previous studies suggested that physical activity (Musich et al., 2022) and sleepiness (Holding et al., 2020) interact with social isolation and could protect against negative health effects of social isolation, and should therefore be explored in future studies that incorporated these outcomes more systematically.
While we see evidence converging on social isolation as a causal risk factor for dementia and cognitive decline, future neuroimaging studies should pay particular attention to questions of temporality in their design to clear up remaining uncertainties. Studies with more numerous timepoints will be of importance to this end and will furthermore allow us to model important aspects like slopes for individual participants (van Doorn et al., 2021). Intervention studies will be the gold standard to provide evidence with regards to the causal role and effect size of social isolation. Multidomain interventions for dementia prevention justifiably become the norm (Stephen et al., 2019), so that effects of reduced social isolation must be investigated as a likely contribution to an aggregate effect.
Illuminating the mechanistic underpinnings of the association should be another focus for future research. Studies might prioritize obtaining reliable proxies for the hypothesized mediators. As elevated cortisol levels, in line with the stress-buffering hypothesis, may exert detrimental effects on cognition and contribute to AD pathology (Ouanes and Popp, 2019), using hair cortisol, a reliable measure of chronic stress (Staufenbiel et al., 2013), could be a promising choice to further investigate this proposed mechanism. In light of the lack of evidence for the stress-buffering hypothesis in our data, alternative mechanistic theories should be pursued, too. The main-effect theory postulates that social relationships foster beneficial health behaviours, affective states and neuroendocrine responses, ultimately protecting neuronal tissue (Kawachi and Berkman, 2001). Others point out that socializing is cognitively demanding and requires engagement with complex environments. In the ‘use-it-or-lose-it’ theory, this is crucial for the maintenance of cognitive function (Hultsch et al., 1999). Promising approaches to answer this research question could be interventions specifically targeting one of the hypothesized detrimental processes in isolated individuals and mediation analyses of multi-wave studies with larger sample sizes. Lastly, reverse causality or simultaneity cannot be completely ruled out yet. However, the observed solid correlations in our healthy sample and the lack of an increase in effect sizes when including participants with dementia or low MMSE scores renders this alternative hypothesis to a causal role of social isolation unlikely.
Moreover, studies investigating social isolation due to lockdown measures and its impact on cognitive and brain health will be of great significance.
In light of the relevance of social isolation for cognitive and general health and well-being (National Academies of Sciences, 2020), its pervasiveness in the elderly population of the global north (Livingston et al., 2020) is alarming. Physical distancing measures have caused an unprecedented rise in the attention to the impact of social isolation but social isolation has been a grave problem before Covid-19 and it will remain a central public health concern thereafter. Existing and future research on reasons for and the role of social isolation in health and disease should provide guidance for the urgently needed development and evaluation of tailored strategies against social isolation and its detrimental effects. These should address social isolation both through intervention strategies on the individual but also societal level, leveraging values like solidarity and communality.
Materials and methods
Study design and preregistration
Request a detailed protocolWe followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) and Committee on Best Practices in Data Analysis and Sharing (COBIDAS) on MRI guidelines in our reporting wherever appropriately applicable.
The study’s preregistration can be found at https://osf.io/8h5v3/. Please refer to it for information on the authors’ previous knowledge of the data and a comprehensive overview of our pre-specified hypotheses and models.
Study population
Request a detailed protocolWe used longitudinal data from the ‘Health Study of the Leipzig Research Centre for Civilization Diseases’ (LIFE). The study was approved by the institutional ethics board of the Medical Faculty of the University of Leipzig and conducted according to the Declaration of Helsinki. The LIFE-Adult-Study is a population-based panel study of around 10,000 randomly selected participants from Leipzig, a major city with 550,000 inhabitants in Germany. A subgroup of around 2600 participants underwent MRI testing at baseline. The baseline examination was conducted from August 2011 to November 2014. Follow-up assessments were performed around 6–7 y after the respective first examinations (Engel et al., 2023). Around 1000 participants of the MRI-subsample returned for follow-up testing.
We included all participants over 50 with MRI data that did not fulfil any of the following exclusion criteria:
Anamnestic history of stroke
Any medical condition (i.e. epilepsy, multiple sclerosis, Parkinson’s disease)/chronic medication use that would compromise cognitive testing (i.e. cancer treatment in the past 12 mo or drugs affecting the central nervous system)
Diagnosed dementia or Mini-Mental State Examination (MMSE)-score < 24
A trained radiologist considered the MRI scans unusable due to brain tumours, or acute ischaemic, haemorrhagic or traumatic lesions
If no MMSE data were available, the participants were excluded if their overall performance in cognitive tests negatively deviated from the wave’s mean by 2 standard deviations (SDs) which is a stricter criterion excluding ~2.6% of the sample compared to ~0.8% excluded based on the MMSE. The exclusion criteria were chosen to reduce the potential of reverse causality, that is, dementia symptoms leading to a loss of social connections, as correlations observed in this cognitively intact sample should not stem from dementia symptoms.
MRI data acquisition, processing, and quality control
Request a detailed protocolWe obtained T1-weighted images on a 3 Tesla Siemens Verio MRI scanner (Siemens Healthcare, Erlangen, Germany) with a 3D MPRAGE protocol and the following parameters: inversion time, 900 ms; repetition time, 2300 ms; echo time, 2.98 ms; flip angle, 9°; field of view, 256 × 240 × 176 mm3; voxel size, 1 × 1 × 1 mm3, shimming: tune-up shim, no fat suppression, whole-brain coverage. We processed the scans with FreeSurfer (FreeSurfer, V5.3.0, RRID:SCR_001847) and the standard cross-sectional pipeline recon-all. FreeSurfer automatically measures hippocampal volume, vertex-wise cortical thickness, and intracranial volume. To ensure high within-subject reliability, we employed FreeSurfer’s longitudinal pipeline on all scans, including those of participants without a follow-up scan. Please see Reuter et al., 2012 for details. Moreover, we smoothed the cortical thickness surfaces with a 10 mm kernel to improve reliability and power (Liem et al., 2015). Different Linux kernels and Ubuntu versions constituted the computational infrastructure during the data acquisition and processing.
Visual quality control was based on the recommendations of Klapwijk et al., 2019. After the baseline data were acquired, our team visually controlled all results of the cross-sectional recon-all pipeline. Additionally, we controlled the outputs of the longitudinal stream of all participants with follow-up data and those whose cross-sectional runs required editing. If we detected errors in the processed scans, we manually edited them (N = 262). We excluded participants from analyses using MRI measures if we deemed the processed scans to be unusable (n = 98).
Variable construction
Social isolation
Request a detailed protocolWe used the standard Lubben Social Network Scale (LSNS) –6 (Lubben et al., 2006) to measure the participants’ social isolation. The questionnaire is a suitable tool to measure social isolation (Valtorta et al., 2016) has a high internal consistency (Cronbach’s α = 0.83), a stable factor structure of the family and non-kin subscale (rotated factor loading comparisons = 0.99) and good convergent validity (correlations with caregiver /emotional support availability and group activity all 0.2–0.46 across multiple sites) (Lubben et al., 2006). In order to make larger scores imply more isolation, we subtracted the actual score from the maximum score of 30.
To quantify changes in social isolation, we subtracted the baseline from the follow-up score. For all baseline observation change in LSNS = 0.
In exploratory analyses testing the standard threshold of 12 points, we converted the continuous scores into a dichotomous categorical variable. Change in LSNS scores for these analyses corresponds to positive or negative category shifts.
Grey matter measures
Request a detailed protocolWe used the hippocampal volume derived from FreeSurfer’s segmentation and averaged it over both hemispheres. Furthermore, we adjusted it for intracranial volume according to the following formula:
where β is the unstandardized regression coefficient of hippocampal volume (HCV) on intracranial volume (ICV) from a linear mixed-effects model (LME) (Jack et al., 1998).
For whole-brain analyses we used the FreeSurfer fsaverage template and cortical thickness as a vertex-wise outcome.
Cognitive functions
Request a detailed protocolWe calculated domain-specific composite scores and calculated them as follows (Beyer et al., 2017):
Executive functions consisted of phonemic and semantic fluency, combined with TMT B/A: executive functions = (z_phonemic fluency + z_semantic fluency + z((TMT B – TMT A)/TMT A))/3.
For the memory score, we defined learning as the sum of three consecutive learning trials of the CERAD word list (10 words), recall as the sum of correctly recalled words after a delay, in which participants performed a nonverbal task, and recognition as the number of correctly recognized words out of a list of 20 presented afterwards: memory = (z_learning + z_recall + z_recognition)/3
Processing speed was defined as the negated z-scored TMT part A score.
Sum-score = z_phonemic fluency + z_semantic fluency + z_sum_learning + z_recall + z_recognition + z((TMT B – TMT A)/TMT A)
Most participants were cognitively tested between 9 a.m. and 1 p.m.
Stress
Request a detailed protocolTrierer Inventar zum chronischen Stress (TICS) is a German questionnaire assessing perceived stress (57 items, six sub-scales, 0–4 points per item). Its sum score is our measure of participants’ chronic stress. The subscales have acceptable to excellent internal consistency (Cronbach’s α = 0.76–0.091) and criterion validity of the work overload sub-scale has been shown by demonstrating a significant correlation with cortisol levels over the course of a work days and its ability to differentiate tinnitus patients from healthy controls (Schulz and Schlotz, 1999).
Control and further variables of interest
Request a detailed protocolMonth and year of birth of the participants and the date of the MRIs were recorded and used to calculate the age to one decimal point. Age = YOM.MOM – YOB.MOB (YOM/MOM = year/month of MRI, YOB/MOB = year/month of birth). If no MRI was available, we used the date of the LSNS.
For follow-up observations, we calculated: change in age = age at follow-up - baseline age. For all baseline observations change in age = 0.
Data on the following variables was only available for the baseline. Henceforth, we used the baseline values of these control variables for both timepoints.
We calculated the body-mass-index (BMI) according to the standard formula: BMI = weight [kg]/(height [m])2
In order to control for hypertension and diabetes, we used dichotomized variables. Participants were categorized as hypertensive if they had a previous diagnosis of hypertension, took antihypertensive medication or had an average systolic blood pressure over 160 mmHg. The systolic blood pressure was measured three times. The first measurement was performed after 5 min of rest and three additional minutes of rest passed between each of the following measurements. Participants were categorized as diabetic if they had a previous diagnosis of diabetes, took antidiabetic medication, or HbA1C measured by turbidimetry was ≥ 6%.
The participants’ education was assessed using an extensive questionnaire (Lampert et al., 2013a) and dichotomously categorized based on prior research on education as a protective factor against dementia (Then et al., 2016). Please see the supplementary text for details.
Participants had to choose their gender in a binary female/male question. Note that the German ‘Geschlecht’ does not differentiate between sex and gender. The lack of a clarification and other options is lamented by the authors.
We used the sum-score of the Center for Epidemiological Studies Depression Scale (CES-D) to measure depressive symptoms (Radloff, 1977).
For a sensitivity analysis, we created a dichotomous variable coded as 1 if participants answered the LSNS questionnaire after March 22, 2020 (first SARS-CoV-2 lockdown in Germany).
For further sensitivity analyses, we used the global Pittsburg Sleep Quality Index (PSQI) score calculated based on the method proposed by the PSQI authors to measure quality of sleep (Buysse et al., 1989) and total physical activity MET-minutes/week as a continuous variable calculated using the International Physical Activity Questionnaire (IPAQ) and its guidelines to obtain a measurement of physical activity (Hagströmer et al., 2006).
To explore general participant characteristics of potential relevance to social isolation, we used data on employment, socioeconomic status, marital status, migration background, and number of persons in the participants’ dwelling. We categorized participants as non-working if they declared not to be gainfully employed due to other reasons than studying, military, or alternative service. We only considered participants to be married if they also lived with their spouses to avoid including separated but not yet divorced persons as this is more appropriate for the topic at hand. Beyond marital status, not discriminating between legally married couples, cohabitees, and other forms of joined living, we used the number of fellow persons living in the dwelling as a continuous variable and also constructed a categorical variable distinguishing those participants that live alone from those living with others. Participants were considered to have a migration background if they stated that they or at least one of their parents was not born in Germany, thus approximating the definition of the Federal Statistical Office of Germany. Socioeconomic status was at baseline calculated as a metric variable according to the guidelines developed at the Robert Koch Institute (Lampert et al., 2013b).
To improve the interpretability of our results, we z-transformed the variables BMI, CESD, TICS, executive function, memory performance, and processing speed by the grand mean and centred the variable baseline age. Additionally, we also centred cognitive performance scores by the baseline mean for a sensitivity analysis.
Outliers and Imputation
Request a detailed protocolWe excluded outliers for our core variables based on a cut-off of 3 SDs (LSNS-score, adjusted hippocampal volume, cognitive functions). Please see Figure 6 (flowchart) for the limited effect of outlier exclusion on sample sizes of the different models. For further details on outlier detection and handling regarding covariates, please see the supplementary text.
To avoid an excessive reduction in sample size due to missing data, we performed imputations for missing predictor variables using the sample mean, distributions based on existing data, or the participant’s mean. Please see the supplementary text for information on our procedures of the respective measures.
Furthermore, we used FIML for analyses using structural equation modelling.
Statistical analyses
Request a detailed protocolAll code can be found at https://github.com/LaurenzLammer/socialisolation, (copy archived at Lammer, 2023). Please see the supplementary text for information on the software used for the analyses.
Statistical modelling
Linear mixed effects models
Request a detailed protocolTo investigate the link between social isolation and our outcomes of interest, we employed LMEs with individual as a random effect.
The general structure of the models in the lme4 syntax was:
Dependent variable ~baseline LSNS + change in LSNS + baseline age + change in age + further control variables + (1|participant).
Please see the supplementary text for explicit formulations of all models. We calculated two models for each hypothesis. In model 1, we included age and gender as control variables. Model 2 additionally included education, hypertension, diabetes, depressive symptoms, and BMI. In model 1, the other risk factors are assumed to mediate the effect of social isolation. In model 2, they are assumed to be confounders (see Appendix 1—figure 3 for a visualization). To measure the effect of ageing, we controlled for baseline age and change in age. Analogously, we differentiated within- and between-subject effects (van de Pol and Wright, 2009) of social isolation. Likewise, we calculated the interaction effect of baseline LSNS and change in LSNS. With this methodology we regressed hippocampal volume, the three cognitive functions, and cortical thickness on baseline LSNS, change in LSNS, and the interaction terms. To measure the overall effect of our predictors of interest, we performed a full-null-model comparison (Bolker et al., 2009). In addition to standard p-values, we calculated Bayes factors (BFs). The relative evidence was measured by dividing the BF for the full model by the BF of the null model (Rouder et al., 2016). This allows us to evaluate the evidence in favour of the full hypothesis compared to the null hypothesis and thus also provide evidence for the absence of an effect (Keysers et al., 2020). We report both measures of significance to offer our readers a comprehensive insight into the data, combining the familiarity of classical frequentist inference with the additional implications of BFs (Keysers et al., 2020).
Sensitivity analyses
Request a detailed protocolFor the first analysis we added whether participants were tested after the start of lockdown measures to all LMEs. In the second analysis we did not exclude participants due to the intake of centrally active or cancer medication and cognitive impairment. To probe the reliability of the coefficients for LSNS_change, we ran an analysis excluding all participants with only one timepoint and used standard mean and within score calculation. Furthermore, we ran two sensitivity analysis testing whether using a hypertension cut-off of 140 mmHg or an MMSE cut-off of <27 as an exclusion criterion would affect our results. Additionally, we checked whether results would differ if cognitive test scores were standardized by the baseline rather than grand mean and whether the inclusion of physical activity or sleep quality as an additional control variable would affect the results. Furthermore, to test for potential differences in the effect of social isolation between women and men, we divided our dataset by gender and recalculated the frequentist LMEs with both resulting datasets. Moreover, we investigated whether the standard LSNS cut-off would be a sensitive measurement to indicate adverse effects of social isolation on our outcomes and thus be helpful for clinical practice. To this end, we ran our models treating social isolation as a dichotomous categorical variable. Additionally, we ran them with an interaction term of the usual variables with a social isolation category variable to explore if we would find evidence for stronger adverse effects of less social contact amongst the participants deemed socially isolated by the standard cut-off.
To explore links of social isolation to general participant characteristics, we ran LMEs with random intercepts and LSNS sum score as dependent variable. We calculated separate models with socioeconomic status, living alone, number of persons sharing the participant’s dwelling, age (differentiated into the two variables baseline age and change in age in one model), employment, gender, chronic stress, migration background, and marital status as independent variables in the full dataset.
Statistical inference
Request a detailed protocolWe report one-sided p-values based on the direction of the predictor/path of interest’s regression coefficient and the direction of our pre-defined hypotheses. To obtain one-sided BFs we sampled 10,000 times from the posterior distribution of our predictor of interest’s effect. Then we multiplied the BF by 2 and the percentage of sampled effects in the direction of our pre-defined hypotheses.
Multiplicity control
Request a detailed protocolOur threshold for significance for all tests was p<0.05. To control for multiple hypothesis testing we FDR-corrected families of tests and each individual whole-brain analysis (see the supplementary text for definition of families).
BFs of 3–10 and BFs of 10–30 are commonly considered to be moderate or strong evidence in favour of a hypothesis. To evaluate these thresholds in light of multiplicity, we conducted two simulation studies described in the supplementary text that revealed that using a BF threshold of 10.75 rather than 3 would keep α below 5% and that this would not substantially decrease power.
Model assumptions
Request a detailed protocolTo ensure that our continuous predictors are normally distributed, we plotted their histograms. We had to log-transform the CES-D, IPAQ, and PSQI scores to obtain a normal distribution.
To rule out major collinearity, we calculated VIFs. The VIFs did not surpass the threshold of 10 (Myers, 1990) in any model.
Furthermore, we tested the stability of our LMEs in R by comparing the estimates obtained from the model based on all data with those obtained from models with the levels of the random effects excluded one at a time. This revealed the models to be fairly stable. Moreover, we visually controlled them for heteroskedasticity with both a histogram and a qq-plot. The qq-plots show a heavy-tailed distribution of the residuals in some models. This is only a minor deficit as the models are not intended to make accurate predictions at specific points (Gelman and Hill, 2006).
Fit indices providing further information on the quality of a model fit using structural equation modelling can be found in Appendix 1—tables 23 and 24 (Schermelleh-Engel et al., 2003). Fit index thresholds were surpassed by multiple mediation models. As the BLCS models are saturated, fit indices are uninformative.
Appendix 1
Outliers
We excluded the datapoints (all measures of the timepoint) of all participants with measures deviating from the mean by 3 SD for our core variables (LSNS-score, adjusted hippocampal volume, cognitive functions). In case of TICS-score deviations by 3 SD we replaced the values with ‘NA’ and hence did not include them in mediation analyses.
Considering confounders, highly implausible values (±4 SD) for CES-D-score or BMI were treated as missing datapoints and we replaced them with values imputed according to our imputation plans listed below in order not to overly reduce the sample size.
All outlier analyses were conducted separately for baseline and follow-up measurements.
Imputation
The data on the control variables education, BMI, diabetes, hypertension, age, and gender were complete or mostly complete. Henceforth, we could impute missing datapoints without inducing severe bias by using the sample mean for continuous variables or values drawn from a distribution determined by the existing data for categorical variables.
However, CES-D-scores were an exception amongst our control variables because the questionnaires often missed a single or a few items. As suggested by Bono et al., 2007, we imputed up to four missing items per participant using the person mean. Similarly, we imputed up to one item in the LSNS and up to six items in the TICS using the person mean.
If results from one of the cognitive tests required to calculate a composite score for a cognitive function was missing, we calculated the score based on the average performance in the remainder of available tests contributing to the composite score, if at least two tests were available.
Appendix 1—figure 4 provides an overview of missingness in relevant variables at different LSNS scores.
Families of tests for multiple comparison correction
The LMEs with hippocampal volume and the cognitive functions as dependent variables form one large family except for models regressing on the interaction of baseline LSNS and change in LSNS. In each family, we separately corrected model 1 and model 2 analyses resulting in two families of 12 tests. Additionally, we FDR-corrected each individual whole-brain analysis using the sided two-stage adaptive FDR-correction in the FreeSurfer-toolbox (Bernal-Rusiel et al., 2013b). All other analyses and the whole-brain analyses were considered to be exploratory and must be evaluated as such.
Education
The participants’ education was assessed using an extensive questionnaire and given a score ranging from 1 (no degree at all) to 7 (A-levels + master’s degree [or equivalent] or promotion) according to prior research (Lampert et al., 2013a). The effects of education and the significance of different degrees are likely to be culture specific. Fortunately, a recent study examined the effects of education in a population of elderly residents of the city of Leipzig. In this study education operationalized as having a tertiary degree or not was found to be a significant predictor of dementia incidence (Then et al., 2016). This is approximated with a cut-off at a score <3.6.
Simulation studies
Although it is sometimes claimed that Bayesian statistics do not require any multiplicity control (Gelman et al., 2012), we do not believe that this is the case in our study. A truly Bayesian approach would require researchers to adjust the priors to all other tests with non-independent hypotheses or datasets (Sjölander and Vansteelandt, 2019). This is hardly feasible and hence, in practice, Bayesian statistics are usually employed without taking all dependencies into account and their results are measured against thresholds similar to those of frequentist statistics. Appendix 1—figure 5 shows how this results in an increasing familywise error rate (FWER) with an increasing number of tests in both Bayesian and frequentist statistics using an example from Keysers et al., 2020. De Jong has provided a solution for this problem for ANOVAs that has been implemented in the JASP software (de Jong, 2019) but there is still a great lack of available tools for researchers using other statistical methods. Henceforth, we decided to conduct a simulation study to find a Bayes factor threshold adjustment that should control our FWER similar to α-adjustments in frequentist statistics.
To find the expected number of false positives for a given number of tests and threshold, we replaced the variables for baseline social isolation and change in social isolation with random normally distributed values with the same SD and kept the original dataset otherwise untouched. Then we calculated our 24 LMEs belonging to the families of tests with the modified dataset and repeated this process 42 times. At a BF threshold of 3, 14 of the 1008 tests were false positives and 881 were detected as true negatives. Appendix 1—figure 6 shows a histogram of the resulting Bayes factors. The study suggests that for the family size of 12 tests in our study a threshold of about 10.75 would ensure a FWER below 5%. Appendix 1—table 25 gives an overview of the false positives and FWERs.
Furthermore, we wanted to see how this threshold adjustment would affect the power of our study. For this simulation study we generated a dataset that closely resembles the actual dataset but has different regression coefficients for baseline social isolation and change in social isolation. Instead of the actual coefficients we set the effect size per point on the LSNS to 0.1, 0.2, or 0.5 y of baseline age. We simulated a dataset and calculated a Bayes factor for each model and each effect size. As we only calculated the LMEs without interaction terms for reasons of simplicity, this resulted in a number of 48 Bayes factors from simulated data for each of our 13 runs totalling 624 tests. While our power for the smallest effect sizes was generally small (<10%), it was 85.6% for baseline social isolation with an effect size of half a year of baseline age. Increasing the threshold to 10.75 would not substantially decrease it (81.7%). Appendix 1—tables 26 and 27 provide an overview of the percentages of false negatives and true positives using the thresholds 3 and 10.75.
Deviations from our preregistration
For the most part, we stuck closely to our preregistered plan in this study but departed from it at some points for different reasons.
We used the function q-value instead of p.adjust for the FDR correction for the simple reason that it provides us with a more comprehensive output. As we set the argument pi to 1, q-value is equivalent to the classic procedure (Storey, 2002).
We originally intended to first perform a full-null model comparison using an ANOVA and only follow this up with the function drop1 in case of a significant value for the respective predictor of interest. Our intention was to avoid any multiplicity problems due to testing all predictors. Using the scope argument of drop1 solved the problem more parsimoniously.
Our plan to exclude participants with two or more lesions in their MRI was the result of an internal equivocation regarding the meaning of an abbreviation. We excluded participants based on the type of lesions but not based on lesion count.
Furthermore, we used FIML for analyses using structural equation modelling. The similar results obtained using our preregistered approach can be found in the pre-print (Lammer et al., 2021). Similarly, results based on R version 3.6.1 can be found there. The final analyses were conducted in R version 4.2.2 because 3.6.1 was discontinued at our institute.
Lastly, we changed from the term sex to gender as it seems more appropriate.
Software
We performed most analyses using R (R Project for Statistical Computing, V4.2.2, RRID:SCR_001905). For the whole-brain analyses we used MATLAB (MATLAB, V9.13 (2022b) RRID:SCR_001622).
We used the package lme4 (R package: lme4, RRID:SCR_015654) to calculate LMEs in R. To obtain reliable p-values, we used the Satterthwaite option from the lmerTest package (R package: lmerTest, RRID:SCR_015656; Kuznetsova et al., 2017). In the whole-brain analyses we employed the MATLAB-toolbox provided by FreeSurfer to calculate vertex-wise LMEs (Bernal-Rusiel et al., 2013a). For mediation analyses and BLCS models we used the sem function from the lavaan package (Rosseel, 2012).
We calculated BFs for all LMEs in R using the BayesFactor package and the functions posterior and generalTestBF with default priors (Rouder and Morey, 2012).
FDR-correction was performed using the q-value function (R package: Qvalue, RRID:SCR_001073) in R and the sided two-stage adaptive FDR-correction in the FreeSurfer-toolbox (Bernal-Rusiel et al., 2013b).
VIFs were calculated using the package car (Fox and Weisberg, 2019).
Reflexivity
Reflexivity, a sensitivity to and acknowledgement of the ways in which scientists shape the collected data and research findings, is an established hallmark of scientific rigour in qualitative research (Mays and Pope, 2000; Sandelowski and Barroso, 2002). The challenges addressed by reflexivity are perhaps more pronounced in but by no means exclusive to qualitative studies. Nevertheless (at least in an openly conducted form), it is largely absent from quantitative studies (Ryan and Golden, 2006). Methodological reforms in quantitative research like preregistrations and registered reports (Nosek et al., 2018; Nosek and Lakens, 2014) are valuable tools to limit the researchers’ potential to make data fit their prior assumptions but their scope is limited. They do not address some of the most fundamental issues in epidemiology: Which analogies are used to make sense of the data, which questions are being raised and answered, and which theories are chosen to explain phenomena (Krieger, 2011)? Disclosing personal characteristics, researchers’ values, and positionality relative to the object of research (Berger, 2015) thus helps readers assess a study and its findings more thoroughly. Additionally, an external evaluation of the presence and prevalence of non-empirical decision vectors (Solomon, 2001) in a field of research can be greatly facilitated. Furthermore, as Stephen J. Gould has put it: “It is dangerous for a scholar even to imagine that he might attain complete neutrality, for then one stops being vigilant about personal preferences and their influences – and then one truly falls victim to the dictates of prejudice” (Gould, 1996).
Henceforth, I, as the first author, want to expand this study by a brief reflection on influences that might have played a role in the formation of this study. I am a medical doctoral student with no prior experience in research and conducted this study as the centrepiece of my planned dissertation. Thus, I entered this project with little prior knowledge. I believe that this both made me more flexible and restricted in my choices. On the one hand I was not dedicated to any specific research programme or topic, but on the other hand my reliance on the advice and support from more senior researchers made me emulate their work and methods in many aspects. Further, my worldview has probably made me tend to epidemiological theories (social epidemiology, eco-social theory) (Berkman et al., 2015; Krieger, 2014) broader than the study of lifestyle factors and hence made me choose social isolation as my research topic. A further characteristic that might be of interest to readers is that during the course of the research, two of my relatives struggled with dementia. Ultimately, this reflexivity is inherently limited, as the use of secondary data precludes me from reflecting on the pivotal processes of data acquisition and participant recruitment.
Explicit equations of all LMEs using the lme4 syntax
Variables in bold are dropped in the null model.
H 1.1 Social isolation is negatively associated with hippocampal volume across individuals.
Model111: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +sex + (1|subject)
Model112: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +sex + hypertension
+diabetes + BMI+CESD + education + (1|subject)
H 1.3 Social isolation is negatively associated with hippocampal volume within individuals.
Model131: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +sex + (1|subject)
Model132: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +sex + hypertension
+diabetes + BMI+CESD + education + (1|subject)
H 1.5 Participants that are socially more isolated at baseline will experience aggravated age- related changes in hippocampal volume over the follow-up period.
Model151: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +
LSNS_bl*age_change +sex + (1|subject)
Model152: HCV ~LSNS_bl +LSNS_change +age_bl +age_change +
LSNS_bl*age_change +sex + hypertension +diabetes + BMI+CES.D +
education + (1|subject)
H 2.1 Social isolation is negatively associated with cognitive functions across individuals.
Model211a: executive function ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
(1|subject)
Model212a: executive function ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
hypertension +diabetes + BMI+CES.D+education + (1|subject)
Model211b: memory performance ~LSNS_bl +LSNS_change +age_bl +age_change +
sex + (1|subject)
Model212b: memory performance ~LSNS_bl +LSNS_change +age_bl +age_change +sex
+ hypertension + diabetes +BMI + CES.D+education + (1|subject)
Model211c: processing speed ~LSNS_bl +LSNS_change +age_bl +age_change +
sex + (1|subject)
Model212c: processing speed ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
hypertension +diabetes + BMI+CES.D+education + (1|subject)
H 2.2 Social isolation is negatively associated with cognitive functions within individuals.
Model221a: executive function ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
(1|subject)
Model222a: executive function ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
hypertension +diabetes + BMI+CES.D+education + (1|subject)
Model221b: memory performance ~LSNS_bl +LSNS_change +age_bl +age_change +sex + (1|subject)
Model222b: memory performance ~LSNS_bl +LSNS_change +age_bl +age_change +sex
+hypertension + diabetes +BMI + CES.D+education + (1|subject)
Model221c: processing speed ~LSNS_bl +LSNS_change +age_bl +age_change +
sex + (1|subject)
Model222c: processing speed ~LSNS_bl +LSNS_change +age_bl +age_change +sex +
hypertension +diabetes + BMI+CES.D+education + (1|subject)
H 2.3 Participants that are socially more baseline will experience aggravated age-related changes in cognitive function over the follow-up period.
Model231a: executive function ~LSNS_bl +age_bl +age_change +LSNS_bl*age_change
+sex + (1|subject)
Model231a: executive function ~LSNS_bl +age_bl +age_change +LSNS_bl*age_change
+sex + hypertension +diabetes + BMI+CES.D+education + (1|subject)
Model231b: memory performance ~LSNS_bl +age_bl +age_change +
LSNS_bl*age_change +sex + (1|subject)
Model231b: memory performance ~LSNS_bl +age_bl +age_change +
LSNS_bl*age_change +sex + hypertension +diabetes + BMI+CES.D+education + (1|subject)
Model231c: processing speed ~LSNS_bl +age_bl +age_change +LSNS_bl*age_change
+sex + (1|subject)
Model231c: processing speed ~LSNS_bl +age_bl +age_change +LSNS_bl*age_change
+sex + hypertension +diabetes + BMI+CES.D+education + (1|subject)
H 5.1 In people who are socially more isolated at baseline, an increase in social isolation
from baseline to follow-up will have a stronger negative association with HCV than in people who are less socially isolated at baseline.
Model511: HCV ~LSNS_bl +LSNS_change +LSNS_bl*LSNS_change +age_bl +
age_change +sex + (1|subject)
Model512: HCV ~LSNS_bl +LSNS_change +LSNS_bl*LSNS_change +
age_bl +age_change +sex + hypertenison +diabetes + BMI+CES.D +
education + (1|subject)
Explicit equations of all LMEs using the FreeSurfer LME syntax
H 1.2 Social isolation is negatively associated with vertex-wise cortical thickness across individuals.
For model 1 we built a matrix consisting of six columns: intercept (all ones), age_bl,
age_change, sex, LSNS_bl and LSNS_change.
The corresponding contrast matrix was [0 0 0 0 1 0].
For model 2 we built a matrix consisting of eleven columns: intercept (all ones), age_bl,
age_change, sex, hypertension, diabetes, education, BMI, CES_D, LSNS_bl and
LSNS_change.
The corresponding contrast matrix was [0 0 0 0 0 0 0 0 0 1 0].
H 1.4 Social isolation is negatively associated with vertex-wise cortical thickness within individuals.
For model 1 we built a matrix consisting of six columns: intercept (all ones), age_bl,
age_change, sex, LSNS_bl and LSNS_change.
The corresponding contrast matrix was [0 0 0 0 0 1].
For model 2 we built a matrix consisting of eleven columns: intercept (all ones), age_bl,
age_change, sex, hypertension, diabetes, education, BMI, CES_D, LSNS_bl and
LSNS_change.
The corresponding contrast matrix was [0 0 0 0 0 0 0 0 0 0 1].
H 1.6 Participants that are socially more isolated at baseline, will experience aggravated age-related changes in cortical thickness over the follow-up period.
For model 1 we built a matrix consisting of seven columns: intercept (all ones), age_bl,
age_change, sex, LSNS_bl, LSNS_change and LSNS_bl*age_change. The last term is an
interaction between baseline LSNS and age_change.
The corresponding contrast matrix was [0 0 0 0 0 0 1].
For model 2 we built a matrix consisting of twelve columns: intercept (all ones),
age_bl, age_change, sex, hypertension, diabetes, education, BMI, CES_D, LSNS_bl,
LSNS_change and LSNS_bl*age_change. The last term is an interaction between baseline
LSNS and age_change.
The corresponding contrast matrix was [0 0 0 0 0 0 0 0 0 0 0 1].
Data availability
This study obtained access to the data from LIFE (Leipziger Forschungszentrum für Zivilisationserkrankungen) under project agreement PV-573. All data will be exclusively shared by LIFE (https://www.uniklinikum-leipzig.de/einrichtungen/life) based on individual's project proposal under data protection rules according to the University of Leipzig and can thus not be shared by the authors directly. All code used for the study is available at https://github.com/LaurenzLammer/socialisolation (copy archived at Lammer, 2023).
References
-
BookSocial EpidemiologyOxford University Press.https://doi.org/10.1093/med/9780195377903.001.0001
-
BookCognitive Aging: Progress in Understanding and Opportunities for ActionNational Academies Press.https://doi.org/10.17226/21693
-
Generalized linear mixed models: a practical guide for Ecology and evolutionTrends in Ecology & Evolution 24:127–135.https://doi.org/10.1016/j.tree.2008.10.008
-
Missing data on the center for epidemiologic studies depression scale: a comparison of 4 imputation techniquesResearch in Social & Administrative Pharmacy 3:1–27.https://doi.org/10.1016/j.sapharm.2006.04.001
-
Memory, navigation and Theta rhythm in the hippocampal-entorhinal systemNature Neuroscience 16:130–138.https://doi.org/10.1038/nn.3304
-
Gray Mattervolume covariance networks, social support, and cognition in older adultsThe Journals of Gerontology. Series B, Psychological Sciences and Social Sciences 75:1219–1229.https://doi.org/10.1093/geronb/gbz023
-
Structural brain correlates of loneliness among older adultsScientific Reports 9:13569.https://doi.org/10.1038/s41598-019-49888-2
-
Cohort profile: the LIFE-adult-studyInternational Journal of Epidemiology 52:e66–e79.https://doi.org/10.1093/ije/dyac114
-
Social isolation and cognitive function in later life: A systematic review and meta-analysisJournal of Alzheimer’s Disease 70:S119–S144.https://doi.org/10.3233/JAD-180501
-
Neuroimaging markers for the prediction and early diagnosis of Alzheimer's disease dementiaTrends in Neurosciences 34:430–442.https://doi.org/10.1016/j.tins.2011.05.005
-
Critical ages in the life course of the adult brain: Nonlinear subcortical agingNeurobiology of Aging 34:2239–2247.https://doi.org/10.1016/j.neurobiolaging.2013.04.006
-
BookAnalytical Methods for Social ResearchCambridge University Press.https://doi.org/10.1017/CBO9780511790942
-
Why we (usually) don't have to worry about multiple comparisonsJournal of Research on Educational Effectiveness 5:189–211.https://doi.org/10.1080/19345747.2011.618213
-
Use it or lose it: engaged Lifestyle as a buffer of cognitive decline in aging?Psychology and Aging 14:245–263.https://doi.org/10.1037//0882-7974.14.2.245
-
Primary prevention of dementia: from Modifiable risk factors to a public brain health agendaSocial Psychiatry and Psychiatric Epidemiology 53:1289–1301.https://doi.org/10.1007/s00127-018-1598-7
-
Association of social engagement with brain volumes assessed by structural MRIJournal of Aging Research 2012:512714.https://doi.org/10.1155/2012/512714
-
Social ties and mental healthJournal of Urban Health 78:458–467.https://doi.org/10.1093/jurban/78.3.458
-
A systematic review and meta-analysis of the effectiveness of acetylcholinesterase inhibitors and Memantine in treating the cognitive symptoms of dementiaDementia and Geriatric Cognitive Disorders 45:131–151.https://doi.org/10.1159/000486546
-
BookEpidemiology and the People’s Health: Theory and ContextOxford University Press.https://doi.org/10.1093/acprof:oso/9780195383874.001.0001
-
Social relationships and cognitive decline: a systematic review and meta-analysis of longitudinal cohort studiesInternational Journal of Epidemiology 45:1169–1206.https://doi.org/10.1093/ije/dyw089
-
lmerTest package: tests in linear mixed effects modelsJournal of Statistical Software 82:1–26.https://doi.org/10.18637/jss.v082.i13
-
SoftwareSocialisolation, version swh:1:rev:96e0b8ba6d3c76a53d84f7a813399f33e5dc13d9Software Heritage.
-
Messung des Sozioökonomischen status in der Studie zur Gesundheit Erwachsener in Deutschland (Degs1)Bundesgesundheitsblatt - Gesundheitsforschung - Gesundheitsschutz 56:631–636.https://doi.org/10.1007/s00103-012-1663-4
-
Messung des Sozioökonomischen status in der Studie „Gesundheit in Deutschland Aktuell" (GEDA)Bundesgesundheitsblatt - Gesundheitsforschung - Gesundheitsschutz 56:131–143.https://doi.org/10.1007/s00103-012-1583-3
-
No strong evidence that social network index is associated with gray matter volume from a data-driven investigationCortex; a Journal Devoted to the Study of the Nervous System and Behavior 125:307–317.https://doi.org/10.1016/j.cortex.2020.01.021
-
Social network diversity and white matter Microstructural integrity in humansSocial Cognitive and Affective Neuroscience 10:1169–1176.https://doi.org/10.1093/scan/nsv001
-
Social Isolation and Loneliness in Older Adults: Opportunities for the Health Care SystemSocial isolation and loneliness in older adults, Social Isolation and Loneliness in Older Adults: Opportunities for the Health Care System, Washington, D.C, The National Academies Press, 10.17226/25663.
-
High Cortisol and the risk of dementia and Alzheimer's disease: A review of the literatureFrontiers in Aging Neuroscience 11:43.https://doi.org/10.3389/fnagi.2019.00043
-
The association between social engagement, loneliness, and risk of dementia: A systematic review and meta-analysisJournal of Alzheimer’s Disease 66:1619–1633.https://doi.org/10.3233/JAD-180439
-
Brain atrophy in Alzheimer's disease and agingAgeing Research Reviews 30:25–48.https://doi.org/10.1016/j.arr.2016.01.002
-
ReportWorld Alzheimer Report 2015. The Global Impact of Dementia. An Analysis of Prevalence, Incidence, Cost and TrendsAlzheimer’s Disease International.
-
The CES-D scale: a self-report depression scale for research in the general populationApplied Psychological Measurement 1:385–401.https://doi.org/10.1177/014662167700100306
-
Social isolation and memory decline in later-lifeThe Journals of Gerontology. Series B, Psychological Sciences and Social Sciences 75:367–376.https://doi.org/10.1093/geronb/gbz152
-
Ecological correlations and the behavior of individualsAmerican Sociological Review 15:351.https://doi.org/10.2307/2087176
-
Do high mental demands at work protect cognitive health in old age via hippocampal volume? results from a community sampleFrontiers in Aging Neuroscience 12:622321.https://doi.org/10.3389/fnagi.2020.622321
-
An R package for structural equation modelingJournal of Statistical Software 48:1–36.https://doi.org/10.18637/jss.v048.i02
-
Default Bayes factors for model selection in regressionMultivariate Behavioral Research 47:877–903.https://doi.org/10.1080/00273171.2012.734737
-
Model comparison in ANOVAPsychonomic Bulletin & Review 23:1779–1786.https://doi.org/10.3758/s13423-016-1026-5
-
Finding the findings in qualitative studiesJournal of Nursing Scholarship 34:213–219.https://doi.org/10.1111/j.1547-5069.2002.00213.x
-
Evaluating the fit of structural equation models: tests of significance and descriptive goodness-of-fit measuresMethods of Psychological Research Online 8:23–74.
-
A time to think: circadian rhythms in human cognitionCognitive Neuropsychology 24:755–789.https://doi.org/10.1080/02643290701754158
-
Variability in brain structure and function reflects lack of peer supportCerebral Cortex 31:4612–4627.https://doi.org/10.1093/cercor/bhab109
-
Frequentist versus Bayesian approaches to multiple testingEuropean Journal of Epidemiology 34:809–821.https://doi.org/10.1007/s10654-019-00517-2
-
The default network of the human brain is associated with perceived social isolationNature Communications 11:6393.https://doi.org/10.1038/s41467-020-20039-w