Assessment of Biolog EcoplateTM method for functional metabolic diversity of aerotolerant pig fecal microbiota

Abstract In the last decades, gut microbiota and its role in mammal host development and health have been increasingly investigated. Metabolites produced by gut microbiota can affect intestinal homeostasis and immune system maturity and activation, and in turn, they can influence the health and growth performance of livestock. Therefore, a better understanding of the functional metabolic capability of the gut microbiota would be appreciated by the scientific community. In this study, the BiologTM Ecoplates technology was applied for studying the metabolic potential of the aerotolerant microbial community of pig fecal samples, evaluating the interference of different storage conditions and cell concentrations. The length of time for which a fecal sample maintained detectable and unchanged microbial metabolic activity was also investigated. Two assays aimed to evaluate differences in the metabolic activities between fresh and snap-frozen fecal samples at different dilutions and at different lengths of times of preservation at −80°C were carried out. The biodiversity and the predicted functionality of the entire bacterial community through a targeted metagenomic approach were also explored. The results highlighted that snap freezing of fecal samples preserved the metabolic activity of the microbial community when compared to fresh feces. Sample storage at −80 °C did not significantly affect the metabolic activity of the microbial community, which was stable for 150 days. Furthermore, the highest metabolic activity was detected with 1:2 to 1:5 dilutions of the stock suspension. BiologTM Ecoplates technology is a rapid and useful method to explore microbial communities’ metabolism in animal fecal samples contributing to investigate host animal physiology. Key points • Freezing of samples can preserve the functional activity of the aerotolerant microbial community for 150 days. • The concentration of microbial cells strongly influences metabolic activity detection. • Sequencing coupled with the BiologTM Ecoplates could be a strategy to evaluate the metabolic potential of the microbiota of the fecal sample. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1007/s00253-021-11449-x.


Introduction
Modern livestock production is based on the improvement of animal efficiency coupled with a reduction of antimicrobial use as keystones for more sustainable production of proteins of animal origin. The gut microbiota plays a key role to achieve these objectives. Indeed, the importance of the interplay between animal and gut microbiota in the regulation of the animals' health has been widely highlighted (Putignani et al. 2014). Generally, microbial communities can contribute to the host homeostasis, affecting the host metabolism (Tremaroli and Bäckhed 2012), influencing processes including nutrient digestion, nutrient synthesis (such as vitamins and fatty acids), and gut epithelial cell turnover (Yeh et al. 2019) as well as influencing the development and activation of the immune system (Trevisi et al. 2018). In monogastric animals, microbial communities have been demonstrated to contribute to the improvement of pig feed efficiency and consequently its productive performances (Kiros et al. 2019).
A. Checcucci and D. Luise contributed equally to this work.
Besides age and host immune competence, diet composition is one of the main factors modulating gut microbial profile and metabolism (David et al. 2014). The gut microbial community, with a huge and diverse concentration of bacteria, viruses, and fungi, contributes to the degradation of feed components during the passage of digesta along the gut intestinal tracts (GIT). Among nutrients, carbon sources and, in particular, complex carbohydrates are largely metabolized by gut microbial communities obtaining oligosaccharides, which are then fermented to short-chain fatty acids (Puertollano et al. 2014). Dietary and endogenous amino acids (AAs) can also be used by the gut microbiota (Dai et al. 2010). In the small intestine, peptides can be metabolized by bacteria resulting in the production of polyamine involved in cell metabolism, growth, and proliferation. In contrast, in the large intestine, AAs can be fermented by the resident microbiota producing detrimental metabolites such as ammonia, indole, and skatole (Bardócz et al. 1995).
Presently, several techniques have been used to investigate the animal gut microbial community composition, in particular high throughput sequencing coupled with bioinformatics analysis. Nevertheless, a proper analysis of the gut microbial community should include information regarding metabolic activity. Microbial metabolism depends on several factors, such as active population in terms of species and quantities, substrate composition, and environmental conditions; thus, meta-transcriptomics and metabolomics have been exploited to investigate the microbial community functions. An additional method to study the metabolism of entire microbial communities is represented by Biolog TM 's platform coupled with Ecoplates (Biolog Inc., Hayward, CA, USA), which is based on the detection and simultaneous measurement of carbon substrate utilization level from the microbial community (Stefanowicz 2006). The method has been widely used for the study of microbial metabolism of environmental sources, including soils and sediments (Gałązka et al. 2017;Lopes et al. 2016), water plants, and grains (Ge et al. 2018). Moreover, Biolog TM 's Ecoplates technique has been applied to evaluate indirectly the functional traits of gut microbial communities of livestock, manure composting, and wastewater (Wang et al. 2018;Gryta et al. 2014). More recently, Phenotype Microarray TM was applied to study the gut and fecal microbial communities of several animal species (Grześkowiak et al. 2020;Yeh et al. 2019), including pigs (Piotrowska et al. 2014), outlining their metabolic profile in certain physiological or environmental conditions (Najdegerami et al. 2012;Pierce et al. 2016).
However, presently, there is no information regarding the effect of the conservation including freezing and length of storage of these complex matrices for their subsequent analysis through the Biolog TM Ecoplates technique. This information could improve and broaden the use of this technique for the evaluation of the microbial metabolic fingerprint in pigs' gut. In fact, in the normal day-by-day sample collection, often samples are maintained frozen before being processed. In this methodological study, the authors investigated, for the first time, the metabolic potential and the functional capacities of the microbiota of pig fecal samples using the Biolog TM platform coupled with Ecoplates system with the following aims: (i) identify the optimal microbial cell concentration for metabolic function description; (ii) compare metabolic activity between fresh and frozen samples; (iii) verify that long storage time (5 months) are not significantly impacting on the correct metabolic analysis of frozen samples. Furthermore, the authors suggest the taxonomical analysis of microbiota as a method that can be coupled with the characterization of microbial community metabolic activity. Indeed, such determination allows for the direct assignment of the metabolic phenotype of the community to specific bacterial groups, which are found to be the most part in pig gut microbiota.

