Skip to main content

Associations between the social environment and early childhood developmental outcomes of Puerto Rican children with prenatal Zika virus exposure: a cross-sectional study

Abstract

Background

Prenatal exposure to the Zika virus can lead to microcephaly and adverse developmental outcomes, even in children without evident birth defects. The social environment plays a crucial role in infant health and developmental trajectories, especially during periods of heightened brain plasticity. The study aimed to assess socioenvironmental factors as predictors of developmental outcomes of 36-month-old children exposed to Zika virus prenatally.

Study design

This cross-sectional study included 53 mothers and 55 children enrolled in the Pediatric Outcomes of Prenatal Zika Exposure cohort study in Puerto Rico. The study performs follow-up developmental assessments of children born to mothers with confirmed and probable Zika virus infection during pregnancy. Mothers completed socioenvironmental questionnaires (e.g., Perceived Neighborhood Scale and US Household Food Insecurity Survey). Children’s developmental outcomes were assessed with the Bayley Scales of Infant and Toddler Development: Third Edition, the Ages and Stages Questionnaires: Third Edition, the Ages and Stages Questionnaire-Socioemotional: Second Edition, and the Child Adjustment and Parent Efficacy Scale.

Results

Linear regression models, adjusting for a child’s sex and age and maternal education, revealed that early life exposure to food insecurity and maternal pregnancy stressors were significantly associated with poorer developmental outcomes in Zika virus-exposed children at 36 months of age. Maternal resilience representation of adaptive ability was associated with the preservation of adequate developmental outcomes in children.

Conclusions

Pregnancy and early childhood are critical life periods for ensuring optimal brain development in children. While the mechanisms in the interaction of children with their environment are complex, the risk and protective factors identified in the study are modifiable through public policy and preventive initiatives. Implementation of comprehensive strategies that improve access to social support programs, educational and nutritional interventions, and mental health services during pregnancy and early childhood can enhance the developmental potential of vulnerable children.

Peer Review reports

Introduction

The environment that surrounds and interacts with children offers opportunities for growth, development, and play that have long-lasting effects on their overall health and well-being. Experts attribute a key role to a stable and responsible caregiving environment in driving normative brain and behavioral development [1]. In fact, relational health in stable, nurturing, and safe environments has been confirmed to mitigate the effects of chronic hardship, which can ensue from disadvantaged social determinants of health, poverty, maternal stress, and food insecurity [2]. In addition, a child’s neighborhood environment is recognized for its critical direct and indirect influences on child health during sensitive life periods from fetal life to early childhood [3]. Protective factors and risks of the early environment interact and contribute to child health, development, and future health disparities through complex mechanisms, including epigenetic, neurohormonal, and relational mechanisms, that merit further study to inform healthcare and public health policy [4].

Congenital infections disrupt neuronal processes through direct infection of developing brain cells, exposure to inflammatory cytokines during maternal infection, and direct damage to the placenta, leading to poor perfusion and malfunction that alters fetal development [5]. Zika virus infection during pregnancy has been associated with congenital Zika syndrome microcephaly and severe brain defects in infants, as well as neurodevelopmental concerns that develop over time in those without detectable birth defects [6,7,8,9]. However, brain plasticity during early childhood provides a window of opportunity for stimulating brain development through interactions with the environment and adaptations based on learning and experiences. These dynamic interactions are influenced by socioenvironmental factors, which have not been thoroughly analyzed as risk or protective factors for developmental outcomes in children with prenatal Zika exposure.

Since the 2016-17 Zika epidemic in Puerto Rico (PR), the Pediatric Outcomes of Prenatal Zika Exposure (POPZE) study systematically followed a cohort of Hispanic children with exposure. Close developmental monitoring and assessments have revealed that at least one-third of the children at 36 months of age have developmental delays in at least one domain (cognitive, language, or motor), with a reduced prevalence of delays noted in some domains over time [6]. Data on developmental delays in Puerto Rican children 0 to 5 years are limited. However, the prevalence of delay in the study children is notably higher than the 2018–2019 US prevalence for emotional, developmental, and behavioral problems of 10.8% among children aged 0–5 years [10]. The occurrence of delay in the language domain decreased from 50.9% at 24 months to 31.9% at 36 months of age. This prevalence is similar to the combined prevalence of 29.7% reported in a recent systematic review by Marbán-Castro et al. [11]. among normocephalic children with prenatal Zika virus exposure. This manuscript provides updated results of developmental outcomes from 55 children at 36 months of age with prenatal Zika virus exposure and explores socioenvironmental factors as predictors of such outcomes.

Methods

Participants in this cross-sectional study include children and mothers enrolled in the POPZE cohort study in Puerto Rico between May 2017 and December 2019 with follow-up developmental assessments at 36 months. The study included professional assessments and parental perspectives on the children’s development and collected socioenvironmental risk and protective factors with an impact on childhood development at the 36-month visit (from September 2019 to March 2021) (see Supplementary Table 1 for detailed information on the timing of assessment and response rate for each instrument and Supplementary Table 2 for psychometric properties of the tools).

The original phase of the study, POPZE I, aimed to characterize the full spectrum of structural and functional abnormalities in children born to mothers infected with ZIKV infection during pregnancy. From May 2017 to June 2019, the team enrolled 114 exposed children and followed 97 of them. With a new funding source starting in August 2019, the team expanded its focus in POPZE II to assess the interaction of socioenvironmental aspects with the exposed children’s development. This phase enrolled the mothers as active participants, with 58 out of 97 (59.8%) re-enrolled for the expanded follow-up visits. The re-consent period took place during the follow-up visits beginning in August 2019. Fifty-five out of 58 (94.8%) participants enrolled in POPZE II completed the assessment of the 36-month visit. Three children were excluded because they did not attend the 36-month scheduled appointment. Participants’ enrollment, inclusion, and exclusion criteria were previously described in a manuscript by Alvarado-Domenech et al. [6]. Briefly, children were born at two academic hospitals in the southern region of Puerto Rico, and their mothers had confirmed (polymerase chain reaction (PCR) positive) or probable (immunoglobulin M (IgM) positive) Zika virus infection. The Ponce Medical School Foundation Institutional Review Board reviewed and approved this study. All parents provided informed written consent for both their own and their child’s participation in the study.

