Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Physiol., 22 September 2025

Sec. Gastrointestinal Sciences

Volume 16 - 2025 | https://doi.org/10.3389/fphys.2025.1643482

Development of a refined neonatal rabbit model of necrotizing enterocolitis with relatively longer survival

Qing-Qing Guo,&#x;Qing-Qing Guo1,2Hao Lin,&#x;Hao Lin2,3Zhi-Xin Wang,Zhi-Xin Wang2,3Jian-Bo WenJian-Bo Wen4Juan-Juan Yang
Juan-Juan Yang5*Shi-Yun Lu,
Shi-Yun Lu2,3*
  • 1Department of Intensive Care Unit, The First Affiliated Hospital, Fujian Medical University, Fuzhou, Fujian, China
  • 2Shengli Clinical Medical College, Fujian Medical University, Fuzhou, Fujian, China
  • 3Department of Gastroenterology, Fuzhou University Affiliated Provincial Hospital, Fuzhou, Fujian, China
  • 4Department of Gastroenterology, Afffliated PingXiang Hospital, Gannan Medical University, Pingxiang, Jiangxi, China
  • 5College of Biological Science and Biotechnology, Fuzhou University, Fuzhou, Fujian, China

Objective: This study aimed to develop an animal model that closely replicates the clinical features of necrotizing enterocolitis (NEC), thereby providing a supportive platform for ongoing NEC research.

Methods: 36 two-week-old neonatal suckling rabbits were randomly assigned to three groups: normal group (Group A, n = 12), traditional group (Group B, n = 12), and modified group (Group C, n = 12). Daily monitoring of general condition, post-modeling survival rate, and survival time was conducted in neonatal rabbits. Histopathological assessment of intestinal morphology was performed using hematoxylin and eosin staining, with scoring conducted to evaluate modeling success. Serum concentrations of C-reactive protein (CRP), tumor necrosis factor-alpha (TNF-α), interleukin-10 (IL-10), occludin, zonula occludens-1 (ZO-1), and zonulin were quantified using ELISA. Protein expression levels of occludin and ZO-1 in intestinal tissues were analyzed via Western blot.

Results: Modified group exhibited significantly higher survival rates and prolonged survival time compared to traditional group (p < 0.05). Both model groups demonstrated elevated serum levels of CRP, TNF-α, IL-10, along with reduced concentrations of occluding, ZO-1 and zonulin when compared to normal group (p < 0.05). Similarly, intestinal expression levels of occludin and ZO-1 were significantly lower in both model groups relative to control (p < 0.05). However, no significant differences were observed between model groups in these parameters (p > 0.05).

Conclusion: Our study successfully established a novel neonatal rabbit model of NEC with prolonged survival time by refining conventional modeling methodologies, thereby providing a new platform for investigating NEC pathogenesis and therapeutic interventions.

1 Introduction

Necrotizing enterocolitis (NEC) is a severe gastrointestinal disorder primarily affecting preterm neonates, with incidence rates inversely correlated with gestational age. It is characterized by an exaggerated inflammatory response accompanied by intestinal edema, hemorrhage, and necrosis, often resulting in intestinal perforation. (Neu and Walker, 2011). Common clinical manifestations include vomiting, abdominal distension, diarrhea, and hematochezia. (Sim et al., 2015). Although the underlying pathogenesis of NEC remains incompletely understood, it is widely accepted to be multifactorial. (Sanchez and Kadrofske, 2019). Contributing factors include prematurity, formula feeding, microbial colonization, intestinal inflammation, and hypoxia-ischemia. (Koike et al., 2020).

Currently conservative treatment options, such as probiotics, antibiotics, and supportive care, have shown limited efficacy. (Wang et al., 2019). In cases refractory to medical management, surgical resection of necrotic intestinal segments may be required. However, surgical intervention is associated with high mortality rates, reported to range from 40% to 50%, as well as increased risk of postoperative complications such as short bowel syndrome, neurodevelopmental delay, and sepsis, with an estimated incidence of 40%–70%. (Carr and Gadepalli, 2019). NEC imposes a substantial burden not only on affected neonates but also on families and healthcare systems due to its high morbidity and associated care costs. (Mowitz et al., 2018). Due to the absence of specific and effective therapeutic options, continued in-depth investigation into the pathophysiology and treatment of NEC remains essential.

The development of reliable and clinically relevant animal models is essential for exploring pathological mechanisms and for evaluating potential therapeutic strategies. Conventional animal models of NEC typically involve exposure to hypertonic formula, hypothermia, and cold stress over a period of 3–5 days. (Ravisankar et al., 2018; Sheng et al., 2015). However, several limitations reduce the translational value of these models. First, the timing of NEC onset in these animal models does not correspond with the typical onset window in human neonates, which generally occurs between 14 and 21 days of life. (Gan et al., 2019). Second, the survival duration in conventional models is short, often ranging from only 2–4 days, restricting the feasibility of extended observation and evaluation. (Ravisankar et al., 2018; Sheng et al., 2015; Roy et al., 2018). Third, simultaneous implementation of modeling and therapeutic interventions introduces confounding factors that compromise the validity of experimental outcomes. (Satoh et al., 2016; Li et al., 2017). Finally, the exclusive use of hypertonic formula for nutritional support during modeling does not reflect the complexity of clinical nutritional management in NEC. (Ravisankar et al., 2018; Roy et al., 2018).

