Growth in early infancy drives optimal brain functional connectivity which predicts cognitive flexibility in later childhood

  1. Chiara Bulgarelli  Is a corresponding author
  2. Anna Blasi
  3. Samantha McCann
  4. Bosiljka Milosavljevic
  5. Giulia Ghillia
  6. Ebrima Mbye
  7. Ebou Touray
  8. Tijan Fadera
  9. Lena Acolatse
  10. Sophie E Moore
  11. Sarah Lloyd-Fox
  12. Clare E Elwell
  13. Adam T Eggebrecht
  14. The BRIGHT Study Team
  1. Centre for Brain and Cognitive Development, Birkbeck, University of London, United Kingdom
  2. Department of Medical Physics and Biomedical Engineering, University College London, United Kingdom
  3. Department of Women and Children’s Health, King's College London, United Kingdom
  4. Medical Research Council Unit The Gambia at the London School of Hygiene and Tropical Medicine, Gambia
  5. Department of Psychology, University of Cambridge, United Kingdom
  6. School of Biological and Experimental Psychology, Queen Mary University of London, United Kingdom
  7. Department of Women and Children’s Health, University of Liverpool, United Kingdom
  8. Nutrition Innovation Centre for Food and Health, School of Biomedical Sciences, Ulster University, Ireland
  9. Mallinckrodt Institute of Radiology, Washington University School of Medicine in St. Louis, United States
9 figures, 8 tables and 1 additional file

Figures

Experimental design.

(A) Measures taken in the BRIGHT project used in this work. The measuring tape represents anthropometric measures, the brain represents fNIRS FC and the test represents the cognitive flexibility assessment. (B) Schematic representation of the spatial layout of the fNIRS array. Sources are marked with red stars, detectors are marked with blue circles, channels are marked with grey lines and numbered with black circles. The channels/optodes used as a reference for the tragus are highlighted in green.

Linear mixed models results showing FC that displayed a statistically significant change with age.

(A) Significant results of the linear mixed model, blue indicates connections that decreased with age, red indicates connections that increased with age. (B) Mean and standard error of the mean (SE) of the functional connections that changed with age (HbO2). Error bars are 1 SE. (C) Violin plot showing the mean ± SD (red circles and lines) and the individual variability (coloured dots) of the FC that showed a change with time (HbO2).

Associations between FC, early growth and later cognitive flexibility.

(A) Significant positive associations are in green, significant negative associations are in orange, and non-significant associations are in blue. * indicates regressions that survived FDR correction for multiple comparisons. (B) Schematic representations of the early FC connections shown to predict cognitive flexibility in preschoolers. Significant positive associations are in green, significant negative associations are in orange (these results did not survive corrections for multiple comparisons).

Appendix 1—figure 1
Schematic representation of the fNIRS array and the functional connections tested.

(A) Each dot represents a channel, colours on the left plot correspond to colours on the right on the baby’s head. Six sections. (B) The 21 connections tested in the linear mixed models. Interhemispheric homotopic connections are in orange (connecting the same regions between hemispheres, that is front left with front right), intrahemispheric connections within section are in green (correlations of channels belonging to the same region), fronto-posterior are in blue (connecting front and middle, middle and back, and front and back regions of the same hemisphere), and crossing interhemispheric connections (interhemispheric non-homotopic, connecting the front and middle, middle and back, and front and back regions of the two hemispheres) are in yellow.

Appendix 1—figure 2
Linear mixed models result showing FC that displayed a statistically significant change with age (fNIRS pre-processing without global signal regression).

(A) Results of the linear mixed model, blue indicates connections that decreased with age, red indicates connections that increase with age. (B) Mean and SE of the functional connections that changed with age. Error bars are 1 SE.

Appendix 1—figure 3
Some examples of scatterplots showing the association between ΔWLZ and FC at 24 months.

The black lines represent the line of best fit and the 95% confidence intervals.

Appendix 1—figure 4
fNIRS preprocessing steps.

The column ‘infants included in the analyses’ in Table 3 refers to those participants whose data survived these preprocessing steps.

Appendix 1—figure 5
The effect of different thresholds for GVTD for motion detection and minimum valid data after pre-processing on data inclusion.

