Skip to main content

Exploring the predictive power of jejunal microbiome composition in clinical and subclinical necrotic enteritis caused by Clostridium perfringens: insights from a broiler chicken model

Abstract

Background

Necrotic enteritis (NE) is a severe intestinal infection that affects both humans and poultry. It is caused by the bacterium Clostridium perfringens (CP), but the precise mechanisms underlying the disease pathogenesis remain elusive. This study aims to develop an NE broiler chicken model, explore the impact of the microbiome on NE pathogenesis, and study the virulence of CP isolates with different toxin gene combinations.

Methods

This study established an animal disease model for NE in broiler chickens. The methodology encompassed inducing abrupt protein changes and immunosuppression in the first experiment, and in the second, challenging chickens with CP isolates containing various toxin genes. NE was evaluated through gross and histopathological scoring of the jejunum. Subsequently, jejunal contents were collected from these birds for microbiome analysis via 16S rRNA amplicon sequencing, followed by sequence analysis to investigate microbial diversity and abundance, employing different bioinformatic approaches.

Results

Our findings reveal that CP infection, combined with an abrupt increase in dietary protein concentration and/or infection with the immunosuppressive variant infectious bursal disease virus (vIBDV), predisposed birds to NE development. We observed a significant decrease (p < 0.0001) in the abundance of Lactobacillus and Romboutsia genera in the jejunum, accompanied by a notable increase (p < 0.0001) in Clostridium and Escherichia. Jejunal microbial dysbiosis and severe NE lesions were particularly evident in birds infected with CP isolates containing cpa, netB, tpeL, and cpb2 toxin genes, compared to CP isolates with other toxin gene combinations. Notably, birds that did not develop clinical or subclinical NE following CP infection exhibited a significantly higher (p < 0.0001) level of Romboutsia. These findings shed light on the complex interplay between CP infection, the gut microbiome, and NE pathogenesis in broiler chickens.

Conclusion

Our study establishes that dysbiosis within the jejunal microbiome serves as a reliable biomarker for detecting subclinical and clinical NE in broiler chicken models. Additionally, we identify the potential of the genera Romboutsia and Lactobacillus as promising candidates for probiotic development, offering effective alternatives to antibiotics in NE prevention and control.

Introduction

Necrotic enteritis (NE) is a severe inflammatory infection that affects the small intestine, specifically the jejunum and ileum. This condition can occur in both humans and poultry. The causative agent of NE is Clostridium perfringens (CP), a gram-positive, rod-shaped, spore-forming bacterium that thrives in various environments like soil, litter, and dust. In chickens, CP can act as an opportunistic pathogen. However, in humans, the primary mode of transmission is through the consumption of contaminated food, including undercooked meat or poultry [1]. CP is divided into seven groups (A to G) based on major toxin producing genes including alpha (cpa), beta (cpb), epsilon (etx), iota (itx), enterotoxin (cpe) and necrotic enteritis B like toxin (netB) [2]. Affected broiler chickens develop NE between 2 and 6 weeks of age [3]. The disease has clinical and subclinical forms; the clinical form includes depression, diarrhea, dehydration, ruffled feathers, low body weight and increased feed conversion ratio (FCR) [4]. Mortality can reach up to 1% per day with a total mortality of 10–40% [5]. The subclinical form of NE is the most common in the field and results in an increased feed conversion ratio (FCR) and reduction in weight gain which leads to overall poor performance and economic losses [6, 7]. Economic losses are expected to rise as many countries as possible [8, 9], including Canada are voluntarily reducing and withdrawing the prophylactic use of antibiotics. Since there is no effective control available against NE, development of alternative prophylactics and therapeutic agents is of utmost importance. In this context, chickens serve as valuable animal models for studying NE as they share some similarities with humans in terms of disease manifestation [10,11,12], host-pathogen interaction [13], genetics [14], intestinal physiology [15], and experimental feasibility[16], thereby facilitating the investigation of potential therapeutic interventions.

The pathogenesis of CP is poorly understood. However, several predisposing factors such as infections with coccidia, immunosuppressive viruses, and abrupt change in protein contents of the diet have been identified as contributing factors for the adhesion and colonization of the bacteria in the intestinal mucosa. In addition, cpa, [17], netB [18] and tpeL [19, 20] toxins play important roles in the pathogenesis of NE. Of these, cpa mainly disrupts the membrane of enterocytes by hydrolysing phosphatidylcholine and sphingomyelin, producing discylglycerol and ceramide respectively [21, 22]. Furthermore, cpa activates the arachidonic acid pathway inducing the production of prostaglandins, tromboxanes, leukotriens and lactate dehydrogenase and result in necrosis of the intestine [23]. NetB is a pore forming toxin [24] and permits an influx of Na+, Cl- and Ca+2 into the cytoplasm and causes osmotic lysis of the cell [25, 26]. Most clinical isolates of CP produce netB [27]. NetB mutants were unable to produce NE in broiler chickens [18]. The co-expression of tpeL and netB by CP can lead to a rapid onset of NE with higher mortality [19, 20]. TpeL induces host cell apoptosis by inhibiting the Ras signalling pathway [28, 29].

The chicken intestinal microbiota plays a crucial role in the growth and health of chickens, contributing towards improved nutrient absorption, improved immune system and inhibition of pathogen colonization by competitive exclusion [30]. Multiple studies demonstrated that Eimeria sp. induced alteration of the cecal microbiome can predispose chickens to the development of NE [31, 32]. However, the impact of the change in the microbial community in the pathogenesis of NE is not clearly understood [32, 33]. The objectives of this study were to; (a) develop subclinical and/or clinical NE broiler chicken animal model by inducing immunosuppression using infectious bursal disease virus (IBDV) and/or by abrupt increase of the protein content in the chicken feed, (b) investigate the impact of the intestinal microbiome in the pathogenesis of NE, and (c) study the virulence of CP isolates containing different combinations (cpa, netB, tpeL and/or cpb2) of toxin genes including their impact on the jejunal microbiome.

Results

The experimental procedures undertaken in this study are outlined in a schematic representation featured in Fig. 1.

Fig. 1
figure 1

Experimental design and sample processing for Clostridium perfringens (CP) animal challenge. This figure provides a diagrammatic representation of the CP animal challenge conducted during both experiments 1 and 2. It delineates the timeline for sample collection and the subsequent processing steps essential for Illumina MiSeq-based next-generation sequencing, targeting the 16S rRNA gene’s V3–V4 region. Image “created with www.BioRender.com” accessed on October 22, 2023

Toxin gene makeup of CP isolates

Among the 57 CP isolates, 44 and 12 isolates were classified as toxinotype G and A, respectively. However, one isolate contained cpa+netB+cpiap+cpb2 and did not fit into the current CP toxinotype classification. Overall, 7 toxin gene combinations were found among the analyzed CP clinical isolates including cpa+netB+cpb2(54.4%), cpa+netB (12.28%), cpa+netB+cpb2+tpeL (10.5%), cpa+cpb2 (8.77%), cpa alone (10.5%), cpa+netB+cpiap+cpb2 (1.7%), and cpa+cpb2+tpeL (1.7%). None of these isolates harbored cpe, cpb or etx toxin genes.

Development of a NE animal model

Clinical signs and mortality

In the developed animal infection model, a total mortality of 10% and 15% (p=0.0173) was observed in the group with vIBDV SK09 + CP by oral gavage (28% protein) and vIBDV SK09 + CP in feed (20% protein) group, respectively (Fig 2A). Death was peracute after the development of clinical signs due to toxemia. The rest of the infected birds were apparently normal until the end of the clinical trial and the only clinical sign observed was diarrhea. At 2 days post-infection (DPI), one bird was found dead in both the vIBDV SK09 + CP by oral gavage (28% protein) group and vIBDV SK09 + CP in feed (20% protein) group. At 3 DPI, one and two birds were found dead in the vIBDV SK09 + CP by oral gavage (28% protein) and vIBDV SK09 + CP in feed (20% protein) groups, respectively (Fig 2A). No mortality was observed in the other experimental groups throughout the study period.

Fig. 2
figure 2

Survival and macroscopic NE lesions of groups of broiler chickens (n = 20/group) post-CP challenge. A Highlights a mortality rate ranging between 10 to 15% in groups of birds pre-exposed to vIBDV SK09 before the CP challenge (p = 0.0173). Key terms: CP: Clostridium perfringens; P: protein; vIBDV: variant infectious bursal disease virus. B Displays a dilated jejunum with patchy, thin, necrotizing areas (arrow) of the intestinal mucosa through the serosal surface of the intestine. C Depicts severe, diffuse, acute necrosis of the intestinal mucosa in the mucosal surface of the intestine

