Skip to main content

Oral administration of Moringa oleifera leaf powder relieves oxidative stress, modulates mucosal immune response and cecal microbiota after exposure to heat stress in New Zealand White rabbits



Heat stress (HS) disrupts the gut barrier allowing the uptake of lipopolysaccharide (LPS) and leads to an inflammatory response and changes in gut microbiota composition. Moringa oleifera leaf powder (MOLP) has been proposed to combat HS, yet its alleviate role is currently under investigation. The current study investigated the effects of chronic HS and MOLP supplementation on changes in redox status and immune response of cecal mucosa along with alteration in cecal microbiota.


A total of 21 young New Zealand White (NZW) rabbits (male) about 32 weeks old (mean body weight of 3318 ± 171 g) reared on a commercial pelleted diet were employed; divided into three groups (n = 7): control (CON, 25 °C), heat stress (HS, 35 °C for 7 h daily), and HS supplemented orally with MOLP (HSM, 35 °C) at 200 mg/kg body weight per day for 4 weeks.


The results demonstrated that MOLP supplementation increased organ index of cecal tissue compared with the HS group (P > 0.05). Levels of malonaldehyde (MDA) and activity of superoxide dismutase (SOD) as well as lactate dehydrogenase (LDH) were reduced in the cecal mucosa of the HSM group compared with the HS group. MOLP downregulated the contents of cecal mucosa LPS, several inflammatory markers (TNF-α/IL-1α/IL-1β), and myeloperoxidase (MPO) in the HSM group (P < 0.05). Secretory immunoglobulin A (SIgA) was increased in the HSM group compared with the HS group (P < 0.05). The transcriptome of cecal mucosa showed that MOLP reduced gene expression relative to several immune factors, including IL-10, IFNG, and RLA, whereas both HS and MOLP increased the gene expression of fat digestion and absorption pathway, including APOA1, FABP1, FABP2, MTTP, and LOC100344166, compared to the CON group (P < 0.001). At the phylum level, the relative abundance of Proteobacteria was increased by HS, while Actinobacteria was significantly increased by HSM compared to other groups (P < 0.05). At genus level, Papillibacter was higher in abundance in HSM groups compared to CON and HS groups (P < 0.05). Higher butyrate concentrations were observed in the HSM group than HS and CON groups (P < 0.05).


In conclusion, HS in growing rabbits resulted in alteration of cecal microbiota at phyla level as well as increased oxidative stress and expression of mucosal inflammatory genes. Whereas, oral MOLP supplementation elevated the relative weight of cecum, affected their immunological and cecal micro-ecosystem function by improving antioxidant status and down-regulating mucosal tissue inflammatory response.


The cecum of rabbits is proportionally the largest of any mammal. It is twice the length of the abdominal cavity and comprises 40–60% of the total volume of the gastrointestinal tract [1]. Moreover, 40% of the gastrointestinal digesta of rabbits is contained in its cecum [2], where the digesta is retained, mixed, and fermented by microbes [3]. As a cecum fermenter, rabbits possess a special re-ingestion system called cecotrophy. This system provides up to 30% of the rabbits’ daily N intake, which is mostly derived from cecal microbes [4, 5]. In addition, the short chain fatty acids (SCFAs) (acetic, formic, propionic, and butyric acids) derived from microbial fermentation in the cecum are actively absorbed through the cecal and colonic walls and utilized by the rabbit as energy sources [1]. Therefore, it is important for rabbits to create and sustain an optimal gut environment regardless of external environmental challenges.

Heat stress (HS) is one of the largest environmental impediments to efficient animal agriculture—especially rabbit husbandry, due to rabbits’ lack of functional sweat glands [6]. Although likely multifactorial, the nexus for many negative consequences of HS might be mediated by hyperthermia’s adverse effects on intestinal epithelial barrier [7]. Heat disrupts the tight junctions between enterocytes lining the intestinal tract and allows the uptake of intestinal pathogens and molecules, such as lipopolysaccharide (LPS) [8]. LPS is a Gram-negative bacterial cell wall component that can stimulate the immune system and, in large doses, cause harmful systemic inflammation [8]. HS suppresses different components of the immune system, thereby enhancing the susceptibility of an animal to various diseases and negatively affecting intestinal mucosa and microbiota composition [9, 10].

Feeding Moringa oleifera under normal conditions has been reported to improve the growth performance, FCR and health status in several livestock species when given at 5% or less of total dry matter intake [11]. Previously, we have also demonstrated the ability of Moringa oleifera leaf powder (MOLP) to improve growth performance, prevent intestinal morphological atrophy effectively, relieve oxidative damage, and downregulate mRNA expression of mucosal inflammatory genes in jejunum of the rabbits under HS conditions [12]. However, whether MOLP can contribute to homeostasis under HS in the cecum of New Zealand White (NZW) rabbits by modulating microbiome and transcriptome response is unknown. Therefore, the present study aimed to evaluate the potential role of MOLP under heat stress; and to explore the link between HS and the gut microbiome and association of HS with antioxidant and immune activity of cecal mucosa by using a high-throughput 16S rRNA sequence and transcriptome analysis.


Experimental design, animals, and housing

A total of 21 male NZW rabbits (32 ± 1 weeks of age) with similar body weight (3318 ± 171 g) were randomly assigned to each of the three treatment groups: control (CON), heat stress (HS), and HS with MOLP supplementation (HSM). Each group had seven replicates/rabbits and each rabbit was retained in individual cage. The nutrient composition of packaged MOLP obtained from commercial sources (Xishuangbanna, Yunnan Province, China) has been given (Supplementary Table 1). All the 21 rabbits were supplied with the pelleted feed ad libitum and fresh water round the clock with exception that HSM group received oral MOLP supplementation administered as gavage syringe at 200 mg/kg body weight once daily in the morning throughout the 28-day experiment. All rabbits in the CON group were maintained at 25 °C throughout the experiment. However, rabbits in the HS and HSM groups underwent thermal stress in thermostatically controlled rooms at 35 ± 1 °C for 7 h daily (09:00–16:00). The temperature of the heat stress houses was brought down to 25 °C for the remaining 17 h of the day after the 7 h of heat wave. Body weight of each rabbit was recorded twice a week and at the end of experiment for the determination of average daily weight gain (ADG). The feed consumption per cage was recorded daily to calculate average daily feed intake (ADFI). Later FCR was calculated by using data of ADG and ADFI.

Sample collection and determination of cecal organ index

The samples were collected towards the end of trial. At day 29, final weighing of rabbits was done after fasting for 6 h. All 21 rabbits were slaughtered between 7:00 and 10:00 in a blinded procedure. The rabbits were humanely slaughtered according to the method described by Nakyinsige et al. [13] by severing carotid artery, jugular vein, trachea, esophagus, and vagus nerve. Their abdomens were incised and after immediate removal of cecum, the weights were measured. Organ indexes were expressed relative to body weight (g of organ per kg of body weight). Next, the cecal digesta samples and the cecal mucosa samples (after digesta removed and tissue rinsed with 0.9% NaCl) were collected into 2 mL Eppendorf tubes and frozen immediately in liquid nitrogen (− 80 °C). The mucosa tissue samples remained in liquid nitrogen while the cecal digesta samples were later kept at − 20 °C until further analysis.

Oxidative stress level in cecal mucosa

Catalase (CAT) in the mucosal tissue of the cecum was measured using an ELISA kit (Angle Gene Biotechnology, Nanjing, China), whereas superoxide dismutase (SOD), glutathione S-transferase (GST), malonaldehyde (MDA), and lactate dehydrogenase (LDH) were measured using another ELISA kits (Nanjing Jiangchang Biotechnology, Nanjing, China) according to the manufacturer’s instructions.

The level of inflammatory markers in cecal mucosa

The LPS, secretory immunoglobulin A (SIgA), tumor necrosis factor alpha (TNF-α), interleukin (IL)-1α, IL-1β and IL-6, and Myeloperoxidase activity (MPO) in the mucosal tissue of the cecum were measured using an ELISA kit (Angle Gene Biotechnology, Nanjing, China).

Transcriptome analysis

