Metagenomic Characterization and Volatile Compounds Determination in Rumen from Saanen Goat Kids Fed Olive Leaves

Simple Summary The aim of this study was to characterize the rumen microbiota of Saanen goat kids fed olive leaves through a high-throughput approach based on 16S rRNA gene sequencing; furthermore, the parallel characterization of rumen volatile profile by solid-phase microextraction coupled with gas chromatography-mass spectrometry has been performed. Twenty goat kids were randomly assigned to two groups. The first group received a basal diet, while in the second one the diet was supplemented with olive leaves. The results showed the dietary supplementation to be able to affect the microbial community in the rumen. Significant differences were specifically observed between the two groups at genera and even family levels characterized by a higher abundance of cellulolytic bacteria in the rumen of goat kids fed olive leaves. In addition, the analysis of volatile compounds at the rumen level has allowed us to highlight differences in relation to the diet and the presence, in the rumen of goat kids fed olive leaves, of compounds indicative of health status. Abstract The accumulation and disposal of by-products deriving from the agro-food industry represents a problem both from an economic and environmental point of view. The use of these matrices in zootechnical nutrition could represent a feasible solution. The aim of the study was to examine the effect of a diet containing olive leaves (OL), a by-product of the olive industry, on the ruminal microbial community of Saanen goat kids and on volatile organic compounds (VOCs) produced during the digestion. Twenty goat kids were randomly divided into two groups of ten goat kids each. The control group (CTR) was fed with a standard diet, while the experimental group (OL+) received a custom-formulated diet containing 10 % OL on a dry matter (DM) basis. After 30 days of trial, genomic DNA was extracted from the rumen liquor and prepared for 16S rRNA-gene sequencing to characterize the rumen microbiota; furthermore, rumen VOCs were also characterized by solid-phase microextraction coupled with gas chromatography-mass spectrometry. The Shannon’s alpha index was not significantly different between the two groups, on the contrary, Bray-Curtis (p < 0.01) and Jaccard (p < 0.01) distances evidenced that feed affected microbial community. Eleven genera were influenced by OL supplementation, with a significant increase (p < 0.05) in Paludibacter, Fibrobacter, Sphaerochaeta Christensenella, Rikenella, Oligosphaera, Candidatus Endomicrobium, Anaerovorax, and Atopobium was observed, while the percentages of Bacteroides and Selenomonas were reduced (p < 0.05). Differences were also observed between the two groups at the family level (p < 0.004). Fibrobacteriaceae, Christensenellaceae, Coriobacteriaceae, Oligosphaeraceae, Candidatus Endomicrobium, and Planctomycetaceae were significantly higher (p < 0.05) in goat kids fed OL diet compared to CTR, while the levels of other identified families, Succinivibrionaceae and Bifidobacteriaceae, were opposite (p < 0.05). Finally, results showed that the main phyla in both groups were Bacteroidetes and Firmicutes; however, no significant differences in the relative abundance of any phyla were observed between the two groups. In addition to what has been reported, the analysis of VOCs at the rumen level showed the ability of the OL integration to induce an increase in hexanoic acid and a parallel decrease in decanal. Furthermore, only in OL+ samples there was the accumulation of α-terpineol to which a wide range of interesting biological properties is attributed. The presence of VOCs associated with health status suggests a favorable role of OL in preserving and improving animal welfare.

microbial population; for that reason, this approach was useful to predict the association between metagenome and metabolic functions using annotations collected in databases of microbial functional genes [22].
To our knowledge, no studies concerning the effect of OL on goat's ruminal bacteria have been previously performed. Therefore, the aim of the present study was to identify possible changes due to OL dietary intake in goat kids rumen microbiota, using a metagenomic approach, and to characterize the production and accumulation of VOCs derived from microbial activity.

Materials and Methods
This study was performed on a commercial farm and animals were handled following the national legislation on animal welfare (Council Directive 2008/119/EC) [23], and then slaughtered in compliance with the Council Regulation 1099/2009 of the European Union [24] on the protection of animals at the time of killing. For the scope of the study, animals did not undergo breeding practices other than those commonly adopted; for this reason, it was not considered necessary to provide further ethical declarations.

