Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Family nurture intervention in the NICU increases autonomic regulation in mothers and children at 4-5 years of age: Follow-up results from a randomized controlled trial

  • Martha G. Welch ,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing

    mgw13@cumc.columbia.edu

    Affiliations Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America, Department of Psychiatry, Columbia University Medical Center, New York, NY, United States of America

  • Joseph L. Barone,

    Roles Data curation, Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America

  • Stephen W. Porges,

    Roles Formal analysis, Methodology, Writing – review & editing

    Affiliations Department of Psychiatry, University of North Carolina at Chapel Hill, Chapel Hill, NC, United States of America, Traumatic Stress Research Consortium, Kinsey Institute, Indiana University, Bloomington, IN, United States of America

  • Amie A. Hane,

    Roles Conceptualization, Formal analysis, Methodology, Writing – review & editing

    Affiliations Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America, Department of Psychology, Williams College, Williamstown, MA, United States of America

  • Katie Y. Kwon,

    Roles Investigation, Methodology, Project administration, Writing – review & editing

    Affiliation Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America

  • Robert J. Ludwig,

    Roles Conceptualization, Investigation, Methodology, Project administration, Writing – review & editing

    Affiliation Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America

  • Raymond I. Stark,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America

  • Amanda L. Surman,

    Roles Data curation, Formal analysis, Writing – original draft

    Affiliation Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America

  • Jacek Kolacz,

    Roles Formal analysis, Writing – review & editing

    Affiliation Traumatic Stress Research Consortium, Kinsey Institute, Indiana University, Bloomington, IN, United States of America

  • Michael M. Myers

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Pediatrics, Columbia University Medical Center, New York, NY, United States of America, Department of Psychiatry, Columbia University Medical Center, New York, NY, United States of America

Abstract

Background

Maturation of multiple neurobehavioral systems, including autonomic regulation, is altered by preterm birth. The purpose of this study was to determine the long-term effects of Family Nurture Intervention (FNI) in the NICU on autonomic regulation of preterm infants and their mothers.

Method

A subset of infants and mothers (48% of infants, 51% of mothers) randomly assigned to either standard are (SC), or SC plus the FNI in the NICU in a prior RCT (ClincalTrials.gov; NCT01439269) returned for follow-up assessments when the children were 4 to 5 years corrected age (CA). ECGs were collected for 10 minutes in mothers and their children while children were in their mothers’ laps. Heart rate, standard deviation for heart rate, respiratory sinus arrhythmia (RSA)–an index of parasympathetic regulation, and a measure of vagal efficiency were quantified.

Results

Both children and mothers in the FNI group had significantly greater levels of RSA compared to the SC group (child: mean difference = 0.60, 95% CI 0.17 to 1.03, p = 0.008; mother: mean difference = 0.64, 95% CI 0.07 to 1.21, p = 0.031). In addition, RSA increased more rapidly in FNI children between infancy and the 4 to 5-year follow-up time point (SC = +3.11±0.16 loge msec2, +3.67±0.19 loge msec2 for FNI, p<0.05). These results show that the rate of increase in RSA from infancy to childhood is more rapid in FNI subjects.

Conclusion

Although these preliminary follow-up results are based on approximately half of subjects originally enrolled in the RCT, they suggest that FNI-NICU led to healthier autonomic regulation in both mother and child, when measured during a brief face-to-face socioemotional interaction. A Pavlovian autonomic co-conditioning mechanism may underly these findings that can be exploited therapeutically.

Introduction

Approximately 10% of infants are born before 36 weeks in the US [1]. Despite remarkable medical progress, these preterm infants are at risk for multiple adverse neurobehavioral outcomes [2, 3]. Paralleling these issues, preterm infants also have altered function of the autonomic nervous system (ANS). Compared with babies born at term age, preterm infants display less vagal regulation of heart rate (HR), as well as alterations in sympathetic modulation of blood pressure [4, 5]. As young adults, preterm infants have lower levels of beat-to-beat heart rate variability (HRV) and slower HR recovery from exercise, all potential markers of unhealthy autonomic function [6].

Maternal factors, epigenetics, and the intrauterine environment, as well as the early postnatal experience, all shape the developing central and autonomic nervous systems. Dysfunction of these systems is implicated in the subsequent development of neuropsychiatric disorders and poorer physical health throughout the lifespan [7, 8]. There is growing evidence that the physical and emotional separation of infants and mothers inherent to care in the Neonatal Intensive Care Unit (NICU) contributes to these long-term adverse autonomic outcomes [911]. As a consequence, some interventions and care practices involving mother/infant physical and emotional contact have been implemented in an effort to improve neurobehavioral outcomes in preterm infants.

One widely practiced NICU intervention is skin-to-skin contact (SSC) between mother and infant. A 2016 Cochrane review of early SSC in the hospital supported the use of early SSC between healthy mothers and infants to promote breastfeeding, but concluded that more rigorously designed studies with larger sample sizes were necessary [12]. However, the review did include two studies of preterm infants that reported infant autonomic regulation outcomes. Infants who received SSC immediately after birth were more likely to achieve cardiorespiratory stabilization in the first hours after birth [13, 14].

Few interventions involving mother/infant interactions in the NICU have rigorously tested autonomic function over time. In one study of preterm infants, it was found that SSC led to more rapid maturation of cardiac vagal tone between 32 and 37 weeks [15]. A follow-up assessment at 10 years of age found significantly increased amplitude of RSA and attenuated stress responses. In addition, mothers’ anxiety in the SSC group was reduced when their preterm infants were 6 to 10 years of age [16].

To address these issues, the Nurture Science Program at Columbia University Medical Center (nurturescienceprogram.org) conducted a randomized controlled trial (RCT) of another NICU-based intervention, Family Nurture Intervention (FNI) comparing standard care alone (SC) with SC plus FNI (ClinicalTrials.gov, NCT01439269). FNI-NICU is aimed at facilitating emotional connection between infants and mothers during their traumatic early separation in the NICU [17]. Specially trained Nurture Specialists, who in this study were former NICU nurses, initiated FNI, on average, 1 week after birth while the infant was still in the incubator. The central focus of FNI is to help mothers and their infants emotionally connect via repeated engagement in “calming sessions.” During these ~1 hour-long sessions, mothers were encouraged to express a full range of feelings and speak emotionally with their infants in their native language. Initially, these sessions involved making contact with their infants through the ports of the incubator, using firm and sustained touch. Later, these sessions occurred outside the incubator with the mother holding their infant while in a comfortable chair. Skin-to-skin holding was encouraged for these sessions but some mothers preferred clothed holding. In addition to calming sessions prior to leaving the unit, mothers were given two small cotton cloths, one to wear in her bra and the other to place under the head of her infant, which they were encouraged to exchange daily. On average, FNI mothers engaged in these activities for 6 hours/week [18].

