Gut Microbiological Disorders Reduce Semen Utilization Rate in Duroc Boars

Although rising evidence suggests that the gut microbiota is closely related to host health, the effects of gut microbiota on male fertility are still rarely explored. This study was to investigate the gut microbiota composition and function, fecal short-chain fatty acids (SCFA), intestinal permeability, and systemic inflammatory status of Duroc boar with high (H group, 100%) and low (L group, <80%) semen utilization rate. Fecal samples, analyzed by 16S ribosomal RNA gene sequencing, displayed taxonomic and functional changes between boars with high and low semen utilization rates. For the gut microbiota composition of the boars, four genera were different between the two groups. The [Ruminococcus] and Sphingobium were enriched in L group boars, then negatively correlated with the semen utilization rate. While RFN20 and Paludibacter were enhanced in the H group, only RFN20 showed a significantly positive correlation with the semen utilization rate of boars. In addition, changes in the metabolic function of the gut microbiota of the two groups were found, including altered branched-chain fatty acid (BCFA) production. Significant increases in plasma endotoxin, zonulin, diamine oxidase, and lipocalin-2 levels were observed in boars with low semen utilization, and also, a similar trend in IL-6 and TNF-α was found. However, the concentration of IL-10 in plasma of boars with high semen utilization rate showed an increasing tendency. These results indicated increased intestinal permeability and systemic inflammation in boars with low semen utilization. Data showed that the composition and functions of gut microbiota varied between boars with high or low semen utilization rates, while the semen utilization rate is notably correlated with the gut microbiota composition, intestinal permeability, and inflammatory status of the boar.

Introduction

Male fertility is of paramount importance for breeding animal herds when artificial insemination is involved (Waberski et al., 2008). Therefore, monitoring and analyzing the quality of boar semen has vast economic meaning for the pig industry (Smital, 2009). Many factors including varicocele, accessory gland infection, immunological factors, and endocrine causes can adversely affect sperm quality (Omu, 2013). Being one of the largest immune, endocrine, and neural organs in the body, the gut may affect male fertility, which has been underestimated. Researches showed that the gut microbiome influences host endocrine and immune functions (Smith et al., 2013; Rastelli et al., 2019). Moreover, emerging data have demonstrated that an aberrant gut microbiota composition is associated with several diseases, including metabolic disorders and inflammatory bowel disorder (Koh et al., 2016). However, the exact effects of gut microbiota on reproduction have not been elucidated.

Healthy gut microbiota plays many beneficial roles in the host gut, such as resistance to colonization of pathogens, immune regulation, and assisted digestion. Simultaneous production of bacterial metabolites, such as short-chain fatty acids (SCFA) and lactic acid, affects the health of the host (Sekirov et al., 2010; Flint et al., 2012), while there is new evidence that dysregulation of gut microbiota can lead to increased intestinal permeability and may act as a part in the development of chronic low-grade inflammation in the host (Cani et al., 2007; Yu et al., 2016). This chronic inflammation may be related to the gut microbiota components like lipopolysaccharide (LPS), which is a chronic innate immune stimulator in pigs (Webel et al., 1997; Gabler et al., 2008; Weber and Kerr, 2008). Once LPS enters the blood circulation through the gut barrier, it may bind to reproductive tract cells or indirectly stimulate immune cells, thereby affecting steroidogenesis and sperm quality (Hedger, 2011).

In this study, we evaluated the compositional, functional, and metabolic differences in the gut microbiota between the Duroc boars with high (n = 29) or low (n = 11) semen utilization rates. We also attempted to investigate the intestinal permeability and inflammatory status of boars of different semen utilization rates. The discoveries will provide new insights into increasing semen quality of the boar during high-intensity production and establish a foundation for further research on the impact of gut microbiota on male reproduction.

Materials and Methods

Animals

A total of 40 Duroc boars aged from12 to 54 months were selected in an artificial insemination center of Yangxiang Joint Stock Company (Guangxi, China). Boar feeding conditions were as our previous study (Wu et al., 2019a). Semen samples were collected by gloved-hand techniques. After collecting, four semen parameters were assessed, including semen volume, sperm concentration, sperm motility, and abnormal sperm rate, according to the method we described earlier (Wu et al., 2019b). According to the cut-off value recommended by previous studies, boar semen with sperm motility below 70% or abnormal rate above 20% is considered unusable for artificial insemination (Rodriguez et al., 2013; Wu et al., 2019a). Since the quality of the semen collected in a single shot is difficult to represent the real situation of the boar, we introduced a new concept called semen utilization rate, which was calculated based on the semen quality of each ejaculation within 3 months (the number of times that semen can be used/the total number of times semen collected). In the experiment, 40 boars were divided separately into a high semen utilization rate group and a low semen utilization rate group. The 11 boars with a semen utilization rate <80% were defined as “low semen utilization rate group” (L group), and those with a 100% semen utilization rate were defined as “high semen utilization rate group” (H group; Table 1). All animal research procedures were conducted in accordance with animal research guidelines issued by the Institutional Animal Care and Use Committee of Huazhong Agricultural University (permit number HZAUSW-2018-014).