Animal Management and Sampling
All twenty Saanen goat kids, homogeneous for age (90 d of age), sex (male), and weight (14.50 ± 1.65 kg), were randomly allotted into two groups (10 animals per group): a control group (CTR) and an experimental group (OL+) whose diet was integrated with OL collected from about 50 olive trees (Olea europaea), more than 50 years of age, belonging to Leccino variety and grown in a limited area of Abruzzo, Italy (Castellalto, Teramo, Italy). The OL chemical composition is reported in a previous study [25]. The trial took place in April and lasted 30 days. The goat kids were housed for the entire trial period in two separate areas of free housing with bunks on straw, a drinking trough, and access to an identical feeding area. The diets were formulated to be isoenergetic; every day, each goat kid of both the group received alfalfa hay administered ad libitum with 1 kg/head/d of a concentrate ( Table 1). The diet of the experimental group was mixed with OL (10 g/100 g on a dry matter (DM) basis of the whole diet). OL harvested during the winter pruning were dried in a ventilated oven at 60 • C until reaching a moisture of 10%, ground, and mixed with the concentrate.  (AOAC, 2000) [26]; neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL) were determined by the detergent procedures of Goering and Van Soest [27].
Within 60-90 min from slaughtering (120 d of age), as reported by Biscarini et al. [22], 500 mL of the whole rumen content (consisting of a mixture of liquid and solid fractions) from the dorsal, central, and ventral of each goat kids was immediately collected, pooled, and filtered through four layers of cheesecloth, and then collected in 50 mL tubes and stored at −20 • C until further analysis (about 1 month).

DNA Extraction and Illumina MiniSeq Sequencing
Total microbial DNA was extracted and purified from the rumen liquid (n = 10 per group) by the Maxwell 16 LEV simplyRNA Blood Kit (Promega, Madison, WI) following the manufacturer's instructions without using the DNaseI. Briefly, 1 mL of rumen fluid was transferred to 2 mL tubes where 0.5 mL of phosphate-buffered saline (PBS) was added. The tubes were then quickly mixed by vortex and centrifuged for 10 min at 300 g to pellet the vegetal matter. A volume of 300 µL supernatant was transferred to the cartridges provided by the kit followed by the automatic extraction/purification workflow on the Maxwell 16 instrument (Promega, Madison, WI). Then, DNA was quantified by Qubit dsDNA HS assay (Thermo Fisher Scientific, Waltham, MA) and diluted to 2 ng/µL for library preparation by means of the Swift Amplicon 16S + ITS panel Kit (Swift Bioscience, Ann Arbor, MI, USA) that targets all the hypervariable regions (V1-V9) of the bacterial 16S rRNA gene. DNA libraries were sequenced on the MiniSeq Illumina platform (Illumina, San Diego, CA, USA) according to the manufacturer's instructions. A 20-strain even mix of bacterial genomic DNA (ATCC MSA-1002) was used as sequencing and bioinformatic analysis control.

Evaluation of Volatile Compounds in Rumen
The characterization of volatile organic compounds (VOCs) in goat kids rumen samples, was performed according to previous research [12,28], by using a divinylbenzenecarboxen-polydimethylsiloxane solid phase microextraction (SPME) fiber. Each sample of rumen was weighed and placed in a test tube. After pressing, 10 mL of rumen fluid was immediately transferred into a glass vial and mixed with 5 mL of saturated NaCl solution (360 g/L) and 10 µL of internal standard solution (4-methyl-2-heptanone; 10 mg/kg in ethanol). The vials were sealed with a polytetrafluoroethylene-silicone septum (Supelco, Bellefonte, PA, USA) and stirred for 30 min at 39 • C to simulate the rumen environment. The VOCs were analyzed by a gas chromatograph (Clarus 580; Perkin Elmer, Waltham, MA, USA) coupled with a mass spectrometer (SQ8S; Perkin Elmer). The gas chromatograph was equipped with an Elite-5MS column (length × internal diameter: 30 × 0.25 mm; film thickness: 0.25 µm; Perkin Elmer). The fiber was thermally desorbed into the GC injector for 30 min in a splitless mode at 250 • C. The column oven initial temperature was 50 • C (held for 1min), followed by a first ramp of 3 • C/min to reach 200 • C (held for 1 min) and a second ramp of 15 • C/min to obtain 250 • C (held for 15 min). Helium was used as the carrier gas with a flow rate of 1 mL/min. Analytes were identified as previously described [29]. The data were expressed as relative abundance, as a percentage of each compound on the sum of the total VOCs detected.
Alpha diversity was used to define the variety and the abundance of species in the rumen. The indices used to calculate alpha diversity were Shannon and Chao1. Bray-Curtis and Jaccard indices were used to analyze similarities or differences among the groups (beta diversity) at the genus, family, and phylum level. The contribution of individual taxa to the overall Bray-Curtis dissimilarity index was computed by SIMPER (similarity percentage, vegan).
Permutational multivariate analysis of variance (Adonis) was used for beta diversity on both Bray-Curtis and Jaccard distances to evaluate quantitative and qualitative dissimilarity between the two groups, respectively. Bray-Curtis dissimilarity index ranges are from 0 to 1; values close to 0 mean experimental groups share the exact same number of each type of species and values close to 1 share none of the same kind of species. Conversely, the Jaccard similarity index varies from 0 to 1 (File S1). The closer to 1, the more similar the two sets of data. Moreover, the Kruskal-Wallis non-parametric ANOVA test was used for differential abundances between the two groups. With specific regard to VOCs analysis, the separation of means between the two groups was assessed by Student's t-test, and differences were considered significant for p values lower than 0.05 and 0.01.