In response to these limitations, the present study aimed to establish a modified NEC model in neonatal suckling rabbits that more accurately simulates clinical disease characteristics. Key modifications included enhancement of nutritional support through parenteral nutrition (PN) to prolong survival and facilitate extended observation periods. Two-week-old suckling rabbits were randomly assigned to three groups: a modified model group, a conventional model group, and a standard breastfeeding control group. NEC induction in both model groups involved hypertonic formula, hypothermia, and hypoxic stress, with the modified group additionally receiving PN. The control group was maintained on maternal breastfeeding. Daily monitoring was conducted to assess general health status. Pathological scores, modeling success rates, and survival outcomes were recorded on days 0, 4, 10, and 17. In addition, serum inflammatory markers and indicators of intestinal barrier integrity were evaluated to characterize systemic and localized effects of the induced NEC models.

2 Materials and methods

2.1 Experimental animals and groupings

A total of 39 two-week-old New Zealand white suckling rabbits (both sexes) were obtained from Songlian Experimental Animal Farm, Songjiang District, Shanghai (license number: SCXK [Shanghai] 2017-0008). Experimental grouping was conducted using a segmented randomized design. Based on body mass, the rabbits were stratified into four sections: Section 1 (160–175 g), Section 2 (176–190 g), Section 3 (191–205 g), and Section 4 (206–220 g). After excluding three animals prior to experimentation, 36 suckling rabbits were randomly allocated into three groups (n = 12 per group) as follows: Group A: (normal group): received breastfeeding only; Group B: (traditional group): subjected to a combination of homemade hypertonic formula, hypoxia, and hypothermia; Group C: (modified group): received PN in addition to the same modeling conditions as Group B. All experimental procedures involving animals were reviewed and approved by the Animal Ethics Committee of Pingxiang Hospital, affiliated with Gannan Medical University (Approval No: 2022243).

2.2 Main reagents and instruments

2.2.1 Drugs and nutritional formulations

Abbott Plus Protein Powder (Abbott, United States), Gimlet Puppy Breast Milk Substitute (PetAg, United States), 20% fat emulsion, 5% compounded amino acids, water-soluble vitamins for injection, fat-soluble vitamins for injection, and Andamax (Huarui Pharmaceutical Co.).

2.2.2 Primary reagents

Enzyme-linked immunosorbent assay (ELISA) kits for C-reactive protein (CRP), tumor necrosis factor-alpha (TNF-α), interleukin-10 (IL-10), and zonulin (Jianglai Biotechnology Co., Shanghai, China); ELISA kits for occludin and zonula occludens-1 (ZO-1) (Thermo Fisher Scientific, United States); primary antibodies for occludin and ZO-1 (Abcam, United Kingdom; ab216327 and ab96587, respectively); β-actin antibody (Immunoway Corporation, United States; YM3028); secondary antibodies for occludin and ZO-1 (Kangwei Century Biological Co., Beijing, China); and secondary antibody for β-actin (Solarbio, Beijing, China).

2.2.3 Main instruments

CY-100 portable digital oxygen meter (Meicheng Electrochemical Analysis Instrument Factory, Jiande, Zhejiang Province); refrigerator maintained at 4 °C (Huamei Cold Chain Technology Co., Zhejiang, China); airtight anoxic chamber and custom-built incubator for neonatal rabbits.

2.3 Induction of animal models

Based on findings from preliminary trials in suckling rabbits, the protocol for inducing the conventional NEC model was refined as follows (Zani et al., 2008):

1. The concentration of Abbott Plus Protein Powder in the hypertonic formula was reduced from 200 g/L to 133 g/L.

2. The frequency of daily gavage was decreased to three times per day.

3. All rabbits were provided with maternal breast milk for 3 days prior to induction for adaptive feeding.

In the modified group, NEC was induced using a combination of homemade hypertonic formula, hypothermia, and hypoxia. To enable prolonged survival for observation under sustained catabolic stress, PN was administered concurrently as the sole energy source throughout the induction phase. In the traditional group, NEC induction was achieved using the same protocol of hypertonic formula, hypothermia, and hypoxia, while receiving isovolumetric saline in place of PN from the outset of the experiment. The normal group received maternal breastfeeding only, without any additional interventions beyond equal-volume milk feedings. The induction phase lasted for 3 days in both model groups. From day four onwards, animals in the modified group were sustained with PN as the sole energy source, while those in the traditional group were provided with an isovolumetric isotonic formula (Figure 1). The normal group continued breastfeeding for the duration of the experimental period.

Figure 1
Flowchart showing an experiment with rabbits starting on postnatal day 14. Two groups:

Figure 1. Establishment method of a modified suckling rabbit model of Necrotizing Enterocolitis with relatively longer survival.

Induction procedures were conducted as follows:

1. Gavage feeding: A hypertonic formula prepared with 10 g Abbott Plus Protein Powder and 75 mL Gimlet milk substitute was administered at 15 mL/(kg·time), three times daily, with intervals of at least 4 h. The gavage volume was adjusted according to the animal’s daily weight.

2. Hypothermic and hypoxic stimulation: Hypothermia was induced by placing animals in a refrigerator maintained at 4 °C for 10 min following the first and third round of gavage each day. Subsequently, hypoxia was induced by placing the animals in a chamber filled with a gas mixture of 95% nitrogen and 5% oxygen at a flow rate of 10 L/min for 10 min. These interventions were applied twice daily (Figure 1).