TABLE 1

Table 1. Semen quality parameters of Duroc boars in different semen utilization rate groups (n = 40).

Sample Collections

Blood samples were harvested by venipuncture from the hindlimb vein of boars in ejaculations. The blood sample was then centrifuged at 3,000 × g for 10 min at 4°C to obtain a plasma sample and stored at −80°C until analysis. The boar’s rectum was massaged to stimulate boar defecation, and then fresh feces were collected and placed in a sterile 5-ml cryopreservation tube. Stool samples were put on the ice first and then transfer it to −80°C for subsequent analysis.

Detection of Intestinal Permeability and Systemic Inflammation Markers in Plasma

Zonulin is a medium well known to modulate gut permeability by modifying intracellular tight junctions. Increased circulating zonulin is considered a potential sign of increased intestinal permeability (Maria Moreno-Navarrete et al., 2012). Simultaneously, endotoxin and diamine oxidase in plasma are also vital indicators for intestinal epithelial permeability and integrity (Luk et al., 1983). The plasma zonulin, endotoxin, diamine oxidase, interleukin-6 (IL-6), IL-10, tumor necrosis factor (TNF)-α, and total lipocalin-2 concentrations were determined by using porcine enzyme-linked immunosorbent assay kits (mlbio, Shanghai, China), according to the manufacturer’s instructions.

Examination of Fecal Short-Chain Fatty Acids and Lactic Acid

The SCFA (acetate, propionate, butyrate, and valerate) concentrations and branched-chain fatty (BCFA; isobutyrate and isovalerate) concentrations in feces of the boars were analyzed by a gas chromatographic method, as described by our previous study (Cheng et al., 2019). Fecal lactic acid concentrations were determined by using porcine enzyme-linked immunosorbent assay kits (mlbio, Shanghai, China), according to the manufacturer’s instructions.

16S rRNA Gene Sequencing and Data Processing

Total genomic DNA was extracted from each fecal sample using QIAamp Fast DNA stool Minikit (Qiagen, Germany), according to the manufacturer’s instructions. The forward primer 341F (5'-ACT CCT ACG GGA GGC AGC AG-3') and the reverse primer 806R (5'-GGA CTA CHV GGG TWT CTA AT-3') were used for amplification of the V3-V4 hypervariable region of a 16S rRNA gene. Thirty cycles of PCR were performed at an annealing temperature of 56°C. PCR amplicons were purified and quantified. The PCR products were used for the construction of the libraries and then paired-end sequenced (2 × 250) on an MiSeq platform (Illumina, United States) at the Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China). All sequencing data were submitted to the NCBI database with accession number PRJNA645685. Sequencing data were processed using QIIME (version 1.8.0). In a few words, the original sequencing reads were sorted to obtain valid sequences. After filtering low-quality sequences, high-quality sequences were clustered into operational taxonomic units (OTUs) with a 97% sequence identity by UCLUST (Edgar, 2010). Then the most abundant sequence in each OTU was selected as the representative sequence. By comparing the representative sequences of OTUs with the template sequence of the Greengenes database (Desantis et al., 2006), the taxonomic information corresponding to each OTU was obtained. To minimize the difference in sequencing depth between samples, the average analysis of 100 evenly resampled OTU subsets under the 90% of the minimum sequencing depth was performed to generate an average, rounded dilution OTU table. Alpha-diversity values of each sample were assessed based on the observed OTUs, Chao1, and Shannon index. Beta-diversity measures dependent on weighted-UniFrac distance were calculated using mothur. LEfSe was conducted to identify bacterial taxa differentially represented between different groups at the genus or higher taxonomy level (Segata et al., 2011). PICRUSt analysis was used to predict the metagenome function of the microbiota (Langille et al., 2013).

Statistical Analyses

The normal distribution of data was confirmed by the Shapiro-Wilk’s test before analysis. Variations in normally distributed semen traits, SCFAs, lactic acid, and markers associated with gut permeability and inflammation were examined by the Student-t test (Statistical Analysis System 9.4; Cary, NC, United States). Significance was reported when the p-value is less than 0.05. The dissimilarities of bacterial relative abundance between two groups and alpha diversity were evaluated using the Mann-Whitney U test. The predicted pathways from PICRUSt were established by using the Welch’s t-test in the STAMP software package (version 2.1.3). Correlations were investigated by using Spearman’s correlation in R 3.6.1 (the R Foundation) with the R Studio 1.2.1335 package and ggplot2 for the heat map. According to the false-discovery rate (FDR) procedure, the relevant results were corrected by FDR analysis, with q of <0.05.

Results

Semen Parameters of Duroc Boars

The semen parameters of the two group boars are given in Table 1. Though there is no significant difference in sperm concentration and total sperm number between the two groups, the low utilization rate group boars showed significantly lower semen volume (p < 0.05), functional sperm number (p < 0.01), sperm motility (p < 0.001), semen utilization rate (p < 0.001) and higher abnormal sperm rate (p < 0.001) compared with the high utilization rate group.