Developmental outcomes

Experienced licensed clinical psychologists administered the Bayley Scales of Infant and Toddler Development, Third Edition (BSID-III) [12] in Spanish to assess the children’s neurodevelopmental functioning based on cognitive, language, and motor domains at 36 months of age. Percentile scores describe the child’s performance in each domain; higher scores imply higher functioning. Parents answered the Ages and Stages Questionnaire, Third Edition (ASQ-3) [13], providing perceptions on their child’s development in communication, gross motor, fine motor, problem solving, and personal-social skills. Higher scores imply better execution on the assessed domain. Parents also completed the Ages and Stages Questionnaire: Social-Emotional, Second Edition (ASQ: SE-2) [14], which assesses children’s responses indicating risk for social or emotional difficulties. Higher scores imply higher social or emotional difficulties. These developmental tests have been widely used to detect early developmental delay or risk of delay in other Zika virus cohorts [6, 7, 15,16,17], allowing comparison with other populations.

The Child Adjustment and Parent Efficacy Scale (CAPES) [18, 19] assesses children’s behavioral and emotional difficulties through the Behavioral and Emotional Adjustment scales, respectively. These scales are combined in the intensity scale, which describes the parent’s perception on the children’s psychological difficulties. The CAPES also assesses parent’s perception of their self-efficacy in managing their child’s behavioral and emotional problems through the Parent Self-Efficacy scale. Some items are reverse scored. Higher scores indicate a higher level of the construct measured.

Socioenvironmental factors

Neighborhood characteristics include the Area Deprivation Index (ADI) [20, 21] and the Perceived Neighborhood Scale (PNS) [22]. The team used the Neighborhood Atlas platform to obtain the participants’ home ADI, an index that ranks the neighborhood based on US Census-identified socioeconomic disadvantage [20, 21]. Participants’ physical addresses were entered on Puerto Rico’s map to obtain the national percentile rank for each participant. A higher percentile score denotes a more disadvantaged neighborhood. Also, the team used a cutoff percentile of 80 to classify families living in higher disadvantage areas. The PNS assesses the parental perception of a neighborhood within four dimensions: social embeddedness, sense of community, neighborhood satisfaction, and perceived crime [22]. Higher scores imply a higher degree of the dimension measured.

Home and family characteristics include the Home Observation Measurement of the Environment – Short Form (HOME-SF) [23], the McMaster Family Assessment Device-Short Form (FAD-SF) [24, 25] and the United States Department of Agriculture (USDA) – Household Food Security Survey Module (HFSSM) [26], all completed by the mothers. The HOME-SF assesses the quality of the cognitive stimulation (e.g., “How often do you read stories to your child?”) and emotional support (e.g., “Does your child see his/her father or father-figure on a daily basis?”) in the child’s home environment [23]. Higher scores indicate a better home environment. The FAD-SF assesses family functioning (e.g., “We can express feelings to each other”) [24]. Higher scores indicate worse functioning [27]. The USDA – HFSSM assesses household food security. Parents reported their access to food of sufficient quantity and quality for the last 12 months. Based on the USDA – HFSSM interpretation scale [26], a score of 1 is considered marginal food insecurity, suggesting a family is vulnerable to food insecurity. Higher scores indicate higher food insecurity.

Maternal characteristics include the following self-report instruments: Zika Virus-Related Prenatal Stress Scale (ZIKV-PSS; developed by researchers; α = 0.89), Beck Depression Inventory (BDI-II) [28], Maternal Resilience Scale (ERESMA) [29], and Social Provisions Scale (SPS) [30]. The ZIKV-PSS assesses stress during pregnancy, including Zika virus-related stressors (e.g., “Receive a Zika diagnosis,” “Not having family support”). To reduce recall bias, mothers were asked to make an effort to place themselves at the stage of pregnancy. The BDI-II explores depression symptomatology (e.g., “Sadness”). ERESMA assesses the resilience of mothers of exceptional children (e.g., “I try to make sure that my child is happy”) [29]. The SPS evaluates perceived social support (e.g., “There are people I can count on in an emergency”) [30]. For each instrument, higher scores indicate a greater degree of the construct assessed.

Statistical analyses

Descriptive statistics were calculated for sociodemographic and socioenvironmental factors and developmental outcomes of the participants in the study. For continuous variables, the mean (SD) was calculated, while percentages were used for categorical variables. Socioenvironmental factors and developmental outcome scores were analyzed as continuous variables. Maternal education was categorized into two groups: (1) high school or lower and (2) higher than high school, representing the highest level of education achieved. Linear regression models were employed, adjusting for child sex, age, and maternal education to examine the associations between each socioenvironmental characteristic and each developmental outcome; P values of < 0.05 denoted statistical significance. Analyses were performed with SAS version 9.4 (SAS Institute).

Results

The sociodemographic and clinical characteristics of the mother and child are presented in Table 1. Fifty-three mothers of 55 children, including two sets of twins, completed the study questionnaires accounting for a 100% response rate. The maternal mean age was 30 years. Of the mothers, 49.1% (26) were unemployed, and 77.4% (41) had an educational level of more than high school, of which 31.7% (13) had technical certificates and 34.2% (14) had associate degrees. Most women were married or cohabitating (n = 40; 75.7%), had public health insurance (n = 41; 77.4%), and had an annual household income <$15,000 (n = 38; 71.7%). Zika virus infection during pregnancy was documented in the first (n = 22; 41.5%), second (n = 18; 33.9%), and third (n = 13; 24.5%) trimesters. Maternal diagnostic tests indicated: PCR positive (n = 16; 30.2%), IgM positive (n = 27; 50.9%), or both (n = 10; 18.9%) assays results. Of the children enrolled in this study phase, none (0%) had microcephaly (head circumference z-score < -2), 2 of 53 (3.8%) had abnormal hearing screening, and 11 of 17 (64.7%) had retinal images abnormalities at birth. The children’s mean age at the 36-month assessment was 36.6 months, and 56.4% (31) were female.

