We investigated effects of neonatal caffeine treatment on brain structure and longitudinal development in a subset of preterm children from a randomized, placebo-controlled trial. Children had MRI at age 11 years (n=118; 64 received caffeine, 54 received placebo), of whom 43 (21 caffeine, 22 placebo) had MRI at term-equivalent age. Global and regional brain volumes, cortical morphometry and white matter microstructure were similar between treatment groups at age 11 years, as was brain development from term-equivalent age to 11 years of age. Any benefits of caffeine on brain structure weaken over time and were not apparent at age 11 years.
1. Schmidt B, Roberts RS, Davis P, et al. Caffeine therapy for apnea of prematurity. N Engl J Med. 2006;354(20):2112-2121.
2. Schmidt B, Roberts RS, Davis P, et al. Long-term effects of caffeine therapy for apnea of prematurity. N Engl J Med. 2007;357(19):1893-1902.
3. Doyle LW, Schmidt B, Anderson PJ, et al. Reduction in developmental coordination disorder with neonatal caffeine therapy. J Pediatr. 2014;165(2):356-359 e352.
4. Schmidt B, Anderson PJ, Doyle LW, et al. Survival without disability to age 5 years after neonatal caffeine therapy for apnea of prematurity. JAMA. 2012;307(3):275-282.
5. Schmidt B, Roberts RS, Anderson PJ, et al. Academic Performance, Motor Function, and Behavior 11 Years After Neonatal Caffeine Citrate Therapy for Apnea of Prematurity: An 11-Year Follow-up of the CAP Randomized Clinical Trial. JAMA Pediatr. 2017;171(6):564-572.
6. Doyle LW, Cheong J, Hunt RW, et al. Caffeine and brain development in very preterm infants. Ann Neurol. 2010;68(5):734-742.
7. Ashburner J and Friston KJ. Unified segmentation. Neuroimage. 2005;26(3):839-851.
8. Dale AM, Fischl B, and Sereno MI. Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179-194.
9. Fischl B, Sereno MI, and Dale AM. Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. Neuroimage. 1999;9(2):195-207.
10. Smith SM, Jenkinson M, Johansen-Berg H, et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage. 2006;31(4):1487-1505.
11. Diedrichsen J, Balsters JH, Flavell J, et al. A probabilistic MR atlas of the human cerebellum. Neuroimage. 2009;46(1):39-46.
12. Patenaude B, Smith SM, Kennedy DN, et al. A Bayesian model of shape and appearance for subcortical brain segmentation. Neuroimage. 2011;56(3):907-922.
13. Desikan RS, Segonne F, Fischl B, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 2006;31(3):968-980.
14. Adamson C, Beare R, Walterfang M, et al. Software pipeline for midsagittal corpus callosum thickness profile processing : automated segmentation, manual editor, thickness profile generator, group-wise statistical comparison and results display. Neuroinformatics. 2014;12(4):595-614.
15. Adamson CL, Wood AG, Chen J, et al. Thickness profile generation for the corpus callosum using Laplace's equation. Hum Brain Mapp. 2011;32(12):2131-2140.
16. Thompson DK, Inder TE, Faggian N, et al. Characterization of the corpus callosum in very preterm and full-term infants utilizing MRI. Neuroimage. 2011;55(2):479-490.
17. Tournier JD, Calamante F, and Connelly A. Improved probabilistic streamlines tractography by 2nd order integration over fibre orientation distributions. International Society for Magnetic Resonance in Medicine (ISMRM). 2010.
18. Beare RJ, Chen J, Kelly CE, et al. Neonatal Brain Tissue Classification with Morphological Adaptation and Unified Segmentation. Front Neuroinform. 2016;10:12.
Figure 1. Global brain volumes at 11 years of age, contrasted between treatment groups.
Data for caffeine and placebo groups are means (standard deviations). Mean differences and P-values are from separate linear regressions adjusted for age at MRI and sex of the child. The plots are visual representations of the mean differences and confidence intervals (CI). Volume units are cm3.
Figure 2. Brain regions and white matter tracts involved in motor function, at 11 years of age, contrasted between treatment groups.
Data for caffeine and placebo groups are means (standard deviations). Mean differences and P-values are from separate linear regressions adjusted for age at MRI and sex of the child, as well as for total intracranial volume for models with volumes. The plots are visual representations of the mean differences and confidence intervals (CI). Units are for volumes cm3 and for diffusivities x10-3 mm2/s. Data are combined across brain hemispheres as there was negligible evidence of group by hemisphere interactions.
Figure 3. The corpus callosum at 11 years, contrasted between treatment groups.
A. Top: Thickness at 100 points along the corpus callosum; bottom: results of t-tests comparing thickness at each point between groups, with colors indicating p-values.
B. Data from subregions. Data for groups are means and standard deviations. Mean differences are from separate linear regressions adjusted for age and sex of the child, as well as for intracranial volume for models with area and thickness. The plots are visual representations of the mean differences and confidence intervals (CI). Units are for areas cm2, thickness cm, and diffusivities x10-3 mm2/s.
Table 1. Longitudinal changes in brain structures, contrasted between treatment groups.
All data are from separate mixed effects models for each MRI measure.
aData for caffeine and placebo groups are regression coefficients (95% confidence intervals (CIs)), indicating the total change in volume between term and 11 years of age.
bThe interaction coefficients (and p-values) are the difference in the longitudinal volume or area change between treatment groups, adjusted for the age difference between scans, and sex. Units are cm3 for volumes and cm2 for areas.
Table 2. Longitudinal changes in diffusion tensor values from white matter regions of interest, contrasted between treatment groups.
All data are from separate mixed effects models for each region and each diffusion measure.
aData are regression coefficients (95% confidence intervals (CIs)), indicating the total change in diffusion values between term and 11.
bThe interaction coefficients (and p-values) are the differences in the longitudinal diffusion change between treatment groups, adjusted for the age difference between scans, and sex. Results are combined across hemispheres as there was negligible evidence of group by hemisphere interactions.
Units are x10-3 mm2/s for diffusivities.