Composition of Gut Microbiota in the Boars

In order to study the composition of gut microbiota in different semen quality boars, 16S ribosomal RNA gene sequencing was used. A total of 1,407,378 high-quality sequencing reads were obtained from 40 samples, ranging from 28,112 to 46,845. Based on 97% species similarity, 50,406 OTUs were obtained from samples of the boars, which involved 7,265 different OTU classifications. The alpha and beta diversities were calculated. Alpha-diversity measures (observed OTUs, Chao1, and Shannon index) of the fecal bacterial community presented no significant differences between the two groups (Figures 1A–C). For beta diversity, the nonmetric multidimensional scaling (NMDS), based on weighted-UniFrac distance, exposed that the gut microbiota showed negligible segregation in the L and H group boars (Figure 1D). Although there is no clear visual separation of groups within the study boars in Figure 1D, the intergroup variation from H to L was higher than the interindividual variation of boars in group H (Figure 1E). The PERMANOVA analysis indicated small but significant differences in the gut microbiota between the L group and the H group (weighted-UniFrac distance, R² = 0.083, p = 0.005). Furthermore, ANOSIM analysis based on weighted-UniFrac distance also confirmed that the microbial composition between the two groups was significantly different (R = 0.170, p = 0.036). As shown in Figure 2, the L and H group boars have different gut microflora composition at the phylum (Figure 2A) and the genus levels (Figure 2B). To identify the key genera in the gut microbiota of the L and H group boars, LDA was conducted to evaluate the effect size of each community that could effectively discriminate different groups (Figure 2C). LDA effect size revealed that in the genus level, two known genera (Sphingobium and [Ruminococcus]) were significantly enriched in the low utilization rate group boars, whereas two genera (Paludibacter and RFN20) were significantly enriched in the high utilization rate group boars (p < 0.05, Wilcoxon rank-sum test; log LDA > 2) (Figure 2C). The relative abundance of these four genera in each sample is shown in Supplementary Figure 1.

FIGURE 1

Figure 1. Alpha and beta diversity of gut microbiota. (A–C) Boxplots of observed operational taxonomic units (OTUs) (A), Chao1 index (B), and Shannon index (C) for the low semen utilization rate (L) and high semen utilization rate (H) groups. Red indicates the H group and blue indicates the L group. Boxes show the medians and the interquartile ranges (IQRs), the whiskers denote the lowest and highest values that were within 1.5 times the IQR from the first and third quartiles, and outliers are shown as individual points. (D) Nonmetric multidimensional scaling (NMDS) plot. The red circles indicate the H group, and the blue triangles indicate the L group. (E) Interindividual variations were determined by average weighted-UniFrac distances between individuals in the L and H group boars. Intraindividual variations were determined by distance-paired H and L group boars. ^***p < 0.001.

FIGURE 2

Figure 2. Microbial profiles at the phylum and genus level. (A) Barplot of relative abundance at the phylum level and (B) genus level for the L and H groups. The y-axis represents the relative abundance of each phylum for the two groups. (C) Linear discriminant analysis (LDA) score for discriminated genera in the L and H groups. The LDA score is calculated by LEfSe. The value suggests that it is increased in the two groups (p < 0.05, Wilcoxon rank-sum test, LDA > 2).

Function Prediction of Gut Microbiota

Based on 16S rRNA sequences, PICRUSt analysis could be used to predict the potential functions of gut microflora. By comparing the KEGG database, it was found that the 49 metabolic pathways in the level 3 KEGG function classes are different between the L and H semen utilization rate group boars (Figure 3). The top three metabolic pathways with the largest differences between the two groups were related to biotin metabolism, tropane, piperidine and pyridine alkaloid biosynthesis, and xylene degradation. Genes for biotin metabolism, and tropane, piperidine, and pyridine alkaloid biosynthesis were significantly depleted in the L group in comparison with the H group. Nevertheless, xylene degradation functions were overrepresented in the L group when compared with the H group. It is worth mentioning that an increase in the function of dioxin degradation was also found in the L group.

FIGURE 3

Figure 3. Analyses of pathways predicted by PICRUSt. The Mann-Whitney U test based on the PICRUSt data set revealed differentially enriched bacterial functions associated either with the L group (blue) or the H group (red).

Analysis of Metabolites of Gut Microbiota

To analyze the metabolites of the gut microbiota in boars of low or high semen utilization rate, we then focused on the fecal SCFAs and lactic acid. The results indicated that the concentration of fecal valerate, isobutyrate, isovalerate, and total acid BCFA from L group boars were higher (p < 0.05) than those from the H group boars, whereas no significant difference in the fecal level of acetate, propionate, butyrate, total SCFAs, and lactic acid was found between the L and H groups (Figure 4).

FIGURE 4