Gross and microscopic lesions

On gross examination, intestines were dilated with gas and patchy necrotic areas were visible through the serosal surface (Fig 2B). Upon opening the intestines, severe, diffuse necrosis was evident in the duodenum, jejunum, and the ileum in the dead or euthanized birds (Fig 2C). Gross NE lesions were observed at a frequency of 5%, 5%, 15%, 10% and 20% in Groups 2 (CP feed, 20% protein), 3 (CP feed, 28% protein), 4 (vIBDV + CP feed [20% protein]), 5 (vIBDV + CP feed [28% protein]) and 6 (vIBDV + CP oral gavage [28% protein]), respectively. No gross lesions were observed in the unchallenged control group. In this NE model, intestinal necrosis was pronounced and prevalent in the jejunum compared to the duodenum or ileum. The histopathological lesions were scored from 0 to 3 (Fig 3A–D). Some of the CP exposed birds had extensive, severe, diffuse necrosis of intestinal villi with the presence of occasional bacterial colonies (Fig 3D) with few birds harboring fibrinonecrotic debris in the lumen. Groups 2, 3, 4, 5 and 6 had histopathological lesions of NE at a frequency of 5%, 55%, 70%, 70% and 55%, respectively (Fig 3E, p < 0.0001). No histologic lesions were observed in the unchallenged negative control group.

Fig. 3
figure 3

Histopathological analysis and scoring of jejunal NE lesions in broiler chickens (n = 20/group). A Represents the normal intestinal mucosa (0 score), B shows mild, focal, acute necrosis of intestinal villi (1 score), C shows moderate, multifocal to coalescing, acute necrosis of intestinal villi (2 score) and D represents severe, diffuse, acute necrosis of intestinal villi (3 score). E percentage microscopic NE lesions in birds (n = 20/group) which were predisposed to abrupt change in protein content in the diet and/or infection with vIBDV followed by exposure to CP (p < 0.0001). (CP: Clostridium. perfringens; P20: 20% protein in the feed; P28: 28% protein in the feed; vIBDV: variant infectious bursal disease virus)

16S rRNA microbiome sequence summary

Next generation 16S rRNA gene sequencing-based microbiome analysis was performed to examine the pattern of dysbiosis (i.e., a shift in the composition and relative abundance of the jejunal microbiome of the intestinal microbial community in CP infected birds as compared to the non-infected control group with healthy intestine). The average number of raw sequences read obtained per sample was 914,290 and it was consistent within and between the groups. After trimming/filtering and removing reads with quality scores of less than 30, an average of 679,150 reads per sample was obtained. The average length of the assembled reads used for taxonomic classification was 465 base pairs (bp). All the FastQ files are deposited in the GenBank with Bio project numbers PRJNA1023687 and PRJNA1024573.

Jejunal relative microbial abundance at the phylum level

The relative abundance of the jejunal microbial composition was presented at the level of phyla, family, and genus (Fig 4A–C). A total of 16 phyla were identified in jejunum with the phylum Firmicutes the most abundant in birds that developed clinical or subclinical NE, and in CP unchallenged negative control groups. However, a significant increase in the level of Firmicutes was noted in birds that developed clinical (89.1%) and subclinical NE as compared to the negative control group (76%, p < 0.0180) (Fig 4A). The second most abundant phylum was Cyanobacteria which constituted 21.8% in the negative control groups. However, the level of Cyanobacteria was significantly reduced in CP challenged birds that developed clinical (3.21%), or subclinical (11.49%) NE as compared to the unchallenged control group (p < 0.0001). The third most abundant phylum was Proteobacteria which significantly increased in birds with clinical NE (7.54%) as compared to the negative control group (0.44%, p < 0.049). Although the relative abundance of Proteobacteria showed a dramatic increase (0.95% to 4.34%) in birds that developed subclinical NE as compared to the negative control group, it was not statistically significant (p < 0.440) (Fig 4A).

Fig. 4
figure 4

Analysis of jejunal microbiome in broiler chickens with subclinical or clinical NE based on 16S amplicon sequencing. A Depicts the relative abundance distribution at the phylum level. B Illustrates the relative abundance distribution at the family level. C Demonstrates the relative abundance distribution at the genus level. D Presents a principal component analysis biplot (n = 8/group) that visually captures the distances and relationships among three distinct groups: C1–C8 (healthy birds), SB1–SB8 (subclinical NE), and CL1–CL8 (clinical NE). Within this biplot, "Loadings" signify the weights of each original variable, while “PC scores” represent the linear combinations of the variables in the dataset

Jejunal relative microbial abundance at the family level

At the family level, the negative control group had predominantly Lactobacillaceae family (64.3%) followed by Peptostreptococcaceae (14.8%), Oscillatoriophycideae (13.4%), and other families (7%) (i.e., others represent families with relative abundance of less than 2%) including Clostridiaceae (1.46%) and Enterobacteriaceae (0.35%) (Fig 4B). Birds with subclinical and clinical NE had a significant decrease in Lactobacillaceae and Oscillatoriophycideae (p<0.0001) as compared to the CP unchallenged control group. There was significant increase in Clostridiaceae in birds with clinical NE (70.4%) followed by subclinical NE (35.3%) as compared to birds not exposed to CP (p<0.0001) (Fig 4B). An increased trend of Enterobacteriaceae was observed as CP infected birds develop subclinical NE (3.40%, p<0.048) and transition to clinical NE (6.42%, p<0.046) (Fig 4B).

Jejunal relative microbial abundance at the genus level

The genus Lactobacillus (54.9%) was the dominant genus in the family Lactobacillaceae (Phylum Firmicutes) in the CP unchallenged control group. In contrast, the genus Lactobacillus was significantly lower in birds with subclinical or clinical NE (p < 0.0001) as compared to the negative control group (Fig. 4C). As expected, the relative abundance of the genus Clostridium (Phylum Firmicutes) was significantly higher in birds that developed subclinical (32.74%) or clinical NE (54.89%) as compared to birds not exposed to CP (1.79%, p < 0.0001) (Fig. 4C). Interestingly, the relative abundance of the genus Romboutsia (Family Peptostreptococcaceae; Phylum Firmicutes) increased in birds that developed subclinical NE (25.57%) as compared to the negative control group (17.73%, p < 0.0323), and significantly reduced in birds that developed clinical NE (16.57%) as compared to birds that developed subclinical NE (p < 0.014). In contrast, there was an increasing trend of the relative abundance of the genus Escherichia (family Enterobacteriaceae, phylum Proteobacteria) as birds develop clinical NE. The proportion of genus Escherichia rose to 4.14% in birds that developed subclinical NE as compared to the CP unchallenged control group (0.31%, p < 0.451). The number of birds that developed clinical NE (8.82%) was statistically significant (p < 0.0216) (Fig. 4C). Based on principal component analysis, the Lactobacillaceae, Enterobacteriaceae and Clostridiaceae families had distinct clustering between the CP unchallenged control group and birds that developed subclinical or clinical NE (Fig. 4D).

Alpha and beta-jejunal microbial diversity of sub-clinical and clinical NE

The microbial diversity within each sample (alpha diversity) of the control birds, birds with subclinical or clinical NE post CP infection was examined by computing richness (observed operational taxonomic units, OTUs), Shannon index, Simpson index and Pielou’s evenness. The box plots generated to visualize the distribution of each metric across the groups displayed a clear distinction in the microbial diversity within the groups (Fig. 5A). The differences in microbial composition between the control, subclinical NE and clinical NE group was assessed by beta diversity analysis by measuring the weighted and unweighted UniFrac distances which were visualized by Principal Co-ordinate Analysis (PCoA) (Fig. 5B). Unweighted UniFrac distances demonstrated statistically significant differences in microbial community structures between the control, subclinical NE and clinical NE groups (p < 0.007; R2 = 0.1226; F-value = 3.4001). However, when considering taxa abundance (weighted UniFrac distances), the differences between the groups were not statistically significant (p < 0.229). The variations in the microbial communities are somewhat influenced by the severity of NE, more notably in the unweighted analysis.

Fig. 5
figure 5

Jejunal microbiome analysis in control and NE-affected birds post-CP infection. A Alpha diversity indices of the jejunal microbiome with box plots indicating distribution across control and NE-affected groups (n = 8/group). B PCoA plots, based on unweighted and weighted UniFrac distances, illustrating microbial community compositions in the same groups

Severity of NE associated with CP isolates containing different combinations of toxin genes

Gross and microscopic lesions

