- 1Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children's Hospital, Boston, MA, United States
- 2Harvard Medical School, Boston, MA, United States
Previously, we reported quantitatively smaller total corpus callosum (CC) and total forebrain size in critically ill term-born and premature patients following complex perioperative critical care for long-gap esophageal atresia (LGEA) that included Foker process repair. We extended our cross-sectional pilot study to determine sub-regional volumes of CC and forebrain using structural brain MRI. Our objective was to evaluate region-specific CC as an in-vivo marker for decreased myelination and/or cortical neural loss of homotopic-like sub-regions of the forebrain. Term-born (n = 13) and premature (n = 13) patients, and healthy naïve controls (n = 21) <1-year corrected age underwent non-sedated MRI using a 3T Siemens scanner, as per IRB approval at Boston Children's Hospital following completion of clinical treatment for Foker process. We used ITK-SNAP (v.3.6) to manually segment six sub-regions of CC and eight sub-regions of forebrain as per previously reported methodology. Group differences were assessed using a general linear model univariate analysis with corrected age at scan as a covariate. Our analysis implicates globally smaller CC and forebrain with sub-region II (viz. rostral body of CC known to connect to pre-motor cortex) to be least affected in comparison to other CC sub-regions in LGEA patients. Our report of smaller subgenual forebrain implicates (mal)adaptation in limbic circuits development in selected group of infant patients following LGEA repair. Future studies should include diffusion tractography studies of CC in further evaluation of what appears to represent global decrease in homotopic-like CC/forebrain size following complex perioperative critical care of infants born with LGEA.
Highlights
- Infants born with long-gap esophageal atresia (LGEA) have decreased forebrain and corpus callosum volumes following Foker process repair
- We report global decrease in homotopic-like CC/forebrain size following complex perioperative critical care of infants born with LGEA
- Decreased non-homotopic subgenual region of the forebrain implicates (mal)adaptation in limbic circuits development in studied group of infant patients.
Introduction
The corpus callosum (CC) represents the major and largest commissural tract in the human brain connecting two cerebral hemispheres. Indeed, callosal axons originate primarily from neurons in layer II/II and layer V of the neocortex (1, 2). Those callosal axons are arranged in topographically ordered manner within the CC, before their projection into their corresponding contralateral cortical areas (3). As such, the CC is divided into a total of six anatomical sub-regions that correspond to topographically ordered sub-regions of forebrain (4), also described as homotopic relationship (5). Development of the CC and its projections in the forebrain occurs rapidly in the third trimester of pregnancy and throughout the first year of life (6, 7). Prematurity and/or critical illness is associated with hindered CC/forebrain development and disruption of subsequent developmental processes (8, 9). As CC myelination has been shown to occur anteriorly to posteriorly (10), we seek to evaluate sub-region-specific CC and forebrain volumes as to identify potential early structural vulnerability that can occur in the context of the complex perioperative critical care in infancy.
Our recent pilot study identified a group of critically ill infants born with long-gap esophageal atresia (LGEA) to have smaller corpus callosum (11) and forebrain (12) volumes following complex perioperative critical care. Unlike short-gap defects that can be repaired by direct anastomosis (13), infants born with long esophageal disconnects (>3 cm) underwent complex perioperative critical care involving tension-induced esophageal growth known as the Foker process (14–16) requiring prolonged sedation ≥5 days (17–19) leading to physical dependence to the drugs of sedation (17, 20, 21). Such complex perioperative critical care spans period of weeks, as described previously [Figure 1 in (19, 22)].
We hypothesized that when compared to healthy infants, both critically ill term- and preterm-born patients following perioperative critical care for LGEA repair with Foker process would exhibit (mal)adaptation in sub-region-specific CC and forebrain volumes. Thus, our novel analysis aimed to quantify homotopic-like CC and forebrain sub-regional volumes using T1-weighted brain MRI in our recent pilot cohort of infants following repair of LGEA (22). Our broader objective was to evaluate sub-region-specific CC volume as a possible in-vivo marker for decreased myelination and/or cortical neural loss of homotopic sub-regions of the forebrain. We used previously published methods of semi-automated segmentation of the CC and forebrain that relied on distinct anatomical landmarks in the infant brain, as described by Witelson (23) and Peterson et al. (24, 25), respectively. Anatomically divided sub-regions of the CC and forebrain were schematically related to each other according to the existing literature and in the spirit of a homotopic relationship (26). Our preliminary data were previously presented only as virtual abstract at a conference (27).
Methods
Study Design and Participants
Our report builds on our previous pilot infant MRI study that received ethical approval from Institutional Review Board as a “no more than minimal risk” study. Specifically, we extended our previous T1-weighted brain MRI analysis of total CC volumes (11) and T2-weighted analysis of total forebrain volumes (12) to the current T1-weighted analysis of sub-regional volumes of both CC, as well as that of the forebrain. Since presented data is from the same infant study subjects, previously described methodological approach (11, 22) for (1) recruitment criteria and (2) MRI scanning process apply to this study as well. Briefly, informed written parental consent was obtained prior to subject participation, in accordance with the Declaration of Helsinki and Good Clinical Practice guidelines. The family of each subject received a $90 gift card for participation in the non-sedated research brain MRI scan.
Patients With LGEA
Eligibility criteria included both term-born [37–42 weeks gestational age (GA) at birth] and preterm-born patients <1 year gestation-corrected age. The preterm-born group included only very preterm (28 to < 32 weeks GA) and moderate or late preterm-born infants (32 to < 37 weeks GA) as defined by The World Health Organization (28). All patients underwent initial LGEA repair at either Boston Children's Hospital or an outside institution in the first 2 months of life. Specifically, Foker process (14–16) allows for tension-induced esophageal growth for long esophageal atresia disconnects (viz. > 3 cm found in LGEA). Such complex perioperative care involves several stages (29): (1) Foker I thoracotomy to place traction sutures onto blind esophageal ends; (2) Esophageal lengthening by continuous traction on the esophagus pouches to induce esophageal tissue growth which requires sedation; (3) Foker II thoracotomy to approximate esophageal ends and perform primary esophageal anastomosis; (4) Post-Foker healing of the anastomosis with weaning of sedation and transition from total parenteral nutrition to enteral feeds prior to hospital discharge. The unique aspect of such complex perioperative care in infancy is not only related to repeated anesthesia and surgery, but to prolonged sedation exposure ≥5 days known to be associated with development of physical dependence to the drugs of sedation (17, 20, 21). Representative timeline illustrating sequence of perioperative critical care for LGEA repair was presented previously (17–19). Associations between individual MRI end-point measures (e.g., number of cranial MRI findings and brain volumes) and the clinical measures of care as to assess the severity of underlying disease in cohort patients will be presented elsewhere.
Exclusion criteria included: (1) extreme prematurity (<28 weeks GA); (2) diagnosis of small for gestational age and/or intrauterine growth restriction (SGA/IUGR) (30, 31); (3) extracorporeal membrane oxygenation exposure; (4) clinically indicated cranial ultrasound findings (e.g., ventricular enlargement with or without gray matter and/or ventricular hemorrhage); (5) neurological disease as documented in clinical record (e.g., seizures); (6) chromosomal abnormalities (e.g., Down's syndrome); (7) prenatal drug exposure to either drugs of abuse or prescription medications; and/or (8) MRI incompatible implants. Indeed, we recruited only those patients born with LGEA that had no clinical evidence of neurological problems at the time of recruitment as per detailed chart review and/or cranial ultrasound findings when available.
Term-Born Controls
Healthy term-born infants (<1 year old) with no prior exposure to surgery, anesthesia, or sedation (i.e., naïve to perioperative treatment) were recruited from a pool of Boston Children's Hospital outpatients and two neighboring newborn centers (Beth Israel Deaconess Medical Center and Brigham and Women's Hospital, Boston, MA). These controls served as a reference baseline for typical sub-regional CC and forebrain size. This study differs from our previous T1-weighted analysis [Table 1 in (11)] by addition of one extra control infant (total n = 21). An updated summary of final group characteristics is shown in Table 1.
MRI Acquisition
All infant patients underwent non-sedated research brain MRI at a point in time just before the hospital discharge following completion of initial surgical repair, or during subsequent admissions for follow up management in the 1st year of life. Thus, subjects were scanned throughout the infancy, depending on the time of their recruitment. A “feed and wrap” approach was used on all infants undergoing a non-sedated research brain MRI scan after the completion of complex perioperative treatment (32–35). Corrected age at scan for all cohort subjects was calculated as follows: postnatal age (weeks) – [40 – GA at birth (weeks)]. Please, refer to our previous publication in regard to individual total CC [Figure 5 in (11)] and total forebrain volumes [Figure 3 in (12)] with respect to corrected GA at scan. Patients were scanned at late evenings or nights using a 3T TrioTim MRI system equipped with 32-channel receive-only head coil and body-transmission (Siemens Healthcare Inc., USA) as per previously described protocol (12, 22). Foam earplugs (Newmatic Medical, Birmingham, AL) and earmuffs (MRI-Safe Neonatal Noise Guards, Universal Medical, Norwood, MA) were used for noise protection. Smaller infants were further supported using beanbags (viz. gentle vacuum bag immobilizer), while infants older than 3 months were allowed to assume a more relaxed position (e.g., arms next to face). Despite undergoing a non-sedated research scan, due to complexity and severity of clinical status, a single physician (DB) continuously monitored all infants for a stable heart rate and oxygenation throughout the MRI acquisition. Structural T1-weighted images were acquired using a MPRAGE sequence (repetition time = 2.52 s; echo time = 1.74 ms; flip angle = 7°; field of view = 192 × 192 mm2; voxel size = 1 × 1 × 1 mm3; 144 sagittal slices). We collected 100% of T1-weighted images for all patients (n = 13/group), and 91% (n = 21/23) for term-born controls (Table 1). Regarding the latter, one infant had partial brain coverage that precluded analysis of total brain volume (n = 20 controls for forebrain analysis) but allowed for CC analysis (n = 21 controls). Although we noted ringing artifact due to motion only in 1/21 controls and 1/13 premature patients, it was mild enough not to require motion correction or obscure tissue segmentation.
Infant Brain Segmentation of T1-Weighted MRI Data
Prior to any tissue segmentation, Freeview (v.2.0) from Freesurfer (the General Hospital Corporation, Boston, MA) was used to correct for any head tilt during MRI acquisition. In other words, preprocessing of T1-weighted images included alignment along the anterior commissure - posterior commissure (AC-PC) line [see Figure 2 in (22)]. A single rater with neuroanatomical expertise blindly performed tissue segmentation, which was subsequently checked by a senior researcher for internal consistency. Specifically, ITK-SNAP software (v.3.6.0; www.itksnap.org) (36) was used for total and sub-regional CC segmentation (Figure 1) and total brain and forebrain segmentation (Figure 2A–C) of T1-weighted images. We also performed surface area analysis, but since finals 2-D data results did not significantly differ from volumetric data, 2D results are not presented.
Figure 1. Sub-regional corpus callosum analysis. Representative T1-weighted corpus callosum (CC) segmentation for a term-born otherwise healthy control (A,A') and a term-born patient (B,B') scanned at 10 and 9 months of age, respectively. Volumes of CC were estimated for a total of 6 CC sub-regions (C). Sub-regional segmentation used the proportion-based parcellation scheme originally proposed by Witelson (23), with slight modification based on Venkatasubramanian et al. (37) (see Methods for details). A, anterior; G, genu; P, posterior.
Figure 2. Sub-regional forebrain analysis. T1-weighted regional forebrain analysis is shown for a term-born otherwise healthy control scanned at 7 months of age. (A) shows T1 image in midsagittal (left) and axial (right) views, noting the key anatomical landmarks used for sub-regional forebrain analysis [the anterior-most aspect of corpus callosum (CC), the anterior commissure (AC), and the posterior commissure (PC)]. Following total brain segmentation (B), the cerebellum (CB) and brainstem (BS) were manually erased to create a mask of total forebrain. Subsequently, forebrain was parcellated into eight sub-regions (C), as per Peterson et al. (24). Orbitofrontal (OF; orange) and dorsal prefrontal (DF; red) sub-regions were combined and analyzed as the prefrontal cortex (PFC) sub-region. Color pattern in the corpus callosum outline (D) schematically illustrates the homotopic-like relationship between different sub-regions of CC and forebrain according to existing literature (26). Specifically, CC region I (genu and rostrum) projects to PFC, region II (rostral body) projects to PM, region III (anterior midbody) projects to SM, region IV (posterior midbody) projects to SM and PO (C, yellow stripes), region V (isthmus) projects to PO, and region VI (splenium) projects to IO and MT. The subgenual (SG) area of forebrain is predominantly comprised of limbic structures, with no major direct projections from the CC. A, anterior; L, left; P, posterior; R, right.
CC Segmentation
We previously described a detailed method for manual tracing of T1-weighted data for total CC volume (cm3) segmentations [Figure 3 in (11)]. The CC sub-regional parcellation scheme used in this study is based on previously established functional topographical division across the white matter tract (38). Furthermore, manual parcellation of the CC sub-regions was done in accordance with previously validated approach by Witelson (23), with slight modifications based on Venkatasubramanian et al. (37). Specifically, we performed sub-regional CC segmentation (Figure 1) as follows:
CC Volumetric Segmentation
Manual CC segmentation (Figure 1A',B') was done in line with previously described protocol for CC segmentation in newborns by Yu et al. (39) to exclude encircling cerebral vasculature and neighboring tissues including the subcallosal area, cingulum, and fornix (columns, body, and crux). Previously, we extensively described and illustrated [Figure 3B–F in (11)] lateral boundaries of CC segmentation: lateral outermost edges of the CC were bounded by the anterior and posterior corona radiata (40).
Sub-regional CC Segmentation
Illustration of parcellation scheme for sub-regional CC division is shown in Figure 1C. As previously described [Figure 2A in (11)], a linear segment running from the most anterior tip of the genu and the most posterior point of the splenium in the mid-sagittal section (Figure 1A,B) was used to create an anterior – posterior (A-P) line that represents CC length in the sagittal view (Figure 1C). ITK-SNAP software (v.3.6.0; www.itksnap.org) (36) was used to measure this A-P length (cm) of the CC. A-P line served for subsequent sub-regional segmentation of CC (Figure 1C). Specifically, regional boundaries were established along five coronal planes corresponding to the crook of the genu (Figure 1C - point “G”), as well as at intervals of one third, one half, two thirds, and four fifths the total length of the A-P line. This resulted in parcellation of 6 CC sub-regions as follows: (I) genu (includes rostrum), (II) rostral body (III) anterior midbody, (IV) posterior midbody, (V) isthmus, and (VI) splenium. In alignment with previous literature (41), rostrum was included as part of the genu segmentation given its very small size and highly variable shape in infancy.
Forebrain Segmentation
We adapted aspects of our method of total brain and forebrain segmentation previously used for T2-weighted analysis (12, 22) to T1-weighted MRI data used in this study. Briefly, we performed the semi-automated approach as follows:
Total Brain and Forebrain Volumetric Segmentation
Total brain volume required the following semi-automated tissue segmentation steps: (i) Skull-stripping of T1-weighted images via manual tracing of the whole brain outline (Figure 2B); and (ii) Partial volume segmentation of cerebrospinal fluid (CSF) using FMRIB's Automated Segmentation Tool (FAST) (42). We used tools in FMRIB Software Library (FSL; v.5.0) to eliminate ventricular system included in the whole brain mask. Specifically, CSF partial volume estimate was (a) thresholded at 99% (eliminating voxels with <99% of their volume comprising CSF), (b) converted to a binary mask, and (c) subtracted from a mask of the whole brain outline in order to generate mask of total brain volume excluding the ventricular system, which underwent additional (d) minor manual editing to draw-in any missing brain tissue, resulting in a finalized total brain mask (Figure 2B). T1-weighted data analysis of total brain volumes in infant with LGEA has been previously published (43) while total brain volume masks served for total and sub-regional forebrain segmentation in this study. Total forebrain tissue mask was created by simply erasing the cerebellum and brainstem (medulla, pons and midbrain) from a total brain mask (Figure 2B,C).
Sub-regional Forebrain Segmentation
Due to ongoing development and myelination throughout a patient's first year of life, there is a decreased gray to white matter tissue contrast in infant brains (44) that poses challenges for automated infant brain segmentation. We selected sub-regional analysis of the forebrain based on previously published scheme by Peterson et al. (24, 25). The reliability and validity of related schemes of brain subdivision have been previously documented (24, 25, 45). Each forebrain mask was divided into eight anatomical sub-regions using one axial plane containing the AC-PC line, and three coronal planes (tangent to genu of CC, the anterior commissure, and the posterior commissure; Figure 2A,C). Thus, the eight forebrain sub-regions, as illustrated in Figure 2C, were segmented: orbitofrontal (OF), dorsal prefrontal (DF), premotor (PM), sensorimotor (SM), parieto-occipital (PO), inferior occipital (IO), midtemporal (MT), and subgenual (SG). Given the small size of the two most anterior sub-regions – the orbitofrontal and the dorsal prefrontal sub-regions were combined and defined as the prefrontal cortex (PFC) sub-region. This resulted in a total of seven forebrain sub-regions included in the final analysis. Note that majority of terminology of each forebrain sub-region is related to the cortical nomenclature, although listed sub-regions encompass the cortex and subcortical tissue of the forebrain.
Structural Quantification
The ITK-SNAP volume estimation tool was used to obtain absolute volumes (cm3) of segmentation masks. Normalization was calculated as part of a whole to allow for better understanding as to how a particular sub-region of interest changes with respect to the whole (viz. total forebrain volume). Specifically, volumetric data of both CC and forebrain sub-regions were normalized as a % total forebrain volume (Figures 3, 4, respectively).
Figure 3. Sub-regional corpus callosum volumes. Graphs display individual normalized sub-regional volumes as a % total forebrain (FB) volume for the three study groups: (1) term-born controls (n = 20; white box), (2) term-born patients (n = 13; black box), and (3) preterm-born patients (n = 13; gray box). As noted in Methods section, forebrain volume used for normalization was not available for 1/21 controls due to incomplete brain coverage, which accounts for final control group n = 20. Normalized volumes as % FB of all CC sub-regions, except for the rostral body, were significantly smaller in both term-born and preterm patients compared to controls (p ≤ 0.002; Table 2), but not between patient groups. Dots represent individual data values, boxes represent the interquartile range (25–75%), and solid horizontal lines represent medians. Roman numerals for each graph also depict classification of CC sub-regions as outlined in the Figure 1.
Figure 4. Sub-regional forebrain volumes. Graphs display individual normalized sub-regional volumes (A–G) as a % forebrain (FB) volume for the three study groups: (1) term-born controls (n = 20; white box), (2) term-born patients (n = 13; black box), and (3) preterm-born patients (n = 13; gray box). As noted in Methods section, forebrain volume was not available for 1/21 controls due to incomplete brain coverage, which accounts for final control group n = 20. We report no group differences for normalized sub-regional forebrain volumes with exception for normalized parieto-occipital (D) and subgenual (G) sub-regions between term-born patients and controls (p ≤ 0.01) and between patient groups (p ≤ 0.04), but not preterm patients and controls (p ≥ 0.42). See also Table 3 for complete statistics. The latter findings may be attributed to higher variability of data in premature patients and should be interpreted with caution. Dots represent individual data values, boxes represent the interquartile range (25–75%), and solid horizontal lines represent medians.
Statistical Analysis
This study used original data from our previous pilot study (22). Since no prior information was available regarding brain findings in the selected cohort of infants with LGEA, we used a convenience sample size of 13 patients/group, based on the anticipated number of eligible infants at our institution and an estimated 50% enrollment rate. Statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS, v.23.0; IBM Corporation, Armonk, NY), and normal distribution of all continuous variables was confirmed using the Shapiro-Wilk test. To account for the potential confounding variable of having subjects scanned at various ages throughout the first year of life, comparisons of sub-regional volumes between the three groups were assessed using a general linear model (GLM) univariate analysis with corrected age at scan as a covariate and Bonferroni adjusted p-values. Statistical significance was assessed at the α <0.05 level after conservative Bonferroni correction that accounted for multiple comparisons. In this report, we present differences between three groups studied. A comprehensive statistical report for absolute and normalized group differences is also condensed in a table format (Tables 2, 3).
Results
In this report, we share quantitative analysis of T1-weighted data for six CC and seven forebrain sub-regional volumes among term- and preterm-born patients following LGEA repair (n = 13/group), and term-born controls (n = 21; Table 1).
Sub-regional Analysis of Corpus Callosum
We previously reported disproportionately smaller total CC volume (as % of total brain volume) in both patient groups (11). Our novel analysis of total CC as a % of total forebrain volume matched these previous findings. Specifically, normalized total CC volumes showed significant group differences [F(2,42) = 16.16, p < 0.001], with both term-born and preterm-born patients having smaller normalized total CC volumes compared to controls, with no differences between patient groups (graph not shown).
To evaluate sub-regional vulnerability of CC, we analyzed volumes of six CC sub-regions as illustrated in Figure 1A',B'. Statistical analysis of absolute values (cm3) between groups is summarized in Table 2. We normalized all six sub-regions of the CC as % total forebrain to see if any of the specific CC sub-regions in patients are smaller when compared to the forebrain as a whole. We report that normalized CC sub-regional volumes (as % forebrain volume; Figure 3) were significantly smaller in both patient groups compared to controls in all sub-regions except for the sub-region II (rostral body). For full statistical details, (see Table 2). As such, results implicate the rostral body, sub-region II of the CC to be least affected.
Sub-regional Analysis of Forebrain
Our current T1-weighted analysis parallels the results of previously published T2-weighted data report (12) that showed proportionally smaller total forebrain (as % of total brain volume). Specifically, no group differences were observed in normalized volumes for total forebrain [F(2,42) = 0.1, p = 0.90; graph not shown].
We extended the study to perform 3-D segmentation of seven forebrain sub-regions according to previous report in literature (24, 25), as illustrated in Figure 2C. Of the parcellated seven forebrain sub-regions, six are known to send homotopic projections via selective sub-regions of CC (26), as schematically color-coded in Figure 2C,D. In the absence of better methodology, presented results should be interpreted with caution due to linear parcellation methodological approach not fully aligned to functional connectivity (see section Study Limitations: Methodological Consideration). Statistical analysis of decreased absolute values (cm3) per sub-region in patients is summarized in Table 3. However, we only found group differences for normalized sub-regional volumes of the forebrain (as % total forebrain volume; Figure 4) for two forebrain sub-regions: parieto-occipital and subgenual sub-regions (term-born patients were different from both controls and premature patients). For full statistical details, (see Table 3). Future studies with more power and narrower clinical criteria for inclusion of premature patients should be considered.
Discussion
Complex perioperative critical care with Foker process (14–16) in infants born with LGEA represents a standard of care. In this manuscript, we bring forth data with intention to bridge the gap in our understanding of such complex care's impact on brain development in infants born with LGEA. This work was built upon previous findings of disproportionally smaller total CC size in both term-born and premature infants following complex LGEA repair (11) and was extended to an investigation of homotopic-like sub-regional CC and forebrain volume analysis. Due to pilot study nature and methodological limitations of the forebrain segmentation in infancy, authors caution the readers against extrapolation of data until future studies confirm presented results.
Homotopic-Like Analysis of Corpus Callosum and Forebrain Volumes
Our current data of sub-regional CC size implicates globally smaller CC in premature and term-born patients. The uniformly decreased normalized CC volume with respect to the whole forebrain [with the exception of the sub-region II of the CC (viz. rostral body)], may implicate delayed and/or impaired maturation of white matter tracts in the setting of LGEA as a congenital anomaly and/or result of LGEA complex perioperative treatment. As such, CC volume analysis may serve as an early marker of brain (mal)adaptations of brain size in the selected group of infants born with LGEA undergoing complex perioperative critical care with Foker process (14–16). Our volumetric analysis showed disproportionally larger normalized volume in the rostral body (sub-region II of the CC) known to connect to the pre-motor cortex, indicating that this region may be least likely affected (schematics summary in Figure 5) in the context of globally smaller brain in selected cohort of infants (22). The anterior callosal fibers such as those projecting from the rostral body have been shown to interconnect to the frontal lobe and relate motor information (6). Larger normalized rostral body findings, and homotopic-like linkages between the rostral body and the premotor area, which was also shown to have potentially altered growth trajectories for the patient groups, may indicate possible protection in this area of the forebrain, an interesting and unexpected finding given the previously established delay in motor development in critically ill infants (46). Conversely, the more adversely affected posterior CC fibers connecting the parietal, temporal, and occipital lobes, the greater likelihood of (mal)adaptations in integration of sensory information (6) might be encountered. To our knowledge, no studies yet have evaluated the neurodevelopmental outcomes in either term- or preterm-born infants following complex perioperative critical care with Foker process for LGEA repair. As such, future clinical evaluations that include motor and sensory assessments should be used to correlate neuroanatomical findings to the neurodevelopmental implications.
Figure 5. Summary schematics of corpus callosum to forebrain homotopic-like projections. Schematic illustrates homotopic projections between corpus callosum and the forebrain sub-regions. In this study, we showed consistent decrease in volumes across sub-regions of CC (Table 2; Figure 3) and forebrain (Table 3; Figure 4) in selected group of infants following complex perioperative critical care for long-gap esophageal repair. Global decrease of the volumes appears least affected in the region II of the CC (viz. rostral body of the CC) known to project to the premotor cortex (green arrow) following esophageal repair with Foker process. Interestingly, a disproportionally smaller normalized volume of subgenual sub-region of the forebrain (not illustrated; see Figure 2C) does not have any direct homotopic connections via CC.
Although no previous studies have investigated sub-regional CC development in the setting of complex perioperative critical care, previous work has shown that the more posterior sub-regions of the CC (viz. anterior midbody, posterior midbody, isthmus, and splenium; Figure 1) are likely to have smaller surface areas, decreased volume of interhemispheric fiber tracts, and reduced fractional anisotropies in very premature babies (41). This is consistent with our findings of decreased volumes in such posterior areas, and possible protective effects in the more anterior rostral body sub-region (sub-region II of CC). Importantly, as CC myelination occurs anteriorly to posteriorly (10), such findings may indicate the posterior sub-regions of the CC are more vulnerable to adverse early life events in the context of the complex perioperative critical care for LGEA repair.
Our current sub-regional analysis of normalized forebrain volumes failed to identify vulnerable sub-region of the forebrain in homotopic-like relation to CC sub-regions. As such, results indicate that uniformly decreased sub-regional forebrain volumes contribute to the global decrease in brain size (22) in the selected pilot cohort of infants following LGEA repair. Interestingly, a disproportionally smaller normalized volume is reported only for the subgenual sub-region of the forebrain (Figure 2C), which does not have any homotopic connections with the CC (Figure 5), and suggest that, possibly - limbic circuits may be adversely impacted following complex perioperative critical care for LGEA repair.
Underlying Brain Mechanisms
The disproportionately smaller CC sub-regional volumes may indicate hypo or incomplete myelination in term-born and premature patients. Whether this is due to decreased oligodendrocyte cells, decreased axonal density, or degradation of the myelin, remains unclear. It should be noted, that previous studies have linked prematurity to focal non-cystic white matter injuries in nearly half of premature newborns with MRI data, and much of this white matter vulnerability has been linked to late oligodendrocyte progenitor cells/ subplate cells (47). Whether this is true of the underlying pathology of term-born infants undergoing complex perioperative care warrants future studies. Previous literature has also reported term-born infants undergoing complex care, including those both with and without major cardiac anomalies, are more likely to have white matter lesions (cardiac and non-cardiac cases) and lower fractional anisotropy in the CC (cardiac cases) (48, 49). Whether the decreased CC size in our sample population of infants with LGEA reflects preexisting brain abnormalities, risk of abnormal development in the setting of perioperative critical care, or a double-hit etiology of brain injury both pre and peri-Foker process, as previously discussed (18), remains unknown. Our most recent preliminary data published as a small case report series (18) are in support of possible double-hit etiology of brain injury (pre-Foker and peri-Foker process) in both premature and term-born infants with LGEA. Future studies in this vulnerable infant population should also include diffusion tensor imaging (DTI) and tractography to shed light on the possible mechanisms affecting white matter microstructural integrity during infancy, a critical age for the white matter maturation and global brain growth.
Neurobehavioral Significance
The exact neurological effects and long-term neurobehavioral impact of complex perioperative care for LGEA repair remain unknown. Previous reports indicate that infants with non-cardiac congenital anomalies that undergo complex perioperative care are at an increased risk for abnormal brain development and poor long-term outcomes (9, 48). More broadly, as the largest white matter tract in the brain that serves to connect and facilitate cross-talk between the two hemispheres, abnormal development of the CC in infancy may lead to altered cognitive/executive function that persists into childhood and adolescence (50). Prematurity, white matter impairment, and decreased CC volume have also all been associated with poor executive functioning and language outcomes (51–53). Since functional implications of reported structural differences remains unknown in infants following complex perioperative critical care for LGEA repair, investigation into clinical correlations and long-term neurodevelopmental outcomes represents the next steps in clinical research investigation.
Study Limitations
Methodological Considerations
As the infant forebrain is undergoing constant rapid myelination, a clear parcellation scheme based on functional significance has yet to be developed. Rather, a previously validated parcellation scheme based on distinct anatomical features rather than functional significance was utilized. Homotopic-like connections between the forebrain and CC were then established according to existing literature (26). In other words, since the parcellation was not performed based on the true connectivity between the corpus callosum and the forebrain, current parcellation scheme does not fully reflect functional relation. Future studies following technical automated segmentation advancement of the forebrain to more defined structures, could be used to assess a more structurally/functionally precise method of understanding homotopic CC-forebrain relationship.
Study Timing
Previous literature has estimated robust linear postnatal growth of cortical gray-matter volume in the first year of life, and comparatively linear—but slower growth of cortical white-matter volume (54). Future studies should include cross-sectional analysis to account for linear growth rates differences between gray and white matter in the first year of life.
Study Size
Small sample size and a lack of subjects scanned at older time points (>8 months of age) pose potential limitations. Age-related changes were taken into consideration by using a general linear model with gestation-corrected age as a covariate, but future studies with larger sample size and a more uniform age distribution are necessary.
Control Group
We only used term-born control infants that served as a comparison to normative size. However, future studies should consider additional controls such as (1) a cohort of non-cardiac LGEA that underwent alternative treatment to the Foker process, (2) a cohort that received only prolonged sedation without surgery, and (3) a cohort of premature infants that received no medical care.
Sex Differences
Previous studies remain inconsistent, with some reporting no sex differences in CC area, and other reporting larger normalized female CC volumes (6, 55). Although we had equal distribution of sex/patient group, our term-born control group consisted predominately of male infants. Future studies should have enough power to allow for sex differences analysis of both CC volumes and its homotopic-like projections to the forebrain.
Pre-existing Findings
Care was taken to exclude infants with previously existing neurological disease or brain abnormalities. However, infants were not given MRI scans prior to treatment, so it is impossible to know if size (mal)adaptations in CC/forebrain size pre-existed to the complex perioperative critical care, as is the case with pre-surgical differences in CC volume of patients with congenital heart defects that has been shown to worsen post-surgery (56).
Conclusions
Our novel sub-regional quantitative analysis strengthens previous observations of qualitative brain atrophy and quantitatively decreased global and regional brain size in a unique pilot cohort of term-born and preterm-born infant patients without previously recognized neurological injury or insult. Future longitudinal studies should include diffusion tractography of the CC and forebrain in further evaluation of what appears to represent posteriorly diminished homotopic-like CC/forebrain volume, and protection of rostral body/premotor regions, following complex treatment of infants born with LGEA.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Boston Children's Hospital Institutional Review Board (IRB-P000007855). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
CM and DB: the conception and manuscript design. MK, CM, and DB: acquisition. MK, DZ, and DB: analysis. MK and DB: drafting the article. All authors: interpretation of data, critical revision for important intellectual content, final approval of the version to be published, and accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding
This work was supported by the NIDA K08 DA035972-01 (to DB) and 2017 Trailblazer Award from the Department of Anesthesiology, Critical Care and Pain Medicine, Boston Children's Hospital (to DB).
Author Disclaimer
The content of this article is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors gratefully acknowledge the infants and their parents who participated in this study. Authors would also like to thank: (i) Russell W. Jennings MD, Dorothy Gallagher RN, and Jean Solodiuk RN, PhD for their help with recruitment; (ii) All MRI technologists for their invaluable help with scanning infants; (iii) colleagues from the Computational Radiology Laboratory, directed by Dr. Simon K. Warfield, PhD for their technical support, as well as (iv) Sonia A. Main, BS Candidate and Madhuri S. Jois, BS Candidate for their help with preprocessing T1-weighted MRI images. Authors express gratitude to both Reviewers for their constructive feedback and invaluable help.
Abbreviations
AC-PC, anterior commissure—posterior commissure; A-P, anterior—posterior; CA, corrected age; CC, corpus callosum; CSF, cerebrospinal fluid; FAST, FMRIB's Automated Segmentation Tool; FSL, FMRIB Software Library; FB, forebrain; DF, dorsal prefrontal sub-region; IO, inferior occipital sub-region; MT, midtemporal sub-region; OF, orbitofrontal sub-region; PFC, prefrontal cortex combined sub-region; PM, premotor sub-region; PO, parieto-occipital sub-region; SG, subgenual sub-region; SM, sensorimotor sub-region; GA, gestational age; GLM, general linear model; IUGR, intrauterine growth restriction; LGEA, long-gap esophageal atresia; MRI, magnetic resonance imaging; SGA, small for gestational age; SPSS, Statistical Package for the Social Sciences; TEF, Tracheo-esophageal fistula.
References
1. Wise SP. The laminar organization of certain afferent and efferent fiber systems in the rat somatosensory cortex. Brain Res. (1975) 90:139–42. doi: 10.1016/0006-8993(75)90688-5
2. Jones EG, Wise SP. Size, laminar and columnar distribution of efferent cells in the sensory-motor cortex of monkeys. J Comp Neurol. (1977) 175:391–438. doi: 10.1002/cne.901750403
3. Zhou J, Wen Y, She L, Sui YN, Liu L, Richards LJ, et al. Axon position within the corpus callosum determines contralateral cortical projection. Proc Natl Acad Sci USA. (2013) 110:E2714–23. doi: 10.1073/pnas.1310233110
4. Goldstein A, Covington BP, Mahabadi N, Mesfin FB. Neuroanatomy, Corpus Callosum. Treasure Island, FL: StatPearls (2020).
5. Wahl M, Strominger Z, Jeremy RJ, Barkovich AJ, Wakahiro M, Sherr EH, et al. Variability of homotopic and heterotopic callosal connectivity in partial agenesis of the corpus callosum: a 3T diffusion tensor imaging and Q-ball tractography study. Am J Neuroradiol. (2009) 30:282–9. doi: 10.3174/ajnr.A1361
6. Tanaka-Arakawa MM, Matsui M, Tanaka C, Uematsu A, Uda S, Miura K, et al. Developmental changes in the corpus callosum from infancy to early adulthood: a structural magnetic resonance imaging study. PLoS ONE. (2015) 10:e0118760. doi: 10.1371/journal.pone.0118760
7. Malinger G, Zakut H. The corpus callosum: normal fetal development as shown by transvaginal sonography. Am J Roentgenol. (1993) 161:1041–3. doi: 10.2214/ajr.161.5.8273605
8. Iwata S, Katayama R, Kinoshita M, Saikusa M, Araki Y, Takashima S, et al. Region-specific growth restriction of brain following preterm birth. Sci Rep. (2016) 6:33995. doi: 10.1038/srep33995
9. Stolwijk LJ, Keunen K, de Vries LS, Groenendaal F, van der Zee DC, van Herwaarden MYA, et al. Neonatal surgery for noncardiac congenital anomalies: neonates at risk of brain injury. J Pediatr. (2017) 182:335–41 e1. doi: 10.1016/j.jpeds.2016.11.080
10. McLaughlin NC, Paul RH, Grieve SM, Williams LM, Laidlaw D, DiCarlo M, et al. Diffusion tensor imaging of the corpus callosum: a cross-sectional study across the lifespan. Int J Dev Neurosci. (2007) 25:215–21. doi: 10.1016/j.ijdevneu.2007.03.008
11. Mongerson CRL, Jaimes C, Zurakowski D, Jennings RW, Bajic D. Infant corpus callosum size after surgery and critical care for long-gap esophageal atresia: qualitative and quantitative MRI. Sci Rep. (2020) 10:6408. doi: 10.1038/s41598-020-63212-3
12. Mongerson CRL, Jennings RW, Zurakowski D, Bajic D. Quantitative MRI study of infant regional brain size following surgery for long-gap esophageal atresia requiring prolonged critical care. Int J Dev Neurosci. (2019) 79:11–20. doi: 10.1016/j.ijdevneu.2019.09.005
13. Hunt RW, Perkins EJ, King S. Peri-operative management of neonates with oesophageal atresia and tracheo-oesophageal fistula. Paediatr Respir Rev. (2016) 19:3–9. doi: 10.1016/j.prrv.2016.01.002
14. Foker JE, Kendall Krosch TC, Catton K, Munro F, Khan KM. Long-gap esophageal atresia treated by growth induction: the biological potential and early follow-up results. Semin Pediatr Surg. (2009) 18:23–9. doi: 10.1053/j.sempedsurg.2008.10.005
15. Kunisaki SM, Foker JE. Surgical advances in the fetus and neonate: esophageal atresia. Clin Perinatol. (2012) 39:349–61. doi: 10.1016/j.clp.2012.04.007
16. Bairdain S, Hamilton TE, Smithers CJ, Manfredi M, Ngo P, Gallagher D, et al. Foker process for the correction of long gap esophageal atresia: primary treatment versus secondary treatment after prior esophageal surgery. J Pediatr Surg. (2015) 50:933–7. doi: 10.1016/j.jpedsurg.2015.03.010
17. Solodiuk JC, Jennings RW, Bajic D. Evaluation of postnatal sedation in full-term infants. Brain Sci. (2019) 9:114–30. doi: 10.3390/brainsci9050114
18. Rudisill SS, Wang JE, Jaimes C, Mongerson CRL, Hansen AR, Jennings RW, et al. Neurologic injury and brain growth in the setting of long-gap esophageal atresia perioperative critical care: a pilot study. Brain Sci. (2019) 9:383–405. doi: 10.3390/brainsci9120383
19. Hodkinson DJ, Mongerson CRL, Jennings RW, Bajic D. Neonatal functional brain maturation in the context of perioperative critical care and pain management: a case report. Heliyon. (2019) 5:e02350. doi: 10.1016/j.heliyon.2019.e02350
20. Anand KJ, Barton BA, McIntosh N, Lagercrantz H, Pelausa E, Young TE, et al. Analgesia and sedation in preterm neonates who require ventilatory support: results from the NOPAIN trial. Neonatal Outcome and Prolonged Analgesia in Neonates. Arch Pediatr Adolesc Med. (1999) 153:331–8. doi: 10.1001/archpedi.153.4.331
21. Dewey WL. Various factors which affect the rate of development of tolerance and physical dependence to abused drugs. NIDA Res Monogr. (1984) 54:39–49. doi: 10.1037/e472472004-001
22. Mongerson CRL, Wilcox SL, Goins SM, Pier DB, Zurakowski D, Jennings RW, et al. Infant brain structural MRI analysis in the context of thoracic noncardiac surgery and critical care. Front Pediatr. (2019) 7:315. doi: 10.3389/fped.2019.00315
23. Witelson SF. Hand and sex differences in the isthmus and genu of the human corpus callosum. A postmortem morphological study. Brain. (1989) 112 (Pt 3):799–835. doi: 10.1093/brain/112.3.799
24. Peterson BS, Vohr B, Staib LH, Cannistraci CJ, Dolberg A, Schneider KC, et al. Regional brain volume abnormalities and long-term cognitive outcome in preterm infants. JAMA. (2000) 284:1939–47. doi: 10.1001/jama.284.15.1939
25. Peterson BS, Anderson AW, Ehrenkranz R, Staib LH, Tageldin M, Colson E, et al. Regional brain volumes and their later neurodevelopmental correlates in term and preterm infants. Pediatrics. (2003) 111:939–48. doi: 10.1542/peds.111.5.939
26. Meyer JW, Makris N, Bates JF, Caviness VS, Kennedy DN. MRI-Based topographic parcellation of human cerebral white matter. Neuroimage. (1999) 9:1–17. doi: 10.1006/nimg.1998.0383
27. Kagan MS, Mongerson CRL, Jois MS, Main SA, Zurakowski D, Jennings RW, et al. Infant regional corpus callosum and forebrain volumes after surgery for long-gap esophageal atresia. In: Organization for Human Brain Mapping (OHBM) Abstract 2020, June 23-July 3, 2020 Annual Meeting, Montreal, Canada (2020).
28. Quinn JA, Munoz FM, Gonik B, Frau L, Cutland C, Mallett-Moore T, et al. Preterm birth: case definition & guidelines for data collection, analysis, and presentation of immunisation safety data. Vaccine. (2016) 34:6047–56. doi: 10.1016/j.vaccine.2016.03.045
29. Liszewski MC, Bairdain S, Buonomo C, Jennings RW, Taylor GA. Imaging of long gap esophageal atresia and the Foker process: expected findings and complications. Pediatr Radiol. (2014) 44:467–75. doi: 10.1007/s00247-013-2847-2
30. Sharma D, Shastri S, Farahbakhsh N, Sharma P. Intrauterine growth restriction - part 1. J Matern Fetal Neonatal Med. (2016) 29:3977–87. doi: 10.3109/14767058.2016.1152249
31. Tsudo M, Uchiyama T, Takatsuki K, Uchino H, Yodoi J. Modulation of Tac antigen on activated human T cells by anti-Tac monoclonal antibody. J Immunol. (1982) 129:592–5.
32. Almli CR, Rivkin MJ, McKinstry RC, Brain Development Cooperative Group. The NIH MRI study of normal brain development (Objective-2): newborns, infants, toddlers, and preschoolers. Neuroimage. (2007) 35:308–25. doi: 10.1016/j.neuroimage.2006.08.058
33. Raschle N, Zuk J, Ortiz-Mantilla S, Sliva DD, Franceschi A, Grant PE, et al. Pediatric neuroimaging in early childhood and infancy: challenges and practical guidelines. Ann NY Acad Sci. (2012) 1252:43–50. doi: 10.1111/j.1749-6632.2012.06457.x
34. Tocchio S, Kline-Fath B, Kanal E, Schmithorst VJ, Panigrahy A. MRI evaluation and safety in the developing brain. Semin Perinatol. (2015) 39:73–104. doi: 10.1053/j.semperi.2015.01.002
35. Paterson SJ, Badridze N, Flax JF, Liu WC, Benasich AA. A Method for Structural MRI Scanning of Non-Sedated Infants. Chicago, IL: Rutgers University (2004).
36. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. (2006) 31:1116–28. doi: 10.1016/j.neuroimage.2006.01.015
37. Venkatasubramanian G, Anthony G, Reddy US, Reddy VV, Jayakumar PN, Benegal V. Corpus callosum abnormalities associated with greater externalizing behaviors in subjects at high risk for alcohol dependence. Psychiatry Res. (2007) 156:209–15. doi: 10.1016/j.pscychresns.2006.12.010
38. Fabri M, Pierpaoli C, Barbaresi P, Polonara G. Functional topography of the corpus callosum investigated by DTI and fMRI. World J Radiol. (2014) 6:895–906. doi: 10.4329/wjr.v6.i12.895
39. Yu X, Zhang Y, Lasky RE, Datta S, Parikh NA, Narayana PA. Comprehensive brain MRI segmentation in high risk preterm newborns. PLoS ONE. (2010) 5:e13874. doi: 10.1371/journal.pone.0013874
40. Wakana S, Jiang H, Nagae-Poetscher LM, van Zijl PC, Mori S. Fiber tract-based atlas of human white matter anatomy. Radiology. (2004) 230:77–87. doi: 10.1148/radiol.2301021640
41. Thompson DK, Inder TE, Faggian N, Johnston L, Warfield SK, Anderson PJ, et al. Characterization of the corpus callosum in very preterm and full-term infants utilizing MRI. Neuroimage. (2011) 55:479–90. doi: 10.1016/j.neuroimage.2010.12.025
42. Zhang Y, Brady M, Smith S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging. (2001) 20:45–57. doi: 10.1109/42.906424
43. Kagan MS, Mongerson CRL, Zurakowski D, Jennings RW, Bajic D. Infant study of hemispheric asymmetry after long-gap esophageal atresia repair. Ann Clin Transl Neurol. (2021) 8:2132–45. doi: 10.1002/acn3.51465
44. Huang H, Zhang J, Wakana S, Zhang W, Ren T, Richards LJ, et al. White and gray matter development in human fetal, newborn and pediatric brains. Neuroimage. (2006) 33:27–38. doi: 10.1016/j.neuroimage.2006.06.009
45. Caviness VS Jr, Meyer J, Makris N, Kennedy DN. MRI-based topographic parcellation of human neocortex: an anatomically specified method with estimate of reliability. J Cogn Neurosci. (1996) 8:566–87. doi: 10.1162/jocn.1996.8.6.566
46. Yang CF, Xue Y, Feng JY, Jia FY, Zhang Y, Li YM. Gross motor developmental dysfunctional outcomes in infantile and toddler pediatric intensive care unit survivors. BMC Pediatr. (2019) 19:508. doi: 10.1186/s12887-019-1893-9
47. Miller SP, Ferriero DM. From selective vulnerability to connectivity: insights from newborn brain imaging. Trends Neurosci. (2009) 32:496–505. doi: 10.1016/j.tins.2009.05.010
48. Stolwijk LJ, Lemmers PM, Harmsen M, Groenendaal F, de Vries LS, van der Zee DC, et al. Neurodevelopmental outcomes after neonatal surgery for major noncardiac anomalies. Pediatrics. (2016) 137:e20151728. doi: 10.1542/peds.2015-1728
49. Mulkey SB, Ou X, Ramakrishnaiah RH, Glasier CM, Swearingen CJ, Melguizo MS, et al. White matter injury in newborns with congenital heart disease: a diffusion tensor imaging study. Pediatr Neurol. (2014) 51:377–83. doi: 10.1016/j.pediatrneurol.2014.04.008
50. D'Antonio F, Pagani G, Familiari A, Khalil A, Sagies TL, Malinger G, et al. Outcomes associated with isolated agenesis of the corpus callosum: a meta-analysis. Pediatrics. (2016) 138:e20160445. doi: 10.1542/peds.2016-0445
51. Nam KW, Castellanos N, Simmons A, Froudist-Walsh S, Allin MP, Walshe M, et al. Alterations in cortical thickness development in preterm-born individuals: implications for high-order cognitive functions. Neuroimage. (2015) 115:64–75. doi: 10.1016/j.neuroimage.2015.04.015
52. Woodward LJ, Clark CA, Bora S, Inder TE. Neonatal white matter abnormalities an important predictor of neurocognitive outcome for very preterm children. PLoS ONE. (2012) 7:e51879. doi: 10.1371/journal.pone.0051879
53. Klebermass-Schrehof K, Aumuller S, Goeral K, Vergesslich-Rothschild K, Fuiko R, Brandstetter S, et al. Biometry of the corpus callosum assessed by 3D ultrasound and its correlation to neurodevelopmental outcome in very low birth weight infants. J Perinatol. (2017) 37:448–53. doi: 10.1038/jp.2016.231
54. Gilmore JH, Knickmeyer RC, Gao W. Imaging structural and functional brain development in early childhood. Nat Rev Neurosci. (2018) 19:123–37. doi: 10.1038/nrn.2018.1
55. Hwang SJ, Ji EK, Lee EK, Kim YM, Shin DY, Cheon YH, et al. Gender differences in the corpus callosum of neonates. Neuroreport. (2004) 15:1029–32. doi: 10.1097/00001756-200404290-00019
Keywords: brain, esophageal atresia, long-gap esophageal atresia, LGEA, morphometric, MRI, pediatric, premature
Citation: Kagan MS, Mongerson CRL, Zurakowski D and Bajic D (2022) Impact of Infant Thoracic Non-cardiac Perioperative Critical Care on Homotopic-Like Corpus Callosum and Forebrain Sub-regional Volumes. Front. Pain Res. 3:788903. doi: 10.3389/fpain.2022.788903
Received: 03 October 2021; Accepted: 09 February 2022;
Published: 07 April 2022.
Edited by:
Emma Duerden, Western University, CanadaReviewed by:
Rod Hunt, Monash University, AustraliaLorenzo Fabrizi, University College London, United Kingdom
Copyright © 2022 Kagan, Mongerson, Zurakowski and Bajic. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dusica Bajic, dusica.bajic@childrens.harvard.edu