Table 1 Maternal and child characteristics and socioenvironmental profile

Socioenvironmental factors

The ADI neighborhood measure revealed a mean percentile rank of 89.4 (SD, 11.1), indicating that most participants (90.6%) live in highly disadvantaged areas (Table 1). However, maternal reports on the PNS revealed mean scores that suggest they perceive an average social embeddedness (25.0; SD, 9.3), a high sense of community (26.4; SD, 6.6) and satisfaction (36.9; SD, 6.6), and a low perceived crime (16.1; SD, 8.1). Regarding the home and family characteristics, parental mean scores on the HOME cognitive stimulation (2.9; SD, 1.0) and emotional support (2.9; SD, 0.98) scales were average (Table 1). Additionally, the mean score (17.3; SD, 6.5) on the FAD-SF scale revealed adequate family functioning. The household food security mean score (1.4; SD, 2.9) suggested that about one-third (21/53; 39.6%) of these families experienced food insecurity vulnerability, a rate disproportionally higher than the prevalence of 13.3% and 6% among the households of US children with and without special healthcare needs, respectively [31]. Regarding maternal characteristics, maternal perceived stress during pregnancy varied (range, 6–59), with a moderate mean score (31.1; SD, 12.2) on the prenatal stress scale (Table 1). Depression symptomatology on the BDI-II mean score was low (8.8; SD, 9.3). Mean scores on the ERESMA scale (211.2; SD, 13.1) suggested that these women manifest high levels of resilience. Additionally, mean scores on the SPS (89.9; SD, 11.8) reflected that they perceive excellent sources of social support (e.g., having a person to rely on when facing difficulties and counting on individuals who recognize and appreciate their skills and abilities).

Developmental outcomes

At 36 months of age, the mean percentile of children’s neurodevelopment scores in the BSID-III indicate average performance in all domains: cognitive (34.5; SD, 17.3), language (29.2; SD, 23.4), and motor (30.8; SD, 21.4) (Table 2). Mean scores on the ASQ-3 were average for communication (44.8; SD, 18.2), gross motor (55.4; SD, 6.6), and fine motor (40.3; SD, 18.1) skills domains. However, parental reports showed below average scores on problem solving (37.2; SD, 15.1) and personal-social (43.2; SD, 14.9) domains, both in the “monitoring zone”, indicating that the children could benefit from further follow-up developmental assessments.

Table 2 Child neurodevelopmental assessment scores and prevalence of developmental delay/risk according to BSID-III, ASQ-3, ASQ: SE-2, and CAPES

Parental reports on the ASQ: SE-2 displayed average scores (54.3; SD, 32.6), suggesting adequate socioemotional development. This was supported by low mean scores on the CAPES intensity (18.9; SD, 11.0), emotional (1.6; SD, 2.1), and behavior (17.3; SD, 9.4) scales. Additionally, mean scores on the CAPES self-efficacy scale (168.9; SD, 27.3) suggested a high level of perceived parental ability to manage their children’s emotional and behavioral problems. On further analyses, no significant differences were found based on trimester of exposure.

Associations between children’s developmental outcomes and socioenvironmental factors

Linear regression models, adjusted by age, child’s sex, and maternal education, were used to explore associations between BSID-III, ASQ-3, ASQ: SE-2, and CAPES outcomes and socioenvironmental factors (Tables 3, 4 and 5). Multiple socioenvironmental factors showed significant associations with child developmental outcomes, but the most salient findings indicated that food insecurity and maternal prenatal stress are linked to poorer developmental outcomes. However, maternal resilience stood out positively, revealing a link with better child outcomes. The following paragraphs illustrate the results of the linear regression analyses by each child’s outcome.

Table 3 Associations between socio-environmental characteristics and BSID-III scores at 36 months (B, SE)*
Table 4 Associations between socio-environmental characteristics and ASQ-3 / ASQ: SE-2 scores at 36 months (B, SE)*

Children who experience higher scores in food insecurity (B = -1.62, SE = 0.77, p = 0.04) and reduced maternal resilience (B = 0.48, SE = 0.19, p = 0.01) had significantly lower scores in the BSID-III cognitive domain (Table 3). Those with higher maternal prenatal stress scores (B = -0.79, SE = 0.23, p = 0.01) and higher ADI percentile ranks, thus living in highly disadvantaged neighborhoods (B = -0.51, SE = 0.26, p = 0.05), had significantly lower scores in the BSID-III language domain. Children with higher maternal prenatal stress scores (B = -0.59, SE = 0.24, p = 0.02) had significantly lower BSID-III motor domain scores.

Consistently, those children with higher scores in food insecurity (B = -1.73, SE = 0.80, p = 0.04) (Table 4) had significantly lower scores in the ASQ-3 communication domain. Those with lower scores in maternal resilience (B = 0.34, SE = 0.17, p = 0.05) had significantly lower scores in the ASQ-3 problem-solving domain. Additionally, children with higher scores in neighborhood perceived crime (B = -0.64, SE = 0.25, p = 0.02), food insecurity (B = -1.63, SE = 0.54, p = 0.005), and maternal prenatal stress (B = -0.36, SE = 0.15, p = 0.02) and lower maternal resilience (B = 0.31, SE = 0.14, p = 0.04) showed significantly lower scores in the ASQ-3 personal-social domain. Children with higher scores for food insecurity (B = 4.28, SE = 1.73, p = 0.02) and maternal prenatal stress (B = 0.73, SE = 0.05, p = 0.04) and lower scores for maternal resilience (B = -1.02, SE = 0.32, p = 0.02) had significantly higher scores for social-emotional difficulties (ASQ: SE-2).

Table 5 Associations between socio-environmental characteristics and CAPES scores at 36 months (B, SE)*