While the intervention encourages some techniques often practiced in the NICU, such as SSC, FNI does not rely solely on SSC and can be initiated while the infants are still confined to the incubator. Indeed, the first deep emotional exchange between infant and mother often took place while the infant was confined to the incubator via comfort touch and maternal emotional expression.

We previously reported many significant results from the FNI-NICU trial [11]. FNI was found to be feasible and safe [19]. FNI mothers had increased maternal sensitivity during normal care-giving activities toward the end of the NICU stay [20], and decreased symptoms of maternal anxiety and depression at 4 months infant corrected age (CA) [21]. Near term, FNI infants showed significantly increased high frequency electroencephalogram (EEG) power [22], altered EEG-coherence [23], and advanced maturation of brain activity [24]. At 18 months of age, FNI infants had improved cognitive and language scores on the Bayley, fewer attention problems on the Child Behavioral Check List, and decreased risk for socioemotional problems as assessed by the Modified Checklist for Autism in Toddlers (MCHAT) [25].

Recently, we reported that FNI infants demonstrated more rapid increases in parasympathetic function from ~35 weeks to ~41 weeks postmenstrual age [26]. In this follow-up study of a the FNI-NICU cohort, we tested the hypothesis that this unique intervention for mothers and infants in the NICU would lead to greater levels of parasympathetic regulation at 4 to 5 years CA. Theoretical considerations and clinical implications are discussed.

Method

Subjects

The data for these analyses were obtained from a subset of the 115 mothers and 150 infants originally recruited into a randomized controlled trial of the Family Nurture Intervention (FNI = 78) or standard NICU care (SC = 72) (Fig 1). This was a single center, parallel-group RCT conducted in the level IV NICU at Morgan Stanley Children’s Hospital of New York. Subjects for this trial were preterm infants born between 26 and 34 weeks of gestation. The trial enrolled 72 infants into the SC group and 78 infants into the intervention group [27].

Of the 56 mothers enrolled into the SC group, 40 gave birth to a single baby and 16 gave birth to twins. Thus, the 72 SC infants were comprised of 40 singletons and 32 were twins. Of the 59 mothers enrolled into the FNI group, 39 gave birth to a single baby and 20 gave birth to twins. For one mother in the FNI group who gave birth to twins, one of twins died prior to enrollment thus, for the FNI group, the 78 infants were comprised of 39 singletons and 39 twins.

The Columbia University Medical Center IRB approved all recruitment, consent, and study procedures on September 24, 2008. Mothers provided written informed consent for themselves and their infants. This study was initiated in January of 2009 with the last enrollment occurring in June of 2012. The 4–5 year follow-up assessments occurred during the period from December of 2013 through January of 2018.

On September 23, 2011 approximately half way through recruitment, the trial was registered (NCT01439269). At the time this trial began, we were unaware of the International Committee of Medical Journal Editor’s (ICMJE) position on this matter. Indeed, the position of the National Institutes of Health, which guided our IRB and Clinical Trials Office, was not clear at the time we initiated the trial. As soon as we became aware of the ICMJE’s position we registered the trial and did not alter the protocol after registration other than to expand the inclusion age range of the infants from 28–32 weeks gestational age (GA) to 26–34 weeks GA, which was reported as required. Thus, although trial registration occurred after initiation of enrollment, the protocol was not changed subsequent to the beginning of recruitment, and outcomes were not affected. The authors confirm that all ongoing and related trials for this intervention are registered.

Family Nurture Intervention and standard NICU care

Study staff met with mothers assigned to the FNI group at the earliest time point after delivery to the intervention procedures. FNI aims to facilitate an emotional connection between mother and infant via calming sessions that include comfort touch, and communication of affect both in the incubator and during clothed contact and SSC. Mothers and infants in the SC group received medical and nursing care that is typical for our NICU; routine contact with medical staff, bedside nursing staff, and access to a psychologist and social worker. Additional study staff who were not responsible for facilitating the intervention met with mothers in both groups once weekly throughout NICU hospitalization to administer questionnaires. See [27] for a comprehensive report of the FNI protocol.

4–5 year assessments

Study subjects were contacted when children were within their 4th and 5th year CA and asked to return to the Nurture Science Laboratory for a follow-up visit. It was not possible to blind either the Nurture Specialists or mothers during the intervention phase of the research. However, the research staff involved in conducting the 4–5 year follow-ups were blinded to group assignment.

ECG acquisition and lap test

Participants were connected to recording equipment shortly after arriving for the visit. Mothers and children were then brought into an empty room that contained only a large, comfortable chair. Mothers were asked to place their child on their lap facing one another. They were instructed to remain facing each other for 10 minutes and asked to “interact as they would normally.” During this period, child and mother ECGs were digitized at 1000 samples/sec using AcqKnowledge software and the MP150 base system with BioNomadix 2Ch Wireless ECG Amplifier from BIOPAC Systems, Inc. (Goleta, CA).

Physiology processing

R-waves from the digitized ECG waveforms were marked using LabChart (ADInstruments, Dunedin, New Zealand). This was done in a 2-step process, first using Labchart’s automatic R-wave detection module, followed by visual inspection of the file to correct misplaced marks. From this software, an inter-beat-interval file was generated. These files were further edited to correct or remove segments of movement artifact using CardioEdit (Brain-Body Center, University of Illinois at Chicago, Chicago, IL). Finally, HR and HRV parameters were analyzed using CardioBatch (Brain-Body Center, University of Illinois at Chicago, Chicago, IL).

Parameters derived by CardioBatch from the corrected R-to-R intervals included HR (bpm), heart period (R-R intervals in msec), the standard deviation of HR (bpm) and respiratory sinus arrhythmia (RSA, in loge msec2). The variance of HR occurring within the frequencies of spontaneous breathing was used as a metric for RSA. This procedure was based on the pioneering research in the electrophysiology of the vagus conducted by Hering [28]. Hering discovered that only the specific vagal fibers that were cardioinhibitory had a respiratory rhythm. Thus, quantifying the variance in HR within the frequencies of spontaneous breathing provides the functional measure of Hering’s observation. These metrics were computed both in 10 second epochs (for estimates of Veff) and in 30 second epochs. The average of the 30 epochs for each child was computed for each of the above parameters. Using the multiple values derived from the 10-second epochs, slopes of the relationship between heart period (i.e., epoch mean of the R-to-R intervals) and RSA for each child and mother were computed as a marker of the efficiency of vagal regulation of dynamic HR [29]. That is, a greater change in heart period per change in RSA, is taken as an index of greater efficiency of vagal modulation of HR.

