Abstract
PURPOSE:
Transcutaneous auricular vagus nerve stimulation (taVNS) is a non-invasive neuromodulation technique that may improve oromotor skills when paired with feeding in at-risk infants, but effects on other motor function and how motor function relates to white matter (WM) microstructure are unknown.
METHODS:
In this prospective study, infants failing oral feeds and slated for gastrostomy tube (G-tube) placement received taVNS paired with bottle feeding daily for 2-3 weeks. The effects of taVNS-paired feeding on general and specific head movements were investigated using the Specific Test of Early infant motor Performance (STEP) and diffusion MRI obtained before and after taVNS treatment. Scores between and within groups (taVNS responders, attained full oral feeds; non-responders, received G-tubes) were compared.
RESULTS:
Performance on head movement items improved significantly in responders but not in non-responders (p < 0.05). Total STEP scores were significantly higher in responders after taVNS treatment than non-responders (p = 0.04). One STEP item, rolling by arm, was associated with significantly greater change in WM tract microstructure (p < 0.05) in the responders.
CONCLUSION:
These results suggest that pairing feeding with taVNS may affect specific head and neck movements to a greater extent in infants who are able to attain full oral feeds.
Abbreviations
Anterior limb of internal capsule
Central nervous system
Cerebral palsy
Diffusional kurtosis imaging
Diffusion MRI
External capsule
Fractional anisotropy
Fidgety movements
Gestational age
Gastrostomy tube
General Movements Assessment
Gross Motor Function Classification System
Inferior frontooccipital fasciculus
Superior longitudinal fasciculus
Specific Test of Early infant motor Performance
Term equivalent age
transcutaneous auricular vagus nerve stimulation
Tract-based spatial statistics
Vagus nerve stimulation
White matter
Background
Non-invasive transcutaneous auricular vagal nerve stimulation (taVNS) paired with infant feeding is a novel therapeutic intervention for infants born preterm who exhibit motor or feeding delays or both [1, 2]. The vagus nerve comprises an extensive synaptic network of afferent and efferent projections integrating the central nervous system (CNS) and autonomic control circuitry of the brain stem to vital internal organs [3]. Vagal afferent projections, including those in the cervical and auricular branches, terminate in numerous nuclei, specifically, the locus coeruleus, raphe nuclei, thalami, and the basal forebrain. Electrical impulses delivered to these nuclei via the vagus nerve in the neck or the auricular branch at the ear result in neurotransmitter release that influences cortical plasticity [4–6].
Multiple preclinical and clinical studies in adults indicate that vagus nerve stimulation (VNS) drives cortical plasticity when combined with a simultaneous sensory or motor activity. Animal studies demonstrate that VNS reorganizes the motor cortex via acetylcholine release in the nucleus basalis, resulting in long-lasting functional performance [4–6]. Data from several recent small and large randomized controlled trials suggest that surgically implanted cervical VNS in conjunction with intensive task practice improves upper-motor function in adult stroke rehabilitation [7–11]. However, implanted cervical VNS is relatively expensive, requires surgery, and may result in adverse events such as wound infection, vocal cord palsy, or hoarseness [8–11]. Non-invasive taVNS targets the largely afferent auricular branch of the vagus nerve, is safe, and improves hand motor function in adult stroke survivors [8, 12]. Further, randomized control trials of VNS and taVNS show greater motor improvements that lasted for at least three months in subacute stroke patients after rehabilitation paired with VNS compared with conventional rehabilitation alone [11, 14].
In a first-of-its-kind infant study, it was demonstrated that taVNS is both feasible and associated with improved oromotor-feeding skills with attainment of full oral feeds in over half of infants failing oral feeds who were gastrostomy tube (G-tube) candidates [1, 2]. Brain stimulation via taVNS was delivered when the infant was actively sucking from a bottle once or twice daily for 2-3 weeks [1, 2]. Most recently, the US Food and Drug Administration designated the taVNS feeding system a Breakthrough Medical Device.
Aside from improved oromotor skills, the authors’ interests lie in evaluating the effect that taVNS treatment might have on other early motor skills. The Specific Test of Early infant motor Performance (STEP) is a novel developmental screening test constructed of ten itemized foundational movements crucial for development of refined motor skills, such as head control against gravity. The ten items on the STEP assess anti-gravity flexion and extension of the head and neck, movement in the arms and legs, and tone in the shoulder girdle and pelvis. Specific items on the STEP reflect the maturation of head movements in the first four months of life and encapsulate the infant’s ability to stabilize the head at midline in supine, prone, side-lying, and supported sitting [15].
To demonstrate the possible effects of this novel taVNS brain stimulation technique on CNS plasticity in a cohort of infants who failed to learn the fundamental motor sequence of feeding, diffusion MRI (dMRI) data and STEP assessments were obtained before and after the treatment course of taVNS-paired feeding. dMRI measures the random thermal movement of water molecules within the underlying tissue microstructure that make it acutely sensitive during the rapid progression, organization and myelination of developing white matter (WM). Of the myriad diffusion metrics, the fractional anisotropy (FA), a dimensionless rotational invariant expressing the water movement’s directional dependence, is the most commonly employed. The FA represents a surrogate biomarker of WM structural integrity and maturation. While studies have shown a correlation between dMRI measures and motor recovery in adult patients with chronic stroke, few studies have directly explored the relationship between FA values and motor learning changes after participation in a short-term, directed intervention [16–19].
It was postulated that circuit activity driven by taVNS-paired feeding would elicit cortical plasticity in and around WM tracts related to feeding (i.e., increased FA values) that correlate to infant motor skills as a whole. This substudy investigated whether this novel therapeutic intervention influenced 1) overall motor skills, 2) specific motor skills of neck and neck control, and 3) WM tracts involved in motor and sensorimotor functions. Moreover, it was of interest whether these skills were different in taVNS responders who attained full oral feeds versus taVNS non-responders who eventually required G-tube placement for oral feeds.
Methods
Design
This study was approved by the Institutional Review Board at the Medical University of South Carolina. Parental consent was obtained prior to enrollment into the taVNS paired-feeding study (Clinicaltrials.gov registration NCT# 04643808). Infants who had failed oral feeding trials and were referred for G-tube placement were recruited sequentially. Infants were either premature (born at <33 weeks gestational age [GA]) or suffered CNS injury such as hypoxic ischemic encephalopathy and were stable without significant respiratory support or unrepaired congenital anomalies that limited feeding volumes. Infants were deemed safe to attempt to take every feed by mouth by occupational or speech therapists.
Intervention
Infants received taVNS paired with bottle feeds once or twice a day for 2-3 weeks. The taVNS treatment was delivered via a left ear electrode with electrical pulses at 25 Hertz, a 500 microseconds pulse width, and current intensity at 0.1 miliAmpere below the infant’s perceptual threshold [1, 20]. The microcurrent stimulation was delivered while the infant was actively sucking and was off during periods of rest or burping; bottle-feeding sessions lasted up to 30 minutes. Infants were categorized as treatment ‘responders’ to taVNS treatment if they achieved full oral feeds (>120 ml/kg/d) with sufficient weight gain for nursery discharge. The clinical results and a detailed description of the stimulation intervention have been reported in prior studies [1, 20].
Outcome measures
STEP
Early motor abilities were measured via the STEP, which was administered by trained pediatric occupational therapists before the start of taVNS intervention (pre-STEP) and upon completion of the intervention (post-STEP) in the neonatal intensive care unit. The STEP’s ten itemized foundational movements can determine if an infant is at risk for developmental delays as early as term and 3 months of age [2]. The STEP takes both qualitative and quantitative aspects of the infant’s movements into consideration and clearly discriminates among preterm infants of different motor abilities [15, 22] (Fig. 1). Excellent sensitivity and specificity of the STEP cut-off scores at term (≤16) and 3 months (≤22) were previously demonstrated to predict the Bayley-III gross motor performance at 12 months [15]. When compared with the Test of Infant Motor Performance, the STEP showed better prediction of delays at 12 months [15].

STEP items related to head and neck movements.
dMRI data were acquired pre- and post-taVNS treatment on a Skyra 3T MRI (Siemens Healthineers, Erlangen, Germany) immediately following feeding with swaddling to induce natural sleep. A diffusional kurtosis imaging (DKI) [23] protocol utilized b-values of b = 0, 1000 and 2000 s/mm2, (3 mm)3 isotropic voxels, echo time = 122 ms, repetition time = 6700 ms with 30 diffusion-encoding directions per non-zero b-value and 40 contiguous axial slices per image volume. Additionally, 9 non-diffusion weighted (b = 0 s/mm2) images were also acquired with matching imaging parameters.
Image volumes and their corresponding diffusion-encoding directions that were corrupted by motion or signal dropout were manually removed prior to image processing. The DESIGNER pipeline [24], reformatted into Python 3 [25], was employed for the full DKI analysis. The pipeline includes Marchenko-Pastur-PCA denoising [26], Gibbs artifact correction [27], eddy correction [28, 29], Gaussian smoothing with a kernel 1.25 times the voxel size [30] and Rician noise bias correction [31]. The output of DESIGNER is tensor-derived diffusion and kurtosis metrics; the FA was specifically focused in on for this investigation.
Data analysis
Change in total STEP scores and STEP item scores between responders and non-responders
Total STEP scores were compared from pre-taVNS to post-taVNS in the responders and non-responders groups by paired t-test. The change in total STEP scores between responder and non-responder groups was compared by independent t-test and normalized for the weeks of development between assessments (ΔSTEP/weeks). To determine specific movements that might be responsive to taVNS-paired feeding intervention, the change in specific STEP item scores within each group was compared by paired t-test and between the two groups by independent t-test, normalized for weeks of development (ΔSTEP/weeks). Statistical analyses were conducted using IBM SPSS (version 21).
FA study template
A neonatal FA study template was generated based on a user-chosen subject FA image (i.e., target) using the optimized method described by Ball et al. [32]. Briefly, all FA images were brain-extracted [33] and registered to the target FA image using a 6 and then 12 degree-of-freedom registration [34]. The resulting affine transformation matrices were concatenated and served as the initialization for non-linear warping of each subject’s FA image to the target FA image [35, 36].
The study-specific FA template resulted from taking the average across all target FA aligned images. This mean FA template was then registered to the Johns Hopkins University neonatal FA template [37] to create the final standard space FA template using similar steps outlined above. The core WM voxels common to all subjects were identified and, using FA <0.15, a binary WM skeleton mask was generated.
The change in FA (ΔFA) after completing taVNS treatment was obtained by taking the difference of the post- and pre-treatment FA maps. Subjects were then dichotomized by total and individual STEP item scores and labeled high (total score >16; item score of 2 or 3) or low performing (total score ≤16; item score of 0 or 1) on the STEP.
A tract-based spatial statistics (TBSS) analysis [38] was then performed for each dichotomized ΔFA set separately while adjusting for the GA difference from pre- to post-taVNS dMRI scans. Individual STEP items tested were head supported sitting, prone, pull to sit, rolling by arm and rolling by leg. A two-sided t-test was employed using the randomize function of FSL, which is a method to perform voxel-wise non-parametric permutation inference testing along with threshold-free cluster enhancement to correct for multiple comparisons [39].
Results
Participants
A total of 27 of the 35 infants in the taVNS feeding study had at least one completed STEP performed before or after taVNS treatment (15 non-responders, 12 responders) (Fig. 2). There were no significant differences between the two groups in GA, birth weight, age at enrollment, or length of time attempting full intended volume for oral feeds prior to receiving taVNS intervention (Table 1). However, there were statistically significant differences in the number of taVNS treatments and days to achieve full oral feeds. Infants in the responder group required fewer taVNS treatments (p = 0.04) and were able to achieve full oral feeds within an average of 10.2±5.6 days. The STEP subgroup included more non-responders than responders, unlike the parent study of 35 infants in which 54% responded to taVNS treatment.

CONSORT flow diagram of the study participants.
Patient demographics and characteristics
Data presented as mean±standard deviation. *Significant at p < 0.05 level; ∧p = 0.005 between pre- and post-STEP in Responders. ap-values from independent t-test or chi-square test between Responders and Non-responders. bData are not applicable. cPatient 11: was transferred patient from another hospital and had >30 days of attempting PO.
Of the total 27 infants, 19 had the STEP performed before and after taVNS treatment (12 non-responders, 7 responders). The average number of days between the pre- and post-STEP tests were similar between the responders (15.3) and non-responders (23.1) (p = 0.22). The STEP scores before taVNS were not significantly different by responder group pre-taVNS treatment (Table 2). Infants who were responders had a significantly greater increase in total STEP score (mean 1.7±2.1) from pre- to post-taVNS treatment than infants who were non-responders and required a G-tube (0.0±0.85, p = 0.04, paired t-test, Fig. 3A, Table 2). Expected change in total STEP score for high performing infants was 1 point over 2 weeks [15]. In addition, responders had significantly greater increase in four specific STEP items directly related to movements in the head and neck compared with non-responders. These included head in supine with visual stimulation (p = 0.01), head in supine with no visual stimulation (p = 0.03), rolling elicited by arm (p = 0.04), and head movements in supported sitting (p = 0.01) (Fig. 3B-D , Table 3). No changes were observed in the non-responder group from before to after treatment in total STEP scores or in any specific STEP item scores (Table 3).
Total STEP scores and specific STEP items between responder and non-responder infants who had both pre- and post-STEP scores; Mean (SD)
Total STEP scores and specific STEP items between responder and non-responder infants who had both pre- and post-STEP scores; Mean (SD)
Data presented as difference in means (post item score - pre item score)±standard deviation. *Significant at P < 0.05 level. a Change in STEP scores divided by weeks between STEP, normalized by weeks of development. Independent t-test was used in all analysis.

Difference in the change in total STEP score (A) and specific STEP items (B, C, D, E) between the taVNS Reponders versus Non-Responders.
Change in total and specific STEP items from pre- to post-taVNS between the Responder and Non-Responder groups, Mean (SD)
*Significant at P < 0.05 level. Paired t-test was used in all analysis.
While taVNS treatment was performed over 2-3 weeks, gains in motor skills or tone might be observed even in this relatively short time period. Analyzing the change in STEP scores normalized by weeks between the pre- and post-STEP assessments (Table 2), there was still a significant difference in the change in total STEP score (p = 0.04) and in three items related to head and neck movement: head in supine with visual stimulation (p = 0.03), head in supine with no visual stimulation (p = 0.03), and head movements in supported sitting (p = 0.01).
Of the total 19 infants with STEP scores, 18 had an MRI performed before and after taVNS treatment (9 non-responders, 9 responders). To determine if the improvements in STEP score were related to WM tract plasticity changes, TBSS analyses of ΔFA between pre- and post-taVNS treatment were performed and dichotomized by high and low performance on total and individual STEP item scores, adjusting for GA between the two MRI scans. There was a significantly greater ΔFA within multiple WM regions (p < 0.05) after taVNS treatment in infants who performed well in rolling by arm compared to infants who performed poorly (Fig. 4). The regions highlighted as different in the rolling by arm analysis correspond to primarily left sided portions of major WM tracts: the anterior limb of internal capsule (ALIC), external capsule (EC), inferior fronto-occipital fasciculus (IFOF), and superior longitudinal fasciculus (SLF) (Fig. 4). The TBSS analysis dichotomized by total STEP performance had no significant findings.

WM tract regions are different following TBSS comparing high and low performing infants in Rolling by Arm (panel A) and Rolling by Leg (panel B) with WM tract atlas regions for refence (panel C). (A) Red-yellow voxels denote the WM tracts with significantly greater positive ΔFA in infants who performed well on STEP item Rolling by Arm item after taVNS treatment vs infants who performed poorly (p < 0.05). The left hemispheric tracts were: anterior and posterior limbs of internal capsule (ALIC and PLIC), external capsule (EC), inferior frontooccipital fasciculus (IFOF), and corona radiata (CR). (B) Regions in red denote WM tracts with significantly greater ΔFA in infants who performed well on STEP item Rolling by Leg item after taVNS treatment vs infants who performed poorly (p < 0.08). The right hemispheric WM clusters correspond to parts of the EC, as well as the ALIC and PLIC, frontal corona radiata, cortico-spinal tract and lower optic radiations. (C) Colored regions from the John’s Hopkins neonatal WM atlas for reference and orientation to the relative location of the IFOF (green), EC (yellow), ALIC (blue), PLIC (light-blue), CR (pink) and optic radiations (red).
A trend to greater ΔFA in infants who performed well on the item Rolling by Leg pre- and post- taVNS treatment (p = 0.1, threshold for TBSS) was observed (Fig. 4). In contrast to rolling by arm, the highlighted WM clusters were heavily right hemispheric with a few clusters of significant voxels along the left optic radiation. The right-sided WM clusters correspond to regions of the EC and ALIC as well as the posterior limb of internal capsule. The SLF, frontal corona radiata and cortico-spinal tract were also included in the trending WM regions.
In this open label study, the possible effects of taVNS paired with oral feeding on the general and specific movements of the head and neck was investigated in a cohort of infants who failed to learn the fundamental motor sequence of feeding. Splitting the group into responders and non-responders allowed comparison of STEP assessments before and after the treatment course of taVNS-paired feeding and investigation of performance on specific STEP items to dMRI data. These results reveal that, although pre-taVNS STEP scores were not different between groups, responders had greater changes in motor function over the short 2-3 week paired taVNS intervention when compared to non-responders.
These intriguing data show that this improvement in overall motor function is largely driven by increases in head and neck control. While bottle-feeding involves coordination of 22 oral-facial muscles for successful intake, the muscles of the head and neck are also activated during feeding and may be affected by the activity of bottle feeding. These results suggest that when feeding is paired with taVNS, these movements may be affected to a greater extent in infants who are able to attain full oral feeds. As STEP total and rolling by arm item scores are not significantly different between responders and non-responders pre-taVNS paired bottle feeding, it is conceivable that taVNS paired with bottle feeding facilitated integration of these other motor circuits. Additionally, from dMRI, FA revealed greater improvements in motor tract WM integrity after taVNS treatment in infants who were high performing on specific STEP items related to head and neck control, as demonstrated during an axial rolling movement. This suggests that tissue anisotropy is increasing more rapidly in high performing infants than in their low performing counterparts, suggesting that increased WM maturation in specific WM tracts is likely reflecting induced neuroplasticity after taVNS-paired bottle feeding.
This is the first study in infants to demonstrate that early taVNS treatment for feeding delays may be associated with improvements in critical early motor movements, especially those related to maturity of antigravity head and neck control. The findings are consistent with the improvement in motor function when VNS was paired with a functional task in adult stroke survivors, demonstrating a generalized, long-lasting effect on motor performance [8, 40]. In addition, studies have shown that VNS paired with intensive motor practice in rats post-stroke promotes reorganization of the primary cortex corresponding to the task footprint and results in greater functional improvement compared to motor practice alone [7]. When VNS was paired with wheel spin activity, cortical representation of the distal forelimb increased, whereas with the lever press activity, cortical representation of more proximal forelimb activity increased compared to rats who practiced the activity alone [41].
While cortical reorganization cannot be directly measured, the maturation of WM tracts involved in these motor activities can be measured via FA. The FA imparts information on the directional dependence of water movement, such that highly aligned, tightly packed fiber bundles will have a large FA (closer to 1) whereas crossing fiber regions or loosely aligned bundles will be representative of lower FA values (closer to 0). The interpretation of FA is complicated by the crossing fiber problem, in which two high FA fibers may cross at right angles resulting in FA = 0, which is not a correct measure of the water movement. However, axonal alignment and progressive rapid myelination during early development has a role in the clear and dramatic increase in FA seen within the WM of the genu and splenium of the corpus callosum. This FA increase stems directly from decreased water mobility across the lipid membrane now sheathing the axons. In healthy development, this FA increase begins to plateau at approximately two to five years of age. In the current study, STEP movements that involve head and neck control were increased in responders versus non-responders after only a brief taVNS paired oral feeding intervention. In addition, greater maturity in the axial rolling items was associated with greater ΔFA values in several major WM tracts involved in coordination of sensorimotor responses (ALIC, EC, IFOF and SLF).
Studies have also shown the importance of post-term axial rolling, in which the infant’s whole body is turned from supine to/from prone and the infant aligns head and neck with the trunk. A new assessment termed the Motor Optimality Score for infants, which is based on the work of Prechtl, includes axial rolling as one of the items [42, 43]. In infants with neurodevelopmental impairment and cerebral palsy (CP), rolling to the side and antigravity movements were decreased compared with typically developing infants [44]. In addition, this study found that movements of the arms and legs to midline, which require antigravity flexion, emerge earlier in neurotypical infants. The item rolling to the side can also help discriminate between grades of CP severity using the Gross Motor Function Classification System (GMFCS), specifically between children classified with GMFCS I-II vs. GMFCS III-V unilateral vs. bilateral CP [44].
Two prior studies have also determined a link between an infant’s feeding abilities and early motor performance using the General Movements Assessment (GMA) [45, 46]. The GMA assesses the quality of spontaneous movements, namely fidgety movements (FM), in the first five months and has been used to predict CP [47–49]. When combined with MRI and neurological examination, the GMA has shown excellent sensitivity and specificity for detecting CP in high-risk infants [50]. Nieuwenhuis et al. [45] showed that abnormal FMs are associated with uncoordinated sucking patterns near term equivalent age (TEA) [45]. Moreover, infants with abnormal GMA scores at term were observed to have a higher risk of developing oromotor-feeding impairment at 12 months [46]. Their data indicated only a weak association of brain abnormalities on structural MRI at TEA and future oromotor-feeding impairment. Nevertheless, taken together, these studies suggest a link between an infant’s early motor performance and early feeding abilities [45, 46].
The current study is limited by the open label design and the relatively small sample size. However, in this first-in-infants study on the impact of taVNS on at-risk infants, multi-modal data was collected on each infant enrolled in the study. This allowed for in-depth analysis of specific early infant motor movements, which have been shown to be significant in other published studies [37, 42]. In addition, the important relationship of dMRI metrics to motor learning changes in infant head control after a brief but focused motor-paired taVNS intervention was explored.
Conclusion
This study provides proof-of-concept that early neuromodulation paired with the motor task of feeding in at-risk infants may have positive and potentially long-lasting effects on motor movements required for head and neck control. When taVNS is delivered while the infant is actively engaging muscles of the head and neck during bottle-feeding, significant changes were observed not just in feeding behaviors but also in other early head and neck movements. The long-term impact of taVNS treatment in these and hopefully future study infants and how improvements in early head and neck movements may impact overall motor skill development later in life will continue to be documented. In addition, future studies should determine the potential long-term impact of taVNS paired with oromotor activity on speech and language domains.
Footnotes
Acknowledgments
Funding for this study was provided by the National Institutes of Health National Center of Neuromodulation for Rehabilitation NIH/NICHD Grant #P2CHD086844; National Institutes of Health Center of Biomedical Research Excellence (COBRE) in Stroke Recovery, Grant #P20GM109040; and NIH NINDS F31NS108623 (H. Moss). We would like to acknowledge the following people who contributed to this study: Sarah Huffman, Morgan Dancy, Will Devries, Sean Thompson and Georgia O’Leary.
Conflict of interest
The authors have no conflicts of interest to report.
