Abstract
Background Research has demonstrated associations between pubertal development and brain maturation. However, existing studies have been limited by small samples, cross-sectional designs, and inconclusive findings regarding coupling directionality and sex differences.
Methods We examined the longitudinal temporal coupling of puberty status assessed using the Pubertal Development Scale (PDS) and magnetic resonance imaging (MRI)-based brain grey and white matter structure. Our sample consisted of 8,896 children and adolescents at baseline (mean age = 9.9) and 6,099 at follow-up (mean age = 11.9) from the Adolescent Brain and Cognitive Development (ABCD) Study.
Results Applying multigroup Bivariate Latent Change Score (BLCS) models, we found that baseline pubertal status predicted the rate of change in cortical thickness among females only, while it predicted rate of change in cortical surface area for both males and females, with a significantly stronger association for females. For cortical surface area, we also found a correlation between baseline pubertal status and area and co-occurring changes over time, with both associations present in males only. Diffusion tensor imaging (DTI) analysis revealed correlated change between pubertal status and fractional anisotropy (FA) for both males and females, but no significant associations for mean diffusivity (MD).
Conclusions Our results suggest that pubertal status in early adolescence predicts cortical maturation, and that the strength of the associations differ between sex. Further research is needed to understand how the associations between puberty and brain maturation are related to environmental and lifestyle factors, and the impact on mental health.
1. Introduction
Adolescence is a transitional stage from childhood to adulthood that is marked by prolonged tissue-specific neurodevelopmental processes. Tethered to this, puberty has been highlighted as a defining period for protracted brain structural maturation during adolescence, with a range of studies linking concurrent trends of puberal and brain development that go beyond the effects of age (Ando et al., 2021; Blakemore et al., 2010; Herting & Sowell, 2017; Vijayakumar et al., 2021). However, existing studies on puberty-brain associations are limited by small sample sizes and cross-sectional designs, in addition to mixed findings, underscoring the need for large-scale longitudinal studies that examine the temporal coupling of pubertal development and brain maturational processes.
Although pubertal development is characterised by progress through the same stages for each individual — involving changes in secondary sex characteristics, rising hormone levels, and emerging reproductive capability — there is a large variation in the timing and tempo of these changes. There are both inter-individual and between-sex differences involved, with the typical age of pubertal onset ranging from 9 to 14 years in human males and 8 to 13 years in human females (Herting et al., 2021; Savin-Williams & Ream, 2006).
The variation in pubertal timing, defined as an individual’s relative pubertal developmental stage compared with same -age and -sex peers (Ullsperger & Nikolas, 2017), has a substantial genetic contribution among both males and females, with heritability estimates ranging from 50 to 80% (Stroud & Davila, 2011). However, a growing body of evidence has also linked variation in pubertal timing to lifestyle and environmental factors, including childhood obesity (Ahmed et al., 2009; Freedman et al., 2002; Herting et al., 2021), socioeconomic status (SES) (Braithwaite et al., 2009; Herting et al., 2021), and nutrition (Parent et al., 2003). An implication of this is that pubertal timing may be partially malleable and, as such, represents a window of opportunity for early intervention that targets youth at risk of health issues related to pubertal timing.
For example, in females, earlier pubertal onset is associated with an increased risk of type II diabetes, breast cancer, and adult obesity (Ahmed et al., 2009), in addition to mental outcomes related to risk-taking (Braams et al., 2015; Collado-Rodriguez et al., 2014), depression and anxiety (Copeland et al., 2019; McNeilly et al., 2022; Mendle et al., 2010; Mendle & Ferrero, 2012), and substance use (Marceau et al., 2019; Patton et al., 2004; Stumper et al., 2019). In males, health risks are better understood by the ‘off-time’ hypothesis (either very early or very late pubertal onset), with worse mental health outcomes (Hayward, 2003), such as depression, anxiety, attention deficit/hyperactivity disorder (ADHD), and risk-taking behaviour being associated with both early (Ge et al., 2001; Huddleston & Ge, 2003; Mendle & Ferrero, 2012) and late (Conley & Rudolph, 2009; Graber et al., 1997) pubertal onset.
In conjunction to pubertal development, the brain shows continued tissue-specific changes that are prone to influence by fluctuating environmental pressures (Ferschmann et al., 2022; Tooley et al., 2021) beyond the complex genetic contribution (Grasby et al., 2020). Furthermore, puberty has been implicated as an influencing factor to the rate of maturation in the brain, and as such recommended as non-negligible when studying adolescent brain maturation (Holm et al., 2022; Wierenga et al., 2018).
Previous cross-sectional (Bramen et al., 2011; Paus et al., 2010; Pfefferbaum et al., 2016) and longitudinal (Herting et al., 2015; Nguyen et al., 2013; Vijayakumar et al., 2021) studies have reported reductions in global grey matter volume and cortical thickness associated with higher pubertal status and testosterone levels, mirroring normative developmental patterns of grey matter during adolescence (Mills et al., 2016; Tamnes et al., 2017). However, when controlling for age, findings have been less consistent, with some studies reporting associations between decreased cortical grey matter and more pubertally developed females (Bramen et al., 2011; Peper, Brouwer, et al., 2009), others reporting positive associations between grey matter volume and testosterone in males (Peper, Brouwer, et al., 2009), and several reporting null findings for both sexes (Bramen et al., 2011; Koolschijn et al., 2014; Peper, Brouwer, et al., 2009; Peper, Schnack, et al., 2009).
For white matter, numerous cross-sectional studies have reported positive associations between pubertal status and white matter volumes (Chavarria et al., 2014; Perrin et al., 2009; Pfefferbaum et al., 2016), with some conflicting findings (Peper, Schnack, et al., 2009). Similarly, studies have reported positive associations between testosterone levels and white matter volume among boys (Paus et al., 2010; Perrin et al., 2008), with inconsistent findings when accounting for age (Peper, Brouwer, et al., 2009).
Diffusion tensor imaging (DTI) provides information regarding white matter architecture and microstructure. DTI studies have generally reported higher fractional anisotropy (FA) and lower mean diffusivity (MD) among more pubertally developed adolescents (Herting et al., 2012; Menzies et al., 2014), mirroring normative developmental patterns (Lebel & Deoni, 2018; Pfefferbaum et al., 2016). However, to date, only one longitudinal study (Herting et al., 2017) has investigated associations between puberty and white matter microstructure using DTI, albeit with a small sample, finding that physical pubertal changes predicted changes in white matter. However, further longitudinal research is necessary to validate and extend previous findings and to examine potential sex differences.
Using a large-scale longitudinal sample (N = 8,896 at baseline; N = 6,099 at follow-up) of children and adolescents between the ages of 9-14 years from the Adolescent Brain Cognitive Development (ABCD) Study cohort, we investigated the developmental associations between puberty and brain structure and tested whether and to what extent these associations differ between males and females. By leveraging a powerful Structural Equation Modeling (SEM) framework, multigroup Bivariate Latent Change Scores (BLCS), we examined how individual and sex differences in pubertal development relate to differences in structural brain maturation. Importantly, BLSC allows us to test the temporal dynamics of the puberty-brain relationship, specifically whether pubertal status and brain features correlate at baseline, brain structure at baseline predicts change in pubertal status, baseline pubertal status predicts change in brain structure, and whether changes in puberty and brain structure co-occur over the course of the study period.
Based on previous research, we expected cross-sectional and longitudinal associations between pubertal status measured with the Pubertal Development Scale (PDS) and global features of cortical structure (surface area, thickness) and white matter microstructure (FA, MD) in the examined age range. Specifically, we hypothesised correlations at baseline across all brain MRI features and for changes in pubertal status and brain structure to co-occur. In terms of coupling effects (i.e., regressions), we hypothesised that baseline pubertal status would predict change in global brain structure and not vice versa. Additionally, we hypothesised sex related effects to be present, with generally stronger puberty-brain associations in females compared to males.
2. Methodology
2.1. Description of sample
The initial sample consisted of children and adolescents that are part of the ongoing longitudinal ABCD Study, including a baseline cohort of ∼11,800 nine- and ten-year-olds to be followed for the course of ten years (Garavan et al., 2018). Data used in the present study were drawn from the ABCD curated annual release 4.0, containing data from baseline up until the second-year visit (https://data-archive.nimh.nih.gov/abcd). All ABCD Study data is stored in the NIMH Data Archive Collection #2573, which is available for registered and authorised users (Request #7474, PI: Westlye). The 4.0 release will be permanently available as a persistent dataset defined in the NDA Study 1299 and has been assigned the DOI 10.15154/1523041.
To minimize confounding effects from complex family-related factors, unrelated participants were chosen by randomly selecting an individual from each family ID, subsequently excluding any siblings (see Supporting Information (SI) Section 1 for details). Additionally, participants were excluded using the ABCD Study exclusion criteria, which included non-English proficiency, general MRI contraindications, a history of a major neurological disorder, traumatic brain injury, extreme premature birth (<28 weeks gestational age), a diagnosis of schizophrenia, intellectual disability, moderate to severe autism spectrum disorder, or substance abuse disorder (Karcher et al., 2018).
Following procedures of quality control (see SI Section 1) for each brain measure (cortical thickness, surface area, FA, and MD), PDS, and our covariates of interest including SES, body-mass index (BMI), and genetic ancestry-derived ethnicity, the final sample consisted of N = 8,896 participants at baseline (4,223 females, 4,673 males) and N = 6,099 participants at follow-up (2,798 females, 3,301 males). Table 1 displays means and standard deviations of each brain MRI metric, PDS, and covariates of interest, split by sex, while the age and PDS distributions in the sample are presented in Figure 1. Ethical approvals for the study are described in SI Section 2.
Descriptive statistics of the sample, including mean ± standard deviation (SD) for each measure.
Age (left) and PDS (right) distribution of the study sample, split by sex and timepoint. Dotted line between plots connects the mean age and PDS values at baseline (timepoint 1) and follow-up (timepoint 2). Age and total PDS scores shown as distributions where green represents females and purple represents males.
2.2. MRI acquisition and processing
Neuroimaging data were acquired at 21 different sites and processed by the ABCD team. A 3-T Siemens Prisma, General Electric 750 or Phillips scanner was used for data acquisition. Protocols used for data acquisition and processing are described elsewhere (Casey et al., 2018; Hagler et al., 2019). In brief, T1-weighted data was acquired by magnetisation-prepared rapid acquisition gradient echo scans with a resolution of 1×1×1 mm3, which was used for generating cortical structural measures, and diffusion-weighted data was obtained by high angular resolution diffusion imaging scans, used for generating white matter microstructural measures.
Two modalities of brain structural measures were used in the present study: grey matter cortical measures and white matter microstructural measures (Hagler et al., 2019). Cortical reconstruction and volumetric segmentation was performed with FreeSurfer 5.3.0 (Dale et al., 1999; Fischl et al., 2002). White matter microstructural measures (DTI) were generated using AtlasTrack, a probabilistic atlas-based method for automated segmentation of white matter fiber tracts (Hagler Jr. et al., 2009). For cortical measures, global measures were generated for cortical surface area and thickness. For DTI, measures of FA and MD were generated over the whole brain. For full pipeline, see SI Section 3.
2.3. Multi-site MRI scanner effects
To adjust for non-biological variance introduced by multi-scanner MRI and acquisition protocols (i.e., scanner effects), a longitudinal harmonisation technique was carried out using the R package LongComBat (Beer et al., 2020). To reduce bias in scaling and spatial resolution of features, our four global brain MRI measures (FA, MD, thickness, area) were harmonised separately. Since the harmonisation is applied to model residuals, the longitudinal ComBat model should ideally match the model in the final analysis. As such, each model was fitted with formulas including age, PDS, BMI, SES, and genetic ancestry-derived ethnicity as covariates. SI Figures 3 and 4 show distribution of each MRI measure before harmonisation and after longitudinal ComBat.
2.4. Pubertal development assessment
Pubertal status was measured using the PDS (Petersen et al., 1988) to examine the perceived development of secondary sex characteristics such as growth spurts, body hair growth, skin changes, breast development and menarche in girls, and voice changes and growth of testes in boys. The PDS is a questionnaire designed to mimic the traditional Tanner staging assessment without the use of reference pictures (Petersen et al., 1988). Due to previous research showing that youth tend to over-report their perceived physical development at younger ages (Schlossberger et al., 1992), the current study utilised caregiver PDS report.
The PDS includes five-items, each rated on a 4-point scale (1 = no development; 2 = development has barely begun; 3 = development is definitely underway; and 4 = development is complete; except menstruation, which is coded 1 = has not begun, 4 = has begun). Thus, higher scores reflect more advanced pubertal development. The PDS has shown high inter-rater reliability between parent and self-rated assessment to clinicians, and correlates highly with plasma levels of gonadal hormones (Carskadon & Acebo, 1993; Koopman-Verhoeff et al., 2020) in addition to the Tanner stages (Koopman-Verhoeff et al., 2020). We used total PDS score for our analyses. When describing our results, PDS and the term pubertal status is used interchangeably for cross-sectional analyses, while pubertal development is used for longitudinal analyses.
2.5. Covariates of interest
Previous research has highlighted several factors that potentially contribute to pubertal timing and development, including BMI, SES, and ethnicity (Ahmed et al., 2009; Braithwaite et al., 2009; Deardorff et al., 2019; Freedman et al., 2002; Herting et al., 2021; Kaplowitz et al., 2001). Therefore, we included BMI, SES, and genetic ancestry-derived ethnicity, in addition to age, as covariates. For full details on how each covariate was measured and derived, see SI Section 4. For figures relating to distribution, frequency, and associations between covariates, see SI Figures 5-7.
2.6. Statistical analysis
The code used in the present study is available in a public repository - Open Science Framework (OSF) - and accessible directly through the OSF webpage https://osf.io/nh649/.
All analyses were carried out using R version 4.0.0 (R Core Team, 2022). A multigroup BLCS model was tested using the lavaan package (Rosseel, 2012) following recommendations provided by (Kievit et al., 2018) for model setup. Four separate models were set up to examine 1) correlations at baseline between pubertal status and MRI measures of global FA, MD, cortical thickness, and surface area; 2) to what extent the baseline pubertal status predicted rate of change in global brain MRI measures; 3) to what extent the baseline score of global brain MRI measure predicted rate of change in the pubertal development; and lastly, 4) whether change co-occurred. To investigate sex differences in the associations, we used multiple group analyses in the framework of SEM, with sex (male, female) as the grouping variable.
We quantified to what extent change in pubertal status (ΔPDS) between baseline and follow-up was a function of brain MRI measure (γ2) and pubertal status (β1) at baseline as following:
And to what extent change in brain MRI measure (MRI) between baseline and follow-up was a function of pubertal status (γ1) and brain MRI measure (β2) at baseline as following:
For an illustrative description of this model, see Figure 2. Due to their relatedness to pubertal development, age, BMI, SES, and genetic ancestry-derived ethnicity were included in each model. Full information for maximum likelihood estimation procedures were used for all BLCS analyses, thereby providing missing data routines that are considered to be adequate (Enders, 2010).
Bivariate latent change score (BLCS) model with two timepoints, examining correlation between pubertal status (PDS) and brain MRI metric (MRI) at baseline (pink), to what extent pubertal status at baseline predicts the rate of change in brain MRI metric (blue: γ1), to what extent brain MRI metric at baseline predicts the rate of change in pubertal status (blue: γ2), and whether changes in PDS and MRI co-occur (yellow) after considering the coupling (γ1, γ2) pathways. Image adapted from Kievit et al. (2018) with permission from author.
Analyses to test sex differences were also carried out, where the model parameters above (1-4) were constrained to be equal across the two groups and tested by means of χ2-difference tests against their respective non-constrained (original) model for each of the four brain MRI metrics. To test sex differences for all associations for all four models, sixteen such tests were carried. The significance threshold was set at p < 0.05, and the results were corrected for multiple comparisons using the false discovery rate (FDR) adjustment (Benjamini & Hochberg, 1995). Both uncorrected (p) and corrected (pcorr) values are reported for these tests.
3. Results
3.1. Descriptive statistics
Descriptive statistics can be found in Table 1. Brain MRI age curves are presented in SI Figure 8. Brain MRI trajectories as a function of PDS are presented in SI Figure 9.
3.2. Multigroup bivariate latent change score (BLCS) models
The main results of the models are summarised in Table 2. For the full models, including analyses testing for sex differences in the associations, see SI Section 5.
Multigroup BLCS model for each brain MRI metric, where numbers 1-4 represent the following effects: (1) correlation at baseline, (2) regression (PDS at baseline related to change in MRI metric), (3) regression (MRI metric at baseline related to change in PDS), (4) correlation of change score.
3.3. Pubertal development and fractional anisotropy
First, a BLCS for pubertal development and FA was modelled and the model fit was good: χ2 (28) = 272.45, CFI = 0.956, RMSEA = 0.044. Inspection of the four parameters of interest, reflecting the four possible puberty-brain relationships (outlined in Section 2.6), showed no statistically significant correlation between pubertal status and FA at baseline for males or females. The results did however show statistically significant associations of correlated change for both males (r = .048, p = .018) and females (r = .079, p < .01). In other words, those with greater increase in pubertal status were, on average, those with greater increase in FA. There were no coupling effects to report. Next, χ2-difference testing for sex differences were carried out. Here, no significant difference between males and females were found for any of the four puberty-brain relationships. For a visual representation of the associations, see Figure 3.
Visualisation of correlated change, where a change in FA is positively related to change in PDS for both males and females. Delta scores calculated as follow-up FA and PDS measures minus baseline FA and PDS measures, respectively. Males represented by lilac triangles and females by green circles.
3.4. Pubertal development and mean diffusivity
Second, a BLCS for pubertal development and MD was modelled and the model fit was good: χ2 (28) = 267.86, CFI = 0.959, RMSEA = 0.043. Inspection of the four parameters of interest revealed no statistically significant correlations or regression coefficients. χ2-difference tests revealed no significant difference between males and females for any of the four puberty-brain relationships.
3.5. Pubertal development and cortical thickness
Third, a BLCS for pubertal development and cortical thickness was modelled and the model fit was good: χ2 (28) = 274.47, CFI = 0.963, RMSEA = 0.044. When examining the four parameters of interest, no statistically significant associations were observed between pubertal status and cortical thickness at baseline for males or females. The results did however reveal statistically significant correlated change for females (r = -.058, p < .01), with no association present for males (r = -.001, p = .94). In terms of coupling effects, pubertal status at baseline predicted rate of change in cortical thickness in females (β = -0.105, p < .01), while no associations were present in males. In other words, having initially higher PDS scores at baseline predicted larger increase between baseline and follow-up cortical thickness among females. χ2-difference tests revealed a significant difference between males and females ((df, 1) = [8.75], p < 0.01. pcorr = 0.049) when comparing model parameters of pubertal status at baseline predicting rate of change in cortical thickness. The three other relationships investigated revealed no significant sex differences after correction for multiple comparisons. Results are visualised in Figure 4.
Visualisation of cortical thickness effects. Figure on the left visualising associations for females, where pubertal status at baseline is associated with change in cortical thickness. Figure on the right visualising correlated change for females.
3.6. Pubertal development and surface area
Lastly, a BLCS for pubertal development and surface area was modelled and the model fit was good: χ2 (28) = 275.74, CFI = 0.982, RMSEA = 0.044. Inspection of the four parameters of interest showed statistically significant (positive) associations for the correlation between pubertal status and surface area at baseline (r = .055, p < .01) and correlated change for males (r = .083, p < .01), with no significant associations among females for either correlated baseline (r = -.001, p = .97) or change (r = .002, p = .91). Coupling effects indicated pubertal status at baseline as being predictive of rate of change in surface area for both males (β = -.078, p < .01) and females (β = -.210, p < .01). In other words, having initially higher PDS scores at baseline predicted larger decrease between baseline and follow-up surface area measures. χ2-difference tests found a statistically significant sex difference ((df, 1) = [12.26], p < 0.01, pcorr < 0.01) when comparing model parameters of pubertal status at baseline predicting rate of change in surface area, thereby showing that the association was significantly stronger for females than males. The three other relationships investigated revealed no significant group differences after correction for multiple comparisons. See Figure 5 for a visual representation for the results.
Visualisation of surface area effects. Top left side figure shows correlation between baseline surface area and PDS, observed in males only. Bottom left side figure shows correlated change, with effects present in males only. Top right-side figure shows association between pubertal status at baseline and change in surface area between the two timepoints, with effects present in both males and females.
4. Discussion
We investigated the coupling pathways of pubertal development and brain maturation in early adolescence using a multi-group BLCS framework in a large-scale longitudinal sample. Our main findings indicate that baseline pubertal status predicted the rate of change in cortical thickness for females only, and that pubertal status predicted rate of change in surface area for both sexes, but more strongly for females than males. Associations were also present for both correlation between pubertal status and surface area at baseline, and co-occurring changes over time, but these associations were present among males only. DTI analyses revealed correlated change over time between pubertal development and FA for both males and females, and no significant associations for MD. Overall, our results suggest sex-related differences in the associations between baseline pubertal status and rate of change in global cortical grey matter features of the brain.
Our study utilised a solid statistical framework to investigate the longitudinal coupling effects of pubertal development and structural brain maturation, providing valuable insight into the temporal dynamics of brain-and-puberty development. Previous research has reported sex-related associations between higher pubertal stage and global cortical grey matter (Herting et al., 2015), with studies generally reporting reduced cortical thickness in females (Bramen et al., 2011; Pfefferbaum et al., 2016) and no associations for males (Bramen et al., 2011; Koolschijn et al., 2014; Peper, Brouwer, et al., 2009). Our findings support previous findings and reveal additional knowledge into the coupling effects, with baseline pubertal status predicting rate of change in cortical thickness, as hypothesised.
For global surface area, baseline levels of pubertal status and surface area were positively related in males, going against our hypothesis. However, male participants who were more pubertally developed at baseline showed less increases in surface area compared to male participants who were less pubertally developed at baseline. Interestingly, our coupling effects provide evidence suggesting that pubertal status at baseline predicts rate of change in surface area (for both males and females), with faster pubertal tempo associated with greater decrease in surface area, supporting our hypothesis. Moreover, despite both sexes showing effects, our findings revealed significantly stronger effects in females than males.
For surface area, our findings can be best understood in the context of the age range in this sample. While the directionality of the coupling effects go in the expected directions (i.e., more pubertally developed predicting surface area decrease), one unexpected finding is the observation that males show the opposite trend for baseline levels of pubertal status and surface area (i.e., positively related). Previous research has reported decreases in surface area associated with more pubertally developed males (Herting et al., 2015). One possible explanation for this could be that many of the males in the present sample are still experiencing increase in surface area (see e.g., Raznahan et al. (2011) and SI Figure 8).
Additionally, as evidenced by Figure 1 and PDS descriptive statistics in Table 1, the majority of the males in our sample are in the early stages of pubertal development compared to females (mid-late stages), and do not progress a lot throughout the study period, with the average PDS moving from 7.5 at baseline to 8.8 at follow-up. These factors highlight a need for future research that includes an age range spanning the entirety pubertal development.
Previous research investigating associations between pubertal development and DTI is sparse. The one existing longitudinal DTI study (Herting et al. 2017) reported that increase in PDS was associated with increase in FA, which was supported in our study using a much larger sample. However, while the prior study reported sex-related effects, our results did not indicate differences between males and females. One potential explanation for this is that Herting and colleagues looked at various regions of the brain while our study observed global FA, meaning if a minority of regional effects (as observed in Herting et al. (2017)) were going the other direction, they may not be evident in the overall effect of increasing FA in the majority of other regions. No evidence of correlation at baseline or coupling effects were observed in either FA or MD, and there were no significant sex differences, which goes against our previously stated hypotheses.
Adolescence is a sensitive period that offers a window of opportunity to shape the fate of the adult brain. Within that period, accelerated pubertal development has been identified as a potential indicator of concern for the developing brain, with earlier pubertal timing being linked to psychological outcomes such as depression, anxiety, and ADHD (Copeland et al., 2019; Ge et al., 2006; Mendle et al., 2010; Mendle & Ferrero, 2012). As such, addressing early pubertal-timing-related risk factors may be crucial in reducing or even halting the development of potentially related psychopathology. Explicitly investigating the brain-puberty relationship with sex differences in mind is also important considering that risk of worse mental health outcomes differs for females (earlier pubertal timing) and males (off-time hypothesis), and this may have ramifications for related brain changes.
With research showing earlier pubertal timing in those facing barriers to positive health and well-being, questions must be raised of how to tackle the challenges posed by environmental factors and disparities in living conditions. Research has observed earlier pubertal timing in concordance with higher child weight status, lower SES, and nutritional needs not being met, and these have disproportionately impacted minority groups more. Although speculative, environmental pressures may be central agents in advanced adaptation to adult roles in the environment, termed ‘adultification’ (Pfeifer & Allen, 2021), where adolescents may be undergoing earlier onset and tempo of puberty, and related neurodevelopmental changes, as a means to adapt to a harsh environment, subsequently increasing their risk of psychopathology. More research is needed to elucidate the potential role of environmental factors on advanced pubertal development and brain maturation.
The current study has several strengths and limitations that warrant further discussion. Firstly, the ABCD Study cohort provides a well-powered sample that is openly available for researchers, improving both transparency and encouraging initiatives testing the robustness of the findings. Second, the longitudinal data and analyses allow for tracking changes over time, albeit while still limited by only one follow-up assessment. Nevertheless, the current multigroup BLCS analysis framework is well tailored to address hypotheses related to both general and fine-grained temporal dynamics over time and is especially powerful for testing cross-domain couplings (Kievit et al., 2018).
While global measures may capture broad effects, we know that pubertal development is associated with regional changes in frontal, temporal and parietal areas (as well as subcortical regions) which provides rationale for looking at these in future research. FA and MD revealed fewer effects compared to grey matter morphology. Future research should utilise advanced diffusion MRI metrics (Beck et al., 2021; Tamnes et al., 2018) that differentiate between intra- and extra-axonal compartments and potentially provides more information than conventional DTI. Lastly, future research would benefit from a combined approach of assessing pubertal development by means of questionnaires as well as hormonal data, given the limitations of PDS (Cheng et al., 2021).
Some general limitations pertaining to ABCD Study and other large-scale volunteer scale databases include biases related to being more healthy, more wealthy, and more educated compared to the general population in their respective countries (Brayne & Moffitt, 2022). This lack of representativeness may mean that those with most adversity and disadvantage will not be represented in this research. Future research should consider more diversified samples to reduce the healthy volunteer bias compromising the generalisability of the findings to more vulnerable populations. On the other hand, the ABCD Study cohort does include more ethnically diverse participants than similar large-scale volunteer databases (Brayne & Moffitt, 2022). Due to participant ages at recruitment, the ABCD Study data does not capture early gonadal processes for many participants (Cheng et al., 2021; Herting et al., 2021), which is more pronounced for females and non-White participants.
In conclusion, our analyses highlight pubertal status in early adolescence as being predictive of longitudinal changes in cortical grey matter structure over time, and that these changes are different for males and females. Further research is needed to understand how puberty impacts the brain, the contribution of environmental and lifestyle factors on pubertal development and brain maturation, and, in turn, the effects on mental health.
Data Availability
Data used in the present study were drawn from the ABCD curated annual release 4.0, containing data from baseline up until the second-year visit (https://data-archive.nimh.nih.gov/abcd). All ABCD Study data is stored in the NIMH Data Archive Collection #2573, which is available for registered and authorised users.
5. Funding
This work was supported by the Research Council of Norway (#223273, #288083, #323951, #276082, #249795, #248238, #298646), and the South-Eastern Norway Regional Health Authority (#2019069, #2021070, #500189, #2014097, #2015073, #2016083, #2018076), KG Jebsen Stiftelsen, and the European Research Council under the European Union’s Horizon 2020 Research and Innovation program (ERC StG, #802998 and RIA #847776).
6. Acknowledgments
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 aged 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.
7. References
