AUTHOR=Rovnaghi Cynthia R. , Rigdon Joseph , Roué Jean-Michel , Ruiz Monica O. , Carrion Victor G. , Anand Kanwaljeet J. S. 

TITLE=Longitudinal Trajectories of Hair Cortisol: Hypothalamic-Pituitary-Adrenal Axis Dysfunction in Early Childhood

JOURNAL=Frontiers in Pediatrics

VOLUME=9

YEAR=2021

URL=https://www.frontiersin.org/journals/pediatrics/articles/10.3389/fped.2021.740343

DOI=10.3389/fped.2021.740343

ISSN=2296-2360

ABSTRACT=<p>The objective of this study was to examine if longitudinal trajectories of hair cortisol concentrations (HCC) measured at two or three yearly time points can identify 1-3 year old children at risk for altered hypothalamic-pituitary-adrenal (HPA)-axis function due to early life stress (ELS). HCC was measured (<italic>N</italic> = 575) in 265 children using a validated enzyme-linked immunosorbent assay. Hair was sampled in Clinic Visits (CV) centered at years 1, 2, and 3 (<italic>n</italic> = 45); 1 and 2 (<italic>n</italic> = 98); 1 and 3 (<italic>n</italic> = 27); 2 and 3 (<italic>n</italic> = 95). Log-transformed HCC values were partitioned using latent class mixed models (LCMM) to minimize the Bayesian Information Criterion. Multivariable linear mixed effects models for ln-HCC as a function of fixed effects for age in months and random effects for participants (to account for repeated measures) were generated to identify the factors associated with class membership. Children in Class 1 (<italic>n</italic> = 69; 9% Black) evidenced declining ln-HCC across early childhood, whereas Class 2 members (<italic>n</italic> = 196; 43% Black) showed mixed trajectories. LCMM with only Class 2 members revealed Class 2A (<italic>n</italic> = 17, 82% Black) with sustained high ln-HCC and Class 2B (<italic>n</italic> = 179, 40% Blacks) with mixed ln-HCC profiles. Another LCMM limited to only Class 2B members revealed Class 2B1 (<italic>n</italic> = 65, 57% Black) with declining ln-HCC values (at higher ranges than Class 1), and Class 2B2 (<italic>n</italic> = 113, 30% Black) with sustained high ln-HCC values. Class 1 may represent hair cortisol trajectories associated with adaptive HPA-axis profiles, whereas 2A, 2B1, and 2B2 may represent allostatic load with dysregulated profiles of HPA-axis function in response to varying exposures to ELS. Sequential longitudinal hair cortisol measurements revealed the allostatic load associated with ELS and the potential for developing maladaptive or dysregulated HPA-axis function in early childhood.</p>