(A) Percentage of data included (left) and seconds of data included (right) per participant. Each line represents an infant, the black line represents the mean value. These graphs are reported from the 12 months sample as example. (B) Number of infants included in the LMM by varying STD threshold for GVTD and minimum length of valid data at different ages.

Appendix 1—figure 6
Strength of correlation between FC in the first and last portion of data.

This graph is reported from the age 12 month sample as example.

Tables

Table 1
FC that significantly changed with age.

Results are displayed in terms of estimated betas, standard errors and p values that survived Bonferroni correction. Regressions that survived Bonferroni correction for multiple comparisons are in bold.

FCFpBaseline(5 months)Betas (SE), p5–8 change Betas (SE), p5–12 change Betas (SE), p5–18 change Betas (SE), p5–24 change Betas (SE), p
HbO2
Frontal interhemispheric11.03<0.0010.28 (0.03),<0.001–0.09 (0.06), 0.12–0.20 (0.05),<0.001–0.25 (0.05),<0.001–0.33 (0.05),<0.001
Left fronto-middle5.8<0.001–0.28 (0.04),<0.0010.06 (0.06), 0.2930.09 (0.05), 0.1020.19 (0.05), <0.0010.23 (0.05), <0.001
Right fronto-middle4.86<0.001–0.70 (0.29), 0.0180.01 (0.05), 0.9820.05 (0.05), 0.2880.13 (0.05), 0.0070.19 (0.05), <0.001
Right frontal-posterior5.52<0.001–0.20 (0.03),<0.0010.16 (0.06), 0.0050.14 (0.05), 0.0060.20 (0.05), <0.0010.22 (0.05), <0.001
HHb
Frontal interhemispheric6.59<0.0010.32 (0.03),<0.001–0.09 (0.06), 0.107–0.22 (0.05),<0.001–0.23 (0.05),<0.001–0.22 (0.05),<0.001
Cross left fronto-right middle5.17<0.001–0.32 (0.04),<0.0010.06 (0.06), 0.3280.07 (0.05), 0.1880.19 (0.05), <0.0010.22 (0.05), <0.001
Right fronto-middle5.55<0.001–0.78 (0.32), 0.0360.01 (0.06), 0.860.05 (0.05), 0.3010.16 (0.05), 0.0020.21 (0.05), <0.001
Right frontal-posterior5.38<0.001–0.24 (0.04),<0.0010.10 (0.06) 0.0730.08 (0.05) 0.1330.23 (0.05), <0.0010.02 (0.05), <0.001
Table 2
Results of the regression analyses of the effect of ΔWLZ on FC at 24 months.

Significant positive associations are in green, significant negative associations are in orange, and non-significant (NS) associations are in blue; * indicates regressions that are still significant after correcting for HCAZ at 7/14 days and ** indicates regressions that are still significant after correcting for neonatal HCAZ and WLZ. Regressions that survived FDR correction for multiple comparisons are in bold.

Frontal interhemispheric FC at 24 months
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthNSF(1,68)=8.69, P=0.004*, R2=0.102F(1,67)=7.65, P=0.007*, R2=0.104F(1,69)=10.4, P=0.002**, R2=0.121F(1,68)=10.98, P=0.001**, R2=0.142F(1,59)=5.40, P=0.024*, R2=0.085
1 monthF(1,77)=5.59, P=0.021**, R2=0.069NSF(1,78)=4.74, P=0.032**, R2=0.058F(1,75)=4.45, P=0.038**, R2=0.057NS
5 monthsNSNSNSNS
8 monthsNSNSNS
12 monthsNSNS
18 monthsNS
24 months
Left fronto-middle FC at 24 months
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthNSNSF(1,67)=4.07, P=0.048**, R2=0.058NSNSF(1,59)=4.24, P=0.044**, R2=0.068
1 monthNSNSNSNSNS
5 monthsNSNSNSNS
8 monthsNSNSNS
12 monthsNSNS
18 monthsNS
24 months
Right fronto-middle FC at 24 months
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthNSNSF(1,67)=4.74, P=0.033, R2=0.067NSF(1,68)=5.29, P=0.024**, R2=0.073NS
1 monthF(1,77)=4.26, P=0.042**, R2=0.053F(1,75)=6.66, P=0.012**, R2=0.083NSF(1,75)=5.33, P=0.024**, R2=0.067NS
5 monthsNSNSNSNS
8 monthsNSNSNS
12 monthsNSNS
18 monthsNS
24 months
Right frontal-posterior FC at 24 months
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthF(1,66)=8.26, P=0.006**, R2=0.113NSNSNSNSNS
1 monthF(1,69)=4.86, P=0.031*, R2=0.067F(1,68)=6.98, P=0.01**, R2=0.094F(1,71)=8.86, P=0.004**, R2=0.112F(1,67)=7.72, P=0.007**, R2=0.105F(1,57)=7.03, P=0.01**, R2=0.112
5 monthsNSNSNSNS
8 monthsNSNSNS
12 monthsNSNS
18 monthsNS
24 months
Table 3
Demographic information per each age.
AgeAge in days (mean ± SD)Sex (M, F)WLZ (mean ± SD)HCZ (mean ± SD)
5 months158.92±9.9145,42–0.24±0.95–0.72±0.92
8 months249.34±17.9728,25–0.24±1.01–0.81±0.90
12 months372.92±14.7542,40–0.58±1.05–0.95±0.95
18 months560.37±24.3848,49–0.85±0.91–0.82±0.93
24 months746.61±23.2549,47–0.44±0.99–0.80±0.97
Appendix 1—table 1
FC that significantly changed with age (fNIRS pre-processing without global signal regression).