Data analyses

Group differences for all variables were determined using analyses of variance. Covariates included factors unrelated to the intervention that might have influenced the child’s and/or mother’s physiology: gestational age at birth, birth weight, age in months at the time of assessment (adjusted for weeks premature), sex, and twin status (Y/N) as covariates.

A second set of analyses were conducted to ensure that significant group differences were not biased due to inclusion of multiple children from a single mother (i.e. twins). Included in these analyses was data from singletons plus data from only one subject from each set of twins. Results presented in this article were not the primary outcome of the above described RCT, rather they were additional measures obtained 4 to 5 years after the primary outcome was obtained. Due to excessive movement artifact not all children provide usable data during this period. For children with more than one follow-up study we included results from the study that was closest to 60 months CA and that had usable ECG data. These numbers reflect inclusion of twins when both twins were assessed. Because twins may not represent independent samples and therefore inflate the degrees of freedom for the analyses, we repeated these analyses when only one twin was included.

Results

Of the original 150 infants from 115 families, 93 children (44 SC, 49 FNI) and their mothers returned for a ~4-year, ~5-year, or both 4-year and 5-year follow-up studies. Data from 71 children (36 SC, 35 FNI) were included in these analyses. Among these were 29 twins (15 SC, 14 FNI). There were 11 sets of twins (6 SC, 5 FNI) for which data from both twins were available and 7 twins (3 SC, 4 FNI) for which data from only 1 twin was available. In one series of analyses of children, all 71 subjects were included.

For the second set of analyses excluding second twins, 30 SC and 30 FNI child values were included. There were 55 maternal assessments; however, 10 of these mothers were assessed twice, once with each twin. For the 11th set of twins, both twins were assessed, but the mother’s data was usable for only one of these tests. Thus, in total 45 mothers were included in these analyses. As for children, two sets of analyses were performed using maternal data. The first included all 55 assessments with 10 mothers included twice, once with each twin. A second set of analyses included data recorded from 45 mothers; 31 mothers (15 SC, 16 FNI) while they were with their singleton child, 11 mothers while they were with twin A, and 3 mothers while with twin B.

The current analyses are based on ECG data obtained during a 10-minute “lap test” (see below). In total, data from 71 children (36 SC and 35 FNI) were included. ECGs of the mothers of these children were also recorded during the lap test. However, due to excessive movement and equipment failure, data from only 45 mothers (23 SC, 22 FNI) were usable. Note that some mothers with twins were assessed twice, once with each twin thus providing a total of 55 maternal recordings (29 SC, 26 FNI). Table 1 provides a summary of infant and maternal demographic characteristics for the subjects included in these analyses vs those who were not. These variables differed only with regard to income. Reported earnings were higher among included FNI families than among those excluded. Among SC families, those included reported a greater proportion of income at the two extremes, while for both study groups, families that were excluded were more likely to have opted not to disclose income information. Table 2 summarizes perinatal clinical characteristics of infants who did and did not return for the follow-up visit(s) at 4 to 5 years corrected age.

thumbnail
Table 1. Demographic variables of mothers, fathers and families who were included (n = 60) or not included in the current analyses (n = 55).

https://doi.org/10.1371/journal.pone.0236930.t001

thumbnail
Table 2. Perinatal clinical characteristics of infants who did (Included N = 71) and did not (n = 57) return for the follow-up visit(s) at 4 to 5 years corrected age, or who returned but did not provide usable data during the lap test (n = 22) (Excluded = 79).

https://doi.org/10.1371/journal.pone.0236930.t002

Effects of FNI on child ANS variables

Fig 2 presents the means (±SE) of RSA for children at the 5-year follow-up during the lap-test. The analyses of covariance showed that FNI infants had significantly higher RSA regardless of which strategy was used (All records: mean difference = 0.60, 95% CI 0.17 to 1.03, F(1,64) = 7.53, p = 0.008; Singleton + Twin “A”: mean difference = 0.50, 95% CI 0.04 to 0.96, F(1,53) = 4.65, p = 0.036).

thumbnail
Fig 2. Means (±SE) of RSA for children at the 5-year follow-up during the lap-test.

The set of bars on the left show results for all children and the set of bars on the right show these results when including only singletons and twin A.

https://doi.org/10.1371/journal.pone.0236930.g002

Table 3 provides the means and SEs for the three other variables analyzed: HR, SD of HR, Veff. Note that the values for HR and SD-HR are based on averages across multiple 30 sec epochs. Values for Veff are derived from 10-second epochs and reflect the slopes of the relationship between heart period (i.e., epoch mean of the R-to-R intervals) and RSA, one value for each subject. Only overall variability in HR (SD-HR) showed a significant group effect with FNI children having increased levels of SD-HR. However, when RSA is included as a covariate in the analysis of SD-HR, group differences are no longer significant (adjusted means and SEs; SC, 7.50±0.23 vs FNI, 7.69±0.23, p = 0.56). This supports the conclusion that it is the high frequency, i.e. RSA, component of total HRV that is affected by the intervention.

thumbnail
Table 3. Means (SE) for heart rate (HR, bpm), standard deviation of HR (SD-HR, bpm), and vagal efficiency (Veff, bpm/in unit of RSA) for standard care (SC) and Family Nurture Intervention (FNI) children assessed at 5 years of age (60.0 months corrected age ± 6.0 months SD).

https://doi.org/10.1371/journal.pone.0236930.t003

Effects of FNI on maternal ANS variables