To balance the cost and to select appropriate number of samples required for transcriptome analysis [14], total RNA from 9 cecal mucosa samples was extracted using the TRIZOL reagent (TaKaRa, Dalian, China). After extraction, the quality of the RNA was assessed by Nanodrop2000 (concentration and OD260/280 detection), agarose gel electrophoresis (RNA integrity detection), and Agilent2100 (RIN value measurement). Only high-quality total RNA samples (≥ 5 μg; ≥ 200 ng/μL; OD260/280 = 1.8 – 2.2) were used to construct sequencing libraries. The extracted samples were processed individually and later used to construct cDNA library following the protocol specified by the Truseq™ RNA sample prep kit (Illumina, San Diego, USA). Later, mRNA isolation and fragmentation (random fragments of 200 bp) were accomplished via Oligo (dT) and metal ions, respectively. mRNA fragments were transcribed into first-strand cDNA using reverse transcriptase and random primers, followed by second-strand cDNA synthesis. The double-stranded cDNA was further performed with end repaired, A base tailed, and indexed adapters were ligated. Then, the ligation product was amplified through polymerase chain reaction (PCR), and Illumina sequencing was performed by Biozerone, Shanghai, China. After an Illumina PE library was constructed, an Illumina platform was adopted to complete the 2 × 150 bp sequencing. Primary sequencing data (raw reads) were then subjected to quality control (QC). Tools such as SeqPrep ( and Sickle ( were also used. After obtaining clean reads, TopHat software (default parameters) was applied to align RNA-Seq reads to mammalian-sized genomes using the ultra-high-throughput short read aligner Bowtie, after which the mapping results were analyzed to identify splice junctions between exons.

A user-friendly software package, RSEM (, was used for quantifying gene and isoform abundances from single-end or paired-end RNA-Seq data. After sequencing, the expression level of each sequence library was standardized as FPKM, and the most differentially expressed genes (DEGs) were selected using EdgeR or deseq2 for further analysis.

Microbiome analysis

The number of animals used in our study for microbiota analysis were decided keeping in view the economic concerns and on the basis of number of animals selected in previous studies [15,16,17,18]. Total genomic DNA of cecal microbiota was extracted from all the 21 samples with 0.3 g of homogenized cecal digesta samples using a bead-beating and phenol-chloroform method [19]. After determining the purity and concentration of DNA, it was finally diluted (1.0 ng/mL). The diluted DNA was used to amplify the V3-V4 hypervariable regions in the 16S rRNA gene amplicons with barcoded primers (341F: 5′-CCTAYGGGRBGCASCAG-3′, 806R: 5′-GGACTACNNGGGTATCTAAT-3′), which is related to the lowest taxonomic assignment error rate for 16S rRNA MiSeq sequencing [20]. PCR was performed in triplicate. Amplicons were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel extraction kit (Axygen Biosciences, Union City, CA, USA) according to the manufacturer’s instructions, after which the amplicons were quantified using QuantiFluor-ST (Promega, Madison, WI, USA). The pooled DNA product was used to construct the Illumina Pair-End library following Illumina’s genomic DNA library preparation procedure. Purified amplicons were pooled in an equimolar and paired-end sequence (2 × 250) on an Illumina Miseq platform (Biozeron Biotechnology, Shanghai, China).

Raw fastq files were trimmed and chimeric sequences were identified and removed from all samples using UCHIME to reduce noise and to obtain Clean Reads. FLASH [21] was used to merge the paired-end reads generated from the DNA fragments. Non-assembled reads were discarded. The Clean Reads were then clustered as operational taxonomic units (OTUs) by scripts of the UPARSE software (version 7.1, http:// with a 97% similarity threshold. The phylogenetic affiliation of each 16S rRNA gene sequence was analyzed by a Ribosomal Database Project Classifier (RRID:SCR_006633) against the Silva (SSU115) 16S rRNA database (RRID:SCR_006423) using a confidence threshold of 70% [22].

Short chain fatty acid (SCFA) analysis

The cecal contents (stored at − 20 °C) from 21 samples were used for analysis of SCFAs using gas chromatography (GC2010 Plus, Shimadzu, Japan) as described previously [23]. Briefly, cecal digesta samples (50 mg) were suspended in 250 μL of ddH2O and were centrifuged at 12,000 × g for 5 min in a microcentrifuge (Microfuge 22R, Beckman Coulter, CA, USA). The SCFA analysis was carried out on the supernatants (200 μL), with crotonic acid as an internal standard. A flame ionization detector and a capillary column (Agilent Technologies, HP-INNOWax, 30 m × 0.25 mm × 0.25 μm, CA, USA) were used, with an injector/detector temperature of 180 °C/180 °C, a column temperature of 130 °C, and a gas flow rate of 30 mL/min.

Statistical analysis

To determine the effect of heat stress and MOLP supplementation, the data for biochemical traits such as anti-oxidative contents and inflammatory markers of cecal mucosa were analyzed by One-way ANOVA using SPSS 21.0 (SPSS Inc., Chicago, IL, USA) software with replicates as experiment units. The results are presented as mean ± standard error of mean (SEM). Significant differences were at P < 0.05. Significant differences due to treatments were separated by Student Newman Keuls multiple range test. Figures 1, 2 and 3 were constructed using the GraphPad Prism (version 8, San Diego, CA, USA).

Fig. 1

Effect of HS and Moringa oleifera leaf powder (MOLP) supplementation on organ index of cecum (n = 7 per group). MOLP: Moringa oleifera leaf powder; CON: control; HS: heat stress; HSM: heat stress with MOLP supplementation. All data is shown as mean values ± standard error of mean (SEM). a,b Means within groups (columns) with different superscript letters are significantly different (P < 0.05)

Fig. 2

Effect of HS and MOLP supplementation on biochemical traits of cecal mucosa in NZW rabbits (n = 7 per group). MOLP: Moringa oleifera leaf powder; CON: control; HS: heat stress; HSM: heat stress with MOLP supplementation; CAT: catalase; SOD: superoxide dismutase; GST: glutathione S-transferase; MDA: malonaldehyde; LDH: lactate dehydrogenase. All data is shown as mean values ± standard error of the mean (SEM). a,b Means within groups (columns) with different superscript letters are significantly different (P < 0.05)

Fig. 3

Effect of HS and MOLP supplementation on LPS, SIgA and markers of inflammation in cecal mucosa of NZW rabbits (n = 7 per group). MOLP: Moringa oleifera leaf powder; CON: control; HS: heat stress; HSM: heat stress with MOLP supplementation; LPS: lipopolysaccharide; sIgA: secretory immunoglobulin A; TNF-α: tumor necrosis factor alpha; IL-1α: interleukin-1α; IL-1β: interleukin-1β; IL-6: interleukin-6; MPO: myeloperoxidase. All data is shown as mean values ± standard error of the mean (SEM). a,b Means within groups (columns) with different superscript letters are significantly different (P < 0.05)

For transcriptomics data analysis, the significance levels of the transcripts were measured using the FDR (False Discovery Rate) control method to rationalize the P-value. The screening criterion to select the differential expressed genes between groups was used when threshold fold change |log2FC| ≥ = 1 & FDR < 0.05 [24, 25] Goatools ( was exploited to identify statistically significantly enriched GO term using Fisher’s exact test. KOBAS 2.0 ( was used to identify a statistically significant enriched pathway using Fisher’s exact test.

For the analysis of microbial data, the community richness estimator (Chao and ACE), diversity indices (Shannon and Simpson), and Good’s coverage as well as relative abundance of phylum and genera among groups were compared using the nonparametric Kruskul-Wallis H test in SPSS 21. Principal component analysis (PCA) was used to compare groups of samples based on unweighted UniFrac distance metrics [26]. Whereas, analysis of LEfSe was carried out using online tool (


MOLP supplementation improved growth parameters and organ index of cecum under HS

The data on growth performance has been given in Supplementary Table 2. Briefly, MOLP supplementation improved growth parameters such as ADG (25.8%, P < 0.05), ADFI (14.6%, P < 0.05) and FCR (12.4%, P < 0.05) compared to HS rabbits. The organ index of cecum to evaluate the effects of HS and Moringa supplementation at the end of trial were also recorded. Compared with the CON group, HS had profoundly reduced organ index of the cecum, whereas MOLP had improved organ index than that of HS group (Fig. 1, P < 0.05).

Effect of HS and MOLP supplementation on cecal oxidative stress

Further, the effect of HS and MOLP on cecal mucosa levels of oxidation and anti-oxidation status was explored. The results showed that MDA levels and SOD activity in cecal mucosa were notably higher in the HS group (Fig. 2, P < 0.05). Activity of Catalase and GST did not differ among the groups. HS significantly elevated LDH activity compared to the CON group, whereas MOLP reduced it (Fig. 2, P < 0.05).

Effect of HS and MOLP supplementation on markers of cecal mucosa inflammation

LPS content was significantly higher in the HS group compared with the CON group, whereas, LPS was restored to the normal level in the HSM group (Fig. 3, P < 0.05). SIgA was significantly lower in the HS group as compared to CON group and HSM groups (Fig. 3, P < 0.05). Inflammatory levels of TNF-α, IL-1α, and IL-1β were also significantly higher in the HS group compared with CON and HSM groups (P < 0.05), while they did not differ between the CON and HSM groups (P > 0.05). Myeloperoxidase (MPO) activity decreased significantly in the HSM group compared to HS group (P < 0.05). The contents of IL-6 did not differ significantly among the groups (P > 0.05).

Transcriptome analysis of cecal mucosa

To better understand changes in mucosa, transcriptome sequencing and differentially expressed gene identification were performed (Supplementary Table 3). The results showed that differentially expressed genes were rich on several KEGG pathways, especially the allograft rejection pathway and the fat digestion and absorption pathway (Fig. 4, P < 0.05). The main characteristic of the cecum exposed to HS is inflammatory immune response indicated by elevated gene expression such as IFNG, IGGC, and CD40LG and an element of the major histocompatibility complex (MHC) in the HS group compared to the CON group Table 1, P < 0.05). MOLP reduced IFNG, PRF1, granzyme B, IL-10, and CD86, and most elements of the MHC were also reduced, in the HSM group compared to the HS group (Table 1, P < 0.05).. Metabolic responses of cecum mucosa were also significantly changed in HS along with increased gene expression of APOA4, APOA1, FABP1, FABP2, MTTP, LOC100344166, LOC100343913, and PLA2G12B, and decreased LOC100339512 and PLA2G2F. MOLP increased FABP2, NPC1L1, ABCG8, and PLA2G1B, and decreased LOC100343913, PLA2G12B, and LOC100359181 (Table 2,  P < 0.05).

Fig. 4

KEGG route enrichment statistics. The X-axis shows the enrichment factor; the Y-axis corresponds to the KEGG pathway. The color of the dot represents the q value, and the dot size represents the number of differentially expressed genes (DEGs) mapped to the reference path. a, b, c represent statistics of the KEGG pathway enrichment of the DEGs of the HS vs. CON, HSM vs. CON, and HSM vs. HS groups, respectively

Table 1 Effect of HS and MOLP supplementation on the log fold change (logFC) of differentially expressed genes (DEGs) in pairwise comparison (HS vs. CON, HSM vs. CON and HSM vs. HS) that encode proteins in allograft rejection pathway
Table 2 The effect of HS and MOLP supplementation on the logfold change (logFC) of differentially expressed genes (DEGs) in pairwise comparison (HS vs. CON, HSM vs. CON and HSM vs. HS) that encode proteins in fat digestion and absorption pathway

Effect of MOLP on cecal microbiota and short chain fatty acids (SCFAs)

The statistical estimates of α-diversity from each sample at a genetic distance of 3%, presented in Supplementary Table 4, revealed no effects on any indices due to HS and dietary intervention, including the number of reads, OTUs, richness estimators (ACE and Chao1), and diversity indices (Shannon and Simpson) (Supplementary Table 4). PCA analysis on microbiota composition showed no distinction among groups (Supplementary Fig. 1). There was also no difference in the relative abundance of microbial genera among groups except for Papillibacter which was significantly higher (P= 0.018) in HSM group as compared to other groups (Supplementary Table 5).

Bacterial taxa (with a relative abundance of > 1% in any of the groups) were subjected to taxonomic composition analysis as previously described [27, 28]. At the phylum level, Firmicutes and Bacteroidetes were the most predominant phyla in the cecal digesta of rabbits, with a total relative abundance of around 90%. In descending order, the predominance of phyla was ranked as follows: Firmicutes > Bacteroidetes > Verrucomicrobia > Proteobacteria > Actinobacteria (Table 3). The relative abundance of Actinobacteria and Proteobacteria were significantly different between the groups (Table 3, P < 0.05), whereas no significant change in relative abundance was observed at the genus level except for Papillibacter (Supplementary Table 5). LDA effect size (LEfSe) analysis further revealed increased relative abundance of Actinobacteria in the HSM group and increased abundance of Proteobacteria and Bacteroidetes in HS group (Fig. 5P < 0.05). Total SCFA concentrations remained unchanged, whereas acetic acid and butyric acid were significantly higher in the HSM group as compared to the CON and HS groups (Table 4, P < 0.05).

Table 3 Effect of HS and MOLP supplementation on microbiota composition in cecal digesta of NZW rabbits at phylum level (relative abundance of top 5 phyla)
Fig. 5

Differential analysis of microbiota community within groups i.e. CON, HS and HSM (n = 7 per group). Linear discrimination analysis (LDA) coupled with effect size (LEfSe) was used to identify the most differentially abundant taxa among the three groups. Only the results meeting an LDA significant threshold of > 3 were shown

Table 4 Effect of HS and MOLP supplementation on the concentrations of short chain fatty acids (SCFAs) in cecal digesta in heat-stressed NZW rabbits (n = 7 per group)


HS is one of the largest impediments to efficient animal agriculture, resulting in hampered growth performance, altered carcass composition, and reduced fertility in rabbits [6, 29,30,31,32]. Traditionally, the negative impacts of HS on animal performance were attributed to associated decreases in feed intake, likely as a strategy to reduce metabolic heat production [7]. Similar results were obtained in our previous study demonstrating significantly reduced ADFI, ADG, and FCR in heat stressed rabbits [12]. Whereas Moringa supplementation significantly improved these traits. As an activated immune system (such that under HS) requires a substantial amount of energy, nutrients (specifically glucose and amino acids) are diverted from anabolic processes to support the immune response and thereby sustain growth and reproduction [33, 34]. In other words, supporting the activated immune system becomes a priority, and the hierarchy of nutrient partitioning is coordinated to ensure survival, whereas the promotion of economically important phenotypes is deemphasized.

Gastrointestinal tract is a major organ affected by heat stress [35]. In actuality, the many negative consequences of HS might be mediated by hyperthermia’s adverse effects on the intestinal epithelial barrier [8]. Liu et al. [36, 37] demonstrated that diurnal heat stress caused significant increase in permeability of small intestine as well as decrease in growth performance of heat stressed pigs compared with controlled ones. Similarly using a 3 day diurnal heat stress model on nursery-grower pig performance, intestinal integrity and endotoxemia; Gabler et al. [35] concluded that HS mediated reduction in ileum integrity and higher circulating LPS, altered blood metabolites and inflammation profiles might result in reduce performance of pigs over summer months. Moreover, dietary antioxidants (vitamin E & Se) have been described to mitigate the impacts of HS on intestinal barrier integrity, associated with a reduction in oxidative stress [37]. Keeping in view the results of current study, we can conclude that disruption of the intestinal epithelial barrier may be the indirect reason for the decreased growth performance and reduced organ index of cecum under HS.

HS causes a diversion of blood flow from the splanchnic bed to the periphery in a thermoregulatory effort to increase radiant heat loss. This altered cardiovascular event reduces blood flow and nutrient delivery to the gastrointestinal tract, resulting in local hypoxia, free radical production, and compromised intestinal architecture [38, 39]. The reactive oxygen species (ROS) readily attack and induce oxidative damage to various biomolecules, including proteins, lipids, lipoproteins, and DNA [40]. Lipid per oxidation (LPO) has been used as a reliable marker of oxidative stress, both in vitro and in vivo [41]. Antioxidant enzymes, such as SOD, catalase, and glutathione peroxidase (GPx), help to reduce oxidative stress in various parts of the cell [42]. These antioxidants may mediate their effect by directly reacting with ROS, quenching them and/or chelating the catalytic metal ions [43]. Herbs like Moringa oleifera contain free radical scavengers, such as polyphenols, flavonoids, and phenolic compounds [44,45,46,47]. Some polyphenol metabolites may have prebiotic, anti-inflammatory, anti-oxidative, anti-carcinogenic, and anti-microbial properties [44,45,46,47]. The present study demonstrated a significant reduction in cecal mucosal MDA content in rabbits fed a Moringa oleifera-supplemented diet. Whereas higher MDA levels and SOD activity in the HS group suggested lipid peroxidation and oxidative stress under HS environment respectively. This reduction could suggest a decrease in LPO and the potential of Moringa oleifera to prevent the formation of excessive free radicals by acting as chain-breaking peroxyl-radical scavengers and to protect low-density lipoproteins from oxidation [48]. Thus, the results of the current study highlight the antioxidant effects rendered by MOLP supplementation under oxidative stress conditions.

HS causes intestinal hyperpermeability—consequently, a dysfunctional intestinal barrier allows for the translocation of both dietary and microbial antigens (e.g., LPS from gram-negative bacteria) from the lumen through the basolateral membrane, eliciting local and systemic inflammatory responses [38, 39]. Accordingly, increased circulatory pro-inflammatory cytokines (e.g., TNFα, IL-1α, IL-1β, IL-6) have been reported in different HS models, presumably linked to increased circulating LPS [49,50,51,52,53]. In support of this, higher cecal mucosa levels of LPS, TNFα, IL-1α, IL-1β, and MPO due to HS, and lower levels due to MOLP, were observed. Further, based on the transcriptome analysis of cecal mucosa, the main characteristic of the cecum exposed to HS is inflammatory immune response indicated by elevated gene expression such as IFNG, IGGC, CD40LG, and an ingredient of the MHC compared to the CON group. MOLP down-regulated IFNG, PRF1, granzyme B, IL-10, and CD86, as well as most of the MHC, compared to the HS group. Up-regulation of IL-10 (secreted by regulatory T-lymphocytes, monocytes and macrophages) in cecal mucosa of HS group seems to inhibit the enhanced production of pro-inflammatory cytokines such as TNF-α and IFNG as reported previously [54]. Under inflammatory conditions, up-regulated IFNG (a potent pro-inflammatory cytokine) mediates both apoptotic and apoptosis-independent loss of intestinal epithelial integrity [55]. Downregulation of INFG in HSM group suggests that the Moringa may be exerting an anti-inflammatory immune response. Results are consistent with other studies indicating down-regulation of pro-inflammatory cytokines (e.g. INFG) with oral application of Lactobacillus-based products in mice [56, 57]. HS affected the functions of immunoglobulins, which play an important role in immune regulation and mucosal defense. IgA protects the mucosal surfaces of the intestine by preventing the entry, binding, and colonization of toxins and pathogens in both animal and human models [58]. The concentration of IgA was significantly lower in the HS group, whereas no difference was observed in the CON and HSM groups. This finding suggests that dietary Moringa oleifera was effective in eliciting a humoral immune response in rabbits under HS. Combined with the findings of our previous study, it can be suggested that HS causes an inflammatory response in both the jejunum [12] and cecum.

Despite marked reductions in feed intake and body weight loss, HS attenuates lipid mobilization in rodents [59] and ruminants [60,61,62,63]. Additionally, increased lipid carcass deposition occurs in heat-stressed growing animals [64], especially when compared to pair-fed thermoneutral controls [65, 66]. In our study, the results showed that both the HS and HSM groups had increased gene expression of fat digestion and absorption pathway, including APOA1, FABP1, FABP2, MTTP, and LOC100344166, compared to the CON group. Moringa has been described to possess anti-hyperlipidemia properties [67]. So, similar expression pattern of genes related to fat digestion and absorption pathway in HS (CON vs. HS) and HSM (CON vs. HSM) group still remains to be made clear for anti-hyperlipidemia effects of Moringa under HS.

HS and Moringa oleifera supplementation had a significant effect on the cecal microbial population at the phylum level assessed in this study. Firmicutes and Bacteroidetes were found to be the two most predominant phyla, followed by Verrucomicrobia, in accordance with previous studies [68], further followed by Proteobacteria and Actinobacteria. It has been reported that Firmicutes, Bacteroidetes, Proteobacteria, and Verrucomicrobia (90.5–97.7%) play an important role in the cecum, duodenum, and jejunum of growing rabbits [69]. In our study, dominant phyla altered by HS were Bacteroidetes (Fig. 5) and Proteobacteria (Table 3, Fig. 5), both being higher in abundance, whereas Moringa oleifera supplementation increased the abundance of Actinobacteria (Table 3). Whereas, at genus level, higher abundance of Papillibacter was found in HSM group (Supplementary Table 3). Shin et al. [70] proposed an increased prevalence of Proteobacteria as a potential diagnostic signature of dysbiosis and intestinal inflammation. Consistent findings have commonly supported an increased prevalence of Proteobacteria to the concept of dysbiosis during metabolic disorders [71, 72]. Likewise, a Chinese pilot study indicated dysbiosis of gut microbiota (significant enrichment of Proteobacteria and Bacilli) along with reduction in short chain fatty acids in patients with Encephalitis [73]. Moreover, gut dysbiosis with an excessive growth of Prevotellaceae and Bacteroidetes has been described in NLRP3 deficient mice under normal conditions [74].

Genus Papillibacter has been found naturally in Shea cake fed anaerobic digester, anaerobic sludge (from the pit of slaughter-houses) and environmental samples [75]. However, the metabolic and functional significance in the ruminal or cecal ecosystem is currently unknown. The only species of this genus ‘Papillibacter cinnaminivorans’ is phylogenetically related to the flavonoid ring-cleaving Clostridium orbiscindens anaerobe which is present in gastrointestinal tract of humans [75]. This may explain the reason why Polyphenol rich Moringa supplementation increased the abundance of Papillibacter in HSM group.

One of the cecum’s main functions is bacterial fermentation of indigestible polysaccharides to produce SCFAs, which can be utilized by the host’s epithelial cells. SCFAs can contribute 20–30% of the caloric requirements in rabbits [76]. Acetate, butyrate, and propionate are end products of the microbial fermentation of macronutrients. As epithelial cells utilize butyrate as their principal energy source, a ready source in the mucus is likely to support early epithelial maturation and growth [77]. SCFAs play an important role in gut physiology. Butyrate, in particular, is widely regarded as health-promoting. As stated, butyrate has been demonstrated as a mediator of anti-inflammatory responses, in the maintenance of the intestinal barrier integrity and in the protection against oxidative stress [78]. It can be stated that the higher butyrate levels indicate self-protection of the gastrointestinal tract under the Moringa oleifera supplementation regimen for HS. Although the results from current study warrant further investigation due to the limited number of rabbits used, yet the present study provided a logical basis for the use of Moringa oleifera in rabbit husbandry under HS conditions.


It can be summarized that HS results in change in redox status while Moringa supplementation reduced oxidative stress in mucosal tissue by reducing protein contents of MDA, SOD and LDH. Heat stress also resulted in dysbiosis of cecal microbiota and reduced integrity of cecal mucosa as evident by higher abundance of Proteobacteria and enhanced contents of LPS respectively, in HS group. Increased inflammatory immune response occurred in HS by elevated expression of pro-inflammatory genes such as IFNG and protein expression of inflammatory markers such as TNF-α, IL-1α, and IL-1β along with higher abundance of putatively pro-inflammatory gram negative Bacteroidetes. Whereas, Moringa supplementation reversed these conditions. Moringa supplementation also improved ADG and organ index of cecum. Moreover, higher abundance of genus Papillibacter and enhanced levels of microbial metabolite ‘Butyrate’ in HSM group demonstrate the ability of polyphenol rich Moringa to favorably module cecal microbial ecosystem. Concluding, Moringa supplementation modulated mucosal immune response and microbial ecosystem of cecum under HS by the mechanism of reduced oxidative stress and inflammatory markers along with reduction in the relative abundance of bacteria linked with dysbiosis.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.



Average daily gain


Average daily feed intake


Feed conversion ratio






Differentially expressed genes


False discovery rate


Heat stress


HS supplemented orally with MOLP


Glutathione S-transferase










Lactate dehydrogenase


LDA effect size






Moringa oleifera leaf powder




New Zealand White


Secretory immunoglobulin A


Short chain fatty acids


Operational taxonomic units


Superoxide dismutase


Tumor necrosis factor alpha


  1. 1.

    Rees Davies R, Rees Davies JAE. Rabbit gastrointestinal physiology. Vet Clin North Am - Exot Anim Pract. 2003;6(1):139–53.

    Article  Google Scholar 

  2. 2.

    Sohn J, Couto MA. Anatomy, Physiology, and Behavior. In The Laboratory Rabbit, Guinea Pig, Hamster, and Other Rodents. Elsevier Inc.; 2012. p. 195–215.

  3. 3.

    Sakaguchi E. Digestive strategies of small hindgut fermenters. Anim Sci J. 2003;74(5):327–37.

    Article  Google Scholar 

  4. 4.

    Belenguer Á, Balcells J, Guada JA, Decoux M, Milne E. Protein recycling in growing rabbits: contribution of microbial lysine to amino acid metabolism. Br J Nutr. 2005;94(5):763–70.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Li X, Min X, Xiao J, Kawasaki K, Ohta N, Sakaguchi E. Utilization of dietary urea nitrogen is stimulated by D-mannitol feeding in rabbits. Anim Sci J. 2012;83(8):605–9.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Marai IFM, Habeeb AAM, Gad AE. Rabbits’ productive, reproductive and physiological performance traits as affected by heat stress: a review. Livest Prod Sci. 2002;78(2):71–90.

    Article  Google Scholar 

  7. 7.

    Baumgard LH, Rhoads RP. Effects of heat stress on Postabsorptive metabolism and energetics. Annu Rev Anim Biosci. 2013;1(1):311–37.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Dokladny K, Zuhl MN, Moseley PL. Intestinal epithelial barrier function and tight junction proteins with heat and exercise. J Appl Physiol. 2016;120(6):692–701.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Aggarwal A, Upadhyay R. Heat stress and immune function. In: Heat stress and animal productivity. New Delhi: Springer India; 2013. p. 113–36.

  10. 10.

    Liu J, Xue C, Sun D, Zhu W, Mao S. Impact of high-grain diet feeding on mucosa-associated bacterial community and gene expression of tight junction proteins in the small intestine of goats. Microbiologyopen. 2019;8:1–12.

    Google Scholar 

  11. 11.

    Falowo AB, Mukumbo FE, Idamokoro EM, Lorenzo JM, Afolayan AJ, Muchenje V. Multi-functional application of Moringa oleifera Lam. in nutrition and animal food products: a review. Food Res Int. 2018;106:317–334., 2018. 

  12. 12.

    Khalid AR, Yasoob TB, Zhang Z, Yu D, Feng J, Zhu X, et al. Supplementation of Moringa oleifera leaf powder orally improved productive performance by enhancing the intestinal health in rabbits under chronic heat stress. J Therm Biol. 2020;93:102680.

  13. 13.

    Nakyinsige K, Sazili AQ, Zulkifli I, Goh YM, Abu Bakar F, Sabow AB. Influence of gas stunning and halal slaughter (no stunning) on rabbits welfare indicators and meat quality. Meat Sci. 2014;98:701–708.

  14. 14.

    Wu S, Liu Y, Duan Y, Wang F, Guo F, Yan F, et al. Intestinal toxicity of deoxynivalenol is limited by supplementation with Lactobacillus plantarum JM113 and consequentially altered gut microbiota in broiler chickens. J Anim Sci Biotechnol. 2018;9:1–13.

  15. 15.

    Alvarenga IC, Aldrich CG, Kohles M. The effect of feed form on diet digestibility and cecal parameters in rabbits. Animals. 2017;7:1–12.

    Article  Google Scholar 

  16. 16.

    Gao K, Pi Y, Mu CL, Peng Y, Huang Z, Zhu WY. Antibiotics-induced modulation of large intestinal microbiota altered aromatic amino acid profile and expression of neurotransmitters in the hypothalamus of piglets. J Neurochem. 2018;146(3):219–34.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Bo Trabi E, Seddik H Eldin, Xie F, Lin L, Mao S. comparison of the rumen bacterial community, rumen fermentation and growth performance of fattening lambs fed low-grain, pelleted or non-pelleted high grain total mixed ration. Anim Feed Sci Technol. 2019;253:1–12. 

  18. 18.

    Yu D, Zhu W, Hang S. Effects of long-term dietary protein restriction on intestinal morphology, digestive enzymes, gut hormones, and colonic microbiota in pigs. Animals. 2019;9:1–14.

    Google Scholar 

  19. 19.

    Dai ZL, Zhang J, Wu G, Zhu WY. Utilization of amino acids by bacteria from the pig small intestine. Amino Acids. 2010;39(5):1201–15.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Amato KR, Yeoman CJ, Kent A, Righini N, Carbonero F, Estrada A, Rex Gaskins H, Stumpf RM, Yildirim S, Torralba M, Gillis M, Wilson BA, Nelson KE, White BA, Leigh SR. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. ISME J. 2013;7(7):1344–53.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Zhao G, Nyman M, Jönsson JÅ. Rapid determination of short-chain fatty acids in colonic contents and faeces of humans and rats by acidified water-extraction and direct-injection gas chromatography. Biomed Chromatogr. 2006;20(8):674–82.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300 Wiley.

    Google Scholar 

  25. 25.

    Dalman MR, Deeter A, Nimishakavi G, Duan ZH. Fold change and p-value cutoffs significantly alter microarray interpretations. BMC Bioinformatics. 2012;13(Suppl 2):1–4.

    Article  Google Scholar 

  26. 26.

    Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Lin L, Xie F, Sun D, Liu J, Zhu W, Mao S. Ruminal microbiome-host crosstalk stimulates the development of the ruminal epithelium in a lamb model. Microbiome. 2019;7:1–16.

    Article  Google Scholar 

  28. 28.

    He J, Guo H, Zheng W, Xue Y, Zhao R, Yao W. Heat stress affects fecal microbial and metabolic alterations of primiparous sows during late gestation. J Anim Sci Biotechnol. 2019;10:1–12.

    Article  Google Scholar 

  29. 29.

    Chiericato GM, Boiti C, Canali C, Rizzi C, Ravarotto L. Effects of heat stress and age on growth performance and endocrine status of male rabbit. World Rabbit Sci. 1995;3:125–31.

    Google Scholar 

  30. 30.

    Dawod A. Impacts of ambient heat stress on growing rabbit performance and carcass traits. RRJoVST. 2015;4:7–13.

    Google Scholar 

  31. 31.

    Marai IFM, Ayyat MS, Abd El-Monem UM. Growth performance and reproductive traits at first parity of New Zealand white female rabbits as affected by heat stress and its alleviation under Egyptian conditions. Trop Anim Health Prod. 2001;33(6):451–62.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    El-Desoky NI, Hashem NM, Elkomy A, Abo-Elezz ZR. Physiological response and semen quality of rabbit bucks supplemented with Moringa leaves ethanolic extract during summer season. Animal. 2017;11(9):1549–57.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Waldron MR, Kulick AE, Bell AW, Overton TR. Acute experimental mastitis is not causal toward the development of energy-related metabolic disorders in early postpartum dairy cows. J dairy Sci. Am Dairy Sci Assoc. 2006;89(2):596–610.

    CAS  Article  Google Scholar 

  34. 34.

    Kvidera SK, Horst EA, Abuajamieh M, Mayorga EJ, Fernandez MVS, Baumgard LH. Glucose requirements of an activated immune system in lactating Holstein cows. J Dairy Sci. 2017;100:2360–74.

  35. 35.

    Gabler NK, Koltes D, Schaumberger S, Murugesan GR, Reisinger N. Diurnal heat stress reduces pig intestinal integrity and increases endotoxin translocation. Transl Anim Sci. 2018;2(1):1–10.

  36. 36.

    Liu F, Celi P, Cottrell JJ, Chauhan SS, Leury BJ, Dunshea FR. Effects of a short-term supranutritional selenium supplementation on redox balance, physiology and insulin-related metabolism in heat-stressed pigs. J Anim Physiol Anim Nutr (Berl). 2018;102:276–85.

  37. 37.

    Liu F, Cottrell JJ, Furness JB, Rivera LR, Kelly FW, Wijesiriwardana U, et al. Selenium and vitamin E together improve intestinal epithelial barrier function and alleviate oxidative stress in heat-stressed pigs. Exp Physiol. 2016;101:801–10.

  38. 38.

    Lambert GP. Stress-induced gastrointestinal barrier dysfunction and its inflammatory effects. J Anim Sci. 2009;87:101–8.

    Article  Google Scholar 

  39. 39.

    Ogden HB, Child RB, Fallowfield JL, Delves SK, Westwood CS, Layden JD. The gastrointestinal exertional heat stroke paradigm: pathophysiology, assessment, severity, aetiology and nutritional countermeasures. Nutrients. 2020;12:1–42.

    Google Scholar 

  40. 40.

    Farber JL. Mechanisms of cell injury by activated oxygen species. JSTOR. 1994;102:17.

  41. 41.

    Sreelatha S, Padma PR. Antioxidant activity and total phenolic content of Moringa oleifera leaves in two stages of maturity. Plant Foods Hum Nutr. 2009;64(4):303–11.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Weydert CJ, Cullen JJ. Measurement of superoxide dismutase, catalase and glutathione peroxidase in cultured cells and tissue. Nat Protoc. 2010;5(1):51–66.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Robak J, Marcinkiewicz E. Scavenging of reactive oxygen species as the mechanism of drug action. Pol J Pharmacol. 1995;47:89–98. Available from:

  44. 44.

    Ma ZF, Ahmad J, Zhang H, Khan I, Muhammad S. Evaluation of phytochemical and medicinal properties of Moringa (Moringa oleifera) as a potential functional food. S Afri J Bot. 2020;129:40–6.

    CAS  Article  Google Scholar 

  45. 45.

    Moyo B, Oyedemi S, Masika PJ, Muchenje V. Polyphenolic content and antioxidant properties of Moringa oleifera leaf extracts and enzymatic activity of liver from goats supplemented with Moringa oleifera leaves/sunflower seed cake. Meat Sci. 2012;91:441–7. 

  46. 46.

    Bholah K, Ramful-Baboolall D, Neergheen-Bhujun VS. Antioxidant activity of Polyphenolic rich Moringa oleifera Lam. extracts in food systems. J Food Biochem. 2015;39(6):733–41.

  47. 47.

    Saini RK, Sivanesan I, Keum YS. Phytochemicals of Moringa oleifera: a review of their nutritional, therapeutic and industrial significance. 3 biotech, vol. 6. Berlin Heidelberg: Springer; 2016. p. 1–14.

    Google Scholar 

  48. 48.

    O’Byrne DJ, Devaraj S, Grundy SM, Jialal I. Comparison of the antioxidant effects of Concord grape juice flavonoids and α-tocopherol on markers of oxidative stress in healthy adults. Am J Clin Nutr. 2002;76(6):1367–74.

    Article  PubMed  Google Scholar 

  49. 49.

    Lambert GP, Gisolfi CV, Berg DJ, Moseley PL, Oberley LW, Kregel KC. Selected contribution: hyperthermia-induced intestinal permeability and the role of oxidative and nitrosative stress. J Appl Physiol. 2002;92:1750–61.

    CAS  Article  Google Scholar 

  50. 50.

    Alhenaky A, Abdelqader A, Abuajamieh M, Al-Fataftah AR. The effect of heat stress on intestinal integrity and Salmonella invasion in broiler birds. J Therm Biol. 2017;70:9–14.

  51. 51.

    Bouchama A, Parhar RS, El-Yazigi A, Sheth K, Al-Sedairy S. Endotoxemia and release of tumor necrosis factor and interleukin 1α in acute heatstroke. J Appl Physiol. 1991;70(6):2640–4.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Bouchama A, Al-Sedairy S, Siddiqui S, Shail E, Rezeig M. Elevated pyrogenic cytokines in heatstroke. Chest. 1993;104(5):1498–502.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Johnson JS, Sapkota A, Lay DC. Rapid cooling after acute hyperthermia alters intestinal morphology and increases the systemic inflammatory response in pigs. J Appl Physiol. 2016;120(10):1249–59.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Gracie JA, Forsey RJ, Chan WL, Gilmour A, Leung BP, Greer MR, et al. A proinflammatory role for IL-18 in rheumatoid arthritis. J Clin Invest. 1999;104(10):1393–401.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Mashukova A, Wald FA, Salas PJ. Tumor necrosis factor alpha and inflammation disrupt the polarity complex in intestinal epithelial cells by a posttranslational mechanism. Mol Cell Biol. 2011;31(4):756–65.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Matsumoto S, Hara T, Hori T, Mitsuyama K, Nagaoka M, Tomiyasu N, et al. Probiotic Lactobacillus-induced improvement in murine chronic inflammatory bowel disease is associated with the down-regulation of pro-inflammatory cytokines in lamina propria mononuclear cells. Clin Exp Immunol. 2005;140(3):417–26.

  57. 57.

    Maassen CBM, Van Holten-Neelen C, Balk F, Heijne Den Bak-Glashouwer MJ, Leer RJ, Laman JD, et al. Strain-dependent induction of cytokine profiles in the gut by orally administered Lactobacillus strains. Vaccine. 2000;18(23):2613–23.

  58. 58.

    Phalipon A, Cardona A, Kraehenbuhl JP, Edelman L, Sansonetti PJ, Corthésy B. Secretory component: a new role in secretory IgA-mediated immune exclusion in vivo. Immunity. 2002;17:107–15.

  59. 59.

    Sanders SR, Cole LC, Flann KL, Baumgard LH, Rhoads RP. Effects of acute heat stress on skeletal muscle gene expression associated with energy metabolism in rats. FASEB J. 2009;23:598–7.

  60. 60.

    Rhoads RP, Sanders SR, Cole L, Skrzypek MVET. Effects of heat stress on glucose homeostasis and metabolic response to an endotoxin challenge in Holstein steers. J Anim Sci. 2009;87:78.

    Google Scholar 

  61. 61.

    Baumgard LH, Wheelock JB, Sanders SR, Moore CE, Green HB, Waldron MR, et al. Postabsorptive carbohydrate adaptations to heat stress and monensin supplementation in lactating Holstein cows. J Dairy Sci; 2011;94:5620–5633. 

  62. 62.

    Wheelock JB, Rhoads RP, VanBaale MJ, Sanders SR, Baumgard LH. Effects of heat stress on energetic metabolism in lactating Holstein cows. J Dairy Sci; 2010;93:644–655. 

  63. 63.

    O’Brien MD, Rhoads RP, Sanders SR, Duff GC, Baumgard LH. Metabolic adaptations to heat stress in growing cattle. Domest Anim Endocrinol. 2010;38(2):86–94.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Geraert PA, Padilha JCF, Guillaumin S. Metabolic and endocrine changes induced by chronic heatexposure in broiler chickens: growth performance, body composition and energy retention. Br J Nutr. 1996;75:195–204.

  65. 65.

    Pearce SC, Gabler NK, Ross JW, Escobar J, Patience JF, Rhoads RP, Baumgard LH. The effects of heat stress and plane of nutrition on metabolism in growing pigs. J Anim Sci. 2013;91(5):2108–18.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Sanz Fernandez MV, Johnson JS, Abuajamieh M, Stoakes SK, Seibert JT, Cox L, et al. Effects of heat stress on carbohydrate and lipid metabolism in growing pigs. Physiol Rep. 2015;3:1–17.

    Article  Google Scholar 

  67. 67.

    Omodanisi EI, Aboua YG, Chegou NN, Oguntibeju OO. Hepatoprotective, antihyperlipidemic, and anti-inflammatory activity of Moringa oleifera in diabetic-induced damage in male wistar rats. Pharm Res. 2017;9:182–7.

    CAS  Google Scholar 

  68. 68.

    Bäuerl C, Collado MC, Zúñiga M, Blas E, Martínez GP. Changes in cecal microbiota and mucosal gene expression revealed new aspects of epizootic rabbit enteropathy. PLoS One. 2014;9(8):e105707.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Wang C, Zhu Y, Li F, Huang L. The effect of Lactobacillus isolates on growth performance, immune response, intestinal bacterial community composition of growing Rex rabbits. J Anim Physiol Anim Nutr (Berl). 2017;101:e1–13.

  70. 70.

    Shin NR, Whon TW, Bae JW. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 2015;33:496–503.

  71. 71.

    Zhang X, Shen D, Fang Z, Jie Z, Qiu X, Zhang C, et al. Human gut microbiota changes reveal the progression of glucose intolerance. PLoS One. 2013;8(8):e71108.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Larsen N, Vogensen FK, Van Den Berg FWJ, Nielsen DS, Andreasen AS, Pedersen BK, et al. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS One. 2010;5(2):e9085.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  73. 73.

    Xu R, Tan C, He Y, Wu Q, Wang H, Yin J. Dysbiosis of gut microbiota and short-chain fatty acids in encephalitis: a Chinese pilot study. Front Immunol. 2020;11:1–13.

    Article  Google Scholar 

  74. 74.

    Henao-Mejia J, Elinav E, Jin C, Hao L, Mehal WZ, Strowig T, et al. Inflammasome-mediated dysbiosis regulates progression of NAFLD and obesity. Nature. 2012;482(7384):179–85.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Patel BKC. Papillibacter. Bergey’s Man Syst Archaea Bact. 2015:1–4.

  76. 76.

    Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70(2):567–90.

    CAS  Article  PubMed  Google Scholar 

  77. 77.

    Polansky O, Sekelova Z, Faldynova M, Sebkova A, Sisak F, Rychlik I. Important metabolic pathways and biological processes expressed by chicken cecal microbiota. Appl Environ Microbiol. 2016;82(5):1569–76.

    CAS  Article  PubMed Central  Google Scholar 

  78. 78.

    SJD O’K. Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol. 2016;13:691–706.

Download references


Not applicable.


This work was supported by the National Key Basic Research Program of China (973 Program) under grant 2013CB127301.

Author information




Talat Bilal Yasoob: Formal analysis, Writing-original draft. Defu Yu: Formal analysis, Writing-original draft. Zhen Zhang: Formal analysis. Abdur Rauf Khalid: Formal analysis. Heba M. Saad: Formal analysis. Xiaofeng Zhu: Formal analysis. Suqin Hang: Formal analysis, Writing-original draft. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Suqin Hang.

Ethics declarations

Ethics approval and consent to participate

The study procedures were adopted after approval by Nanjing Agricultural University (LicenseNo.SYXK-2017-0027).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1: Table S1.

Nutrient composition per 100 g of Moringa oleifera leaf powder (MOLP).

Additional file 2: Table S2.

Effect of heat stress and MOLP supplementation on productive performance of New Zealand White rabbits.

Additional file 3: Table S3.

Summary of the sample comparison results of differentially expressed genes (DEGs).

Additional file 4: Table S4.

Summary of Illumina Miseq sequence data and statistical analysis of the bacterial diversity in the cecum of CON, HS, and HSM groups (n = 7 per group).

Additional file 5: Table S5.

Effect of HS and MOLP supplementation on the relative abundance of microbial genera (with percentages greater than 0.5% in any of groups) in cecum of NZW rabbits (n = 7 per group).

Additional file 6: Figure S1.

Unweighted UniFrac Principal component analysis (PCA) of microbiota composition in cecum based on operational taxonomic unit (OUT) data (n = 7 per gorup). CON: control group; HS: heat stress group; HSM: heat stress with MOLP supplementation group; MOLP: Moringa oleifera leaf powder.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yasoob, T.B., Yu, D., Khalid, A.R. et al. Oral administration of Moringa oleifera leaf powder relieves oxidative stress, modulates mucosal immune response and cecal microbiota after exposure to heat stress in New Zealand White rabbits. J Animal Sci Biotechnol 12, 66 (2021).

Download citation


  • Cecal microbiota
  • Heat stress
  • Moringa oleifera leaf powder
  • Rabbits
  • Short chain fatty acids