Results are displayed in terms of estimated betas, standard errors, and p values.

FCFpBaseline(5 months)Betas (SE), p5–8 change Betas (SE), p5–12 change Betas (SE), p5–18 change Betas (SE), p5–24 change Betas (SE), p
HbO2
Left fronto-middle10.4<0.001–0.55 (0.29),<0.0630.05 (0.05), 0.2390.07 (0.04), 0.1250.16 (0.04),<0.0010.26 (0.04),<0.001
Right fronto-middle4.82<0.001–0.24 (0.28), 0.3980.03 (0.05), 0.5350.05 (0.04), 0.2460.08 (0.04), 0.0510.18 (0.04),<0.001
HHb
Frontal interhemispheric3.50.0020.18 (0.33) 0.591–0.07 (0.06), 0.260–0.07 (0.05), 0.214–0.17 (0.05), 0.002–0.17 (0.05), 0.002
Appendix 1—table 2
Results of the regression analyses of the effect of FC on cognitive flexibility in younger and older preschoolers.

Significant positive associations are in green, significant negative associations are in orange, and non-significant (NS) associations are in blue.

Cognitive flexibility in younger preschoolersCognitive flexibility in older preschoolers
Frontal interhemispheric FC5 monthsF(1,38)=4.21, p=0.047, R2=0.102NS
8 monthsNSNS
12 monthsNSNS
18 monthsNSF(1,45)=4.82, p=0.030, R2=0.115
24 monthsNSNS
Left fronto-middle FC5 monthsNSNS
8 monthsNSNS
12 monthsF(1,33)=7.86, p=0.009, R2=0.197NS
18 monthsNSF(1,45)=5.72, p=0.021, R2=0.115
24 monthsNSNS
Right fronto-middle FC5 monthsNSNS
8 monthsNSNS
12 monthsF(1,38)=4.82, p=0.034, R2=0.115NS
18 monthsNSNS
24 monthsNSNS
Right frontal-posterior FC5 monthsNSNS
8 monthsF(1,32)=5.6, p=0.024, R2=0.153NS
12 monthsNSNS
18 moNSNS
24 moNSF(1,49)=4.85, P=0.032, R2=0.092
Appendix 1—table 3
Results of the correlational analyses between changes in growth (ΔWLZ) and cognitive flexibility in younger and older preschoolers.
Cognitive flexibility in younger preschoolers
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthr(54)=−0.169, p=0.214r(51)=−0.167, p=0.232r(51)=−0.249, p=0.072r(53)=−0.198, p=0.147r(51)=−0.151, p=0.280r(41)=−0.241, p=0.102
1 mor(64)=0.031, p=0.805r(63)=0.001, p=0.993r(66)=0.013, p=0.916r(62)=0.084, p=0.511r(48)=−0.021, p=0.885
5 mor(64)=−0.075, p=0.548r(66)=−0.032, p=0.798r(63)=0.008, p=0.947r(50)=−0.091, p=0.520
8 mor(66)=−0.044, p=0.719r(62)=0.067, p=0.600r(50)=−0.041, p=0.775
12 mor(64)=0.097, p=0.439r(51)=−0.012, p=0.932
18 mor(49)=−0.107, p=0.454
24 mo
Cognitive flexibility in older preschoolers
ΔWLZbirth1 month5 months8 months12 months18 months24 months
birthr(55)=0.036, p=0.789r(58)=−0.122, p=0.353r(56)=−0.138, p=0.303r(56)=−0.153, p=0.252r(55)=−0.149, p=0.268r(53)=−0.137, p=0.319
1 monthr(64)=−0.138, p=0.268r(61)=−0.113, p=0.379r(62)=−0.113, p=0.294r(59)=−0.013, p=0.919r(59)=−0.069, p=0.597
5 monthsr(72)=−0.002, p=0.984r(73)=−0.025, p=0.829r(68)=0,064, p=0.601r(70)=0.046, p=0.704
8 monthsr(72)=0.026, p=0.827r(68)=0.098, p=0.421r(69)=0.034, p=0.777
12 monthsr(68)=0.091, p=0.456r(70)=−0.057, p=0.634
18 monthsr(65)=−0.103, p=0.406
24 months
Appendix 1—table 4
Results of the correlational analyses between changes in functional connectivity between 5 and 24 months and cognitive flexibility in younger and older preschoolers.
Cognitive flexibility
ΔFC between 5 and 24 monthsYounger preschoolersOlder preschoolers
Frontal interhemispheric connectivityr(12)=−0.368, p=0.196r(20)=−0.188, p=0.401
Left fronto-middle FCr(10)=0.253, p=0.428r(17)=0.270, p=0.264
Right fronto-middle FCr(12)=0.248, p=0.392r(20)=0.192, p=0.392
Right frontal-posterior FCr(12)=0.188, p=0.520r(20)=0.164, p=0.465
Appendix 1—table 5
Characteristics of included and excluded participants and seconds of data included in the analyses at each age.