Acquisition of Metagenomic Data and Analysis
We obtained a mean of 236,831 and 277,576 reads for the CTR and OL+ group, respectively. From 75-97% starting reads matched primers. Over 4.1 Mb non-chimeric reads, 2.2 Mb (80%) from OL integration feed, and 1.9 Mb (83%) from conventional feed were used in generating consensus sequences for analysis ( Table 2). All the taxa of the ATCC MSA-1002 were correctly classified at the genus level. Data are reported as the mean of the ten records per group ± standard deviation.

Effects of Diet on Microbial Community Composition
Alpha diversity (Shannon and Chao1) was calculated at the genus level between two feed groups OL+ and CTR ( Figure 1). Although OL+ appears to be slightly higher overall than CTR, the differences were not statistically significant (p > 0.05).
Beta-diversity analysis showed OL+ and CTR are significantly different (Adonis tests of both Bray-Curtis and Jaccard indices, p ≤ 0.001). In this study, Bray-Curtis and Jaccard indices evidenced that OL integration in feed affects the rumen microbial community. The contribution of individual genera to the overall Bray-Curtis dissimilarity is reported in Table 3. The analysis of the sequences identified a total of 157 genera in the two groups; however, the genera with a relative abundance greater than 1% were 16 for CTR and 14 for OL+ (Table S1). Percentage contribution to Bray-Curtis dissimilarity showed that 11 genera were influenced by OL integration (Figure 2). Among them, in OL+, a higher number of reads (p < 0.05) was observed in Paludibacter (259 counts vs 1061.5 counts in CTR and OL+, respectively), Fibrobacter (849.2 counts vs 1502.7 counts in CTR and OL+, respectively), Sphaerochaeta (142.3 counts vs 386.2 counts in CTR and OL+, respectively), Christensenella (76.1 counts vs 192.9 counts in CTR and OL+, respectively), Rikenella (44.4 counts vs 165.7 counts in CTR and OL+, respectively), Oligosphaera (49.7 counts vs 137 counts in CTR and OL+, respectively), Candidatus Endomicrobium (3.4 counts vs 77.6 counts in CTR and OL+, respectively), Anaerovorax (12.3 counts vs 71.3 counts in CTR and OL+, respectively), Atopobium (12.7 counts vs 40.6 counts in CTR and OL+, respectively). Conversely, OL integration reduced the percentage of Bacteroides (233.2 counts vs 67.3 counts in CTR and OL+, respectively) and Selenomonas (43 counts vs 11.6 counts in CTR and OL+, respectively).
Beta diversity was then calculated at family taxonomic rank with statistical significance slightly lower (p < 0.004) than genus-level tests. At this taxon level, the sequences could be assigned to thirty-two different families (Table 4).

Identification of Volatile Compounds
The analysis of VOCs in the rumen samples obtained from the goat kids made it possible to identify a total of 24 compounds belonging to different chemical families: acids, aldehydes, alcohols, ketones, esters, and hydrocarbons. As shown in Table 5, the OL introduction into the animal's diet showed effectiveness in inducing some variations compared to the control group. Specifically, in OL+ samples, a higher concentration of hexanoic acid (7.07 ± 0.76% vs 9.35 ± 0.89% in CTR and OL+, respectively, p < 0.05) was found, while significantly lower values were found for the decanal (1.89 ± 0.17% vs 0.85 ± 0.10% in CTR and OL+, respectively, p < 0.01). Furthermore, the finding concerning the presence of α-terpineol (3.06 ± 0.26%) only in samples obtained following the OL dietary supplementation is peculiar.  1 Volatile compounds (VOCs) are expressed as relative mean percentages (%) of total detected compounds ± standard deviation (SD). nd = not detectable; ns = not significant; * p < 0.05; ** p < 0.01.

