{'fraction_non_alphanumeric': 0.030705394190871368, 'fraction_numerical': 0.015991062879029683, 'mean_word_length': 4.940652256351915, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 3, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Amplitude of low-frequency fluctuation may be an early predictor of delayed motor development due to neonatal hyperbilirubinemia: a fMRI study

Background Acute bilirubin encephalopathy or kernicterus is the worst consequence of brain damage caused by the elevation of total unbound serum bilirubin (TSB) in neonates. The present study aimed to visualize the characteristic brain regions of neonates with hyperbilirubinemia (HB) using functional magnetic resonance imaging (fMRI) and to measure the amplitude of low-frequency fluctuation (ALFF) values. Methods This was a prospective cohort study, which included newborns with HB who were hospitalized at the Children’s Hospital of Fudan University. The control group included neonates admitted with neonatal simple wet lung or pneumonia without neurological disease or brain injury. Newborns were divided into a severe hyperbilirubinemia group (SHB), moderate HB group, and control group based on TSB levels. The newborns completed routine MRI combined with fMRI scans and brainstem auditory evoked potentials (BAEPs) during their hospitalization. Results A total of 251 newborns were included in this study. There were 45 patients in the SHB group, 65 in the HB group, and 141 in the control group. The average ALFF value in the basal ganglia region in the SHB group was the highest, which was greater than that in the HB and control groups (P<0.001). The ALFF increased with an increase in TSB concentration. Based on the results of the Bayley Scales of infant development assessment, we further found that the most significant difference in ALFF remained in the basal ganglia region between newborns with motor development scores above 70 (including 70) and below 70. Correlation analysis revealed a strong negative correlation between motor development scores and ALFF (r=−0.691, P<0.001). When ALFF alone was used to predict motor development, the sensitivity was 89%. When ALFF was combined with TSB and BEAP results, the area under the ROC curve was the largest (AUC =0.85). The sensitivity and specificity of the model were 67.86% and 90.77%, respectively. Conclusions The ALFF value may be able to serve as an early imaging biomarker and has a greater sensitivity than TSB or BAEP results in predicting long-term motor development (18 m) in HB.


Introduction
Acute bilirubin encephalopathy or kernicterus is the worst consequence of brain damage caused by the elevation of total serum bilirubin (TSB) in neonates (1), which may lead to severe neurological sequelae in the later stage (2)(3)(4). The TSB concentration is markedly elevated (over 20 mg/dL) over a short period of time in newborns, which may result in clinical evidence of brain damage, ranging from subtle neurologic abnormalities to severe encephalopathy or permanent neurologic damage to death (5). Therefore, in neonates with hyperbilirubinemia (HB), the TSB concentration is an important routine indicator to determine the severity of the disease and requires regular monitoring of the concentration.
At present, brain injury caused by a high concentration of TSB usually needs to be confirmed by routine magnetic resonance imaging (MRI). Studies have reported that T1 signal increases appeared in the bilateral regions of the pallidum and sub-thalamic nucleus on routine MRI within the first seven days after birth in neonates with acute bilirubin encephalopathy (6), and T2 hyperintensity was found in kernicterus or chronic bilirubin encephalopathy (7). However, these manifestations are still not considered the gold standard for the diagnosis of bilirubin encephalopathy. Neuropathological studies have shown that large amounts of microglia can be found in basal ganglia areas in vitro and in animal models of kernicterus (8,9). However, the significance of these neuronal proliferations is unclear. Brainstem auditory evoked potentials (BAEPs) (3,(10)(11)(12)(13)(14) and routine MRI provide only qualitative findings and do not have a high specificity in the development of long-term bilirubin encephalopathy or poor neurological prognosis. In the diagnosis of bilirubin encephalopathy, there is an urgent need for quantitative biological parameters that allow rapid assessment of neuroimaging changes at an early stage.
Compared to routine MRI, functional magnetic resonance imaging (fMRI) is also a non-invasive, nonradioactive technique for observing brain activity and has been widely used in neuroscience and cognitive science (15,16). Amplitude of low-frequency fluctuation (ALFF) is a relatively stable indicator in fMRI sequences that can respond to the low-frequency electrical activity of local neurons in the resting state. Most studies using ALFF as a biomarker have focused on specific brain region findings and functional connectivity in patients with psychiatric and neurodegenerative disorders (17)(18)(19)(20)(21)(22)(23). ALFF has the potential to be a special tool in the study of psychiatric or neurological disorders to assess the electrical activity of neurons in the resting state early in the course of the disease (24,25). The present study was designed to visualize the characteristic brain regions of hyperbilirubinemic neonates using fMRI and to measure the ALFF values in specific brain regions in the resting state. We will use ALFF values in combination with serum total bilirubin values and BAEP results to predict neurobehavioral developmental outcomes in neonates at 18 months of age. We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/tp-20-447).

Study design
This was a prospective cohort study. The study selected ALFF as a neonatal exposure factor with the purpose of observing the long-term neurodevelopmental status (18 months old) of neonates with HB at different postmenstrual ages (PMA). The study included newborns with HB who were hospitalized in the Department of Neonates, Children's Hospital of Fudan University between November 2016 and December 2017. These neonates were hospitalized for HB. The control group included neonates admitted with neonatal simple wet lung or pneumonia without HB during hospitalization, neurological disease, or brain injury. The newborns underwent routine MRI and fMRI scans during hospitalization. Newborns with serum total bilirubin above the exchange line were included in the severe hyperbilirubinemia group (SHB); those with serum total bilirubin values above the phototherapy line but below the exchange line were included in the moderate HB group, and the control group consisted of newborns with simple wet lung or pneumonia (control).
The inclusion criteria were as follows: (I) gestational age of 37 weeks or more for full-term babies; (II) diagnosis of neonatal HB based on postnatal age with reference to the bilirubin reference curve of Bhutani et al. (1,(26)(27)(28)(29); (III) the guardian of the newborn signed the informed consent for examination and clinical treatment.
The exclusion criteria were as follows: (I) newborns who were more than 28 days old or who had already been treated in another hospital. (II) Newborns suffering from other neurological diseases: (i) suspected or diagnosed congenital structural malformation of the central nervous system; (ii) the diseases that cause acute significant brain injury, including brain tumors, cerebral hemorrhage, hydrocephalus, hypoxic-ischemic encephalopathy, meningitis, and other central nervous system infections; d. traumatic brain injury. (III) The time from the discovery of HB to completion of the MR scan was longer than seven days. (IV) The data quality does not meet the requirements of research and analysis: (i) sudden change in the condition; (ii) during the scan, the newborn woke up and produced a large number of motor artifacts, which could not be used for reporting and analysis.
The study protocols were registered at www.clinicaltrial. gov (NCT02544100) and were approved by the Medical Ethics Committee of Children's Hospital of Fudan University (2020-247). Written informed consent was obtained from the parents of the infant participants involved in our study. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Clinical baseline and laboratory data
Baseline data included sex, PMA, birth weight, hypertension during pregnancy, hyperglycemia during pregnancy, thyroid dysfunction during pregnancy, parental hepatitis, 5 min Apgar scores, 10 min Apgar scores, routine MR test results (abnormal signals in the basal ganglia or subthalamic region), neonatal hemolytic disease, TSB, and BAEP results. All demographic and clinical information were obtained from the HIS clinical information system.
Neonates with HB were admitted to the hospital for phototherapy, and serum bilirubin levels were monitored. Transcutaneous bilirubin was monitored every 12 h for the first 96 h after phototherapy. If transcutaneous bilirubin exceeded the phototherapy line, we continued phototherapy and monitored the TSB every 24 hours until the TSB dropped below the phototherapy line and did not fluctuate.

Phototherapy treatment
According to the phototherapy reference curve (27,30) and phototherapy method (31), all neonates with HB received standard intensity phototherapy [8-10 μW/(cm 2 ·nm)] when their serum total bilirubin values exceeded the phototherapy threshold. When the serum total bilirubin value reached or exceeded the threshold of blood exchange, they received continuous strong phototherapy for 4-6 h [30 μW/(cm 2 ·nm)], and when the serum total bilirubin decreased to the blood exchange threshold of 50 μmol/L (3 mg/dL), the strong phototherapy mode was adjusted to the standard phototherapy mode. While the neonate was receiving phototherapy, we applied light-shielding eye shields to cover both eyes, covered the perineum with a diaper, and exposed as much skin as possible in other areas. Strong phototherapy was first administered to neonates while preparing for blood exchange, and blood should be exchanged immediately if TSB levels did not decrease or even continued to increase, or for neonates with immune hemolysis if TSB did not decrease to 34-50 μmol/L (2-3 mg/dL) after phototherapy. However, in this study, considering the principle of the fMRI sequence, the effect of blood exchange on functional MRI results is unknown (blood exchange may affect cerebral tissue blood flow), and we did not include neonates who received blood exchange in this study to control for population bias.

Image acquisition
Cranial MR scans of the newborn were performed within seven days of phototherapy (blue light). Specifically, MR examinations were conducted using GE Discovery MR750 3.0 T scanners (GE Healthcare, Milwaukee, WI, USA) with an eight-channel phased-array head coil with foam padding to reduce head motion. The neonate was sedated with phenobarbital (5 mg/kg) half an hour before the cranial MR examination, so that he/she could remain stationary during the examination to reduce artifacts and ensure image quality. Routine MRI included axial and sagittal spin echo T1-weighted sequences [time of repetition (TR) 400 ms, time of echo (TE) 10 ms], axial fast spin echo T2-weighted sequence (TR 3983 ms, TE 85 ms), axial T2 FLAIR sequence (TR 9,000 ms, TE 90 ms), and axial diffusionweighted imaging sequence (TR 5,000 ms, TE 75 ms, B 800 s/mm 2 ).
Subsequently, an echo-planar imaging sequence sensitive to T2* blood oxygenation level dependent (BOLD) contrast material was used to acquire fMRI images with the following settings: TR 2,000 ms; TE 30 ms; voxel size 3.1×3.1×3.0 mm 3 ; flip angle 90°; field of view (FOV) 200 mm; and matrix size 64×64. Each fMRI run included 160 volumes. A minimum of one run (5 min 20 s) was obtained for each infant. During MRI scanning, we strictly controlled the image artifacts of head movement. When the data in this sequence had image artifacts, we rescanned the sequence again until the quality of the data reached the analysis requirements.

The assessment of BAEPs
The recording of BAEPs was carried out with a Nicolet Viking IIIP machine (Nicolet Biomedical, Inc., Madison, WI, USA) at the bedside or in the neuroelectrophysiological room for neonates at the corrected age of two months. Tests were conducted by single-channel differential recording with an active electrode attached to the high forehead and a reference electrode attached to the ipsilateral mastoids. A neonatologist was responsible for BAEP testing, and BAEPs were examined with reference to previous studies that have been reported (12,32,33). The results of BAEPs were interpreted by otolaryngologists and verified by pediatric neurologists. We classified the results as normal or abnormal with brainstem impairment according to the waveform data of BAEPs (12,34). A test result was defined as normal when the repeatable waves V were detectable at less than 40 dB nHL, and latencies and interpeak latencies appeared within 2 SD of the normal for appropriate age groups at our laboratory. Abnormal BAEPs were described as waves V above 40 dB nHL, and latencies of waves I, III, and V exceeded 2 SD of the normal IP of I, III, and V and IPL of I-III, III-V, and I-V (32).

The assessment of the neonatal long-term developmental outcomes
We used the third version of the Bayley Scales of Infant and Toddler Development (BSID-III) to assess the developmental status of the infants at 18 months of age (18 months). The evaluation was divided into three areas: cognition, language, and body movement to systematically evaluate the development level of the subjects (35). It was confirmed that the baby had good sleep and adequate diet before the evaluation. The development evaluations of all subjects were carried out by the same evaluator, who was qualified to conduct BSID assessments. In view of the young age of the subjects, each subject was accompanied by a parent at the time of the evaluation, and the parents were advised not to give any hints to the subjects during the evaluation period. Before the evaluation, the evaluators had a simple communication with the subjects to increase their affinity, which helped to eliminate the tension of infants and give full play to the best level. If three items were passed in a row, the test was continued; otherwise, the previous starting point was selected, and the evaluation was terminated if five items failed in a row. A total of 324 test items were evaluated, including 91 items of cognition, 72 items of gross exercise, 66 items of fine exercise, 48 items of receptive communication, and 47 items of expressive communication. Among them, each item passed to get 1 point, did not pass to get 0 points, the final scores were accumulated to obtain the original score, and we calculated the composite score equivalents according to the scaled scores.

Image processing and extraction of the ALFF value
D a t a p r e p r o c e s s i n g w a s p e r f o r m e d u s i n g S P M 8 (http://www.fil.ion.ecl.ac.uk/spm). ALFF analysis was performed using the in-house software DPARSF (Version 4.3, http:// www.restfmri.net /forum/DPARSF) (36,37). The time course of each voxel was first converted to the frequency domain using a fast Fourier transformation. Then, the square root of the power spectrum was computed and averaged throughout the bandpass (0.01-0.08 Hz) at each voxel, that is, the ALFF. The ALFF of each voxel was divided by the global ALFF value for standardization purposes. After the color was drawn in significantly different brain areas, the standardized ALFF value reflected the degree of its raw ALFF value relative to the average ALFF value of the brain regions with the most significant difference.

Statistical analysis
Analyses were conducted using Stata version 12.0 software (Stata Corp LP, College Station, TX, USA). Differences in demographic variables were analyzed by one-way analysis of variance (ANOVA) for continuous variables and Chi-square test for categorical variables. The imaging statistical group analysis was conducted with a one-way ANOVA test to reveal the difference in the mean ALFF among the groups of three TSB concentrations. Voxel-based comparison of ALFF across the SHB, HB, and control groups was performed using one-way ANOVA in SPM8. All wholebrain voxel-wise analyses were corrected for multiple comparisons to achieve a whole-brain false-positive rate of 0.05. The ALFF difference between the BSID normal group and the BSID abnormal group was compared using the independent samples t-test. Receiver operating characteristic (ROC) curves of TSB, ALFF, and BAEPs were plotted to conduct the area under the curve (AUC) measurements, which were used to predict the outcomes of BSID (18 m) in neonates.

Participants
Based on the inclusion and exclusion criteria, 251 newborns were enrolled in this study (Figure 1), including 45, 65, and 141 patients in the SHB, HB, and control groups, respectively. The mean gestational age was greater than 39 weeks, and the mean birth weight was greater than 3,000 g (P>0.05). There was no difference in the sex composition of the three groups of newborns (P>0.05). Detailed baseline data are shown in Table 1. Of the three groups, six in the SHB group and nine in the HB group had hemolytic syndrome (P=0.001). All neonates with HB underwent a phototherapy intervention. In the SHB group, 16 (35.5%) had abnormal results on conventional MRI; in the HB group, only 5 (7.7%) had abnormal results (P<0.001). In the SHB group, 20 (44.4%) had abnormal BAEP results; in the HB group, 6 had abnormal BAEP

ALFF of over responses to neonates with HB
We compared the ALFF values of fMRI in neonates under three conditions of HB: SHB, HB, and control groups. Using one-way ANOVA analysis, we observed significant increases in ALFF in the bilateral thalamus in neonates with SHB ( Figure 2A). The area at the red-yellow border is the brain region with the most significant ALFF differences between the three groups, shown in the figure (Figure 2A) as the basal ganglia region, and is involved bilaterally. In addition, the more yellow the highlighted area, the greater the ALFF difference. The most significant difference in voxel ALFF values in this region was at the location of the left thalamus. Indeed, the ALFF increased with an increase in TSB concentration under HB. A clear upward trend can be observed in the figure. Pearson correlation analysis showed that TSB was remarkably correlated with ALFF on both sides of the thalamus (r thalamus =0.514, P<0.001, Figure 2B). Among the three groups, the average value of ALFF in the SHB group was the highest, which was greater than that in the HB and control groups (P<0.001, Figure 2C).

Comparison of ALFF values for different findings on routine MRI
We compared the ALFF values of the different findings on routine MRI ( Table 2). In the SHB group, routine MRI showed abnormal signals in the basal ganglia or hypothalamus in 16 neonates, while there were no abnormalities in these local brain regions in the other 29 cases. In the HB group, routine MRI showed abnormal signals in the basal ganglia or hypothalamus in five neonates, while the other 60 showed no abnormalities in these local brain regions. No statistically significant differences were found when comparing the ALFF values between the two groups with normal and abnormal routine MRI findings (P>0.05).

ALFF of over responses to neonates with low motor scores of BSID
A total of 174 newborns completed the follow-up for BSID assessment at 18 months of age. Table 3 shows the neonates' scores on each scale after the BSID assessment. The results showed significant differences in motor development scale scores at 18 months of age among the three groups at different TSB levels (P<0.001). There were no statistically significant differences in cognitive and language development scale scores among the three groups (P>0.05). Based on the results obtained in Table 3, we further compared the differences in ALFF values between newborns with motor development scores above 70 (including 70) and below 70, and analyzed the correlation between motor development scores and ALFF in 18-month-old neonates. By comparing the fMRI images of the two groups, we found that the most significant difference remained in the basal ganglia region ( Figure 3A). The most significant areas of ALFF differences were involved in the basal ganglia region bilaterally, and the left side involved more volume than the right. In addition, the more yellow the highlighted area, the greater the ALFF difference. The most significant difference in voxel ALFF values in this region was also at the location of the left thalamus. Correlation analysis revealed a strong negative correlation between motor development scores and ALFF (r=−0.691, P<0.001) ( Figure 3B). Meanwhile, the level of ALFF was significantly lower in the normal group than in the abnormal group on the BSID assessment ( Figure 3C).

Prediction of the short-term motor scores in neonatal HB
We used four variables, including ALFF, TSB, BAEPs, and routine MRI results of the basal ganglia, to predict the incidence of motor outcomes of BSID using the ROC curve. As shown in Figure 4, when ALFF was combined with TSB, BEAPs, and routine MRI results, the area under the ROC curve was the largest (AUC =0.85). The sensitivity and specificity of the model were 67.86% and 90.77%, respectively. However, when only using the ALFF model for prediction, the sensitivity can increase to 89%. The specific parameters of each ROC curve are shown in Figure 4.

Discussion
Functional MRI is based on the BOLD principle (38,39). Oxygen in the brain is transported to neuronal cells via hemoglobin. When neural activity increases, oxygen consumption also increases, which causes an increase in blood flow to the local area. When hemoglobin  binds to oxygen molecules to form oxyhemoglobin, it is antimagnetic. When oxygen molecules are free from hemoglobin, they are paramagnetic. This magnetic difference then causes a small change in the magnetic resonance signal, which is the basic principle of fMRI.
In the resting state of the human body, researchers have found low-frequency electrical activity in neuronal cells even after excluding physiological noise disturbances such as respiration and heart rate (38). With the continuous development of new fMRI techniques, spontaneous neural activity in the resting state can be captured and analyzed by a number of specific methods, including ALFF, which will The location of significant ALFF difference across two groups under t-test analysis. The red and yellow highlighted area represents the area with the most significant difference in ALFF between the two groups, with the voxel with the most significant difference in ALFF located on the left thalamus. The brighter the area, the more significant the difference. (B) Correlation analysis of motor scores of BSID and ALFF value extracted from the ROI of basal ganglia. The y-axis represents the motor scores of BSID and the x-axis represents the ALFF value. be discussed below. Since Biswal et al. (38) found that spontaneous lowfrequency (nearly 0.01-0.1 Hz) fluctuation (LFF) is highly synchronous among different cortices and deep gray matter, fMRI technology has been developed as a new branch of neuroimaging. ALFF (38) can be used to measure the amplitude of the spontaneous activity of local neurons. The average value of low frequency amplitude (i.e., lowfrequency signal of BOLD) during spontaneous activity is calculated to describe the spontaneous activity of neuronal cells under voxels, which reflects the level of spontaneous activity of each voxel in the resting state from the energy perspective and receives feedback from the regional activity during the resting state.
A 2013 study in JAMA Neurology reported the relationship between ALFF and cognitive outcomes in patients with severe traumatic axonal injury. The results showed that patients had greater ALFF in the frontal regions, which was correlated with cognitive performance (40). Some scholars have concluded that abnormalities in spontaneous brain activity in acute mild traumatic brain injury patients exist in the frontal lobe as well as in distributed brain regions associated with integrative, sensory, and emotional roles, and that abnormal spontaneous neuronal activity in different brain regions could be better detected by ALFF (41). Another study revealed a significant positive correlation between ALFF in the cingulate gyrus and clinical cognitive dysfunction (42) in patients with mild traumatic brain injury. In addition, longitudinal follow-up of the recovery of brain-injured patients revealed significant localized changes in ALFF compared to healthy individuals (43,44). However, there is no relevant application of ALFF in the neonatal brain injury population.
Recent studies in vitro and neonatal animal models have found that a large number of small glial cells around the injured nerve cells are activated, or calcium deposition occurs in neurons, which may be associated with TSB induction (8,9,45). The response of TSB concentration to brain injury in children is indirect, which mainly affects the nervous tissue through complex neuro-metabolic pathways. Therefore, ALFF can be used as a direct and quantitative biomarker to reflect the effect of TSB on the brain in the neonatal population. Therefore, we demonstrated that the occurrence of HB would induce an increase in ALFF in the basal ganglia, and this association would increase with the concentration of TSB. ALFF can be used to observe the effect of elevated TSB at or near the basal ganglia in neonates and quantitatively observe the changes in spontaneous neuronal activity in the gray matter nuclei around the basal ganglia.
MRI is often used clinically to observe brain injury in children with HB and to predict long-term prognosis (46)(47)(48)(49). Due to the deposition of TSB, the increased signal intensity of T2-weighted imaging can be found in the basal ganglia of the affected newborn (50). The injury characteristics of near-term and full-term infants are similar, and injury to the basal ganglia area is the main target of neonatal bilirubin encephalopathy. This conclusion was also confirmed in our study. Thus, the increase in TSB may induce the signal change of basal ganglia area in neuro-images, both the neurotoxicity of unconjugated TSB in the brain and these key brain regions should be considered in the clinical diagnosis.
At present, the T1-or T2-weighted signal changes of the globus pallidus caused by the increase in TSB are the most recognized neuroimaging features in the academic community (51,52). Sometimes, we can see in the routine MRI reports of newborns with HB that there are symmetrical abnormal signals in the pallidum and thalamic nuclei of the brain by T1-or T2-weighted imaging. In clinical practice, however, pediatricians found that neonates with T1 or T2 signal changes in the basal ganglia do not have neurodevelopmental disorders in most cases (7), and this is consistent with the results of the present study. Based on the results in Figure 2, we found that routine MRI can ignore neonates who have no obvious changes in the T1 and T2 signals. Furthermore, the routine sequence cannot quantify brain injury in the basal ganglia. Therefore, fMRI may compensate for the shortcomings of routine MRI sequences in these aspects and quantitatively demonstrate that the significantly affected brain areas were in the basal ganglia, especially the subthalamus.
Some studies have suggested that ALFF is associated with cognitive function following brain injury (24,40,53). It is well known that Bayley scores can be used for the long-term assessment of neonatal neurological behavior and development. Immediately following, we analyzed the difference between the normal and abnormal groups of Bayley motor scores using the ALFF value, and found a very strong difference between these two groups. The functional connectivity of certain brain regions in the injured population is enhanced, and the ALFF values are also increased, which may be related to the compensation of nerve cells (40). In infants with low motor development scores, the changes in ALFF in the basal ganglia area of the brain are most significant. Based on ROC curves, ALFF can be a good predictor of motor development outcomes. When using ALFF combined with TSB and BAEP results to predict motor development outcomes, the AUC and specificity were the largest among all models. Therefore, when we need to synthetically predict the long-term motor development of children, we can use ALFF combined with TSB and BAEP results for better specificity. Compared with TSB or any other variables to predict motor development in neonates with HB, the sensitivity was greatly improved by ALFF. We later found that there were patients with high ALFF who had negative MR or BAEP, as shown in Figure 2. In patients with altered ALFF, routine MR does not necessarily report positive results, and the addition of routine MR or BAEP may instead reduce the true positive rate. This also means that ALFF has a higher sensitivity compared to other examinations and can detect earlier cases reported as negative by routine MR or BAEPs, which should be the significance of ALFF.

Limitations of the study
Since the monitoring of TSB is still not dynamically available in real time, the highest bilirubin values obtained in this study may be underestimated. It is possible that some of the neonates with SHB may have been placed in the milder group. However, we still found significant differences in the ALFF values between the three groups. For the neonates to better perform the MR examination, we used a small amount of phenobarbital prior to the examination to help sedate the neonates and reduce head movement and artifacts. Although we ensured image quality, sedation may affect the baseline levels of quantitative measures such as ALFF in the fMRI sequence. However, even if it was effective, the overall trend of ALFF changes with the change in bilirubin will not change. We considered immature brain development in preterm infants, which may be an important confounding factor affecting ALFF levels; therefore, the study did not include the preterm population. Therefore, the results of the present study cannot be generalized to preterm infants. In addition, the scan time for routine MRI (10 min) plus an fMRI sequence (15 min) takes at least 25 min and will be longer if head movement or artifacts are present. This may not be suitable for newborns with unstable conditions. We hope that more research will focus on this field and use a larger sample to verify our results.

Conclusions
The basal ganglia region in HB has a high ALFF value. ALFF may be an early imaging biomarker and has greater sensitivity than TSB or BAEP results in predicting longterm motor development (18 m) from HB. When combined with TSB and BAEP, results can be obtained with higher specificity to predict long-term motor developmental outcomes in newborns with HB. In the field of pediatrics, especially neonatology, fMRI and ALFF are still in the preliminary stage of research. There is still much research on the understanding and prognosis of bilirubin encephalopathy.