3. PN: In the modified group, central venous access was established via the internal jugular vein 1 day prior to induction. PN was initiated on the day of modeling and maintained throughout the experiment. The PN formulation was adapted from previously published protocols, with a reduced glucose content to account for the energy derived from formula feeding during the induction phase (see Table 1 for formulation details). (Wu et al., 2007; Cai et al., 2013) During the initial three-day induction period, PN was administered at 4 mL/kg/h for approximately 12 h per day, resulting in a daily total of 50 mL/kg and 40 kcal/kg. From day four onwards, the infusion rate was increased to 5 mL/kg/h for over 16 h per day, totaling approximately 100 mL/kg and 80 kcal/kg daily.

4. Volume-matched intervention in the traditional group: During the induction phase, animals in this group received saline in volumes equivalent to those of PN. From day 4 onward, they were provided with an equivalent volume of isotonic formula for nutritional support.

5. Analgesia: To alleviate potential pain and distress associated with the procedures, all animals in the two intervention groups (traditional and modified) received subcutaneous injections of carprofen at a dose of 4 mg/kg once daily throughout the modeling period.

6. Housing conditions: Outside of hypothermic or hypoxic stimulation periods, all animals were housed in incubators maintained at a temperature of 30°C–32 °C and relative humidity of 60%–80%.

Table 1
www.frontiersin.org

Table 1. Detailed formula of nutrient solution for parenteral nutrition.

2.4 General condition monitoring of suckling rabbits

The overall condition of the suckling rabbits was assessed every 6 h from the onset of the modeling phase. Observations included mental status, feeding behavior, activity levels, stool characteristics, abdominal distension, and signs of respiratory distress. Body weight was recorded daily at a consistent time point in the morning. Animals exhibiting moribund conditions were humanely euthanized in accordance with ethical guidelines to facilitate timely specimen collection.

2.5 Timing of specimen collection

Prior to the initiation of modeling, three suckling rabbits were euthanized to obtain baseline blood and small intestinal tissue samples. Subsequently, three animals from each group were euthanized at specified time points: Day 1, Week 1, and Week 2 following the modeling procedure (corresponding to experimental Days 4, 11, and 17, respectively). Blood and intestinal tissue samples were collected by trained personnel using aseptic technique. Samples were transported on ice, promptly aliquoted, and stored at −80 °C. In addition, specimens were collected immediately upon identification of moribund animals to preserve sample integrity.

2.6 Histopathological evaluation of the intestine

2.6.1 Gross examination

Euthanized animals were positioned on a sterile dissection platform, and the abdominal cavity was opened following serum sample collection. The entire gastrointestinal tract, including the stomach and intestinal segments, was excised, rinsed with sterile saline in a Petri dish, and transferred to a clean tray for macroscopic examination. Gross pathological features of the intestines were documented and photographed.

2.6.2 Light microscopic examination

Approximately 2–3 cm of distal ileum tissue (from a blinded intestinal section) was isolated and subjected to hematoxylin and eosin (H&E) staining. Histopathological scoring of the intestinal tissue was performed using the Nadler scale as follows (Nadler et al., 2000):

Score 0: Intact intestinal villi and epithelium with normal histological appearance.

Score 1: Mild separation and swelling in the submucosa or lamina propria.

Score 2: Moderate separation and swelling in the submucosa or lamina propria, with edema in submucosal or muscular layers.

Score 3: Severe separation and swelling in the submucosa or lamina propria, with marked edema in the submucosal or muscular layers, and localized villous loss.

Score 4: Complete villous loss with evidence of intestinal necrosis.

Histopathological evaluations were conducted by two independent investigators under blinded conditions. In cases of significant discrepancy, a third pathologist was consulted to reach a consensus. A final pathology score of ≥2 was defined as indicative of NEC and was used to determine the modeling success rate in each group. Representative pathological images were captured and archived.

2.7 Survival analysis and modeling success rate

Mean survival time and survival rates were calculated at different time points (prior to modeling, during induction, and following the modeling period), based on the survival duration of individual animals. The modeling success rates of the two model groups were determined using histopathological scores, enabling comparison of differences in modeling outcomes between the groups.

2.8 ELISA for serum inflammatory factors and intestinal mucosal injury markers

Whole blood samples were allowed to clot at room temperature for 2 h, followed by centrifugation at 4 °C for 20 min at 2000 rpm. The serum was subsequently collected and stored at −20 °C. Concentrations of CRP, TNF-α, IL-10, occludin, ZO-1, and zonulin were quantified using ELISA kits, according to the manufacturers’ instructions.

2.9 Western blot analysis of small intestinal mucosal injury markers

Total cytoplasmic proteins were extracted from small intestinal tissues stored at −80 °C and quantified using the bicinchoninic acid method. A 50 μg protein sample from each extract was subjected to sodium dodecyl sulfate-polyacrylamide gel electrophoresis, followed by transfer to a polyvinylidene fluoride membrane. Membranes were blocked with 5% skimmed milk prior to antibody incubation.

Primary antibodies were incubated overnight at 4 °C. After washing, membranes were incubated with secondary antibodies at room temperature for 1 h with gentle agitation. Protein expression levels of occludin and ZO-1 were detected using the ECL method. All procedures were carried out in strict accordance with the respective reagent protocols.

2.10 Statistical analysis

All statistical analyses were performed using SPSS software version 22.0 (SPSS Inc., Chicago, IL, United States). Quantitative data are presented as mean ± standard deviation (x̄ ± s). One way analysis of variance or non-parametric tests were employed for intergroup comparisons, while categorical data were analyzed using the chi-squared test. Survival analyses were conducted using the Kaplan-Meier method. A p-value of <0.05 was considered statistically significant.

