1. Introduction
Pet food production was 34.96 million tons all over the world in 2023, and it was still increasing while most of other animal feed was decreased [
1]. Companion animals are always provided the best by their owners, whose demand for high-quality pet food is very high [
2]. One key aspect that pet food owners and manufacturers consider is the protein source and content [
3]. The content of protein in pet food is high, which in canine food ranges from 17.3% to 36.6% [
4]. To satisfy the standard of pet food as human food, some companies and some owners choose not to use or feed products containing by-product meal [
5]. Some pet foods do not contain plant protein, and some pet owners have expressed concern that gluten in grains may be a source of allergies in dogs [
6]. So, it is a big cost of animal protein originally supplied to humans in pet food. Attention to the environment, animal welfare, and climate change are encouraging institutions and individuals to seek alternatives to conventional animal proteins.
Along with the 32,000-year history of the parallel evolution between dogs and humans and adapting to agricultural-based living conditions, dogs evolved from carnivores to omnivores due to large changes in their food source [
7]. Both animal ingredients and plant ingredients containing large amounts of protein and starch can be digested and absorbed by dogs [
8,
9]. To save animal protein-sourced protein and protect the environment, it would be a better way to use more plant-sourced protein and recycled protein to replace part of animal-sourced protein in pet food.
Traditional animal protein, such as chicken, fish, and meat by-products, and plant proteins, such as soybean meal and corn gluten meal, have been the major protein ingredients used in the formulation of commercial pet foods [
10,
11,
12]. In recent years, more novel sustainable ingredients with high protein, such as insect meal and single-cell protein, have also entered the pet food market [
11,
13]. Yeast extract is the water-soluble extract produced from yeast waste streams, such as
Saccharomyces cerevisiae, and separated from inner yeast cells. It could be a functional source of nutrients, since yeast extract is rich in proteins, amino acids, nucleotides, sugars, and a variety of trace elements [
14]. Lin et al. showed that yeast products may be beneficial to adult dogs by positively altering the gut microbiota, enhancing immune capacity, and reducing inflammation [
15]. Because of insect proteins’ low land use, lower greenhouse gas emissions, and low water pollution, they may contribute to sustainable food production as an alternative source of animal protein. Mealworm meal is a kind of by-product after yellow mealworms (
Tenebrio molitor) larvae defatted, which has a high quality and quantity of protein and amino acids [
16,
17]. Pet food with insect-based ingredients was poorly accepted for human consumption to feed their pets. Insect-based pet foods proved to be attractive for purchase only when consumers were well informed about the product’s properties in terms of sustainability and healthiness for their pets [
18]. So, it is necessary to analyze different protein ingredients in pet food comprehensively and systematically to help pet food companies and pet owners know the utilization of these environment-friendly protein sources.
Knowledge of the energy values and digestibility of ingredients is important to correctly balance pet food [
19,
20]. Current research on the effective energy value of pet food typically recommends the modified Atwater equation or predictive equations based primarily on fixed energy values and digestibility coefficients associated with the chemical composition of diets to estimate the metabolizable energy (ME) content of pet foods [
20,
21,
22,
23]. However, the equations do not apply to all ingredients and may overestimate the food energy of animals [
24,
25]. The effective energy value of pet food is based on the ME energy system. In opposition to ME, net energy (NE) is a more precise evaluation of the true energy value of the feed [
26], because it takes the heat increment (HI) from the digestion and metabolism of feeds into account [
27,
28]. We can find the rule of energy metabolism of different protein ingredients by total heat production (THP) and HI, which would be useful for losing weight in pets and patients during nutrition recovery.
The most accurate method to assess the effective energy value of feed is to evaluate the animal’s real digestive and metabolic conditions in vivo. The difference method is more suitable for determining the nutrient digestibility and the effective energy value of single ingredients in vivo [
29,
30,
31,
32]. A 30% substitution ratio in the difference method in vivo has been shown to effectively assess the energy content of poultry by-product meal for beagles in our previous study [
20]. Traditional protein sources in pet foods include poultry by-product meal, fish meal (FM), meat and bone meal (MBM), corn gluten meal (CGM), and soybean meal (SBM) [
33,
34,
35]. Recently, high-quality sustainable protein resources such as insect meals and yeast products have also been used in pet foods [
36,
37,
38].
This study aimed to determine the effective energy values of FM, MBM, CGM, SBM, mealworm meal (MM), and yeast extract (YE) by using the difference method, measure the nitrogen metabolism and heat production, and assess the feces score for beagles. Through stepwise regression analysis of the measured energy value and chemical composition of ingredients fed to beagles, we also derived predictive equations for the effective energy value of protein ingredients.
2. Materials and Methods
2.1. Diets and Ingredients
The basic diet (BD) was formulated to satisfy the nutritional needs of adult canines [
21].
Table 1 shows the seven test diets, including the BD and diets involving 30% substitution of the BD, each replaced with fish meal (FM), meat and bone meal (MBM), corn gluten meal (CGM), soybean meal (SM), mealworm meal (MM), and yeast extract (YE), respectively. All diets were mixed uniformly in powder form.
2.2. Animals, Housing, and Experimental Design
The experiment took place at the companion animal testing center of the Special Animal and Plant Sciences of the Chinese Academy of Agriculture Sciences (Changchun, China). Beagles were kept in indoor enclosures covering floor space, adhering to prescribed light cycles, temperatures, and sanitation practices following the Animal Welfare Act guidelines. Before the experiment, all dogs had undergone deworming and vaccination, and no medications were administered throughout the study [
20].
Throughout the experiment, except for the fasting period, dogs were provided with two meals of equal size at 09:00 and 14:00, with unrestricted access to fresh water. And the daily food intake of each beagle was recorded. All diets were provided as a mixture blended with water; the ratio between powder and water was 1:2.
The average weight of the six healthy adult female beagles was 15.07 ± 2.15 kg, and their body condition score (BCS) ranged from 4.5/9 to 5.5/9 [
39]. The six dogs were each fed one of the seven diet treatments, according to a 7 × 6 incomplete Latin square design.
The beagles were individually housed in respiration chambers with a volume of 0.42 m
3 [
40]. Indirect calorimetry was performed as described by Zhang et al. [
20] and conducted for seven periods. Each experiment period lasted for 10 days, including a 3-day adaptation period followed by a 7-day testing period (including a 4-day feeding period and a 3-day fasting period). Between each experiment period, have a 7-day washing period fed on BD. The beagles were weighed at the start of the feeding period and at the start and end of the fasting period. At 09:00 a.m. on d 0 of each experiment period, beagles were transferred to the chamber to adapt. Each dog was changed into a living chamber in each experiment period in proper order. Throughout the 7-day testing period, O
2 consumption and CO
2 production volumes were measured continuously for 4 consecutive days to assess total heat production (THP) and 3 consecutive days to assess fasting heat production (FHP), employing the Brouwer equation [
41].
2.3. Fecal Score
During the feeding period, the fecal samples of each dog were scored every day. Fecal score was used using the following 5-point system: 1 = very hard, dry pellets. 2 = hard, formed, remains firm and soft; 3 = soft, formed, retains shape; 4 = unformed stool, pasty and slushy; and 5 = watery diarrhea [
42,
43]. The ideal fecal score was 2 to 3, indicating well-formed stools that were convenient to collect without being excessively dry [
44].
2.4. Sample Collection
During the feeding period, total feces from each dog were weighed and collected once daily for 4 days. All fecal samples were stored at −20 °C. At the end of each experiment period, feces samples from each dog during each feeding period were uniformly mixed and dried at 65 °C, then smashed and sifted with a 1 mm screen before chemical analysis.
Urine was collected daily during the 7-day testing period and then mixed with 10 mL of 10% sulfuric acid and measured for volume. Urine samples were separately mixed at the end of the feeding and fasting periods for each dog with each testing period, and then stored at −20 °C until analysis.
2.5. Chemical Analyses
Diets, ingredients, and feces were analyzed for dry matter (DM) (AOAC method 934.01 [
45]). Nitrogen in all the diets, ingredients, feces, and urine samples was determined using the standard procedure (AOAC method 984.13 [
45]), and crude protein (CP) was calculated as nitrogen × 6.25. The ether extract (EE), ash, calcium (Ca), phosphorus (P), and amino acids (AAs) contents in the diets, ingredients, and fecal samples were analyzed with AOAC 920.39, 967.05, 968.08, 985.01, and 994.12 [
45]. The gross energy (GE) in the diets, ingredients, feces, and urine samples was determined using an adiabatic bomb calorimeter (IKA C2000, Staufer, Germany), with benzoic acid employed as the standard. The aflatoxin B1 and vomitoxin contents of the ingredients were determined by the use of ELISA kits (Sinobestbio Co., Ltd., Shanghai, China).
2.6. Calculations
The apparent total tract digestibility (ATTD) of energy and nutrients of test diets was calculated using the following equation:
The ATTD and effective energy value of test ingredients were calculated as previously described by Adeola [
46]:
where ID, TDD, and BDD were the apparent digestibility of the ingredients, test diets, and BD, respectively (%); IE, TDE, and BDE were the energy value of the ingredients, test diets and BD, respectively, (MJ/kg DM); and X was the substitution ratio of the ingredients.
The values of DE, ME, and NE in the diet were calculated as follows [
26]:
The THP and HI of beagles were determined using the following equations [
41]:
where VO
2 was O
2 consumption, and VCO
2 was CO
2 production. To account for the effect of body weight on energy metabolism and respiration between animals, the data were converted to metabolic weight [
20].
2.7. Statistical Analyses
The data were presented in the format of the mean ± SEM and analyzed by using one-way ANOVA for energy value, nitrogen balance, O2 consumption, and CO2 production. Distinctions among diets or ingredients were assessed through Duncan’s multiple range test, with a significance level set at p < 0.05. Pearson’s correlation analysis was conducted to explore associations among various nutrients, energy values of ingredients, and equations. The estimation of equations was conducted using multiple linear regression through the stepwise method in SPSS 25.0 (SPSS Inc., Chicago, IL, USA). A graphical representation of correlation coefficients was generated using GraphPad Prism 9.0 software.
3. Results
3.1. Nutrient Composition of Test Ingredients
The analyzed chemical composition of ingredients (DM basis) is shown in
Table 2. The analyzed content of CP in the six test ingredients is listed in decreasing order as MM, FM, CGM, MBM, SBM, and YE, and all the test ingredients had a protein level greater than 40%. The concentrations of ash 34.83%, Ca 12.77%, and P 5.43% were found to be greater in MBM than in the other ingredients. Compared with CGM, SBM, and YE, MM, MBM, and FM had a greater EE content.
CGM had the highest gross energy content of 22.76 MJ/kg, while MBM had the lowest at 16.13 MJ/kg. Among the test ingredients, MM contained the highest levels of cysteine, threonine, arginine, valine, and leucine; FM was higher in lysine, histidine, and isoleucine; and CGM had the highest methionine, tyrosine, and phenylalanine content, which matched the higher CP content of the ingredients.
3.2. The Energy Values and the ATTD of GE and Nutrients of Diets
The ATTD of CP in the MBM diet was significantly lower than that of the BD, FM, and CGM diets (
p < 0.05) (
Table 3). The ATTD of CF in beagles fed FM and SBM diets was lower when compared with other diets (
p < 0.05). The ATTD of DM in the MBM diet was significantly lower than in other diets (
p < 0.05). Moreover, the ATTD of organic matter (OM) and GE in MM was the lowest among the diets.
In terms of the energy value content of the test dietary diets, the gross energy of the MM diet was higher than that of other diets (p < 0.05). The FE values of the BD, FM, and CGM diets were significantly lower than those of other diets (p < 0.05). The UE of the CGM diet was the highest at 0.95 MJ/kg, significantly higher than the BD, FM, and MBM diets (p < 0.05).
The MM diet had the highest levels of DE and ME at 18.46 MJ/kg and 17.80 MJ/kg, while the MBM diet had the lowest at 13.31 MJ/kg and 12.61 MJ/kg. No significant variations were observed in NE between the BD, FM, and MM diets (p > 0.05).
The energy conversion efficiency of the ME:GE ratio of the FM diet was significantly greater than the MBM, SBM, and YE diets (p < 0.05). There were no significant differences seen for the ME:DE and NE:ME ratios (p > 0.05). The ME:DE ratio ranged from 94.63 to 97.48% among the seven diets, while the range of NE:ME is 75.47% to 86.07%. The ratios of NE:ME of the BD, FM, and MM diets were all above 80%.
3.3. Nitrogen Balance and Heat Production for Different Diets in Beagles
The data on the effects of test diets on nitrogen balance and heat production in beagles are presented in
Table 4. No significant effect was observed for ME intake among the diets (
p > 0.05). THP and HI were unaffected by the diets (
p > 0.05). The HI of the diets listed in descending order as the YE, CGM, MBM, SBM, FM, and MM diets, and BD as the lowest one. There were no effects of NI, UN, RN, NPU, or PBV among the diets (
p > 0.05). The FN of BD was significantly lower than the FM, MBM, CGM, MM, and YE diets (
p < 0.05).
3.4. The Energy Values and the ATTD of Nutrients of the Test Ingredients
The ATTD of nutrients, as well as the DE, ME, and NE content of test ingredients, are shown in
Table 5. Beagles fed FM, CGM, SBM, and YE had greater ATTD of DM and OM compared to those fed MM (
p < 0.05). No distinctions were observed in the ATTD of CP and GE between the MBM and MM (
p > 0.05), but they were lower compared to the other four ingredients (
p < 0.05). The ATTD of CF in SBM was significantly lower than MBM, CGM, MM, and YE (
p < 0.05). Overall, the ATTD of nutrients among the six ingredients was the lowest for MBM and MM and the highest for FM and CGM.
The energy value of the six ingredients was significantly different (p < 0.05). The DE values (MJ/kg DM) in descending order were MM at 22.95, CGM at 17.46, FM at 16.48, SBM at 15.36, YE at 15.11, and MBM at 6.73, and MM was significantly higher in comparison to the remaining five ingredients (p < 0.05). The ME content of MBM was significantly lower in comparison to the other five ingredients (p < 0.05). The NE of the FM, CGM, and MM were higher than that of the MBM (p < 0.05). In terms of energy utilization efficiency for the test ingredients, the ME:DE ratio ranged from 68.85% to 97.25%, with MM being significantly lower compared to the other ingredients (p < 0.05), and the NE:ME ratio ranged from 60.86% to 94.42%.
3.5. The Prediction Equations of the Energy Values of Poultry By-Product Meal, Fish Meal, Meat and Bone Meal, Corn Gluten Meal, Soybean Meal, Mealworm Meal, and Yeast Extract for Beagles
Combined with the findings of earlier research conducted by our team [
20], the correlation between nutrient composition and DE, ME, and NE content of PBM and the six protein ingredients tested is presented in
Figure 1. The ash content exhibited a negative correlation with OM (
p < 0.01). A negative correlation (
p < 0.01) was observed between the content of carbohydrate and CP (
p < 0.01) and EE content. The content of CF exhibited a negative correlation with EE (
p < 0.01) and a positive correlation with carbohydrates (
p < 0.01).
The DE value demonstrated a positive correlation with ME, NE, and CP (p < 0.01) and a negative correlation with ash content (p < 0.01). The ME value showed a positive correlation with NE and OM (p < 0.01) and a negative correlation with ash (p < 0.01). The NE content was positively correlated with the CP content (p < 0.01).
Based on energy values and chemical composition, stepwise regression analysis was conducted to establish predictive equations for the effective energy, such as DE, ME, and NE (MJ/kg DM), of the seven ingredients, as shown in
Table 6. The GE was the first predictor of DE content with R
2 = 0.889 and RSD = 1.487 (
p < 0.001), however, the precision of the equation was enhanced when CF was involved in the predictive equation with R
2 = 0.964 and RSD = 0.845 (
p < 0.001). The DE content had a strong correlation with the ME content, so it could be used as the only predictor in the ME prediction equations, where R
2 = 0.799 and RSD = 0.117 (
p < 0.001).
Protein and fiber content can serve as predictors of the effective energy value content of the ingredients. The prediction equations for DE, ME, and NE of the seven diets were: DE = 26.991 − 0.521ash − 0.143CHO − 0.446CF + 0.266EE where R2 = 0.964 and RSD = 0.845; ME = 16.521 − 0.267ash − 0.319GE − 0.287CF + 0.16CP where R2 = 0.919 and RSD = 0.899; and NE = 0.303 + 0.212CP − 0.154EE − 0.146CF where R2 = 0.930 and RSD = 0.582.
The deviations between the calculated values using the prediction equation in this study and the measured values of ME for MBM, MM, YE, CGM, FM, and SBM were all less than 10%, as shown in
Table 7. When calculated using the NRC recommended equation, the deviations in ME for FM, CGM, SBM, MM, and YE compared to the measured were below 10%, with MBM reaching 67.13%. The differences between the calculated values using the prediction equation and the measured values of NE for six ingredients were all less than 10%.
3.6. Fecal Characteristics
Fecal characteristics, including fecal moisture content and fecal score, are shown in
Table 8. Dogs fed the SBM diet had a higher fecal moisture content of 73.49% contrasted with the other diets (
p < 0.05), and the MBM diet had the lowest fecal moisture content of 56.59% among the diets. The feces of the MBM diet contained more than 40% ash, and the other diets contained less than 20%. All fecal scores were within an acceptable range using the 5-point scale referenced previously. The YE diet had the highest fecal score of 3.19 compared with BD at 2.47 (
p < 0.05) and did not differ among the FM, MBM, CGM, SBM, and MM diets (
p > 0.05).