Fig 3 presents the means (SE) of RSA for mothers at their child’s 4 to 5 year follow-up visit during the lap-test. Results from the analyses of covariance showed that mothers of FNI infants had significantly higher RSA for both sets of analyses (All records: mean difference = 0.64, 95% CI 0.07 to 1.21, F(1,48 = 4.94, p = 0.031; Singleton + Twin “A”: mean difference = 0.69, 95% CI 0.05 to 1.33, F(1,38) = 4.62, p = 0.038).

thumbnail
Fig 3. Means (±SE) of RSA for mothers at their child’s 4 to 5 year follow-up during the lap-test.

The set of bars on the left show results for all children. These studies of mothers reflect multiple measures from some mothers tested with both of their twins. The set of bars on the right show results when including only singletons and twin A.

https://doi.org/10.1371/journal.pone.0236930.g003

Table 4 provides the means and SEs for the three other maternal variables analyzed: HR, SD of HR, Veff. As was the case for child variables, only the SD of HR showed a significant group difference in these three variables; FNI higher children had higher SDs than SC.

thumbnail
Table 4. Means (±SE) for heart rate (HR, bpm), standard deviation of HR (SD-HR, bpm), and vagal efficiency (Veff, bpm/ loge units of RSA) for mothers of standard care (SC) and Family Nurture Intervention (FNI) children assessed at 5 years of age (60.0 months corrected age ± 6.0 months SD).

https://doi.org/10.1371/journal.pone.0236930.t004

Relationships between child and mother variables

In addition to analyses of mean values we also computed the relationships between child and mother variables. The partial correlations, adjusting for the same covariates as for the analyses of covariance above, are shown in Table 5.

thumbnail
Table 5. Partial correlations between child and mother lap test autonomic variables at the 4–5 year follow-up visit.

Correlations were adjusted for effects of sex, gestational age, body weight, corrected age at testing, twin status on both the child and mother values.

https://doi.org/10.1371/journal.pone.0236930.t005

Relationships between infant and child autonomic variables

In a prior publication [26] we reported on RSA and Veff values obtained at ~35 weeks and at term age (~40)weeks postmenstrual age, as well for the changes in these measures between the two time-points. The results showed that early life developmental changes in RSA and Veff were significantly greater for FNI infants. Table 6 summarizes correlations between these prior measures of RSA and Veff and the RSA measures obtained at 4–5 years of age. These analyses showed that, across all subjects, RSA at either ~35 or ~40 weeks and in either of two sleep states (quiet sleep, QS; active sleep, AS) were significantly correlated with values obtained at 4 to 5 years of age while the children were awake and on their mother’s lap. While these correlations were positive for both SC and FNI, the correlations within groups reached significance only for FNI subjects.

thumbnail
Table 6. Correlations between RSA and vagal efficiency in infancy (~35 and ~40 weeks postmenstrual age) and RSA at 4 to 5 years of age for all subjects combined and for Standard Care (SC) and Family Nurture Intervention (FNI) subjects separately.

https://doi.org/10.1371/journal.pone.0236930.t006

In another set of analyses, we determined whether the change in RSA from ~40 weeks to 4 and 5 years was significantly different for SC and FNI subjects (Table 6). There were 27 SC and 21 FNI subjects that had RSA at ~40 weeks and again at 4–5 years. Repeated measures ANOVAs for RSA from the infant to child time period showed significant age x group interactions (p<0.05) for data obtained during infancy in both quiet and AS states, after adjusting the ~40 week RSA values for effects of gestational age at birth, birth weight, sex, twin status, and the postmenstrual age at the term age measurement, and adjusting the 4–5 year RSA values for gestational age at birth, birth weight, sex, twin status, and the CA at the 4–5 year time point. The mean (±SE) increases in RSA during QS from ~40 weeks to those when awake at 4–5 years were +3.11±0.16 loge msec2 for SC and +3.67±0.19 loge msec2 for FNI. These results show that the rate of increase in RSA from infancy to childhood is more rapid in FNI subjects.

Discussion

Over the past several decades many interventions have been designed to ameliorate the adverse consequences of premature birth. Some of these have focused primarily on reducing stress, trauma, and anxiety in mothers and fathers of these infants. In some cases, these interventions fall under the general category of Family Centered Care. The goals of these interventions is to help parents and fathers manage the challenges associated with taking care of these more vulnerable infants and, in turn, provide the infants a more supportive environment. Reviews of these interventions, support the value of such approaches for reducing parental distress and, in some cases infant outcomes of weight gain and length of stay, although effect sizes are generally modest (for review of these approaches see [3033]. Other interventions are directed more toward the infants’ experiences in the NICU and regulation of the NICU environment. Included among these are NICCAP [34], sound reduction [35], massage [36], and skin-to-skin or kangaroo care [12]. As shown in the cited reviews, there is some evidence to support use of these interventions as ways to improve growth and neurodevelopmental outcomes of preterm infants, though again, results vary and with few exceptions (e.g. [16]) have not followed infants long term nor have jointly focused on infant and mother outcomes.

In this current report we provide evidence that a new NICU-based intervention, (FNI) alters long-term physiological outcomes of premature infants and their mothers. Specifically, we show that FNI in the NICU increased parasympathetic activity of both mothers and children at 4- to 5-years CA, when compared with mothers and children who received SC alone at the same time point. In other analyses, we combined current data with previous data analyses and found that the rate of change in parasympathetic regulation from ~41 weeks postmenstrual age to 4 to 5 years CA as measured by RSA was greater in FNI children [26]. Both of these findings extend our results at the earlier time point, which showed that FNI infants exhibit more rapid maturation of infant parasympathetic regulation from ~35 to ~41 weeks postmenstrual age. Together, these results indicate that infants who received FNI in the NICU have accelerated development in autonomic regulation from infancy well into early childhood.

In some respects, the finding that a few hours of facilitated mother-infant co-regulatory calming interactions in the NICU can lead to such long lasting changes in physiology may seem surprising. However, extensive animal literature describes the powerful effects of early life experiences on later life behavior and physiology [3740]. Indeed, in studies of mother-pup rearing interactions among normal laboratory rats we found that: 1)”certain types of interactions between mothers and their pups may contribute to individual differences in cardiovascular system development”, 2) “…the cardiovascular system may be shaped by experiences of early life, which are embedded in mother/young behavioral interactions”, and 3)… “that some aspects of maternal behavior may change as a function of pup characteristics” [41, 42]. In other words, evidence has shown for some time that 5-year findings of HRV in the FNI group, which followed a relatively brief calming cycle intervention between mothers and preterm infants in the NICU, are consistent with findings by Feldman et al, who studied the effect of 14 days of SSC in the NICU on autonomic function into childhood.