3 Results

3.1 Alterations in the general status of suckling rabbits

3.1.1 Changes in general condition

Suckling rabbits in the traditional group began exhibiting signs of lethargy, diminished reactivity, and reduced activity approximately 36 h after the initiation of the experiment. After 48 h, these rabbits developed profuse diarrhea, including instances of hematochezia, accompanied by a high mortality rate. In the modified group, signs of lethargy, decreased activity, abdominal bloating, and diarrhea were observed starting on the seventh day after modeling (i.e., the 10th day of the experiment). Most of these rabbits experienced diarrhea, bloody stools, and death between the eighth- and ninth-days following modeling (i.e., the 11th and 12th days of the experiment). Mortality data across the three groups at various time points are presented in Table 2.

Table 2
www.frontiersin.org

Table 2. Mortality in the three groups of suckling rabbits at each time point.

3.1.2 Changes in body mass

Body mass in the normal group increased significantly during the experimental period, whereas the model groups demonstrated significant decreases. Notably, both model groups demonstrated significantly different body mass changes compared to the control group (normal group vs. traditional group: 91.20 ± 22.11 g vs. 69.42 ± 43.67 g; normal group vs. modified group: 91.20 ± 22.11 g vs. 30.58 ± 12.49 g; p < 0.001) (Figure 2a). Additionally, body mass loss was significantly greater in the modified group compared to the traditional group (p < 0.05).

Figure 2
A set of images presenting scientific data: (a) A bar graph showing weight changes in three groups: normal, traditional, and modified. Significant differences are noted with asterisks. (b) Images of intestinal sections from the normal, traditional, and modified groups highlighting morphological differences. (c) Microscopic views of intestinal tissues from the three groups, showing varying levels of structural integrity. (d) A bar graph depicting NEC pathological scores for the groups, with significant differences marked by asterisks.

Figure 2. Alterations in general status and intestinal histopathological changes of suckling rabbits in each group. (a) Changes in body mass. (b) Macroscopic presentations of the intestine. (c) Comparison of pathological sections under light microscope: Sections were stained with hematoxylin and eosin (HE) using standard techniques and photographed at ×10 magnification. (d) NEC pathological scores. Note: *p < 0.05, **p < 0.001.

3.2 Intestinal histopathological changes

3.2.1 Gross specimen observations

In the normal group, the intestines of suckling rabbits exhibited yellowish coloration with good elasticity and no observable signs of congestion or distention. In contrast, the intestines from both model groups demonstrated varying degrees of intestinal dilation, congestion, and discoloration, ranging to a purplish hue. No significant differences in gross intestinal pathology were identified between the two model groups (Figure 2b).

3.2.2 Light microscopic observations

Histological analysis of intestinal tissue in the normal group revealed normal structural features, including intact epithelial linings and villi, orderly glandular architecture, and absence of mucosal, submucosal, or lamina propria congestion, edema, or separation. In contrast, both model groups exhibited extensive histopathological damage. Noted features included villous degeneration, glandular disorganization or loss, tissue fragmentation, necrosis, and detachment of epithelial structures. Edema and separation of the submucosa and lamina propria were also observed, along with thinning or rupture of the muscularis layer (Figure 2c).

3.2.3 Pathological scores

Mean histopathological scores were 0.33 ± 0.49 for the normal group, 2.09 ± 0.66 for the traditional group, and 3.07 ± 1.36 for the modified group. Both model groups exhibited significantly higher scores compared to the normal group (p < 0.001), indicating more severe intestinal injury. Although the modified group demonstrated slightly higher mean scores than the traditional group, the difference was not statistically significant (p > 0.05) (Figure 2d).

3.3 Survival time, survival rate and modeling success rate of suckling rabbits

3.3.1 Survival time

The mean survival time in the normal group was 17.00 ± 0.00 days. In the traditional group, the mean survival time was significantly shorter, at 2.55 ± 1.21 days (p < 0.001 vs. normal). The modified group exhibited an intermediate survival duration, with a mean of 9.83 ± 4.26 days, which was also significantly shorter than the normal group (p < 0.001) but significantly longer than the traditional group (p < 0.001) (Figure 3a).

Figure 3
Bar charts and a survival curve compare normal, traditional, and modified groups in terms of survival time, survival rates at different days, and modeling success rate. The normal group consistently shows higher values across metrics, with significant differences noted by asterisks. The survival curve further illustrates these differences over eighteen days, with the normal group maintaining a 100% survival rate, while traditional and modified groups decrease rapidly.

Figure 3. Comparison of survival time, survival rate and modeling success rate in suckling rabbits. (a) Survival time. (b) 4-day survival rate. (c) 11-day survival rate. (d) 15-day survival rate. (e) Modeling success rates among model groups. (f) Kaplan-Meier survival curve. Note: *p < 0.05, **p < 0.001.

3.3.2 Survival rates

1. At 4 days post-experiment, survival rates were 100% in the normal group, 75.0% in the modified group, and 8.3% in the traditional group. Both model groups exhibited significantly lower survival rates compared to the normal group; however, the difference between the modified group and the normal group was not statistically significant. The modified group had a significantly higher survival rate than the traditional group (p < 0.05) (Figure 3b).

2. By day 11, survival rates were 100% for the normal group, 50.0% for the modified group, and 0% for the traditional group. Both model groups had significantly reduced survival compared to the normal group (p < 0.05), and the survival rate in the modified group was significantly higher than that in the traditional group (p < 0.05) (Figure 3c).