Materials and methods
Animal housing, sample collection, preparation, and studies design Healthy post-weaning piglets (weaned at 28 days of age) were included in the experiment and housed in pens at the University of Bologna. Room temperature was kept controlled at 30°C at the beginning and 25°C at the end of the experiment, with a 1°C decrease every 3 days. Infrared lamps were located above the piglets for the first 7 days post-weaning. Piglets were fed a standard commercial diet without antimicrobials and based on barley, wheat, corn, milk whey powder, and soy protein concentrate (diet composition: 17.9% crude protein and 5.0% crude fibers). Piglets had free access to feed and water throughout the experimental period; feed was ad libitum supplied in a dry feeder.
On day 35 of age, fresh feces samples were collected into a sterile collection tube of 120 mL immediately after defecation from healthy piglets. Then, to limit the environmental contamination, approximately 5g of feces was taken from the non-superficial part of the feces, homogenized, and divided into aliquots of 500 mg each. For experiment 1, two aliquots were prepared from one homogenized fecal sample: the first aliquot was stored at room temperature for 4 h (fresh sample) and the second aliquot was immediately frozen in liquid nitrogen for 5 min and then stored at −80°C (frozen sample). Both samples were analyzed after 4 h. For experiment 2, one homogenized fecal sample was divided into four aliquots that were immediately frozen in liquid nitrogen for 5 min and then stored at −80°C.

Experiment 1: comparison between fresh and frozen feces
Fresh and frozen samples (500 mg each) were first suspended in 50 mL sterile PBS (pH 7.2) and vortexed: these suspensions were named stock suspensions (SSUSPs). SSUSPs were shaken for 30 min at 21°C and then allowed to rest for 15 min at 4°C , according to Ecoplates' sample preparation protocol used by Gryta et al. (2014) and Núñez-Díaz et al. (2017). For each SSUSP, two different dilutions were prepared: 1:10 (approximately, 9 × 10 4 bacterial cells/g) and 1:100 (approximately, 9 × 10 3 bacterial cells/g). The bacterial cell count of SSUSP samples was performed on Luria Bertani (LB) plates, after an incubation of 24h at 39°C. The SSUSP and the diluted samples were then analyzed using Biolog TM Ecoplates. Dilution showing the highest metabolic activity was the basis for choosing the dilutions of experiment 2.