WD = withdrawn, D = deceased, MV = missed visit, DD = developmental delay, NIRS not undertaken = the participant was assessed but did not want or could not perform the NIRS assessments, FC not undertaken = the participant was assessed with other NIRS task, but not FC, Fussed out = the participant wore the headband and the FC acquisition had started but the participant showed signs of fussiness soon after the start of the acquisition, MP = missing pictures of the headband placement, EM = missing event markers, TI = technical issues during the NIRS testing session. The proportion of children included in the analysis was computed based on the infants with FC data.

AgeN-%Not testedNIRS not undertakenFC not undertakenInfants with FC dataFussed outExperimental errorsHeadband PlacementToo many channels excludedNot enough data after pre-processingInfants included in the analysesSeconds of data included in the analyses (mean ± SD)Inclusion rate (from the 204 infants recruited)
WDDMVDDMPEMTI
5 monthsN2123107179161133754787382.68±92.7842%
%0.980.490.981.474.903.4387.748.935.391.671.675.582.7926.2548.6
8 monthsN7153181415676012366053372.93±80.6625%
%3.430.492.451.478.826.8676.474.483.8400.6414.73.8438.4633.97
12 monthsN9143171315741021025682372.93±80.6640%
%4.410.491.961.478.336.3776.962.540.6301.276.361.2735.6652.22
18 monthsN8115312516040011643897388.85±81.6347%
%3.920.497.351.475.882.4578.432.5000.62102.523.7560.62
24 monthsN131291341530200644596399.56±79.8647%
%6.370.4914.20.491.471.967501.30003.922.6129.4162.74

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Chiara Bulgarelli
  2. Anna Blasi
  3. Samantha McCann
  4. Bosiljka Milosavljevic
  5. Giulia Ghillia
  6. Ebrima Mbye
  7. Ebou Touray
  8. Tijan Fadera
  9. Lena Acolatse
  10. Sophie E Moore
  11. Sarah Lloyd-Fox
  12. Clare E Elwell
  13. Adam T Eggebrecht
  14. The BRIGHT Study Team
(2026)
Growth in early infancy drives optimal brain functional connectivity which predicts cognitive flexibility in later childhood
eLife 13:RP94194.
https://doi.org/10.7554/eLife.94194.4