Article Text
Abstract
Objective To assess the relationship between the 5 min Apgar score and developmental vulnerability at 5 years of age.
Design Population-based retrospective cohort study.
Setting Manitoba, Canada.
Participants All children born between 1999 and 2006 at term gestation, with a documented 5 min Apgar score.
Exposure 5 min Apgar score.
Main outcome measures Childhood development at 5 years of age, expressed as vulnerability (absent vs present) on five domains of the Early Development Instrument: physical health, social competence, emotional maturity, language and cognitive development, and communication skills.
Results Of the 33 883 children in the study, most (82%) had an Apgar score of 9; 1% of children had a score <7 and 5.6% had a score of 10. Children with Apgar scores <10 had higher odds of vulnerability on the physical domain at age 5 years compared with children with a score of 10 (eg, adjusted OR (aOR) for Apgar 9=1.23, 95% CI 1.05 to 1.44). Similarly, children with Apgar scores of <10 were more vulnerable on the emotional domain (eg, aOR for Apgar 9=1.20, 95% CI 1.03 to 1.41). Nevertheless, the Apgar-based prognostic model had a poor sensitivity for physical vulnerability (19%, 95% CI 18% to 20%). Although the Apgar score-based prognostic model had reasonable calibration ability and risk-stratification accuracy for identifying developmentally vulnerable children, classification accuracy was poor.
Conclusions The risk of developmental vulnerability at 5 years of age is inversely associated with the 5 min Apgar score across its entire range, and the score can serve as a population-level indicator of developmental risk.
- Child Psychology
- Comm Child Health
- Epidemiology
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
What is already known on this topic
Low 5 min Apgar scores are associated with adverse short-term and long-term cognitive outcomes and developmental impairment.
A comprehensive view of the relation between Apgar scores and the child's early developmental health has been largely overlooked in the existing literature.
What this study adds
The 5 min Apgar score is inversely associated with a risk of developmental vulnerability in early childhood across its entire range.
There is an increased risk of developmental vulnerability among children with an Apgar score of 9, compared with children with a score of 10.
The Apgar score based prognostic model created reasonably acceptable categories of low and high-risk children for developing an early intervention program.
Introduction
Introduced by Virginia Apgar in 1952, the Apgar score measures the health status of newborn children based on heart rate, respiration, colour, muscle tone and reflex irritability.1 The Apgar score at birth and 1 min after birth was initially developed to assess the immediate condition of the newborn and the potential need for resuscitation. Subsequent studies have demonstrated the 5 min Apgar score's ability to predict longer-term cognitive outcomes such as a reduced IQ, a lower numeracy and literacy score,2 neurological disability and low cognitive function at age 19 years.3 In school-age children, low Apgar scores have been linked to minor motor, language, speech and developmental impairments.4–6
Most follow-up studies examining the long-term prognostic value of the Apgar score on child outcomes have focused on the association between extremely low 5 min Apgar scores, such as those linked to severe birth asphyxia, and distinct neurocognitive profiles such as IQ and academic ability. Virtually, no prior study has examined the developmental correlates of the Apgar score as a continuous variable, across the entire spectrum of recorded scores.2 ,7 ,8 In addition, there are no studies that have assessed the prognostic performance of the 5 min Apgar score for identifying developmental vulnerability among school-age children. A programme of early diagnosis could identify high-risk children who would potentially benefit from interventions designed to mitigate the developmental problems at school-going age. We, therefore, carried out a population-based study to comprehensively examine the relationship between the 5 min Apgar score and developmental vulnerability at 5 years of age.
Methods
Study design
All children born in Manitoba, Canada, between 1 January 1999 and 31 December 2006, with a gestational age of 37 weeks or greater and a documented 5 min Apgar score, as well as a completed Early Development Instrument (EDI) assessment at age 5 years between 2005 and 2011, were included in the study cohort. Data were obtained from the Manitoba Population Health Research Data Repository, which collates information on health and social services utilisation for all residents of Manitoba.9 Two data sources were used to determine the socioeconomic status (SES) of children in the study: the provincial Employment and Income Assistance data (identifying those requiring social assistance or income supplementation) and Census data (providing mean household income by area of residence). Finally, EDI10 data, which provided information on early childhood developmental outcomes and school readiness, were accessed through linkage with the Healthy Child Manitoba Office. All data files used in this study were anonymised, and linkage at the individual level was performed using a scrambled personal health identification number. The reliability and validity of these data sources have been previously documented.11 ,12
Prognostic variables and outcomes
The primary prognostic variable of interest was the routinely collected 5 min Apgar score (hereafter referred to as the ‘Apgar score’) from the hospital record. Apgar scores were analysed in several ways: (a) using conventional categories: Apgar values of 0–6, 7, 8 and 9–10; (b) revised categories: Apgar values of 0–3, 4, 5, 6, 7, 8, 9 and 10 (values from 0 to 3 were grouped together because of small numbers); (c) as a continuous variable and (d) as a continuous variable with a quadratic term to evaluate the potential for a non-linear relationship with developmental outcomes.
Other independent variables included in the prognostic model were infant sex (male vs female), birth weight-for-gestational age, age of the child in years at the time of EDI assessment, gestational age at birth in completed weeks (37, 38, 39, 40, 41 and ≥42), breastfeeding initiation (yes vs no) and SES. Birth weight-for-gestational age was categorised as: small (<10th percentile), appropriate (10th–90th percentile) and large (>90th percentile).13 Each child's family income was derived from the mean household income in the child's residential area (based on postal code) obtained from the 2006 Canadian Census data.14–16 This was complemented by parental receipt of income assistance at any time from the child's birth to the EDI assessment. Family SES was based on income quintile and receipt of income assistance, and was categorised as low (quintile 1 or income assistance recipients),12 medium (quintiles 2 and 3) and high (quintiles 4 and 5).
The primary outcome of interest was childhood development, as measured by the EDI. The EDI was routinely administered biennially in all 37 public school divisions in Manitoba beginning in 2005/06. Teachers—remunerated by their school district—completed the EDI for each child in their kindergarten class (age range 5–7 years) midway through the school year. The EDI consists of 104 binary and Likert-scale items designed to tap five core areas of early childhood development:10 ,17 physical health and well-being; social competence, emotional maturity, language and cognitive development and communication skills and general knowledge. Children were considered vulnerable on a domain if their scores fell below the 10th percentile value18 based on the national EDI scores.19
Statistical analyses
The frequency of each Apgar-score value was calculated within categories of maternal and infant characteristics. Variables chosen as prognostic determinants were selected because of presumed association with developmental vulnerability based on the literature,18 ,20 as well as availability in our data source. Multivariable logistic regression was used to examine the association between each independent variable and vulnerability on each domain of the EDI. Results were expressed as ORs with 95% CIs. Models were evaluated for goodness of fit using the Akaike information criterion (AIC).
In addition to modelling vulnerability on the EDI, we also assessed the performance of the prognostic model for identifying any developmental vulnerability. The model-predicted probability of vulnerability on each EDI domain was categorised as identifying low risk (<10%), medium risk (10%–19%) and high risk (≥20%) of an adverse developmental outcome. The model-predicted probability of vulnerability on any of the EDI domains was categorised as identifying low risk (<30%), medium risk (30%–44%) and high risk (≥45%) of any adverse developmental outcome. These cut-off values were chosen based on the provincial prevalence of vulnerability, since approximately 11% of children are vulnerable on an EDI domain, and 28% of children are vulnerable on one or more domains.21 We assessed performance of the prognostic model in terms of calibration ability, stratification capacity and classification accuracy22 ,23 to identify children at high risk for developmental vulnerability.
A two-sided p value <0.05 was used to define statistical significance. No adjustment was made for multiple comparisons in this exploratory study. The University of Manitoba Health Research Ethics Board sanctioned the study, and the Manitoba Health Information Privacy Committee approved data access. Analyses were performed using SAS V.9.2 (SAS Institute, Cary, North Carolina, USA).
Results
There were 33 883 children (mean age=5.7 years) with a gestational age of ≥37 weeks and complete Apgar and EDI data who were included in the study. The majority (82%) of children had an Apgar score of 9, and 1% had a score <7. Only 0.2% of children had an Apgar score ≤4 (table 1). As expected, low Apgar scores were also more common among males, small-for-gestational age live births, children of mothers who did not initiate breast feeding and those with a low SES.
The prevalence of vulnerability within each of the five EDI domains by Apgar score is shown in figure 1. Overall, the prevalence of vulnerability in one or more domains of the EDI was 28%, with physical and language domains having the highest rates of vulnerability at 12% and 12%, respectively. There was a graded, decreasing trend in the rate of vulnerability on all domains of the EDI with increasing Apgar score. Vulnerability rates for the physical (23%), social (18%) and emotional (16%) domains were highest among children with an Apgar score of 5.
The prognostic model for categorising children into low (<10%), medium (10%–19%) and high risk groupings (≥20%) for vulnerability on the physical domain showed reasonable calibration ability; 7.6% and 25% of low-risk versus high-risk children were found to be vulnerable on the physical domain (table 2). Similar results were obtained with regard to calibration ability on the other EDI domains. The model also successfully stratified children into reasonably sized risk categories, with a large proportion of children classified as low or medium risk, and approximately 5%–12% of children identified as being at high risk for vulnerability in each of the five EDI domains. However, the classification accuracy of the model was poor; the model had a detection rate (sensitivity) between 11% and 26% for vulnerability on the five different domains of the EDI. Similar results were obtained for the prognostic model identifying vulnerability on any of the EDI domains (table 2).
Using the conventional categorisation of Apgar scores (table 3), we observed that compared with those with an Apgar score of 9–10, the odds of vulnerability on the physical domain of the EDI were significantly higher among children with an Apgar score of 0–6 (OR 1.55, 95% CI 1.17 to 2.06). A graded excess in odds of vulnerability was evident with increasing Apgar score though rates of developmental vulnerability among those with Apgar scores of 7 and 8 were not significantly different from those with an Apgar score of 9–10.
Analyses with modelling Apgar scores using the revised categories, however, revealed that the odds of vulnerability on the physical domain increased in a graded fashion across the full range of scores (table 4). Children with an Apgar score of 7, 8 or 9 had significantly higher odds of vulnerability on the physical domain compared with children with an Apgar score of 10 (OR 1.45, 1.24 and 1.23 for Apgar 7, 8 and 9, respectively; table 4). Similarly, children with an Apgar score of 7, 8 or 9 had higher odds of vulnerability on the emotional domain compared with children with an Apgar score of 10 (OR 1.60, 1.24 and 1.20 for Apgar 7, 8 and 9, respectively). Children with scores <4 had 2.71 times higher odds of vulnerability in the language domain compared with children with an Apgar score of 10 (OR 2.71, 95% CI 1.22 to 6.01).
Finally, analyses with Apgar score modelled as a continuous variable (see online supplementary appendix table 1) showed that unit increase in the Apgar score was associated with a significant decrease in vulnerability on the physical (OR 0.91, 95% CI 0.86 to 0.95), social (OR 0.94, 95% CI 0.89 to 0.99) and emotional domains (OR 0.93, 95% CI 0.89 to 0.98). Non-significant decreases in vulnerability were observed on the language and communication domain. Modelling Apgar scores as a continuous variable or as a continuous variable with a quadratic term did not improve model fit (as assessed by the AIC) or change the pattern of the findings.
Discussion
Our population-based study showed an association between the Apgar score and early childhood developmental vulnerability at 5 years of age, with lower 5 min Apgar scores being associated, in a graded continuous manner, with an increased risk of vulnerability on the physical health and well-being, emotional and social domains of the EDI. The Apgar score-based prognostic model created reasonably acceptable categories of low-risk and high-risk children for creating an early intervention programme, but failed to identify a significant proportion of children who were subsequently found to be vulnerable on one or more developmental domains. On the other hand, our study showed an increased risk of vulnerability on the physical and emotional domains even among children with an Apgar score of 9, compared with children with a score of 10.
Although the literature suggests that low Apgar scores are of concern in the short and long term, our results showed that the negative association between Apgar scores and developmental adversity extended in a linear, graded fashion across the entire range of scores. Such a finding—revealing a continuous linear association with a conventionally categorical predictor variable—parallels the identification of SES as a graded, continuous, population-level predictor of morbidity, whereas previously only outright poverty was viewed as important for health and well-being.24 Our findings are consistent with previous studies showing an increased need for special education among children with Apgar scores of 7 or 8 at 5 min after birth compared with children with Apgar score 9 or 10 and the need for extra resources in kindergarten.2 ,4 ,8 ,20
A neurobiological account for why Apgar scores might be significant, continuous-level predictors of developmental outcomes cannot be provided by the data available in this study. Many of the observable physiological components of the Apgar score are mediated by shifts in autonomic nervous system (ANS) activation, and autonomic reactivity to emotional and physical stressors (stressors such as labour and birth) has been linked to chronic disease processes in adults25 and to both developmental26 and health outcomes27 in children. Mothers’ prenatal adversity, moreover, has been associated with trajectories of postnatal development of autonomic reactivity,28 suggesting that even emotional stressors antecedent to birth can influence ANS responses. Thus, the present finding—that the full range of Apgar scores predicts developmental vulnerability at kindergarten—could represent the influence of heightened stress reactivity (indexed via ANS responses to the stress of birth) on long-term developmental endpoints.
Our prognostic model's findings are in line with previous studies that showed that the specificity of the Apgar score as a prognostic tool was far better than its sensitivity.3 ,29 ,30 Further, risk assessment research also shows that biological criteria alone fail to identify 65% of children who, subsequently, experience high rates of poor health and disappointing educational outcomes.31 Although our findings indicate that Apgar scores cannot be effectively used at the individual level to determine prognosis or trigger preventive intervention, our findings could have substantial public health significance as a population's Apgar profile could serve as an indicator of the burden of adverse developmental outcomes in children. Additionally, documentation of differences in 5 min Apgar profiles between populations (eg, a lower frequency of infants with a 5 min Apgar score of 9) could lead to an aetiological search for causes and interventions that improve developmental outcomes in children.
The strengths of our study included the ability to access comprehensive health and education-related databases at the population level. The EDI assessment by teachers avoided reliance on parental or self-report of developmental outcomes. Nevertheless, there may have been some individual differences in teachers’ ability to evaluate developmental outcomes.19 The SES of children was determined by their SES in the first 5 years after birth, although SES at birth would have been preferable for a prognostic model. However, SES as constructed in our study is unlikely to have changed substantially during the study period. Further, our study was restricted to the comparatively healthy subset of all term live births, as children with severe chronic illnesses may not have enrolled in kindergarten or may have enrolled in special needs schools. This was not a serious limitation as our objective was to examine the role of the Apgar score as a marker of childhood development in normal children.
In summary, our study showed that the risk of developmental vulnerability at 5 years of age was inversely associated with the 5 min Apgar score and that this relation extended across the entire range of scores. Thus, even a minor degree of physiological dysfunction soon after birth, reflected in the Apgar score, may indicate a slightly higher risk for developmental vulnerability in later childhood. Although the low sensitivity of the prognostic model based on the 5 min Apgar score limits the clinical utility of this model for early intervention, the Apgar score might well serve as a population-level indicator of developmental risk.
Acknowledgments
Appreciation is extended to the Manitoba Health Policy Centre (MCHP) staff for providing the data under project (HIPC #2012/2013-07) and specifically to Charles Burchill at MCHP for his continuous support and guidance throughout the linkage and data analyses process.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
- Data supplement 1 - Online supplement
Footnotes
Contributors NR: conceptualised and designed the study, carried out all the data analysis, drafted the initial manuscript and contributed to the acquisition of data. WTB: conceptualisation and design of the study, advised in the interpretation of data and critically reviewed and revised the manuscript. MB, DJ and HT assisted in interpretation of the data and critically reviewed the manuscript. RAM contributed to the acquisition of the data, interpretation of the data and critically reviewed the manuscript. KSJ supervised the conceptualisation and design of the study, advised in the interpretation of data and critically reviewed and revised the manuscript. All authors approved the final manuscript as submitted.
Funding This work was supported by the Canadian Institute of Health Research (MOP-119393, PI Helen Tremlett). NR is funded by the Multiple Sclerosis Society of Canada (PhD Research Studentships). She has received travel grants from the Canadian Institutes of Health Research (CIHR) and the endMS Research and Training Network/Multiple Sclerosis Society of Canada. WTB is funded by the Canadian Institute for Advanced Research; the BC Leadership Chair in Child Development also supported this work. MB is the recipient of the Manitoba Centre for Health Policy Population-Based Child Health Research award funded by the Government of Manitoba. KSJ is the recipient of a Chair in maternal, fetal and infant health services research from the CIHR. His work is also supported by the Child and Family Research Institute. HT is funded by the Multiple Sclerosis Society of Canada (Don Paty Career Development Award), Michael Smith Foundation for Health Research, and is the Canada Research Chair for Neuroepidemiology and Multiple Sclerosis. She has received research support from the US National Multiple Sclerosis Society, CIHR, and UK MS Trust; speaker honoraria and/or travel expenses to attend conferences from the Consortium of MS Centres, US National MS Society, the University of British Columbia Multiple Sclerosis Research Program, Bayer Pharmaceutical (speaker, 2010, honoraria declined), Teva Pharmaceuticals (speaker 2011), ECTRIMS (2011, 2012, 2013), UK MS Trust (2011), the Chesapeake Health Education Program, US Veterans Affairs (2012, honorarium declined), Novartis (2012), the Consortium of MS Centers (2013) and Biogen (honorarium declined). Unless otherwise stated, all speaker honoraria are either donated to an MS charity or to an unrestricted grant for use by her research group. RAM receives research funding from CIHR, Multiple Sclerosis Society of Canada, Multiple Sclerosis Scientific Research Foundation, National Multiple Sclerosis Society, Rx & D Health Research Foundation, and has conducted clinical trials funded by Sanofi-Aventis.
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.