Biolog TM phenotype assays
The carbon source utilization profile of pig fecal microbiota was assessed by using the fully automated Biolog Omnilog TM (Biolog, Hayward, CA, USA) System through Ecoplates Phenotype MicroArray (PM). These 96-well microplates contain three replicate sets of 31 lyophilized relevant C substrates, along with a tetrazolium redox dye (Insam 1997). C substrates were grouped by chemical class (eight carbohydrates, eight carboxylic acids, four polymers, six amino acids, two amines, and three miscellaneous substrates; Table 1). Substrate utilization by the microbial community is evidenced by a color change in each plate well, as the redox dye is reduced to tetrazolium violet. Redox technology used by the PM system measures cell energy (NADH) production following substrate utilization as a universal reporter. More specifically, if the microbial cells have an active metabolism, there would be a flow of electrons to NADH, which determines the reduction of a tetrazolium dye and the consequent production of a purple color.
For the study design, samples were diluted according to the experimental design previously described and 150 μl of sample supernatant was dispensed into each well of the Ecoplates (Piotrowska et al. 2014).
The Ecoplates were then incubated at 39°C in aerobic conditions in the dark. Substrate's consumption was recorded by the automated Omnilog TM System every 15 min for 96 h, to allow reaching plateau phase when the substrates are degraded. Microbial activity in each microplate well was expressed as metabolic activity value (MAV), which is the Biolog TM index based on the ability to use different compounds by a bacterial community.

Microbial community profiling
Total bacterial DNA for microbiota profile analysis was extracted from the frozen fecal samples used in experiment 2 by using FastDNA TM Spin Kit for Soil (MP Biomedicals Europe, LLC). DNA quantity and quality were evaluated using NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE, USA). The V3-V4 hypervariable regions of the 16S rRNA gene were then amplified   and dilution were included as factors. Data of experiment 2 were fitted using a linear model in which the time of storage, dilution, and their interaction were included as factors. ANOVA and Tukey's test analysis were then carried out. The packages "car" (version 3.0.3), "lsmeans" (version 2.30.0), and "multcomp" (version 1.4.10) within the R software (R Core Team, 2019) were used to carry out the statistical analysis and "ggplot2" (version 3.3.0), "plyr" version (1.8.6), "tidyr" (version1.0.2) packages were used to create figures. Results were considered significant at P ≤ 0.05 and tendencies at 0.05 ≤ P ≤ 0.10.

Next-generation sequencing analysis
Microbial sequence data were analyzed using the DADA2 pipeline (Callahan et al. 2016), and taxonomic categories were assigned by using the Silva database (release 132) as reference (Quast et al. 2012). To explore the metabolic potential of the microbial communities, a functional prediction based on the 16S rRNA gene sequences was performed using Tax4Fun (Aßhauer et al. 2015) that provides a good approximation to functional profiles obtained from metagenomic shotgun sequencing approaches. Tax4fun produces a KEGG Orthology (KO) abundance table, whereas KO identifies an ortholog of an experimentally characterized gene/protein. Subsequently, the KO enrichment analysis was carried out using Microbiome Analyst online tool (Dhariwal, et al. 2017). Enrichment analysis was performed via over-representation analysis (ORA), using the hypergeometric test to test if a particular group of KOs are represented more than expected by chance within the uploaded list of genes. P values were corrected for multiple comparisons using Benjamini and Hochberg's false discovery rate (FDR). Results were considered significant at FDR < 0.05.

