Abstract
Close friendships are important for mental health and cognition in late childhood. However, whether the more close friends the better, and the underlying neurobiological mechanisms are unknown. Using the Adolescent Brain Cognitive Developmental study, we identified nonlinear associations between the number of close friends, mental health, cognition, and brain structure. Although few close friends were associated with poor mental health, low cognitive functions and small areas of the social brain (e.g., the orbitofrontal cortex, the anterior cingulate cortex, the anterior insula and the temporo-parietal junction), increasing the number of close friends beyond a level (around 5) was no longer associated with better mental health and larger cortical areas, and was even related to lower cognition. In children having no more than 5 close friends, the cortical areas related to the number of close friends revealed correlations with the density of μ-opioid receptors and the expression of OPRM1 and OPRK1, and could partly mediate the association between the number of close friends, ADHD symptoms, and crystalized intelligence. Longitudinal analyses showed that both too few and too many close friends at baseline were associated with more ADHD symptoms and lower crystalized intelligence 2 years later. Additionally, we found friendship network size was nonlinearly associated with well-being and academic performance in an independent social network dataset of middle school students. These findings challenge the traditional idea of “the more, the better”, and provide insights into potential brain and molecular mechanisms.
Introduction
Late childhood and its transition towards adolescence is a period marked by decreasing parental influence alongside increasing peer influence. It is a period critical for social interaction, during which friendships are especially important1. During this period, the social brain is still undergoing significant development, in parallel with changes in social cognition2. Evidence suggests that psychiatric disorders often have an onset in adolescence3, which may be partly influenced by the concurrent changes in the social environment and brain4. Therefore, understanding the relationship between friendship, mental health, and cognition during this period, and the underlying brain mechanisms, is of considerable clinical and public health importance.
It has been well-established that positive social relationships such as close friendships are essential for mental health and cognition in children and adolescents5–7. However, it remains unclear whether having more close friends is necessarily better, and if not, whether there is an optimal quantity of close friendships. Cognitive constraints and time resources limit the number of close social ties that an individual can maintain simultaneously8. The innermost layer of the friendship group with the highest emotional closeness is around 5 close friends (the so-called Dunbar’s number)9. For now, only a few empirical studies have examined the nonlinear association between social relationships, mental health, and cognition. A national investigation reported that adolescents with either too many or too few friends had higher levels of depressive symptoms10. Two large-scale studies found that the benefits of social interactions for well-being were nearly negligible once the quantity reached a moderate level11,12. A significant U-shaped effect was detected between positive relations with others and cognitive performance13. Overall, the assumption of linearity still dominates studies of social relationships, and the effect of the friendship network size at the high end remains largely unexplored.
The social brain hypothesis suggests that complex social selection pressures drive the evolution of brain size: for example, relative neocortex volume correlates with social group size in primates14. Both experimental manipulation and free-ranging macaque studies found that social network size could predict mid-superior temporal sulcus volume15,16, a region in which neurons respond to socially relevant stimuli such as face expression and head movement to make or break social contact17. In humans, a series of studies demonstrated that the medial prefrontal cortex (mPFC, i.e. orbitofrontal [OFC] and anterior cingulate cortex [ACC]), the cortex in the superior temporal sulcus (STS), the temporo-parietal junction (TPJ), amygdala and the anterior insula were involved in social cognitive processes18,19. The μ-opioid receptor that is widely distributed in the brain especially in regions implicated in social pain such as the ACC and anterior insula20, has been recognized to be crucial to the formation and maintenance of friendships8. Variation in the μ-opioid receptor gene (OPRM1) is associated with individual differences in rejection sensitivity21. A range of other neurotransmitters including dopamine, serotonin, GABA and noradrenaline may interact with the opioids, and is involved in social affiliation and social behaviour22. Dysregulation of the social brain and neurotransmitter systems has also been proposed in the pathophysiology of major psychiatric disorders23. While it has been proposed that the association between social relationships and mental health is mediated through changes in the social brain24, there is a lack of empirical evidence in late childhood and adolescence.
In this study, based on the existing literature, we hypothesized a nonlinear association between the number of close friends, mental health, and cognitive outcomes. We used data from the Adolescent Brain Cognitive Developmental (ABCD) study25, and another independent social network dataset26 (N > 23,000 in total; Figure 1a). Two different analytic approaches were used to evaluate the nonlinear relation between friendship quantity (predictor) and mental health and cognition (outcome). A significant quadratic term would indicate the presence of nonlinearity, and then a two-lines test27 was conducted to estimate an interrupted regression and to identify the breakpoint (Figure 1b). To explore the underlying neurobiological mechanisms, we further tested the nonlinear association between the number of close friends and brain structure, and then correlated the related brain differences with the density of 8 neurotransmitter systems, and the expression of OPRM1 and κ-opioid receptor (OPRK1) (Figure 1c). Finally, longitudinal and mediation analyses were conducted to uncover the direction and direct association between the number of close friends, mental health, cognition, and brain structure (Figure 1d). We hypothesized that the number of close friends was nonlinearly related to mental health, cognition, and the social brain; and that the nonlinear association between the number of close friends and behavioral measures may be mediated by brain and molecular mechanisms.
Results
Demographic characteristics
In the ABCD study, 7,512 participants (3,625 [48.3%] female, aged 9.91 ± 0.62 years) provided self-reported number of close friends, a broad range of mental health and cognitive measures, and quality-controlled MRI data at baseline, and 4,290 of them (2,044 [47.7%] female, aged 11.49 ± 0.66 years) had two-year follow-up data available. Detailed population characteristics are shown in Table 1. In the social network dataset, 16,065 subjects from 48 middle schools (8,065 [50.3%] female, aged 12.00 ± 1.03 years) who had complete key variables were included (Table S1).
Nonlinear association between the number of close friends, mental health and cognition
The number of close friends was significantly associated with 12 out of 20 mental health measures, and 7 out of 10 cognitive scores at baseline (the total F-value of the linear and quadratic terms, p < 0.05/30; Figure 2a-2g). For these 19 outcomes except the withdrawn/depressed, all quadratic terms reached significance after Bonferroni corrections (p < 0.05/60), and all quadratic models provided a significantly better fit than the corresponding linear models (F = [13.25, 55.53], all p < 0.001). For mental health, the greatest effect sizes of the quadratic terms were observed for social problems (β = 0.06, t = 5.92, p = 3.3⨉10−9, ΔR2 = 0.43%) and attention problems (β = 0.08, t = 5.83, p = 5.8⨉10−9, ΔR2 = 0.42%). For cognition, the greatest effect sizes of the quadratic terms were observed for total intelligence (β = -0.06, t = -7.45, p = 1.0⨉10−13, ΔR2 = 0.50%) and crystalized intelligence (β = -0.05, t = -6.87, p = 6.7⨉10−12, ΔR2 = 0.43%) (Table S2 and Figure S1). The findings were robust with respect to random choice of the siblings (Figure S2).
The average breakpoint of the number of close friends for the mental health and cognitive outcomes with significant quadratic terms was 4.89 ± 0.66 (Figure 2h). The ideal number of close friends was 5 and the closer to that number, the better for participants’ mental health. Similarly, for cognition the optimal number of close friends was 5 and the further from that number the poorer the participants’ cognition. These nonlinear associations were consistent in males and females (Figure S3). However, the number of same-sex close friends, but not of opposite-sex close friends, was significantly related to mental health and cognition (27 out of 30 measures with a significant F-value after Bonferroni correction), and children with 4.00 ± 0.60 same-sex close friends had the best mental health and cognitive functions (Figure S4). Finally, same analyses were performed using the cross-sectional data collected at 2 years later, and the nonlinear associations of the number of close friends with ADHD symptoms and crystalized intelligence remained significant, with an average breakpoint of 4.83 ± 0.75 close friends (Figure S5).
The number of close friends was quadratically associated with brain structure
In the ABCD study, the number of close friends was significantly associated with the total cortical area (F = 6.29, p = 1.0×10−3; Figure 3d), and the total cortical volume (F = 5.80, p = 3.1×10−3). No significant relationship between the number of close friends and mean cortical thickness (F = 0.62, p = 0.54), and total subcortical volume (F = 3.94, p = 0.02) was found.
After FDR correction (q < 0.05), the significant cortical areas associated with the number of close friends were mainly located in the OFC, insula, the ACC, the anterior temporal cortex, and the TPJ (Figure 3a and Table S3). The brain region with the largest effect size for the linear term was the OFC (left medial OFC [area 11l and 13l] and lateral OFC [area 47m and 47s]). The quadratic terms of the number of close friends for all these regions were significant (Figure 3b), and the greatest effect sizes were observed in the temporal pole (left STGa: β = -0.001, t = -4.20, p = 2.8⨉10−5, ΔR2 = 0.18%, Figure 3e; right TGd: β = -0.007, t = -4.11, p = 4.0⨉10−5, ΔR2 = 0.18%, Figure 3f). These findings were robust for random choice of the siblings (Figure S6). Similar findings were found for cortical volumes (Table S4 and Figure S7). As the correlation of cortical area and cortical volume with the number of close friends is high (r = 0.78, p = 3.3×10−76) and cortical area and volume themselves are highly correlated (r = 0.92, p = 9.0×10−151; Figure S8), we focused on cortical area in the following analyses.
Further, two-lines tests suggested that participants with around 5 close friends (breakpoint = 5.30 ± 0.85) had the largest areas in these cortical regions (Figure S9). To illustrate the patterns of nonlinear relationships, we performed linear regression models in participants with ≤ 5 and > 5 close friends, respectively. Similar regions to those found with quadratic models including the OFC, insula, the ACC, and temporal cortex were significant after FDR correction in the ≤ 5 group (Figure S10a and S10b), and the largest effect size was observed in the OFC (Figure S10c). However, the number of close friends was not related to cortical area in the > 5 group (Figure S10d).
The differences in cortical area related to the number of close friends in the two groups were not correlated (r = -0.02, p = 0.78; Figure S10e).
Relationship to molecular architecture
As the number of close friends was nonlinearly associated with cortical area and the significant regions were only found in participants with no more than 5 close friends, we focused on the brain difference pattern for the number of close friends in the ≤ 5 group. We found that the correlations between the spatial pattern of cortical area related to the number of close friends and densities of neurotransmitters were not significant except for the μ-opioid receptor (Spearman’s rho = 0.44, Bonferroni corrected pperm = 0.02; Figure 4a and 4b). Transcriptomic analyses showed that OPRM1 (Spearman’s rho = 0.45, pperm = 0.001; Figure 4c) and OPRK1 (Spearman’s rho = 0.46, pperm = 0.002; Figure 4d) were highly expressed in regions related to the number of close friends.
Longitudinal and mediation results
As the nonlinear association between the number of close friends and ADHD symptoms is relatively strong and robust, and for cognitive outcomes, only crystalized intelligence was collected at 2-year follow-up in the ABCD study, we focused on these two measures in longitudinal and mediation analyses. The cross-lagged panel model (CLPM) revealed that participants having closer to 5 close friends had fewer ADHD symptoms 2 years later (β = 0.04, p < 0.001; Figure 5a). CLPMs in separate groups confirmed that more close friends contributed to fewer ADHD symptoms in ≤ 5 group (β = -0.04, p = 0.003; Figure S11a), but the effect reversed in the > 5 group (β = 0.05, p = 0.019; Figure S11b). The relationship between the absolute difference of close friend number to 5 and crystalized intelligence was bidirectional (Figure 5b). Only in the ≤ 5 group was a significant negative correlation found between crystalized intelligence at baseline and the number of close friends at 2-year follow-up (β = -0.06, p = 0.001; Figure S11c and S11d).
Mediation analyses were used to determine whether and the extent to which the association between the number of close friends, ADHD symptoms, and crystalized intelligence could be explained by the identified cortical areas in the ≤ 5 group. The total identified cortical area partly mediated the association between the number of close friends and ADHD symptoms (6.5%, 95%CI [3.4%, 14%]; path a*b: -0.005, 95%CI [-0.008, -0.003]; Figure 5c), and the mediation effects of individual significant regions ranged from 1.52% to 4.68% (Figure 5d). Similarly, the association between the number of close friends and crystalized intelligence was partly mediated by the total identified cortical area (13.5%, 95%CI [8.1%, 27%]; path a*b: 0.008, 95%CI [0.005, 0.01]; Figure 5e), ranging from 2.58% to 8.52% for each significant region (Figure 5f).
Findings in an independent social network dataset
Utilizing the social network dataset allowed us to validate and extend findings in the ABCD study, as it is an independent and large dataset, the friendship network was generated by nomination, and different measures of mental health and cognition were collected (i.e., well-being and grade point average [GPA]). Three indicators of friendship network size (i.e., outdegree, indegree, and reciprocal degree; Figure S12) were significantly related to well-being (indegree: F = 38.63, p = 1.8×10−17; outdegree: F = 33.55, p = 2.9×10−15; reciprocal degree: F = 53.87, p = 4.8×10−24; Figure S13a) and GPA (indegree: F = 28.08, p = 6.7 × 10−13; outdegree: F = 46.66, p = 6.2 × 10−21; reciprocal degree: F = 192.65, p = 2.1×10−83; Figure S13b). Specifically, for well-being, all linear terms were significant, but only the quadratic term of outdegree was significant after Bonferroni correction (β = -2.9⨉10−4, t = -3.67, p = 2.4⨉10−4, ΔR2 = 0.07%; Figure S13c). For GPA, the quadratic terms of all three indicators were significant, and the greatest effect size was observed in the outdegree (β = -0.001, t = - 6.02, p = 1.8⨉10−9, ΔR2 = 0.17%; Figure S13c). The two-lines tests revealed that the positive association of outdegree with well-being and GPA diminished once the outward nomination reached 7 or 8 (Figure S13d). The results confirmed that friendship network size especially outdegree was nonlinearly related to mental health and cognitive outcomes.
Discussion
The present study showed that close friendship quantity was associated with better mental health and higher cognitive functions in late childhood, however, the beneficial effects diminished or reversed when increasing the number of close friends beyond a moderate level. Consistent with the hypothesis, a quadratic association between the number of close friends and the areas of social brain regions was found, including the orbitofrontal cortex (OFC), the anterior cingulate cortex (ACC), insula, anterior temporal cortex and temporo-parietal junction (TPJ), and these regions mediated the nonlinear association with behavior. In addition, the brain differences related to the number of close friends were correlated with measures of the endogenous opioid involvement of the brain regions.
Applying two different analytic approaches in two independent large-scale datasets, we provide compelling evidence that the number of close friends was nonlinearly associated with various mental health and cognitive outcomes. Too large a social network size or too frequent social contacts were not beneficial for well-being and were even detrimental for physical and mental health10–12,28. One explanation is that the quantity of relationships an individual can maintain at the same time is limited by a combination of cognitive capacity and time8. There is a trade-off between the quantity and quality of friendships, suggesting that an increasing number of close friends may be associated with less intimacy. Meanwhile, as time resources are finite, spending too much time on social activities may lead to insufficient time for study and thereby to lower academic performance. Indeed, adolescents are susceptible to peer influence29. The presence of a peer may increase risk-taking behaviors which can be detrimental to mental health30, and reduce cognitive performance31. Having more close friends may increase the possibility of this kind of influence. In addition, ADHD symptoms or novelty seeking may lead to searching for many superficial friendships. Further, in the ABCD study, we identified an optimal number of 5 close friends, which is consistent with the Dunbar’s number9. People devote about 40% of their total social efforts (e.g., time and emotional capital) to just their 5 most important people32. In a phone-call dataset of almost 35 million users and 6 billion calls, a layered structure was found with the innermost layer of an average of 4.1 people33. Our study extends previous findings that there may be an optimal size of close friendship network in late childhood in terms of its relationship to mental health, cognition, and brain structure. The number of 5 close relationships appears to be consistent with the social network size limited by cognitive and time resources found in anthropological studies, and greatly extends previous research by relating this to brain systems and transmitter systems in the brain.
Although friendships are recognized to be especially significant for psychosocial development in late childhood and adolescence, we know relatively little about how brain development relates to close friendships34. Consistent with the behavioral findings, we found that the number of close friends was nonlinearly related to the cortical areas of social brain regions. Children with more close friends had larger areas of social brain regions, but the positive relationship only held for the group with no more than 5 close friends. The areas of social brain regions partly mediated the relationship of the number of close friends with ADHD symptoms and crystalized intelligence. There are two major systems in the brain related to social behavior: an affective system of the OFC, ACC and the anterior insula, and a mentalizing system including the TPJ35,36. The dorsal ACC and anterior insula play an important role in social pain (i.e., painful feelings associated with social disconnection)37. The OFC receives information about socially relevant stimuli such as face expression and gesture from the cortex in the superior temporal sulcus17,38, and is involved in social behavior by representing social stimuli in terms of their reward value35,39,40. The OFC can then influence behavior by a number of routes including the ACC35,41–43. The OFC volume is associated with social network size, which is partly mediated by mentalizing competence44. Previously published meta-analysis studies reported an overlap in brain activation between all mentalizing tasks in the mPFC and posterior TPJ45. Stimulation of the TPJ improves self-other representations46. Animal studies confirmed a casual effect of social relationships on brain development. Adolescent rodents with deprivation of peer contacts showed brain level changes including reduced synaptic pruning in the prefrontal cortex47.
It is known that the endogenous opioid system has an important role in social affiliative processes22. In vivo human positron emission tomography (PET) studies found that μ-opioid receptor regulation in brain regions such as the amygdala, anterior insula, and the ACC may preserve and promote emotional well-being in the social environment48. Variation in the μ-opioid receptor gene (OPRM1) was associated with individual differences in rejection sensitivity, which was mediated by dorsal ACC activity in social rejection21. OPRM1 variation was also related to social hedonic capacity49. Pain tolerance, which is associated with activation of the μ-opioid receptor, was correlated with social network size in humans50. Social behaviors like social laughter and social touch increased pleasurable sensations and triggered endogenous opioid release to maintain social relationships51–53. The opioid system is associated with major psychiatric disorders especially depression54.
Several issues should be taken into account when considering our findings. First, as an association study, no causal conclusion should be made in this study. Whether the number of close friends drives the social brain development, or children with larger social brains have more close friends, is unclear. A bidirectional relationship has been reported in the literature14. Second, the egocentric network size used in the ABCD study is not able to differentiate reciprocal and unilateral friendships. It is possible that the findings may differ using different indicators55,56. For example, in the social network dataset, we showed that outdegree but not indegree and reciprocal degree was quadratically associated with the GPA. Third, the quality of close friendships has not been considered in this study. It has been reported that the relationship between having more friends and fewer depressive symptoms in adolescence was mediated by a sense of belonging57. Although current findings on the relative importance of friendship quantity and quality are inconsistent58,59, it is necessary to include their interactions in future studies. Finally, although we found children with 5 close friends had the best mental health and cognitive functions, and the largest social brain areas, the optimal size cannot be generalized to populations of other ages and cultures. Extension of the findings in other cultural and age groups would be useful.
In conclusion, this study provides new evidence going beyond previous research that a larger number of close friends up to a moderate level in late childhood is associated with better mental health and higher cognitive functions, and that this can be partly explained by the size of the social brain including the OFC and ACC, and the endogenous opioid system. This study may have implications for targeted friendship interventions in the transition from late childhood to early adolescence.
Materials and Methods
Participants and behavioral measures
The ABCD study
The Adolescent Brain Cognitive Development (ABCD) Study is tracking the brain development and health of a nationally-representative sample of children aged 9 to 11 years from 21 centers throughout the United States (https://abcdstudy.org). Parents’ full written informed consent and all children’s assent were obtained by each center, and research procedures and ethical guidelines were followed in accordance with the Institutional Review Boards60. The current study was conducted on the ABCD Data Release 4.0. At baseline, 8,835 individuals from 7,512 families (6,225 [82.9%] with a child, 1,252 [16.7%] with 2 children, 34 [0.5%] with 3 children, and 1 [0.01%] with 4 children) had complete behavioral and structural MRI data. To avoid the influence of family relatedness, we randomly picked only one child in each family, finally resulting in 7,512 children, of whom 4,290 had two-year follow-up data.
Close friendships are characterized by enjoying spending time together, having fun, and trust. Participants were asked how many close friends that are boys and girls they have, respectively. Mental health problems were rated by the parent using the Child Behavior Checklist (CBCL), a validated and widely-used assessment of childhood behavior 61. The CBCL contains 20 empirically based subscales spanning emotional, social and behavioral domains. Raw scores were used in analyses, higher scores indicating more severe problems. Cognitive functions were assessed by the NIH Toolbox®62, which has good reliability and validity in children63. The toolbox consists of seven different tasks covering episodic memory, executive function, attention, working memory, processing speed, and language abilities, and also provides three composites of crystalized, fluid, and total intelligence64. Uncorrected standard scores were used in analyses. All ten cognitive scores were available at baseline, but only crystalized intelligence was collected two years later.
Social network dataset
In order to validate and generalize the findings in the ABCD study, we utilized a publicly available dataset of a social network experiment, conducted among students in 56 middle schools in New Jersey, the United States26 (https://www.icpsr.umich.edu/web/civicleads/studies/37070). All parents and students provided informed consent for the survey, and the research protocol was approved by the Princeton University Institutional Review Board. Participants were asked to report which other students (up to ten) in their school they chose to spend time with in the last few weeks. We were able to generate a directed friendship network within each school, and therefore identified three kinds of network measures for each individual: outdegree (i.e., outward nomination), indegree (i.e., inward nomination), and reciprocal degree (i.e., reciprocal nomination). Well-being was assessed by three questions: “I feel like I belong at this school”, “I have stayed home from school because of problems with other students”, and “During the past month, I have often been bothered by feeling sad and down”11. Cognitive function was measured by the GPA on a 4.0 scale, obtained from school administrative records.
Structural MRI data
In the ABCD study, 3D T1- and T2-weighted structural images were collected using 3T scanners at 21 data collecting sites65. The detailed preprocessing pipeline has been described elsewhere66. In brief, we used FreeSurfer v6.0 to preprocess the minimal preprocessed T1- and T2-weighted images downloaded from the ABCD study, including cortical surface reconstruction, subcortical segmentation, smoothed by a Gaussian kernel (FWHM = 10 mm), and estimation of morphometric measures (i.e., cortical area, thickness, and volume). Then, the cortical surface of each subject was registered to a standard fsaverage space and parcellated into 180 cortical regions per hemisphere as defined in the Human Connectome Project multimodal parcellation (HCP-MMP) atlas67. Volumetric reconstructions of subcortical structures were also obtained based on the Aseg atlas68.
Neurochemical data
Fourteen receptors and transporters across eight different neurotransmitter systems (serotonin: 5HT1a, 5HT1b, 5HT2a, 5HT4, and 5HTT; dopamine: D1, D2, and DAT; GABA: GABAa; glutamate: mGluR5; norepinephrine: NAT; cannabinoid: CB1; opioid: MOR; acetylcholine: VAChT) were investigated. Density estimates were derived from average group maps of healthy volunteers scanned in prior PET and SPECT studies (Table S5). All density maps were downloaded online (https://github.com/juryxy/JuSpace/tree/JuSpace_v1.3/JuSpace_v1.3/PETatlas), which had been registered and normalized into the Montreal Neurological Institute (MNI) space, and linearly rescaled to 0-10069. For comparability, the HCP atlas in fsaverage space was converted to individual surface space (“mri_surf2surf”) of the MNI brain template ch270 which was preprocessed by Freesurfer (“recon-all”), and then was projected to volume (“mri_label2vol”). The density maps were parcellated into the 360 cortical regions as the structural MRI data according to the volume-based HCP-MMP atlas. Specifically, for the μ-opioid receptor, occipital cortex served as the reference region71, and was therefore excluded in analysis.
Transcriptomic data
Gene expression data was from six neurotypical adult brains in the Allen Human Brain Atlas72. We focused on the opioid receptor genes (i.e., OPRM1 and OPRK1). The preprocessed transcriptomic data were imported from et al73 (https://doi.org/10.6084/m9.figshare.6852911), including probe-to-gene re-annotation, intensity-based data filtering, and probe selection using RNA-seq data as a reference. Then, samples were assigned to brain regions according to the volume-based HCP-MMP atlas, and expression values were averaged within each region. Since right hemisphere data were only available for two donors, analyses were conducted on the left hemisphere only, finally resulting in 177 brain regions.
Statistical analysis
Nonlinear association analysis
The nonlinear associations of close friendship quantity with mental health, cognition, and brain structure were investigated. The quantity of close friendship was log-transformed [log10(x+1)] in analyses as it has a skewed distribution55 (Figure S14). Two different analytic approaches were used to robustly evaluate the nonlinear relationships. First, we fitted a quadratic regression model (y = bx2 + ax + c) with close friendship quantity as the independent variable. Close friendship quantity was mean-centered to ensure that the linear (a) and quadratic (b) terms were orthogonal. Three statistical parameters were of interest: a total F-value of linear and quadratic terms, reflecting the association between close friendship quantity and the measures of interest74; the quadratic term, indicating the presence of a nonlinear association; and the linear term. The effect size of the quadratic term was calculated by the change in the overall proportion of variance (adjusted R2) between the quadratic model and the corresponding linear model, and the effect size of the linear term was the ΔR2 between the linear model and the model with only covariates. The model fits of quadratic and linear models were compared by analysis of variance. Although quadratic regression is widely used in psychosocial studies to detect the presence of nonlinearity11,75, simulation studies showed that this approach for testing a U-shaped effect has a high false positive rate27. Therefore, we conducted a two-lines test27 once a significant quadratic term was found, which could estimate a data-driven breakpoint. We then split the data accordingly to fit two linear models, respectively. If the segment slopes have opposite signs and both of them are significant, a U-shaped relation exists. Same analytic approaches were used in behavioral and neuroimaging analyses. In the ABCD study, sex, age, parent education level, household income, ethnicity, puberty, BMI and site were used as covariates of no interest for the behavioral analyses. For the neuroimaging analyses, we additionally controlled for handedness, head motion and MRI manufacturer. In the social network dataset, we controlled for sex, age, grade, whether the subject was or was not new to the school, and whether or not most friends went to this school. Bonferroni correction was used in behavioral analyses, and FDR correction was used in neuroimaging analyses.
Several sensitivity analyses were performed. To examine the potential sex influence, we conducted nonlinear association analyses in male and female, respectively. The effect of the sex of close friends was tested by separating close friends into same-sex and opposite-sex ones. To validate the findings from data at baseline, we replicated the same analyses using the cross-sectional data collected at 2 years later. For neuroimaging analyses, if significant nonlinear associations were detected, we also conducted linear regression models in two groups split by the average breakpoint, respectively.
Spatial correlation with neurotransmitter density and gene expression
Unthresholded t-statistic maps of brain structure associated with close friendship quantity in two groups (i.e., split by the average breakpoint) were used to correlate with neurotransmitter density and gene expression level by Spearman’s rank correlation. Bootstrapping was performed to ensure the robustness, and the significance was tested by 5,000 times permutation, in which the correlation was re-computed using null t-statistic maps obtained by label shuffling for close friendship quantity76.
Cross-lagged panel analysis
Longitudinal relationships of close friendship quantity with ADHD symptoms (i.e., cbcl_scr_syn_attention) and crystalized intelligence was investigated using a classic two-wave cross-lagged panel model (CLPM) implemented by Mplus 7.0. Firstly, we conducted CLPMs using the absolute value of the difference between close friendship quantity and the breakpoint, and then established CLPMs for participants with the quantity of close friendship ≤ breakpoint and > breakpoint at baseline, respectively. We controlled for several stable (i.e., sex, parent education level, ethnicity, and site) and time-variant variables (i.e., age, household income, and puberty) in these models. The model parameters were estimated by the full information maximum likelihood method77, and the model fit was interpreted using common thresholds of good fit78.
Mediation analysis
We used the baseline data in the ABCD study to test the associations between close friendship quantity, ADHD symptoms, crystalized intelligence, and brain structure. The total area of the significant brain regions was used as the mediator. Variables were normalized and then entered into the model. Sex, age, parent education level, household income, ethnicity, pubertal status, BMI, handedness, head motion, MRI manufacturer and site were used as covariates of no interest. In addition to the total area, the mediation effect of individual significant regions was also evaluated, the p-values of the mediation effect corrected by FDR correction. Total, direct and indirect associations were estimated by bootstrapping 10,000 times, and the 95% bias-corrected and accelerated confidence interval (CI) was reported. Analyses were performed using the R mediation package.
Data Availability
The ABCD data used in this study is Data Release 4.0 (http://dx.doi.org/10.15154/1523041). The data are available by request from the NIMH Data Archive (https://data-archive.nimh.nih.gov/abcd). The social network dataset is conducted among students in 56 middle schools in New Jersey, the United States. This dataset is publicly available (https://www.icpsr.umich.edu/web/civicleads/studies/37070).
Acknowledgements
We thank the children and families whose ongoing participation made this study possible. Data used in the preparation of this article were obtained from the Adolescent Brain Cognitive DevelopmentSM (ABCD) Study (https://abcdstudy.org), held in the NIMH Data Archive (NDA). This is a multisite, longitudinal study designed to recruit more than 10,000 children age 9-10 and follow them over 10 years into early adulthood. The ABCD Study® is supported by the National Institutes of Health and additional federal partners under award numbers U01DA041048, U01DA050989, U01DA051016, U01DA041022, U01DA051018, U01DA051037, U01DA050987, U01DA041174, U01DA041106, U01DA041117, U01DA041028, U01DA041134, U01DA050988, U01DA051039, U01DA041156, U01DA041025, U01DA041120, U01DA051038, U01DA041148, U01DA041093, U01DA041089, U24DA041123, U24DA041147. A full list of supporters is available at https://abcdstudy.org/federal-partners.html. A listing of participating sites and a complete listing of the study investigators can be found at https://abcdstudy.org/consortium_members/. ABCD consortium investigators designed and implemented the study and/or provided data but did not necessarily participate in the analysis or writing of this report. This manuscript reflects the views of the authors and may not reflect the opinions or views of the NIH or ABCD consortium investigators. C.S. is supported by grants from the National Natural Sciences Foundation of China (No. 82101617). W.C. is supported by grants from the National Natural Sciences Foundation of China (No. 82071997), the Shanghai Rising-Star Program (No. 21QA1408700). J.F. is supported by National Key R&D Program of China (No. 2018YFC1312904 and No. 2019YFA0709502), Shanghai Municipal Science and Technology Major Project (No. 2018SHZDZX01), ZJ Lab, Shanghai Center for Brain Science and Brain-Inspired Technology, and the 111 Project (No. B18015).
Footnotes
Competing Interest Statement: The authors declare no competing interests.
References
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