The association of the presence of different combinations of toxin genes in CP isolates and the severity of NE in infected birds was evaluated in the CP broiler chicken infection model. Group 2 (infected with CP containing cpa, netB, cpb2 and tpeL genes), Group 3 (infected with CP containing cpa, cpb2 and netB genes); and Group 4 (infected with CP containing cpa and cpb2 genes) had gross NE lesions at the rate of 15%, 5% and 10%, respectively. No gross lesions were noted in group 1 (unchallenged control group). Groups 2, 3, and 4 had histopathological lesions of NE at the rate of 90%, 20% and 10%, respectively. No histopathologic lesions were observed in the unchallenged negative control group (Fig. 6A).

Fig. 6
figure 6

Impact of different CP Toxin gene combinations on NE histologic lesions and jejunal microbiome composition in broiler chickens. A Percentage of broiler chickens (n = 20/group) manifesting NE histopathological lesions following exposure to different CP toxin gene combinations (p = 0.0001). B Relative abundance (n = 20/group, n = 10 for control) of the jejunal microbiome at the phylum level post-CP exposure. C Composition of the jejunal microbiome at the family level. D Detailed breakdown at the genus level of the jejunal microbiome. CP: Clostridium perfringens

Jejunal relative microbial abundance at the phylum level

The jejunal microbial diversity was presented as relative abundance at the phyla, family, and genus levels (Fig 6B–D). A total of 16 phyla were identified. Like the first experiment, the phylum Firmicutes was the most abundant in all the experimental groups. However, the proportion of Firmicutes increased in birds exposed to cpa+netB+tpeL+cpb2 (86.45%, p < 0.0017), cpa+netB +cpb2 (85.55%, p = 0.0053), or cpa+cpb2 (94.34%, p < 0.0001) as compared to unchallenged control birds (77.52%) (Fig 6B). Cyanobacteria was the second most abundant phylum with its proportion significantly reduced in all challenged groups as compared to the unchallenged negative control group (p < 0.0001). The third most abundant phylum was Proteobacteria whose proportion was 0.95% in the unchallenged control group. In contrast, the abundance of Proteobacteria increased in birds exposed to cpa+netB+tpeL+cpb2 (5.94%), cpa+netB+cpb2 (4.50%), or cpa+cpb2 (1.23%) (p < 0.2960, Fig 6B).

Jejunal relative microbial abundance at the family level

Microbial abundance at the family level in the unchallenged negative control group was dominated by Lactobacillaceae (61.99%) followed by Peptostreptococcaceae (16.02%), Oscillatoriophycideae (10.78%) and other families (5.40%) (i.e., including Clostridiaceae [1.21%] and Enterobacteriaceae [0.31%]) (Fig 6C). Birds exposed to CP isolates containing cpa+netB+tpeL+cpb2, cpa+netB+cpb2 and cpa+cpb2 had a significant decrease in Lactobacillaceae (17.42%, 18.37%, 13.02%, respectively) as compared to the unchallenged control group (p < 0.0001). Similarly, birds exposed to cpa+netB+tpeL+cpb2, cpa+netB+cpb2 and cpa+cpb2 had a decreased proportion of Oscillatoriophycideae (i.e., 3.10%, 7.23%, and 3.87%, respectively) compared to the unchallenged control group (p<0.3250). There was a significant increase in Clostridiaceae in birds exposed to cpa+netB+tpeL+cpb2 (p < 0.0001) compared to the unchallenged control group. An increase in the proportion of Clostridiaceae was also noted in cpa+netB+cpb2 (p = 0.3540) and cpa+cpb2 (p = 0.4400) compared to the unchallenged negative control group. Similarly, a significant increase in the abundance of Peptostreptococcaceae was observed in birds exposed to cpa+netB+cpb2 (54.32%) and cpa+cpb2 (71.80%) compared to birds not exposed to CP (16.02%, p < 0.0001) and in birds challenged with CP containing cpa+netB+tpeL+cpb2 (11.70%, p < 0.774) (Fig 6C).

Jejunal relative microbial abundance at the genus level

At the genus level, Lactobacillus, Oscillatoria, Romboutsia, Clostridium, Escherichia and Weissella were identified with a relative abundance of 57.31%, 18.81%, 16.40%, 1.55%, 0.38% and 0.16%, respectively in the unchallenged control group (Fig. 6D). In contrast, there was a significant increase in Clostridium (49.88%) in birds challenged with CP containing cpa + netB + tpeL + cpb2 compared to the negative control group (p < 0.0001). In the same group, there was an increase in the proportion of Escherichia (5.34%), Romboutsia (17.67%) and Weissella (0.23%), and a significant decrease in Lactobacillus (13.28%, p < 0.0001). The genus Oscillatoria was undetectable in the group. The birds exposed to the CP isolate containing cpa + netB + cpb2 genes had a significantly increased proportion of Romboutsia (53.32%) and a significantly reduced level of Lactobacillus (16.18%, p < 0.0001) as compared to the other groups (p < 0.0001). A rise in the proportion of Oscillatoria (14.62%, p < 0.4690), Clostridium (8.52%, p < 0.0760) and Escherichia (3.49%, p = 0.7050) were also observed. In addition, birds exposed to CP isolate containing cpa + cpb2 had Romboutsia (71.23%), Lactobacillus (11.79%), Clostridium (7.88%), Oscillatoria (3.30%), and Escherichia (0.96%) (Fig. 6D).

Alpha and beta-jejunal microbial diversity based on CP toxin genes

The microbial diversity within each sample (alpha diversity) of the group of birds infected with CP containing cpa, netB, cpb2 and tpeL genes (CP 4T), cpa, cpb2 and netB genes (CP 3T), cpa and cpb2 genes (CP 2T) or non-infected control birds was examined by computing richness (observed OTUs), Shannon index, Simpson index and Pielou’s evenness. Based on the distribution of each metric across the groups, a clear distinction in the microbial diversity within the groups (Fig. 7A) was observed. The differences in microbial composition between the CP 4T, CP 3T, CP 2T and non-infected control groups was evaluated by beta diversity analysis. The measurement of the weighted and unweighted UniFrac distances was visualized by PCoA (Fig. 7B). Unweighted UniFrac distances demonstrated statistically significant differences in microbial community structures between the CP 4T, CP 3T, CP 2T and non-infected control groups (p < 0.001; R2 = 0.4932; F-value = 21.085). In addition, weighted UniFrac distances (taxa abundance), between the CP 4T, CP 3T, CP 2T and non-infected control groups were statistically significant (p < 0.001; R2 = 0.81816; F-value = 97.485). The high R2 values suggest that the toxin types in CP isolates are significant contributors to the observed variations in the microbial communities.

Fig. 7
figure 7

Jejunal microbiome diversity and composition in control and CP toxin-exposed broiler chickens. A Alpha diversity indices with box plots showing distribution across control and different CP toxin combinations infected groups (n = 20/group, n = 10 for control). B PCoA plots of microbial community compositions, using unweighted and weighted UniFrac distances. Toxin labels: CP 4T includes cpa, netB, cpb2, tpeL genes; CP 3T has cpa, cpb2, netB genes; CP 2T contains cpa and cpb2 genes

Differential microbial abundance in NE affected birds.

DESeq2 results have discerned distinct bacterial abundance patterns tailored to specific toxin combinations. The Additional file 2: Fig. S1 and Additional file 3: Fig. S2 provide the statistical analysis of microbial abundance across different comparison groups. The 4T group was characterized by an overabundance of taxa such as Kurthia, Escherichia, and Staphylococcus, primarily within the Firmicutes and Proteobacteria phyla (p < 0.01) (Fig. 8A). In contrast, the 3T group exhibited a decreased abundance of bacterial classifications like Ralstonia, Eubacteriaceae, and Variovorax, mainly from the Proteobacteria and Firmicutes phyla (p < 0.01). Notably, bacterial taxa including Oscillatoriaceae, Clostridium, and Peptostreptococcaceae, primarily from the Cyanobacteria and Firmicutes phyla, maintained consistent abundance across all toxin combinations, inclusive of the 2T group (p ≤ 0.01).

Fig. 8
figure 8

Venn diagram illustrating differential abundance values in the microbiome of broiler chickens. This figure displays a Venn diagram heatmap showing LogFC values of microbiome abundance for five conditions: A “CP 4T”, “CP 3T”, “CP 2T”, and B “Clinical” and “Subclinical” conditions. Each cell in the heatmap represents a LogFC value, with non-significant values in white. The x-axis represents LogFC value in the clinical conditions, and the y-axis represents microbiome classifications