The CAPES intensity scale combines the emotional adjustment and behavior scales, thus reflecting emotional and behavioral problems. Children with lower scores on HOME emotional support (B = -2.99, SE = 1.39, p = 0.04) and maternal resilience (B = -0.36, SE = 0.11, p = 0.002) and higher scores on food insecurity (B = 1.37, SE = 0.46, p = 0.005), maternal prenatal stress (B = 0.36, SE = 0.12, p = 0.003) and depression symptomatology (B = 0.48, SE = 0.15, p = 0.002) had significantly higher scores on the intensity scale (Table 5). Those with higher scores in social embeddedness (B = 1.30, SE = 0.38, p = 0.001), sense of community (B = 1.62, SE = 0.52, p = 0.003), neighborhood satisfaction (B = 1.27, SE = 0.60, p = 0.04), social support (B = 1.16, SE = 0.30, p = 0.0003), and maternal resilience (B = 1.21, SE = 0.28, p < 0.0001) had significantly higher maternal self-efficacy scores. Conversely, higher scores in family functioning (indicating worst functioning) (B = -2.32, SE = 0.58, p < 0.0001) and maternal prenatal stress (B = -0.70, SE = 0.33, p = 0.04) and depression symptomatology (B = -2.08, SE = 0.32, p < 0.0001) were significantly linked with lower maternal self-efficacy scores.

Discussion

This study explored socioenvironmental variables as risk and protective factors in the development of 55 Hispanic/Latino children from Puerto Rico whose neurodevelopment is biologically vulnerable from prenatal Zika exposure. The findings presented significant associations between early life exposure to food insecurity and maternal pregnancy stressors in the context of social vulnerabilities with developmental risks in these children, highlighting two environmental impediments to the national goal of eliminating disparities and promoting healthy development and well-being across all life stages for all [32]. However, maternal resilience stood out as a remarkable adaptive ability linked to better children’s developmental outcomes.

Social determinants of health significantly impact overall health, well-being, and quality of life [33]. However, beyond addressing each of the vulnerabilities individually, we must consider the impact of the interrelated experiences of material hardships (food insecurity, housing, medical, etc.) [34] that families face. Nevertheless, participants’ characteristics highlight socioeconomic vulnerabilities strongly associated in other studies with the risk of gross motor [35], fine motor [36], communication [37], cognitive [38, 39], and language [40] delays and emotional problems [38, 41,42,43]. Participants’ vulnerabilities include high unemployment rates, limited household income, public health insurance, and living in disadvantaged areas. These occur in a critical context of years of economic recession in Puerto Rico. Currently, 55% of Puerto Rican families with children have a household income below the poverty level [44], which is disproportionately higher than the rate in the continental U.S. (14.1%) [45], denoting an increased risk for children not attaining developmental potential due to extreme poverty [46, 47]. In addition, there are inherent characteristics of children with developmental disabilities that put them at risk for disadvantage. For example, children with mental, behavioral, and emotional difficulties more often live in the lowest income categories [48].

Food insecurity may be influenced by income and employment status [49] and can limit a child’s intake of essential nutrients, with implications for poor brain development [50]. Despite the compensation of the Nutrition Assistance for Puerto Rico program, the inflation and elevated cost of living represents challenges to Puerto Rican families and jeopardizes access to adequate food for family consumption. Similar to studies linking child malnutrition and lower cognitive scores [51], children in this study who experienced food insecurity had lower performance in BSID cognitive skills. A 24-month longitudinal study by Milner et al. [52]. determined that food insecurity’s timing, intensity, and duration are associated with lower ASQ communication, social, and gross motor scores. While the food insecurity assessment was from the past 12 months, we also found associations with lower ASQ-3 social and communication scores but not with motor scores. We identified other associations with increased behavioral and emotional problems, consistent with other studies of preschool children [53, 54]. Evidence shows that children with special healthcare needs are twice as likely to be food insecure as children without healthcare special needs [31], possibly due to the higher expenditures on services required by this population. Participants in this study showed a greater vulnerability to food insecurity than children with special healthcare needs, with a prevalence of 41.8% versus 13.3% [31].

The findings highlight food insecurity as a noteworthy risk factor for early childhood development within this vulnerable population, even after adjusting for important confounders. Indeed, Hobbs et al. [51] explored the associations between food insecurity and behavioral and cognitive outcomes in 5-year-old children from disadvantaged families. The analysis involved examining these associations across percentiles of food insecurity distribution (10th, 25th, 50th, 75th, and 90th), while controlling for potential confounders such as maternal depression, parenting stress, and material hardship. They found that both externalizing and internalizing behavioral problems were linked to food insecurity across all percentiles. Despite differences in the methodology used to explore associations, this finding is consistent with the association between food insecurity and behavioral and emotional difficulties in POPZE II children, as measured by the CAPES intensity scale and ASQ: SE-2. Regarding cognitive skills, Hobbs et al. [51] reported variability in the evaluated outcomes and the food insecurity percentiles where associations were identified. Food insecurity was associated with lower receptive vocabulary scores in children beyond the 50th percentile. Conversely, low scores in letter-word identification skills were associated with food insecurity only in children at the 10th percentile. The variability in cognitive outcomes might be due to differences in the measures used to assess cognition. Despite the variability, these findings align with POPZE II findings supporting a negative association between food insecurity and cognitive skills. Difficulties identified could potentially have repercussions for psychosocial problems later in adolescence [55] and adulthood [56] in already developmentally vulnerable children.

Cohorts describing birth and follow-up consequences of prenatal Zika have not extensively explored socioenvironmental factors in association to developmental outcomes to allow for comparisons with the Puerto Rican experience. Nevertheless, poor socioenvironmental conditions are known to increase the health risks of vulnerable populations during epidemics. Inadequate home conditions can facilitate Zika virus transmission, and malnutrition and co-infections in the context of socioeconomic disadvantage can affect immune status and the response to infections [57]. In this sense, the study setting serves to highlight the socioenvironmental experiences of many families and children in different countries impacted by the Zika epidemic.