3. At day 15, survival rates were 100% in the normal group, 8.3% in the modified group, and 0% in the traditional group. Both model groups had significantly lower survival rates than the normal group (p < 0.05), but no statistically significant difference was observed between the two model groups (Figure 3d).

3.3.3 Modelling success rate

The modeling success rate was 41.7% in the traditional group and 75.0% in the modified group. Although the modified group exhibited a higher success rate, the difference was not statistically significant (p < 0.001) (Figure 3e).

3.3.4 Survival curves

Survival curve comparisons indicated a significantly lower survival rate in the traditional group relative to both the normal and modified model groups (p < 0.001). Additionally, the modified group demonstrated a significantly lower survival rate than the normal group (p < 0.001) (Figure 3f).

3.4 ELISA for serum inflammatory markers: CRP, TNF-α, and IL-10

3.4.1 CRP comparison

Serum CRP levels were significantly elevated in both the traditional and modified groups compared to the normal group (p < 0.001). Although the CRP level in the modified group was slightly higher than in the traditional group, this difference did not reach statistical significance (p > 0.05) (Figure 4a).

Figure 4
Graphs and a Western blot image depict the expression levels of various proteins and markers in three groups: normal, traditional, and modified. Panels (a) to (f) show protein levels of CRP, IL-10, TNF-α, Occludin, ZO-1, and Zonulin, respectively, with significant differences indicated by asterisks. Panels (g) and (h) present relative expression of Occludin and ZO-1 normalized to β-actin. Panel (i) shows a Western blot for ZO-1, Occludin, and β-actin. Statistical significance is denoted as * for p<0.05, ** for p<0.01, with

Figure 4. Comparison of serum inflammatory factors and intestinal mucosal damage markers among all groups. (a) Serum CRP level. (b) Serum IL-10 level. (c) Serum TNF-α levels. (d) Serum occludin level. (e) Serum ZO-1 level. (F) Serum zonulin protein level. (g) Intestinal occludin level. (h) Intestinal ZO-1 level. (i) Comparison of intestinal tight junction protein expression. Note: *p < 0.05, **p < 0.001.

3.4.2 IL-10 comparison

IL-10 expression was significantly increased in both model groups relative to the normal group, with statistically significant differences observed (p < 0.001). While IL-10 levels were slightly higher in the modified group than in the traditional group, the difference was not statistically significant (p > 0.05) (Figure 4b).

3.4.3 TNF-α comparison

Serum TNF-α levels were significantly elevated in both model groups compared to the normal group (p < 0.001). A slightly higher TNF-α concentration was observed in the modified group compared to the traditional group, but the difference lacked statistical significance (p > 0.05) (Figure 4c).

3.5 ELISA for serum occludin, ZO-1, and zonulin

3.5.1 Occludin comparison

Occludin expression levels were significantly reduced in both model groups compared to the normal group (p < 0.001). Although the occludin level in the modified group was higher than that in the traditional group, the difference was not statistically significant (p > 0.05) (Figure 4d).

3.5.2 ZO-1 comparison

ZO-1 expression levels were significantly decreased in both model groups relative to the normal group, with statistical significance (p < 0.001). The modified group demonstrated a slightly higher ZO-1 level compared to the traditional group, but this difference did not reach statistical significance (p > 0.05) (Figure 4e).

3.5.3 Zonulin protein comparison

Zonulin expression levels were significantly elevated in both model groups compared to the control group (p < 0.05). The zonulin level in the modified model group was slightly higher than that in the general model group, although this difference was not statistically significant (p > 0.05) (Figure 4f).

3.6 Western blot assays for occludin and ZO-1

3.6.1 Occludin comparison

Occludin expression was significantly reduced in both model groups compared to the normal group, with no significant difference observed between the two model groups (Figures 4g,i).

3.6.2 ZO-1 comparison

Similarly, ZO-1 expression was significantly decreased in both model groups relative to the normal group, with no significant difference observed between the two model groups (Figures 4h,i).

4 Discussion

The findings of this study indicate that, relative to the conventional modeling approach, the modified NEC modeling method applied to suckling rabbits confers distinct advantages, including prolonged survival time and disease characteristics that more closely resemble the pathogenesis observed in clinical NEC. Consequently, the modified method represents a superior approach for generating animal models of NEC.

Due to the absence of specific clinical symptoms and physical signs, histopathological changes in the intestine remain the definitive standard for clinical diagnosis of NEC. (Shao, 2019). These changes primarily encompass neutrophil infiltration in the intestinal wall, disruption of villous structures, submucosal edema, and hemorrhage. (Papillon et al., 2013). Furthermore, pathological evaluation constitutes the gold standard for confirming successful NEC model establishment. (Nadler et al., 2000). In the present study, both model groups demonstrated marked disruption of intestinal tissue architecture, including villous destruction and submucosal congestion and edema, as observed under light microscopy. These histopathological alterations conformed to the NADLER scoring criteria for NEC diagnosis and were consistent with pathological findings reported in pediatric NEC cases.

In addition to microscopic pathology, gross morphological changes observed in intestines affected by NEC in pediatric patients include dilatation, thinning of the intestinal wall, congestion with a purple-black discoloration, necrosis, and perforation. (Ballanc et al., 1990). The gross intestinal alterations identified in the model groups in this study predominantly involved diminished elasticity, severe dilatation, congestion, and necrosis, closely mirroring those seen in clinical NEC.