Discussion
The rumen microbial community structure is affected by animal species, breed, age, parity, and lactation stage. However, diet is the major determinant of ruminal microbial composition, and changes in the feeding strategy can lead to rapid and dramatic changes in gut microorganisms [9][10][11]. Rumen microorganisms metabolize polysaccharides, proteins, and lipids in VFAs, microbial proteins, and vitamins that are used by animals [34]. To date, several metagenomic studies have been focused on the effect of different dietary protocols on the rumen microbiome. However, knowledge of microbial ecology is still limited and, in particular, the effects of dietary OL supplementation on the ruminal goat microbiome have not been reported in the literature.
Independently from the diet, in the present study, the main phyla identified in the rumen were Bacteroidetes, followed by Firmicutes, whereas subdominant phyla were Proteobacteria, Verrumicrobia, Fibrobacteres, and Actinobacteria. This result is consistent with previous studies published by Cremonesi et al. [35], who reported that the rumen composition of Alpine dairy goats was dominated by Bacteroidetes with about 61.2% relative abundance, followed by Firmicutes, Proteobacteria, and Verrucomicrobia at 24.2%, 4.1%, and 3.3% relative abundances, respectively. Conversely, in free-ranging Moxotó breed goats, a lower relative abundance of Bacteroidetes and a high abundance of Firmicutes was observed in the rumen microbiome [36]. This aspect, on the one hand, may depend on the different genetics of the animals involved in these two experiments; on the other, it is certainly influenced by the diet of the animals, supporting the fact that the presence of different substrates at the rumen level is effective in changing the microbiota composition.
Members of phylum Bacteroidetes are the primary degraders of complex polysaccharides in the plant cell wall, and they have higher mean glycoside hydrolase (GH) enzymes and polysaccharide lipases (PLs) genes per genome as well as signal peptide-containing GH and PLs compared with the members of the phylum Firmicutes which have a lower ability for polysaccharide degradation and seem to possess more cellulose hydrolysis capacity. Moreover, they are known for their production of butyric acid or any other bacterial phyla in the gut intestinal tract [37,38].
The data obtained in the present study showed a low (0.4:1) Firmicutes/Bacteroidetes ratio in both groups. Similarly, in Vietnamese goat rumen, a dominance of Bacteroides compared to Firmicutes and a low Firmicutes/Bacteroidetes ratio correlated to a high diversity of cellulolytic GH enzymes originating from Bacteroidetes was observed [39]. The ratio between Firmicutes and Bacteroidetes is influenced by many factors such as breed, age, and diet [38]. During the life cycle of the goat, this ratio gradually decreases, changing from 2.1:1 at 80 days old to 1.7:1 at 100 days old and to 0.3:1 at 110 days old [40]. This decrease has been correlated to the adaptation of rumen bacteria to a plant diet. Recently, in goats fed a high-grain diet, it has been observed an increase in ruminal acidity and a very high Firmicutes/Bacteroidetes ratio (~3:1), which are thought to be unhealthy changes [41]. The data obtained in this study suggest that the percentage of OL integration did not interfere with the normal rumen ecosystem.
Regarding the taxonomy at the family level, in our samples, the highest percentage contributions to Bray-Curtis dissimilarity were Prevotellaceae, Subdivision5_genera_incertae_sedis, Succinivibrionaceae, and Ruminococcaceae. However, OL integration induced significant increases in different families: Fibrobacteraceae, Christensenellaceae, Coriobacteriaceae, Oligosphaeraceae, Candidatus Endomicrobium, and Planctomycetaceae; conversely, Succinovibrionaceae and Bifidobacteriaceae decreased. Among them, even if not significant, a higher abundance of Ruminococcaceae was detected in OL+ with respect to CTR. Ruminococcaceae are fibrolytic bacteria that may degrade numerous polysaccharides such as starch, cellulose, and xylan, and produce short-chain fatty acids (SCFAs). Moreover, they play a fundamental role in the biohydrogenation (BH) of dietary unsaturated fatty acids (UFAs) [42,43], but at the same time, they are susceptible to this type of fatty acids (FAs) [44]. A number of observations suggest that Ruminococcaceae and Christensenellaceae are potentially beneficial bacteria because they participate in the positive regulation of the intestinal environment and are linked to immunomodulation and healthy homeostasis [45].
In dairy goat, Fibrobacteriaceae has been correlated with FAs involved in the BH pathway of α-linolenic acid (C18:3n3-ALA) [35]. The increase in Fibrobacteriaceae in the OL+ group may be due to a higher intake of C18:3n3 which represents the major FA of OL [46]. In the rumen, different bacteria families are responsible for the BH of dietary UFAs. The BH process identified two important microbial transformations of fats in the rumen: lipolysis and BH. Bacteria belonging to the Acidaminococcaceae family, such as Anaerovibrio lipolytica, hydrolyze the ester bond of tri-and di-glycerides; on the contrary, Butyrivibrio (Lachnospiraceae family) lipases hydrolyze phospholipids. The Butyrivibrio group, including the genus Butyrivibrio and Pseudobutyrivibrio, and the species B. proteoclasticus, are able to hydrogenate C18:3 and form cis9, trans11-conjugated linoleic acid (C18:2, c9, t11-CLA) and vaccenic acid (C18:1, t11-VA) [47][48][49]. Previous studies carried out on different species of ruminants have reported that diets supplemented with agroindustrial by-products containing considerable amounts of plant secondary metabolites such as different kinds of polyphenols, saponins, and essential oils, are able to modulate rumen BH [50,51]. However, the effects on the rumen microbial population differ in relation to several factors such as the rumen passage rate, interaction with basal diet, different amounts of lipids in the diet, and the specific composition of phenolic substances. In this trial, no significant differences between the two groups were observed in Lachnospirace and Acidaminococcacea families, and the relative abundance of Butyrivibrio, Pseudobutirivibrio, and Anaerovibrio genera did not change, suggesting that OL dietary integration was not able to affect ruminal BH. However, further analysis is required to better understand the possible influence of OL intake on lipid metabolism.
At the genera level, the results showed a higher abundance of Paludibacter, Fibrobacter, Sphaerochaeta Christensenella, Rikenella, Oligosphaera, Candidatus Endomicrobium, Anaerovorax, and Atopobium in the rumen of goat kids fed with OL; on the contrary, Bacteroides and Selenomonas were negatively influenced. In both groups, Prevotella (phylum Bacteroidetes) was the most dominant genus, in accordance with studies in goats and in other ruminants [52,53]. Prevotella comprises many species involved in the degradation of different substrates such as hemicelluloses, pectin, starch, protein, and simple sugars as energy sources to produce succinate as the major fermentation end product. In the present study, no significant differences between the two groups were observed in the abundance of Prevotella. These data were in accordance with previous studies where it has been demonstrated that the abundance of Prevotella did not always change in relation to different diets. Liu et al. [53] observed a higher abundance, although not significant, of Prevotella in goats fed a complete feed diet with respect to goats fed an all-forage diet. Bekele et al. [52] showed an increasing trend of Prevotella in the concentrate-fed group compared to the roughage group in the rumen of sheep, although the abundance did not differ significantly in the various ration groups. On the contrary, Huo et al. [54] found that the abundance of Prevotella in the rumen from hay-fed goats was higher than in concentrate-fed animals.
Among the main fibrolytic bacteria (Fibrobacter, Ruminoccoccus, Butyrivibrio, and Clostridium spp.), significant variations were observed only in Fibrobacter, even if Ruminococcus spp. and Clostridium spp. were present at a higher relative abundance in the OL+ group with respect to the CTR. Fibrobacter is the major cellulose-degrading bacterial species in the rumen of herbivorous animals and is present in high-fiber diets [55]. Compared with other fibrolytic species such as Ruminococcus, it digests fiber faster. Different studies have demonstrated that in animals fed a tannin-rich diet, the fiber-degrading bacteria, Fibrobacter and Ruminococcus, were reduced [56]. In this study, a specific evaluation of tannins content in the used OL was not performed; however, this matrix was demonstrated to be relatively rich in these compounds [57]; therefore, it must be taken into account that the method of administration, the proportion of OL in the diet, and the time of the trial used in our study could have affected the observed effects on fibrolytic bacteria.
Recently, it has been highlighted that bioactive plant compounds such as tannin and polyphenols can modify rumen methane production. Moreover, in vitro studies have shown a reduction in methane production in the rumen combined with different parts of olive [19]. Many methanogens belong to the phylum Euryarcheota (Archea domain) and the most common are from the genera Methanobrevibacter, Methanimicrococcus, and Methanobacterium, which convert H 2 and CO 2 into methane. Less abundant are Methylotrophs, Methanosarcinales, Methanosphaera, and Methanomassiliicoccaceae which produce methane using methylamines and methanol [11].
In our study, methanogens were scarcely represented, accounting for only 9.9 of Euryarcheota and 14.1 of the number of observed OTUs in CTR and OL+, respectively, with a relative abundance of about 0.01% in both groups. At deeper taxonomy levels, among methanogens, three different genera were detected in both groups: Methanimicroccus (Methanosarcinaceae family), Methanobrevibacter (Methanobacteriaces family), and Methanomassiliicoccus (Methanomassiloliicoccaceae family); however, their relative abundance was lower than 0.01%.
VOCs are produced by the metabolism of microorganisms that are able to ferment feed components. The diet can directly affect the rumen microbiota, and subsequently the production of VOCs. The specific composition of VOCs produced from a given bacterial community will depend on the diversity and metabolic types/groups of species present. The detection of these compounds in the rumen is a useful tool to evaluate physiological changes due to animal nutrition and numerous analytical approaches have been developed to accurately characterize and measure VOCs.
In the present study, the headspace sampling with solid-phase microextraction coupled with gas chromatography-mass spectrometry was an effective approach that allowed to highlight 24 compounds belonging to different chemical families: acids, aldehydes, alcohols, ketones, esters, and hydrocarbons. First of all, the dietary supplementation seems to have induced an increase in hexanoic acid and a decrease in decanal at the rumen level. This condition assumes particular relevance when compared with what was previously reported by Spinhirne et al. [58], who exploited the solid-phase microextraction followed by gas chromatography-mass spectrometry for the analysis of volatile organic compounds in bovine breath. In particular, it was highlighted that decanal was more associated with clinically morbid steers, while the presence of hexanoic acid was symptomatic of healthy steers. As far as our study is concerned, it is therefore plausible that dietary supplementation with OL has somehow had positive effects on animal health. In support of what has been reported, the presence only in OL+ samples of α-terpineol, a monocyclic monoterpene tertiary alcohol naturally present in vegetable species, should also be highlighted. In fact, this compound is associated with a wide range of interesting biological properties for its antioxidant, anticancer, anticonvulsant, antiulcer, and antihypertensive potential; this compound was even credited with insecticidal properties [59].
Overall, what has been reported therefore lays the basis for hypothesizing an effect of dietary OL supplementation in improving the welfare conditions of the animals involved in the experimentation. In fact, the reduction in the straight-chain aldehydes derived from the oxidative degradation of UFAs is generally associated in mammals with a condition of reduced oxidative stress. Furthermore, the concomitance with the accumulation of hexanoic acid confirms this aspect as it indicates the predominance of triglycerides degradation following enzymatic events [60]. However, from this point of view, it is necessary to conduct further and more in-depth evaluations, which can clarify the metabolomic aspects associated with these findings, also in light of the changes observed at the metagenomic level.

Conclusions
Our data suggest the capability of OL to modulate the ruminal microbial community in goat kids. The OL addition seemed to improve fiber degradation modulating cellulolytic bacteria. The changes observed at the metagenomic level were associated with a different rumen volatile profile, and the presence of VOCs associated with health status suggests a favorable role of OL in preserving and improving animal welfare. Further investigations are needed to improve knowledge of the functional genomic framework of olive leaves degradation in the rumen and to examine their impact on lipid and biohydrogenation metabolism.
Supplementary Materials: The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/vetsci9090452/s1, File S1: Beta diversity; Table S1: Relative abundance (%) of the main genera identified in the control group (CTR) and experimental group (OL+).  Institutional Review Board Statement: Ethical review and approval were waived for this study due to the fact that this research used animals and data from commercial farms which were handled following the national legislation on animal welfare (DL n.126, 07/07/2011, EC Directive 2008/119/EC), and then slaughtered complying with the EU Regulation 1099/2009 on the protection of animals at the time of killing. The ruminal sampling was carried out following the killing of the animals, therefore were not introduced practices different than those normally envisaged.
Informed Consent Statement: Written consent was obtained from all the animal owners enrolled in the study.
Data Availability Statement: All data are available following reasonable request to the corresponding author.