Results
Experiment 1: effect of sample storage and cell concentration on functional metabolic activity The effect of the storage conditions (fresh sample vs frozen sample), as well as of the different bacterial cells' concentration (9 × 10 5 cells/g (SSUSP sample), 9 x 10 4 cells/g (1:10 dilution), and 9 × 10 3 cells/g (1:100 dilution)) on the functional metabolic activity (expressed as MAVs) of fecal samples was evaluated.
Supplemental Fig. S1 shows the functional metabolic diversity of fresh and frozen samples at each dilution tested for all the time points during the 96-h incubation time on each carbon source category. Results highlighted the utilization of the complex carbon sources, carboxylic acids, amino acids, and amines in both fresh and frozen sample dilutions except for dilution 1:100 in frozen condition, in which MAVs remained close to zero. Generally, lower MAVs were observed for the fresh sample diluted at 1:100 dilution compared with SSUSP and 1:10 dilution. In wells filled with carbohydrates, complex carbon sources, and carboxylic acids, the MAVs reached by the fresh sample at the dilution 1:10 at 96 h were numerically higher in respect to the MAVs of the frozen sample at the same dilutions. In wells filled with amines and amino acids, a different trend was observed: the SSUSP sample reached the highest MAVs already after 25 h of incubation (Supplemental Fig. S1). Table 2 shows the mean of MAVs for the different carbon source classes and compounds and the effect of sample conservation (fresh vs frozen) and sample dilution after 96 h of incubation time. At the end of incubation, the results indicated that, besides the sample storage and bacterial cell concentration, the microbial community in fecal samples could not utilize or poorly utilized (< 40 MAVs) the following substrates as a carbon source: erythritol (except for fresh 1:10 sample), 2-hydroxy benzoic acid, 4-hydroxy benzoic acid, itaconic acid, and keto butyric acid. On the other hand, the sources that had higher values were pyruvic acid methyl ester, D-xylose, D-galacturonic acid, L-asparagine, N-acetyl-D-glucosamine, glycogen, D-cellobiose, D-xylose, and D-lactose ( Table 2).
Experiment 2: preservation of the functional metabolic activity in samples with different bacterial cell concentrations and time of storage Experiment 2 allowed to verify the impact on metabolic activity of low temperature storage at different time points (T1 = 1 day; T2 = 5 days; T3 = 45 days; T4 = 150 days) using three different bacterial concentration (SSUSP and dilutions 1:2 and 1:5). Figure 1 and Supplemental Fig. S2A and S2B show the functional metabolic diversity (represented as MAVs), during the incubation time, of 1:2 and 1:5 dilution of SSUSP sample, respectively. The storage time (T1, T2, T3, and T4) for every carbon source was evaluated. The results highlighted the utilization of carbon source in all dilutions and storage time of the samples (Fig. 1). At the end of 96 h of incubation, for both the sample dilutions and timing of storage, the MAVs of carbohydrates and complex carbon sources were the highest among the biochemical categories while those related to amines were the lowest (Fig. 2). Table 3 reports the MAVs' mean of the different carbon sources after 96 h of incubation, the effect of sample dilution, length of storage, and the interaction between dilution and length of storage factors.
To investigate the microbial functional diversity, the sequences obtained from 16S rRNA sequencing were annotated to the KEGG database using Tax4fun. The KO table produced was used to perform the enrichment analysis, to identify pathways that are represented more than expected by chance. The results of the enrichment analysis are reported in Supplemental Table S1. The most enriched pathway resulted in the biosynthesis of amino acids, carbon metabolism, methane metabolism, porphyrin and chlorophyll metabolism, carbon fixation pathways, pyruvate, propanoate, glycine, serine and threonine metabolism, peptidoglycan biosynthesis, and glycolysis/gluconeogenesis (P FDR < 0.01).