The maintenance of intestinal mucosal barrier integrity relies on tight junction proteins such as ZO-1 and occludin, which are critical for epithelial barrier function. (Slifer and Blikslager, 2020). Zonulin, a human protein that modulates intestinal permeability through interactions with tight junction proteins, serves as a valuable biomarker for intestinal mucosal injury. (Barbaro et al., 2020).

Ravisankar et al. reported that disruption of tight junction proteins compromises gut barrier function, permitting translocation of microbial products and toxins into systemic circulation, thereby contributing to the pathogenesis of diseases including NEC and sepsis. (Ravisankar et al., 2018). Similarly, Bein et al. observed that expressions of occludin and ZO-1 are significantly reduced in neonates with NEC compared to healthy controls, implicating increased intestinal permeability and gut tissue necrosis in affected patients. (Bein et al., 2018). In alignment with these observations, the current study found significantly elevated serum zonulin levels along with decreased serum and intestinal expressions of tight junction proteins in both model groups relative to controls. These results collectively indicate pronounced intestinal mucosal injury, necrosis, and increased permeability in the modeled suckling rabbits.

It is now recognized that NEC is associated with increased localized expression of various inflammatory factors, including interleukin-6 (IL-6), interleukin-8 (IL-8), IL-10, transforming growth factor, and TNF, which serve as nonspecific markers for NEC. (Cho et al., 2016). In the present study, serum levels of inflammatory factors—CRP, TNF-α, and IL-10—were significantly elevated in the model groups compared to controls, indirectly confirming the successful establishment of the NEC model.

The construction of reliable animal models constitutes the foundation of experimental research; consequently, the validity of related investigations depends on whether these models accurately replicate clinical disease characteristics such as age of onset, survival time, and the capacity to simulate clinical treatment processes. Clinically, most neonates with NEC develop the disease between 11 and 24 days of age. (Gan et al., 2019). However, conventional NEC models often utilize neonatal rats that are sacrificed after only 2–3 days of induction, which does not align with the typical age of disease onset in human neonates. (Ravisankar et al., 2018; Sheng et al., 2015).

In contrast, this study employed 2-week-old suckling rabbits, with disease induction occurring after 3 days, thereby more closely reflecting the clinical timing of NEC onset. The primary clinical manifestations of NEC in neonates include hematochezia, abdominal distension, vomiting, and hypoactive bowel sounds, which correspond to the symptoms observed in the experimental model groups. (Shao, 2019). The current clinical mortality rate for neonates with NEC within 5 days of admission is approximately 12%, and survival rates for Bell stage I and II patients can exceed 80% with conservative medical management. (Shao, 2019; Luig et al., 2005). A significant increase in both survival time and survival rate in the modified model group was demonstrated in this study, thereby approximating the survival outcomes observed in clinical pediatric populations.

Many existing experimental NEC studies initiate treatment interventions concurrently with disease induction, a practice that may compromise experimental accuracy. (Satoh et al., 2016; Li et al., 2017). In contrast, the extended survival time and increased survival rate observed in the modified model group allowed for the administration of interventions only after the modeling phase, minimizing potential confounding effects and enhancing the precision of experimental results.

Clinically, neonates with NEC often require intensive maintenance therapy, including empirical antibiotic treatment, gastrointestinal decompression, PN, fluid resuscitation, correction of hematological and metabolic imbalances, and support for circulatory and respiratory functions. (Knell et al., 2019). While experimental constraints preclude full replication of these complex clinical treatment conditions, this study enhanced the energy supply to suckling rabbits through simulated clinical PN, which contributed to correcting water-electrolyte imbalances and thus more closely approximated the clinical management of NEC.

This improved experimental protocol may have broader applicability beyond the construction of a prolonged NEC model. Animal models of severe acute pancreatitis, severe ulcerative colitis, and severe cholangitis are often associated with high mortality rates due to extensive fluid loss, elevated energy demands, and complications such as infectious shock. (Antoniou et al., 2012; Xu et al., 2010; Hu et al., 2014). Enhancing energy supply and facilitating fluid resuscitation through PN via jugular vein catheterization during the modeling period may prolong survival, reduce mortality, and provide an extended window for subsequent experimental interventions. Additionally, the potential application of this modified modeling approach to other inflammatory bowel disease models, including ulcerative colitis and Crohn’s disease, warrants further investigation and represents a future research focus of the present study team.

Nevertheless, this study has several limitations. First, the relatively small sample size in this study may have influenced the reliability of the results. Second, while the provision of supportive care via the central venous access is the most plausible explanation for the prolonged survival, the precise physiological and molecular mechanisms behind this improvement remain unexplored in this study. Finally, it is important to acknowledge that all current NEC models, including the modified one presented here, are developed based on incomplete understanding of the disease’s exact etiology and pathophysiology. Therefore, future studies are necessary not only to validate these findings with larger sample sizes but also to utilize this model to investigate the specific mechanisms of survival improvement and to continue refining modeling techniques as our understanding of NEC deepens.

In conclusion, beyond prolonged survival, our modified NEC model offers several key advantages: (1) Enhanced clinical relevance by more closely mimicking the human disease course; (2) A prolonged therapeutic window that facilitates the investigation of disease mechanisms and the evaluation of potential treatments; (3) Elimination of confounding factors by temporally separating the disease induction phase from subsequent experimental interventions, thereby increasing the accuracy of experimental outcomes; and (4) The versatility of the central venous access, which not only provides essential nutritional and fluid support but also serves as a platform for administering various therapeutics (e.g., antibiotics, anticoagulants), making the experimental process more aligned with clinical practice. Notably, this ‘life line’ platform is applicable beyond NEC research to other severe disease models such as severe acute pancreatitis and sepsis.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The animal study was approved by ethics Committee of Affiliated PingXiang Hospital (No. 2022243). The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