Figure 4. Content of (A) short-chain fatty acids (SCFA) and (B) lactic acid in the feces of different semen utilization rate group boars. n = 11 boars in the L group and 29 boars in the H group. ^*p < 0.05.

Intestinal Permeability and Systemic Inflammation Levels of Boars

To determine the intestinal permeability of boars, we assessed endotoxin, zonulin, and diamine oxidase in the plasma. As shown in Figure 5, the plasma endotoxin level of boars with low semen utilization was significantly increased compared to boars with high semen utilization (Figure 5A). Similarly, the change trends of plasma zonulin and diamine oxidase concentration was consistent with endotoxin (Figures 5B,C). We further four parameters related to the immune activation and system inflammation. The results showed that the concentrations of plasma lipocalin-2 were significantly lower (p < 0.05) in the high semen utilization rate group boars; a similar trend was observed in the proinflammatory cytokines IL-6 (p = 0.054) and TNF-α (p = 0.077). However, different from the trend of proinflammatory cytokines, the concentration of anti-inflammatory mediator IL-10 (p = 0.068) in the plasma of boars with high semen utilization rate showed an increasing tendency (Figures 5D–G).

FIGURE 5

Figure 5. Plasma concentrations of (A) endotoxin, (B) zonulin, (C) diamine oxidase, (D) Lipocalin-2, (E) interleukin (IL)-6, (F) tumor necrosis factor (TNF)-α, and (G) IL-10 of different semen quality utilization rate group boars. Data are presented as mean ± SEM (n = 10 boars in the L group and 27 boars in the H group). ^*p < 0.05, ^**p < 0.01.

Correlations Between Parameters Related to Boar Health

Spearman correlation analysis revealed that Sphingobium could be an important genus contributing to the intestinal permeability defect. When combining all the data together, the analysis showed a positive, and yet statistically not significant, association trend between plasma endotoxin and the abundance of Sphingobium (Figure 6A). However, after focusing on those boars’ samples with the abundance of Sphingobium > 0% (The relative abundance of Sphingobium in each sample is shown in Supplementary Figure 1), the results revealed a strong positive correlation (n = 19, r = 0.519, p = 0.023; Figure 6B), suggesting that high abundance of both Sphingobium or the other unidentified genera within Proteobacteria plays critical roles in regulating intestinal permeability. The plasma endotoxin concentration showed a positive correlation with plasma IL-6 (n = 36, r = 0.873, p < 0.001) (Figure 6C). Additionally, there was a weak positive correlation trend between the abundance of Sphingobium and fecal BCFAs (Figure 6D). Spearman correlation analysis was used to analyze the potential connections between clinical status and the semen parameters of boars (Figure 6E). The semen utilization rate was positively correlated with the genus RFN20 (p < 0.05) but negatively correlated with the genera [Ruminococcus] (p < 0.05) and Sphingobium (p < 0.05), plasma endotoxin levels (p < 0.05), plasma diamine oxidase levels (p < 0.05), plasma IL-6 levels (p < 0.05), and plasma lipocalin-2 levels (p < 0.05). The abnormal sperm rate was positively correlated with the plasma diamine oxidase levels (p < 0.05). In addition, the functional sperm number was positively correlated with the genus Paludibacter.

FIGURE 6

Figure 6. Spearman correlation analysis. (A) Correlation analysis between plasma endotoxin and the abundance of Sphingobium (n = 36). (B) Correlation analysis between plasma endotoxin and the abundance of Sphingobium (n = 19). (C) Correlation between plasma endotoxin and plasma IL-6. (D) Correlation between the abundance of Sphingobium and fecal BCFAs. (E) Heat map of the Spearman r correlations between the clinical status and semen parameters of boars. ^*p < 0.05.

Discussion

The association between gut microbiota and host health has been widely studied. However, the effects of gut microbiota on semen quality were still less reported. In this study, we present evidence of such relationships in boars. The most significant findings of the present study are: (i) taxonomic and functional changes in the gut microbiome with boars in different semen utilization rate, and (ii) gut permeability dysfunction and altered immune status in low semen utilization rate boars. These findings suggested a novel relationship between gut microbiota and semen utilization rate in boar.

Using 16S ribosomal RNA gene sequencing analysis, we investigated both phylum and genus level changes in the boar gut microbiota. The results showed that different semen utilization rates of Duroc boars have diverse intestinal microbiota composition. The Firmicutes and Bacteroidetes phyla were the maximum abundant in total sequences, being similar to previous findings in the rectal microbiota of pigs (Wang et al., 2019). In addition, we also observed that Sphingobium and [Ruminococcus] were enriched in the L group, while Paludibacter and RFN20 were enhanced in the H group boars. We hypothesized that analyzing the correlation of semen utilization rate and taxa abundance would be a better, more precise way to distinguish bacteria that had significant impacts on semen quality. Spearman correlation analysis showed that the genus Sphingobium and [Ruminococcus] were negatively correlated with the semen utilization rate, while the genus RFN20 was positively correlated with it. Sphingobium can degrade many toxic organic pollutants, such as polycyclic aromatic compounds (Wittich et al., 2007). We found that the relative abundance of Sphingobium in the L group increased, and at the same time, the expression of metabolizing pollutants signaling pathways also increased. This means that boars in the L group may be exposed to more environmental pollutants, and pollutants may be a potential risk factor for the deterioration of boar semen quality (Bolden et al., 2015). [Ruminococcus], Paludibacter and RFN20 are capable of degrading dietary fiber and fermenting complex carbohydrates (Ze et al., 2012; Fernandes et al., 2014; Cui et al., 2019). SCFAs are mainly fermented products that escape the absorption of carbohydrates and proteins in the small intestine during digestion (Rios-Covian et al., 2016). SCFAs, especially butyrate, is beneficial in regulating the gut barrier (Huang et al., 2015) and also protects against inflammation (Flint et al., 2012). However, in this study, alterations in gut permeability and inflammatory status in the L group boars may not be related to SCFAs because the data showed that there was no significant difference in fecal concentrations of SCFAs between the L and H groups. However, levels of fecal BCFAs in the L group were higher than those of the H group boars. BCFAs can be considered as markers of proteolysis in the colon. It is believed that protein fermentation can result in toxic or potentially toxic metabolites such as ammonia, BCFAs (e.g., isobutyrate and isovalerate), and phenolic compounds (Macfarlane and Macfarlane, 2012; Nyangale et al., 2012). The increase in BCFA content in the L group may be owing to the enrichment of Proteobacteria. Bacteria of Proteobacteria are mainly responsible for the use of amino acids (Booijink, 2009), and in our study, we found a weak positive correlation tendency between the abundance of Sphingobium and fecal BCFA.

Likewise, we found an increased abundance of Gram-negative Proteobacteria (Alphaproteobacteria, Sphingomonadales, Sphingomonadaceae, and Sphingobium) in the L group boars. Proteobacteria is a minor component in the balanced gut-associated microbial community (Eckburg et al., 2005). Expansive of Proteobacteria is observed during conditions of low-level intestinal inflammation, including irritable bowel syndrome and metabolic syndrome (Morgan et al., 2012). Numerous studies have shown that the abnormal expansion of Proteobacteria is a potential diagnostic microbial characteristic of gut microbiota imbalance and epithelial dysfunction (Shin et al., 2015; Litvak et al., 2017). Then, we evaluated boar intestinal permeability and systemic inflammatory response. We found that the contents of the three indicators of intestinal permeability of boars with low semen utilization rates were increased. Furthermore, consistent with what was mentioned before, we found a positive correlation between plasma endotoxin and the abundance of genus Sphingobium.

Circulating LPS activates the immune system via toll-like receptors, which result in proinflammatory cytokines like IL-6 or TNF-α and acute-phase protein production (Kimball et al., 2003; de La Serre et al., 2010; Rakhshandeh and De Lange, 2012). In this study, we found that the plasma endotoxin concentration of boars with low semen utilization was significantly increased, accompanied by an increase in the level of systemic inflammation. Spearman analysis showed that the concentration of plasma endotoxin was positively correlated with IL-6. The above findings are consistent with previous studies. Proinflammatory conditions may lead to increased ROS formation, and sperm membrane lipid peroxidation decreased motility and increased DNA damage (Urata et al., 2001; Kasturi et al., 2008). Besides, male hypogonadism may be associated with the chronic inflammation that inhibits testosterone synthesis (Kalyani and Dobs, 2007). Both LPS and proinflammatory cytokines can affect the production of sex steroids. Nevertheless, decreased testosterone is associated with many adverse reactions, including reduced sperm production, libido, and sexual function (Matsumoto, 2002; Bassil and Morley, 2010).

Additionally, we also explored the function of the gut microbiota of boars with different semen utilization rates. Among these functional pathways, we noticed that xylene degradation and dioxin degradation pathways are increasingly expressed in L group boars. This may be related to the increased abundance of Sphingobium. Members of this genus are mostly found in saprophytic soils and water, and many isolates are involved in the biodegradation of many toxic organic pollutants such as nonylphenol (Ushiba et al., 2003), phenanthrene (Prakash and Lal, 2006), and polycyclic aromatic compounds (Wittich et al., 2007). Environmental pollutants of xylene and dioxins have carcinogenic effects and reproductive toxicity on animals (Bolden et al., 2015; Johnson et al., 2020). These indicate that the boar semen quality is affected by environmental pollutants. The above analysis further confirms that the gut microbial structure and microbial metabolism are shaped in boars with low semen utilization.