In the FNI study, the intervention was delivered by Nurture Specialists who were trained to assess emotional connection between mother and infant and to facilitate a connection even before SSC was possible while the infant was still confined to the incubator. Mothers in the SC group were able to engage in nurturing activities of their choosing, which included SSC, non-SSC holding and communication. FNI mothers were asked to engage in non-SSC activities including exchange of scent cloths, comfort touch, vocal soothing, emotional exchange, and SSC when feasible. With regard to the improved scores on the Bayley assessments at 18 months of age, we found effects of FNI were not dependent on the amount of SSC [25]. Therefore, while the autonomic findings are similar in both studies, the differences in intervention protocol raise some intriguing questions about the mechanism of change.

FNI was adapted for the NICU environment from an intervention previously developed by Martha G. Welch MD for older children [43, 44] and was designed to facilitate the emotional connection between mothers and infants in the NICU [17]. The intervention is based on calming cycle theory, which predicts that socioemotional autonomic function can be changed over time via repeated calming sessions between mother and infant [4547]. A central tenet of calming cycle theory is that autonomic Pavlovian co-conditioning between mother and fetus during gestation and the postnatal period promotes co-regulation of autonomic states and emotional connection, which underlies improved socioemotional function over time [17, 45]. Calming cycle theory is consistent with the physiological postulates of polyvagal theory [48, 49], which describe how the autonomic states of the mother and infant can be measured via HRV, cardiac vagal tone (i.e., RSA) and vagal efficiency (Veff) to determine the dyad’s relational health.

Our results show increased RSA in both mother and child during face-to-face contact, which we believe supports the idea that FNI-NICU led to a healthier parasympathetic calming cardiac reflex 4 to 5 years after the intervention. Thus, our findings are consistent with two theories: polyvagal theory, which elucidates the function of the mammalian vagal signaling system and how ANS states provide a physiological scaffold for emotion regulation. Our findings are also consistent with calming cycle theory, which describes a Pavlovian learning mechanism that can be exploited to change autonomic states and socioemotional relationship between mother and infant over time through repeated calming cycles.

In this study we observed higher RSA in both FNI infants and FNI mothers, but neither infants nor mothers in the FNI group had lower HRs than those in the SC group. One possible explanation for this might be related to the level of engagement between FNI mothers and infants. Our results were obtained during a brief 10-min face-to-face “lap test” at the beginning of a battery of assessments, when the pair was instructed to engage with one another as they normally would. Although we have not coded the behavioral interactions at this 4 to 5-year time point, behavioral analysis of mother/infant interactions at 4-months CA showed differences between FNI mothers and their infants [17]. For example, FNI mothers engaged more often in sustained positive touch, followed by brief moments of negative touch, followed again with positive touch. As well, FNI infants exhibited more sustained positive vocal affect and were more likely to transition sooner from negative to positive vocal affect. If greater engagement and reciprocal “give and take” characterize the interactions between FNI mothers and children, one could reasonably conclude that FNI dyads display greater variability in motor activity, when compared with SC dyads. The higher overall variability in HR (SD-HR across 30 sec epochs) found in the current study would support that hypothesis. In as much as greater motor activity would be expected to lead to higher HRs, this might counteract the effect of greater parasympathetic activity on lowering HR.

RSA is a well-established index of efferent parasympathetic input via vagal pathways to the heart (e.g., [50]). However, measurement of ANS function provides more than a good biomarker of early intervention. Our study of FNI-NICU shows that RSA may provide a valuable on-going measure of parent-child socioemotional health as well. Our findings also provide strong support for the hypothesis that there is abundant plasticity in early ANS development that makes it a prime candidate for tracking effectiveness of early interventions via RSA. Prior studies of RSA have shown developmental stability of this measure during early development, for example, from 9 months to 3 years of age [51].The current results extend these developmental correlations to both earlier (~35 weeks postmenstrual age) and later (4 to 5 years corrected) ages. Thus, our data shows RSA to be a predictable measure of autonomic function by demonstrating both responsiveness to early environmental influences and long term individual difference in stability.

Reciprocity in mother and child motor activity might contribute to the positive correlations between mothers and children with regard to HR and standard deviation of HR. However, the correlations between mothers and children for RSA were not significant despite the significant effects of the intervention on both mother and child. This supports the conclusion that, at the group level, FNI has significant effects on both mothers and children. However, at the individual dyad level, in some cases the magnitude of the effect of the intervention on RSA was stronger for the mother than for the child while for other cases the effects were stronger for the child than the mother thus leading to a lack of mother-child correlation. Nonetheless, both calming cycle theory and polyvagal theory predict that the multiple positive prosocial, attentional, and executive function effects reported in the FNI-NICU study should be associated with increases in cardiac vagal tone (RSA) in both mother and child [52, 53]. Indeed, enhanced vagal regulation over the life span should result in multiple psychological and physiological benefits [54]. Such physiological benefits could be quite significant for the mother as well. There is considerable literature to suggest that mothers are at increased risk for developing cardiovascular disease after delivery of a preterm infant [55, 56]. Given the potentially protective influence of cardiac vagal activity on the development of cardiovascular disease [57, 58], we hypothesize that facilitating emotional connection between mothers and preterm infants will provide a buffer to maternal risk for cardiovascular disease.

Limitations

The current findings support the conclusion that FNI during the NICU stay has effects on cardiac vagal regulation lasting for years after discharge; yet, there are important caveats to these findings. One is the fact that the lap test physiological data obtained were from subsets of 55 of 115 (48%) mothers and 71 of 150 (47%) children originally enrolled in the RCT. The most recent follow-up visit prior to 4 years was at 18 months. Given the inter-visit gap of two and a half years, a near 50% follow-up rate in a city and country where relocation is common was not surprising. However, this does raise the issue of whether the FNI and SC mothers and children who returned for this visit represent special groups of the subjects enrolled in the RCT. It is also possible that FNI increases resilience against effects of a stressful environment and group differences would be similar in the lost to follow-up subjects. Accordingly, it is not possible to answer this question with the currently available data.

Indeed, there are many possible events and circumstances that could alter the course of RSA trajectories and, in turn, the effects of FNI observed. Another of these is the degree to which the recommendation made to FNI subjects at discharge to continue calming cycle sessions at home was adhered. It is also possible that some SC mothers and infants achieved the desired emotional connection after discharge on their own, which would introduce additional variance in our results. However, there were no home visits or home intervention aspects to this trial or documentation of childcare activities within the home. Thus, we cannot determine the extent to which these factors influenced our outcomes.

The design of the FNI-NICU RCT included preterm infants with a wide range of gestational ages. Because the study had this wide range and the associated heterogeneity in maternal, fetal and neonatal medical conditions, and morbidities we included gestational age at birth and birth weight as covariates in our analyses. Nonetheless, we cannot know whether the effects of this intervention would be equal for all preterm infants.