Q-QG: Project administration, Formal Analysis, Writing – original draft, Data curation, Conceptualization. HL: Data curation, Project administration, Conceptualization, Writing – original draft, Methodology. Z-XW: Writing – original draft, Visualization, Data curation, Resources, Project administration, Validation. J-BW: Conceptualization, Formal Analysis, Methodology, Project administration, Writing – original draft. J-JY: Writing – review and editing, Supervision, Software, Investigation, Validation. S-YL: Writing – review and editing, Conceptualization, Methodology, Project administration, Formal Analysis.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. Natural Science Fundation of Fujian Province (No.2023J01557). Natural Science Fundation of Fujian Province (No.2022J01560). The Leading Project Foundation of Science and Technology, Fujian Province (No. 2023Y0053). Joint Funds for Innovation of Science and Technology, Fujian Province (No. 2023Y9325).

Acknowledgments

We would like to acknowledge the hard and dedicated work of all the staff that implemented the intervention and evaluation components of the study.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Abbreviations

NEC, necrotizing enterocolitis; PN, parenteral nutrition; HE, hematoxylin and eosin; CRP, C-reactive protein; TNF-α, Tumor Necrosis Factor -α; IL-10, Interleukin-10; ZO-1, Zonula Occludens-1; BCA, Bicinchoninic Acid; PVDF, Polyvinylidene fluoride; ECL, Enhanced Chemiluminescence; ANOVA, analysis of variance; IL-6, Interleukin-6; IL-8, Interleukin-8; TGF, Transforming Growth Factor; TNF, Tumor Necrosis Factor.

References

Antoniou E., Margonis G. A., Angelou A., Pikouli A., Argiri P., Karavokyros I., et al. (2012)2016). The TNBS-induced colitis animal model: an overview. Ann. Med. Surg. 11 (11), 9–15. doi:10.1016/j.amsu.2016.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Ballance W. A., Dahms B. B., Shenker N., Kliegman R. M. (1990). Pathology of neonatal necrotizing enterocolitis: a ten-year experience. J. Pediatr. 117 (1 Pt 2), 6–13. doi:10.1016/s0022-3476(05)81124-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbaro M. R., Cremon C., Morselli-Labate A. M., Di Sabatino A., Giuffrida P., Corazza G. R., et al. (2020). Serum zonulin and its diagnostic performance in non-coeliac gluten sensitivity. Gut 69 (11), 1966–1974. doi:10.1136/gutjnl-2019-319281

PubMed Abstract | CrossRef Full Text | Google Scholar

Bein A., Eventov-Friedman S., Arbell D., Schwartz B. (2018). Intestinal tight junctions are severely altered in NEC preterm neonates. Pediatr. Neonatol. 59 (5), 464–473. doi:10.1016/j.pedneo.2017.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai W., Tang Q. Y., Wang Y., Feng Y., Wu J., Qian L. X., et al. (2013). Guidelines for the clinical use of neonatal nutrition support in China. J. Clin. Pediatr. 30 (12), 1177–1182.

Google Scholar

Carr B. D., Gadepalli S. K. (2019). Does surgical management alter outcome in necrotizing enterocolitis? Clin. perinatology 46 (1), 89–100. doi:10.1016/j.clp.2018.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho S. X., Berger P. J., Nold-Petry C. A., Nold M. F. (2016). The immunological landscape in necrotising enterocolitis. Expert Rev. Mol. Med. 18, e12. doi:10.1017/erm.2016.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Gan X., Mao J., Li J. (2019). Clinical features and outcomes of neonatal necrotizing enterocolitis. Chin. J. Contemp. Pediatr. 21 (10), 1028–1032. doi:10.7499/j.issn.1008-8830.2019.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu Z. J., Jia C. K., Fan Z. N. (2014). Establishment of pig models of acute cholangitis of severe type under endoscope. J. Chin. Pract. Diagnosis Ther. 28 (7), 652–654.

Google Scholar

Knell J., Han S. M., Jaksic T., Modi B. P. (2019). Current status of necrotizing enterocolitis. Curr. problems Surg. 56 (1), 11–38. doi:10.1067/j.cpsurg.2018.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Koike Y., Li B., Ganji N., Zhu H., Miyake H., Chen Y., et al. (2020). Remote ischemic conditioning counteracts the intestinal damage of necrotizing enterocolitis by improving intestinal microcirculation. Nat. Commun. 11 (1), 4950–50. doi:10.1038/s41467-020-18750-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li X. M., Li X. W., Shang Q. J., Gao Z., Hao F., Guo H., et al. (2017). Fecal microbiota transplantation (FMT) could reverse the severity of experimental necrotizing enterocolitis (NEC) via oxidative stress modulation. Free Radic. Biol. and Med. 108, 32–43. doi:10.1016/j.freeradbiomed.2017.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Luig M., Lui K.NSW and ANG (2005). Epidemiology of necrotizing enterocolitis--Part I: changing regional trends in extremely preterm infants over 14 years. J. Paediatr. child health 41 (4), 169–173. doi:10.1111/j.1440-1754.2005.00582.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mowitz M. E., Dukhovny D., Zupancic J. A. F. (2018). The cost of necrotizing enterocolitis in premature infants. Seminars fetal and neonatal Med. 23 (6), 416–419. doi:10.1016/j.siny.2018.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nadler E. P., Dickinson E., Knisely A., Zhang X. R., Boyle P., Beer-Stolz D., et al. (2000). Expression of inducible nitric oxide synthase and interleukin-12 in experimental necrotizing enterocolitis. J. Surg. Res. 92 (1), 71–77. doi:10.1006/jsre.2000.5877