Further dissecting the clinical manifestations of NE, profound variations in bacterial abundance emerged. Birds with clinically severe symptoms showed a decline in several bacteria from the "Firmicutes" phylum, such as Clostridioides and Lysinibacillus, and also from the "Proteobacteria" phylum, like Shigella and Rhizobium (p < 0.01), (Fig. 8B). Conversely, birds manifesting subclinical conditions exhibited reduced abundances of bacteria primarily from the "Firmicutes" phylum (p < 0.01). Specifically, among the clinically severe birds, there was an elevated abundance of the ‘Clostridiales’ (p < 0.01).

Discussion

NE is a complex disease and a variety of predisposing factors like Eimeria sp. infections, immunosuppressive viral diseases such as IBDV, and poor management factors such as high stocking density, abrupt change in chicken feed ingredients and poor litter quality play a critical role in the occurrence and severity of the disease [34,35,36]. The degree of virulence of CP also plays a major role in the pathogenesis of NE. CP isolates expressing the plasmid encoded netB toxin have been associated with cases of NE [37], and the netB gene is infrequently detected in CP isolated from healthy chickens compared to CP isolates linked to NE outbreaks [18, 27, 38]. In addition, CP isolates expressing tpeL toxin (large clostridial toxin [LCT] family) [39] along with netB toxin produce a more severe form of NE [20]. Apart from netB and tpeL toxins, the role other CP toxins play in the pathogenesis of the NE is not well understood. Here, we investigated the potential interdependece of cpa, netB, cbp2 and tpeL toxins produced by CP in the pathogenesis of NE. To best of our knowledge it is the first study demonstrating the microbiome analysis based on toxin genes of different CP isolates in broiler chickens and their role in the development of NE. We aslo examined the jejunal microbial dysbiosis associated with subclinical and clinical NE induced by CP isolates containing dfferent combination of toxin genes.

Out of the 57 clinical isolates of CP characterized under this study, the majority contained cpa and netB toxin genes followed by CPA as a major toxin. A small minority of CP isolates contained tpeL gene. This is in line with previous reports indicating that NE outbreaks are mostly associated with CP expressing netB toxin [18, 27, 38]. Under this project, we successfully developed a subclinical and/or clinical NE broiler chicken animal model based on vIBDV-SK09 [40] as an immunosuppresant and/or by abrupt increase in the proportion of plant based protein in the chicken diet. In the past, Eimeria sp. infection and fish meal based NE animal models have been reported [31, 32, 41, 42]. Since in vitro growth and maintainance of Eimeria sp. is difficult and time consuming, we believe that our animal model is an easy fit for studying the pathogenesis of NE. Our results demonstrated that immunosupression with vIBDV-SK09 follwowed by CP infection produced increased mortality and severe gross and histolopathologic lesions suggesting immunosuppression as one of the most improtant predisposing factors to NE development. In addition, an abrupt increase in the concentration of protein in the broiler feed effectively precipiateted NE in CP exposed chickens. High concentration of protein in the feed tend to induce CP to express an increased amount of cpa toxin. Since CP lacks enzymes for amino acid biosynthesis, high protein diet also acts as a source of amino acids and a substrate for the proliferation of CP [5, 43].

In our animal model, both subclinical and clinical NE were consistently observed in the jejunum. Therefore, our microbiome based microbial diversity study was focused on jejunal samples collected from healthy, subclinical or clinical NE cases. In contrast, most past microbiome studies were focused on the cecum despite the fact that NE lesions occur in the small intestine [31, 32, 44, 45]. These studies also used exposure of birds to Eimeria species as a predisposing factor which could directly cause pathologic lesions in the intestine. To the best of our knowledge, this is the first study which demonstrated changes in the patterns of the jejunal microbiome at the different stages of NE development post CP infection. In this study, Lactobacillus (57.11%), Oscillatoria (18.81%), Romboutsia (16.40%) and Clostridium (1.35%) were identified as the common bacterial genera consistently identified in the jejunum of healthy broiler chickens. The dominance of Lactobacillus in the jejunal microbial community of healthy chickens was also reported in previous studies [46]. Birds that developed subclinical or clinical NE after exposure to CP containing different combinations of toxin genes had significantly reduced proportion of Lactobacillus and significantly increased level of Clostridium. Similar phenomenon were reported in the ileum [47] and jejunum [46] of birds after experimental induction of NE. The competitive exclusion of Lactobacillus might be essential for CP’s successful colonization of the jejunum and induction of the development of NE. The genus Lactobacillus is known for its probiotic properties and maintain intestinal health by regulating the intestinal microbiota [48] and by modulating the immune system [49] Lactobacillus can also produce bacteriostatic bacteriocin-like compound against pathogenic bacteria [50]. Moreover, the bacteria produces lactic acid upon fermentation of carbohydrates and reduces the intestinal pH making an inhabitable environment for acid sensitive pathogens [51, 52]. Hence, suppression of Lactobacillus was considered beneficial for the colonization and growth of enteric pathogens [46]. Reduction in the proportion of other lactic acid producing bacteria like the genus Weissella were reported in the jejunum and ceca of birds that developed NE [32, 46]. Interestingly, the level of E. coli in the jejunum increased dramatically in the subclinical NE cases as compared to the unchallenged group. The proportion of the genus Escherichia continued to rise significantly as the CP challenged birds developed clinical NE which are consistent with previous reports [53, 54]. Similar enrichement of Escherichia coli was also reported in the ileum and cecum of mice after infection with porcine CP [55]. It was proposed that E. coli may counteract the serum endotoxin secreted by CP [56]. E. coli and CP might also have a synergestic effect in the pathogenesis of NE. But further studies are required to clearly elucidate the role E. coli plays in the development of NE after CP infection.

In addition, our study indicates that the change in the microbial diversity and reduction in the proportion of Lactobacillus sp. was dependent upon the CP toxins. The presence of netB and tpeL toxin genes caused the most dysbiosis as compared to a combination of cpa and cpb2 toxin genes alone. A distinct pattern in the jejunal microbial population was observed between unchallenged negative control birds and birds that developed subclinical, or clinical NE after challenge with CP containing cpa+netB+tpeL+cpb2 genes. This indicates that dysbiosis of the normal intestinal microbial community is an essential step in the pathogenesis of NE. Only a few broilers exposed to CP isolates containing either cpa+netB+cpb2 or cpa+cpb2 developed NE. Interestingly, birds which did not develop either subclinical or clinical NE had a change in the jejunal microbial community as compared to the negative control group (data not shown). More specifically, these birds had a statstically significant increase in the abundance of the genus Rombutsia as compared to the challeneged groups and the unchallenged control groups. This suggests that the genus Rombutsia may contribute to the inhibition of the proliferaiton of CP and progression of the infection to subclinical or clinical NE. Previous studies allude that Rombutsia metabolizes monosaccharides and disacharides and produces acetic, lactic and fomic acids [57]. Buffered formic acid have shown to improve the the performace of broiler chickens with improved feed conversion ratio [58]. Likewise, supplementation of Astragalus polysaccharide to NE affected birds led to the proliferation of Romboutsia in the ileum and was correlated with a decreased abundance of CP in the cecum [59]. Altogether, the genus Rombutsia might play a significant role in the competitive exclusion of CP and helps in the prevention/control of NE during exposure to CP containing cpa+netB+cpb2 or cpa+cpb2 genes. Our results provide evidence that infection of birds with CP containing cpa+netB+tpeL+cpb2 genes can effectively dampen the proliferation Romboutsia and Lactobacillus and produce NE.

Beta diversity analyses, comparing microbial composition across groups, revealed significant variations linked to NE severity, influenced by toxin gene combinations. CP 3T (cpa, cpb2, netB genes) and CP 2T (cpa and cpb2 genes) exhibit unique beta diversity patterns, shedding light on the interplay between bacterial composition, toxin genes, and NE development in broiler chickens.

Differential abundance analysis of the microbiome in broiler chickens unveils the intricate microbial dynamics underlying disease pathology. In this study, DESeq analysis revealed a stable presence of certain bacterial groups across varying toxin combinations may indicate their inherent resilience or pivotal role in the context of NE in chickens. Intriguingly, both clinical conditions showcased decreased abundance of certain bacteria from the "Firmicutes" and "Proteobacteria" phyla, albeit in distinct ratios. This unveils the nuanced microbial shifts associated with varied clinical presentations of NE, deepening our understanding of the disease's pathology. Collectively, these findings emphasize the intricate dynamics between bacterial abundance and its potential implications in shaping the progression and pathogenesis of NE.

