Exploring the role of gut microbiota in host feeding behavior among breeds in swine

The interplay between the gut microbiota and feeding behavior has consequences for host metabolism and health. The present study aimed to explore gut microbiota overall influence on feeding behavior traits and to identify specific microbes associated with the traits in three commercial swine breeds at three growth stages. Feeding behavior measures were obtained from 651 pigs of three breeds (Duroc, Landrace, and Large White) from an average 73 to 163 days of age. Seven feeding behavior traits covered the information of feed intake, feeder occupation time, feeding rate, and the number of visits to the feeder. Rectal swabs were collected from each pig at 73 ± 3, 123 ± 4, and 158 ± 4 days of age. DNA was extracted and subjected to 16 S rRNA gene sequencing. Differences in feeding behavior traits among breeds during each period were found. The proportion of phenotypic variances of feeding behavior explained by the gut microbial composition was small to moderate (ranged from 0.09 to 0.31). A total of 21, 10, and 35 amplicon sequence variants were found to be significantly (q-value < 0.05) associated with feeding behavior traits for Duroc, Landrace, and Large White across the three sampling time points. The identified amplicon sequence variants were annotated to five phyla, with Firmicutes being the most abundant. Those amplicon sequence variants were assigned to 28 genera, mainly including Christensenellaceae_R-7_group, Ruminococcaceae_UCG-004, Dorea, Ruminococcaceae_UCG-014, and Marvinbryantia. This study demonstrated the importance of the gut microbial composition in interacting with the host feeding behavior and identified multiple archaea and bacteria associated with feeding behavior measures in pigs from either Duroc, Landrace, or Large White breeds at three growth stages. Our study provides insight into the interaction between gut microbiota and feeding behavior and highlights the genetic background and age effects in swine microbial studies.