In this study, mothers experienced prenatal Zika virus-related biopsychosocial stressors in the context of social determinants of health vulnerabilities and disadvantages at the individual, family, and neighborhood levels (e.g., low education attainment and household income, food insecurity, and living in disadvantaged areas), factors that have been strongly associated with developmental and health risks in children [58, 59]. Although prenatal stress data were explored retrospectively, introducing recall bias and possible minimization of symptomatology, the findings support previous literature connecting prenatal stress and suboptimal language [60, 61], motor [62], personal-social [63], social-emotional [64, 65], and behavioral [61, 66] abilities in children. The results can also be linked to accumulated evidence that repeated adverse experiences cause “toxic” stress-induced dysregulated stress responses in the mother that can be transmitted intergenerationally to the offspring with long-term consequences in health, development, and well-being [67, 68].

From a positive perspective, despite biopsychosocial vulnerabilities and their stressors, mothers report the ability to adapt and overcome. Maternal resilience was a protective factor in these children’s cognitive, problem solving, personal-social, social-emotional, and behavioral skills. The buffering effect of maternal resilience on child development in the face of adversity is supported by improved family functioning skills and self-efficacy in handling child behavior. Of note, resilience emerges from the dynamic interplay of multiple processes occurring across and between systems [69]. Despite living in disadvantaged areas, families perceive a high sense of community and neighborhood satisfaction. Individuals tend to be more resilient when they have family, social, and community support [70]. In this context, the interaction of children with these systems is critical in learning and developing resilience skills. The role of families and caregivers and the impact of Zika virus exposure on children and families is considerable and should not be overlooked [71]. A positive environment, modeling from significant others, nurturing caregivers, and responsive caregiving promotes positive adjustment and better developmental outcomes [3, 72, 73].

Study limitations include the cross-sectional design and limited sample size, which reduces the generalization of the results. Additionally, data on maternal prenatal stress were collected retrospectively, introducing a recall bias. The lack of a control group precludes the study team from determining whether identified risks and protective factors are shared among groups or are specific to children exposed to Zika virus prenatally. Furthermore, the attrition from the original cohort, mostly due to natural disasters, might have impacted the capacity to retain children with the more severe outcomes in the cohort. We also acknowledge that the psychometric assessment of the newly developed Zika Virus-Related Prenatal Stress Scale requires further validation analyses.

Despite study limitations, to our knowledge, this study is among the first to contribute to assessing multiple socioenvironmental factors as predictors of developmental outcomes of children exposed to Zika virus prenatally. Through a comprehensive assessment implementing standardized procedures with valid and reliable instruments, the findings also contribute to a greater understanding of the factors that may buffer the adverse impact of the Zika virus on the exposed children’s brain development. Immediate family resilience and community supportive influences can be plausible explanations for improved developmental outcomes in some domains observed over time [6] and may add to the understanding of the absence of statistically significant differences in developmental outcomes between exposed and unexposed children reported by prenatal Zika studies of 18- to 30-month-old children who were asymptomatic at birth [16, 74]. The study will publish results of school readiness assessments of exposed and unexposed children from the same community, thereby addressing whether the exposure status determines significant differences in cognitive functioning and other developmental outcomes in the context of relevant socioenvironmental factors.

Conclusions

Pregnancy and the first years of a child’s life are vital for optimal brain development, as they are periods of great sensitivity and vulnerability for establishing the structural and functional abilities of the brain. Most of the risk and protective factors identified in this study are modifiable through public policy that supports the implementation of comprehensive strategies such as increasing access to social support programs, enhancing educational and nutritional interventions, and providing mental health support to pregnant women and families. Addressing social determinants of health requires transdisciplinary collaboration from policymakers, community partners, and healthcare providers to ensure maternal and child well-being and better health outcomes. Timely provisions of resources and interventions can reduce inequities and enhance the development potential of children and families.

Data availability

The datasets used and analyzed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

ABR:

Automated Auditory Brainstem

ADI:

Area deprivation index

ASQ-3:

Ages and Stages Questionnaire, Third Edition

ASQ: SE-2:

Ages and Stages Questionnaire: Social-Emotional, Second Edition

BDI-II:

Beck Depression Inventory – Second Edition

BSID-III:

Bayley Scales of Infant and Toddler Development, Third Edition

CAPES:

Child Adjustment and Parent Efficacy Scale

HC:

Head circumference

HOME-SF:

Home Observation Measurement of the Environment – Short Form

HFSSM:

Household Food Security Survey Module

IgM:

Immunoglobulin M

PCR:

Polymerase chain reaction

RetCam:

Retinal imaging with RetCam 3 Visualization System

SPS:

Social Provisions Scale

USDA:

United States Department of Agriculture

ZIKV-PSS:

Zika Virus-Related Prenatal Stress Scale