Discussion
Compared to meta-transcriptomics, metabolomics, and other tools generally applied to investigate the microbial community functions, the Phenotype Microarray TM is a relatively fast and cost-effective method to examine microbial metabolic functions and phenotypes. In particular, the Ecoplates technique enables the study of whole microbial communities and investigating their carbon substrate utilization, outlining a community metabolic fingerprinting.
The present study evaluated the inoculum density and the storage conditions needed to preserve a detectable metabolic activity of the microbial community present in piglet fecal specimens.
A significant effect of the inoculum density on the consumption of the carbon substrates was observed. Indeed, considering the results of experiment 1, only the dilutions of the original fecal sample, which contained about 9 × 10 5 bacterial cells/g and the dilution containing 9 × 10 4 bacterial cells/g showed a detectable functional activity. Furthermore, differently diluted samples showed differences in the metabolic performances for carboxylic acids (higher in 9 × 10 4 bacterial cells/g dilution in respect to SSUSP sample) and amines consumption (lower in 9 × 10 4 bacterial cells/g dilution in respect to SSUSP sample). These results suggested that not only the composition of the microbial population, as already described by Szabó et al. (2007), but also the total cell concentration of the microbial community and the growth condition in the presence of O 2 can affect its general functions and metabolism (for example, its ability to degrade organic matter). Indeed, it is possible that a large number of cells (bacterial, but not only) can generate a premature saturation of the Ecoplates system, rapidly consuming the substrates in the wells. This may generate a pH variation counterproductive to microbial growth and the Biolog TM staining measurement (Biolog, personal communication). Therefore, it is usual that undiluted inoculum shows the worst growth performances than diluted ones. Microbiologically, dilution cultures have proven advantageous in the study of oligotrophic bacteria, since inoculum dilution maintains numerically abundant organisms while eliminating fast-growing competitors that typically prevail in enrichment cultures (De Fede et al. 2001).
A potential more uniform response for all the categories of compounds was found with an inoculum density of 9 × 10 5 (SSUSP) and 9 × 10 4 (1:10) cells/g; however, the dilution 1:10 showed a lower MAV trend for amines and amino acids compared to SSUSP samples that in the other side had a lower trend for carbohydrates; therefore, two intermediate dilutions T1, 1 day after sampling; T2, 15 days after sampling; T3, 45 days after sampling; T4, 150 days after sampling 1:5 (1.8 × 10 5 cells/g) and 1:2 (4.5 × 10 5 cells/g) were utilized in experiment 2 to optimize the metabolic activity detection.
Recently, the study of Metzler-Zebeli et al. (2016) showed that, in pig feces, the freezing process can influence the DNA yield and bacterial abundances, while the snap freezing and storage temperature had only little effect on the bacterial abundance. Therefore, this study evaluated the type of storage and the preservation of the bacterial community activity for 5-month storage at −80°C.
Considering the optimal dilutions of the samples, the results showed no difference in the metabolic activities of the pig fecal microbial community when the same inoculum density was compared in both fresh and frozen samples. Therefore, the immediate freezing and storage at −80°C of pig fecal samples did not affect the microbial metabolic potential and substrates utilization after 96 h of incubation. Results also suggested that sample freezing, when performed immediately after the sample collection, can ensure a robust and error-minimizing analysis of the microbial metabolism. Therefore, there is no need to process a sample immediately after collection, as reported in several studies (Piotrowska et al. 2014;Yeh et al. 2019), resulting in higher analysis flexibility and a broader application of this technique.
Furthermore, the results obtained in experiment 2 confirmed that the sample storage at −80°C for 150 days did not affect the substrate consumption by the microbial community. The utilization rate of carbon sources showed a similar trend during the incubation time at all tested microbial concentrations, and it was not influenced by the time of storage of the samples. Also, when considering the single compound utilization after 96 h of incubation, overall, the time of storage had no statistically significant effect.
Overall, in agreement with Grześkowiak et al. (2020) and Kong et al. (2013), carbohydrates were the most consumed substrate category, while the least used carbon sources were the ammines, that vary depending on the microbial community tested (Grześkowiak et al. 2020;Kong et al. 2013). The utilization of carbohydrate sources by the pig fecal microbial community was confirmed by the enrichment analysis on the list of KEGG genes predicted from the taxonomic composition of the microbial community of our sample. Indeed, our results evidenced that the major enriched pathways were related to carbohydrate. According to our Biolog TM results, the most consumed carbohydrates, on which the MAVs reached after 96 h of incubation remained almost unaltered in all dilution tested, were glucose-1-phosphate, α-D-lactose, D-xylose, D-cellobiose, and N-acetyl-D-glucosamine. In accordance with the observation of Grześkowiak et al. (2020), our data confirmed that the primary substrate for the pig gut microbiota is glucose in its different forms. Furthermore, the always elevated consumption of lactose can be due to the high abundance of Lactobacillus in the fecal microbial community of our sample. This genus, which is typical of the gut microbiota of young pigs, tends to increase according to the quantity of dietary lactose that is high in the diet of young piglets (Daly et al. 2014). Accordingly, Lactobacillus abundance undergoes a strong reduction in growing and adult pigs (Han et al. 2018). The taxonomic composition of the microbial community observed in the present study agrees with the common fecal bacterial profile of the weaned pigs (Chen et al., 2017;Motta et al., 2019). Interestingly, the high consumption of D-cellobiose should allow hypothesizing unaltered metabolic activity of cellulolytic bacteria (Gupta et al. 2012), such as those belonging to the family of Ruminococcaceae, classically associated with ruminants' and non-ruminants' gut, such as horses and pigs, largely present in our samples. Another highly abundant genus in our fecal samples was Prevotella. Prevotella abundancy in the gut of post-weaned pigs is generally connected to the intake of starch and fiber-rich diet which increases in the  (Ivarsson et al. 2014). However, its abundance is positively correlated with the content of xylose of the pigs' diet (Ivarsson et al. 2014) and it is known that some species of Prevotella can utilize monosaccharides including glucose (Takahashi and Yamada 2000).
Apart from carbohydrates, the modern diets for piglets have a concentration of highly digestible crude protein of around 17% and include free amino acids to satisfy piglet nutrient requirements. Human and animal gut microbes may extensively utilize amino acids for the synthesis of proteins in a species-dependent manner. In our study, L-asparagine and L-serine were the most used amino acids. Previous studies related to the metabolism of the amino acids in microbes in the pig intestinal content demonstrated that asparagine, together with aspartate, is one of the most rapidly fermented amino acids by ileal and jejunal mixed bacteria (Dai et al. 2010). This can explain the high consumption rate of asparagine by the fecal microbes of our samples. The results of the prediction of the pig feces microbiome functions obtained from 16S rRNA sequencing revealed that the pathway of serine and threonine metabolism was one of the most represented pathways identified by Tax4fun, which confirmed recent data on the functions of the pigs' cecum and ileum microbiome (Umu et al. 2020). In supporting this, it is known that in the large intestine of pigs a large variety of bacteria can use mucus and endogenous sources that are rich in threonine and serine (Le Floc'h et al. 2018). That information can explain the always high utilization of L-serine in all the samples and dilutions observed in the present study.
Our results showed a relatively low consumption of amines as carbon sources, which (i) further decreased at lower bacterial cell concentration (1:100 dilution) and (ii) numerically decreased with sample storage (from T1 to T4 in dilution 1:5). As for amines, the Ecoplates includes phenylethylamine and putrescine. Phenylethylamine is the decarboxylation product of phenylalanine, and it is catalyzed by the enzyme tyrosine decarboxylase (Kong et al. 2013) as a result of microbial fermentation. Recently, it has been investigated as a potential antimicrobial compound in cross-contamination of food products, influencing bacterial growth and biofilm formation through the interference with bacterial signal transduction pathways (Lynnes et al. 2014). Therefore, we can hypothesize that a similar inhibitory effect of phenylethylamine was exerted on piglets' fecal microbiota, in particular at the lower bacterial cell concentration. Putrescine is one of the main polyamines formed in the intestinal tract from bacteria by decarboxylation of lysine and ornithine. The pathways for bacterial degradation of putrescine are related to those of γ-aminobutyric acid (GABA), arginine, and ornithine. The pathways for putrescine degradation are transaminase pathways and the glutamylation pathways (Kurihara et al. 2008), which can degrade putrescine as a sole carbon source at ≤30°C , but not at higher temperatures, as confirmed by the results of our experiments where the incubation temperature was set at 39°C to mimic the pig internal body temperature. As previously reported, the less consumed carbon sources vary depending on the microbial community. In our study, low utilization of amines (phenylethylamine and putrescine) by the piglet fecal microbial community was observed. The main pathways of polyamine catabolism are the oxidative deamination by Cu 2+ diamine oxidase enzyme (DAO,EC 1.4.3.22) and by the monoamine oxidase enzyme (N-acetylputrescine oxidase, MAO, EC 1.4.3.4). Only a few strains of lactic acid bacteria (LAB) have been proven to degrade amines (Callejón et al. 2014) and in particularly stringent pH conditions. Thus, this may explain the results obtained in the metabolic profile analysis of our samples. Indeed, we tested fecal samples of 35-day-old pigs where the microbiota profile highlighted a high relative abundance of genera belonging to the lactic acid-producing bacteria.
The present study is the first pilot methodological study aimed to characterize the phenotype of pig fecal samples' microbial community. The results were comparable with the metagenomics data. Nevertheless, more additional aspects need to be considered. Firstly, in addition to bacteria, the microbial community of pig fecal samples is characterized by other microorganisms, such as yeasts and fungi (Ramayo-Caldas et al., 2020). This may have led to discrepancies between the experimental results of functional metabolism detected by the Ecoplates and the prediction analysis on 16S rRNA performed on the bacterial microbiota. Secondly, the aerobic incubation of sample cultures during Ecoplates experiments may have excluded the strict anaerobic microorganisms, which generally are represented in high percentage in pigs' fecal samples and whose metabolism plays a fundamental role among gut microbes (Crespo-Piazuelo et al., 2018). Thus, further Ecoplates TM experiments on fecal samples under anaerobic conditions must integrate these results and are desirable to deeper investigate the metabolic activity of microbial community simulating intestinal conditions of pigs.
In conclusion, the results obtained in this study evidenced that the functional metabolic activity of the aerotolerant microbes in pig's fecal samples can be preserved without statistically significant variation until 150 days after a snap freezing and storage at −80°C. Therefore, the sample's microbial community remained representative within 5 months from the sampling. Furthermore, the bacterial cell concentration, if between 4.5 × 10 5 and 1.8 × 10 5 bacterial cells/g, does not influence the microbial community metabolism in different carbon sources tested. These results suggested high flexibility of the analysis and broader application of this technique that can represent a rapid and representative method to explore the metabolic capabilities of the microbial community in animal samples which, coupled with the investigation of species distribution in the samples, could contribute to explain the health status and the physiology of the animals. Furthermore, the utilization of prediction methods for the functional gene family's abundance present in the fecal microbial communities through the taxonomical characterization of bacterial community is of importance for the prediction and the comparison of the functional potential of fecal samples belonging to different individuals. The application of Biolog Ecoplates coupled with sequencing prediction method may contribute to deeply understand the interplay between gut microbiota and host, helping in moving from a sterile description of the gut microbiota to a more accurate exploitation of its potential in livestock animals.
Nevertheless, the combination of the PM technique with 16S rRNA sequencing data can present some limitations (neither can represent the totality of the metabolic capabilities of feces microbial community), it is necessary to consider that no condition (dilutions, states of conservation, temperatures, oxygenation conditions) is suitable for all microbes present in the samples. Therefore, the connection reported in this study can still represent the first attempt for a better description of gut microbiota, contributing to give a further perspective regarding the key role of gut microbiota in pigs.
In the future, for a complete assessment of metabolic activity, it will be necessary to perform Phenotype Microarray TM with Ecoplates setting up both aerobic and anaerobiobic condition, coupled with a targeted and untargeted metagenomics sequencing for the taxonomical and functional characterization of such entire fecal microbial community.
Author contribution AC: investigation, formal analysis, writing -original draft, writing -review and editing; DL: formal analysis, data curation, writing -original draft, writing -review and editing; MM: formal analysis, writing -original draft; CF: formal analysis, writing -original draft; BP: investigation, conceptualization; MP: investigation, conceptualization, project administration, writing -review and editing; PT: investigation, conceptualization, project administration, writing -review and editing Funding Open access funding provided by Alma Mater Studiorum -Università di Bologna within the CRUI-CARE Agreement.
Data availability All data generated or analyzed during this study are included in this published article (and its supplementary materials file).
Code availability The datasets on microbiota data supporting the conclusions of this article are available in the PiGutNet area of the European Nucleotide Database (ENA) repository with the number PRJEB40788.

Declarations
Ethics approval Fecal samples were collected after defecation. No direct procedures were specifically carried out on the animals for this study, and therefore, this work did not require ethical approval under Italian law (Decreto Legislativo 4 March 2014, n. 26).

Conflict of interest The authors declare no competing interests.
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/.