PubMed Abstract | CrossRef Full Text | Google Scholar

Neu J., Walker W. A. (2011). Necrotizing enterocolitis. N. Engl. J. Med. 364 (3), 255–264. doi:10.1056/NEJMra1005408

PubMed Abstract | CrossRef Full Text | Google Scholar

Papillon S., Castle S. L., Gayer C. P., Ford H. R. (2013). Necrotizing enterocolitis: contemporary management and outcomes. Adv. Pediatr. 60 (1), 263–279. doi:10.1016/j.yapd.2013.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravisankar S., Tatum R., Garg P. M., Herco M., Shekhawat P. S., Chen Y. H. (2018). Necrotizing enterocolitis leads to disruption of tight junctions and increase in gut permeability in a mouse model. BMC Pediatr. 18 (1), 372. doi:10.1186/s12887-018-1346-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy S. K., Meng Q. H., Sadowitz B. D., Kollisch-Singule M., Yepuri N., Satalin J., et al. (2018). Enteral administration of bacteria fermented formula in newborn piglets: a high fidelity model for necrotizing enterocolitis (NEC). PloS one 13 (7), e0201172. doi:10.1371/journal.pone.0201172

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez J. B., Kadrofske M. (2019). Necrotizing enterocolitis. Neurogastroenterol. Motil. 31 (3), e13569. doi:10.1111/nmo.13569

PubMed Abstract | CrossRef Full Text | Google Scholar

Satoh T., Izumi H., Iwabuchi N., Odamaki T., Namba K., Abe F., et al. (2016). Bifidobacterium breve prevents necrotising enterocolitis by suppressing inflammatory responses in a preterm rat model. Benef. microbes 7 (1), 75–82. doi:10.3920/BM2015.0035

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao X. M. (2019). Practical neonatology. 5th Edition. Beijing: People's Medical Publishing House.

Google Scholar

Sheng Q. F., Lv Z. B., Cai W., Song H. L., Qian L. X., Huaibo Mu H. B., et al. (2015). Effects of recombinant human beta-defensin-3 in neonatal rats with necrotizing enterocolitis. Chin. J. Pediatr. Surg. 36 (2), 109–112. doi:10.1038/pr.2014.93

CrossRef Full Text | Google Scholar

Sim K., Shaw A. G., Randell P., Cox M. J., McClure Z. E., Li M. S., et al. (2015). Dysbiosis anticipating necrotizing enterocolitis in very premature infants. Clin. Infect. Dis. 60 (3), 389–397. doi:10.1093/cid/ciu822

PubMed Abstract | CrossRef Full Text | Google Scholar

Slifer Z. M., Blikslager A. T. (2020). The integral role of tight junction proteins in the repair of injured intestinal epithelium. Int. J. Mol. Sci. 21 (3), 972. doi:10.3390/ijms21030972

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang K. W., Tao G. Z., Sylvester K. G. (2019). Recent advances in prevention and therapies for clinical or experimental necrotizing enterocolitis. Dig. Dis. Sci. 64 (11), 3078–3085. doi:10.1007/s10620-019-05618-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu J., Hong L., Cai W., Tang Q., Shi C. (2007). Glutamine attenuates TPN-associated liver injury in infant rabbits. Eur. J. Pediatr. 166 (6), 601–606. doi:10.1007/s00431-006-0294-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu J. J., Liu J., Hou X. H. (2010). Relationship between different moulding methods and intestinal bacterial translocation in an animal model of acute pancreatitis. Chin. J. Pancreatol. 10 (3), 225–226.

Google Scholar

Zani A., Cordischi L., Cananzi M., De Coppi P., Smith V. V., Eaton S., et al. (2008). Assessment of a neonatal rat model of necrotizing enterocolitis. Z. fur Kinderchir. 18 (6), 423–426. doi:10.1055/s-2008-1038951

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: animal model, intestinal tissue, necrotizing enterocolitis, New Zealand rabbits, protein expression

Citation: Guo Q-Q, Lin H, Wang Z-X, Wen J-B, Yang J-J and Lu S-Y (2025) Development of a refined neonatal rabbit model of necrotizing enterocolitis with relatively longer survival. Front. Physiol. 16:1643482. doi: 10.3389/fphys.2025.1643482

Received: 09 June 2025; Accepted: 05 September 2025;
Published: 22 September 2025.

Edited by:

Minesh Khashu, University Hospitals Dorset NHS Foundation Trust, United Kingdom

Reviewed by:

Amilcare Barca, University of Salento, Italy
Dhirendra Kumar Singh, University of North Carolina at Chapel Hill, United States

Copyright © 2025 Guo, Lin, Wang, Wen, Yang and Lu. 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: Shi-yun Lu, bHVzaGl5dW5AZmptdS5lZHUuY24=; Juan-juan Yang, anVhbnlhbmdAZnp1LmVkdS5jbg==

These authors have contributed equally to this work

Disclaimer: 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.