Next Article in Journal
Transcriptome Analysis Reveals the Genes Involved in Bifidobacterium Longum FGSZY16M3 Biofilm Formation
Previous Article in Journal
Explorative Meta-Analysis of 417 Extant Archaeal Genomes to Predict Their Contribution to the Total Microbiome Functionality
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 2
10.3390/microorganisms9020382
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance

by
Yuxuan Hong
1,2,
Yang Cheng
1,2,
Leluo Guan
3,
Zutao Zhou
1,2,
Xiaowen Li
1,2,
Deshi Shi
1,2 and
Yuncai Xiao
1,2,*
1
College of Veterinary Medicine, Huazhong Agricultural University, Wuhan 430070, China
2
State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan 430070, China
3
Department of Agricultural, Food & Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(2), 382; https://doi.org/10.3390/microorganisms9020382
Submission received: 6 January 2021 / Revised: 4 February 2021 / Accepted: 10 February 2021 / Published: 13 February 2021
(This article belongs to the Section Gut Microbiota)
Browse Figures
Review Reports Versions Notes

Abstract

:
Bacillus amyloliquefaciens TL promotes broiler chicken performance by improving nutrient absorption and utilization and reducing intestinal inflammation. In this study, RNA-sequencing (RNA-seq)-based transcriptomes of ileal tissues collected from probiotic-fed and control broiler chickens were analyzed to elucidate the effects of the probiotic B. amyloliquefaciens TL, as a feed additive, on the gut immune function. In total, 475 genes were significantly differentially expressed between the ileum of probiotic-fed and control birds. The expression of genes encoding pyruvate kinase, prothymosin-α, and heat stress proteins was high in the ileum of probiotic-fed birds (FPKM > 500), but not in the control group. The gene ontology functional enrichment and pathway enrichment analyses revealed that the uniquely expressed genes in the control group were mostly involved in immune responses, whereas those in the probiotic group were involved in fibroblast growth factor receptor signaling pathways and positive regulation of cell proliferation. Bacillus amyloliquefaciens TL downregulated the expression of certain proinflammatory factors and affected the cytokine–cytokine receptor interaction pathway. Furthermore, B. amyloliquefaciens TL in broiler diets altered the expression of genes involved in immune functions in the ileum. Thus, it might contribute to improved broiler growth by regulating the immune system and reducing intestinal damage in broilers.

1. Introduction

Over the past 20 years, the broiler chicken industry has made substantial advances to improve bird productivity through breeding and nutrition management. Various feed additives such as antimicrobials, probiotics, and prebiotics have been widely applied to improve the growth performance of birds [1]. Among them, probiotics have been shown to improve immunity, gut function, health, and growth parameters in broiler chickens [2,3]. Moreover, supplementing the diet with probiotics has been reported to enhance nutrient digestibility and improve the cecal microflora composition in broiler chickens [4]. It has been demonstrated that dietary supplementation of Bacillus spp. can alter the intestinal microbiota to help prevent infectious diseases and increase productivity [5,6,7,8,9]. It has also been shown that some Bacillus subtilis strains can promote growth in chickens [10]. However, the effects of probiotics on gastrointestinal tract physiology and function in chickens are still poorly understood.
In birds, the intestine, especially the small intestine, is the main site of absorption and metabolism of nutrients derived from feed; it is also a primary site for pathogen entry [11,12]. It has important immunological, endocrine, and regulatory functions [13] that substantially affect animal health [12,14]. The ileum, the longest segment of the small intestine, harbors the majority of gut-associated lymphoid tissue [15], has a higher abundance of microbiota and stronger fermenting ability than the duodenum and jejunum, and it produces various microbial metabolites [16,17]. The gut microbiome can affect host animal health in several aspects, including the provision of nutrients and vitamins, protection against pathogens, development of the immune system, and homeostasis of epithelial mucosa [18]. Therefore, the establishment and maintenance of beneficial interactions between the host and its associated gut microbiota is a key health requirement [19]. To date, most studies have mainly focused on the cecum and reported that probiotics can inhibit the negative effects of intestinal disorders in broiler chickens by reducing cecal pathogens and regulating stress reactions [20,21]. However, there is still a gap in our understanding regarding the ileum and its functional changes in response to probiotics in chicken.
Studies have revealed the effects of immunity on broiler growth performance [22,23,24]. Klasing and Austic showed that immune system stimulation leads to a decrease in the growth performance of growing chicks [25]. In general, immune system activation is energetically costly, and long-term stimulation can have negative effects on the host [24]. Our previous study showed that dietary supplementation with Bacillus amyloliquefaciens TL improves the growth of broilers [26]. We speculated that this phenomenon was due to B. amyloliquefaciens TL, which can alter the gut microbiota and subsequently lower intestinal inflammation and the systemic immune response, diverting more energy toward growth. Therefore, in this study, a transcriptome analysis of the ileum of broiler chickens was performed to elucidate the functional changes in this tissue in response to the probiotic B. amyloliquefaciens TL. This study will provide a theoretical basis for future studies on the regulatory and molecular mechanisms underlying the effects of probiotics on animal performance.

2. Materials and Methods

2.1. Laboratory Animals

The present study was performed in strict accordance with the Guide for the Care and Use of Laboratory Animals Monitoring Committee of Hubei Province, China, and the protocol was approved by the Committee on the Ethics of Animal Experiments of the College of Veterinary Medicine, Huazhong Agricultural University (NO. HZAUCH-2017-012). For this study, chickens were reared in cages. In total, 120 1-day-old male Cobb broiler chickens were randomly divided into two groups—control and probiotic—with each group consisting of four pens and 15 chickens per pen (n = 60 chickens per group in a total of 120). The control group chickens were fed a basal diet (no drugs or additives), and the probiotic group chickens were fed grain with B. amyloliquefaciens TL 4 × 106 cfu/g) from 1 to 35 days of age. The B. amyloliquefaciens TL strain was obtained from Hubei Huada Real Technology Co. (Wuhan, China). B. amyloliquefaciens TL were added to the feed in the form of dry powder at a concentration of 200 g/ton, with an effective viable number of 2.0 × 1010 cfu/g; the dose was advised by the product description. Previous studies in this laboratory have shown that the effect of probiotics at a concentration of 4×106 cfu/g is similar to the feed additive chlortetracycline at 50 g/ton [26].
The chicken coop (which housed all pens) and the surrounding areas were disinfected with potassium permanganate and formalin before the trial. The temperature of the chicken coop was maintained at approximately 33 °C until the chickens were 7 days old; then, it was gradually reduced to 23 °C when the chickens were between 7 and 21 days of age and maintained at 23 °C thereafter. To keep the chickens healthy, the coop was cleared of manure daily, and the chickens had access to free ventilation; water and feed were provided ad libitum. Feed consumption by chickens in both groups and residual feed were recorded daily. Body weights of the chickens were measured on days 1, 7, 14, 21, and 35 (n = 60 chickens per group in a total of 120). During the experiment, the control group and the probiotic group each had two chickens die over 1–7 days due to individual weakness. Thereafter, the weights of the dead broilers and the remaining feed weight were recorded to correct the growth performance data.

2.2. Tissue Sampling

Six chickens were sampled from the control and probiotic groups on day 21 (n = 6 samples per group in a total of 12 samples), and the pen number and weight of each individual before slaughter were recorded. The method of selecting samples is: select the individuals in the four pens that are closest to the average weight of each pen. In order to expand the biological replicate sample, we selected the two closest to the average weight of each group from all individuals in each group. We selected the most representative individuals for sampling as much as possible to sequence the gut transcriptome and used these representative data to predict the entire population. The basic information of the sampled individuals is shown in Supplementary Table S1. Following their respective fasting periods, the birds were euthanized by cervical dislocation, and the intestinal tissues were collected for further analyses. All tissues were washed with ice-cold phosphate buffer solution (Sigma-Aldrich, St. Louis, MO, USA), snap-frozen in liquid nitrogen, and stored at −80 °C until mRNA expression analysis.

2.3. RNA Preparation and Sequencing

The total RNA was extracted from each ileal sample using the RNATM.iso PLUS Kit (Takara Biotechnology, Shiga, Japan). RNA degradation and purity were assessed by 1% agarose gel electrophoresis, and purity was further examined using a NanoPhotometer® spectrophotometer (Implen, Westlake Village, CA, USA) [27]. RNA was quantified using the Qubit® RNA Assay Kit and Qubit® 2.0 Fluorometer (Life Technologies, San Diego, CA, USA) [27]. RNA integrity was assessed using the RNA Nano 6000 Assay Kit with a Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA) [27]. The samples that passed quality control standards (specifically, the target band of RNA electrophoresis gel is bright and clear, there is no diffusion area in the swimming lane, and no protein and DNA contamination; RNA Integrity Number (RIN) value is close to 10; 28S/18S is greater than or equal to 1.5; 1.7 < OD260/280 < 2.0; OD260/230 ratio is 2.5) were submitted to Novogene Co., Ltd. (Beijing, China) for library preparation and sequencing using the Illumina HiSeq platform, and 125–/150 bp paired-end reads were generated (n = 6 samples per group in a total of 12 samples) [28].

