Growth of Very Preterm Infants in a Low-Resourced Rural Setting after Affiliation with a Human Milk Bank.

The extrauterine growth restriction (EUGR) of very preterm infants has been associated with long-term complications and neurodevelopmental problems. EUGR has been reported at higher rates in low resource settings. There is limited research investigating how metropolitan human milk banks contribute to the growth outcomes of very preterm infants cared in rural areas. The setting of this study is located at a rural county in Taiwan and affiliated with the Taiwan Southern Human Milk Bank. Donor human milk was provided through a novel supplemental system. A renewal nutritional protocol was initiated as a quality improvement project after the affiliated program. This study aimed to compare the clinical morbidities and growth outcome at term equivalent age (TEA) of preterm infants less than 33 weeks of gestational age before (Epoch-I, July 2015-June 2018, n = 40) and after the new implementation (Epoch-II, July 2018-December 2020, n = 42). The Epoch-II group significantly increased in bodyweight z-score at TEA ((-0.02 ± 1.00) versus Epoch-I group (-0.84 ± 1.08), p = 0.002). In multivariate regression models, the statistical difference between two epochs in bodyweight z-score changes from birth to TEA was still noted. Modern human milk banks may facilitate the nutritional protocol renewal in rural areas and improve the growth outcomes of very preterm infants cared for. Establishing more distribution sites of milk banks should be encouraged.

1. Introduction

Extrauterine growth restriction (EUGR), defined as an anthropometric measure that is lower than the 10th percentile, is prevalent and occurs in the majority of very preterm infants [1]. Preterm infants with EUGR at discharge have been shown to develop short stature, poor neurodevelopmental outcomes, and cardiovascular and metabolic problems after they grow up [2,3,4,5,6]. EUGR has been linked to intrauterine growth restriction, maternal preeclampsia, low gestational age (GA), feeding intolerance, delay in first enteral feeding, and inadequate nutritional supplement [1,7,8,9,10,11].

To minimize the prevalence of EUGR, standardized nutritional protocols have been proposed in recent years. Early trophic feeding (ETF) and optimal parenteral nutrition (PN) supplements have been recommended to be implemented in every neonatal intensive care unit (NICU) by several studies [12,13,14,15,16,17]. However, in practice, discrepancies between city and rural institutions have been largely overlooked despite the well-known benefits of ETF reported [18,19,20]. In rural hospitals, enteral feeding could be delayed due to the fear of necrotizing enterocolitis (NEC) and the shortage of pediatric surgeons [13,15,21]. The accessibility of a human milk bank (HMB) could be an important factor associated with initiating ETF. The benefits of implementation of an HMB in NICUs have been reported by several studies (Table S1) [22,23,24,25,26]. To date, seldom studies report the influence on growth outcome of these very preterm infants after the implementation of an HMB in NICUs.

Our institution is located at a low-resource area in Taiwan, and the accessibility of donor human milk (DHM) had been difficult for years. However, since July 2018, we have been the first distribution site of the Taiwan Southern Human Milk Bank (TSHMB) (Table S2). Under a novel cooperation model (see methodology Section 2.3), DHM was readily provided without any cost. This cooperation also provides an opportunity to implement a renewal nutrition protocol to fit the recommendations of the ETF guideline of the Taiwan Society of Neonatology [27].

This study primarily aimed to review the growth outcome of very preterm infants at term equivalent age (TEA, post-menstrual age 37 to 42 weeks) before and after the quality improvement (QI) of the renewal nutritional protocol related to the affiliation with an HMB [28]. This study also investigated the effectiveness and safety of the implementation as the secondary outcome.

2. Materials and Methods

2.1. Study Design

We conducted this retrospective study and included all preterm infants less than 33 weeks of GA born between July 2015 and December 2020. Infants with congenital disease or major congenital anomalies involving the cardiovascular or gastrointestinal system were excluded.

Considering the basis of the nutritional protocol before and after the QI project implementation, enrolled infants were separated into two groups for historical comparison: Epoch-I group (July 2015–June 2018) and Epoch-II group (July 2018–December 2020). The study was approved by the Institutional Review Board of MacKay Memorial Hospital (IRB approval number: 21MMHIS079e).

2.2. Study Setting and Affiliation with Human Milk Bank

This study was conducted at a seven-bed NICU at Taitung MacKay Memorial Hospital in Taitung County, Taiwan. Taitung is a rural county located in the east-south of Taiwan Island.