Our findings supported the changes in gut microbiota composition and function of the different semen utilization rate boars. However, the reasons for these changes are still unclear and needs further investigation. Additionally, the semen utilization rate is affected by many factors. Through this study, we were unable to determine the specific mechanism by which gut microbiota changed semen parameters. Nonetheless, given the overall importance of gut microbiota to male fertility, there is no doubt that further research is needed to clarify these potential mechanisms. Our research exhibited that the structure and function of the gut microbiota of boars with different semen utilization rates have changed significantly. Specifically, the accumulation of Sphingobium and other Proteobacteria bacteria in the gut of the L group boars may contribute to the increase in plasma endotoxins, which in turn will cause systemic inflammation and reduce the quality of boar semen. This study suggested the potential correlation of the semen utilization rate of boars with gut microbiota, intestinal permeability, and inflammation status stress. This investigation provided some novel insights into the differences between the gut microbiota and microbial function of boars with different semen utilization rates. We reported first the vital effects of gut microbiota on boar semen quality and provided a new perspective for understanding male fertility.

Data Availability Statement

All sequencing data were uploaded to the NCBI database with accession number PRJNA645685.

Ethics Statement

The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Huazhong Agricultural University (permit number HZAUSW-2018-014). Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

JP conceived and designed the experiments, and wrote and revised the manuscript. SJ conceived and designed the experiments, and revised the manuscript. LG performed the experiments, analyzed the data, and wrote part of the manuscript. YW performed the experiments and took part in the data analysis. CW and HW analyzed the data. JT and HS collected the samples. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key Research and Development Project of China (2017YFD0502004), China Agriculture Research System (CARS-36), and Hubei Province Science and Technology Innovation Major Project (2019ABA081).

Conflict of Interest

JT and HS were employed by the company YangXiang Joint Stock.

The remaining 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.

Acknowledgments

We thank the Yangxiang Joint Stock Company for providing boar samples.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.581926/full#supplementary-material

References

Bassil, N., and Morley, J. E. (2010). Late-life onset Hypogonadism: a review. Clin. Geriatr. Med. 26, 197–222. doi: 10.1016/j.cger.2010.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolden, A. L., Kwiatkowski, C. F., and Colborn, T. (2015). New look at BTEX: are ambient levels a problem? Environ. Sci. Technol. 49, 5261–5276. doi: 10.1021/es505316f

CrossRef Full Text | Google Scholar

Booijink, C. C. (2009). Analysis of diversity and function of the human small intestinal microbiota. Wageningen Universiteit.

Google Scholar

Cani, P. D., Amar, J., Iglesias, M. A., Poggi, M., Knauf, C., Bastelica, D., et al. (2007). Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes 56, 1761–1772. doi: 10.2337/db06-1491

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, C. S., Wei, H. K., Wang, P., Yu, H. C., Zhang, X. M., Jiang, S. W., et al. (2019). Early intervention with faecal microbiota transplantation: an effective means to improve growth performance and the intestinal development of suckling piglets. Animal 13, 533–541. doi: 10.1017/S1751731118001611

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, J., Mai, G., Wang, Z., Liu, Q., Zhou, V., Ma, Y., et al. (2019). Metagenomic insights into a cellulose-rich niche reveal microbial cooperation in cellulose degradation. Front. Microbiol. 10:618. doi: 10.3389/fmicb.2019.00618

PubMed Abstract | CrossRef Full Text | Google Scholar

De La Serre, C. B., Ellis, C. L., Lee, J., Hartman, A. L., Rutledge, J. C., and Raybould, H. E. (2010). Propensity to high-fat diet-induced obesity in rats is associated with changes in the gut microbiota and gut inflammation. Am. J. Physiol. Gastrointest. Liver Physiol. 299, G440–G448. doi: 10.1152/ajpgi.00098.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635–1638. doi: 10.1126/science.1110591

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandes, K. A., Kittelmann, S., Rogers, C. W., Gee, E. K., Bolwell, C. F., Bermingham, E. N., et al. (2014). Faecal microbiota of forage-fed horses in New Zealand and the population dynamics of microbial communities following dietary change. PLoS One 9:11. doi: 10.1371/journal.pone.0112846

PubMed Abstract | CrossRef Full Text | Google Scholar

Flint, H. J., Scott, K. P., Louis, P., and Duncan, S. H. (2012). The role of the gut microbiota in nutrition and health. Nat. Rev. Gastroenterol. Hepatol. 9, 577–589. doi: 10.1038/nrgastro.2012.156

PubMed Abstract | CrossRef Full Text | Google Scholar

Gabler, N. K., Spencer, J. D., Webel, D. M., and Spurlock, M. E. (2008). n-3 PUFA attenuate lipopolysaccharide-induced down-regulation of toll-like receptor 4 expression in porcine adipose tissue but does not alter the expression of other immune modulators. J. Nutr. Biochem. 19, 8–15. doi: 10.1016/j.jnutbio.2006.11.014

CrossRef Full Text | Google Scholar