2.4. Quality Control

Transcriptome assembly and annotation protocols were provided by Novogene Co., Ltd. (Beijing, China). To obtain high-quality sequences, the raw reads were filtered by removing the following: (1) reads containing adapters; (2) reads containing poly-N; (3) low-quality reads [29]. Once the clean data were obtained, Q20, Q30, and GC content were calculated. The subsequent analyses were performed using the high-quality clean data of [27].

2.5. Mapping Reads to the Reference Genome

Reference genome and gene model annotation files were downloaded directly from the genome website (https://www.ncbi.nlm.nih.gov (accessed on 13 February 2021) [10/2017]). An index was built for the reference genome using Bowtie v2.2.3, and paired-end clean reads were aligned to the reference genome using TopHat v2.0.12 [27]. We selected TopHat as the mapping tool because it can generate a database of splice junctions based on the gene model annotation file, yielding the best mapping results compared to all nonsplice mapping tools [30].

2.6. Quantification of Gene Expression

HTSeq v0.6.1 was used to count the read numbers mapped to each gene. The fragments per kilobase of transcript per million fragments mapped (FPKM) of each gene were then calculated based on the length of the gene and the number of reads mapped to it. The expected FPKM number simultaneously considers the effect of sequencing depth and gene length on the number of reads and is currently the most commonly used method for estimating gene expression levels [31].

2.7. Differential Expression Analysis

Differential expression analysis was performed based on two conditions/groups (two biological replicates per condition) using the DESeq R package (1.18.0). DESeq provides statistical routines for detecting differential expression in the digital gene expression data using a model that is based on a negative binomial distribution. The resulting p-values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate [32]. Genes with adjusted p-values < 0.05 were detected by DESeq and were assigned as differentially expressed fragments (i.e., per kilobase of transcript sequence per million base pairs sequenced).

2.8. Gene Ontology (GO) and Pathway Enrichment

Based on the DEG data after screening, pathway enrichment with the Kyoto Encyclopedia of Genes and Genomes (KEGG) and GO was performed using DAVID (https://david.ncifcrf.gov/ (accessed on 13 February 2021) [10/2017]), with Gallus gallus as the reference.

2.9. Quantitative Real-Time Polymerase Chain Reaction Analysis

For quantitative real-time polymerase chain reaction (qRT-PCR) validation experiments, six genes were randomly selected to assess the RNA-sequencing (RNA-Seq) data (Table 1). For this analysis, 1 μg of RNA was reverse-transcribed into cDNA using the PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Dalian, China) according to the manufacturer’s instructions. cDNA was diluted 10-fold and used for real-time PCR analyses with a Bio-Rad CFX96TM System and signal detection protocols in accordance with the manufacturer’s instructions (Takara). qRT-PCR experiments were performed in three technical replicates. β-Actin was used as the endogenous control. The expression of individual genes was normalized to that of β-Actin, a house-keeping gene [33]. Primers used for qRT-PCR are shown in Table 1. Data were analyzed using GraphPad Prism v 5.0 Software (GraphPad, Inc., La Jolla, CA, USA).

2.10. Statistical Analyses and Data Records

Statistical evaluations were performed using SPSS for Windows, version 22 (SPSS, Inc., Chicago, IL, USA). Graphs were plotted using GraphPad Prism 5 software (GraphPad, Inc.). The results are presented as the mean ± SEM. The significance levels for all analyses were set as p < 0.05 (*) and p < 0.01 (**).

3. Results

3.1. Growth Performance

The average daily gain (ADG) and overall feed conversion ratio (FCR) are shown in Table 2. The average daily gain of the probiotic group was significantly higher than that of the control group at 0–7 (p < 0.05) and 14–21 days of age, and the difference was the most significant at 14–21 days (p < 0.01). At 14–21 days of age, the average weight and average daily gain of broilers in the probiotic group were 25.33% and 51.47% higher than those in the control group, respectively. Additionally, the overall FCR of broilers in the probiotic group was lower than that in the control group (Table 2).

3.2. Transcriptome Analysis

The transcriptome data from the 12 ileum samples are presented in Supplementary Table S2 (n = 6 samples per group in a total of 12 samples). After removing the low-quality data, 6.83–10.38 Gb of clean data were obtained for further analyses, with >88% of the filtered reads having a Phred quality score of >30 (base Q30). As shown in Figure 1a, the Pearson correlation (R2) of the biological repeats in the control and probiotic groups was high, indicating that the sequencing data were suitable for subsequent analyses. The transcriptome profiles of the ileal tissues from the probiotic group were separated from those of the control group in the Principal Component Analysis (PCA) plot (Figure 1b). The expression of 14,033 and 14,446 genes was detected in the ileal tissues of birds from the control and probiotic groups, respectively, with 13,759 genes expressed in both groups (defined as the core transcriptome) (Figure 1c), based on the FPKM, which was >1 and was detected in >50% of the samples in each group.

3.3. Functional Analysis of Broiler Chicken Transcriptomes

We performed a functional enrichment analysis of all genes expressed in the ileal tissues of broilers in the two groups (n = 6 samples per group in a total of 12 samples). The top 60 functionally enriched GO terms (the top 20 terms each in the biological process (BP), cellular component (CC), and molecular function (MF) categories) in the ileal tissue of the control and probiotic groups are shown in Figure 2a,b, respectively. Among them, 54 GO terms were enriched in both groups, six were enriched in the control group only, and six were enriched only in the probiotic group (Figure 2c and Supplementary Table S3). In addition to normal cellular activities such as “cell division”, “DNA replication”, and “intracellular protein transport”, the GO terms related to intestinal function that were enriched in both groups included “positive regulation of I-kappaB kinase/NF-kappaB signaling”, “cellular response to mechanical stimulus”, “positive regulation of defense response to virus by host”, and “xenophagy”.
We performed functional enrichment analysis of the uniquely expressed genes in the ileal tissues of the two groups of broilers (n = 6 samples per group in a total of 12 samples). In total, 274 genes were only expressed in the control group and 687 were expressed only in the probiotic group. The following GO terms were enriched for gene function in the control group: “immune response”, “chemotaxis”, “regulation of inflammatory response”, “thrombin receptor signaling pathway”, and “somatic hypermutation of immunoglobulin genes” (Figure 3a). The following GO terms were enriched in the uniquely expressed genes of the probiotic group: “fibroblast growth factor receptor signaling pathway”, “positive regulation of cell proliferation”, “integral component of plasma membrane”, and “G-protein coupled serotonin receptor activity” (Figure 3b). The significantly enriched KEGG pathways involving the uniquely expressed genes in the control group included “intestinal immune network for IgA production” and “cytokine–cytokine receptor interaction” (Figure 3c). The uniquely expressed genes in the probiotic group were mainly enriched in “regulation of actin cytoskeleton” and “neuroactive ligand–receptor interaction” (Figure 3d). There were 243 and 229 highly expressed genes (FPKM > 500) in the control and probiotic groups, respectively (Supplementary Table S4a,b). Genes encoding pyruvate kinase M (PKM), prothymosin-α, and heat stress proteins (HSPs) were highly expressed in the ileum of the probiotic group but not in that the control group (Supplementary Table S4c).

3.4. Differentially Expressed Genes in the Ileal Tissue and Their Functions

In total, 475 genes were found to be differentially expressed in the ileum (p < 0.05) (n = 6 samples per group in a total of 12 samples) with 321 upregulated (67.58%) and 154 downregulated genes (32.42%) in the probiotic group based on the cut-off for differentially expressed genes (DEGs) with p-adj < 0.05 and detected in >50% of samples within each group (Supplementary Table S5). The GO enrichment analysis showed that the functions of DEGs could be classified as BP, CC, and MF categories. The enriched GO terms of the DEGs are shown in Figure 4a; a high proportion of the DEGs upregulated in the probiotic group were involved in BPs including “regulation of cell migration”, “multicellular organism development”, “neural crest cell migration”, and “coronary vasculature development” categories (Figure 4a and Supplementary Table S6). The enriched GO terms of the downregulated DEGs in the probiotic group were in the “establishment of T cell polarity” and “immunoglobulin production” categories (Figure 4b and Supplementary Table S6). In addition, KEGG pathways were used to predict gene functions. The pathway analysis of downregulated DEGs between the two groups highlighted “cytokine–cytokine receptor interaction” as the most significantly influenced pathway. The upregulated DEGs were significantly enriched in the “calcium signaling pathway”, “focal adhesion”, and “insulin signaling pathway” (Table 3 and Figure 4c,d).

3.5. Confirmation of Gene Expression Data by Quantitative Reverse Transcription Polymerase Chain Reaction

To validate the identified ileal DEGs between the control and probiotic groups determined by RNA-Seq, six DEGs (CXCR5, CXCL13, IL-22, ACACB, LEPR, and LAMA5) were selected for qRT-PCR analysis. The results of high-throughput RNA-seq revealed that ACACB, LEPR, and LAMA5 were upregulated, whereas CXCR5, CXCL13, and IL-22 were downregulated. The qRT-PCR results showed similar trends of upregulation and downregulation (Figure 5), validating the high-throughput RNA-seq results.

4. Discussion

Our study showed that poultry feed supplemented with B. amyloliquefaciens TL significantly increased the weight of broilers at 21 days (p < 0.05) and reduced the feed:gain ratio. The results indicate that B. amyloliquefaciens TL as a probiotic promotes growth and improves feed efficiency. A previous study has shown that although dietary supplementation with the probiotic B. subtilis does not have any effect on broiler growth performance between 1 to 21 days, broilers fed probiotics presented higher ADG and average daily feed intake and a lower FCR than broilers with no probiotic supplementation (p < 0.05) at 22–42 days [34]. Furthermore, Willie et al. reported that broiler chicks supplemented with commercial probiotic cultures (Lactobacillus acidophilus, Lactobacillus casei, Bifidobacterium, and Enterococcus faecalis) showed the highest weight gain of 0.62 kg at 3 weeks [35]. Furthermore, the addition of probiotics (L. acidophilus, B. subtilis, and Clostridium butyricum) to broiler diets resulted in the highest weight gain at 22–35 days, a 7.91% increase over that of the control group (p < 0.05) [36]. Although performance improvement varied in these trials, it is evident that probiotic supplementation can lead to higher productivity of broilers at different growth stages. However, these studies did not use B. amyloliquefaciens TL, and its growth-promoting effect in broilers was particularly evident at 14–21 days of age, which is different from the results of previous studies. The present study results suggest that the growth-promoting effect of B. amyloliquefaciens TL is better during the broiler brooding period (1–21 days) than the fattening period (22–35 days).
This is the first study to report the transcriptome and its functions in the ileum of broiler chickens at 21 days of age and identify whether B. amyloliquefaciens TL could affect tissue function at the molecular level. The transcriptome profiles of the ileal mucosal genes in two fast-growing chicken hybrids (HA and HB) at 43 days have been reported [12]. The body weight, daily weight gain, and daily feed intake of HB broilers were higher than those of the HA broilers, but HA broilers showed earlier development characteristics. The ileal mucosal transcriptome of Zampiga et al. showed that a high percentage of gene sets involved in cell energy metabolism, mitochondria structure and functionality, ribosome structure, protein synthesis, cell structure and integrity, and antioxidant and detox mechanisms could be observed in the HA group. Gene sets involved in immune system activation, signal transduction and cell signal transduction, DNA remodeling and replication-chromatin/histone modification, cell activation, migration and adhesion, inflammation, and bone remodeling were detected in the HB group [12]. Here, the functional analysis of the core transcriptome revealed that in addition to some normal cellular activities such as cell division and protein transport, which were similar to the results of Zampiga et al., the following functions were enriched: “positive regulation of I-kappaB kinase/NF-kappaB signaling”, “cellular response to mechanical stimulus”, “positive regulation of defense response to virus by host”, and “xenophagy”. This indicates that the ileum is involved in positively regulating the activation of the immune system to resist bacterial and viral invasions in the intestine. Xenophagy is a selective form of autophagy and is also specifically used to inhibit cellular activities of intracellular bacteria [37]. These results concur with the transcription profile of HB broilers [12]. The ileum transcriptome is involved in the regulation of some immune signaling pathways and defense responses to viruses and bacteria, suggesting that these immune functions are among the core functions of the chicken ileum.
In the present study, the uniquely expressed genes in the ileal tissue of the control group were mostly related to immune responses (e.g., CD28, CXCR5, VPREB3, and CCR6), with "immune response" and "chemotaxis" GO terms and "intestinal immune network for IgA production" and "cytokine–cytokine receptor interaction" pathways. CD28 is a secondary signaling receptor that activates T cells by interacting with costimulatory molecules on antigen-presenting cells, such as CD80/86 [38]. In chickens, CD28 is expressed in most αβ T cells, and it has similar functional properties to mammalian CD28 [39]. CCR6 is a specific receptor for CCL20, and both have been cloned and expressed in chickens. Chicken CCR6 is mainly expressed in bone marrow, secondary lymphoid organs, and on the mucosal surface [40]. The expression of chicken CCR6 suggests immature dendritic cells [41]. Chicken VpreB3 is expressed in immature and/or mature IgM-positive cells, binds to free IgLC, and negatively regulates its maturation and secretion [42]. Studies have shown that after the injection of lipopolysaccharides, the CXCR5 level in the bursa decreases, and the CXCR5 mRNA level in the cecum tonsil increases. Chicken CXCR5 might be involved in the migration of bursal immune cells [43]. Here, these genes were expressed only in the control group, indicating that the immune activity of the intestines of 21-day-old broilers is more active in the control than in the probiotic group. Studies have shown that probiotics can effectively reduce the risk of intestinal infection and colonization of multiple pathogens [44,45]. We speculate that in the ileum of birds receiving B. amyloliquefaciens, there could be fewer pathogens. However, future studies are needed to quantify the pathogen populations between the two groups to support such speculations.
Here, FGF19, FGF16, FGF10, TRIM71, and AREG were only expressed in the probiotic group, and they are mainly associated with the GO terms “fibroblast growth factor receptor signaling pathway” and “positive regulation of cell proliferation.” Fibroblast growth factors (FGFs) are signaling proteins that were originally discovered as stimuli for the growth of fibroblasts or epithelial cells and have multiple functions in development, metabolism, and nerve function [46]. Chicken FGF10 is expressed predominantly in abdominal fat and is involved in preadipocyte differentiation and adipocyte maturation [47]. FGF19 is mainly secreted by ileal enterocytes in response to bile acid activation of the nuclear receptor FXR, and then, FGF19 acts on the cell surface receptor complex in hepatocytes to inhibit bile acid synthesis and gluconeogenesis and stimulate glycogen protein synthesis [48]. Studies on rodents claim that FGF16 also has a regulatory effect on bile acid metabolism [49]. However, in birds, the role of FGF19 in retinal development and that of FGF16 in inner ear development have been reported [50,51], whereas research on the functions of FGF19 and FGF16 in the gut is relatively limited. Considering these findings, the unique expression of the three genes of the FGF family in the probiotic group suggests that glycogen and protein syntheses and fat differentiation were more active in these broilers. AREG is a ligand for epidermal growth factor receptor (EGFR) and is involved in a wide range of physiological processes, including breast development, blastocyst implantation, keratinocyte proliferation, nerve and bone tissue development, and axon growth [52,53]. TRIM71 belongs to the TRIM-NHL protein family. The TRIM proteins control important cellular processes, such as intracellular signaling, innate immunity, transcription, autophagy, carcinogenesis, and antiretroviral action. For example, TRIM22 can inhibit human immunodeficiency virus 1 long-term and repeat-driven transcription according to the RING domain, and TRIM62 plays an important role in restricting the replication of avian leukemia virus subgroup J. In this study, TRIM71 was uniquely expressed in the probiotic group. TRIM71 plays a conservative role in regulating early development and differentiation and has the function of inhibiting tumorigenesis [54], but there is no related function in birds. The unique expression of these genes in the probiotic group indicates that this group has higher disease resistance, which would enable the rapid growth of broilers at this stage.
The ileal transcriptome showed that the PKM, prothymosin-α, and HSPA8 genes were highly expressed in the probiotic group, but not in the control group. PKM activation might increase glucose metabolism flux by catalyzing the conversion of phosphoenolpyruvate to ATP and pyruvate, which is used as a raw material for synthesizing triglycerides from acetyl-CoA [55]. Therefore, the higher PKM expression in the ileal tissues of the probiotic group indicates that the broilers in this group have higher metabolism and energy requirements. Prothymosin-α is an immunomodulator with antiviral, antifungal, and anticancer functions, with several features similar to IL-1α [56]. The high expression of prothymosin-α in the probiotic group might lead to improved immune function and disease resistance, which could be one of the reasons for the better growth of B. amyloliquefaciens TL-fed chickens. HSPA8 encodes a member of the heat shock protein 70 family, which plays a major role in maintaining a stable cellular environment under high heat stress conditions in animals [57,58]. One study showed that Hsp70 inhibition enhances heat stress injury and apoptosis in chicken primary cardiomyocytes, and Hsp70 has a cytoprotective effect during heat stress in chickens [59]. In addition, the HSP70 family has an immunomodulatory effect, which can suppress the immune phenotype by abrogating immune cell interactions, thereby inhibiting the inflammatory pathway [58]. The high expression of HSP70 in the probiotic group also indicates that feeding broilers the probiotic B. amyloliquefaciens TL enables them to respond more effectively to heat stress by regulating immunity. We did not perform a heat stress challenge in this study, and this conclusion will be verified by subsequent experiments.
Among the DEGs, the immune-related DEGs (such as IL-8, CXCL13L2, CXCR5, TNFRSF13C, and IFNGR1) showed significantly different expressions. IL-8, a downregulated DEG in the probiotic group, is an important immunomodulatory factor. IL-8 was the first CXC chemokine identified in chickens [60], and it plays an important role in the mucosa by attracting chemokines and neutrophils to cause an inflammatory response [61]. IL-8 mRNA was shown to increase in the intestinal tissues and livers of birds after infection with Salmonella [62]. According to Ateya et al., supplementing probiotics, acidifiers, and symbiotics to the diet of broilers infected with Escherichia coli can reduce the expression of ileal IL-8, which regulates intestinal inflammation in response to E. coli infection and minimizes inflammation-induced damage, thereby improving growth performance [63]. Here, the downregulation of IL-8 expression in the ileum of the probiotic group indicates that the supplementation with TL reduces chronic inflammation of the intestine. It might also reduce energy consumption caused by intestinal inflammation pressure.
CXCL13 and its receptor CXCR5 (which was downregulated in the probiotic group) also play important roles in infection, inflammation, and the immune response [64]. Studies have shown that chicken CXCR5 is expressed on the surface of all B cell and T cell subpopulations, and CXCL13 effectively binds to cells expressing CXCR5. CXCL13 is a B cell chemoattractant and plays an important role in lymphoid neoplasia. It is also extensively involved in the pathogenesis of several autoimmune diseases, inflammatory conditions, and lymphoproliferative disorders [65]. The chemokine-receptor pair (CXCL13-CXCR5) is also conserved between mammalian and avian species [66]. Considering the important role of this complex in inflammatory diseases, its downregulation indicates a lower level of ileal inflammation in broiler chickens in the probiotic group.
TNFRSF13C encodes the BAFF receptor (BAFFR). Chicken BAFF is synthesized by bursal stromal cells or immature B cells [67]. It is expressed mainly in the bursa of Fabricius. The receptor has also been detected in peripheral B lymphocytes in organs such as the cecal tonsil, and it can bind to BAFF [68]. BAFF treatment increases the number of B cells in the cecal tonsil and the expression of BAFFR [68]. Chicken BAFF and BAFFR play important roles in early B cell survival and mature B cell proliferation in chickens [68]. Downregulation of TNFRSF13C indicates that B. amyloliquefaciens TL might reduce peripheral B lymphocyte numbers in the ileum and intestine-specific immune response.
Avian IFNGR1 is one of the receptors of IFN-γ, which is mainly expressed in chicken lymphoid tissues and organs; it is also expressed in tissue cells such as the small intestine and breast muscle cells [69]. In birds, IFN-γ is a key cytokine for I-type T helper (Th1) cell responses and is essential for controlling infection by intracellular pathogens [70]. Avian IFN-γ can be produced by macrophages, dendritic cells, and CD4+ and CD8+ T cells, as in mammals [71]. It has been reported that IFN-γ in mammals can prevent the infiltration and destruction of virions in cells during virus infection, as well as the formation of virions or germination of viruses. It can further stimulate cells of the immune system to protect the human body [72]. Kidane et al. reported that in vitro attenuated Histomonas meleagridis vaccination results in high levels of IFN-γ in the cecum of turkeys, which might involve important mechanisms of immune responses against avian diseases [71]. In the present study, the decrease in IFNGR1 expression can be attributed to the lower stimulation of the ileal immune system by viruses and pathogens and a lower Th1 immune response in the ileum of broilers in the probiotic group. Similar studies have shown that a decrease in IFN-γ production by spleen cells in Lactobacillus-fed chickens reflects a selective decrease in Th1 cell activation [73].
CCL19 could act as a lymphocyte chemoattractant factor via CCR7 in some virus infections or inflammatory disorders [74]. For example, studies have shown that CCL19 played an important role in T cell migration into bursae during IBDV infection [75,76]. In short, the role of CCL19 in chickens is summarized as follows: on the one hand, the main role of CCL19 is to recruit immune cells into the lesions to limit virus infection and inflammatory responses; on the other hand, CCL19 could facilitate virus invasion [76]. The downregulation of the CCL19 gene in the ileum of broiler chickens in the probiotic treatment group showed that TL reduced the chemotaxis of T cells, thereby reducing the energy consumption of the immune response.
The GO analysis of DEGs revealed that the downregulated genes were predominantly involved in the “establishment of T cell polarity” and “immunoglobulin production”. The KEGG analysis of the downregulated genes indicated that they were mainly involved in the “cytokine–cytokine receptor interaction”, and most DEGs enriched in this pathway are related to the immune response. Figure 6 shows the proposed regulatory and molecular mechanisms of B. amyloliquefaciens TL. We conclude that the effect of B. amyloliquefaciens TL mainly involves the downregulation of the expression of some immune-related genes in the ileum. In the probiotic group, the expression of some important inflammatory factors and their receptors at 21 days were lower than those in the control group, thereby reducing the activity of immune cells and the production of immunoglobulins. This indicates the following benefits of supplementing probiotics in the feed: a reduction in intestinal inflammation levels and intestinal damage.
In recent years, increasing attention has been paid to the use of probiotics as an alternative feed additive to antibiotics. Previous studies have confirmed the ability of probiotics to enhance the gut microbial balance and consequently the natural defense system of animals against pathogenic bacteria. During production, broilers are susceptible to environmental and pathogenic stresses, which are detrimental to growth. Probiotics act as potential immune modulators and protect these birds during growth. Gene expression analysis in the cecal tonsils of chickens treated with probiotics and challenged with Salmonella indicated that probiotics inhibit IL-12 and IFN-γ secretion [77]. In the present study, the reduced immune activity of the probiotic group might have assisted broilers in diverting nutritional energy to promote body growth. Excessive immune stress in broilers can affect growth performance. Therefore, the industry should avoid unnecessary innate immune activity, such as chronic inflammation. The probiotic B. amyloliquefaciens TL can reduce immune stress to some extent, which explains the weight gain in our probiotic group. The present study provides evidence that B. amyloliquefaciens TL positively contributes to the health of the host, indicating its probiotic potential. Future research must focus on the mechanism of the interaction between the intestinal flora and specific host intestinal epithelial cells. The effectiveness of probiotics depends on the bacterial strain and host health status. Therefore, it is important to determine the effect of probiotics in different hosts in different health states.

5. Conclusions

In conclusion, the results of this study indicate that B. amyloliquefaciens TL alters ileal gene expression. This probiotic reduced the expression of genes involved in the inflammatory response, intestinal inflammatory factors, and receptors in the ileal tissue of 21-day-old broilers. This suggests that probiotics might reduce the damage to epithelial cells from inflammation, and thereby reduce the burden of pathogens and energy consumption caused by immune stimulation and activity, concomitantly improving broiler growth. Our study provides new insights into the role of probiotics in the diet. The findings might help optimize the use of probiotics as feed additives, instead of antibiotics, to achieve high productivity in the broiler industry.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/9/2/382/s1, Table S1: Information of sampled individuals; Table S2: Summary of sequencing data for each sample; Table S3a: GO terms enriched by group C and P; Table S3b: GO terms uniquely enriched by group C; Table S3c: GO terms uniquely enriched by group P; Table S4a: Highly abundant genes in group C; Table S4b: Highly abundant genes in group P; Table S4c: Genes highly expressed in group P but not abundantly expressed in group C; Table S5a: ALL DEGs between group C and group P; Table S5b: UP DEGs between group C and group P; Table S5c: DOWN DEGs between group C and group P; Table S6: GO terms enriched with DEGs.

Author Contributions

Conceptualization, Y.X. and L.G.; methodology, Y.X., L.G., and Y.C.; resources, Z.Z.; data curation, Y.H.; writing—original draft preparation, Y.H. and Y.C.; writing—review and editing, L.G., Y.H., and Y.X.; supervision, Z.Z., D.S. and X.L.; project administration, Y.X.; funding acquisition, Y.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2017YFD05010001, and the Innovative Job Funds of Agricultural Science and Technology of Hubei Province, grant number 2019-620-000-001-30.

Institutional Review Board Statement

The experimental animal procedures were approved by the Committee on the Ethics of Animal Experiments of the College of Veterinary Medicine, Huazhong Agricultural University, Wuhan, China, under permit No. HZAUCH-2017-012 (2017).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

We are grateful to Shuaifei Feng (Huazhong Agricultural University) for help in genome sequence analysis. We thank Novogene Co., Ltd. (Beijing, China; https://www.novogene.com/ (accessed on 13 February 2021) [02/2021]) for the technical support in genome sequencing, assembly, and annotation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jin, L.; Ho, Y.; Abdullah, N.; Jalaludin, S. Probiotics in poultry: Modes of action. World Poult. Sci. J. 1997, 53, 351–368. [Google Scholar] [CrossRef]
  2. Al-Khalaifa, H.; Al-Nasser, A.; Al-Surayee, T.; Al-Kandari, S.; Al-Enzi, N.; Al-Sharrah, T.; Ragheb, G.; Al-Qalaf, S.; Mohammed, A. Effect of dietary probiotics and prebiotics on the performance of broiler chickens. Poult. Sci. 2019, 98, 4465–4479. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, Y.; Wang, B.; Zeng, Z.; Liu, R.; Tang, L.; Gong, L.; Li, W. Effects of probiotics Lactobacillus plantarum 16 and Paenibacillus polymyxa 10 on intestinal barrier function, antioxidative capacity, apoptosis, immune response, and biochemical parameters in broilers. Poult. Sci. 2019, 98, 5028–5039. [Google Scholar] [CrossRef] [PubMed]
  4. Majidi-Mosleh, A.; Sadeghi, A.; Mousavi, S.; Chamani, M.; Zarei, A. Ileal MUC2 gene expression and microbial population, but not growth performance and immune response, are influenced by in ovo injection of probiotics in broiler chickens. Br. Poult. Sci. 2017, 58, 40–45. [Google Scholar] [CrossRef] [PubMed]
  5. Cheng, Y.H.; Zhang, N.; Han, J.C.; Chang, C.W.; Hsiao, F.S.H.; Yu, Y.H. Optimization of surfactin production from Bacillus subtilis in fermentation and its effects on Clostridium perfringens-induced necrotic enteritis and growth performance in broilers. J. Anim. Physiol. Anim. Nutr. 2018, 102, 1232–1244. [Google Scholar] [CrossRef] [PubMed]
  6. Li, C.L.; Wang, J.; Zhang, H.J.; Wu, S.G.; Hui, Q.R.; Yang, C.B.; Fang, R.J.; Qi, G.H. Intestinal Morphologic and Microbiota Responses to Dietary Bacillus spp. in a Broiler Chicken Model. Front. Physiol. 2018, 9, 1968. [Google Scholar] [CrossRef]
  7. Grant, A.; Gay, C.G.; Lillehoj, H.S. Bacillus spp. as direct-fed microbial antibiotic alternatives to enhance growth, immunity, and gut health in poultry. Avian Pathol. 2018, 47, 339–351. [Google Scholar] [CrossRef] [PubMed]
  8. Whelan, R.A.; Doranalli, K.; Rinttilä, T.; Vienola, K.; Jurgens, G.; Apajalahti, J. The impact of Bacillus subtilis DSM 32315 on the pathology, performance, and intestinal microbiome of broiler chickens in a necrotic enteritis challenge. Poult. Sci. 2019, 98, 3450–3463. [Google Scholar] [CrossRef]
  9. Ramlucken, U.; Ramchuran, S.O.; Moonsamy, G.; Lalloo, R.; Thantsha, M.S.; Jansen van Rensburg, C. A novel Bacillus based multi-strain probiotic improves growth performance and intestinal properties of Clostridium perfringens challenged broilers. Poult. Sci. 2020, 99, 331–341. [Google Scholar] [CrossRef] [PubMed]
  10. Opalinski, M.; Maiorka, A.; Dahlke, F.; Cunha, F.; Vargas, F.; Cardozo, E. On the use of a probiotic (Bacillus subtilis-strain DSM 17299) as growth promoter in broiler diets. Braz. J. Poult. Sci. 2007, 9, 99–103. [Google Scholar] [CrossRef] [Green Version]
  11. Abreu, M.T. Toll-like receptor signalling in the intestinal epithelium: How bacterial recognition shapes intestinal function. Nat. Rev. Immunol. 2010, 10, 131. [Google Scholar] [CrossRef] [PubMed]
  12. Zampiga, M.; Bertocchi, M.; Bosi, P.; Trevisi, P.; Meluzzi, A.; Sirri, F. Differences in productive performance and intestinal transcriptomic profile in two modern fast-growing chicken hybrids. J. Anim. Physiol. Anim. Nutr. 2019, 103, 125–134. [Google Scholar] [CrossRef] [Green Version]
  13. Scanes, C.G.; Pierzchala-Koziec, K. Biology of the gastro-intestinal tract in poultry. Avian Biol. Res. 2014, 7, 193–222. [Google Scholar] [CrossRef]
  14. Sugiharto, S. Role of nutraceuticals in gut health and growth performance of poultry. J. Saudi Soc. Agric. Sci. 2016, 15, 99–111. [Google Scholar] [CrossRef] [Green Version]
  15. Barr, T.; Sureshchandra, S.; Ruegger, P.; Zhang, J.; Ma, W.; Borneman, J.; Grant, K.; Messaoudi, I. Concurrent gut transcriptome and microbiota profiling following chronic ethanol consumption in nonhuman primates. Gut Microbes 2018, 9, 338–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. O’Hara, A.M.; Shanahan, F. The gut flora as a forgotten organ. EMBO Rep. 2006, 7, 688–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhang, C.; Liu, Y.; Chen, S.; Qiao, Y.; Zheng, Y.; Xu, M.; Wang, Z.; Hou, J.; Wang, J.; Fan, H. Effects of Intranasal Pseudorabies Virus AH02LA Infection on Microbial Community and Immune Status in the Ileum and Colon of Piglets. Viruses 2019, 11, 518. [Google Scholar] [CrossRef] [Green Version]
  18. Lin, L.; Zhang, J. Role of intestinal microbiota and metabolites on gut homeostasis and human diseases. BMC Immunol. 2017, 18, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Sommer, F.; Bäckhed, F. The gut microbiota—Masters of host development and physiology. Nat. Rev. Microbiol. 2013, 11, 227. [Google Scholar] [CrossRef]
  20. Mohammed, A.A.; Jiang, S.; Jacobs, J.A.; Cheng, H.W. Effect of a synbiotic supplement on cecal microbial ecology, antioxidant status, and immune response of broiler chickens reared under heat stress. Poult. Sci. 2019, 98, 4408–4415. [Google Scholar] [CrossRef]
  21. Hu, Y.; Wang, L.; Shao, D.; Wang, Q.; Wu, Y.; Han, Y.; Shi, S. Selectived and Reshaped Early Dominant Microbial Community in the Cecum with Similar Proportions and Better Homogenization and Species Diversity due to Organic Acids as AGP Alternatives Mediate Their Effects on Broilers Growth. Front. Microbiol. 2019, 10, 2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Yang, X.; Guo, Y.; He, X.; Yuan, J.; Yang, Y.; Wang, Z. Growth performance and immune responses in chickens after challenge with lipopolysaccharide and modulation by dietary different oils. Animal 2008, 2, 216–223. [Google Scholar] [CrossRef] [Green Version]
  23. Lai, H.T.; Nieuwland, M.G.; Kemp, B.; Aarnink, A.J.; Parmentier, H.K. Effects of repeated intratracheally administered lipopolysaccharide on primary and secondary specific antibody responses and on body weight gain of broilers. Poult. Sci. 2011, 90, 337–351. [Google Scholar] [CrossRef]
  24. Yang, X.; Li, W.; Feng, Y.; Yao, J. Effects of immune stress on growth performance, immunity, and cecal microflora in chickens. Poult. Sci. 2011, 90, 2740–2746. [Google Scholar] [CrossRef]
  25. Klasing, K.; Austic, R. Changes in protein degradation in chickens due to an inflammatory challenge. Proc. Soc. Exp. Biol. Med. 1984, 176, 292–296. [Google Scholar] [CrossRef] [PubMed]
  26. Hong, Y.; Cheng, Y.; Li, Y.; Li, X.; Zhou, Z.; Shi, D.; Li, Z.; Xiao, Y. Preliminary Study on the Effect of Bacillus amyloliquefaciens TL on Cecal Bacterial Community Structure of Broiler Chickens. BioMed. Res. Int. 2019, 2019, 5431354. [Google Scholar] [CrossRef] [Green Version]
  27. Liu, W.; Qiu, X.; Song, C.; Sun, Y.; Meng, C.; Liao, Y.; Tan, L.; Ding, Z.; Liu, X.; Ding, C. Deep Sequencing-Based Transcriptome Profiling Reveals Avian Interferon-Stimulated Genes and Provides Comprehensive Insight into Newcastle Disease Virus-Induced Host Responses. Viruses 2018, 10, 162. [Google Scholar] [CrossRef] [Green Version]
  28. Li, X.; Jia, L.; Chen, X.; Dong, Y.; Ren, X.; Dong, Y.; Chen, Y.; Xie, L.; Liu, M.; Shiota, C.; et al. Islet α-cell Inflammation Induced by NF-κB inducing kinase (NIK) Leads to Hypoglycemia, Pancreatitis, Growth Retardation, and Postnatal Death in Mice. Theranostics 2018, 8, 5960–5971. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, Y.; Miao, G.; Wu, Q.; Lin, F.; You, C.; Wang, S.; Aweya, J.J.; Ma, H. Transcriptome sequencing and molecular markers discovery in the gonads of Portunus sanguinolentus. Sci. Data 2018, 5, 180131. [Google Scholar] [CrossRef] [PubMed]
  30. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  31. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511. [Google Scholar] [CrossRef] [Green Version]
  32. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  33. Wang, Y.; Gao, M.; Zhu, F.; Li, X.; Yang, Y.; Yan, Q.; Jia, L.; Xie, L.; Chen, Z. METTL3 is essential for postnatal development of brown adipose tissue and energy expenditure in mice. Nat. Commun. 2020, 11, 1648. [Google Scholar] [CrossRef] [Green Version]
  34. Bai, K.; Huang, Q.; Zhang, J.; He, J.; Wang, T. Supplemental effects of probiotic Bacillus subtilis fmbJ on growth performance, antioxidant capacity, and meat quality of broiler chickens. Poult. Sci. 2016, 96, pew246. [Google Scholar] [CrossRef]
  35. Willis, W.L.; Isikhuemhen, O.S.; Ibrahim, S.A. Performance assessment of broiler chickens given mushroom extract alone or in combination with probiotics. Poult. Sci. 2007, 86, 1856–1860. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, Z.F.; Kim, I.H. Effects of multistrain probiotics on growth performance, apparent ileal nutrient digestibility, blood characteristics, cecal microbial shedding, and excreta odor contents in broilers. Poult. Sci. 2014, 93, 364–370. [Google Scholar] [CrossRef] [PubMed]
  37. Al Azzaz, J.; Rieu, A.; Aires, V.; Delmas, D.; Chluba, J.; Winckler, P.; Bringer, M.A.; Lamarche, J.; Vervandier-Fasseur, D.; Dalle, F.; et al. Resveratrol-Induced Xenophagy Promotes Intracellular Bacteria Clearance in Intestinal Epithelial Cells and Macrophages. Front. Immunol. 2018, 9, 3149. [Google Scholar] [CrossRef]
  38. Lee, I.K.; Bae, S.; Gu, M.J.; You, S.J.; Kim, G.; Park, S.-M.; Jeung, W.-H.; Ko, K.H.; Cho, K.J.; Kang, J.S. H9N2-specific IgG and CD4+ CD25+ T cells in broilers fed a diet supplemented with organic acids. Poult. Sci. 2017, 96, 1063–1070. [Google Scholar] [CrossRef] [PubMed]
  39. Koskela, K.; Arstila, T.P.; Lassila, O. Costimulatory function of CD28 in avian gammadelta T cells is evolutionarily conserved. Scand. J. Immunol. 1998, 48, 635–641. [Google Scholar] [CrossRef] [PubMed]
  40. Munoz, I.; Berges, M.; Bonsergent, C.; Cormier-Aline, F.; Quéré, P.; Sibille, P. Cloning, expression and functional characterization of chicken CCR6 and its ligand CCL20. Mol. Immunol. 2009, 47, 551–559. [Google Scholar] [CrossRef] [PubMed]
  41. Wu, Z.; Hu, T.; Kaiser, P. Chicken CCR6 and CCR7 are markers for immature and mature dendritic cells respectively. Dev. Comp. Immunol. 2011, 35, 563–567. [Google Scholar] [CrossRef] [PubMed]
  42. Rosnet, O.; Blanco-Betancourt, C.; Grivel, K.; Richter, K.; Schiff, C. Binding of free immunoglobulin light chains to VpreB3 inhibits their maturation and secretion in chicken B cells. J. Biol. Chem. 2004, 279, 10228–10236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Annamalai, T.; Selvaraj, R.K. Chemokine receptor CCR7 and CXCR5 mRNA in chickens following inflammation or vaccination. Poult. Sci. 2011, 90, 1695–1700. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, G.; Zhang, Y.; Song, X.; Xia, Y.; Lai, P.F.; Ai, L. Lactobacillus casei LC2W can inhibit the colonization of Escherichia coli O157:H7 in vivo and reduce the severity of colitis. Food Funct. 2019, 10, 5843–5852. [Google Scholar] [CrossRef]
  45. Redweik, G.A.J.; Stromberg, Z.R.; Van Goor, A.; Mellata, M. Protection against avian pathogenic Escherichia coli and Salmonella Kentucky exhibited in chickens given both probiotics and live Salmonella vaccine. Poult. Sci. 2020, 99, 752–762. [Google Scholar] [CrossRef] [PubMed]
  46. Itoh, N.; Ornitz, D.M. Fibroblast growth factors: From molecular evolution to roles in development, metabolism and disease. J. Biochem. 2011, 149, 121–130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Matsubara, Y.; Aoki, M.; Endo, T.; Sato, K. Characterization of the expression profiles of adipogenesis-related factors, ZNF423, KLFs and FGF10, during preadipocyte differentiation and abdominal adipose tissue development in chickens. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2013, 165, 189–195. [Google Scholar] [CrossRef] [PubMed]
  48. Kliewer, S.A.; Mangelsdorf, D.J. Bile Acids as Hormones: The FXR-FGF15/19 Pathway. Dig. Dis. 2015, 33, 327–331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Rulifson, I.C.; Collins, P.; Miao, L.; Nojima, D.; Lee, K.J.; Hardy, M.; Gupte, J.; Hensley, K.; Samayoa, K.; Cam, C.; et al. In Vitro and In Vivo Analyses Reveal Profound Effects of Fibroblast Growth Factor 16 as a Metabolic Regulator. J. Biol. Chem. 2017, 292, 1951–1969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Okamoto, M.; Bito, T.; Noji, S.; Ohuchi, H. Subtype-specific expression of Fgf19 during horizontal cell development of the chicken retina. Gene Expr. Patterns 2009, 9, 306–313. [Google Scholar] [CrossRef]
  51. Olaya-Sánchez, D.; Sánchez-Guardado, L.; Ohta, S.; Chapman, S.C.; Schoenwolf, G.C.; Puelles, L.; Hidalgo-Sánchez, M. Fgf3 and Fgf16 expression patterns define spatial and temporal domains in the developing chick inner ear. Brain Struct. Funct. 2017, 222, 131–149. [Google Scholar] [CrossRef] [Green Version]
  52. Huang, Y.; Tsai, H.; Wang, S.; Kuo, S.; Fong, Y.; Tang, C. Amphiregulin promotes vascular endothelial growth factor-C expression and lymphangiogenesis through STAT3 activation in human chondrosarcoma cells. Cell Physiol. Biochem. 2019, 52, 1–15. [Google Scholar]
  53. Busser, B.; Sancey, L.; Brambilla, E.; Coll, J.L.; Hurbin, A. The multiple roles of amphiregulin in human cancer. Biochim. Biophys. Acta 2011, 1816, 119–131. [Google Scholar] [CrossRef] [PubMed]
  54. Yin, J.; Kim, T.-H.; Park, N.; Shin, D.; Choi, H.I.; Cho, S.; Park, J.B.; Kim, J.H. TRIM71 suppresses tumorigenesis via modulation of Lin28B-let-7-HMGA2 signaling. Oncotarget 2016, 7, 79854. [Google Scholar] [CrossRef]
  55. Qi, W.; Keenan, H.A.; Li, Q.; Ishikado, A.; Kannt, A.; Sadowski, T.; Yorek, M.A.; Wu, I.-H.; Lockhart, S.; Coppey, L.J. Pyruvate kinase M2 activation may protect against the progression of diabetic glomerular pathology and mitochondrial dysfunction. Nat. Med. 2017, 23, 753. [Google Scholar] [CrossRef] [PubMed]
  56. Mosoian, A. Intracellular and extracellular cytokine-like functions of prothymosin α: Implications for the development of immunotherapies. Future Med. Chem. 2011, 3, 1199–1208. [Google Scholar] [CrossRef] [PubMed]
  57. Calderwood, S.K.; Khaleque, M.A.; Sawyer, D.B.; Ciocca, D.R. Heat shock proteins in cancer: Chaperones of tumorigenesis. Trends Biochem. Sci. 2006, 31, 164–172. [Google Scholar] [CrossRef]
  58. Hassan, F.-U.; Nawaz, A.; Rehman, M.S.; Ali, M.A.; Dilshad, S.M.; Yang, C. Prospects of HSP70 as a genetic marker for thermo-tolerance and immuno-modulation in animals under climate change scenario. Anim. Nutr. 2019, 5, 340–350. [Google Scholar] [CrossRef] [PubMed]
  59. Xu, J.; Tang, S.; Song, E.; Yin, B.; Bao, E. Inhibition of heat shock protein 70 intensifies heat-stressed damage and apoptosis of chicken primary myocardial cells in vitro. Mol. Med. Rep. 2017, 15, 2881–2889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Parcells, M.S.; Lin, S.-F.; Dienglewicz, R.L.; Majerciak, V.; Robinson, D.R.; Chen, H.-C.; Wu, Z.; Dubyak, G.R.; Brunovskis, P.; Hunt, H.D. Marek’s disease virus (MDV) encodes an interleukin-8 homolog (vIL-8): Characterization of the vIL-8 protein and a vIL-8 deletion mutant MDV. J. Virol. 2001, 75, 5159–5173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Wu, Y.F.; Shien, J.H.; Yin, H.H.; Chiow, S.H.; Lee, L.H. Structural and functional homology among chicken, duck, goose, turkey and pigeon interleukin-8 proteins. Vet. Immunol. Immunopathol. 2008, 125, 205–215. [Google Scholar] [CrossRef] [PubMed]
  62. Withanage, G.; Kaiser, P.; Wigley, P.; Powers, C.; Mastroeni, P.; Brooks, H.; Barrow, P.; Smith, A.; Maskell, D.; McConnell, I. Rapid expression of chemokines and proinflammatory cytokines in newly hatched chickens infected with Salmonella enterica serovar typhimurium. Infect. Immun. 2004, 72, 2152–2159. [Google Scholar] [CrossRef] [Green Version]
  63. Ateya, A.I.; Arafat, N.; Saleh, R.M.; Ghanem, H.M.; Naguib, D.; Radwan, H.A.; Elseady, Y. Intestinal gene expressions in broiler chickens infected with Escherichia coli and dietary supplemented with probiotic, acidifier and synbiotic. Vet. Res. Commun. 2019, 43, 131–142. [Google Scholar] [CrossRef] [PubMed]
  64. Ohandjo, A.Q.; Liu, Z.; Dammer, E.B.; Dill, C.D.; Griffen, T.L.; Carey, K.M.; Hinton, D.E.; Meller, R.; Lillard, J.W. Transcriptome network Analysis Identifies CXCL13-CXCR5 Signaling Modules in the prostate tumor immune Microenvironment. Sci. Rep. 2019, 9, 1–13. [Google Scholar] [CrossRef] [PubMed]
  65. Kazanietz, M.G.; Durando, M.; Cooke, M. CXCL13 and Its Receptor CXCR5 in Cancer: Inflammation, Immune Response, and Beyond. Front. Endocrinol. 2019, 10, 471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Haertle, S.; Alzuheir, I.; Busalt, F.; Waters, V.; Kaiser, P.; Kaufer, B.B. Identification of the receptor and cellular ortholog of the Marek’s Disease Virus (MDV) CXC Chemokine. Front. Microbiol. 2017, 8, 2543. [Google Scholar] [CrossRef] [Green Version]
  67. Schneider, K.; Kothlow, S.; Schneider, P.; Tardivel, A.; Göbel, T.; Kaspers, B.; Staeheli, P. Chicken BAFF—A highly conserved cytokine that mediates B cell survival. Int. Immunol. 2004, 16, 139–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Kothlow, S.; Morgenroth, I.; Graef, Y.; Schneider, K.; Riehl, I.; Staeheli, P.; Schneider, P.; Kaspers, B. Unique and conserved functions of B cell-activating factor of the TNF family (BAFF) in the chicken. Int. Immunol. 2007, 19, 203–215. [Google Scholar] [CrossRef] [Green Version]
  69. Han, X.; Chen, T.; Wang, M. Molecular cloning and characterization of chicken interferon-gamma receptor alpha-chain. J. Interferon Cytokine Res. 2008, 28, 445–454. [Google Scholar] [CrossRef] [PubMed]
  70. Kaiser, P.; Stäheli, P. Avian cytokines and chemokines. In Avian Immunology; Elsevier: Amsterdam, The Netherlands, 2014; pp. 189–204. [Google Scholar]
  71. Kidane, F.A.; Mitra, T.; Wernsdorf, P.; Hess, M.; Liebhart, D. Allocation of Interferon Gamma mRNA Positive Cells in Caecum Hallmarks a Protective Trait Against Histomonosis. Front. Immunol. 2018, 9, 1164. [Google Scholar] [CrossRef] [Green Version]
  72. Bandurska, K.; Krol, I.; Myga-Nowak, M. Interferons: Between structure and function. Postepy Hig. Med. Dosw. (Online) 2014, 68, 428–440. [Google Scholar] [CrossRef] [PubMed]
  73. Brisbin, J.T.; Gong, J.; Orouji, S.; Esufali, J.; Mallick, A.I.; Parvizi, P.; Shewen, P.E.; Sharif, S. Oral treatment of chickens with lactobacilli influences elicitation of immune responses. Clin. Vaccine Immunol. 2011, 18, 1447–1455. [Google Scholar] [CrossRef]
  74. Förster, R.; Davalos-Misslitz, A.C.; Rot, A. CCR7 and its ligands: Balancing immunity and tolerance. Nat. Rev. Immunol. 2008, 8, 362–371. [Google Scholar] [CrossRef] [PubMed]
  75. Ou, C.; Wang, Q.; Zhang, Y.; Kong, W.; Zhang, S.; Yu, Y.; Ma, J.; Liu, X.; Kong, X. Transcription profiles of the responses of chicken bursae of Fabricius to IBDV in different timing phases. Virol. J. 2017, 14, 93. [Google Scholar] [CrossRef] [PubMed]
  76. Wang, Q.; Ou, C.; Wei, X.; Yu, Y.; Jiang, J.; Zhang, Y.; Ma, J.; Liu, X.; Zhang, G. CC chemokine ligand 19 might act as the main bursal T cell chemoattractant factor during IBDV infection. Poult. Sci. 2019, 98, 688–694. [Google Scholar] [CrossRef] [PubMed]
  77. Haghighi, H.R.; Abdul-Careem, M.F.; Dara, R.A.; Chambers, J.R.; Sharif, S. Cytokine gene expression in chicken cecal tonsils following treatment with probiotics and Salmonella infection. Vet. Microbiol. 2008, 126, 225–233. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 09 00382 g001 550
Figure 1. Transcriptome profiles of ileal tissues from 21-day-old broilers. (a) Correlation coefficient heat map between individual samples. The value is the square of the correlation coefficient of Pearson (R2) (p < 0.01). (b) Principal component analysis of transcriptomes from control group and probiotics group. (c) Venn diagram of the number of expressed genes detected in control group and probiotics group. “P” = Bacillus amyloliquefaciens TL group; “C” = control group.
Figure 1. Transcriptome profiles of ileal tissues from 21-day-old broilers. (a) Correlation coefficient heat map between individual samples. The value is the square of the correlation coefficient of Pearson (R2) (p < 0.01). (b) Principal component analysis of transcriptomes from control group and probiotics group. (c) Venn diagram of the number of expressed genes detected in control group and probiotics group. “P” = Bacillus amyloliquefaciens TL group; “C” = control group.
Microorganisms 09 00382 g001
Microorganisms 09 00382 g002 550
Figure 2. Gene Ontology (GO) terms of ileal transcriptomes. (a) GO enrichment analysis of the all expressed genes scattered in the C group (top 60). (b) GO enrichment analysis of all expressed genes scattered in the P group (top 60). (c) Venn diagram of the first 60 GO terms of Groups C and P. “BP” = biological process, “CC” = cellular component, “MF” = molecular function, “C” = control group, and “P” = Bacillus amyloliquefaciens TL group.
Figure 2. Gene Ontology (GO) terms of ileal transcriptomes. (a) GO enrichment analysis of the all expressed genes scattered in the C group (top 60). (b) GO enrichment analysis of all expressed genes scattered in the P group (top 60). (c) Venn diagram of the first 60 GO terms of Groups C and P. “BP” = biological process, “CC” = cellular component, “MF” = molecular function, “C” = control group, and “P” = Bacillus amyloliquefaciens TL group.
Microorganisms 09 00382 g002
Microorganisms 09 00382 g003 550
Figure 3. Functional enrichment analysis of uniquely expressed genes (a) Gene ontology (GO) enrichment analysis of the uniquely expressed genes scattered in the C group. (b) GO enrichment analysis of the uniquely expressed genes scattered in the P group. (c) Enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed genes (DEGs) in the C group. (d) Enriched KEGG pathways of DEGs in the P group. “BP” = biological process, “CC” = cellular component, “MF” = molecular function, “C” = control group, and “P” = Bacillus amyloliquefaciens TL group.
Figure 3. Functional enrichment analysis of uniquely expressed genes (a) Gene ontology (GO) enrichment analysis of the uniquely expressed genes scattered in the C group. (b) GO enrichment analysis of the uniquely expressed genes scattered in the P group. (c) Enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed genes (DEGs) in the C group. (d) Enriched KEGG pathways of DEGs in the P group. “BP” = biological process, “CC” = cellular component, “MF” = molecular function, “C” = control group, and “P” = Bacillus amyloliquefaciens TL group.
Microorganisms 09 00382 g003
Microorganisms 09 00382 g004 550
Figure 4. Functional enrichment analysis of differentially expressed genes (DEGs). (a) Gene ontology (GO) enrichment analysis of the upregulated DEGs between C group and p group. (b) GO enrichment analysis of the downregulated DEGs between C group and P group. “BP” = biological process, “CC” = cellular component, and “MF” = molecular function. (c) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of upregulated DEGs. (d) KEGG pathway enrichment of downregulated DEGs. The vertical axis represents the name of the pathway, and the horizontal axis represents the Rich factor. The size of the point indicates the number of DEGs in the pathway, and the color of the point corresponds to a different Q value ranges. “P” = Bacillus amyloliquefaciens TL group; “C” = control group.
Figure 4. Functional enrichment analysis of differentially expressed genes (DEGs). (a) Gene ontology (GO) enrichment analysis of the upregulated DEGs between C group and p group. (b) GO enrichment analysis of the downregulated DEGs between C group and P group. “BP” = biological process, “CC” = cellular component, and “MF” = molecular function. (c) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of upregulated DEGs. (d) KEGG pathway enrichment of downregulated DEGs. The vertical axis represents the name of the pathway, and the horizontal axis represents the Rich factor. The size of the point indicates the number of DEGs in the pathway, and the color of the point corresponds to a different Q value ranges. “P” = Bacillus amyloliquefaciens TL group; “C” = control group.
Microorganisms 09 00382 g004
Microorganisms 09 00382 g005 550
Figure 5. The fold expression changes in six randomly selected genes to validate the RNA-seq results.
Figure 5. The fold expression changes in six randomly selected genes to validate the RNA-seq results.
Microorganisms 09 00382 g005
Microorganisms 09 00382 g006 550
Figure 6. B. amyloliquefaciens TL acts on broiler ileum tissues and significantly downregulates immune-related genes in the “cytokine–cytokine receptor interaction” pathway. (a) “cytokine–cytokine receptor interaction” pathway. The genes in the green boxes represent genes expressed in birds, and the genes in the purple box represent the downregulated genes related to immune activity obtained in this study. (b) Schematic diagram of the function of the cytokine or cytokine receptor encoded by the downregulated gene. The interaction between the chemokine CXCL13 and its sole receptor CXCR5 guides B cells to B cell follicles and participates in the zoning of germinal centers. It also attracts T cells to B cell areas and mediates the interaction between B and T cells [66]; the interaction of BAFF and its receptor BAFFR (encoded by TNFSF13C) plays a central mediator in the homeostasis of mature B cell subsets, mainly through the modulation of cell survival. Unlike mammals, BAFF is also involved in early B cell development in birds [68]; the interaction of IL-8 with IL8RA, a G-protein-coupled receptor is responsible for the activation of neutrophils. An IL-8-activated neutrophil is characterized by enhanced migration, phagocytosis, superoxide generation, and granule release [60]; IFN-γ interacts with its receptor IFNGR1, activates type I macrophages, enhances the activation of both MHC-I and MHC-II molecules, and assists in antigen presentation and processing, as well as removal of intracellular pathogens, and prevents virus replication. Moreover, IFN-γ activates T-helper I-type immune responses [70]; the chemokine CCL19 interacts with its receptor CCR7 and plays a pivotal role in T cell and dendritic cell trafficking into Bursa of Fabricius and other lymphoid tissues [76]. B. amyloliquefaciens TL acts directly or indirectly (indicated by the pink dotted arrow in the figure) on immune cells to reduce the expression of corresponding cytokines or cytokine receptors, thereby downregulating the above-mentioned immune activity and reducing the energy consumption of immune activation, which is beneficial for the accumulation of energy in broilers.
Figure 6. B. amyloliquefaciens TL acts on broiler ileum tissues and significantly downregulates immune-related genes in the “cytokine–cytokine receptor interaction” pathway. (a) “cytokine–cytokine receptor interaction” pathway. The genes in the green boxes represent genes expressed in birds, and the genes in the purple box represent the downregulated genes related to immune activity obtained in this study. (b) Schematic diagram of the function of the cytokine or cytokine receptor encoded by the downregulated gene. The interaction between the chemokine CXCL13 and its sole receptor CXCR5 guides B cells to B cell follicles and participates in the zoning of germinal centers. It also attracts T cells to B cell areas and mediates the interaction between B and T cells [66]; the interaction of BAFF and its receptor BAFFR (encoded by TNFSF13C) plays a central mediator in the homeostasis of mature B cell subsets, mainly through the modulation of cell survival. Unlike mammals, BAFF is also involved in early B cell development in birds [68]; the interaction of IL-8 with IL8RA, a G-protein-coupled receptor is responsible for the activation of neutrophils. An IL-8-activated neutrophil is characterized by enhanced migration, phagocytosis, superoxide generation, and granule release [60]; IFN-γ interacts with its receptor IFNGR1, activates type I macrophages, enhances the activation of both MHC-I and MHC-II molecules, and assists in antigen presentation and processing, as well as removal of intracellular pathogens, and prevents virus replication. Moreover, IFN-γ activates T-helper I-type immune responses [70]; the chemokine CCL19 interacts with its receptor CCR7 and plays a pivotal role in T cell and dendritic cell trafficking into Bursa of Fabricius and other lymphoid tissues [76]. B. amyloliquefaciens TL acts directly or indirectly (indicated by the pink dotted arrow in the figure) on immune cells to reduce the expression of corresponding cytokines or cytokine receptors, thereby downregulating the above-mentioned immune activity and reducing the energy consumption of immune activation, which is beneficial for the accumulation of energy in broilers.
Microorganisms 09 00382 g006
Table
Table 1. Primers used in this study.
Table 1. Primers used in this study.
Name.Sequence (5′–3′)
β-actinF: TATTGCTGCGCTCGTTGTTG
R: TGGCCCATACCAACCATCAC
CXCR5F: CCTATGACTTGAGCCTGGTGG
R: TCCCAGCACGAACATAAGCAG
LAMA5F: GCACATCCCATAACGAACGC
R: TCAGGACGAGGGGAATTTGC
ACACBF: TTCGGGACTTCAACCGTGAG
R: GGCTGCTTAAAATCCCGCAG
CXCL13F: CAGCCATCCTGGAAGCCAAC
R: GGATCCACACAGATCCTCTCG
LEPRF: AGTGCAAACATGCAAGCGAG
R: CAGCTTGCCTTCAACCCAAC
Table
Table 2. Effects of Bacillus amyloliquefaciens TL on the growth performance of broilers.
Table 2. Effects of Bacillus amyloliquefaciens TL on the growth performance of broilers.
GroupControlProbioticsp-Value
1 days40.84 ± 1.4839.58 ± 1.190.065
Average weight7 days147.48 ± 11.08157.35 ± 8.580.053
(g)14 days499.75 ± 20.67532.50 ± 24.090.056
21 days858.62 ± 30.74 A1076.12 ± 98.48 B0.000
-35 days1901.38 ± 76.36 a2156.12 ± 198.61 b0.017
1–7 days15.23 ± 1.38 a16.82 ± 1.08 b0.022
Average daily 8–14 days50.32 ± 1.5153.59 ± 2.330.074
gain (g)15–21 days51.27 ± 1.74 A77.66 ± 10.74 B0.000
22–35 days80.21 ± 3.7283.08 ± 9.230.608
Feed conversion ratio2.021.62
Note: The result is expressed as Mean ± SEM (n = 60 chickens per group in a total of 120). Different capital letters between columns indicate extremely significant differences (p < 0.01); different lowercase letters between columns indicate significant differences (p < 0.05); the same lowercase letters or no letters indicate nonsignificant differences (p > 0.05).
Table
Table 3. RNA-sequencing (RNA-seq) analysis in Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway.
Table 3. RNA-sequencing (RNA-seq) analysis in Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway.
KEGG Pathwayp-ValueGene NumberMolecules
Up group
Calcium signaling pathway0.0598665PHKA1, ADCY3, P2RX1, HTR6, ITPR3
Focal adhesion0.0883085COL27A1, LAMA5, RAPGEF1, KDR, COL5A1
Insulin signaling pathway0.0890464TSC2, PHKA1, ACACB, RAPGEF1
Neuroactive ligand–receptor interaction0.0915556SSTR1, GRIK4, GRM2, CHRM4, P2RX1, HTR6
Down group
Cytokine–cytokine receptor interaction0.0012606TNFRSF13C, IFNGR1, CXCL13L2, IL20RA, CCL19, CXCR5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hong, Y.; Cheng, Y.; Guan, L.; Zhou, Z.; Li, X.; Shi, D.; Xiao, Y. Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance. Microorganisms 2021, 9, 382. https://doi.org/10.3390/microorganisms9020382

AMA Style

Hong Y, Cheng Y, Guan L, Zhou Z, Li X, Shi D, Xiao Y. Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance. Microorganisms. 2021; 9(2):382. https://doi.org/10.3390/microorganisms9020382

Chicago/Turabian Style

Hong, Yuxuan, Yang Cheng, Leluo Guan, Zutao Zhou, Xiaowen Li, Deshi Shi, and Yuncai Xiao. 2021. "Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance" Microorganisms 9, no. 2: 382. https://doi.org/10.3390/microorganisms9020382

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop