- 1Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Ministry of Education, College of Animal Science and Veterinary Medicine, Huazhong Agricultural University, Wuhan, Hubei, China
- 2Animal Production Department, Faculty of Agriculture, Moshtohor, Benha University, Mushthar, Egypt
- 3Institute of Animal Science, Guangdong Academy of Agricultural Sciences, State Key Laboratory of Livestock and Poultry Breeding, Key Laboratory of Animal Nutrition and Feed Science in South China, Ministry of Agriculture, Key Laboratory of Poultry Genetics and Breeding, Ministry of Agriculture, Guangdong Key Laboratory of Animal Breeding and Nutrition, Guangdong Public Laboratory of Animal Breeding and Nutrition, Guangzhou, China
- 4Animal and Poultry Production Department, Faculty of Technology and Development, Zagazig University, Zagazig, Egypt
- 5Department of Poultry Production, Faculty of Agriculture, Assiut University, Assiut, Egypt
- 6Department of Animal Breeding and Genetics, Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur, Pakistan
- 7Department of Animal Production, National Research Centre, Giza, Egypt
- 8Department of Nutrition and Clinical Nutrition, Faculty of Veterinary Medicine, Cairo University, Giza, Egypt
Starting phase of laying chicken life is the building stone for rearing and production stages. Since, fecal microbial transplantation (FMT) regulates the gut microbial diversity and affects the productive performance of the bird. The aim of this study is to evaluate the effect of FMT from feed-efficient broiler chicken could program the diversity of gut microbiota and growth of recipient native slow growing egg-laying chicks. For this, a total of 150 (one-day-old) Jing Hong chicks were randomly assigned into two groups, each group consisted of 5 replicates (n = 15 bird/ replicate). The control group (CON) and FMT recipient birds (FMT) fed on basal diet, the FMT group received an oral daily dose of FMT prepared from Cobb-500 chickens. The FMT performed from the 1d to 28d of age, through the experimental period, feed intake and body weight were recorded weekly. At the end of a 28-day trial, carcass traits were assessed and cecal samples were collected for microbiome assessment via 16S rRNA-based metagenomic analysis to characterize the diversity and functions of microbial communities. The data were statistically analyzed using R software. Body weight and body weight gain increased, and FCR decreased (p = 0.01) in FMT group. The relative abundance of Firmicutes and the Firmicutes/Bacteroidetes (F/B) ratio were increased due to FMT administration (p = 0.01). A higher relative abundance of Lactobacillus, Lactococcus, and Bifidobacterium were presented in the FMT group. Meanwhile, Enterococcus, Helicobacter, and Bacteroides were more abundant in the CON group (p < 0.01). Kyoto encyclopedia of genes and genomes (KEGG) pathways for microbial functions regarding amino acid metabolism, secondary metabolites biosynthesis, carbohydrate metabolism, energy metabolism, and enzyme families, cofactors, and vitamins were significantly annotated in the FMT group. Overall, FMT administration from the donor of highly feed-efficient broilers improved weight gain by reshaping a distinct gut microbiome, which may be related to the metabolism and health in the recipients laying chicks, providing new insight on the application of the FMT technique for early life programming of laying chickens.
Introduction
The early maturation of laying chicks and early body weight gain was found to be a key step to achieve the flock uniformity and egg production efficiency (Underwood et al., 2021). One key factor regulating the birds body weight is the microbiota that colonize their gut, the close symbiotic relationship between a host and its gut microbiota affects energy harvest from diet, regulates the host well-being, and promotes the innate and the acquired immunity (Deweerdt, 2014). Additionally, gut microbiota secures the permeability of the intestinal mucosa and regulates the fermentation and absorption of dietary nutrients, which may explain the microbiota importance in the regulation of energy and growth activities (Sonnenburg and Bäckhed, 2016). In our previous studies, we investigated some factors that affect gut microbiota development in chicken, including performance type (Elokil et al., 2020a), medicine (Elokil et al., 2020b), fertility (Elokil et al., 2020a), and diet (Abouelezz et al., 2019). The FMT has been reported as the most effective way for whole microbial community changes, whereas probiotic and antibiotic tools usually alter up to three orders of the 100 trillion native microorganisms of the gut (Borody and Khoruts, 2012; Cammarota et al., 2014). The FMT, also known as stool/fecal transplantation or fecal bacterio-therapy, is the infusion of a liquid filtrate of stools from a donor individual into the gut of a recipient individual to modulate the microbial community for therapeutic or medical purposes (Zhang et al., 2012) and to establish a durable alteration of the recipient’s gut microbiota (Grehan et al., 2010).
The microbial Firmicutes/Bacteroidetes (F/B) ratio is important with positive correlation of somatic growth, as these taxa include bacteria that ferment indigestible carbohydrates providing an extra source of energy for the host, which is eventually stored in skeletal muscle and adipose tissue and used for somatic growth (Cornejo-Pareja et al., 2019). Body weight selection was shown to result in quantitative genetic-correlated responses in gut microbiota and in significant differences in heritability of the F/B ratio, and the genetic correlations between high and low body weight chickens were estimated (Meng et al., 2014). Moreover, the higher F/B ratio of gut microbiota has been associated with overweight in humans (Karvonen et al., 2019), mice (Vadder de and Mithieux, 2018), and rabbits (Zeng et al., 2015). The gut microbiota of broiler chicks, which have a high F/B ratio, were found to have higher ability to ferment diet ingredients to produce volatile fatty acids than that of laying chicks (Walugembe et al., 2015). This suggests that a high F/B ratio in the gut microbial diversity and/or richness is associated with body weight gain. On the other side, the intentional change of gut microbiota through fecal microbiota transplantation (FMT), probiotics, and antibiotics has been used for several purposes including improvement of bird’s growth performance.
Jing Hong chicken was selected for this study as it is a prominent native breed in China, which is widely distributed in Hubei Province as one of the important egg varieties with high economic value of egg production. It has a high resistance for endemic diseases, rough feeding conditions, and extensive breeding environment. Therefore, the aim of the present study was to use FMT as a tool to induce early life programming of laying chicks via modulation of their gut microbiota, leading to early body weight maturation. Thereby, providing experimental evidence to understand how the intentional programming of gut microbiome affects the growth performance of laying chicks.
Materials and methods
Experimental design, diet, and bird management
In a 28-day trial. A total of 150-day-old unsexed of newly hatched chicks from Jing Hong laying strain were randomly allocated into two treatment groups (n = 75/group) as follows: (1) birds fed on the basal diet (CON), and (2) birds fed on basal diet and received fecal microbiota transplantation (FMT). Each dietary treatment consisted of 5 replicates (15 birds/ replicate). The chicks were housed in growth metal cages in an enclose poultry farm of Huazhong Agricultural University throughout the study period, each cage was equipped with one manual feeder and drinker to ensure the animals’ nutritional requirements. A tray with a parchment paper was placed under each cage for fecal sample collection. The diet was formulated following the guidelines of Chinese Feeding Standards using the Feed database in China (Beijing, 2004), diet ingredients and nutrient composition is provided in Table 1. All birds were fed on the same iso-nitrogenous and iso-caloric diet and water was provided ad-libitum.
Preparation and inoculation of the FMT material
Seventy-five newly hatched Cobb-500 broiler chicks were used daily as a stool donor. The FMT was prepared following the methods (Matsuoka et al., 2014). Briefly, 10 g of fresh fecal samples were collected daily in the morning into the sterile tube (50 ml). The white part of the excreta was removed immediately, because it mainly comprises uric acid, and the fecal part was mixed with 0.75% saline in 1:6 ratios (6 ml of 0.75% saline for each gram of feces) and span at low speed (800 × g) for 3 min at 4°C to separate undigested feed and particulate material from the microbial fraction. Keeping the mixture on ice until precipitates were fully settled down, the supernatant was collected and filtered with the sterile gauze to get fecal suspension.
Performing FMT in laying chickens
A total of 150 females of the recipients Jing Hong chicks similar age with the donor Cobb-500 (one-day-old) were randomly assigned into two groups (FMT and CON), each group consisted from 75 chicks into 5 replicates (n = 15 bird/ replicate). Each laying chick in FMT group was provided with an oral fresh dosage of 1 ml of FMT solution daily from day 1 to day 28 of age, at early morning and before having access to feed and water. The transplant was orally administered using 1-ml oral dosage via oral-syringes with soft flexible tips that were positioned at the back of the chicks’ tongues. The chicks were then supervised to ensure that they had swallowed the material and feed was with-held for 60 min after administration.
Growth performance and carcass traits
Live body weight was recorded weekly and used to determine the weekly body weight gain and feed intake was recorded weekly, FCR was then calculated. At 28 days of age, 11 healthy chicks of group average BW were selected and euthanized. The chicks’ breasts (pectoralis major) and thigh (drumstick) muscles were skinned, de-boned and weighted; their livers weight, and the lengths (cm) of their duodenum, jejunum, ileum, and ceca were measured similar to Abdelatty et al. (2021).
Microbial DNA extraction
To obtain samples for microbial genes analysis, fresh samples from the cecum (n = 11/ group) were aseptically collected 5–10 min after slaughter under anaerobic conditions and placed into sterile 2-ml cryo-tubes (Sarstedt, Nümbrecht, Germany), then, samples were immediately stored at −80°C until total microbial DNA extraction. Then, all the sample tubes were instantaneously snap-frozen in liquid nitrogen, and subsequently stored at −80°C for extraction of total microbial DNA later.
The microbial DNA was extracted from 300-μl cecal samples using a DNA stool mini kit (QIAmp DNA Stool mini Kit; QIAGEN, Hilden, Germany) following the manufacturer’s instructions. A Qubit 2.0 Fluorometer (Thermo Fisher Scientific Co., Ltd., Shanghai, China) with a Qubit dsDNA HS assay kit (Life Technologies) was used to quantify DNA concentration.
PCR amplification and 16S rRNA sequencing
An aliquot of each of the extracted DNA samples was used as a template for PCR amplification. Additionally, 16S ribosomal RNA, library preparation, and DNA sequencing were performed by a commercial provider (Personalbio Co., Ltd., Shanghai, China). To generate an approximate amplicon size of 570 bp, primers (forward 5′ CCTAYGGGRBGCASCAG GNG 3′, reverse 5′ GGACTACNNGGGTATCTAAT 3′) targeting V3-V4 region amplicons were used for amplification and PCR products were purified (Busato et al., 2022; Tawfik et al., 2022). The PCR conditions were as follows: initial denaturation, annealing, and extension were carried out and repeated at 94°C for 4 min, 94°C for 30 s, 50°C for 45 s, and 72°C for 30 s for 25 cycles.
Sequence quality control and calculation of operational taxonomic units (OTUs)
For each library, amplified libraries of V3-V4 region amplicons were pooled and sequenced using the Illumina Miseq 2000 platform sequencer, including 250-bp paired-end reads generated with a 7-cycle index read (sequences with an overlap longer than 10 bp without any mismatch were assembled). The resultant overlapping paired-end reads were stitched and quality filtered using Microsynth after removing barcode and primer sequences. To perform sequence quality control, the pre-filtered and stitched reads were processed using the QIIME software (Quantitative Insights into Microbial Ecology, v1.8.0).1 Briefly, the raw tags were merged based on the overlap of two reads. Then, clean tags were created by pretreating the raw tags and removing the chimeric sequences to generate effective tags. Finally, quantitative analysis of the sequences was performed using the QIIME software, and the chimeric sequences were eliminated using the USEARCH software (v5.2.236).2
Bioinformatics analysis
Quality sequences were counted for each sample after the removal of the chimeric sequences (shorter than 160 bp), and statistical estimations were created for the distribution of sequence length using the R software to characterize the length distribution of the sequences contained in each sample. Then, sequences were clustered into OTUs with 97% similarity using open-reference OTU picking and UCLUST (Edgar, 2010). In addition, Specaccum analysis was applied to check if all sample sizes and the OTU abundance matrix were sufficient to estimate community richness (Supplementary Figure S1). Finally, the taxonomy of OTUs was determined using the Greengenes default database in QIIME (DeSantis et al., 2006). The most abundant OTUs between groups (CON and FMT) were then classified against the NCBI nucleotide database using BLASTN for taxonomic classification, targeting the 16S rRNA marker. For the analysis of bacteria and archaea, the databases for the 16S rRNA gene, Greengenes (Release 13.8)3 and RDP (Ribosomal Database Project, Release 11.1),4 were used by default to identify the OTU diversity among the samples and between the groups. A Venn diagram was constructed to calculate the total number of OTUs per sample (i.e., group) using the R software.
Annotation of microbial composition
A rarefaction depth of 12,000 sequences was used to perform α- and β-diversity analyses among samples; the α-diversity was determined using the unweighted and weighted UniFrac distance (Lozupone et al., 2007). The analysis of α-diversity (Shannon, Simpson, Chao1, and ACE indices) was applied to find the drivers of variation in the microbial community structure in the samples (Shannon, 1948; Simpson, 1949; Chao and Shen, 2004). In addition, rarefaction curves of β-diversity (PCA, PCoA, NMDS, and UPGMA cluster analyses) in the samples were calculated using a maximum rarefaction depth of 12,000 sequences and the observed OTU index (Ramett, 2007). The orders of phylogenetic tree construction (MEGAN and KRONA) were used for microbial interactive visualization and to visually display the results of species annotation between groups based on OTU tags. Then, a heat map was constructed according to the distribution of the top 50 abundances and the degree of similarity between the samples. The order analyses of PLS-DA, LDA, PERMANOVA, and ANOSIM were performed to determine the variation in community structure between groups for screening key species.
Annotation of microbial function prediction
To predict the metabolic functioning of bacteria based on total genome sequences from the 16S rRNA gene, Kyoto Encyclopedia of Genes and Genomes (KEGG) function annotation of the sequence was carried out based on Tax4Fun, then visualized using the statistical analysis of metagenomic profiles (STAMP) software package (Langille et al., 2013). KEGG orthologies (KOs) were categorized into KEGG level 2 pathways with the removal of the non-microbial pathways. The classified genes were divided into six categories, namely, Metabolism, Genetic Information Processing, Environmental Information Processing, Cellular Processes, Organismal Systems, and Human Diseases, each of which was further divided into multiple levels. The nearest sequenced taxon index (NSTI) was used to provide a measure of the availability of sequenced reference genomes for each OTU in a sample. It calculates the average branch length for each sample between an OTU and the nearest sequenced reference genome based on the Greengenes reference phylogeny, accounting for OUT abundance. Finally, the R software was used for the cluster analysis of the top 50 most abundant functional predictions in each sample, presented by the heat map (Langille et al., 2013).
Statistical analysis
The analysis of β-diversity included PCA and PCoA, which were calculated using weighted UniFrac and unweighted UniFrac. One-way ANOVA and the LM procedure of the R software version 3.2.2, R Core were applied to estimate the different bacterial taxa in CON and FMT groups (R Development Core Team, 2015). Each chick individual was considered an experimental unit, and FMT was included as a fixed effect in the statistical model. All differences were considered significant at p < 0.05 and were considered trends when p < 0.10. Pairwise comparison was performed using T tests.
Results
Growth performance and carcass traits
The differences in body weight, daily weight gain, internal organs weight, and length of the digestive tract sections between CON and FMT groups of laying chickens are shown in Table 2. FMT increased final body weight (p < 0.01) by 2.5% and average daily gain increased (p < 0.01) by 2.3%. Feed conversion ratio was lower in FMT group compared with the CON group (p = 0.01), breast muscle weight increased by 9.7% and thigh muscle weight increased by 17.63% in FMT chickens compared with the control counterpart (p = 0.01). Liver weight increased by 8.5% (p = 0.01). On the other side, duodenum, ilium, and cecum length was shorter in MFT group compared with the CON one (p < 0.05).
Table 2. Effect of oral fecal microbiota transplantation from the donor commercial broiler chickens to the recipient native Jing Hong laying chickens on body weight and some carcass traits1.
16S rRNA sequencing metrics and quality
High-throughput sequencing generated from 22 chickens (n = 11/ dietary treatment) yielded a total of 775,945 reads (average of 35,270 and range of 32,836–38,767 reads per sample; Supplementary Figure S2). Average read length was 160 bp, and the distributions of sequence lengths shown in OTUs were generated based on the greengene database using quantitative insights into microbial ecology (QIIME) and characterized for different taxonomic levels including domain, phylum, class, order, family, and genus. The taxa considered as common in the samples were used in further analysis. Statistical number of OTUs at each classification level among groups of control and antibiotics-exposed chicks are presented in Supplementary Table S1. A total of 12 phyla, 23 classes, 35 orders, 60 families, 77 genera, and 31 species were identified in these samples Supplementary Table S2.
Alpha and beta diversities
Total observed OTU count, α indexes (ACE and Chao1), and diversity indicators (Simpson and Shannon) between groups are summarized in Table 3. The total observed OTUs was significantly increased (p < 0.01) in the FMT group than in the CON group. In contrast, the means of ACE and Chao1 were significantly decreased (p < 0.01) in the FMT group than in the control group. Similarly, the Simpson and Shannon indices were significantly decreased in the FMT group than in the CON group. In addition, β-diversity indicators, principal component analysis (PCA), non-metric multidimensional scaling (NMDS), and principal coordinate analysis (PCoA) were obtained to measure the between-group and inter-group distance of the CON and FMT groups.
Table 3. Observed operational taxonomic units (OTUs) and α-diversity measures of bacterial communities between chicks from the control and fecal microbiota transplantation (FMT) groups.
The PCA based on unweighted UniFrac distance is presented in Figure 1. The PCA plot revealed that the distance between samples of the CON group were smaller than that of the FMT group, indicating that FMT had a more diverse gut microbiota. The R value of the analysis of similarities (ANOSIM; 0.38) was correlated among samples within the CON and FMT groups (p < 0.01). Similarly, the weighted and unweighted values of NMDS indicated that samples of the CON group formed a tighter cluster than those of the FMT group. Regarding the R value obtained from the mutational multivariate analysis of variance (PERMANOVA), there were significant correlations in weighted (R = 0.25; p < 0.05) and unweighted (R = 0.19; p < 0.01) UniFrac distances among samples and between groups of the CON and FMT groups. In addition, the lowest group distance was significantly estimated (p < 0.05) within groups in the comparative analysis boxplot.
Figure 1. The PCA (principal component analysis) based on unweighted UniFrac distance. (A) the PCA, (B) non-metric multidimensional scaling (NMDS) of weighted and unweighted of unifrac distance, and (C) boxplot for comparative analysis of inter-group/group differences in unifrac distance between groups of control and FMT chicks. Two independent cluster showed significant correlated by ANOSIM (R = 0.38***) and PERMANOVA of weighted (R = 0.25**) and unweighted (R = 0.19***) of unifrac distance between control and FMT of laying-chicks.
Microbial structure and composition
Figure 2 illustrates the effect of FMT on gut microbiota structure of Jing Hong chickens. The number of OTUs was significantly increased (p < 0.01) in the FMT group, with a high abundance of unique OTUs when compared with the CON group. Additionally, a heat map of the top 50 most abundant genera in the microbiome community combined with a cluster analysis revealed similar microbiome composition among samples of the CON and FMT groups. The differences between the CON and FMT groups regarding the relative abundance (percent of total sequences) of each phylum and family of taxa from the microbiome are presented in Figure 2 and Table 4.
Figure 2. Impact of FMT on the Chickens Gut Microbiota Structure. (A) Venn diagram of OUT number between groups (Control and FMT) of common and unique units. (B) Heat map showing the genera with significant differences of relative abundances between groups. (C) Taxonomic composition analysis between groups at each classification level of phylum and family control and FMT groups.
Table 4. Relative abundance (percent of total sequences) of taxa in the cecal microbiota of chicks from the control and fecal microbiota transplantation (FMT) groups estimated using metagenomic analysis1.
At the phylum level, the relative abundances of Firmicutes increased and that of Bacteroidetes decreased (p < 0.01) in the FMT group, and the Firmicutes/Bacteroidetes ratio (F/B ratio) increased in the FMT group than in the CON group (p < 0.01; Table 4). The differences between the CON and FMT groups regarding the relative abundance (percent of total sequences) of each phylum and family of taxa from the microbiome are presented in Figure 2 and Table 4. At the phylum level, the relative abundances of Firmicutes was increased (FMT vs. CON, 89.05 vs. 75.39%) and that of Bacteroidetes was decreased (FMT vs. CON, 2.51 vs. 5.80%) in the FMT group (p < 0.05). Consequently the Firmicutes/Bacteroidetes ratio (F/B ratio) increased (FMT vs. CON, 35.45 vs. 12.99%) in the FMT group than in the CON group (p < 0.01; Table 4). On the other hand, the relative abundance of Proteobacteria (FMT vs. CON, 3.52 vs. 7.98%) and Campilobacterota (FMT vs. CON, 0.164 vs. 0.811%) were decreased in the FMT group chickens (p < 0.01; Table 4). Similarly, there was decreased (p < 0.05) in the relative abundance of phyla of Verrucomicrobia (FMT vs. CON, 0.740 vs. 1.36%) and TM7 (FMT vs. CON, 0.005 vs. 0.042%) in the FMT group than in the CON group as presented in Figure 2 and Table 4.
At the genus level, there were significant increased the relative abundances of some symbiotic taxa such as Lactobacillus (FMT vs. CON, 17.87 vs. 9.12%), Lactococcus (FMT vs. CON, 2.51 vs. 5.80%) and Bifidobacterium (FMT vs. CON, 0.064 vs. 0.051%) in the FMT group compared with the CON group (p < 0.01; Figure 2; Table 4). On the other hand, there were significant decreased the relative abundances of some pathogenic taxa such as Enterococcus (FMT vs. CON, 0.55 vs. 0.86%), Helicobacter (FMT vs. CON, 1.027 vs. 1.589%) and Bacteroides (FMT vs. CON, 0.054 vs. 1.30%) in the FMT group in compared with the CON group (p < 0.01; Figure 2; Table 4).
Reshaping gut microbial community by FMT in laying chicks
The results of the LEfSe analysis based on the linear discriminant analysis (LDA) of the specific bacterial taxa that are associated with FMT are presented in Supplementary Figure S3. The LEfSe cladogram representative of the structure of the host–microbiota axis shows a significant shift of microbiota between groups, including a total of 64 bacterial taxa that were significantly different between the CON and FMT groups, with only seven taxa (Lactobacillallus, Lactobacillaceae, Lactobacillales, Bacilli, Dietzia, Dietziaceae, and Aeriscardovia) in the FMT group. The LDA score plot presented group-enriched taxa that were significant between the two groups (p < 0.05).
Figure 3 illustrates the comparison of abundance among groups at the levels of phylum and genus were performed using a Metastats analysis and indicated an increase of the phyla Firmicutes and Acidobacteria and of the genera Dietzia and Ruminococcus in the FMT group. The similarity between samples was evaluated by unweighted pair group method with arithmetic mean (UPGMA) analysis, which showed that the samples of the CON and FMT groups were mostly similar and distributed into a different cluster of hierarchical trees (Supplementary Figure S4).
Figure 3. Metastats comparison statistics test for significant differences of taxa abundance between groups of control and FMT chicks. (A) six phylum (B) twenty genus presented higher F/B ratio in FMT compare with control group.
Microbial function prediction
The relative abundance comparison between the CON and FMT groups based on the second-level of KEGG pathways analysis is shown in Table 5 and Figure 4. A total of six categories of cellular processes, environmental information processing, genetic information processing, immune information processing, metabolism processing, and organismal systems processing were revealed between the CON and FMT groups.
Table 5. Relative abundance of the functional microbiome predicted by KEGG second-level pathways between chicks from the control and FMT groups estimated using the statistical analysis of metagenomic profiles.
Figure 4. Effect of FMT on Functional Prediction of KEGG Second-level Analysis. (A) KEGG orthologous gene cluster abundance heat map combined with cluster analysis based on the similarity of the functional group abundance distribution between samples (red represents a functional group with higher abundance in the corresponding sample, and green represents a functional group with lower abundance); venn diagram of the common and unique functional group, (B) PLS-DA discriminant analysis plot showed samples belonging to the same group are closer to each other, and the distance between the points of different groups is farther.
In the category of cellular processes, the levels of cell growth and death, and transport and catabolism had a higher (p < 0.01) relative abundance in the FMT group than in the CON group. The level of membrane transport pathway in the category of environmental information processing was significantly increased (p < 0.01) in the FMT group than in the CON group (Table 5). Prediction of functional microbiome in the category of genetic information processing showed that three pathway levels of DNA replication and repair, transcription, and translation were significantly increased (p < 0.01) in FMT group. One level of infectious diseases pathway under the category of immune information processing was increased (p < 0.05) in the FMT group when compared with the CON group. Regarding the category of metabolism processing, the differences between groups (FMT and CON) in the pathway levels regarding amino acid metabolism (p < 0.01), secondary metabolites biosynthesis (p < 0.01), carbohydrate metabolism (p < 0.01), energy metabolism (p < 0.05), enzyme families (p < 0.05), and cofactors and vitamins (p < 0.05) were higher in the FMT groups than in the CON group. The pathway level of endocrine system was significantly increased (p < 0.01) in the FMT group when compared with the CON group.
The top 50 most abundant KEGG orthologous genes were clustered in the heat map combined with the analysis between the CON and FMT groups; samples of the control group formed clusters based on their similar microbial compositions. The Venn diagram of the common and unique functional groups revealed that the FMT group had more functional groups than the CON group. PLS discriminant analysis (PLS-DA) plot showed that samples belonging to the same group were more similar to each other than those from different groups.
Discussion
There is increasing interest in early life programming of chickens for better health condition and performance. The recent reports showed promising results on the application of FMT in early life of birds (Siegerstetter et al., 2018; Yan et al., 2021) to modulate gut microbiota. Therefore, we investigated the safety of daily dosage of FMT from broiler donor chicks and their ability to alter the composition, structure, and function of gut microbiota and improve metabolic outcomes in the early life of native egg laying chicks with low performance. The gut microbiota of the Cobb-500 broiler (donors) has been previously characterized in several studies, including a distinct microbial community regarding improvement in feed-efficient and weight gain (Siegerstetter et al., 2018; Richards et al., 2019).
In the present study, the dominant phyla were Firmicutes and Proteobacteria followed by Bacteroidetes, Cyanobacteria, and Actinobacteria; however, the proportion of each phylum fluctuated and was affected by the FMT treatment. We herein observed that Firmicutes and the F/B ratio were increased in the FMT group, which was also consistent with an increase in body weight gain of chicks. We therefore speculate that this may be because Firmicutes and a high F/B ratio in the cecum microbiota allow a more efficient use of feed energy, since, previous studies have shown that the F/B ratio is associated with mammalian weight gain in different species, such as human (Karvonen et al., 2019), mice (Everard et al., 2014), and rabbit (Zeng et al., 2015). Additionally, the difference between the control group and the FMT group may be attributed to the latter having obtained abundant microbiota from the broiler chicks, which have a higher abundance of microbes in their digestive tract and can thus intake more nutrients.
Moreover, most α-diversity indices (Ace, Chao, Shannon, Simpson) were higher in the control group than in the FMT group, suggesting that the gut microbiota in the control group was more diverse than in the FMT group, which appeared to be more homogeneous based on the overabundance of Firmicutes and the F/B ratio. Additionally, in the current study, Lactobacillus, Lactococcus, and Bifidobacterium were more abundant in chicks received FMT. Lactobacillus is known to be involved in carbohydrate (Goh and Klaenhammer, 2014) and lipid metabolism (Falcinelli et al., 2015) and in immune system metabolism (Xin et al., 2014). Lactobacillus is the most important bacterial genus as it facilitates nutrient absorption, enhances host immunity, and prevents intestinal inflammatory responses (Wang et al., 2018). Collectively, Lactobacilli and Bifidobacteria are purported beneficial for gut physiology and body weight gain.
Gut microbiota composition is one of the most internal factors affecting physiological responses of the host through host–microbiome interaction, and a normal gut microbiota community allows gut permeability, nutrient digestibility, and anti-inflammatory reactions (Elokil et al., 2020c). Our results suggest that the shift in microbe distribution resulting from FMT successfully improved body weight gain. Therefore, the gut microbiota of the FMT group can digest complex and simple carbohydrates and produce more nutrients, such as volatile fatty acids, microbial proteins, and vitamins, in order for the host to gain weight in comparison with the control group.
Nowadays, CRISPR-Cas technologies, which were developed from the prokaryotic immune system, have recently made it possible for scientists to examine and modify organisms with an unparalleled ease and effectiveness (Vercoe et al., 2013; Ahmed et al., 2018). The delivery of CRISPR-Cas systems to bacterial populations is another option. These systems can be designed to target and destroy particular microbiome components (Ramachandran and Bikard, 2019). Together, these techniques offer fascinating chances to explore the intricate relationships that exist between the microbiome’s constituent parts and human bodies. They also open up fresh possibilities for the creation of medications that specifically target the microbiome (Stout et al., 2018). Most bacteria have CRISPR-Cas systems, which are present in 40% of them. In certain circumstances, it may be possible to use endogenous CRISPR-Cas systems, while in other cases, it may be possible to introduce designed CRISPR-Cas systems into the target bacteria (Hullahalli et al., 2017). The genomes of probiotic or microbiome-associated yeast, bacteria, and bacteriophages can all be altered using these methods. Additionally, based on their sequence, they can be employed to eradicate particular strains without affecting the remainder of the microbiome (Jiang et al., 2013).
Our findings revealed that FMT treatment promoted early body weight gain in Jing Hong chickens, speculating the physical growth will be completed early and will increase the length of laying time in the hens of FMT group than control group (Choi and Cho, 2016). Predictably, cocks in FMT group will have good physical and external morphometric traits such as comb height and color, wattle length, earlobe width, wingspan, girth circumference, drumstick length and testicular weight (Metzler-Zebeli et al., 2016). Since, the body weight- age relationship of laying chicks is directly related to the egg weight and number of egg produced for the laying chickens (Di Masso et al., 1998), and egg weight attains proper flock uniformity which is a key criterion for the success of egg-laying performance (Lacin et al., 2008). Which was achieved in the current study through FMT that regulated gut microbiome uniformity, and increased.
Conclusion
In conclusion, our findings showed that a distinct microbial community were colonized the gut due to FMT administration from the donor of highly feed-efficient broilers to laying chicks, which improved intestinal functions and antimicrobial pathogenicity compared with the control group. In addition, early life programming of slow growing chicks breed (Jing Hong) by FMT application was successful via increasing the cecal abundance of Firmicutes, Lactobacillus, Lactococcus, and Bifidobacterium, which may be related to increase the carbohydrate metabolism and health in laying chicks. FMT could be a potential strategy to improve animal growth performance.
Data availability statement
The data presented in the study are deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2021) in National Genomics Data Center (Nucleic Acids Res 2022), China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences repository, accession number GSA: CRA007995 that are publicly accessible at https://ngdc.cncb.ac.cn/gsa.
Ethics statement
The animal study was reviewed and approved by The experimental procedures used in this study met the guidelines of the Care and Use of Laboratory Animals of the Standing Committee of Hubei People’s Congress (No. 5) and approved by the Biological Studies Animal Care Committee of Hubei Province P.R. China, and the Ethics Committee of Huazhong Agricultural University, PR China, with approval number HZAUCH-19006. In this study, all efforts were made to minimize the animal suffering and was carried out in compliance with the ARRIVE guidelines.
Author contributions
AbE, H-ZL, and SL contributed to conception and design of the study and performed the experiment. WC, KA, ME, and KM organized the database. H-ZI, AhE, MM, and AA performed the statistical analysis. SL and WC were partial funding collaborators in this project. All authors equally wrote and revised the manuscript, contributed to the article, and approved the submitted version.
Funding
This work was supported by Natural Science Foundation of China (No.32072707); Science and Technology Major Project of Hubei Province (2020ABA016); the grants from the National Key R&D Program of China (SQ2021YFD1300002); The Special Fund for Scientific Innovation Strategy-Construction of High Level Academy of Agriculture Science (R2020PY-JX008 and 202106TD); Egyptian Ministry of Scientific Research, The Science, Technology & Innovation Funding Authority (39382). The funders had no role in the design, analysis, interpretation and writing of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1022783/full#supplementary-material
Footnotes
References
Abdelatty, A. M., Mandouh, M. I., Mohamed, S. A., Busato, S., Badr, O. A. M., Bionaz, M., et al. (2021). Azolla leaf meal at 5% of the diet improves growth performance, intestinal morphology and p70S6K1 activation, and affects cecal microbiota in broiler chicken. Animal 15:100362. doi: 10.1016/j.animal.2021.100362
Abouelezz, K., Abou-Hadied, M., Yuan, J., Elokil, A. A., Wang, G., Wang, S., et al. (2019). Nutritional impacts of dietary oregano and Enviva essential oils on the performance, gut microbiota and blood biochemicals of growing ducks. Animal 13, 2216–2222. doi: 10.1017/S1751731119000508
Ahmed, W., Hafeez, M. A., Ahmad, R., and Mahmood, S. (2018). CRISPR-Cas system in regulation of immunity and virulence of bacterial pathogens. Gene Reports 13, 151–157. doi: 10.1016/j.genrep.2018.10.004
Borody, T. J., and Khoruts, A. (2012). Fecal microbiota transplantation and emerging applications. Nat. Rev. Gastroenterol. Hepatol. 9, 88–96. doi: 10.1038/nrgastro.2011.244
Busato, S., Ford, H. R., Abdelatty, A. M., Estill, C. T., and Bionaz, M. (2022). Peroxisome proliferator activated receptor activation in precision-cut bovine liver slices reveals novel putative PPAR targets in Periparturient dairy cows. Front. Vet. Sci. 9:931264. doi: 10.3389/fvets.2022.931264
Cammarota, G., Ianiro, G., Bibbo, S., and Gasbarrini, A. (2014). Gut microbiota modulation: probiotics, antibiotics or fecal microbiota transplantation. Intern. Emerg. Med. 9, 365–373. doi: 10.1007/s11739-014-1069-4
Chao, A., and Shen, T. J. (2004). Nonparametric prediction in species sampling. J. Agric. Biol. Environ. Stat. 9, 253–269. doi: 10.1198/108571104X3262
Choi, H. H., and Cho, Y.-S. (2016). Fecal microbiota transplantation: current applications, effectiveness, and future perspectives. Clin. Endosc. 49, 257–265. doi: 10.5946/ce.2015.117
Cornejo-Pareja, I., Munoz-Garach, A., Clemente-Postigo, M., and Tinahones, F. J. (2019). Importance of gut microbiota in obesity. Eur. J. Clin. Nutr. 72, 26–37. doi: 10.1038/s41430-018-0306-8
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Deweerdt, S. (2014). Microbiome: A complicated relationship status. Nature 508, S61–S63. doi: 10.1038/508S61a
Di Masso, R. J., Dottavio, A. M., Canet, Z. E., and Font, M. T. (1998). Body weight and egg weight dynamics in layers. Poult. Sci. 77, 791–796. doi: 10.1093/ps/77.6.791
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Elokil, A. A., Abouelezz, K., Adetula, A. A., Ahmad, H. I., Mo, C., Sun, C., et al. (2020a). Investigation of the impact of gut microbiotas on fertility of stored sperm by types of hens. Poult. Sci. 99, 1174–1184. doi: 10.1016/j.psj.2019.10.048
Elokil, A. A., Abouelezz, K. F. M., Ahmad, H. I., Pan, Y., and Li, S. (2020b). Investigation of the impacts of antibiotic exposure on the diversity of the gut microbiota in chicks. Animals 10:896. doi: 10.3390/ani10050896
Elokil, A. A., Magdy, M., Melak, S., Ishfaq, H., Bhuiyan, A., Cui, L., et al. (2020c). Faecal microbiome sequences in relation to the egg-laying performance of hens using amplicon-based metagenomic association analysis. Animal 14, 706–715. doi: 10.1017/S1751731119002428
Everard, A., Lazarevic, V. N., Johansson, M., Ståhlman, M., Backhed, F., Delzenne, N. M., et al. (2014). Microbiome of prebiotic-treated mice reveals novel targets involved in host response during obesity. ISME J. 8, 2116–2130. doi: 10.1038/ismej.2014.45
Falcinelli, S., Picchietti, S., Rodiles, A., Cossignani, L., Merrifield, D. L., Taddei, A. R., et al. (2015). Lactobacillus rhamnosus lowers zebrafish lipid content by changing gut microbiota and host transcription of genes involved in lipid metabolism. Sci. Rep. 5, 1–11. doi: 10.1038/srep09336
Goh, Y. J., and Klaenhammer, T. R. (2014). Insights into glycogen metabolism in Lactobacillus acidophilus: impact on carbohydrate metabolism, stress tolerance and gut retention. Microb. Cell Factories 13:94. doi: 10.1186/s12934-014-0094-3
Grehan, M. J., Borody, T. J., Leis, S. M., Campbell, J., Mitchell, H., and Wettstein, A. (2010). Durable alteration of the colonic microbiota by the administration of donor fecal flora. J. Clin. Gastroenterol. 44, 551–561. doi: 10.1097/MCG.0b013e3181e5d06b
Hullahalli, K., Rodrigues, M., and Palmer, K. L. (2017). Exploiting CRISPR-Cas to manipulate enterococcus faecalis populations. elife 6:e26664. doi: 10.7554/eLife.26664
Jiang, W., Bikard, D., Cox, D., Zhang, F., and Marraffini, L. A. (2013). RNA-guided editing of bacterial genomes using CRISPR-Cas systems. Nat. Biotechnol. 31, 233–239. doi: 10.1038/nbt.2508
Karvonen, A. M., Sordillo, J. E., Gold, D. R., Bacharier, L. B., O’Connor, G. T., Zeiger, R. S., et al. (2019). Gut microbiota and overweight in 3-year old children. Int. J. Obes. 43, 713–723. doi: 10.1038/s41366-018-0290-z
Lacin, E., Yildiz, A., Esenbuga, N., and Macit, M. (2008). Effects of differences in the initial body weight of groups on laying performance and egg quality parameters of Lohmann laying hens. Czech J. Anim. Sci. 53, 466–471. doi: 10.17221/341-CJAS
Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lozupone, C. A., Hamady, M., Kelley, S. T., and Knight, R. (2007). Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 73, 1576–1585. doi: 10.1128/AEM.01996-06
Matsuoka, K., Mizuno, S., Hayashi, A., Hisamatsu, T., Naganuma, M., and Kanai, T. (2014). Fecal microbiota transplantation for gastrointestinal diseases. Keio J. Med. 63, 69–74. doi: 10.2302/kjm.2014-0006-RE
Meng, H., Zhang, Y., Zhao, L., Zhao, W., He, C., Honaker, C. F., et al. (2014). Body weight selection affects quantitative genetic correlated responses in gut microbiota. Plo S one 9:e89862. doi: 10.1371/journal.pone.0089862
Metzler-Zebeli, B. U., Lawlor, P. G., Magowan, E., and Zebeli, Q. (2016). Effect of freezing conditions on fecal bacterial composition in pigs. Animals 6:18. doi: 10.3390/ani6030018
Ramachandran, G., and Bikard, D. (2019). Editing the microbiome the CRISPR way. Philos. Trans. R. Soc. B 374:20180103. doi: 10.1098/rstb.2018.0103
Ramett, A. (2007). Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 62, 142–160. doi: 10.1111/j.1574-6941.2007.00375.x
R Development Core Team (2015). R: a language and environment for statistical computing. Vienna: R foundation for statistical computing. Available at: http://www.r-project.org/
Richards, P., Fothergill, J., Bernardeau, M., and Wigley, P. (2019). Development of the caecal microbiota in three broiler breeds. Frontiers in veterinary science 6:201. doi: 10.3389/fvets.2019.00201
Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Siegerstetter, S. C., Petri, R. M., Magowan, E., Lawlor, P. G., Zebeli, Q., O’Connell, N. E., et al. (2018). Fecal microbiota transplant from highly feed-efficient donors shows little effect on age-related changes in feed-efficiency-associated fecal microbiota from chickens. Appl. Environ. Microbiol. 84, e02330–e02317. doi: 10.1128/AEM.02330-17
Sonnenburg, J. L., and Bäckhed, F. (2016). Diet–microbiota interactions as moderators of human metabolism. Nature 535, 56–64. doi: 10.1038/nature18846
Stout, E. A., Sanozky-Dawes, R., Goh, Y. J., Crawley, A. B., Klaenhammer, T. R., and Barrangou, R. (2018). Deletion-based escape of CRISPR-Cas9 targeting in lactobacillus gasseri. Microbiology 164, 1098–1111. doi: 10.1099/mic.0.000689
Tawfik, A., Hasanan, K., Abdullah, M., Badr, O. A., Awad, H. M., Elsamadony, M. E.-D., et al. (2022). Graphene enhanced detoxification of wastewater rich 4-nitrophenol in multistage anaerobic reactor followed by baffled high-rate algal pond. J. Hazard. Mater. 424:127395. doi: 10.1016/j.jhazmat.2021.127395
Underwood, G., Andrews, D., Phung, T., and Edwards, L. E. (2021). Incubation, hatchery practice and the welfare of layer hens. Animal Production Science 61:867. doi: 10.1071/AN20391
Vadder de, F., and Mithieux, G. (2018). Gut-brain signaling in energy homeostasis: the unexpected role of microbiota-derived succinate. J. Endocrinol. 236, R105–R108. doi: 10.1530/JOE-17-0542
Vercoe, R. B., Chang, J. T., Dy, R. L., Taylor, C., Gristwood, T., Clulow, J. S., et al. (2013). Cytotoxic chromosomal targeting by CRISPR/Cas systems can reshape bacterial genomes and expel or remodel pathogenicity islands. PLoS Genet. 9:e1003454. doi: 10.1371/journal.pgen.1003454
Walugembe, M., Hsieh, J. C. F., Koszewski, N. J., Lamont, S. J., Persia, M. E., and Rothschild, M. F. (2015). Effects of dietary fiber on cecal short-chain fatty acid and cecal microbiota of broiler and laying-hen chicks. Poult. Sci. 94, 2351–2359. doi: 10.3382/ps/pev242
Wang, J., Ji, H., Wang, S., Liu, H., Zhang, W., Zhang, D., et al. (2018). Probiotic Lactobacillus plantarum promotes intestinal barrier function by strengthening the epithelium and modulating gut microbiota. Front. Microbiol. 9:1953. doi: 10.3389/fmicb.2018.01953
Xin, J., Zeng, D., Wang, H., Ni, X., Yi, D., Pan, K., et al. (2014). Preventing non-alcoholic fatty liver disease through Lactobacillus johnsonii BS15 by attenuating inflammation and mitochondrial injury and improving gut environment in obese mice. Appl. Microbiol. Biotechnol. 98, 6817–6829. doi: 10.1007/s00253-014-5752-1
Yan, C., Xiao, J., Li, Z., Liu, H., Zhao, X., Liu, J., et al. (2021). Exogenous fecal microbial transplantation alters fearfulness, intestinal morphology, and gut microbiota in broilers. Frontiers in Veterinary Science 8:987. doi: 10.3389/fvets.2021.706987
Zeng, B., Han, S., Wang, P., Wen, B., Jian, W., Guo, W., et al. (2015). The bacterial communities associated with fecal types and body weight of rex rabbits. Sci. Rep. 5, 1–8. doi: 10.1038/srep09342
Keywords: microbial transplantation, early life programming, 16S rRNA gene, body weight gain, laying chicken
Citation: Elokil AA, Chen W, Mahrose K, Elattrouny MM, Abouelezz KFM, Ahmad HI, Liu H-Z, Elolimy AA, Mandouh MI, Abdelatty AM and Li S (2022) Early life microbiota transplantation from highly feed-efficient broiler improved weight gain by reshaping the gut microbiota in laying chicken. Front. Microbiol. 13:1022783. doi: 10.3389/fmicb.2022.1022783
Edited by:
Xihong Zhou, Institute of Subtropical Agriculture (CAS), ChinaReviewed by:
Waqas Ahmed, The University of Tennessee, United StatesAqeel Ahmad, University of Florida, United States
Copyright © 2022 Elokil, Chen, Mahrose, Elattrouny, Abouelezz, Ahmad, Liu, Elolimy, Mandouh, Abdelatty and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Chen, cwei010230@163.com; Shijun Li, lishijun@mail.hzau.edu.cn
†ORCID: Abdelmotaleb A. Elokil https://orcid.org/0000-0001-7175-1408
Khaled F. M. Abouelezz https://orcid.org/0000-0002-6141-1095