Hedger, M. P. (2011). Toll-like receptors and signalling in spermatogenesis and testicular responses to inflammation-a perspective. J. Reprod. Immunol. 88, 130–141. doi: 10.1016/j.jri.2011.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Song, P. X., Fan, P. X., Hou, C. L., Thacker, P., and Ma, X. (2015). Dietary sodium butyrate decreases postweaning diarrhea by modulating intestinal permeability and changing the bacterial communities in weaned piglets. J. Nutr. 145, 2774–2780. doi: 10.3945/jn.115.217406

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, K. J., Passage, J., Lin, H., Sriram, S., and Budinsky, R. A. (2020). Dioxin male rat reproductive toxicity mode of action and relative potency of 2,3,7,8-tetrachlorodibenzo-p-dioxin and 2,3,7,8-tetrachlorodibenzofuran characterized by fetal pituitary and testis transcriptome profiling. Reprod. Toxicol. 94, 146–162. doi: 10.1016/j.reprotox.2020.02.008

CrossRef Full Text | Google Scholar

Kalyani, R. R., and Dobs, A. S. (2007). Androgen deficiency, diabetes, and the metabolic syndrome in men. Curr. Opin. Endocrinol. Diabetes Obes. 14, 226–234. doi: 10.1097/MED.0b013e32814db856

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasturi, S. S., Tannir, J., and Brannigan, R. E. (2008). The metabolic syndrome and male infertility. J. Androl. 29, 251–259. doi: 10.2164/jandrol.107.003731

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimball, S. R., Orellana, R. A., O’connor, P. M. J., Suryawan, A., Bush, J. A., Nguyen, H. V., et al. (2003). Endotoxin induces differential regulation of mTOR-dependent signaling in skeletal muscle and liver of neonatal pigs. Am. J. Physio. Endocrinol. Metab. 285, E637–E644. doi: 10.1152/ajpendo.00340.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Koh, A., De Vadder, F., Kovatcheva-Datchary, P., and Backhed, F. (2016). From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell 165, 1332–1345. doi: 10.1016/j.cell.2016.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Langille, M. G. I., Zaneveld, J., Caporaso, J. G., Mcdonald, D., Knights, D., Reyes, J., 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

CrossRef Full Text | Google Scholar

Litvak, Y., Byndloss, M. X., Tsolis, R. M., and Baumler, A. J. (2017). Dysbiotic proteobacteria expansion: a microbial signature of epithelial dysfunction. Curr. Opin. Microbiol. 39, 1–6. doi: 10.1016/j.mib.2017.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Luk, G. D., Bayless, T. M., and Baylin, S. B. (1983). Plasma postheparin diamine oxidase. Sensitive provocative test for quantitating length of acute intestinal mucosal injury in the rat. J. Clin. Invest. 71, 1308–1315. doi: 10.1172/jci110881

PubMed Abstract | CrossRef Full Text | Google Scholar

Macfarlane, G. T., and Macfarlane, S. (2012). Bacteria, colonic fermentation, and gastrointestinal health. J. AOAC Int. 95, 50–60. doi: 10.5740/jaoacint.SGE_Macfarlane

PubMed Abstract | CrossRef Full Text | Google Scholar

Maria Moreno-Navarrete, J., Sabater, M., Ortega, F., Ricart, W., and Manuel Fernandez-Real, J. (2012). Circulating zonulin, a marker of intestinal permeability, is increased in association with obesity-associated insulin resistance. PLoS One 7:e37160. doi: 10.1371/journal.pone.0037160

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumoto, A. M. (2002). Andropause: clinical implications of the decline in serum testosterone levels with aging in men. J. Gerontol A Biol. Sci. Medi. Sci. 57, M76–M99. doi: 10.1093/gerona/57.2.M76

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, X. C., Tickle, T. L., Sokol, H., Gevers, D., Devaney, K. L., Ward, D. V., et al. (2012). Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 13:R79. doi: 10.1186/gb-2012-13-9-r79

PubMed Abstract | CrossRef Full Text | Google Scholar

Nyangale, E. P., Mottram, D. S., and Gibson, G. R. (2012). Gut microbial activity, implications for health and disease: the potential role of metabolite analysis. J. Proteome Res. 11, 5573–5585. doi: 10.1021/pr300637d

CrossRef Full Text | Google Scholar

Omu, A. E. (2013). Sperm parameters: paradigmatic index of good health and longevity. Med. Princ. Pract. 22, 30–42. doi: 10.1159/000354208

PubMed Abstract | CrossRef Full Text | Google Scholar

Prakash, O., and Lal, R. (2006). Description of Sphingobium fuliginis sp. nov., a phenanthrene-degrading bacterium from a fly ash dumping site, and reclassification of Sphingomonas cloacae as Sphingobium cloacae comb. nov. Int. J. Syst. Evol. Microbiol. 56, 2147–2152. doi: 10.1099/ijs.0.64080-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rakhshandeh, A., and De Lange, C. F. M. (2012). Evaluation of chronic immune system stimulation models in growing pigs. Animal 6, 305–310. doi: 10.1017/s1751731111001522

PubMed Abstract | CrossRef Full Text | Google Scholar

Rastelli, M., Cani, P. D., and Knauf, C. (2019). The gut microbiome influences host endocrine functions. Endocr. Rev. 40, 1271–1284. doi: 10.1210/er.2018-00280

PubMed Abstract | CrossRef Full Text | Google Scholar