Approximately 70 neonates were treated at the unit yearly, including approximately 15 infants with a very low birth weight. One neonatologist, one resident, and one nurse practitioner were regularly in charge of the patients. The unit was equipped with a milk preparation room.

Since July 2018, this unit began affiliating with TSHMB, which is 200 km from the unit. The first small distribution site of pasteurized DHM was set up. DHM is always ready to be used for newborns whose mothers’ own milk (MOM) is insufficient or unavailable. The families are requested to sign consent for using DHM and the cost is supported by a government project of the Ministry of Health and Welfare, Taiwan. When there is a shortage of DHM in the unit, a new supply will arrive within 24 h after placing an order—orders can be conveniently made through Line®, a social network service (Table S2). The DHM delivery are through a convenient cold chain network operated by a Taiwan private company (SEVEN-ELEVEN ®, Taiwan), which also owned more than 6000 stores by 2021. A shipment usually includes an ice bucket (Dometic® WCI-33, Sweden) with around 50 bottles (100 cc DHM/bottle).

2.3. Nutrition Protocol

The brief comparisons between the two epochs are listed in Appendix A Table A1. The details of the implemented protocol are explained in the following subsection.

Table A1

Comparison of nutritional practices in two epochs.

	Epoch-I(n = 39) 1	Epoch-II(n = 41) 1	p-Value
Age of initial nutrition management, days
PN	1.13 ± 0.74	0.22 ± 0.42	<0.001
First enteral feeding	3.27 ± 3.75	0.02 ± 0.16	<0.001
Milk type at first feeding 2			<0.001
MOM	5 (13.5%)	1 (2.5%)
DHM	0 (0%)	32 (80%)
Preterm formula	32 (86.5%)	7 (17.5%)
Fortification initiating 3
Fortification with exclusive MOM or DHM	5 (13.5%) 2	31 (77.5%) 2	<0.001
Post-menstrual age, weeks	31.8 ± 2.4	30.0 ± 2.0	0.008
Postnatal age, days	19.1 ± 16.1	10.9 ± 4.0	0.003
Length of PN use, days	28.4 ± 14.6	8.0 ± 4.9	<0.001
Feeding amount (with or without fortification) reached 120 mL/kg/day, days	30.1 ± 14.0	10.6 ± 3.7	<0.001
PICC insertion	29 (74.4%)	10 (25.6%)	<0.001

Data are presented as the mean ± standard deviation or number (percentage). PN: parental nutrition; MOM: mother’s own milk; DHM: donor human milk; PICC: peripherally inserted central catheter. 1 One patient died on the same day admitted to NICU in each group was excluded from analysis. 2 Two patients in Epoch-I and one patient in Epoch-II died before starting enteral feeding were 3 Fortifications mean shifting 20 kcal/cc to higher caloric density. Statistical significance was assumed for p < 0.05 (indicated in bold).

2.3.1. Parental Nutrition

The standardized nutritional protocol, modified based on recommendations from the Taiwan Society of Neonatology, started the PN soon after the patient was born, targeting 3–3.5 g per kilogram per day (g/kg/d) of protein by the first day of life (DOL) and increasing to 4 g/kg/d. Lipids were administrated at 1 g/kg/d by first or second DOL and increased to 3 g/kg/d with an increment rate of 0.5–1 g/kg/d. The dextrose concentrations were started at a glucose infusion rate of 4.8 milligrams per kilogram per minute (mg/kg/min) and then increased by daily fluid requirement.

The calcium and phosphorus content in PN was adjusted to the proper 1.3:1 molar ratio to reach the recommendation of the American Academy of Pediatrics.

PN was ceased once the feeding amount reached 120 milliliters per kilogram per day (mL/kg/d) if there was no further indication for intravenous fluid administration.

2.3.2. Enteral Nutrition