One of the limitations of this study is that the microbiome study was based on 16S rRNA gene sequencing and the microbial classififcation was only possible upto the genus level. Whole genome metagenomic studies will need to be performed in the future for microbiome analysis at the species level. The significance of the chicken intestinal mycobiome and virome in the pathogenesis of NE should also be addressed in the future. Our study also focused on the jejunum as it is the prinicpal site for NE development. But the impact of NE in the microbial community present in the different parts of the small and large intestine needs to be investigated. As E. coli was found in high proportions in CP infected birds which developed subclinical or clinical NE, the role E. coli plays in the pathogenesis of NE requires further investigation. In addition, the specific Rombutsia and Lactobacillus species which maintain the health of the jejunum of birds should be identified and tested as potential probiotics for the control and prevention of NE.

In summary, our first objective successfully established a broiler chicken model for NE through an abrupt increase in plant-based protein in the diet and the induction of immunosuppression via vIBDV-SK09 infection. CP isolates carrying the cpa+netB+tpeL+cpb2 genes induced microscopic lesions. Our second objective consistently revealed distinct differences in the jejunal microbial community between CP-infected birds and healthy ones, serving as a potential marker for identifying subclinical and clinical NE in broiler chickens. The proliferation of E. coli in the jejunum of CP-infected birds may play a crucial role in exacerbating NE pathogenesis, although further confirmation is needed. Additionally, our third objective indicates the importance of netB and tpeL genes in NE pathogenesis, with up to 90% of CP-infected birds carrying cpa+netB+tpeL+cpb2 genes developing microscopic lesions in the jejunum and significantly impacting jejunal microbial diversity. As the poultry industry moves away from prophylactic antimicrobial use, the resurgence of NE is a growing concern. Our study provides evidence that the Romboutsia and Lactobacillus genera could be promising candidates for developing probiotic products as effective alternatives for preventing and controlling NE.

Materials and methods

CP isolates, and virotyping by polymerase chain reaction (PCR)

Fifty-seven CP isolates were characterized by polymerase chain reaction (PCR) for the presence of toxin genes. The 50 isolates were kindly provided by Dr. Patrick Boerlin, Department of Pathobiology, Ontario Veterinary College (OVC), University of Guelph, Canada. The remaining isolates were provided by Prairie Diagnostic Services (PDS) Inc., University of Saskatchewan. All the CP were isolated from field cases of NE in chickens across the provinces of Alberta, Saskatchewan and Ontario (unpublished). Each CP isolate was individually cultured in 3 mL of brain heart infusion (BHI) broth incubated at 37 °C for 18 hr under anaerobic conditions. Total genomic DNA was extracted using DNeasy Blood and Tissue kit (Qiagen) according to the manufacturer’s instructions. The quality and quantity of the extracted genomic DNA was determined by using the NanoDrop™ spectrophotometer (Thermo Scientific™). The targeted toxin genes of CP (i.e., major toxin genes; including cpa, cpb, etx and iap, as well as minor toxin genes including cpe, netB, cpb2 and tpeL) were examined by PCR as previously described [60]. Briefly, all the PCR reactions were carried out in a final reaction volume of 50 µL using a template DNA of 100 ng and 2 units of Taq DNA polymerase (Thermo Fisher Scientific). The PCR program consisted of 35 cycles of denaturation (95°C for 30 sec), annealing (To varied depending on the gene amplified and performed for 30 sec), and elongation (72°C for 1 min). The primer sequences, annealing To and fragment size of each toxin gene amplicon is given in Additional file 1: Table S1. Reference CP strains (ATCC 13124, Clostridium perfringens; Strain S 107, contain cpa gene; ATCC 3626, Clostridium perfringens; type b, contains cpa, cpb, etx, cpb2, and pfoA genes and an in house CP 21 strains which contains cpa, netB, cpb2, and tpeL genes) were included as controls. The amplified PCR products were analyzed by gel electrophoresis.

Preparation of CP for broiler chicken challenge

Clinical isolates of CP containing cpa, netB, cpb2 and tpeL (toxin type G) toxin genes were selected from our culture collection. To prepare the challenge inoculum, each CP isolate was streaked on separate 5% Columbia sheep blood agar plates in duplicate (Thermo Scientific, Canada) and incubated at 37 ºC for 24 hr, under anaerobic conditions. Single colonies were individually transferred to cooked meat broth media (CMM, Sigma-Aldrich, Canada) and incubated at 37 ºC for 24 hr, under anaerobic conditions. 10 mL of cultured CMM broth was diluted with 490 mL of CMM (dilution factor 50; v/v) and incubated for 24 hr at 37 °C. Subsequently, the CMM culture was used to inoculate fluid thioglycollate media (FTG, Sigma-Aldrich, Canada) (3% [v/v]) and incubated anaerobically at 37 ºC for 15–16 h and 1×109 colony forming units (CFU) per mL was used to challenge broilers at 20, 21 and 22 days of age.

NE animal model development protocol

Day old broiler chickens (n=120, Fraser Valley Chick Sales, Abbotsford, BC) were placed in the animal care unit (ACU)’s Animal Biosafety level (ABSL)-2 facility, Western College of Veterinary Medicine (WCVM), University of Saskatchewan, Canada. The birds were randomly allocated to 6 groups (n=20/group) as follows: (1) no vIBDV-SK09 and no CP challenge; (2) CP challenge in-feed with 20% protein ; (3) CP challenge in-feed with 28% protein; (4) vIBDV-SK09 challenge and CP challenge in-feed with 20% protein; (5) vIBDV-SK09 challenge and CP challenge in-feed with 28% protein; (6) vIBDV-SK09 challenge and CP challenge orally (1 mL of CP culture twice daily for three days) with 28% protein feed. Each group was kept in a separate pen in the same room. Induction of immunosuppression after infection with vIBDV-SK09 followed by an abrupt increase in protein content in the feed was used to predispose broiler chickens to CP infection and NE development. The CP isolate used for this experiment contained cpa, netB, cpb2 and tpeL genes. Until 18 days of age, the diet of day-old broiler chickens contained 20% of raised without antibiotics (RWA) poultry starter (Farm ChoiceTM RWA, MasterFeeds, Canada). Birds were then orally infected with 1x103 TCID50/bird (50% tissue culture infectious dose) of vIBDV-SK09 [40] at 17 days of age to induce immunosuppression. vIBDV-SK09 was used since it causes severe immunosuppression by depleting 90% of B-cells in the peripheral blood [61]. Later, feed was withdrawn at 19 days of age for 10-12 hr followed by feeding with 28% RWA feed. The 28% RWA feed was prepared by mixing 25% RWA turkey/gamebird starter crumble (MasterFeeds, Canada, Canada) with 38% layer/grower supplement (MasterFeeds, Canada) at a 10:3 ratio. Broilers were fed with the 28% protein ration with added FTG-grown CP culture in 1:1 (v/w) ratio for 3 consecutive days (20 to 22 days of age). Clinical signs were scored three times per day as follows; 0 = normal; 1= ruffled feathers, reluctant to stand, may peck some feed, may have loose droppings and depressed but still responsive to the environment; 2 = ruffled feathers, sitting posture, hesitates to move, may show brown foamy diarrhea, do not forage, severely depressed; 3 = bird found dead without any pre-monetary signs. Birds were euthanized at 23 days of age and post-mortem examinations were conducted. A diagramatic representation of the experimental procedures is provided (Fig 1). All the animal experiments were approved by the Animal Research Ethics Board, University of Saskatchewan and the experimental procedures were adhered to the Canadian Council on Animal Care guidelines for humane animal use.

Sample collection protocol for microbiome study

At 23 days of age, the broilers were euthanized, and post-mortem examination conducted. Gross lesions were recorded, and sections of the duodenum, jejunum and ileum were collected for histopathology as described earlier. At the end of the trial, contents of the jejunum (n=8/group) were collected to study the effect of CP challenge on the intestinal microbiome. Samples from broilers with no CP challenge, broilers with macroscopic lesions of NE, and broilers with only microscopic lesions of NE were investigated for the composition of the jejunal microbiome.

Animal challenge experiments with CP isolates containing different toxin genes

To investigate the role of different toxin combinations in the pathogenesis of NE, challenge experiments were conducted using different CP isolates on the animal-infection model described above. A total of 80 day-old broiler chicks (Prairie Pride, SK, Canada) were randomly divided into four groups (n=20/group): (1) no CP challenge; (2) CP strain containing cpa, netB, cpb2 and tpeL genes; (3) a CP strain containing cpa, cpb2 and netB genes; and (4) CP strain containing cpa and cpb2 genes. The CP culture preparation and animal challenge were similar as mentioned before. Broilers were monitored for development of clinical signs and mortality three times per day until the end of the trial. The experiment was terminated at 3 days post challenge and necropsy was performed. Gross lesions were recorded, and sections of intestines were collected for histopathology. The jejunal contents were collected and stored at -80 °C freezer for DNA extraction and 16S amplicon sequencing.

Gross and histopathological scoring of the intestine

All dead and euthanized broilers were necropsied (n=20/group) and intestines were carefully opened to record the gross NE lesions. The duodenum, jejunum, and ileum of broilers with and without gross NE lesions were collected in 10% neutral buffered formalin (pH=7.4) for histologic examination. Processing of the fixed intsetinal tissues and hematoxylin and eosin (H & E) staining was performed as described previously [62]. Briefly, fixed intestines were transferred to graded alcohol to xylene solutions for dehydration, followed by clearing. The tissues were then embedded in paraffin wax for sectioning. The tissue sections with a thickness of 3-5 μm were stained with conventional H & E and examined under a light microscope. Microscopic lesions were scored as: 0 = no lesions/ healthy mucosa; 1 = focal necrosis of intestinal villi, acute; 2 = necrosis of intestinal villi, multifocal to coalescing, acute; 3 = diffuse, necrosis of intestinal villi, acute, severe.

Genomic DNA extraction and 16S amplicon sequencing

Total DNA was extracted from jejunal contents (200 mg/sample) using a QIAamp® Fast DNA Stool Mini kit (Qiagen Inc.). The DNA concentration was measured by a NanoDrop spectrometer (Thermo Scientific™, USA). The 16S amplicon PCR (targeting the V3-4 region) was performed as per the Illumina 16S metagenomic sequencing library preparation protocol. The PCR products were purified using AMPure XP magnetic beads (Beckman Coulter, Inc, CA, USA). The size of the amplicons was confirmed by gel electrophoresis. Sequence libraries were prepared using the Nextera XT DNA Library Preparation Kit (Illumina, San Diego, USA) with 24 or 96 indexes as per the manufacturer’s instructions. The libraries were cleaned up with AMPure XP magnetic beads and analyzed by a bioanalyzer using a high sensitivity DNA kit (Agilent, Santa Clara, CA). The concentration of the amplicon sequence libraries was measured using the Qubit dsDNA BR (broad range) assay kit (Thermo Fisher Scientific, Waltham, USA). The normalized libraries were pooled, denatured with 0.2N NaOH, diluted using H1 hybridization buffer and loaded onto the MiSeq V3 (600 cycles) cartridge. Sequencing was performed using the Illumina MiSeq platform (Illumina, San Diego, US) as per the company’s protocol. The sequence quality was examined in real time using the Sequencing Analysis Viewer (SAV) V1.8.

Amplicon sequence assembly and taxonomic classification

The FASTQ paired sequence data was extracted from the MiSeq and the data quality metrices were examined by the Illumina SAV V1.8 software. The FASTQ files were then imported into Geneious Prime 2022.01 software (Dotmatics Inc.) and microbiome analysis was performed using the 16S Biodiversity analysis tool. The paired end reads were quality trimmed using BBDuk plugin with the minimum quality score and minimum length set at 30 and 100, respectively. The trimmed sequences were merged using the BBMerge tool. The contigs were then assembled and clustered into operational taxonomic units (OTUs) using the Geneious de novo assembler. The representative sequences were BLASTed to a curated 16S rRNA gene reference sequence database from NCBI to assign taxonomy. Finally, the BLAST results were used as a targeted taxonomic database to quantify the biodiversity in the full data set using the Geneious Prime sequence classifier plugin.

Microbiome diversity and abundance analysis

The jejunal microbial community composition data collated from both data sets were used to perform microbiome diversity analysis. The first dataset contained the relative abundance of microbial taxa for three different groups (n = 8/group; total = 24 samples) which included birds that developed subclinical, and clinical NE post CP infection and non-infected control group. The second dataset contained the relative abundance of jejunal microbial taxa for 4 different groups (20 birds/group except group 1 [n = 10]) including birds which were infected with a CP strain containing cpa, netB, cpb2 and tpeL genes (CP 4 T), a CP strain containing cpa, cpb2 and netB genes (CP 3 T), CP strain containing cpa and cpb2 genes (CP 2 T), and a non-infected control group. To evaluate the diversity within individual samples (alpha diversity), we calculated metrics including richness (observed OTUs), the Shannon index, the Simpson index, and Pielou’s evenness. This analysis was performed using the skbio.diversity.alpha module within the scikit-bio library in Python [63]. The differences in microbial composition between groups (beta diversity) were analyzed by computing the weighted and unweighted UniFrac Distances. These phylogenetic distance metrics were computed using the phyloseq and ape packages in R. Using the computed beta diversity metrices, PCoA was conducted to visualize the differences in microbial community structures among the samples in two-dimensional space. Permutational Multivariate Analysis of Variance (PERMANOVA) was employed to statistically test the differences in microbial compositions between the groups. This was executed using the `adonis2` function from the `vegan` package in R [64]. The test was applied on both unweighted (presence/absence) and weighted (abundance) UniFrac distances. Additionally, DESeq2 analysis was performed to identify taxa that exhibited significant differential abundance between experimental conditions or groups in microbiome datasets derived from sequencing data. This analysis was performed using the DESeq2 package in conjunction with the phyloseq [65] package in the R program.

Statistical analysis

Survival, gross and histopathological data were analyzed using GraphPad Prism V6.0 software (Dotmatics Inc.) with a significance level of p < 0.05. Survival or mortality data were compared between the groups using the long-rank (Mantel-cox) and chi-square tests, respectively. Microscopic lesions of the intestine were analyzed between the groups using non-parametric Kruskal-Wallis one-way ANOVA and Tukey’s multiple comparison (p < 0.0001). Further comparison of microscopic lesion scores between two groups was analysed with two-tailed Mann Whitney test (p < 0.0001). Family, phylum, and genus level relative microbial diversity were analyzed by Ordinary two-way ANOVA (p < 0.0001). Further, Turkey’s multiple comparison was performed to identify variation between the different experimental groups at 95% confidence interval. Individual and group sample variance at the family level was examined by principal component analysis using GraphPad Prism V6.0 software (Dotmatics Inc.). Alpha diversity metrics across test groups were assessed using the Kruskal-Wallis and Dunn’s post hoc tests. Beta diversity was analyzed using weighted and unweighted UniFrac distances. Differential abundance, based on DESeq2 logFC values (± 1.5 fold), was tested using the Wald test. A p-value of p ≤ 0.01 was considered statistically significant.

Availability of data and materials

The data supporting the findings of this study are available from the authors upon request.