Background
Feeding behavior is a collection of activities that reflect the feelings of hunger or satiety of animals at different physiological or developmental stages [1]. The feeding behavior of pigs can be evaluated by feed intake, time spent eating, number of meals per day, length of and interval between meals, and feeding rate [2,3].
Numerous studies have addressed the significant relationship between feeding behavior and growth performance [4], carcass quality [4], and feed efficiency [5,6] in swine. Feeding behavior conveys essential information for producers to optimize feed management and identify sick animals to achieve production goals.
The regulation of feeding behavior results from a vast network of interactions between the gastrointestinal (GI) tract and the brain. Miroschnikow et al. [7] found evidence of gut-derived chemical or neural signals functioning on the central nervous system to regulate various physiological activities, including feeding behavior. The diverse microorganisms colonizing the GI tract (hereafter referred to as gut microbiota) are extensively involved in nutrient absorption, energy metabolism, immune activities, and disease resistance through profound symbiotic interactions with the host in humans and animals [8,9]. In human and rodent models, studies suggest that the gut microbiota modifies brain function and modulates appetite and eating behavior [1,10,11]. Meanwhile, diet and feeding patterns have been found to induce rapid shifts in the gut microbial composition [12,13]. This intimate relationship between the gut microbiota and host feeding behavior relies on energy exchange. The microbiota relies on the host for nutrition to maintain its growth and population, while the host benefits from the metabolic activity of microbiota, partly digesting and assimilating the produced metabolites [1].
Host genetic background has been found to influence both feeding behavior [14] and gut microbiota composition [15] in pigs. Additionally, both feeding behavior [2] and gut microbial community [16][17][18] are dynamic over time. The gut microbiota matures with age and relatively stabilizes in adults, with the main changes occurring during the post-natal and peri-weaning periods, which are mainly associated with maternal, nutritional, and environmental factors [16,[19][20][21]. Longitudinal investigations in pigs have revealed several enterotypes of the gut microbiota and their relationships with growth performance at different stages during the growing-finishing phase [17,22]. These findings indicate that the interaction between the gut microbiota and host feeding behavior may depend on factors such as host genetics and age.
However, limited information on this interaction is available in different breeds of pigs during the growingfinishing period. Our previous study characterized differences in the gut microbial composition and investigated its associations with feed efficiency at three time points across the three breeds [23]. Based on that, this study aims (i) to compare feeding behavior traits among three commercial breeds; (ii) to investigate the role of the gut microbiota in explaining the variations of phenotypic traits; (iii) to identify the specific microbes that are associated with feeding behavior at different stages during the growing-finishing stage in swine.

Experimental animals and sample collection
Data used in this study were provided by Smithfield Premium Genetics (SPG; Rose Hill, NC, USA). Procedures for rectal swab collection were approved by the Institutional Animal Care and Use Committee, North Carolina State University, and National Pork Board. Microbial and feed efficiency data from this experiment are described in details in [23] and will therefore not be discussed in the present study. This study included 651 boars of either Duroc (DR; n =205), Landrace (LR; n = 226), or Large White (LW; n = 220) breed. Boars were the progeny of the mating of 28 sires and 124 dams for DR, 27 sires and 154 dams for LR, and 45 sires and 161 dams for LW. The number of siblings per family ranged from one to five for DR and from one to three for LR and LW. Most of the families (84% for DR, 90% for LR, and 96% for LW) had one or two siblings. The distributions of the family relationships were similar across the three breeds. The family relationship among individuals is provided in Additional file (1) Pigs weaned at the mean age of 24 days were kept in single-breed groups and housed in pens with an average of nine, eleven, and ten pigs per pen (12.4 m 2 ) for DR, LR, and LW, respectively. There were 64 such pens located across eight rooms (8 pens per room). Each room had all three breeds with headcounts as presented in Additional file (2) The experiment was conducted concurrently for three breeds in each room from May to December 2017. Each pen was installed with a single-space FIRE feeder (Osborne Industries Inc., Osborne, KS), which allows one pig to visit at a time. The feed consumption, feeder occupation time, bodyweight, and animal identifier were recorded whenever a pig visited the feeder. Pigs had 24-hour access to the feeder, and they were provided the same pelleted feed over the entire study. Details of the diet formulation are provided in Additional file (3) Feeding events were recorded on three breeds from 73 ± 3 to 163 ± 6 days of age. Rectal swabs were collected from each pig repeatedly at three time points (73 ± 3, 123 ± 4, and 158 ± 4 days of age).

DNA extraction and 16 S rRNA gene sequencing
The procedures to extract DNA (gDNA) from each rectal swab for Illumina library preparation are described in Lu et al. [22]. After DNA extraction, 100 µL of crude DNA was purified using a QIAquick 96 PCR purification kit (Qiagen, MD, USA). Purification procedures were carried out according to the manufacturer's protocol with a few minor modifications described in Lu et al. [22]. Purified DNA concentration was quantified using the Qubit 2.0 Fluorometer (Thermo Fisher Scientific Inc., MA, USA).
The V4 region of the 16 S rRNA gene from genomic DNA was amplified using the primers 515 F (5-GTG CCA GCMGCC GCG GTAA-3) and 806R (5-GGA CTA CHVGGG TWT CTAAT-3) to generate indexed libraries, as described in [22,24]. Equal amounts of amplicon libraries quantified using the Qubit dsDNA assay kit (Thermo Fisher Scientific Inc., MA, USA) were pooled together. The pooled amplicons were purified using the Agencourt AMPure XP beads (Beckman Coulter Inc., IN, USA), following the manufacturer's protocol. Purified pools were supplemented with 5-10% PhiX control DNA and were sequenced with Illumina MiSeq paired-end 2 × 250 bp + 13 bp index reactions using the 600v3 kit. Negative controls of microbial DNA-free water were included to confirm the absence of bacterial contamination during the Illumina library preparation. ATCC mock community standards (20-strain staggered mix, MSA-1003; ATCC, VA, USA) were used as positive controls to confirm the absence of lineage-specific biases in the upstream methods. Samples were sequenced at the DNA Sequencing Innovation Lab at the Center for Genome Sciences and Systems Biology at Washington University (St. Louis, MO, USA).

Data editing Phenotypic data
Data quality control was performed on the feeding events recorded by the FIRE feeder to identify and remove errors and outliers due to feeder malfunctions or animal-feeder interactions. The quality control categories and procedures on FIRE feeder data used in this study are described by Casey et al. [25]. After quality control, a total of 470,414 feeding records remained in the dataset. Corresponding to the three time points of rectal swab collection (T1: 73 ± 3 days of age; T2: 123 ± 4 days of age; T3: 158 ± 4 days of age), feeding records were split into three periods (P1: 73 ± 3 to 98 days of age; P2: 99 to 140 days of age; P3: 141 to 163 ± 6 days of age). The breakpoints at 98 and 140 days of age are the midpoints between T1 and T2 and T2 and T3 of rectal swab collection. For each period, seven feeding behavior measures, including the average daily feed intake (ADFI), average daily feeder occupation time (AOTD), average daily feeding rate (ADFR), average daily number of visits to the feeder (ANVD), average feed intake per visit (AFIV), average feeder occupation time per visit (AOTV), and average feeding rate per visit (AFRV), were calculated for each pig similar to what described by Lu et al. [6].

Sequencing data
The raw sequencing reads obtained from the Illumina platform were converted into fastq format files using MiSeq Reporter. Paired sequences were merged into a single sequencing file using FLASH (v1.2.11), with the minimum overlap length was set to 100 bp and the maximum length was set to 250 bp [26]. Sequences were trimmed to remove the primers and were demultiplexed using the Quantitative Insights Into Microbial Ecology (QIIME2) [27]. The amplicon sequence variant (ASV) feature table was constructed and denoised using the Divisive Amplicon Denoising Algorithms 2 (DADA2) with default settings and no truncation [28]. Taxonomic information was annotated using the Ribosomal Database Project (RDP) Classifier (v2.4) based on the SILVA reference database (v132) [29] with a confidence cutoff of 0.8. The ASVs were filtered by a minimum total count of 1000 across all the samples and were removed with a prevalence rate less than 0.05 at each sampling time point. A total of 750, 724, and 824 ASVs remained after the quality control for time point 1 (T1; 73 days of age), time point 2 (T2; 123 days of age), and time point 3 (T3; 158 days of age), respectively.

Effects of breed on feeding behavior traits
Before analysis, phenotypes normality assumptions were verified, and traits ADFR during P1, AFRV during P1 and P2, and ANVD, AFIV, and AOTV during all three periods were log-transformed in the R environment (v4.0.2) [30]. To explore the relationship between feeding behavior traits, pairwise Spearman's rank correlation coefficient (r s ) was estimated for each breed group during each period due to the normality assumption's violation in some traits [31]. Statistical significance for the correlations was considered present at α ≤ 0.05.
To estimate the effect of breed and to account for the environmental impact on feeding behavior traits, a linear mixed model was constructed using the PROC GLIM-MIX in SAS (v9.4, SAS Institute, Carry, NC, USA). The model was conducted on each trait during each period separately, with Breed and Room fitted as the fixed effects and Sire and Pen fitted as the random effects. The least-squares means for fixed effects were obtained using the LSMEANS statement of the PROC GLIMMIX, and they were compared among breeds during each period. A P-value less than or equal to 0.05 was considered significantly different.

Microbial parameter of host feeding behavior traits
We applied the following two models to investigate the proportion of variance in host phenotypes due to differences in the gut microbial composition among individuals for each breed. Microbial measurements at T1 (start of P1), T2 (midpoint of P2), and T3 (end of P3) were used in the following models to estimate the proportion of variation in feeding behavior during the respective period. First, we defined a baseline model as: where y ijkl is the feeding behavior trait, µ is the overall intercept, Room i is the fixed effect of i th room, sire j is the random effect of j th sire, pen k is the random effect of k th pen, and e ijkl is the residual. The effects of sire, pen, and residual were assumed to be distributed as N(0, Iσ 2 p ), N(0, Iσ 2 p ), and N(0, Iσ 2 p ), respectively; where I was an identity matrix, and σ 2 2 , σ 2 p , and σ 2 e were the variances for the respective effects.
To estimate the variance in feeding behavior that is contributed by the gut microbiota of pigs at each growth stage, we then defined a univariate mixed model, in which microbial information was added to the baseline model as: where the model parameters are the same described in the baseline model, with the addition of m l , which contains the random microbial effect for animal l at the given time point with m ~ N (0, Mσ 2 m ), where M is the microbial relationship matrix and σ 2 m is the microbial variance. The matrix M was created for each time point separately as follows: We first constructed a matrix V based on the ASVs feature table with a dimension n × N, where n is the number of animals and N is the number of ASVs for each time point. Each element V ij represented relative abundance of ASV j in animal i for a given time point. A value of 0.01 was added to each V ij to guarantee positive definiteness, as described by Camarinha-Silva et al. [32]. Then, the matrix V was used to calculate the elements of matrix X with the same dimension as: Thus, the matrix X contained log-transformed, centered, and scaled relative abundance of ASV. The microbial relationship matrix M was then computed as M = 1 N XX T , representing the covariance between individuals based on their microbial composition resemblance.
We fitted the mixed models using a package BGLR [33] in R environment (v4.0.2) [30]. A Markov chain Monte Carlo (MCMC) algorithm was used to perform the analyses with 150,000 iterations, 50,000 iterations discarded as burn-in, and 10 iterations set as the thinning interval. Convergence of the models was visually checked by looking at trace plots of each parameter posterior distribution. The total phenotypic variance was obtained as the sum of σ 2 s , σ 2 p , σ 2 m , and σ 2 e . The proportion of phenotypic variance explained by the microbiota was calculated as the σ 2 m over the phenotypic variance.
This proportion was defined with the term "microbiability" (m 2 ) by Difford et al. [34,35] and represented the overall interaction of the microbiota with the host phenotype. ANOVA was performed using the PROC GLM in SAS (v9.4, SAS Institute, Carry, NC, USA) to compare the microbiability estimates among feeding behavior traits, breeds, and time points, as well as their interactions. The least-squares means for fixed effects and contrasts were obtained with the LSMEANS statement of the PROC GLM. A P-value less than or equal to 0.05 was considered significant.

Association between ASVs and feeding behavior traits
To investigate the association between individual ASVs and feeding behavior traits, the following linear mixed model was constructed using the PROC GLIMMIX in SAS (v9.4, SAS Institute, Cary, NC, USA): where y ijklmn is the feeding behavior trait at each time point, µ is the overall intercept, Breed i is the fixed effect of i th breed, Room j is the fixed effect of j th room, asv k is the covariate of the k th ASV, asv * Breed (ik) is the fixed effect of the (ik) th ASV and breed interaction, Sire l is the random effect of l th sire within breed, Pen m is the random effect of m th pen within room, and e ijklmn is the residual. The marginal effects of individual ASVs were evaluated, and the orthogonal contrasts were performed among breed groups at each time point. The Benjamini-Hochberg method was used to adjust the P-values to control the false discovery rate (FDR). An FDR of 0.05 was set as the cutoff threshold to declare the significance of the association between ASV and trait.

Breed feeding behavior characteristics at three time points
Phenotypic correlations for individual breed groups by period are depicted in Fig. 1. Generally, correlations between feeding behavior traits were similar across breeds during each period. Strong correlations were found between AFIV and AOTV (r > 0.75, P < 0.001), AFRV and ADFR (r > 0.85, P < 0.001), ANVD and AFIV (r < -0.80, P < 0.001), as well as ANVD and AOTV (r < -0.70, P < 0.001). As expected, the feeding rate (ADFR and AFRV) had negative relationships with the feeder occupation time (AOTD and AOTV).
The phenotypic differences of average daily feed intake (ADFI) among breeds have been reported in our previous study [23]. Thus, we excluded it and reported the results of the remaining six feeding behavior traits in the present study. Least squares means and 95% confidence intervals of feeding behavior traits for the effect of breed during each period are presented in Fig. 2. Significant differences between breeds were found for every feeding behavior trait during each period, except for feeding rate traits (ADFR and AFRV) during P3. The LW pigs spent more time y ijklmn = + Breed i + Room j + asv k + asv * Breed (ik) + sire l + pen m + e ijklmn in the feeder daily (AOTD) than the other two breeds across periods. The DR pigs visited the feeder less frequently (ANVD) but consumed more feed (AFIV) and spent more time visiting the feeder (AOTV) compared to the other breeds, especially during P2 and P3. The LW pigs had smaller feeding rates (ADFR and AFRV) than DR and LR pigs during P1 and P2, with no difference noticed during P3 among breeds. Differences

Microbiability of host feeding behavior traits
Relative abundance of the gut microbial taxa of DR, LR, and LW pigs at T1, T2 and T3 is described in our previous study [23]. The relative proportion of the phenotypic variation in host feeding behavior traits attributed to gut microbial composition, defined as "microbiability" (m 2 ), is depicted by breed and time point in Fig. 3. The m 2 estimates ranged from 0.10 (SE = 0.04) to 0.28 (SE = 0.11) for DR pigs, 0.09 (SE = 0.04) to 0.23 (SE = 0.09) for LR pigs, and 0.10 (SE = 0.04) to 0.31 (SE = 0.09) for LW pigs across the three time points. Results of differences in m 2 estimates across breeds, sampling time points, or feeding behavior traits are summarized in Additional file 6. The m 2 estimates were similar between DR and LR (P = 0.42), as well as between LR and LW (P = 0.32), but they were greater for DR than LW (P = 0.02). When the m 2 estimates were compared across the three time points, they were similar between T1 and T2 (P = 0.34) and between T2 and T3 (P = 0.28), while they were significantly higher at T3 than T1 (P = 0.01). Furthermore, compared to other feeding behavior traits, the gut microbial composition explained a greater proportion of the total variation in ADFI. There were no significant interaction effects between breed, time point, and feeding behavior trait. Estimates of variance components for each trait obtained from two models are depicted in Additional file 7. Compared to the baseline model, the inclusion of microbial information in the model helped reduce the proportion of variance previously included in the residual component for most feeding behavior traits.

Specific microbes associated with feeding behavior traits
A false discovery rate of 0.05 was considered the threshold to identify ASVs significantly associated with feeding behavior traits. The relationship between the gut microbiota and ADFI was reported in our previous study [23]. Thus, we only reported the results of the other feeding behavior traits in this section. Within each breed group, we identified 4, 6, and 11 ASVs for DR pigs, 2, 0, and 8 ASVs for LR pigs, and 3, 5, and 27 ASVs for LW pigs at T1, T2, and T3 associated with feeding behavior traits, respectively ( Table 1). The majority of the identified ASVs were associated with AOTD for all breed groups at T3. However, at the threshold of FDR  Venn diagrams of the number of common and unique identified ASVs among breed groups and time points are presented in Fig. 5. Identified ASVs were breed-specific at each time point with a minor overlapping among breed groups. There was only one ASV in common between LR and LW pigs at T1, and one ASV was shared between DR and LW pigs at T2 (Fig. 5 A). Three ASVs were shared by LW pigs with either DR or LR pigs, and no overlapping was found between DR and LR pigs at T3 (Figure 5 A). Similarly, when comparing identified ASVs within the breed group across time points, only two ASVs were found the same between T2 and T3 in LW pigs, and no shared ASV was found across time points in DR and LR pigs (Fig. 5B). The scaled coefficient estimates of identified ASVs by breed group at three time points are depicted in Fig. 6. A total of 55 unique ASVs were either positively or negatively associated with feeding behavior traits. Nine ASVs were positively associated with AFRV (n = 5), AOTD (n = 3), and AFIV (n = 1) in three breed groups at T1, and most of them were annotated to Lachnospiraceae family belonging to Firmicutes phylum (Fig. 6 A). Five ASVs negatively associated with AFIV (n = 1), AFRV (n = 1), and AOTD (n = 2) for DR pigs and AOTD (n = 1) for LW pigs at T2 were annotated to Methanobrevibacter, Weissella, and Ruminococcaceae_UCG-004 at the genus level. Besides, six ASVs positively associated with AOTD for DR and LW pigs were annotated to Methanobrevibacter (Euryarchaeota phylum), Dorea, Marvinbryantia, and Ruminococcaceae (Firmicutes phylum) at the genus level ( Fig. 6B). At T3, 29 ASVs with positive associations across breeds were annotated to 16 genus levels and mainly classified to either Blautia, Dorea, Ruminococcus_2, and Subdol. Most of the 17 ASVs with negative associations in three breeds were annotated to Christensenellaceae_R-7_ group, Family_XIII_AD3011_group, and Ruminococ-caceae_UCG-004 genus in Firmicutes phylum.
Among the 55 ASVs, eight were associated with feeding behavior traits across several breeds or time points (Fig. 6). For example, the ASV10 assigned to Ruminococ-caceae_UCG-004 genus was adversely associated with AOTD in DR and LR at T2 and in LR and LW at T3. Similarly, the ASV43 assigned to Christensenellaceae_R-7_ group was negatively associated with AOTD in DR and LW at T3. Additionally, the ASV15 belonging to Marvinbryantia genus and the ASV16 belonging to Dorea genus were positively associated with AOTD in LR and LW at T3, and in DR at T2 and LW at T3, respectively. Furthermore, several ASVs were associated with more than one feeding behavior trait; ASV6 (unassigned genus) with AFIV in LR and AFRV in LW at T1; ASV9 (Methanobrevibacter genus) with AFRV and AOTD in DR at T2; ASV42 (unassigned genus) with ANVD and AFIV in DR at T3; and ASV50 (Family_XIII_AD3011_group genus) with ANVD and AOTD in LW at T3. The contrasts of the coefficient estimate for the most significant ASV (least FDR) in each trait-breed-time point category are depicted in Additional file 8.

Discussion
Phenotypic correlations between feeding behavior traits in this study are similar to the results previously reported by Fernández et al. [14]. In agreement with Rauw et al. [36] and Carcò et al. [37], we observed an increase in both the feeding rate and feed consumption of pigs with age, most likely due to their increased physical capacity. As expected, we observed different feeding patterns between breeds. We found that Large White pigs had a lower feeding rate than Duroc and Landrace pigs during the first two periods. In contrast, Fernández et al. [14] found that Duroc had lowest feeding rate of the three breeds, defining them as "slow eaters". Variations in feeding behaviors among rooms or pens for growing pigs were noticed in this study. It is expected that with only one Feeder equipped per pen, social ranks in grouphoused pigs may lead to competition at feeding and result in a discrepancy in feeding behavior across pens [38].
First proposed by Difford et al. [34], estimation of microbiability has recently become a popular approach to quantify the proportion of variation in host phenotypes explained by the microbial composition in agricultural animal species. Difford et al. [34] reported that the ruminal microbiota composition was responsible for 15% of the total phenotypic variation in CH 4 emissions of dairy cows. Wen et al. [39] indicated that the duodenal and cecal microbiota explained 24% and 21% of fat deposition variation in chicken, respectively. Moreover, Vollmar et at. [40] reported that the variations in the gut microbiota were associated with 9%, 18%, and 27% of variations in the feed intake, daily gain, feed per gain ratio in Japanese Quail, respectively. Our microbiability estimates ranged from 0.09 to 0.31, which suggested that the gut microbiota was associated with a small to moderate amount of variation in host feeding behavior among animals. In this study, breed effects on microbiability were discovered. Differences in the microbiability estimates between Duroc and Large White pigs suggest that the gut microbial composition may differ in predicting the host phenotype in populations with diverse genetic backgrounds. Furthermore, our findings show that the gut microbial composition of pigs at the finishing stage (T3) is more informative in explaining the variations in feeding behavior than at the start of the growth period (T1). The differences in the microbiability estimates across time points could be explained by several factors. For example, if all microbial samples had been collected at a uniform time point (i.e., the midpoint) within each feeding behavior measurement period, we might have seen fewer variances in the microbiability. In addition, transitions in the diet and environment prior to our first sampling may have contributed to differences in the microbiability between growth stages. We also observed variations in the microbiability estimates across feeding behavior traits, with the highest values for average daily feed intake, ranging from 0.17 to 0.31, which are similar to the estimate of 0.16 for feed intake of Pietrain sows at the slaughter age estimated by Camarinha-Silva et al. [32]. To the best of our knowledge, our study is the first one to estimate the microbiability for feeder occupation time, feeding rate, and number of visits to the feeder in multiple breeds at three time points during the growth in swine. These findings suggest that host genetics and age may have an impact on microbiability for different traits in pigs, but more research is needed to confirm the findings in different populations.
Previous studies have suggested that the gut microbial composition varies among individuals due to breed and age in pigs [15,16,23]. In this study, the identified ASVs associated with feeding behavior were also breed-specific and stage-specific during the growing to finishing period, with a small part of ASVs shared among breeds or sampling time points. At the last sampling time point (158 days of age), the greatest number of ASVs were found associated with feeding behavior across breeds compared to the previous two sampling time points (73 and 123 days of age). This might be due to the relative higher alpha diversity measured by the Shannon Index and greater number of clusters in the gut microbiota at 158 days compared to early time points in all three breeds, as described in our previous study [23].
Among the 40 ASVs positively associated with feeding behavior, more than half of the ones associated with feed intake, feeding rate, and feeder occupation time were assigned to Lachnospiraceae family in Firmicutes phylum. Similarly, Cox et al. [41] identified four ASVs in Lachnospiraceae family that exhibited greater abundances in people with good appetite than people with poor appetite by comparing their gut microbial composition. At the genus level, two ASVs robust in positive associations with the daily feeder occupation time in multiple breeds and time points were annotated to Marvinbryantia and Dorea. The remaining identified ASVs positively associated with the daily feeder occupation time were mainly annotated to Blautia and Rumi-nococcaceae_UCG-014. Interestingly, these microbes have been found correlated to complex host behaviors, such as stressor-induced behavior and depressive-like behavior in human and animal models [42][43][44][45]. Shifts in the abundance of Blautia [43], Dorea [42], and Marvinbryantia [45] were observed in stressed and depressed mice. Similarly, reduction in the abundance of Rumino-coccaceae_UCG-014 was linked to the development of anxiety-like behaviors in humans [44]. The genus Prevotella positively associated with the daily feeder occupation time of Duroc pigs in our study were also found positively correlated with the appetite of Duroc pigs in a previous study conducted by Yang et al. [46]. One ASV belonging to Lactobacillus genus was found to associate with the daily feeding rate in Duroc pigs at 73 days of age in this study. Lactobacillus strains, commonly used as probiotics promoting overall intestinal health, have been found to influence the host eating behavior through various pathways [47,48].
The present study also identified 17 ASVs that had negative associations with host feeding behavior traits. At the family level, the majority of these ASVs were classified into Ruminococcaceae, Christensenellaceae, or Family_XIII. Numerous studies have indicated that the intestinal abundance of Christensenellaceae is negatively related to host body mass index and fat mass in human studies [49][50][51]. Similarly, a significant number of microbes belonging to Ruminococcaceae family colonized in the gut have been associated with the lower weight gain and lean phenotype in humans [52]. In our study, several microbes from the Ruminococcaceae and Christensenellaceae families were negatively associated with daily feeder occupation time or the number of visits to the feeder. It is possible that the reduced feeding behavior might be linked to the low fatness deposit and weight gain through specific gut taxa, but further studies are needed in this sense. Moreover, both positive and negative associations were identified between microbes from Ruminococcaceae family and feeding behavior measures, including daily feeder occupation time, feed intake per visit to the feeder, and feeding rate across breeds and time points in this study. Similarly, Yang et al. [46] reported one positive and two negative associations between Ruminococcaceae family members and feed intake in pigs. These findings suggest that Ruminococcaceae family may contain a diverse range of members in the gut, each with distinct functions. Further investigations on the specific individuals belonging to Ruminococcaceae family are needed to clarify this aspect. At the genus level, three ASVs involved in multiple negative associations with the daily feeder occupation time and number of visits to the feeder in several breeds and time points were annotated to Ruminococcaceae_UCG-004, Christensenellaceae_R-7_ group, and Family_XIII_AD3011_group. Among these microbes, Christensenellaceae_R-7_group [53], Rumi-nococcaceae_UCG-004 [54] are known to produce short-chain fatty acids, which are essential metabolites in regulating host energy homeostasis and appetite through gut-brain communication in the gut [55]. Accumulating research has indicated that elevated SCFAs production in the gut could suppress the host appetite [55][56][57], thus modifying feeding behavior. We identified an ASV belonging to Methanobrevibacter genus in a negative association with the average feeding rate per visit and also in a positive association with the average daily feeder occupation time in Duroc pigs at 123 days in the present study. Interestingly, several studies have shown that people with anorexia nervosa have higher levels of Methanobrevibacter in their gut [58,59]. The microbes found in this study could serve as candidates for further clarifications on the relationship between the gut microbiota and host feeding behavior. However, differences in the environment [60], experimental design [61], the choice of sequencing platform [62], and 16 S target regions [63] limit the comparisons of results from various microbiome studies. Thus, more studies are still needed to validate the results in different populations with various conditions and further establish the causality between specific microbes and feeding behavior in swine.

Conclusions
In this study, we observed different feeding patterns recorded by the electronic feeder among Duroc, Landrace, and Large White pigs at three time points during the growing-finishing period. By modeling the microbial relationship among animal individuals and estimating the microbiability, we quantified a low to moderate amount of variation in host feeding behavior associated with the gut microbial composition in each breed group across the three time points. Moreover, we identified a total of 55 unique microbes differently associated with feed intake, feeder occupation time, number of visits to the feeder, and feeding rate in three breed groups, which might provide information for future studies aiming at modulating feeding disorders and improving production in growing pigs. Our study provides insights to better understand the importance of the gut microbiota in interacting with host feeding behavior and highlights the influences of host genetic background and age in porcine microbiome studies.
Abbreviations ASV: Amplicon sequence variant; DR: Duroc; LR: Landrace; LW: Large White; ADFI: Average daily feed intake; AOTD: Average daily feeder occupation time; ADFR: Average daily feeding rate; ANVD: Average daily number of visits to the feeder; AFIV: Average feed intake per visit; AOTV: Average feeder occupation time per visit; AFRV: Average feeding rate per visit.