References

  1. Bick J, Nelson CA. Early adverse experiences and the developing brain. Neuropsychopharmacol. 2016;41(1):177–96.

    Article  CAS  Google Scholar 

  2. Garner A, Yogman M, Committee on Psychosocial Aspects of Child and Family Health, Section on Developmental, and Behavioral Pediatrics, Council on Early Childhood. Preventing childhood toxic stress: partnering with families and communities to promote relational health. Pediatrics. 2021;148(2):e2021052582.

    Article  PubMed  Google Scholar 

  3. National Scientific Council on the Developing Child. Place matters: The environment we create shapes the foundations of healthy development: Working paper No. 16. 2023. https://developingchild.harvard.edu/place-matters-the-environment-we-create-shapes-the-foundations-of-healthy-development/. Accessed 8 August 2023.

  4. Shonkoff JP, Boyce WT, McEwen BS. Neuroscience, molecular biology, and the childhood roots of health disparities: building a new framework for health promotion and disease prevention. JAMA. 2009;301(21):2252–9.

    Article  CAS  PubMed  Google Scholar 

  5. Fortin O, Mulkey SB. Neurodevelopmental outcomes in congenital and perinatal infections. Curr Opin Infect Dis. 2023;36(5):405–13.

    Article  CAS  PubMed  Google Scholar 

  6. Alvarado-Domenech LI, Rivera-Amill V, Appleton AA, Rosario-Villafañe V, Repollet-Carrer I, Borges-Rodríguez M, et al. Early childhood neurodevelopmental outcomes in children with prenatal Zika virus exposure: a cohort study in Puerto Rico. J Pediatr. 2022;247:38–e455.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Faiçal AV, de Oliveira JC, Oliveira JVV, de Almeida BL, Agra IA, Alcantara LCJ, et al. Neurodevelopmental delay in normocephalic children with in utero exposure to Zika virus. BMJ Paediatr Open. 2019;3(1):e000486.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Mulkey SB, Arroyave-Wessel M, Peyton C, Bulas DI, Fourzali Y, et al. Neurodevelopmental abnormalities in children with in utero Zika virus exposure without congenital Zika syndrome. JAMA Pediatr. 2020;174(3):269–76.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Nielsen-Saines K, Brasil P, Kerin T, Vasconcelos Z, Gabaglia CR, et al. Delayed childhood neurodevelopment and neurosensory alterations in the second year of life in a prospective cohort of ZIKV-exposed children. Nat Med. 2019;25(8):1213–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 2018–2019. National Survey of Children’s Health (NSCH) data query. https://www.childhealthdata.org/browse/survey/results?q=7965&r=1&g=812. Accessed 8 August 2023.

  11. Marbán-Castro E, Goncé A, Fumadó V, Romero-Acevedo L, Bardají A. Zika virus infection in pregnant women and their children: a review. Eur J Obstet Gynecol Reprod Biol. 2021;265:162–8.

    Article  PubMed  Google Scholar 

  12. Bayley N. Bayley Scales of Infant and Toddler Development: Administration Manual. 3rd ed. United States of America: The Psychological Corporation; 2006.

    Google Scholar 

  13. Squires J, Twombley E, Bricker D, Potter L. Ages and stages questionnaires (ASQ): a parent-completed child monitoring system. 3rd ed. Baltimore, MD: Brookes; 2009.

    Google Scholar 

  14. Squires J, Briker D, Twombly E. The ages and stages questionnaires: Social-Emotional (ASQ:SE) user’s guide. Baltimore: Paul H; 2002.

    Google Scholar 

  15. Lopes Moreira ME, Nielsen-Saines K, Brasil P, Kerin T, Damasceno L, Pone M, et al. Neurodevelopment in infants exposed to Zika virus in utero. N Engl J Med. 2018;379(24):2377–9.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gerzson LR, de Almeida CS, Silva JHD, Feitosa MMA, de Oliveira LN, Schuler-Faccini L. Neurodevelopment of nonmicrocephalic children, after 18 months of life, exposed prenatally to Zika virus. J Child Neurol. 2020;35(4).

  17. Alger J, Cafferata ML, López R, Wiggins LD, Callejas A, Castillo M, et al. Neurodevelopmental assessment of normocephalic children born to Zika virus exposed and unexposed pregnant people. Pediatr Res. 2024;95(2):566–72.

    Article  CAS  PubMed  Google Scholar 

  18. Morawska A, Sanders M, Haslam D, Filus A, Fletcher R. Child Adjustment and parent efficacy scale (CAPES): development and initial validation of a parent report measure. Aust Psychol. 2014;44:241–52.

    Article  Google Scholar 

  19. Seijo D, Tomé D, Sanmarco J, Morawska A, Fariña F. Spanish adaptation and validation of the Child Adjustment and parent efficacy scale. Sustainability. 2021;13(9):4647.

    Article  Google Scholar 

  20. Kind AJH, Buckingham WR. Making Neighborhood-Disadvantage Metrics Accessible - The Neighborhood Atlas. N Engl J Med. 2018;378(26):2456–8.

    Article  PubMed  PubMed Central  Google Scholar 

  21. University of Wisconsin School of Medicine and Public Health. 2015 Area Deprivation Index v3.2. https://www.neighborhoodatlas.medicine.wisc.edu/. Accessed 8 August 2023.

  22. Martinez ML, Black M, Starr RH. Factorial structure of the Perceived Neighborhood Scale (PNS): a test of longitudinal invariance. J Commun Psychol. 2002;30(1):23–43.

    Article  Google Scholar 

  23. Bradley RH, Caldwell BM. Home observation for measurement of the environment: a revision of the preschool scale. Am J Ment Defic. 1979;84(3):235–44.

    CAS  PubMed  Google Scholar 

  24. Epstein N, Bishop D, Levin S. The McMaster Model of Family Functioning. J Marital Fam Ther. 1978;4:19–31.

    Article  Google Scholar 

  25. Barroilhet S, Cano-Prous A, Cervera-Enguix S, Forjaz MJ, Guillén-Grima F. A Spanish version of the Family Assessment device. Soc Psychiatry Psychiatr Epidemiol. 2009;44(12):1051–65.

    Article  PubMed  Google Scholar 

  26. Bickel G, Nord M, Price C, Hamilton W, Cook J. Guide to measuring household food security. Alexandria, VA, USA: US Department of Agriculture Food and Nutrition Service; 2000.

    Google Scholar 

  27. Miller IW, Ryan CE, Keitner GI, Bishop DS, Epstein NB. The McMaster Approach to families: theory, assessment, treatment and research. J Family Ther. 2000;22(2):168–89.

    Article  Google Scholar 

  28. Beck A, Steer R, Brown G. BDI-II Manual. San Antonio: The Psychological Corporation; 1996.

    Google Scholar 

  29. Roque MP, Acle G, García M. Escala de resiliencia materna: Un estudio de validación en una muestra de madres con niños especiales. Rev Iberoam Diagn Ev. 2009;1(27):107–32.

    Google Scholar 

  30. Cutrona C, Russel D. The provisions of social relationships and adaptation to stress. Adv Personal Relationships. 1987;1(1):37–67.

    Google Scholar 

  31. Balistreri KS. Food insufficiency and children with special healthcare needs. Public Health. 2019;167:55–61.

    Article  CAS  PubMed  Google Scholar 

  32. Office of Disease Prevention and Health Promotion. Healthy People 2030. In: U.S. Department of Health and Human Services. https://health.gov/healthypeople. Accessed 8 August 2023.

  33. Poulain T, Vogel M, Kiess W. Review on the role of socioeconomic status in child health and development. Curr Opin Pediatr. 2020;32(2):308–14.

    Article  PubMed  Google Scholar 

  34. Schenck-Fontaine A, Ryan RM. Poverty, material hardship, and children’s outcomes: A nuanced understanding of material hardship in childhood. Child (Basel). 2022,9(7).

  35. Veldman SL, Jones RA, Chandler P, Robinson LE, Okely AD. Prevalence and risk factors of gross motor delay in pre-schoolers. J Paediatr Child Health. 2020;56(4):571–6.

    Article  PubMed  Google Scholar 

  36. Camden C, Héguy L, Casoli M, Roy M, Rivard L, Berbari J, Ppop Research Team CM. Preschoolers’ Developmental Profiles and School-Readiness in a Low-Income Canadian City: A Cross-Sectional Survey. Int J Environ Res Public Health. 2020;17(7).

  37. Hammer CS, Morgan P, Farkas G, Hillemeier M, Bitetti D, Maczuga S. Late talkers: a Population-based study of risk factors and School Readiness consequences. J Speech Lang Hear Res. 2017;60(3):607–26.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Dea C, Gauvin L, Fournier M, Goldfeld S. Does Place Matter? An International comparison of early childhood development outcomes between the metropolitan areas of Melbourne, Australia and Montreal, Canada. Int J Environ Res Public Health. 2019;16(16):2915.

    Article  PubMed  PubMed Central  Google Scholar 

  39. González L, Cortés-Sancho R, Murcia M, Ballester F, Rebagliato M, Rodríguez-Bernal CL. The role of parental social class, education and unemployment on child cognitive development. Gac Sanit. 2020;34(1):51–60.

    Article  PubMed  Google Scholar 

  40. Schneider T, Linberg T. Development of socio-economic gaps in children’s language skills in Germany. Longit Life Course Stud. 2021;13(1):87–120.

    Article  PubMed  Google Scholar 

  41. Goldfeld S, O’Connor M, Chong S, Gray S, O’Connor E, Woolfenden S, et al. The impact of multidimensional disadvantage over childhood on developmental outcomes in Australia. Int J Epidemiol. 2018;47(5):1485–96.

    Article  PubMed  Google Scholar 

  42. Teymoori A, Côté SM, Jones BL, Nagin DS, Boivin M, Vitaro F, Orri M, Tremblay RE. Risk factors associated with boys’ and girls’ developmental trajectories of physical aggression from early childhood through early adolescence. JAMA Netw Open. 2018;1(8):e186364.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Herrmann J, Vogel M, Pietzner D, Kroll E, Wagner O, Schwarz S, et al. Factors associated with the emotional health of children: high family income as a protective factor. Eur Child Adolesc Psychiatry. 2018;27(3):319–28.

    Article  PubMed  Google Scholar 

  44. Puerto Rico. Deficient on child and youth well-being by maintaining a 55% poverty level. https://www.juventudpr.org/blog/deficiente-puerto-rico-en-el-bienestar-de-la-ninez-y-la-juventud-al-mantener-un-nivel-de-pobreza-de-55. Accessed 10 August 2023.

  45. Food Stamps / Supplemental Nutrition Assistance Program (SNAP). https://data.census.gov/table?q=unemployment&t=Income+and+Poverty&tid=ACSST1Y2021.S1702. Accessed 10 August 2023.

  46. Black MM, Walker SP, Fernald LCH, Andersen CT, DiGirolamo AM, Lu C, et al. Early childhood development coming of age: science through the life course. Lancet. 2017;389(10064):77–90.

    Article  PubMed  Google Scholar 

  47. Gil JD, Ewerling F, Ferreira LZ, Barros AJ. Early childhood suspected developmental delay in 63 low- and middle-income countries: large within- and between-country inequalities documented using national health surveys. J Glob Health. 2020;10(1):010427.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Cree RA, Bitsko RH, Robinson LR, Holbrook JR, Danielson ML, Smith C, et al. Health care, family, and community factors associated with mental, behavioral, and developmental disorders and poverty among children aged 2–8 years - United States, 2016. MMWR Morb Mortal Wkly Rep. 2018;67(50):1377–83.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Food Security in the U.S. https://www.ers.usda.gov/topics/food-nutrition-assistance/food-security-in-the-u-s/. Accessed 10 August 2023.

  50. Jun S, Zeh MJ, Eicher-Miller HA, Bailey RL. Children’s dietary quality and micronutrient adequacy by food security in the household and among household children. Nutrients. 2019;11(5):965.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Hobbs S, King C. The unequal impact of food insecurity on cognitive and behavioral outcomes among 5-year-old urban children. J Nutr Educ Behav. 2018;50(7):687–94.

    Article  PubMed  Google Scholar 

  52. Milner EM, Fiorella KJ, Mattah BJ, Bukusi E, Fernald LCH. Timing, intensity, and duration of household food insecurity are associated with early childhood development in Kenya. Matern Child Nutr. 2018;14(2):e12543.

    Article  PubMed  Google Scholar 

  53. Dean G, Vitolins MZ, Skelton JA, Ip EH, Lucas CB, Brown CL. The association of food insecurity with mental health in preschool-aged children and their parents. Pediatr Res. 2023;94(1):290–5.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Thomas MMC, Miller DP, Morrissey TW. Food insecurity and child health. Pediatrics. 2019;144(4).

  55. Paquin V, Muckle G, Bolanis D, Courtemanche Y, Castellanos-Ryan N, Boivin M, et al. Longitudinal trajectories of food insecurity in childhood and their associations with mental health and functioning in adolescence. JAMA Netw Open. 2021;4(12):e2140085.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Cook JT, Black M, Chilton M, Cutts D, Ettinger de Cuba S, Heeren TC, et al. Are food insecurity’s health impacts underestimated in the U.S. population? Marginal food security also predicts adverse health outcomes in young U.S. children and mothers. Adv Nutr. 2013;4(1):51–61.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Barbeito-Andrés J, Schuler-Faccini L, Garcez PP. Why is congenital Zika syndrome asymmetrically distributed among human populations? PLoS Biol. 2018;16(8):e2006592.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Entringer S, Buss C, Wadhwa PD. Prenatal stress, development, health and disease risk: a psychobiological perspective-2015 Curt Richter Award Paper. Psychoneuroendocrino. 2015;62:366–75.

    Article  Google Scholar 

  59. Bruner CACE. Place, race, and poverty: building hope for children. Acad Pediatr. 2017;17(7S):S123–9.

    Article  PubMed  Google Scholar 

  60. Laplante DP, Brunet A, Schmitz N, Ciampi A, King S. Project Ice Storm: prenatal maternal stress affects cognitive and linguistic functioning in 5 1/2-year-old children. J Am Acad Child Adolesc Psychiatry. 2008;47(9):1063–72.

    Article  PubMed  Google Scholar 

  61. D’Souza S, Crawford CN, Buckley J, Underwood L, Peterson ER, Bird A, et al. Antenatal determinants of early childhood talking delay and behavioural difficulties. Infant Behav Dev. 2019;57:101388.

    Article  PubMed  Google Scholar 

  62. Haselbeck C, Niederberger U, Kulle A, Wache K, Brauner E, Gutermuth M, et al. Prenatal maternal distress seems to be associated with the infant’s temperament and motor development: an explorative study. J Neural Transm (Vienna). 2017;124(7):881–90.

    Article  CAS  PubMed  Google Scholar 

  63. Wei Q, Zou J, Ma X, Xiao X, Zhang Y, Shi H. Prospective associations between various prenatal exposures to maternal psychological stress and neurodevelopment in children within 24 months after birth. J Affect Disord. 2023;327:101–10.

    Article  PubMed  Google Scholar 

  64. Wu Y, Espinosa KM, Barnett SD, Kapse A, Quistorff JL, Lopez C, et al. Association of elevated maternal psychological distress, altered fetal brain, and offspring cognitive and social-emotional outcomes at 18 months. JAMA Netw Open. 2022;5(4):e229244.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Duguay G, Garon-Bissonnette J, Lemieux R, Dubois-Comtois K, Mayrand K, Berthelot N. Socioemotional development in infants of pregnant women during the COVID-19 pandemic: the role of prenatal and postnatal maternal distress. Child Adolesc Psychiatry Ment Health. 2022;16(1):28.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Tuovinen S, Lahti-Pulkkinen M, Girchenko P, Heinonen K, Lahti J, Reynolds RM, et al. Maternal antenatal stress and mental and behavioral disorders in their children. J Affect Disord. 2021;278:57–65.

    Article  PubMed  Google Scholar 

  67. Ahmad SI, Shih EW, LeWinn KZ, Rivera L, Graff JC, Mason WA, et al. Intergenerational transmission of effects of women’s stressors during pregnancy: child psychopathology and the protective role of parenting. Front Psychiatry. 2022;13:838535.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Foss S, Petty CR, Howell C, Mendonca J, Bosse A, Waber DP, et al. Associations among maternal lifetime trauma, psychological symptoms in pregnancy, and infant stress reactivity and regulation. Dev Psychopathol. 2023;35(4):1714–31.

    Article  PubMed  Google Scholar 

  69. Masten AS. Resilience in developing systems: progress and promise as the fourth wave rises. Dev Psychopathol. 2007;19(3):921–30.

    Article  PubMed  Google Scholar 

  70. Southwick SM, Bonanno GA, Masten AS, Panter-Brick C, Yehuda R. Resilience definitions, theory, and challenges: interdisciplinary perspectives. Eur J Psychotraumatol. 2014;1;5.

  71. Kapogiannis BG, Chakhtoura N, Hazra R, Spong CY. Bridging knowledge gaps to Understand how Zika Virus exposure and infection affect Child Development. JAMA Pediatr. 2017;171(5):478–85.

    Article  PubMed  Google Scholar 

  72. Bethell CD, Gombojav N, Whitaker RC. Family resilience and connection promote flourishing among US children, even amid adversity. Health Aff (Millwood). 2019;38(5):729–37.

    Article  PubMed  Google Scholar 

  73. Salinas J, O’Donnell A, Kojis DJ, Pase MP, DeCarli C, Rentz DM, et al. Association of social support with brain volume and cognition. JAMA Netw Open. 2021;4(8):e2121122.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Blackmon K, Evans R, Fernandes M, Landon B, Noel T, Macpherson C, et al. Neurodevelopment in normocephalic children with and without prenatal Zika virus exposure. Arch Dis Child. 2022;107(3):244–50.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The team thanks the children and families for participating in this study.