The FNI trial design also allowed for enrollment of mothers with twins. Although it is not always clear how independent mothers with twin “A” are from mothers with twin “B”, in general, it is usually best to not consider twins as independent subjects. There are multiple strategies for dealing with this issue including regression modelling with adjustment for clustering/nesting of children within mothers or using a mixed effects model. We opted to present the results two ways, one treating twins as independent samples and a second, more conservative approach of including data from only one twin.

Another limitation of this report not related to the findings per se is our understanding of how FNI can bring about such long-term changes in autonomic regulation. As noted above there are many NICU interventions tested and other new ones have been designed [59]. For all of these approaches we believe there is need to propose and test theories of change and mechanisms that would account for their effectiveness. While we have proposed such a theory of change for FNI, which includes experiencing mother infant physiological co-regulation during calming sessions, establishment of mother/infant emotional connection, and counter conditioning of adverse NICU experiences [60], we did not test this theory in the first RCT of FNI. Thus, explicating underlying mechanisms for early interventions, be they behavioral or physiological, remains a gap in the field of early intervention and preventative medicine approaches to improving outcomes in infants born too early. We are addressing this need in our ongoing replication trial of FNI that includes repeated acquisition of mother and infant physiology to better quantify mother/infant co-regulation and changes in physiological synchrony over time as an index of conditioning. In addition, we will assess emotional connection using the Welch Emotional Connection Screen [61, 62] to test the hypothesis that emotional connection underlies effectiveness of the intervention.

Conclusions

This follow-up study of our previous FNI-NICU trial cohort at 4 to 5 year suggest that FNI may result in healthier autonomic regulation in both mother and children compared to SC. The results are consistent with other clinical and preclinical findings with respect to early postnatal mother/infant interactions and extend the research to longitudinal follow-up in childhood. These results are also consistent with polyvagal theory, which provides the evolutionary and ontogenetic underpinnings of mammalian vagal physiology in prosocial behavior, and calming cycle theory, which describes the mother and infant/child autonomic co-conditioning mechanism through which socioemotional autonomic regulation can be changed from maladaptive to adaptive. Findings support the conclusion that facilitating mother/infant autonomic co-regulation through repeated calming sessions in the NICU is a feasible and promising way to improve neurobehavioral outcomes in preterm infants.

Supporting information

S1 Checklist. CONSORT 2010 checklist of information to include when reporting a randomised trial*.

https://doi.org/10.1371/journal.pone.0236930.s001

(DOC)

Acknowledgments

The authors thank all the participants in this RCT and the NICU staff at Morgan Stanley Children’s Hospital.