The protocol emphasized the use of human milk, including MOM and DHM, and the practice of ETF. Owing to the limited supplement of DHM, patients who were born at less than 30 weeks of gestational age had the top priority to access DHM if MOM was unavailable. Those that were born between 31 to 33 weeks of gestational age whose MOM was unavailable had the next priority. Only those who had unstable clinical condition delayed the initiation of enteral feeding. The volume of ETF was set at 15 mL/kg/d. The length of ETF was set according to the birth weight of the patient. Five-day ETF was set for those with birth weights less than 700 gm, 3-day for those with birth weight between 701 gm and 1000 gm, and 2-day for those with birth weights between 1001 gm and 1250 gm. After ETF was accomplished, the rate of incensement was set at 20 mL/kg/d per day. The feedings were interrupted, or workups were implemented if one of the following was present: (1) bilious or bloody residuals, (2) existing peritoneal signs on physical examination, (3) respiratory problems probably related to feeding. Two-staged fortification was applied, and we fortified human milk with human milk fortifier or concentrated preterm formula [29]. First (22 Kcal/oz) and second (24 Kacl/oz) fortifications were performed after the feeding amount reached 100–120 and 150 mL/kg/d. The increment of feeding amount would be held for 3 days after the first fortification.

2.4. Outcomes and Data Collection

The variables collected from maternal and neonatal electronic chart included perinatal history, the day of nutrition intervention after admission, length of hospital stay, anthropometrical measurement at TEA, medication history, and neonatal morbidities during hospitalization. The GA of a fetus was calculated from the established date of delivery according to the last menstrual period.

The primary outcome data were the Δz-score of body weight (BW) at TEA, the difference of bodyweight z-score between birth and TEA. We calculated the z-scores and percentiles on the basis of reference data on the website [30]. The secondary outcomes were the nutritional practice and the major neonatal morbidity during hospitalization including bronchopulmonary dysplasia (BPD), treated retinopathy of prematurity (ROP), NEC ≥ stage IIA, and late-onset sepsis (LOS) [31,32,33].

2.5. Statistics and Analysis

All analyses were performed using SPSS (Version 26, IBM, Armonk, NY, USA). A p-value of 0.05 was considered significant. We compared continuous variables with Student’s t-test or Mann–Whitney’s U test, whereas we compared categorical data with the chi-square test or Fisher’s exact test, where applicable. The outcome of infants was compared between the study periods utilizing the analysis of covariance. In analysis models, covariates were selected with a priori for clinical consideration.

3. Results

Eighty-four preterm infants were included during the study period, and excluded infants included one extensive Klippel–Trenaunay syndrome and one congenital complex heart disease. (Figure 1). Among these 82 neonates enrolled in the analysis, 40 were in the Epoch-I group and 42 were in the Epoch-II group.

Figure 1

Flow chart of patient enrollment. n: case number.

The Epoch-II group had earlier administration of PN, first enteral feeding, milk fortification, and higher exclusive MOM or DHM feeding before fortification than the Epoch-I group. The length of PN administration was shorter in the Epoch-II group (Appendix A Table A1).

The GA was 28.7 ± 2.4 and 28.1 ± 2.6 (weeks, p = 0.229) and birth weight was 1120 ± 363 and 1132 ± 360 (grams, p = 0.912) in Epoch-I and Epoch-II groups, respectively. The Epoch-I group had lower z-scores of birth weight (−0.32 ± 0.93 vs. 0.07 ± 0.71, p = 0.033), higher pH of 1st blood gas (7.29 ± 0.96 vs. 7.24 ± 0.91, p = 0.008), and higher rate of small for gestational age (20.0% vs. 2.4%, p = 0.007), maternal absent or reverse-end diastolic velocity (20.0% vs. 2.4%, p = 0.007), and usage of umbilical line (20.0% vs. 4.8%, p = 0.035). Otherwise, there was no significant difference between the two groups in baseline demographic characteristics (Table 1).

Table 1

Baseline demographic characteristics.

	n for Analysis	Epoch-I(n = 40)	Epoch-II(n = 42)	p-Value
Perinatal variables
Gestational age, weeks	82	28.7 ± 2.4	28.1 ± 2.6	0.229
Birth weight, grams	82	1120 ± 363	1132 ± 360	0.912
Z-score of birth weight	82	−0.32 ± 0.93	0.07 ± 0.71	0.033
Small for gestational age	82	8 (20.0%)	1 (2.4%)	0.007
Gender: male	82	26 (65.0%)	26 (61.9%)	0.771
Cesarean delivery	82	21 (52.5%)	20 (47.6%)	0.659
AEDV/REDV	82	8 (20.0%)	1 (2.4%)	0.007
Multi-gestational pregnancy	82	4 (10.0%)	10 (23.8%)	0.097
Antenatal steroid	82	34 (85.0%)	36 (85.7%)	0.927
Neonatal variablesAPGAR score
1st minute, median (range)	82	7 (2–9)	7 (1–10)	0.167
5th minute, median (range)	82	9 (5–10)	9 (2–10)	0.111
pH of 1st blood gas analysis	82	7.29 ± 0.96	7.24 ± 0.91	0.008
Surfactant administration	82	16 (40%)	11 (26.2%)	0.183
Umbilical line	82	8 (20%)	2 (4.8%)	0.035
hsPDA requiring surgery	82	7 (17.5%)	3 (7.1%)	0.185
Late onset sepsis	80 1	12 (30.8%)	1 (2.4%)	<0.001
Necrotizing enterocolitis ≥ IIA	80 1	3 (7.7%)	3 (7.3%)	1.000
Treated retinopathy of prematurity	73 2	2 (5.3%)	6 (16.2%)	0.279
Bronchopulmonary dysplasia	73 2	15 (41.6%)	9 (24.3%)	0.115
Death	82	4 (10%)	5 (11.9%)	0.938
Discharge PMA, weeks	70 3	38.5 ± 3.0	38.5 ± 5.1	0.987
Length of hospital stay, days	70 3	67.1 ± 30.0	72.1 ± 47.1	0.600
Z-score of body weight at TEA	68 4	−0.84 ± 1.08	−0.02 ± 1.00	0.002
EUGR	68 4	11 (33%)	4 (11.4%)	0.029
Persistent MOM feeding after discharge (with or without mixing with PDF)	68 4	12 (36.4%)	16 (45.7%)	0.434

Data are presented as the mean ± standard deviation or number (percentage). n: case number; AEDV: absence of end diastolic velocity; REDV: reverse of end diastolic velocity; hsPDA: hemodynamically significant patent ductus arteriosus; PMA: post-menstrual age; TEA: term equivalent age; EUGR: extrauterine growth restriction; MOM: mother’s own milk; PDF: post-discharge formula. 1 One patient died at the same day admitted to NICU in each group was excluded from analysis. 2 Deaths were excluded from analysis. 3 Seven patients in Epoch-I and 5 patients in Epoch-II died or transferred to other hospital were excluded from analysis. 4 Seven patients in Epoch-I and 7 patients in Epoch-II died, transferred to other hospital, or lost to follow up were excluded from analysis. Statistical significance was assumed for p < 0.05 (indicated in bold).

Seven patients in each group were excluded from the growth outcome analysis because of inaccessibility of BW at TEA (Epoch-I: four died, two were transferred to other institutions, one was lost to follow up; Epoch-II: five died, two were lost to follow up). The mean follow-up rate was 82.9%.

Although both groups shared similar post-menstrual ages at discharge and lengths of hospital stay, the Epoch-I group had lower BWs (3217 ± 552 vs. 3572 ± 448, grams, p = 0.006) and z-scores for BW (−0.84 ± 1.08 vs. −0.02 ± 1.00, p = 0.002) at TEA.

Factors which may be associated with the growth outcomes were listed and analyzed (Table 2). Infants with higher gestational age were related to better Δz-score of body weight (mean coefficient = 0.198; 95% confidence interval [CI]: 0.098–0.298, p < 0.001). Neonatal morbidities had significantly negative association with Δz-score of body weight at TEA, including PDA ligation (mean coefficient = −0.957, 95% CI: −1.705–−0.209, p = 0.013) and bronchopulmonary dysplasia (mean coefficient = −0.659, 95% CI: −1.156–−0.162, p = 0.010).

Table 2

Dependence of Δz-score of body weight at term equivalent age on clinical variables: univariate analysis.

	Univariate
	B	95%CI (LB, UB)	p-Value
Gestational age, week	0.198	(0.098, 0.298)	<0.001
Small for gestational age	0.292	(−0.544, 1.129)	0.488
Z-score of birth weight	−0.234	(−0.527, 0.060)	0.117
Sex	0.331	(−0.146, 0.808)	0.171
Cesarean section	0.024	(−0.453, 0.510)	0.920
Antenatal steroid	−0.106	(−0.778, 0.567)	0.755
Multi-gestational pregnancy	0.587	(−0.079, 1.244)	0.079
PDA ligation	−0.957	(−1.705, −0.209)	0.013
Bronchopulmonary dysplasia	−0.659	(−1.156, −0.162)	0.010
Late onset sepsis	−0.580	(−1.217, 0.056)	0.073
Epoch (ref. Epoch I)	0.347	(−0.122, 0.816)	0.145

Linear regression was performed for each independent variable. Statistical significance was assumed for p < 0.05 (indicated in bold). PDA: patent ductus arteriosus; B: mean of coefficients; CI: confidence interval; LB: lower border; UB: upper border.

Although there was no statistical significance between both epochs in the univariate model in regards to the Δz-score of BW (mean coefficient 0.347, 95% CI: −0.122–0.816, p = 0.145) (Table 2), the differences were noted significantly after multivariate linear regression (Table 3, models 1 and 2). In the model 3, epoch had borderline significance in a positive association with Δz-score of BW, which was adjusted with one of neonatal morbidities, PDA ligation (mean coefficient = 0.429, 95% CI: −0.016–0.874, p = 0.058) (Table 3). In the model 3, the p value of epoch was more significant than the p values of PDA ligation, which may hint the effect of epoch may be a more important factor than neonatal morbidities related to Δz-score of BW. Similar results were observed when selecting late onset sepsis or bronchopulmonary dysplasia into multivariate models (Table S5).

Table 3

Dependence of Δz-score of body weight at term equivalent age on clinical variables: multivariate analyses.

	Multivariate Model 1			Multivariate Model 2			Multivariate Model 3
	B	95%CI (LB, UB)	p	B	95%CI (LB, UB)	p	B	95%CI (LB, UB)	p
Gestational age	0.204	(0.102, 0.306)	<0.001	0.200	(0.099, 0.301)	<0.001	0.189	(0.084, 0.295)	0.001
Z-score of birth weight	−0.123	(−0.406, 0.161)	0.392	−0.130	(−0.411, 0.150)	0.356	−0.122	(−0.404, 0.160)	0.389
Sex	0.373	(−0.047, 0.794)	0.081	0.339	(−0.078, 0.756)	0.109	0.320	(−0.101, 0.742)	0.133
Multi-gestational Pregnancy				0.481	(−0.095, 1.057)	0.100	0.456	(−0.126, 1.037)	0.122
PDA ligation							−0.276	(−1.014, 0.461)	0.457
Epoch (ref. Epoch I)	0.501	(0.073, 0.929)	0.022	0.479	(0.056, 0.902)	0.027	0.429	(−0.016, 0.874)	0.058

Multivariate linear regression was performed for each dependent variable. Statistical significance was assumed for p < 0.05 (indicated in bold). PDA: patent ductus arteriosus; B: mean of coefficients; CI: confidence interval; LB: lower border; UB: upper border.

We further tested using absolute BW and z-score of BW at TEA separately as primary growth outcomes. Epoch played significantly important roles in both dependent variables after multivariate linear regression (Tables S3 and S4).

For secondary outcomes, no statistical difference was noted between the two epochs on the aspect of death, NEC ≥ IIA, treated ROP, and BPD. The prevalence of LOS was higher in the Epoch-I group (31% vs. 2.4%, p < 0.001) (Table 1).

4. Discussion

In this QI project coordinated with an HMB, we found a significant decreasing prevalence of EUGR and Δz-score of BW of very preterm infants at TEA after implementing a standardized nutritional protocol with intensive nutritional supplementation and ETF. We also found a markedly decreased prevalence of LOS without an increase in the death rate or NEC for those receiving the new nutritional protocol. To our best knowledge, the present study is the first that discusses the effectiveness and safety of the standardized nutritional protocol that is held in a rural neonatal intensive care unit with a novel utilization of HMB.

4.1. The Influence of DHM Accessibility on Feeding Policy

The increased accessibility of DHM after affiliation with HMB may impact the feeding policy [34]. The first impact may be related to age of initiation of feeding (Appendix A Table A1). During Epoch-II, most of the preterm babies were routinely fed on the day of their birth unless the clinical condition was unstable based on the new nutritional protocol. Only one patient who received postnatal resuscitation due to perinatal asphyxia did not meet the criteria for early enteral feeding. Physicians fed the preterm infants earlier in Epoch-II —3.3 days earlier than Epoch-I. The very preterm infants of Epoch-II reached full enteral feeding earlier—19.5 days earlier than those in Epoch-I. Our findings were consistent with previous studies reporting that the nutritional practices directly related to EUGR included age at first enteral feeding and full enteral feeding, and early feeding is related to early achievement of full enteral feeding [35].

Higher rates of exclusive MOM or DHM feeding before fortification were noted in the Epoch-II group. This resulted from the cooperation between our institute and TSHMB since July 2018. Since then, almost all very low birth weight infants have been fed with their MOM or DHM in our unit. Previous studies reported that using MOM and DHM decreased the incidence of feeding intolerance and the risk of NEC [36]. Although we did not observe a decrease in the rate of NEC between two epochs, the use of MOM and DHM did facilitate our medical staff, including attending physicians, pediatric residents, and nurses, to enforce the new feeding protocol.

4.2. Renewal of Nutritional Protocol, Growth Outcomes, and Neonatal Morbidities

In this study, z-score of BW, Δz-score of BW, and presence of EUGR at TEA reached statistical difference between the two groups. Although the significance of the presence of EUGR disappeared after logistic regression analysis, the statistical differences of z-score of BW, Δz-score of BW at TEA persisted after multivariate linear regression analysis. The results were compatible with previous studies that showed a standardized nutritional protocol could reduce the incidence of EUGR in preterm infants [12,15]. The lower Δz-score of BW at TEA could explain why the patients in Epoch-II had a higher z-score of BW at TEA while comparing with those in Epoch-I.

In addition, we did not interrupt the feeding based on the gastric residual volume, since current evidence showed that routine checking of the gastric residual volume would delay the time of reaching full enteral feeding without benefits for preventing NEC [37]. This may also be a reason why the time to reach full enteral feeding decreased so dramatically after implementing the standardized nutritional protocol. Along with a shorter time to reach full enteral feeding, the length of PN administration became shorter, and the rate of central line usage and LOS became lower.

4.3. The Benefits, Cost, and Challenge of HMB during COVID-19 Pandemic

NEC has been related to high medical cost [38]. Fortunately, HMB has been reported to increase human milk availability in reginal or community NICUs and has been associated with a decrease in the NEC rate [34]. However, those studies seldom reported the growth outcome after the increased availability of human milk from milk banks. Our study has linked the HMB affiliation and a renewal of nutrition protocol. In addition, our study has shown the improvement of growth outcome in very preterm infants without increasing NEC rate. Our beneficial findings may help to promote the establishment of affiliation sites of HMB for the debate of the high cost of DHM [39]. After further analysis (not shown in tables), the cost of DHM in this study was 2–3 US$ per ounce, which was much lower than 3–5 US$ per ounce reported in other studies [39,40]. Importantly, the rate of persistent MOM feeding after discharge non-significant increased from 36.4% to 45.7% after affiliation with the HMB (Table 1), which has been reported in the work of Delfosse [41].

Moreover, it is worth noting that the study period of Epoch-II was during the COVID-19 pandemic. Despite the persistent use of MOM after discharge, which did yield statistical significance between the two groups, we did observe a higher portion of MOM use in Epoch-II after the patient was discharged. It has been reported that DHM supplementation could increase the rate of exclusive breastfeeding [22,42]. Considering the benefits of MOM to long-term neurodevelopmental outcomes, maintaining an adequate supplement of safe DHM during the pandemic became a major issue [43].

4.4. Limitations

This study had several limitations. First, as with many comparisons with historical cohorts, the main limitation of the study is the possibility that changes in management of infants other than HDM may influence the outcomes. For example, the initial PN was started earlier in Epoch-II than in Epoch-I. Under limited case number in this rural hospital, we only could investigate with historical cohorts to explain the effect of implanting “a bundle” of changes after affiliation of a human milk bank. Changes compatible with the implantation were also observed, including shorter duration of total parenteral nutrition, earlier age of initial feeding and full feeding, more infants received fortified human milk not fortified formula, less umbilical lines, less peripheral inserted lines, and less catheter related late onset sepsis.

Second, the small sample size of preterm infants enrolled in the trial was relatively small; we could not control more potential confounders in the multivariate analysis to see if the Δz-score of body weight TEA is associated with the epoch. However, the epoch could explain significant associations with the growth in some results of the many multivariate models we tested. To overcome these disadvantages, a large population-based, multi-institution prospective study should be conducted in the future.

Third, there was a lack of data regarding the exact daily caloric intake and other anthropometrical outcomes, such as head circumference and body length at TEA.

5. Conclusions

A renewed nutrition protocol coordinated with an HMB may improve the growth outcomes of very preterm infants born in rural areas. Future research should continue to explore the effects of distribution sites of the modern HMB in multi-hospitals.