References

  1. Hatheway CL. Toxigenic clostridia. Clin Microbiol Rev. 1990;3:66–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Rood JI, Adams V, Lacey J, Lyras D, McClane BA, Melville SB, et al. Expansion of the Clostridium perfringens toxin-based typing scheme. Anaerobe. 2018;53:5–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Cooper KK, Songer JG. Virulence of Clostridium perfringens in an experimental model of poultry necrotic enteritis. Vet Microbiol. 2010;142:323–8.

    Article  PubMed  Google Scholar 

  4. Kaldhusdal M, Hofshagen M. Barley inclusion and avoparcin supplementation in broiler diets.: 2. clinical, pathological, and bacteriological findings in a mild form of necrotic enteritis. Poult Sci. 1992;71:1145–53.

    Article  CAS  PubMed  Google Scholar 

  5. McDevitt R, Brooker J, Acamovic T, Sparks N. Necrotic enteritis; a continuing challenge for the poultry industry. Worlds Poult Sci J. 2006;62:221–47.

    Article  Google Scholar 

  6. Kaldhusdal M, Hofshagen M, Løvland A, Langstrand H, Redhead K. Necrotic enteritis challenge models with broiler chickens raised on litter: evaluation of preconditions, Clostridium perfringens strains and outcome variables. FEMS Immunol Med Microbiol. 1999;24:337–43.

    Article  CAS  PubMed  Google Scholar 

  7. Shojadoost B, Vince AR, Prescott JF. The successful experimental induction of necrotic enteritis in chickens by Clostridium perfringens: a critical review. Vet Res. 2012;43:74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Casewell M, Friis C, Marco E, McMullin P, Phillips I. The European ban on growth-promoting antibiotics and emerging consequences for human and animal health. J Antimicrob Chemother. 2003;52:159–61.

    Article  CAS  PubMed  Google Scholar 

  9. Pulcini C, Bush K, Craig WA, Frimodt-Møller N, Grayson ML, Mouton JW, et al. Forgotten antibiotics: an inventory in Europe, the United States, Canada, and Australia. Clin Infect Dis. 2012;54:268–74.

    Article  PubMed  Google Scholar 

  10. Fathima S, Hakeem W, Shanmugasundaram R, Selvaraj R. Necrotic enteritis in broiler chickens: a review on the pathogen, pathogenesis, and prevention. Microorganisms. 2022;10(10):1958.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Bazacliu C, Neu J. Necrotizing enterocolitis: long term complications. Curr Pediatr Rev. 2019;15:115–24.

    Article  PubMed  Google Scholar 

  12. Bazacliu C, Neu J. Pathophysiology of necrotizing enterocolitis: an update. Curr Pediatr Rev. 2019;15:68–87.

    Article  PubMed  Google Scholar 

  13. Deblais L, Kathayat D, Helmy YA, Closs G, Rajashekara G. Translating ‘big data’: better understanding of host-pathogen interactions to control bacterial foodborne pathogens in poultry. Anim Health Res Rev. 2020;21:15–35.

    Article  PubMed  Google Scholar 

  14. Beacon TH, Davie JR. The chicken model organism for epigenomic research. Genome. 2021;64:476–89.

    Article  CAS  PubMed  Google Scholar 

  15. Uzal FA, McClane BA, Cheung JK, Theoret J, Garcia JP, Moore RJ, et al. Animal models to study the pathogenesis of human and animal Clostridium perfringens infections. Vet Microbiol. 2015;179:23–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Garcia P, Wang Y, Viallet J, Macek JZ. The chicken embryo model: a novel and relevant model for immune-based studies. Front Immunol. 2021;12: 791081.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Keyburn AL, Sheedy SA, Ford ME, Williamson MM, Awad MM, Rood JI, et al. Alpha-toxin of Clostridium perfringens is not an essential virulence factor in necrotic enteritis in chickens. Infect Immun. 2006;74:6496–500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Keyburn AL, Boyce JD, Vaz P, Bannam TL, Ford ME, Parker D, et al. NetB, a new toxin that is associated with avian necrotic enteritis caused by Clostridium perfringens. PLoS Pathog. 2008;4: e26.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Chalmers G, Bruce H, Hunter D, Parreira V, Kulkarni R, Jiang Y-F, et al. Multilocus sequence typing analysis of Clostridium perfringens isolates from necrotic enteritis outbreaks in broiler chicken populations. J Clin Microbiol. 2008;46:3957–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Coursodon CF, Glock RD, Moore KL, Cooper KK, Songer JG. TpeL-producing strains of Clostridium perfringens type A are highly virulent for broiler chicks. Anaerobe. 2012;18:117–21.

    Article  CAS  PubMed  Google Scholar 

  21. Oda M, Terao Y, Sakurai J, Nagahama M. Membrane-binding mechanism of Clostridium perfringens alpha-toxin. Toxins. 2015;7:5268–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Oda M, Matsuno T, Shiihara R, Ochi S, Yamauchi R, Saito Y, et al. The relationship between the metabolism of sphingomyelin species and the hemolysis of sheep erythrocytes induced by Clostridium perfringens α-toxins. J Lipid Res. 2008;49:1039–47.

    Article  CAS  PubMed  Google Scholar 

  23. Navarro MA, McClane BA, Uzal FA. Mechanisms of action and cell death associated with Clostridium perfringens toxins. Toxins. 2018;10:212.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lepp D, Roxas B, Parreira VR, Marri PR, Rosey EL, Gong J, et al. Identification of novel pathogenicity loci in Clostridium perfringens strains that cause avian necrotic enteritis. PLoS ONE. 2010;5: e10795.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Prescott JF, Parreira VR, Mehdizadeh Gohari I, Lepp D, Gong J. The pathogenesis of necrotic enteritis in chickens: what we know and what we need to know: a review. Avian Pathol. 2016;45:288–94.

    Article  CAS  PubMed  Google Scholar 

  26. Rood JI, Keyburn AL, Moore RJ. NetB and necrotic enteritis: the hole movable story. Avian Pathol. 2016;45:295–301.

    Article  CAS  PubMed  Google Scholar 

  27. Martin TG, Smyth JA. Prevalence of netB among some clinical isolates of Clostridium perfringens from animals in the United States. Vet Microbiol. 2009;136:202–5.

    Article  CAS  PubMed  Google Scholar 

  28. Gu C, Lillehoj HS, Sun Z, Lee Y, Zhao H, Xianyu Z, et al. Characterization of virulent netB+/tpeL+ Clostridium perfringens strains from necrotic enteritis—affected broiler chicken farms. Avian Dis. 2019;63:461–7.

    Article  PubMed  Google Scholar 

  29. Guttenberg G, Hornei S, Jank T, Schwan C, Lü W, Einsle O, et al. Molecular characteristics of Clostridium perfringens TpeL toxin and consequences of mono-O-GlcNAcylation of Ras in living cells. J Biol Chem. 2012;287:24929–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Xiao Y, Xiang Y, Zhou W, Chen J, Li K, Yang H. Microbial community mapping in intestinal tract of broiler chicken. Poult Sci. 2017;96:1387–93.

    Article  CAS  PubMed  Google Scholar 

  31. Stanley D, Keyburn AL, Denman SE, Moore RJ. Changes in the caecal microflora of chickens following Clostridium perfringens challenge to induce necrotic enteritis. Vet Microbiol. 2012;159:155–62.

    Article  PubMed  Google Scholar 

  32. Stanley D, Wu S-B, Rodgers N, Swick RA, Moore RJ. Differential responses of cecal microbiota to fishmeal, Eimeria and Clostridium perfringens in a necrotic enteritis challenge model in chickens. PLoS ONE. 2014;9: e104739.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Macdonald SE, Nolan MJ, Harman K, Boulton K, Hume DA, Tomley FM, et al. Effects of Eimeria tenella infection on chicken caecal microbiome diversity, exploring variation associated with severity of pathology. PLoS ONE. 2017;12: e0184890.

    Article  PubMed  PubMed Central  Google Scholar 

  34. McReynolds J, Byrd J, Anderson R, Moore R, Edrington T, Genovese K, et al. Evaluation of immunosuppressants and dietary mechanisms in an experimental disease model for necrotic enteritis. Poult Sci. 2004;83:1948–52.

    Article  CAS  PubMed  Google Scholar 

  35. Timbermont L, Haesebrouck F, Ducatelle R, Van Immerseel F. Necrotic enteritis in broilers: an updated review on the pathogenesis. Avian Pathol. 2011;40:341–7.

    Article  CAS  PubMed  Google Scholar 

  36. Tsiouris V, Georgopoulou I, Batzios C, Pappaioannou N, Ducatelle R, Fortomaris P. High stocking density as a predisposing factor for necrotic enteritis in broiler chicks. Avian Pathol. 2015;44:59–66.

    Article  CAS  PubMed  Google Scholar 

  37. Keyburn AL, Bannam TL, Moore RJ, Rood JI. NetB, a pore-forming toxin from necrotic enteritis strains of Clostridium perfringens. Toxins. 2010;2:1913–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Tolooe A, Shojadoost B, Peighambar S, Tamaddon Y. Prevalence of netB gene among Clostridium perfringens isolates obtained from healthy and diseased chickens. J Anim Vet Adv. 2011;10:106–10.

    Article  CAS  Google Scholar 

  39. Nakano V, Ignacio A, Llanco L, Bueris V, Sircili MP, Avila-Campos M. Multilocus sequence typing analyses of Clostridium perfringens type A strains harboring tpeL and netB genes. Anaerobe. 2017;44:99–105.

    Article  CAS  PubMed  Google Scholar 

  40. Kurukulsuriya S, Ahmed KA, Ojkic D, Gunawardana T, Gupta A, Goonewardene K, et al. Circulating strains of variant infectious bursal disease virus may pose a challenge for antibiotic-free chicken farming in Canada. Res Vet Sci. 2016;108:54–9.

    Article  PubMed  Google Scholar 

  41. Lu M, Li RW, Zhao H, Yan X, Lillehoj HS, Sun Z, et al. Effects of Eimeria maxima and Clostridium perfringens infections on cecal microbial composition and the possible correlation with body weight gain in broiler chickens. Res Vet Sci. 2020;132:142–9.

    Article  CAS  PubMed  Google Scholar 

  42. Wu S-B, Stanley D, Rodgers N, Swick RA, Moore RJ. Two necrotic enteritis predisposing factors, dietary fishmeal and Eimeria infection, induce large changes in the caecal microbiota of broiler chickens. Vet Microbiol. 2014;169:188–97.

    Article  PubMed  Google Scholar 

  43. Lan Y, Verstegen M, Tamminga S, Williams B. The role of the commensal gut microbial community in broiler chickens. World’s Poult Sci J. 2005;61:95–104.

    Article  Google Scholar 

  44. Akerele G, Al Hakeem WG, Lourenco J, Selvaraj RK. The effect of necrotic enteritis challenge on production performance, cecal microbiome, and cecal tonsil transcriptome in broilers. Pathogens. 2022;11:839.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Tang Y, Zhang X, Wang Y, Guo Y, Zhu P, Li G, et al. Dietary ellagic acid ameliorated Clostridium perfringens-induced subclinical necrotic enteritis in broilers via regulating inflammation and cecal microbiota. J Anim Sci Biotechnol. 2022;13:1–18.

    Google Scholar 

  46. Yang W-Y, Lee Y, Lu H, Chou C-H, Wang C. Analysis of gut microbiota and the effect of lauric acid against necrotic enteritis in Clostridium perfringens and Eimeria side-by-side challenge model. PLoS ONE. 2019;14: e0205784.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Yang Q, Liu J, Wang X, Robinson K, Whitmore MA, Stewart SN, et al. Identification of an intestinal microbiota signature associated with the severity of necrotic enteritis. Front Microbiol. 2021;12: 703693.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Fusco V, Quero GM, Cho G-S, Kabisch J, Meske D, Neve H, et al. The genus Weissella: taxonomy, ecology and biotechnological potential. Front Microbiol. 2015;6:155.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Lin W, Yu B, Lin C, Hwang W, Tsen H. Immune effect of heat-killed multistrain of Lactobacillus acidophilus against Salmonella typhimurium invasion to mice. J Appl Microbiol. 2007;102:22–31.

    Article  CAS  PubMed  Google Scholar 

  50. Chateau N, Castellanos I, Deschamps A. Distribution of pathogen inhibition in the Lactobacillus isolates of a commercial probiotic consortium. J Appl Bacteriol. 1993;74:36–40.

    Article  CAS  PubMed  Google Scholar 

  51. Belenguer A, Duncan SH, Holtrop G, Anderson SE, Lobley GE, Flint HJ. Impact of pH on lactate formation and utilization by human fecal microbial communities. Appl Environ Microbiol. 2007;73:6526–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Sengupta R, Altermann E, Anderson RC, McNabb WC, Moughan PJ, Roy NC. The role of cell surface architecture of lactobacilli in host-microbe interactions in the gastrointestinal tract. Mediators Inflamm. 2013;2013: 237921.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Du E, Gan L, Li Z, Wang W, Liu D, Guo Y. In vitro antibacterial activity of thymol and carvacrol and their effects on broiler chickens challenged with Clostridium perfringens. J Anim Sci Biotechnol. 2015;6:1–12.

    Article  CAS  Google Scholar 

  54. Du E, Wang W, Gan L, Li Z, Guo S, Guo Y. Effects of thymol and carvacrol supplementation on intestinal integrity and immune responses of broiler chickens challenged with Clostridium perfringens. J Anim Sci Biotechnol. 2016;7:1–10.

    Article  Google Scholar 

  55. Jiang Z, Su W, Wen C, Li W, Zhang Y, Gong T, et al. Effect of porcine Clostridium perfringens on intestinal barrier, immunity, and quantitative analysis of intestinal bacterial communities in mice. Front Vetx Sci. 2022;9: 881878.

    Article  Google Scholar 

  56. Li Z, Wang W, Liu D, Guo Y. Effects of Lactobacillus acidophilus on gut microbiota composition in broilers challenged with Clostridium perfringens. PLoS ONE. 2017;12: e0188634.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Gerritsen J, Fuentes S, Grievink W, van Niftrik L, Tindall BJ, Timmerman HM, et al. Characterization of Romboutsia ilealis gen. nov., sp. nov., isolated from the gastro-intestinal tract of a rat, and proposal for the reclassification of five closely related members of the genus Clostridium into the genera Romboutsia gen. nov., Intestinibacter gen. nov., Terrisporobacter gen. nov. and Asaccharospora gen. nov. Int J Syst Evol Microbiol. 2014;64:1600–16.

    Article  CAS  PubMed  Google Scholar 

  58. Gharib-Naseri K, Kheravii SK, Li L, Wu S-B. Buffered formic acid and a monoglyceride blend coordinately alleviate subclinical necrotic enteritis impact in broiler chickens. Poult Sci. 2021;100: 101214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Song B, Li P, Yan S, Liu Y, Gao M, Lv H, et al. Effects of dietary astragalus polysaccharide supplementation on the Th17/Treg balance and the gut microbiota of broiler chickens challenged with Necrotic Enteritis. Front Immunol. 2022;13: 781934.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Timbermont L, Lanckriet A, Gholamiandehkordi AR, Pasmans F, Martel A, Haesebrouck F, et al. Origin of Clostridium perfringens isolates determines the ability to induce necrotic enteritis in broilers. Comp Immunol Microbiol Infect Dis. 2009;32:503–12.

    Article  PubMed  Google Scholar 

  61. Kurukulasuriya S, Ahmed KA, Ojkic D, Gunawardana T, Goonewardene K, Gupta A, et al. Evaluation of five circulating strains of variant infectious bursal disease virus (varIBDV) for their immunogenicity as broiler breeder vaccines and protective efficacy in neonatal broiler chicks. Vaccine. 2022;40:5608–14.

    Article  CAS  PubMed  Google Scholar 

  62. Feldman AT, Wolfe D. Tissue processing and hematoxylin and eosin staining. Histopathology. 2014;1180:31–43.

    Article  CAS  Google Scholar 

  63. Virtanen P, Gommers R, Oliphant TE, Haberland M, Reddy T, Cournapeau D, et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods. 2020;17:261–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Costantini MS, Medeiros MC, Crampton LH, Reed FA. Wild gut microbiomes reveal individuals, species, and location as drivers of variation in two critically endangered Hawaiian honeycreepers. PeerJ. 2021;9: e12291.

    Article  PubMed  PubMed Central  Google Scholar 

  65. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8: e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are thankful to Dr. Philip Griebel (Vaccinology and Immunotherapy, School of Public Health, University of Saskatchewan) and Dr. Ngeleka Musangu (Prairie Diagnistic Services Inc., University of Saskatchewan) for their expertise in all aspects of our study and for their help in reviewing the manuscript. The authors are grateful to the animal health technicians at the Animal Care Unit, Western College of Veterinary Medicine, University of Saskatchewan. Authors sincerely thank Dr. Babajan Banagapalli (Faculty of Medicine, King Abdulaziz University, Saudi Arabia) for his commendable help in computational programming steps and in data visualization.