References

  1. 1. Murphy SL, Mathews TJ, Martin JA, Minkovitz CS, Strobino DM. Annual Summary of Vital Statistics: 2013–2014. Pediatrics. 2017;139(6):e20163239. pmid:28814547.
  2. 2. Barfield WD. Public Health Implications of Very Preterm Birth. Clin Perinatol. 2018;45(3):565–77. pmid:30144856 PubMed Central PMCID: PMC6193262.
  3. 3. Agrawal S, Rao SC, Bulsara MK, Patole SK. Prevalence of Autism Spectrum Disorder in Preterm Infants: A Meta-analysis. Pediatrics. 2018;142(3):e20180134. pmid:30076190.
  4. 4. Patural H, Barthelemy JC, Pichot V, Mazzocchi C, Teyssier G, Damon G, et al. Birth prematurity determines prolonged autonomic nervous system immaturity. Clin Auton Res. 2004;14(6):391–5. pmid:15666067.
  5. 5. Yiallourou SR, Witcombe NB, Sands SA, Walker AM, Horne RS. The development of autonomic cardiovascular control is altered by preterm birth. Early Hum Dev. 2013;89(3):145–52. pmid:23058299.
  6. 6. Haraldsdottir K, Watson AM, Beshish AG, Pegelow DF, Palta M, Tetri LH, et al. Heart rate recovery after maximal exercise is impaired in healthy young adults born preterm. Eur J Appl Physiol. 2019;119(4):857–66. pmid:30635708.
  7. 7. Mulkey SB, du Plessis AJ. Autonomic nervous system development and its impact on neuropsychiatric outcome. Pediatr Res. 2019;85(2):120–6. pmid:30166644; PubMed Central PMCID: PMC6353676.
  8. 8. Kolacz J, Kovacic KK, Porges SW. Traumatic stress and the autonomic brain-gut connection in development: Polyvagal Theory as an integrative framework for psychosocial and gastrointestinal pathology. Dev Psychobiol. 2019. pmid:30953358.
  9. 9. Feldman R. Mutual influences between child emotion regulation and parent-child reciprocity support development across the first 10 years of life: Implications for developmental psychopathology. Dev Psychopathol. 2015;27(4 Pt 1):1007–23. pmid:26439059.
  10. 10. Korja R, Latva R, Lehtonen L. The effects of preterm birth on mother-infant interaction and attachment during the infant's first two years. Acta Obstet Gynecol Scand. 2012;91(2):164–73. pmid:22007730.
  11. 11. Welch MG, Myers MM. Advances in family-based interventions in the neonatal ICU. Curr Opin Pediatr. 2016;28(2):163–9. pmid:26886783.
  12. 12. Moore ER, Bergman N, Anderson GC, Medley N. Early skin-to-skin contact for mothers and their healthy newborn infants. Cochrane Database Syst Rev. 2016;11:CD003519. pmid:27885658.
  13. 13. Bergman NJ, Linley LL, Fawcus SR. Randomized controlled trial of skin-to-skin contact from birth versus conventional incubator for physiological stabilization in 1200- to 2199-gram newborns. Acta Paediatr. 2004;93(6):779–85. pmid:15244227.
  14. 14. Chi Luong K, Long Nguyen T, Huynh Thi DH, Carrara HP, Bergman NJ. Newly born low birthweight infants stabilise better in skin-to-skin contact than when separated from their mothers: a randomised controlled trial. Acta Paediatr. 2016;105(4):381–90. pmid:26303808.
  15. 15. Feldman R, Eidelman AI. Maternal postpartum behavior and the emergence of infant-mother and infant-father synchrony in preterm and full-term infants: the role of neonatal vagal tone. Dev Psychobiol. 2007;49(3):290–302. pmid:17380505.
  16. 16. Feldman R, Rosenthal Z, Eidelman AI. Maternal-preterm skin-to-skin contact enhances child physiologic organization and cognitive control across the first 10 years of life. Biol Psychiatry. 2014;75(1):56–64. pmid:24094511.
  17. 17. Hane AA, LaCoursiere JN, Mitsuyama M, Wieman S, Ludwig RJ, Kwon KY, et al. The Welch Emotional Connection Screen: validation of a brief mother-infant relational health screen. Acta Paediatr. 2019;108(4):615–25. pmid:29959878.
  18. 18. Welch MG, Hofer MA, Stark RI, Andrews HF, Austin J, Glickstein SB, et al. Randomized controlled trial of Family Nurture Intervention in the NICU: assessments of length of stay, feasibility and safety. BMC Pediatr. 2013;13:148. Epub 2013/09/26. pmid:24063360; PubMed Central PMCID: PMC3851000.
  19. 19. Welch MG, Hofer MA, Stark RI, Andrews HF, Austin J, Glickstein SB, et al. Randomized controlled trial of Family Nurture Intervention in the NICU: assessments of length of stay, feasibility and safety. BMC pediatrics. 2013;13:148. pmid:24063360; PubMed Central PMCID: PMC3851000.
  20. 20. Hane AA, Myers MM, Hofer MA, Ludwig RJ, Halperin MS, Austin J, et al. Family nurture intervention improves the quality of maternal caregiving in the neonatal intensive care unit: evidence from a randomized controlled trial. Journal of developmental and behavioral pediatrics: JDBP. 2015;36(3):188–96. Epub 2015/03/11. pmid:25757070.
  21. 21. Welch MG, Halperin MS, Austin J, Stark RI, Hofer MA, Hane AA, et al. Depression and anxiety symptoms of mothers of preterm infants are decreased at 4 months corrected age with Family Nurture Intervention in the NICU. Arch Womens Ment Health. 2016;19(1):51–61. pmid:25724391.
  22. 22. Welch MG, Myers MM, Grieve PG, Isler JR, Fifer WP, Sahni R, et al. Electroencephalographic activity of preterm infants is increased by Family Nurture Intervention: a randomized controlled trial in the NICU. Clinical neurophysiology: official journal of the International Federation of Clinical Neurophysiology. 2014;125(4):675–84. pmid:24140072
  23. 23. Myers MM, Grieve PG, Stark RI, Isler JR, Hofer MA, Yang J, et al. Family Nurture Intervention in preterm infants alters frontal cortical functional connectivity assessed by EEG coherence. Acta paediatrica (Oslo, Norway: 1992). 2015;104(7):670–7. Epub 2015/03/17. pmid:25772627.
  24. 24. Welch MG, Stark RI, Grieve PG, Ludwig RJ, Isler JR, Barone JL, et al. Family nurture intervention in preterm infants increases early development of cortical activity and independence of regional power trajectories. Acta Paediatr. 2017;106(12):1952–60. pmid:28850710.
  25. 25. Welch MG, Firestein MR, Austin J, Hane AA, Stark RI, Hofer MA, et al. Family Nurture Intervention in the Neonatal Intensive Care Unit improves social-relatedness, attention, and neurodevelopment of preterm infants at 18 months in a randomized controlled trial. Journal of child psychology and psychiatry, and allied disciplines. 2015;56(11):1202–11. Epub 2015/03/13. pmid:25763525.
  26. 26. Porges SW, Davila MI, Lewis GF, Kolacz J, Okonmah-Obazee S, Hane AA, et al. Autonomic regulation of preterm infants is enhanced by Family Nurture Intervention. Dev Psychobiol. 2019. pmid:30868570.
  27. 27. Welch MG, Hofer MA, Brunelli SA, Stark RI, Andrews HF, Austin J, et al. Family nurture intervention (FNI): methods and treatment protocol of a randomized controlled trial in the NICU. BMC pediatrics. 2012;12:14. pmid:22314029; PubMed Central PMCID: PMC3394087.
  28. 28. Hering HE. Functional test of the heart vagi in man. Munch Med Wochenschr. 1910;1910(57):1931–3.
  29. 29. Porges SW, Doussard-Roosevelt JA, Stifter CA, McClenny BD, Riniolo TC. Sleep state and vagal regulation of heart period patterns in the human newborn: an extension of the polyvagal theory. Psychophysiology. 1999;36(1):14–21. pmid:10098376.
  30. 30. Sabnis A, Fojo S, Nayak SS, Lopez E, Tarn DM, Zeltzer L. Reducing parental trauma and stress in neonatal intensive care: systematic review and meta-analysis of hospital interventions. J Perinatol. 2019;39(3):375–86. pmid:30659239.
  31. 31. Ding X, Zhu L, Zhang R, Wang L, Wang TT, Latour JM. Effects of family-centred care interventions on preterm infants and parents in neonatal intensive care units: A systematic review and meta-analysis of randomised controlled trials. Aust Crit Care. 2019;32(1):63–75. pmid:30554939.
  32. 32. Benzies KM, Magill-Evans JE, Hayden KA, Ballantyne M. Key components of early intervention programs for preterm infants and their parents: a systematic review and meta-analysis. BMC Pregnancy Childbirth. 2013;13 Suppl 1:S10. pmid:23445560; PubMed Central PMCID: PMC3561170.
  33. 33. Chertok IR, McCrone S, Parker D, Leslie N. Review of interventions to reduce stress among mothers of infants in the NICU. Adv Neonatal Care. 2014;14(1):30–7. pmid:24472886.
  34. 34. Ohlsson A, Jacobs SE. NIDCAP: a systematic review and meta-analyses of randomized controlled trials. Pediatrics. 2013;131(3):e881–93. pmid:23420913.
  35. 35. Almadhoob A, Ohlsson A. Sound reduction management in the neonatal intensive care unit for preterm or very low birth weight infants. Cochrane Database Syst Rev. 2015;1:CD010333. pmid:25633155.
  36. 36. Alvarez MJ, Fernandez D, Gomez-Salgado J, Rodriguez-Gonzalez D, Roson M, Lapena S. The effects of massage therapy in hospitalized preterm neonates: A systematic review. Int J Nurs Stud. 2017;69:119–36. pmid:28235686.
  37. 37. Denenberg VH. Early experience and emotional development. Sci Am. 1963;208:138–46. pmid:14026976.
  38. 38. Levine S. Developmental determinants of sensitivity and resistance to stress. Psychoneuroendocrinology. 2005;30(10):939–46. Epub 2005/06/17. S0306-4530(05)00090-9 [pii] pmid:15958281.
  39. 39. Francis DD, Caldji C, Champagne F, Plotsky PM, Meaney MJ. The role of corticotropin-releasing factor—norepinephrine systems in mediating the effects of early experience on the development of behavioral and endocrine responses to stress. Biol Psychiatry. 1999;46(9):1153–66. pmid:10560022.
  40. 40. McEwen BS. Early life influences on life-long patterns of behavior and health. Ment Retard Dev Disabil Res Rev. 2003;9(3):149–54. pmid:12953293.
  41. 41. Myers MM, Brunelli SA, Shair HN, Squire JM, Hofer MA. Relationships between maternal behavior of SHR and WKY dams and adult blood pressures of cross-fostered F1 pups. Dev Psychobiol. 1989;22(1):55–67. Epub 1989/01/01. pmid:2912813.
  42. 42. Myers MM, Brunelli SA, Squire JM, Shindeldecker RD, Hofer MA. Maternal behavior of SHR rats and its relationship to offspring blood pressures. Dev Psychobiol. 1989;22(1):29–53. Epub 1989/01/01. pmid:2912812.
  43. 43. Welch MG. Holding time: How to eliminate conflict, temper tantrums, and sibling rivalry and raise happy, loving, sucessful children. New York, N.Y.: Simon and Schuster; 1988.
  44. 44. Welch MG, Northrup RS, Welch-Horan TB, Ludwig RJ, Austin CL, Jacobson JS. Outcomes of Prolonged Parent-Child Embrace Therapy among 102 children with behavioral disorders. Complement Ther Clin Pract. 2006;12(1):3–12. pmid:16401524.
  45. 45. Welch MG. Calming cycle theory: the role of visceral/autonomic learning in early mother and infant/child behaviour and development. Acta Paediatr. 2016;105(11):1266–74. pmid:27536908.
  46. 46. Welch MG, Ludwig RJ. Mother/Infant Emotional Communication Through the Lens of Visceral/Autonomic Learning. Filippa M, Kuhn P, Westrup B, editors. New York, NY: Springer International Publishing; 2017. 271–94 p.
  47. 47. Welch MG, Ludwig RJ. Calming Cycle Theory and the Co-Regulation of Oxytocin. Psychodyn Psychiatry. 2017;45(4):519–40. pmid:29244620.
  48. 48. Porges SW. The polyvagal perspective. Biol Psychol. 2007;74(2):116–43. pmid:17049418; PubMed Central PMCID: PMC1868418.
  49. 49. Porges SW. The polyvagal theory: new insights into adaptive reactions of the autonomic nervous system. Cleve Clin J Med. 2009;76 Suppl 2:S86–90. pmid:19376991; PubMed Central PMCID: PMC3108032.
  50. 50. Lewis GF, Furman SA, McCool MF, Porges SW. Statistical strategies to quantify respiratory sinus arrhythmia: are commonly used metrics equivalent? Biol Psychol. 2012;89(2):349–64. pmid:22138367; PubMed Central PMCID: PMC3269511.
  51. 51. Porges SW, Doussard-Roosevelt JA, Portales AL, Suess PE. Cardiac vagal tone: stability and relation to difficultness in infants and 3-year-olds. Dev Psychobiol. 1994;27(5):289–300. pmid:7926281.
  52. 52. Porges SW. Making the world safe for our children: Down-regulating defense and up-regulating social engagement to 'optimise' the human experience. Children Australia. 2015;40(2):114.
  53. 53. Porges SW, Furman SA. The Early Development of the Autonomic Nervous System Provides a Neural Platform for Social Behavior: A Polyvagal Perspective. Infant Child Dev. 2011;20(1):106–18. pmid:21516219; PubMed Central PMCID: PMC3079208.
  54. 54. Kemp AH, Koenig J, Thayer JF. From psychological moments to mortality: A multidisciplinary synthesis on heart rate variability spanning the continuum of time. Neurosci Biobehav Rev. 2017;83:547–67. pmid:28888535.
  55. 55. Heida KY, Bots ML, de Groot CJ, van Dunne FM, Hammoud NM, Hoek A, et al. Cardiovascular risk management after reproductive and pregnancy-related disorders: A Dutch multidisciplinary evidence-based guideline. Eur J Prev Cardiol. 2016;23(17):1863–79. pmid:27432836.
  56. 56. Tanz LJ, Stuart JJ, Williams PL, Rimm EB, Missmer SA, Rexrode KM, et al. Preterm Delivery and Maternal Cardiovascular Disease in Young and Middle-Aged Adult Women. Circulation. 2017;135(6):578–89. pmid:28153993; PubMed Central PMCID: PMC5308562.
  57. 57. Thayer JF, Yamamoto SS, Brosschot JF. The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int J Cardiol. 2010;141(2):122–31. pmid:19910061.
  58. 58. Kubota Y, Chen LY, Whitsel EA, Folsom AR. Heart rate variability and lifetime risk of cardiovascular disease: the Atherosclerosis Risk in Communities Study. Ann Epidemiol. 2017;27(10):619–25 e2. pmid:29033120; PubMed Central PMCID: PMC5821272.
  59. 59. Silveira RC, Mendes EW, Fuentefria RN, Valentini NC, Procianoy RS. Early intervention program for very low birth weight preterm infants and their parents: a study protocol. BMC pediatrics. 2018;18(1):268. pmid:30092772; PubMed Central PMCID: PMC6085617.
  60. 60. Welch MG. Calming Cycle Theory: The Role of Visceral/Autonomic Learning in Early Mother and Infant/Child Behavior and Development. Acta paediatrica (Oslo, Norway: 1992). 2016. pmid:27536908.
  61. 61. Frosch CA, Fagan MA, Lopez MA, Middlemiss W, Chang M, Hane AA, et al. Validation study showed that ratings on the Welch Emotional Connection Screen at infant age six months are associated with child behavioural problems at age three years. Acta paediatrica (Oslo, Norway: 1992). 2019;108(5):889–95. pmid:30702768.
  62. 62. Hane AA, LaCoursiere JN, Mitsuyama M, Wieman S, Ludwig RJ, Kwon KY, et al. The Welch Emotional Connection Screen: validation of a brief mother-infant relational health screen. Acta paediatrica (Oslo, Norway: 1992). 2019;108(4):615–25. pmid:29959878.