Rios-Covian, D., Ruas-Madiedo, P., Margolles, A., Gueimonde, M., De Los Reyes-Gavilan, C. G., and Salazar, N. (2016). Intestinal short chain fatty acids and their link with diet and human health. Front. Microbiol. 7:185. doi: 10.3389/fmicb.2016.00185

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez, A. L., Rijsselaere, T., Beek, J., Vyt, P., Van Soom, A., and Maes, D. (2013). Boar seminal plasma components and their relation with semen quality. Syst. Biol. Reprod. Med. 59, 5–12. doi: 10.3109/19396368.2012.725120

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/Gb-2011-12-6-R60

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekirov, I., Russell, S. L., Antunes, L. C., and Finlay, B. B. (2010). Gut microbiota in health and disease. Physiol. Rev. 90, 859–904. doi: 10.1152/physrev.00045.2009

CrossRef Full Text | Google Scholar

Shin, N. -R., Whon, T. W., and Bae, J. -W. (2015). Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33, 496–503. doi: 10.1016/j.tibtech.2015.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Smital, J. (2009). Effects influencing boar semen. Anim. Reprod. Sci. 110, 335–346. doi: 10.1016/j.anireprosci.2008.01.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, P. M., Howitt, M. R., Panikov, N., Michaud, M., Gallini, C. A., Bohlooly, Y. M., et al. (2013). The microbial metabolites, short-chain fatty acids, regulate colonic treg cell homeostasis. Science 341, 569–573. doi: 10.1126/science.1241165

PubMed Abstract | CrossRef Full Text | Google Scholar

Urata, K., Narahara, H., Tanaka, Y., Egashira, T., Takayama, F., and Miyakawa, I. (2001). Effect of endotoxin-induced reactive oxygen species on sperm motility. Fertil. Steril. 76, 163–166. doi: 10.1016/S0015-0282(01)01850-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ushiba, Y., Takahara, Y., and Ohta, H. (2003). Sphingobium amiense sp. nov., a novel nonylphenol-degrading bacterium isolated from a river sediment. Int. J. Syst. Evol. Microbiol. 53, 2045–2048. doi: 10.1099/ijs.0.02581-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Waberski, D., Petrunkina, A. M., and Toepfer-Petersen, E. (2008). Can external quality control improve pig AI efficiency? Theriogenology 70, 1346–1351. doi: 10.1016/j.theriogenology.2008.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Tsai, T., Deng, F., Wei, X., Chai, J., Knapp, J., et al. (2019). Longitudinal investigation of the swine gut microbiome from birth to market reveals stage and growth performance associated bacteria. Microbiome 7:109. doi: 10.1186/s40168-019-0721-7

CrossRef Full Text | Google Scholar

Webel, D. M., Finck, B. N., Baker, D. H., and Johnson, R. W. (1997). Time course of increased plasma cytokines, cortisol, and urea nitrogen in pigs following intraperitoneal injection of lipopolysaccharide. J. Anim. Sci. 75, 1514–1520. doi: 10.2527/1997.7561514x

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, T. E., and Kerr, B. J. (2008). Effect of sodium butyrate on growth performance and response to lipopolysaccharide in weanling pigs. J. Anim. Sci. 86, 442–450. doi: 10.2527/jas.2007-0499

PubMed Abstract | CrossRef Full Text | Google Scholar

Wittich, R. M., Busse, H. J., Kampfer, P., Tiirola, M., Wieser, M., Macedo, A. J., et al. (2007). Sphingobium aromaticiconvertens sp nov., a xenobiotic-compound-degrading bacterium from polluted river sediment. Int. J. Syst. Evol. Microbiol. 57, 306–310. doi: 10.1099/ijs.0.64433-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y. H., Lai, W., Liu, Z. H., Wei, H. K., Zhou, Y. F., Tan, J. J., et al. (2019a). Serum and seminal plasma element concentrations in relation to semen quality in Duroc boars. Biol. Trace Elem. Res. 189, 85–94. doi: 10.1007/s12011-018-1459-y

CrossRef Full Text | Google Scholar

Wu, Y. H., Guo, L. L., Liu, Z. H., Wei, H. K., Zhou, Y. F., Tan, J. J., et al. (2019b). Microelements in seminal and serum plasma are associated with fresh semen quality in Yorkshire boars. Theriogenology 132, 88–94. doi: 10.1016/j.theriogenology.2019.04.002

CrossRef Full Text | Google Scholar

Yu, Y., Lu, L., Sun, J., Petrof, E. O., and Claud, E. C. (2016). Preterm infant gut microbiota affects intestinal epithelial development in a humanized microbiome gnotobiotic mouse model. Am. J. Physiol. Gastrointest. Liver Physiol. 311, G521–G532. doi: 10.1152/ajpgi.00022.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Ze, X., Duncan, S. H., Louis, P., and Flint, H. J. (2012). Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J. 6, 1535–1543. doi: 10.1038/ismej.2012.4

PubMed Abstract | CrossRef Full Text | Google Scholar