- Research
- Open access
- Published:
Multi-omics analysis reveals gut microbiota-ovary axis contributed to the follicular development difference between Meishan and Landrace × Yorkshire sows
Journal of Animal Science and Biotechnology volume 14, Article number: 68 (2023)
Abstract
Background
The mechanism by which Meishan (MS) sows are superior to white crossbred sows in ovarian follicle development remains unclear. Given gut microbiota could regulate female ovarian function and reproductive capacity, this study aimed to determine the role of gut microbiota-ovary axis on follicular development in sows.
Methods
We compared the ovarian follicular development, gut microbiota, plasma metabolome, and follicular fluid metabolome between MS and Landrace × Yorkshire (L × Y) sows. A H2O2-induced cell apoptosis model was used to evaluate the effects of multi-omics identified metabolites on the apoptosis of porcine ovarian granulosa cells in vitro.
Results
Compared with L × Y sows, MS sows have greater ovary weight and improved follicular development, including the greater counts of large follicles of diameter ≥ 5 mm, secondary follicles, and antral follicles, but lesser atretic follicles. The ovarian granulosa cells in MS sows had alleviated apoptosis, which was indicated by the increased BCL-2, decreased caspases-3, and decreased cleaved caspases-3 than in L × Y sows. The ovarian follicular fluid of MS sows had higher concentrations of estradiol, progesterone, follicle-stimulating hormone, luteinizing hormone, and insulin like growth factor 1 than L × Y sows. Gut microbiota of MS sows formed a distinct cluster and had improved alpha diversity, including increased Shannon and decreased Simpson than those of L × Y sows. Corresponding to the enhanced function of carbohydrate metabolism and elevated short-chain fatty acids (SCFAs) in feces, the differential metabolites in plasma between MS and L × Y sows are also mainly enriched in pathways of fatty acid metabolism. There were significant correlations among SCFAs with follicular development, ovarian granulosa cells apoptosis, and follicular fluid hormones, respectively. Noteworthily, compared with L × Y sows, MS sows had higher follicular fluid SCFAs concentrations which could ameliorate H2O2-induced porcine granulosa cells apoptosis in vitro.
Conclusion
MS sows have more secondary and antral follicles, but fewer atretic follicles and apoptotic ovarian granulosa cells, as well as harbored a distinctive gut microbiota than L × Y sows. Gut microbiota may participate in regulating ovarian follicular development via SCFAs affecting granulosa cells apoptosis in sows.
Background
Ovulation rate, the number of oocytes ovulated from matured ovarian follicles per estrus cycle, is the decisive factor affecting the litter size in sows [1, 2]. More than 99% of the developing ovarian follicles undergo atretic degeneration before ovulation in mammalian [3, 4]. Earlier comparative study showed that high fertility Meishan (MS) sows had a higher number of antral/pre-ovulatory follicles than white crossbred sows [5, 6]. Miller et al. found that follicular fluid estradiol concentrations and total estradiol synthesis by granulosa cells were higher in MS sows than in Large-White hybrid sows, which may account for the superiority of MS sows in follicular development [5]. However, the mechanism by which MS sows maintain a superiority in follicular development during the follicular phase remains unclear.
The gut microbiome, as part of hologenome (the collective genomic content of a host and its microbiome), has reported to affect ovarian function and reproductive capacity [7–11]. For sows, dietary prebiotics and probiotics supplementation could improve reproductive performance in association the gut microbiota [12–14]. Moreover, large-scale comparative research between high- and low-reproductive performance sows from different commercial farms suggested that litter size of sows was likely associated with the compositions of the gut microbiota [15]. This close correlation between litter size and gut microbiota had been also found in MS sows [16]. The decrease in litter size of MS sows induced by commercial feeds could be improved by crude fiber supplementation in association with gut microbiota [16]. Gut microbiota had also been proved to mediate the benefits of dietary inulin and cellulose supplementation on ovarian oocyte mature in sows [17]. Overall, gut microbiota may involve in regulation of the litter size of sows by affecting the ovarian follicular development and subsequent ovulation rate.
Given that MS sows can have 11.3 more pre-ovulatory follicles and 10 more released oocytes in an estrus cycle than Large-White hybrid sows [5], the aim of this study was to investigate the role of gut microbiota-ovary axis on follicular development by comparing ovarian follicular development, gut microbiota, plasma metabolome, and follicular fluid metabolome in MS and L × Y sows under the same feeding and housing conditions.
Methods
Experimental design and animals
This animal experiment was carried out in accordance with the recommendations of the Guide for the Care and Use of Laboratory Animals Monitoring Committee of Hubei Province, China, and the protocol was approved by the Scientific Ethics Committee of Huazhong Agricultural University (approval numbers HZAUSW-2018-015).
Eight MS sows (age 95 ± 7 d and weight 26.94 ± 3.48 kg) from the National Conservation Farm of MS Pigs (Jiading, Shanghai) and 7 L × Y sows (age 90 ± 4 d and weight 38.98 ± 4.45 kg) from a commercial farm were housed in the same pig-house and had ad libitum access to water and diets. The diets (Table S1) were prepared according to the National Research Council [18]. At 24 h after detection of the third estrus (the average estrous cycle is 21 d), 7 MS sows and 7 L × Y sows were sacrificed for sampling (one MS sow was not reached puberty and was removed). The blood samples were obtained via anterior vena cava puncture into heparinized vacutainer tubes and then centrifuged (4 °C, 3,500 × g, for 10 min) for plasma. Immediately after slaughter, feces from rectum were collected and placed in liquid nitrogen. The weight and large follicle count (diameter ≥ 5 mm) of left and right ovaries were determined, respectively. Half of the left ovary was fixed in 4% paraformaldehyde solution for measuring ovarian morphology and the rest were placed in liquid nitrogen. The contents of follicles (diameter ≥ 3 mm) from the right ovary were pooled and collected in a tube and allowed to settle for 5 min [19]. The supernatant was removed and centrifuged at 1,900 × g at 4 °C for 30 min to separate granular cells from the follicular fluid [19]. The follicular fluid and granular cells were stored at −80 °C until further analysis.
Reagents
Paraformaldehyde (#P6148, Sigma-Aldrich, St. Louis, MO, USA), DMEM (#D6429, Sigma-Aldrich, St. Louis, MO, USA), FBS (#SH30406, HyClone, Logan, Utah, USA), penicillin/streptomycin (#15140122, Gibco, Grand Island, NY, USA), H2O2 (#323381, Sigma-Aldrich, St. Louis, MO, USA), BCL-2 (#3498S, CST, Danvers, MA, USA), Caspase-3 (#9662S, CST, Danvers, MA, USA), and Cleaved Caspase-3 (#9664S, CST, Danvers, MA, USA).
Hematoxylin and eosin (H&E) staining
H&E staining was conducted as reported [20]. Briefly, after fixation in 4% paraformaldehyde (#P6148, Sigma-Aldrich, St. Louis, MO, USA) overnight at 4 °C, ovary tissues (n = 7) were dehydrated in a graded series of alcohols and then embedded in paraffin wax. The section (10 μm) of embedded ovary tissue was dried at 56 °C for 24 h. Next, slides were rehydrated in a graded series of ethanol. After a brief wash in distilled water, we incubated the slides with hematoxylin solution. The sections were then washed with running tap water to remove excess hematoxylin. Then, we differentiated the sections in 1% acid alcohol for 30 s and then washed them with running tap water for 1 min. This step was followed by an incubation in the eosin counterstain, subsequent dehydration in a graded series of ethanol, and immersion in xylene.
Ovarian follicle counting
Follicle counting was performed as previously described with minor modifications [21]. Briefly, primordial follicles were characterized as having a compact oocyte surrounded by a single layer of flattened granular cells, while primary follicles were identified by the presence of an enlarged oocyte surrounded by a single layer of cuboidal granular cells. Secondary follicles were defined as having an enlarged oocyte surrounded by at least a partial or complete second layer of cuboidal granular cells but no more than four complete layers of cuboidal granular cells. Antral follicles were characterized by the presence of areas of follicular fluid (antrum) or a single large antral space. Follicles at the primordial, primary, and preantral stages of development were deemed atretic if the oocyte was degenerating (convoluted and condensed or fragmented) or absent. To avoid repeat counting a follicle, the numbers of primordial and primary follicles were counted in the randomly selected ten areas of the same size (4 mm2), and the numbers of secondary, antral, and atretic follicles were counted in the whole section of H&E-stained ovary (n = 7).
Culture and treatment of sow granular cells
Fresh sow ovaries (20 selected ovaries from 13 sows at 2–3 estrous cycle) were obtained from a commercial slaughterhouse. Sow ovarian granulosa cells were collected from follicles (about 200 follicles with diameter 3–6 mm in total) using a syringe [22, 23]. The granulosa cells were washed and centrifuged (1,000 × g, 5 min), and then inoculated in a 6-cm dish. Cells were cultured in a complete culture medium containing 89% DMEM (#D6429, Sigma-Aldrich, St. Louis, MO, USA), 10% FBS (#SH30406, HyClone, Logan, Utah, USA), and 1% penicillin-streptomycin (#15140122, Gibco, Grand Island, NY, USA) at 37 °C under 5% CO2 throughout this trial. When the confluence of the cells reached approximately 90%, the cells were transferred to 12-well plate for the next experiment. H2O2-induced apoptosis model of sow granulosa cells was constructed as previously described with minor modifications [24]. Briefly, granular cells were pre-treated with SCFAs (acetate 95.72 μmol/L, propionate 8.65 μmol/L, butyrate 12.56 μmol/L, isobutyrate 6.38 μmol/L, valerate 2.21 μmol/L, and isovalerate 14.91 μmol/L) or the vehicle PBS for 4 h. Then, the granulosa cells were exposed to H2O2 (#323381, Sigma-Aldrich, St. Louis, MO, USA) (200 μmol/L) [24] or the vehicle PBS for another 24 h. Subsequently, we conducted apoptosis analysis of granular cells with the terminal terminal-deoxynucleoitidyl transferase mediated nick end labeling (TUNEL) assay (n = 7).
Apoptosis analysis of granular cells
The TUNEL assay was used to mark apoptotic ovarian granular cells in paraffin section [25]. For paraffin ovary slides (10 μm), after deparaffinize, rehydrate, antigen retrieval, and permeabilization, slides were treated with appropriate amount mix of TDT enzyme, dUTP and buffer (1:5:50) in the tunel kit according to the number of slices and tissue size. These slides were then placed in a flat wet box, incubate at 37 °C for 2 h. Follow DAPI counterstain in nucleus, slides were washed three times with PBS (pH 7.4) in a Rocker device, 5 min each. For granular cells in vitro, after fixation in 4% paraformaldehyde and permeabilization with 0.3% Triton X-100 for 5 min, cells were treated for TUNEL assay in the dark. In the negative control group, recombinant TdT Enzyme was replaced with ddH2O. The cells with red fluorescence were defined as apoptotic cells, and the ratio of TUNEL positive granular cells were determined using computer-assisted morphometry (Nikon ECLIPSE Ti2, Tokyo, Japan).
Western blot analysis
The granular cells were homogenized in RIPA buffer and then on ice for 30 min. After centrifugation at 12,000 × g, 4 °C for 15 min, the supernatants were collected, and their protein concentrations were determined using BCA kit (Thermo Scientific, 23250). An equivalent amount of protein (12 μg) was separated by polyacrylamide gel electrophoresis and transferred to a polyvinylidene fluoride membrane. After TBST cleaning and sealing with 5% skim milk, the membrane was incubated with rabbit polyclonal antibodies BCL-2 (CST, #3498S, 1:1,000), Caspase-3 (CST, #9662S, 1:1,000), and Cleaved Caspase-3 (CST, #9664S, 1:1,000) at 4 °C overnight. The membrane was then incubated with secondary antibody (1:10,000) at room temperature for 1.5 h. Bands were measured by densitometry using image J software and relative protein expression levels were standardized with β-actin.
Quantitative enzyme-linked immunosorbent assay (ELISA)
The concentrations of follicle-stimulating hormone (FSH, #CSB-E06791p, sensitivity 10 mIU/mL, and inter/intra assays CV < 15%), estradiol (#CSB-E06844p, sensitivity 40 pg/mL, and inter/intra assays CV < 15%), progesterone(#CSB-E12869p, sensitivity 0.2 ng/mL, and inter/intra assays CV < 15%), and insulin-like growth factor 1(IGF-1, #CSB-E06829p, sensitivity 1.56 ng/mL, intra assays CV < 8%, and inter assays CV < 10%) in follicular fluid were determined using commercial pig ELISA kits (Cusabio, Wuhan, China) according to the manufacturer’s instructions.
Gut microbiota profiling
After thoroughly homogenizing with zirconia beads, the total genomic DNA of a single sow’s fecal microbiota was extracted using the TGuide S96 Magnetic Soil/Stool DNA kit (TIANGEN, Beijing, China) according to manufacturer’s instructions. Concentration and purity of the DNA were determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Rockland, DE, USA). V4 hypervariable regions of 16S rRNA genes (515F, 5′-GTGYCAGCMGCCGCGGTAA-3′; 806R2, 5′-GGACTACNVGGGTWTCTAAT-3′) was used for PCR amplification. Purified amplicon products were pooled in equimolar and paired end sequenced (2 × 250) on an Novaseq 6000 platform (Illumina, San Diego, USA) at Biomarker Technologies Co., Ltd. (Beijing, China).
The sequencing raw data was analyzed using the QIIME2 platform (version 2021.6) [26]. USEARCH [27] (version 10.0) was employed to cluster sequences into operational taxonomic units (OTUs) with similarity of over 97%. The taxonomy of each OTU representative sequence was determined by RDP Classifier against the 16S rRNA gene database Silva [28] (Release132) and Greengenes [29] (version 13.5). Principal co-ordinates analysis (PCoA) analysis and permutational multivariate analysis of variance (PERMANOVA, with 999 Monte Carlo permutations) were conducted based on Bray–Curtis distances using the package “vegan” in R software (version 4.3.1) [30]. The different genera were identified using linear discriminant analysis (LDA) effect size (LEfSe) analysis [31].
Analysis of gut microbial function and metabolites
The function of gut microbiota was analyzed using PICRUSt2 based on 16S rDNA data and reference genome [32]. The gut microbial metabolites short-chain fatty acid (SCFAs) of acetate, propionate, butyrate, valerate, isobutyrate, and isovalerate in rectal feces were determined using gas chromatography with a modification of the previous method [20]. In brief, 1 g of feces samples was weighed into a 2-mL centrifuge tube with 1 mL of methanol added. After being vortexed for 30 s, the sample was centrifuged (12,000 × g) at 4 °C for 10 min. The supernatant was transferred into a new centrifuge tube (2-mL) and mixed with 0.2 mL 25% metaphosphoric acid. After 30 min at 4 °C, the tube was centrifuged (12,000 × g) again at 4 °C for 10 min. To quantify SCFAs, a calibration curve for SCFAs was constructed.
Untargeted metabolomics profiling
Plasma metabolites were extracted as previously described [20]. Metabolomics analysis of plasma sample was performed by the liquid chromatography-mass spectrometry (LC-MS) system followed machine orders. First, all chromatographic separations were performed using an ultra-performance liquid chromatography system (Waters, Milford, MA, USA). For the reversed phase separation, plasma sample used an ACQUITY UPLC BEH C18 column (2.1 mm ×100 mm, 1.7 μm) and liver tissue sample used an ACQUITY UPLC HSS T3 column (2.1 mm × 100 mm, 1.8 μm). A high-resolution tandem mass spectrometer Xevo G2 XS QTOF (Waters) was used to detect metabolites eluted from the column. The Q-TOF was operated in both positive and negative ion modes. In order to test the stability of the LC-MS during the whole acquisition, a quality control sample (Pool of all samples) was gained after every 10 samples. Peak extraction is mainly achieved through the commercial software Progenesis QI (version 2.2), including peak alignment, peak extraction, normalization, deconvolution, and compound identification. Based on quality control sample information, local polynomial regression fitting signal correction (Quality control–based robust LOESS signal correction, QC-RSC) is performed on the real sample signal [33].
SCFAs-targeted absolute quantification metabolomics
The absolute concentrations of SCFAs in ovarian follicular fluid (n = 7) were determined using liquid chromatography-tandem mass spectrometry (LC-MS/MS, Waters I-class–AB Sciex QTRAP 6500) with multiple reaction monitoring. This system allows accurate quantification of low abundance metabolites (up to pg level). Briefly, sample preprocess is as follows. 1) Preparation of standard curve: Take 7 kinds of SCFAs mixed standards (acetate, propionate, butyrate, valerate, isobutyrate, and isovalerate) and carry out gradient dilution; 2) Take 20 μL sample and standard curve sample, add 60 μL cold MeOH/ACN (2:1) of standard song, shake for 5min, precipitate at −20 °C for 4 h, centrifuge at 20,000 × g at 4°C for 15 min, and take 40 μL supernatant in an EP tube; 3) Derivatization: Add 20 μL 200mmol/L 3-NPH (solvent 50% ACN) and 20 μL 120 mmol/L EDC-6% pyridine mixture (solvent 50% ACN) into the EP tube containing 40 μL supernatant respectively, and place it in a metal bath incubate with shaking at 40 °C for 30 min. Cool to room temperature on ice and centrifuge briefly. Add 80 μL 1000D internal standard (solvent is 10% ACN) to the derivatized EP tube, and mix well; 4) Add 90 μL H2O and 90 μL derivatized sample to the filter plate, centrifuge at 3,000 × g at 4 °C for 5 min, take 90 μL filtered liquid in a new 96-well plate, and load 10 μL on the LC-MS. In MultiQuant software (SCIEX, USA), the default parameters are used for automatic identification and integration of each MRM transition (ion pair), and manual inspection is assisted. The content of SCFAs (ng/mL) were obtained by substituting the integrated peak area of the target index in the sample into the standard curve.
Statistical analysis
Experimental data were analyzed by one-way analysis of variance tests, followed by Fisher’s least significant difference and the Duncan multiple comparison test with GraphPad 8.0 software. Results were presented as mean ± SEM. Significance was presented as *P < 0.05, **P < 0.01, and ***P < 0.001.
Results
Ovary weight and ovarian follicles development between MS and L × Y sows
To determine the role of gut microbiota-ovary axis on follicular development in sows, MS and L × Y sows were housed in the same pig-house and fed with the same diet until the third estrus cycle (Fig. 1A). The absolute and relative ovary weights including the left and right of MS sows were greater than those of L × Y sows (P < 0.05) (Fig. 1B and C). MS sows have a greater count of the follicles with diameter ≥ 5 mm (P < 0.01) (Fig. 1D). Further morphological analysis of H&E-stained ovary sections showed that MS sows have greater counts of secondary (P < 0.01) and antral follicles (P < 0.001), but have a smaller count of atretic follicle (P < 0.001) (Fig. 1F).
The apoptosis level of ovarian granulosa cells between MS and L × Y sows
TUNEL-stained ovary sections showed that granulosa cells from MS sows have fewer TUNEL-positive signal than from L × Y sows (P < 0.001) (Fig. 2A and B). In the negative control group, false positive reactivity was not observed. Correspondingly, ovarian granulosa cells of MS sows have higher level of anti-apoptosis protein BCL-2 (P < 0.001), whereas have lower levels of pro-apoptosis proteins caspases-3 (P < 0.01) and cleaved caspases-3 (P < 0.001) than those of L × Y sows (Fig. 2C–G).
The concentrations of hormones in ovarian follicular fluid between MS and L × Y sows
Follicular fluid was collected from fresh ovaries (Fig. 3A). The ovarian follicular fluid of MS sows has higher concentrations of estradiol (P < 0.001, Fig. 3B), progesterone (P < 0.01, Fig. 3C), FSH (P < 0.05, Fig. 3D), LH (P < 0.01, Fig. 3E), and IGF-1 (P < 0.001, Fig. 3F) than L × Y sows.
The diversity and structure of gut microbiota between MS and L × Y sows
Gut microbiota analysis of rectum feces showed that there are 87 unique OTUs belong to MS sows, 126 unique OTUs belong to L × Y sows, and 873 OTUs shared in both MS and L × Y sows (Fig. 4A). Beta diversity presented by PCoA showed MS sows formed a distinct cluster markedly away from L × Y sows (P = 0.0008, PC1 and PC2 explained 40.77% and 29.86% of the variation, respectively) (Fig. 4B). Compared with L × Y sows, MS sows have improved alpha diversity of gut microbiota including increased Shannon (P < 0.05) and decreased Simpson (P < 0.05). Linked-bar plots showed that MS sows have an obvious different composition of gut microbiota at the levels of phylum (Fig. 4D), family (Fig. 4E), and genus (Fig. 4F), respectively.
The marker bacteria, function, and metabolites of gut microbiota between MS and L × Y sows
As shown in Fig. 5A and B, cladogram of LEfSe and histogram of LDA score showed that Turicibacteraceae, Turicibacterales, Clostridiaceae, Fibrobacter, Oscillospira, and Bacteroidaceae were the marker bacteria of MS sows; Erysipelotrichaceae, Erysipelotrichales, Erysipelotrichi, Veillonellaceae, Lachnospiraceae, Streptococcaceae, Bacilli, and Lactobacillales were the marker bacteria of L × Y sows. Compared with the L × Y sows, gut microbiota of MS sows has an increased function of carbohydrate metabolism (P = 0.001) but have decreased functions of translation (P = 0.041) and membrane transport (P = 0.049) (Fig. 5C). Correspondingly, MS sows have increased concentrations of SCFAs including propionate (P < 0.01), butyrate (P < 0.01), isobutyrate (P < 0.05), valerate (P < 0.05), and isovalerate (P < 0.001) in rectal feces than L × Y sows (Fig. 5D).
Metabolomic profiling of plasma between MS and L × Y sows
Heatmaps of untargeted metabolomic profiling showed that MS sows had an obviously different plasma metabolomic profiling compared with L × Y sows (Fig. 6A and B). PCA plots based on the ion signal intensity showed that the plasma metabolomic profiling of MS sows was differed significantly from those of L × Y sows, including positive ion mode (P = 0.004, Fig. 6C) and negative ion mode (P = 0.0025, Fig. 6D), respectively. As shown in Fig. 6E, these different metabolites in plasma enriched in pathways of fatty acid metabolism, steroid hormone biosynthesis, alpha-linolenic acid metabolism, butanoate metabolism, primary bile acid biosynthesis, fatty acid elongation, estrogen signaling pathway, biosynthesis of amino acids, ABC transporters, and 2-oxocarboxylic acid metabolism, etc. Correspondingly, MS sows have increased concentrations of plasma SCFAs including propionate (P < 0.01), butyrate (P < 0.05), isobutyrate (P < 0.05), valerate (P < 0.05), and isovalerate (P < 0.01) than L × Y sows (Fig. 6F).
Correlation analysis among gut microbiome and follicular development indices
As shown in Fig. 7, the differential gut microbial function were correlated with follicular development indices, especially carbohydrate metabolism was significant (P < 0.01) positively associated with follicles of diameter ≥ 5 mm, antral follicles, estradiol, progesterone, LH, IGF-1, and BCL-2; but was negatively associated with apoptosis, caspases-3, and cleaved caspases-3, respectively. Correspondingly, for SCFAs in feces, acetate was negatively associated with atretic follicles (P < 0.05); propionate was significant (P < 0.05) positively associated with antral follicles, LH, progesterone, estradiol, and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis; butyrate was significant (P < 0.05) positively associated with antral follicles, progesterone, estradiol, and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis; isobutyrate was significant (P < 0.05) positively associated with antral follicles, secondary follicles, IGF-1, FSH, progesterone, and estradiol, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis; valerate was significant (P < 0.001) positively associated with antral follicles, but was significant (P < 0.05) negatively associated with atretic follicles and cleaved caspases-3; isovalerate was significant (P < 0.05) positively associated with antral follicles, secondary follicles, follicles of diameter ≥ 5 mm, IGF-1, LH, progesterone, estradiol, and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis. For SCFAs in plasma, acetate was significant (P < 0.01) positively associated with primary follicles; propionate was significant (P < 0.05) positively associated with antral follicles and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, and caspases-3; butyrate was significant (P < 0.05) positively associated with antral follicles, FSH, progesterone, estradiol, and BCL-2, but was significant (P < 0.05) negatively associated with caspases-3, and apoptosis; isobutyrate was significant (P < 0.05) positively associated with antral follicles, follicles of diameter ≥ 5 mm, and estradiol, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis; valerate was significant (P < 0.001) positively associated with antral follicles, secondary follicles, and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, and apoptosis; isovalerate was significant (P < 0.001) positively associated with antral follicles, secondary follicles, follicles of diameter ≥ 5 mm, IGF-1, progesterone, estradiol, and BCL-2, but was significant (P < 0.05) negatively associated with atretic follicles, cleaved caspases-3, caspases-3, and apoptosis.
Effects of physiological concentration of SCFAs on granulosa cells apoptosis in vitro
As shown in Fig. 8A, SCFAs-targeted absolute quantitative metabolomics analysis was conducted to determine the physiological concentrations of SCFAs in the ovarian follicular fluid collected from fresh ovarian follicles of MS and L × Y sows, respectively. The results of SCFAs-targeted absolute quantitative metabolomics showed that ovarian follicular fluid of MS sows had increased concentrations of propionate (P < 0.001), butyrate (P < 0.01), isobutyrate (P < 0.001), valerate (P < 0.01), isovalerate (P < 0.001), and total SCFAs (P < 0.001) than L × Y sows (Fig. 8B). Porcine ovarian granulosa cells were collected from fresh ovaries and then cultured in vitro (Fig. 8C). Fig. 8D and E showed that H2O2 treatment markedly increased apoptosis of porcine ovarian granulosa cells compared with the Ctrl (P < 0.001), which was ameliorated by SCFAs treatment at the concentrations in follicular fluid (P < 0.01). In the negative control group, false positive reactivity was not observed.
Discussion
Here, we found that MS sows have greater ovary weight and have more follicles with diameter ≥ 5 mm than L × Y sows. Further morphological analysis indicated that MS sows have more secondary and antral follicles but have less atretic follicle than L × Y sows. These results are comparable with an earlier comparative study which showed that hyper-fertility MS sows even had 11.3 more pre-ovulatory follicles and ovulate 10 more occytes in an estrus cycle than Large-White hybrid sows [5]. Ovarian granulosa cell apoptosis-induced follicular atresia results in over 99% of the follicles undergo atretic degeneration before ovulation [3, 4]. For the first time, we demonstrated MS sows had significantly fewer apoptotic granulosa cells in antral follicles than L × Y sows with TUNEL-staining in situ. This extraordinary phenotype was further supported by the attenuated apoptosis signal pathway (increased BCL-2, decreased caspases-3 and cleaved caspases-3) of ovarian granulosa cell in MS sows than in L × Y sows. In sows, ovarian granulosa cells can produce steroid hormones (estradiol and progesterone) and provide an appropriate microenvironment for oocytes to meiosis and maturation [19, 34, 35]. Given the determinative role of granulosa cells in follicular development [34, 35], the high-quality granulosa cells may account for the hyper-ovulation rate of MS sows. Granulosa cells in antral follicles synthesize plenty of FSH receptors and LH receptors on each cell, making them extremely sensitive to FSH and LH [35]. Correspondingly, we found MS sows had increased concentrations of estradiol, progesterone, FSH, LH, and IGF-1 in ovarian follicular fluid than L × Y sows. Given that granulosa cells produce estradiol [34, 35], the reduced apoptotic granulosa cells may account for the higher estradiol concentrations in MS follicular fluid than in hybrid Large White pigs [5]. In addition, FSH together with LH can stimulate granulosa cells synthesize estradiol and progesterone, and also coordinate follicular maturation and ovulation [36]. On the other hand, IGF-1 can synergize with FSH to stimulate granulosa cell proliferation and steroidogenesis [37]. Plasma concentrations of both FSH and LH are also positively associated with the ovulation rate in crossbreed sows [38]. Taken the above evidence together, high-quality ovarian granulosa cells may account for the superiority of MS sows in follicular development and ovulation rate than L × Y sows although the underlying mechanism remains unclear.
Gut microbiota, known as the second genome and a substantial endocrine organ, has been reported to affect host reproductive hormones production and metabolism [7, 39–43]. A recent study that recruited 18 fertile female subjects and 18 patients with infertility matched by age showed marked differences in the abundance of gut microbiota between fertile and infertile females [44]. Moreover, partially hydrolyzed guar gum supplementation improved gut microbial dysbiosis and the success of pregnancy in females with infertility [44]. Jiang et al. [16] found that fecal progesterone and estradiol concentrations were increased in MS sows whose litter sizes and gut microbiota had been improved by dietary crude fiber supplementation. Hence, it could be speculated that gut microbiota may contribute to the differences of follicular hormone concentrations between MS and L × Y sows. Under the same feed and housing conditions, MS sows developed a distinctive gut microbiota of improved structure and diversity than L × Y sows. Turicibacteraceae, Turicibacterales, Clostridiaceae, Fibrobacter, Oscillospira, and Bacteroidaceae were the marker bacteria of MS sows, whereas Erysipelotrichaceae, Erysipelotrichales, Erysipelotrichi, Veillonellaceae, Lachnospiraceae, Streptococcaceae, Bacilli ,and Lactobacillales were the marker bacteria of L × Y sows. Compared with the L × Y sows, gut microbiota of MS sows has an increased function of carbohydrate metabolism but have decreased functions of translation and membrane transport. Correspondingly, MS sows have increased concentrations of SCFAs including propionate, butyrate, isobutyrate, valerate, and isovalerate in rectal feces than L × Y sows. Compared well with our results, the alpha diversity of gut microbiota, the relative abundance of Fibrobacter, and fecal SCFAs were increased in MS sows whose litter sizes had been improved by dietary crude fiber supplementation [16].
To explore the connection between the gut microbiota and the distal ovarian follicular development, we profiled the plasma metabolome features using untargeted metabolomic and demonstrated a markedly distinguishing plasma metabolomics profiling between MS and L × Y sows. These different metabolites in plasma enriched in pathways of fatty acid metabolism, steroid hormone biosynthesis, alpha-linolenic acid metabolism, butanoate metabolism, primary bile acid biosynthesis, fatty acid elongation, estrogen signaling pathway, biosynthesis of amino acids, ABC transporters, and 2-Oxocarboxylic acid metabolism, etc. Correspondingly, MS sows have increased concentrations of plasma SCFAs including propionate, butyrate, isobutyrate, valerate, and isovalerate than L × Y sows. In addition, correlation analysis showed that the concentrations of SCFAs in both feces and plasma were significant negatively associated with the apoptosis of ovarian granulosa cells, indicating that the increased SCFAs may contribute to the porcine granulosa cells apoptosis amelioration observed in MS sows. In proximal tissues such as the intestine, adipose, and pancreas, SCFAs can regulate host immunity, inflammation, and metabolism by activating G protein-coupled receptors (GPCRs) [45]. However, the activation of GPCRs by SCFAs is receptor-preferred and concentration-dependent [46]. In vivo, butyrate regulates porcine granulosa cells to secrete progesterone and estradiol via cAMP signaling pathway in a dose-dependent manner [47]. The effective dose of butyrate on the secretion of progesterone was over 0.05 mmol/L, and the effective dose for estradiol secretion was over 5 mmol/L [47]. These evidences indicated that the in vitro levels of SCFAs may not fully reflect the physiological levels and function in vivo [47]. Earlier study using proton nuclear magnetic resonance spectroscopy showed that the acetate concentrations were 0.34–1.5 mmol/L in ovarian follicular fluids of pig [48]. However, to our knowledge, the physiological concentrations of the others SCFAs in porcine ovarian follicular fluid remain unreported because of the concentrations are too low to be determined using a conventional method [14, 47]. Targeted metabolomics has a pg-level sensitivity for accurate quantification of low-abundance metabolites. With targeted bile acid metabolomics, Gu et al. [49] determined bile acids in human plasma at the nmol/L level. In this study, with SCFAs-targeted absolute quantitative metabolomics, for the first time, we reported that ovarian follicular fluid of MS sows had higher SCFAs concentrations than those of L × Y sows including propionate (8.65 ± 1.55 vs. 4.15 ± 1.29 μmol/L), butyrate (12.56 ± 3.43 vs. 6.56 ± 1.57 μmol/L), isobutyrate (6.38 ± 1.34 vs. 3.54 ± 0.63 μmol/L), valerate (2.21 ± 0.42 vs. 1.56 ± 0.29 μmol/L), isovalerate(14.91 ± 3.60 vs. 6.11 ± 1.90 μmol/L), and total SCFAs (140.42 ± 30.50 vs. 100.00 ± 16.50 μmol/L), respectively. Inconsistent with Gosden et al., we found that the acetate concentration in ovarian follicular fluid was 86.91 ± 18.42 μmol/L, which was markedly lower than 0.34–1.5 mmol/L [48]. This variation may be because of the different detection method. Given that the effect of SCFAs is concentration-dependent [46, 47], we determined whether the SCFAs were sufficient to ameliorate porcine granulosa cells apoptosis in vitro. Not surprisingly, these physiological concentrations of SCFAs in porcine ovarian follicular fluid could ameliorate H2O2-induced porcine granulosa cells apoptosis in vitro. Evidence in vitro suggested that acetate and propionate can activate GPR43 at low concentration (1 mmol/L) to inhibit the apoptosis of cultured cells in vitro; while at high concentration (10 mmol/L), they promote the apoptosis of colorectal cancer cells [50, 51]. In vitro experiments showed that butyric acid could promote the expression of GPCRs genes in porcine ovarian granulosa cells only when it reached 10 mmol/L [47]. It is generally believed that the physiological concentration of SCFAs in porcine follicles cannot reach the threshold of activating GPCRs [14, 52]. In this study, the physiological concentrations of SCFAs in porcine ovarian follicular fluid are also not reach the threshold of activating GPCRs [14, 52]. Considering SCFAs share the GPCRs and have competitive antagonistic/additive relationship against GPCRs [45], the role of GPCRs in SCFAs ameliorating porcine granulosa cells apoptosis remains uncertain. On the other hand, recent studies have shown that affecting the epigenetic modification of cells in the intestine and even the distant organs is a new mechanism by which gut microbial metabolites SCFAs can also regulate biological processes such as cell cycle, cell proliferation, and apoptosis [45, 53–56]. Such as SCFAs can regulate histone acetylation modification by providing acetylation substrates or affecting the activities of histone deacetylases and histone acetylases, thereby regulating cell physiology [45]. In addition, SCFAs have the potential to affect novel modifications of histones-propionylation and butyrylation [54, 56]. Although the mechanism by which physiological concentrations of SCFAs ameliorates porcine granulosa cells apoptosis is complicated and remains to be identified, our results suggested that modulating gut microbiota to produce SCFAs or supplementation with SCFAs directly may be a novel approach to ameliorate ovarian dysfunction induced by granulosa cells apoptosis. For example, accumulating evidence indicated that high-energy diets could decrease oocyte quality and ovarian function, increase granulosa cells apoptosis and follicular atresia, thus potentially decreasing reproductive performance and shortening the reproductive span [14, 57, 58]. For sows, Zhuo et al. [59] found that high-fat (high-energy) diets could increase granulosa cells apoptosis and follicular atresia, reduce SCFAs-producing gut microbes, whereas dietary fibre supplementation restored this loss and increased the relative abundance of representative SCFAs-producing gut microbes. The authors suggested that dietary fibre supplementation protected against high fat feeding-induced granulosa cells apoptosis and follicular atresia in an indirectly approach via gut microbiota-SCFAs-serotonin/melatonin-ovary axis [59]. On the other hand, according to our results, dietary fibre supplementation may decrease high fat feeding-induced follicular atresia directly by gut microbiota-derived SCFAs ameliorating ovarian granulosa cells apoptosis in sows. Overall, these results provide insight into preventing granulosa cells apoptosis-related ovarian dysfunction.
Conclusions
This study provides insight into the regulating effects of gut microbiota-ovary axis on follicular development in sows by a comprehensively comparative study between hyper-ovulation rate MS pigs and normal-ovulation rate L × Y pigs. The high levels of estradiol, progesterone, FSH, LH, and IGF-1, together with the less apoptotic granulosa cells in ovarian follicles may account for the hyper-ovulation rate of MS pigs. MS pigs have higher SCFAs concentrations in porcine ovarian follicular fluid, which may be because of their enhanced carbohydrate metabolism function of gut microbiota than L × Y pigs. These physiological concentrations of SCFAs in porcine ovarian follicular fluid of MS pigs could ameliorate porcine granulosa cells apoptosis in vitro. Therefore, gut microbiota participate in regulating ovarian follicular development via SCFAs affecting granulosa cells apoptosis in sows.
Availability of data and materials
All data generated or analyzed during this study can be made available by the corresponding author upon reasonable request.
Abbreviations
- BCL-2:
-
B-cell lymphoma 2
- ELISA:
-
Enzyme-linked immunosorbent assay
- FSH:
-
Follicle-stimulating hormone
- GPCRs:
-
G protein-coupled receptors
- H&E:
-
Hematoxylin and eosin
- IGF-1:
-
Insulin-like growth factor 1
- LC-MS:
-
Liquid chromatography-mass spectrometry
- LDA:
-
Linear discriminant analysis
- LEfSe:
-
Linear discriminant analysis effect size analysis
- LH:
-
Luteinizing hormone
- L × Y:
-
Landrace × Yorkshire
- MS:
-
Meishan
- OTU:
-
Operational taxonomic units
- PCoA:
-
Principal co-ordinates analysis
- PICRUSt2:
-
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2
- SCFAs:
-
Short chain fatty acids
- TUNEL:
-
Terminal terminal-deoxynucleoitidyl transferase mediated nick end labeling
References
Hunter MG, Picton HM. Effect of hcg administration at the onset of estrus on early embryo survival and development in meishan gilts. Anim Reprod Sci. 1995;38(3):231–8.
Rosendo A, Druet T, Gogué J, Canario L, Bidanel JP. Correlated responses for litter traits to six generations of selection for ovulation rate or prenatal survival in French Large White pigs. J Anim Sci. 2007;85(7):1615–24.
Shen M, Jiang Y, Guan Z, Cao Y, Li L, Liu H, et al. Protective mechanism of FSH against oxidative damage in mouse ovarian granulosa cells by repressing autophagy. Autophagy. 2017;13(8):1364–85.
Manabe N, Goto Y, Matsuda-Minehata F, Inoue N, Maeda A, Sakamaki K, et al. Regulation mechanism of selective atresia in porcine follicles: regulation of granulosa cell apoptosis during atresia. J Reprod Dev. 2004;50(5):493–514.
Miller A, Picton H, Craigon J, Hunter MG. Follicle dynamics and aromatase activity in high-ovulating Meishan sows and in Large-White hybrid contemporaries. Biol Reprod. 1998;58(6):1372–8.
Dufour JJ, Mariana JC. Comparative follicular development in Meishan and Large White gilts during prepubertal periods and its relation to hormonal stimulation. Biol Reprod. 1993;48(5):1020–5.
Moran NA, Sloan DB. The hologenome concept: helpful or hollow? PLoS Biol. 2015;13(12): e1002311.
Koren O, Goodrich Julia K, Cullender Tyler C, Spor A, Laitinen K, Kling Bäckhed H, et al. Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell. 2012;150(3):470–80.
Qi X, Yun C, Sun L, Xia J, Wu Q, Wang Y, et al. Gut microbiota–bile acid–interleukin-22 axis orchestrates polycystic ovary syndrome. Nat Med. 2019;25(8):1225–33.
Shimizu K, Muranaka Y, Fujimura R, Ishida H, Tazume S, Shimamura T. Normalization of reproductive function in germfree mice following bacterial contamination. Exp Anim. 1998;47(3):151–8.
Williams CL, Ybarra AR, Meredith AN, Durrant BS, Tubbs CW. Gut microbiota and phytoestrogen-associated infertility in Southern White Rhinoceros. mBio. 2019;10(2):e00311–19.
Cao M, Li Y, Wu QJ, Zhang P, Li WT, Mao ZY, et al. Effects of dietary clostridium butyricum addition to sows in late gestation and lactation on reproductive performance and intestinal microbiota. J Anim Sci. 2019;97(8):3426–39.
Kritas SK, Marubashi T, Filioussis G, Petridou E, Christodoulopoulos G, Burriel AR, et al. Reproductive performance of sows was improved by administration of a sporing bacillary probiotic (Bacillus subtilis C-3102). J Anim Sci. 2015;93(1):405–13.
Zhuo Y, Cao M, Gong Y, Tang L, Jiang X, Li Y, et al. Gut microbial metabolism of dietary fibre protects against high energy feeding induced ovarian follicular atresia in a pig model. Br J Nutr. 2021;125(1):38–49.
Uryu H, Tsukahara T, Ishikawa H, Oi M, Otake S, Yamane I, et al. Comparison of productivity and fecal microbiotas of sows in commercial farms. Microorganisms. 2020;8(10):1469.
Jiang X, Lu N, Xue Y, Liu S, Lei H, Tu W, et al. Crude fiber modulates the fecal microbiome and steroid hormones in pregnant Meishan sows. Gen Comp Endocrinol. 2019;277:141–7.
Men Z, Cao M, Gong Y, Hua L, Zhang R, Zhu X, et al. Microbial and metabolomic mechanisms mediating the effects of dietary inulin and cellulose supplementation on porcine oocyte and uterine development. J Anim Sci Biotechnol. 2022;13(1):14.
National Research Council. Nutrient requirements of swine: eleventh. revised. Washington, DC: The National Academies Press; 2012.
Costermans NGJ, Soede NM, van Tricht F, Blokland M, Kemp B, Keijer J, et al. Follicular fluid steroid profile in sows: relationship to follicle size and oocyte quality. Biol Reprod. 2020;102(3):740–9.
Xu B, Qin W, Yan Y, Tang Y, Zhou S, Huang J, et al. Gut microbiota contributes to the development of endometrial glands in gilts during the ovary-dependent period. J Anim Sci Biotechnol. 2021;12(1):57.
Wang S, Zheng Y, Li J, Yu Y, Zhang W, Song M, et al. Single-cell transcriptomic atlas of primate ovarian aging. Cell. 2020;180(3):585–600.
Zheng Y, Ma L, Liu N, Tang X, Guo S, Zhang B, et al. Autophagy and apoptosis of porcine ovarian granulosa cells during follicular development. Animals (Basel). 2019;9(12):1111.
Zhang X, Tao Q, Shang J, Xu Y, Zhang L, Ma Y, et al. MiR-26a promotes apoptosis of porcine granulosa cells by targeting the 3β-hydroxysteroid-Δ24-reductase gene. Asian-Australas J Anim Sci. 2020;33(4):547–55.
Zhang J, Ren Q, Chen J, Gao B, Wang X, Zhang Z, et al. Autophagy contributes to oxidative stress-induced apoptosis in porcine granulosa cells. Reprod Sci. 2021;28(8):2147–60.
Gavrieli Y, Sherman Y, Ben-Sasson SA. Identification of programmed cell death in situ via specific labeling of nuclear DNA fragmentation. J Cell Biol. 1992;119(3):493–501.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72.
Wang X, Tsai T, Deng F, Wei X, Chai J, Knapp J, et al. Longitudinal investigation of the swine gut microbiome from birth to market reveals stage and growth performance associated bacteria. Microbiome. 2019;7(1):109.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. 2020;38(6):685–8.
Dunn WB, Broadhurst D, Begley P, Zelena E, Francis-McIntyre S, Anderson N, et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc. 2011;6(7):1060–83.
Palma GA, Argañaraz ME, Barrera AD, Rodler D, Mutto AÁ, Sinowatz F. Biology and biotechnology of follicle development. Sci World J. 2012;2012:938138.
Knox RV. Follicle development in pigs: State of the art. Mol Reprod Dev. 2022. https://doi.org/10.1002/mrd.23576.
Pierce JG, Parsons TF. Glycoprotein hormones: dtructure and function. Annu Rev Biochem. 1981;50(1):465–95.
Liu J, Koenigsfeld AT, Cantley TC, Boyd CK, Kobayashi Y, Lucy MC. Growth and the initiation of steroidogenesis in porcine follicles are associated with unique patterns of gene expression for individual componentsof the ovarian insulin-like growth factor system. Biol Reprod. 2000;63(3):942–52.
Knox R, Vatzias G, Naber C, Zimmerman D. Plasma gonadotropins and ovarian hormones during the estrous cycle in high compared to low ovulation rate gilts. J Anim Sci. 2003;81(1):249–60.
Evans JM, Morris LS, Marchesi JR. The gut microbiome: the role of a virtual organ in the endocrinology of the host. J Endocrinol. 2013;218(3):R37–47.
Clarke G, Stilling RM, Kennedy PJ, Stanton C, Cryan JF, Dinan TG. Minireview: gut microbiota: the neglected endocrine organ. Mol Endocrinol. 2014;28(8):1221–38.
Qi X, Yun C, Pang Y, Qiao J. The impact of the gut microbiota on the reproductive and metabolic endocrine system. Gut Microbes. 2021;13(1):1894070.
Han Q, Wang J, Li W, Chen Z-J, Du Y. Androgen-induced gut dysbiosis disrupts glucolipid metabolism and endocrinal functions in polycystic ovary syndrome. Microbiome. 2021;9(1):101.
Silva MSB, Giacobini P. Don’t trust your gut: when gut microbiota disrupt fertility. Cell Metab. 2019;30(4):616–8.
Komiya S, Naito Y, Okada H, Matsuo Y, Hirota K, Takagi T, et al. Characterizing the gut microbiota in females with infertility and preliminary results of a water-soluble dietary fiber intervention study. J Clin Biochem Nutr. 2020;67(1):105–11.
van der Hee B, Wells JM. Microbial regulation of host physiology by short-chain fatty acids. Trends Microbiol. 2021;29(8):700–12.
Le Poul E, Loison C, Struyf S, Springael J-Y, Lannoy V, Decobecq M-E, et al. Functional characterization of human receptors for short chain fatty acids and their role in polymorphonuclear cell activation. J Biol Chem. 2003;278(28):25481–9.
Lu N, Li M, Lei H, Jiang X, Tu W, Lu Y, et al. Butyric acid regulates progesterone and estradiol secretion via cAMP signaling pathway in porcine granulosa cells. J Steroid Biochem Mol Biol. 2017;172:89–97.
Gosden RG, Sadler IH, Reed D, Hunter RH. Characterization of ovarian follicular fluids of sheep, pigs and cows using proton nuclear magnetic resonance spectroscopy. Experientia. 1990;46(10):1012–5.
Gu Y, Wang X, Li J, Zhang Y, Zhong H, Liu R, et al. Analyses of gut microbiota and plasma bile acids enable stratification of patients for antidiabetic treatment. Nat Commun. 2017;8(1):1785.
Jan G, Belzacq AS, Haouzi D, Rouault A, Métivier D, Kroemer G, et al. Propionibacteria induce apoptosis of colorectal carcinoma cells via short-chain fatty acids acting on mitochondria. Cell Death Differ. 2002;9(2):179–88.
Pingitore A, Gonzalez-Abuin N, Ruz-Maldonado I, Huang GC, Frost G, Persaud SJ. Short chain fatty acids stimulate insulin secretion and reduce apoptosis in mouse and human islets in vitro: role of free fatty acid receptor 2. Diabetes Obes Metab. 2019;21(2):330–9.
Koh A, De Vadder F, Kovatcheva-Datchary P, Bäckhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. 2016;165(6):1332–45.
Thomas SP, Denu JM. Short-chain fatty acids activate acetyltransferase p300. eLife. 2021;10:e72171.
Kebede AF, Nieborak A, Shahidian LZ, Le Gras S, Richter F, Gómez DA, et al. Histone propionylation is a mark of active chromatin. Nat Struct Mol Biol. 2017;24(12):1048–56.
Krautkramer KA, Kreznar JH, Romano KA, Vivas EI, Barrett-Wilt GA, Rabaglia ME, et al. Diet-microbiota interactions mediate global epigenetic programming in multiple host tissues. Mol Cell. 2016;64(5):982–92.
Chen Y, Sprung R, Tang Y, Ball H, Sangras B, Kim SC, et al. Lysine propionylation and butyrylation are novel post-translational modifications in histones. Mol Cell Proteomics. 2007;6(5):812–9.
Skaznik-Wikiel ME, Swindle DC, Allshouse AA, Polotsky AJ, McManaman JL. High-fat diet causes subfertility and compromised ovarian function independent of obesity in mice. Biol Reprod. 2016;94(5):108.
Reynolds KA, Boudoures AL, Chi MM, Wang Q, Moley KH. Adverse effects of obesity and/or high-fat diet on oocyte quality and metabolism are not reversible with resumption of regular diet in mice. Reprod Fertil Dev. 2015;27(4):716–24.
Zhou P, Zhao Y, Zhang P, Li Y, Gui T, Wang J, et al. Microbial mechanistic insight into the role of inulin in improving maternal health in a pregnant sow model. Front Microbiol. 2017;8:2242.
Funding
This work was supported by the National Natural Science Foundation of China (32230099 and 31925037) and Hubei Hongshan Laboratory (2021hszd018).
Author information
Authors and Affiliations
Contributions
BX, LM, and XY designed the research. BX, YT, SZ, and JH performed the research. BX, WQ, YW, YT, SZ, and XY analyzed the data. BX and XY wrote the paper with the help of all authors. All authors read and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The animal model and experimental procedures used in this experiment were approved by the Scientific Ethics Committee of Huazhong Agricultural University (approval numbers HZAUSW-2018-015).
Consent for publication
Not applicable.
Competing of interest
The authors declare that they have no competing interests.
Supplementary Information
Additional file1:
Table S1. Ingredients and nutrients composition of the diet used in this trial.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Xu, B., Qin, W., Chen, Y. et al. Multi-omics analysis reveals gut microbiota-ovary axis contributed to the follicular development difference between Meishan and Landrace × Yorkshire sows. J Animal Sci Biotechnol 14, 68 (2023). https://doi.org/10.1186/s40104-023-00865-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40104-023-00865-w