Funding

Financial support was provided by grants from the Chicken Farmers of Saskatchewan and Poultry ScienceCluster III (Canadian Poultry Research Council and Agriculture and Agri-Food Canada).

Author information

Authors and Affiliations

Authors

Contributions

HG conducted the acquisition, analysis, interpretation of the data and drafted the manuscript under supervision of LA and SG. HG, LA, IS, SP, BCL, SG, AD and KAA helped in acquiring the experimental material and supported the animal experiments. HG, LA, NS verified the underlying data and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Susantha Gomis.

Ethics declarations

Ethics approval and consent to participate

All the animal experiments were approved by the Animal Research Ethics Board, University of Saskatchewan.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

The primer sequences and the annealing temperature used to amplify CP toxin genes including the expected size of the amplicons.

Additional file 2: Figure S1.

Differential gene expression in clinical and subclinical conditions. A Volcano plots showcase the balance between gene expression magnitude and statistical significance. B MA plots illustrate the correlation of average gene expression with log2 fold changes. C Box plots highlight the top 5 classifications based on their significant LogFC and p-values.

Additional file 3: Figure S2

. Detailed analysis of differential gene expression. A Volcano plots provide a holistic representation, emphasizing classifications with pronounced differential expressions. B MA plots detail the data, showing average expression versus log2 fold changes. C Box plots display the top 5 classifications chosen for their pronounced LogFC and compelling p-values.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gautam, H., Ayalew, L.E., Shaik, N.A. et al. Exploring the predictive power of jejunal microbiome composition in clinical and subclinical necrotic enteritis caused by Clostridium perfringens: insights from a broiler chicken model. J Transl Med 22, 80 (2024). https://doi.org/10.1186/s12967-023-04728-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04728-w

Keywords