Funding

Funding for this study was provided by the National Institutes of Health (NIH), National Institute of Minority Health and Health Disparities (NIMHD) (Grant U54MD007579). The NIH-NIMHD had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; or the decision to submit the manuscript for publication. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Author information

Authors and Affiliations

Authors

Contributions

MR-R and LIA-D conceptualized, designed, implemented, supervised, acquired funding, conducted data curation and analysis, and wrote the original draft. AAA analyzed the data and contributed to the original draft. VR-V, IR-C, and MB-R collected and performed data curation. LD-P and KM-R contributed to data curation. MG, PV-G, CM-O, and OO-R collected clinical data. VR-A contributed statistical resources. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Mary Rodríguez-Rabassa.

Ethics declarations

Ethics approval and consent to participate

The Ponce Medical School Foundation Institutional Review Board reviewed and approved this study (IRB approval number 1903008631; May 20, 2019). All parents provided informed written consent for both their own and their child’s participation in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rodríguez-Rabassa, M., Appleton, A.A., Rosario-Villafañe, V. et al. Associations between the social environment and early childhood developmental outcomes of Puerto Rican children with prenatal Zika virus exposure: a cross-sectional study. BMC Pediatr 24, 342 (2024). https://doi.org/10.1186/s12887-024-04806-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12887-024-04806-y

Keywords