Obese dogs exhibit different fecal microbiome and specic microbial networks compared with normal weight dogs

Page created by Willard Ramsey
 
CONTINUE READING
Obese dogs exhibit different fecal microbiome and specic microbial networks compared with normal weight dogs
Obese dogs exhibit different fecal microbiome and
specific microbial networks compared with normal
weight dogs
Hanbeen Kim
 Life and Industry Convergence Research Institute, Pusan National University
Jakyeom Seo
 Life and Industry Convergence Research Institute, Pusan National University
Tansol Park
 Chung-Ang University
Kangmin Seo
 National Institute of Animal Science
Hyunwoo Cho
 National Institute of Animal Science
Julan Chun
 National Institute of Animal Science
Kihyun Kim (  kihyun@korea.kr )
 National Institute of Animal Science

Article

Keywords: Body condition score, Obesity, Fecal microbiome, Specific network, Dog

Posted Date: June 6th, 2022

DOI: https://doi.org/10.21203/rs.3.rs-1616770/v2

License:   This work is licensed under a Creative Commons Attribution 4.0 International License.
Read Full License

                                               Page 1/25
Obese dogs exhibit different fecal microbiome and specic microbial networks compared with normal weight dogs
Abstract
Background
Canine obesity is a major health concern that predisposes dogs to various disorders. The fecal
microbiome has been attracting attention because of their impact on energy efficiency and metabolic
disorders of host. However, little is known about specific microbial interactions, and how these may be
affected by obesity in dogs. The objective of this study was to investigate the differences in fecal
microbiome and specific microbial networks between obese and normal dogs. A total of 20 beagle dogs
(males = 12, body weight [BW]: 10.5 ± 1.08 kg; females = 8, BW: 11.3 ± 1.71 kg; all 2-year-old) were fed to
meet the maintenance energy requirements for 18 weeks. Then, 12 beagle dogs were selected based on
body condition score (BCS) and divided into two groups: high BCS group (HBCS; BCS range: 7–9, males
= 4, females = 2) and normal BCS group (NBCS; BCS range: 4–6, males = 4, females = 2). In the final week
of the experiment, fecal samples were collected directly from the rectum, before breakfast, for analyzing
the fecal microbiome using 16S rRNA gene amplicon sequencing.

Results
The HBCS group had a significantly higher final BW than the NBCS group (P < 0.01). The relative
abundances of Faecalibacterium, Phascolarctobacterium, Megamonas, Bacteroides, Mucispirillum, and
an unclassified genus within Ruminococcaceae were significantly higher in the HBCS group than those in
the NBCS group (P < 0.05). Furthermore, some Kyoto Encyclopedia of Genes and Genomes (KEGG)
modules related to amino acid biosynthesis and B vitamins biosynthesis were enriched in the HBCS
group (P < 0.10), whereas those related to carbohydrate metabolism were enriched in the NBCS group (P <
0.10). Microbial network analysis revealed distinct co-occurrence and mutually exclusive interactions
between the HBCS and NBCS groups.

Conclusions
In conclusion, several genera related to short-chain fatty acid production were enriched in the HBCS
group. The enriched KEGG modules in the HBCS group enhanced energy efficiency through cross-feeding
between auxotrophs and prototrophs. However, further studies are needed to investigate how specific
networks can be interpreted in the context of fermentation characteristics in the lower gut and obesity in
dogs.

Background
Owing to the growing interest in dog welfare, more people have started to take an interest in the health
problems of dogs [1]. Among them, canine obesity is considered a major metabolic disease because it
can predispose dogs to a variety of disorders, such as cardiovascular disease [2], expiratory airway

                                                  Page 2/25
Obese dogs exhibit different fecal microbiome and specic microbial networks compared with normal weight dogs
dysfunction [3], and diabetes mellitus [4], and has detrimental effects on longevity [5]. The main cause of
obesity in dogs is an imbalance between energy intake and expenditure [6]. Furthermore, some studies
suggest that breed, feed quality, and the owner’s interest in maintaining pet’s health are also possible
causes [7–9]. Therefore, the veterinary community is paying increasing attention to the prevention of
obesity in dogs.

Previously, companion dogs have been evaluated for being overweight and obese based on their body
weight (BW), body condition score (BCS), and body composition (fat mass and lean body mass) [10–12].
These methods have been recognized to be efficient in terms of cost and convenience. However, given the
limited data on the ideal weight for each breed and the subjective opinions of various veterinarians, these
methods need to be improved for better reliability [6]. Accordingly, recent advances in omics technology
are expected to replace traditional methods and provide new methods for investigating obesity indicators.

The gastrointestinal tract of dogs, especially the colon, is habituated by a complex microbiota that
contributes to energy source absorption, metabolic functions, and immunogenic responses [13]. Several
studies have investigated changes in the fecal microbiota in response to the modulation of dietary
components, including carbohydrates [14], proteins [15], and fats [16]. In addition, some studies have
reported differences in the fecal microbiota of obese dogs after a weight-loss program [17–20]. Taken
together, these studies highlight the importance of the fecal microbiota in sustaining canine health and
controlling obesity.

Previous studies have also reported differences in the gut microbiota among obese and lean individuals
in humans and mice [21–25]. These changes are characterized by a reduced diversity of gut microbiota
[21, 22], increased Firmicutes/Bacteroidota ratio [23, 24], and better energy harvesting [25] in the obese
group.

However, only a few studies have investigated the differences in the fecal microbiota between obese and
lean dogs [17, 20, 26]. The composition of the gut microbiota in dogs can differ based on several factors,
including diet, age, breed, and disorders [27–29]. Previous studies had several limitations in interpreting
the relationship between obesity and gut microbiota, including the lack of unification in terms of breeds
or ages [17, 20, 26]. Therefore, the present study investigated the effects of obesity on the gut microbiota
and microbial functional features of dogs of a similar age and breed. We hypothesized that fecal
microbiota, microbial functional profiles, and networks among microbiota would differ between obese
and normal dogs.

Therefore, this study aimed to investigate differences in fecal microbiota composition, microbial
functional profiles, and specific networks between obese and normal BCS groups after maintenance
energy requirements were met.

Methods
Animals and experimental design
                                                  Page 3/25
This study was conducted using neutered beagles owned by the National Institute of Animal Science in
Korea. This experiment was conducted in accordance with the methods approved by the Animal Care and
Use Committee National Institute of Animal (NIAS-2019-370). All methods for animal care are performed
in accordance with the ARRIVE guidelines. Twenty beagle dogs (males = 12, initial body weight [BW]: 10.5
± 1.08 kg; females = 8, initial BW: 11.3 ± 1.71 kg; all 2-year-old) were enrolled in this study. The beagles
lived individually in an indoor space (1.8 m × 2.6 m) maintained at a temperature of 22–24 °C and
relative humidity from 60% to 80% using automatic temperature control and continuous ventilation. All
dogs were fed twice a day (at 10:00 and 16:00), had access to drinking water ad libitum, and engaged in
approximately 3 h of outdoor activity every day. The dogs were fed an extruded pellet-type commercial
dry food (Iskhan All-life33®; Wooriwa Ltd., Korea) that met the maintenance energy requirements
estimated based on the formula suggested by the Association of American Feed Control Officials
(AAFCO; metabolizable energy requirement = 132 kcal × kg BW0.75/day) [30] (Table 1). The experiment
was conducted for 18 weeks. In the first and final day of the experiment, the BCS was measured on a
nine-point scale according to the criteria of Laflamme D [12], and 12 neutered beagles (eight males and
four females, all 2-years-old) were selected based on the final BCS 24 h before collecting feces and
divided into two groups, each consisting of four males and two females. The normal BCS group (NBCS)
was composed of six healthy dogs with a normal BCS (4–6), and the high BCS group (HBCS) was
composed of six dogs with a high BCS (7–9) (Table 2).

Fecal sample collection and DNA extraction

Twenty-four hours before collecting feces, feeding was restricted. Fecal samples were collected directly
from the rectum before breakfast (9:00 to 10:00), immersed in liquid nitrogen, and stored in a deep freezer
at -80 °C until analysis. DNA was extracted from the fecal samples using the NucleoSpin DNA Stool kit
(Macherey-Nagel, Germany) according to the manufacturer’s instructions. Briefly, fecal samples (180–220
mg) were lysed using the MN Bead Tube Holder, and the genomic DNA was harvested through
NucleoSpinDNA stool column centrifugation at 13,000 × g for 1 min. In the final step, genomic DNA was
eluted into the elution buffer after centrifugation at 13,000 × g for 1 min and stored at -20 ℃.

16S rRNA gene amplicon sequencing and data processing

For PCR amplification targeting the V3 and V4 variable regions of 16S rRNA, the following primer
sets were used [31]: forward primer 5′-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG CCT ACG
GGN GGC WGC AG-3′ and reverse primer 5′-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GGA
CTA CHV GGG TAT CTA ATC C-3′. Illumina adapter overhang nucleotide sequences were added to the
gene-specific sequences. The locus-specific sequences were as follows: forward overhang 5′-TCG TCG
GCA GCG TCA GAT GTG TAT AAG AGA CAG-3′ and reverse overhang 5′-GTC TCG TGG GCT CGG AGA TGT
GTA TAA GAG ACA G-3′. To purify the PCR products, KAPA HiFi HotStart ReadyMix (KAPA Biosystems,
USA) and Agencourt AMPure XP system (Beckman Coulter Genomics, USA) were used. The libraries were
sequenced using the Illumina MiSeq instrument (2 × 300 paired-end sequencing; Illumina, CA, USA).

                                                  Page 4/25
The barcode sequences were trimmed using Cutadapt (version 1.11) [32]. Then, the raw amplicon
sequences were analyzed using built-in plugins within Quantitative Insights into Microbial Ecology 2
(QIIME2) (version 2019.04) [33]. First, the DADA2 plugin was used to remove primer sequences, denoise
low-quality reads (Q < 25), merge forward and reverse reads, and remove chimeric sequences. The raw
16S rRNA gene sequences are available in the National Center for Biotechnology Information sequence
read archive (BioProject accession: PRJNA818097). The amplicon sequence variants (ASVs) were
taxonomically classified using the Silva 16S rRNA gene database (version SSU138) [34]. Taxonomically
unassigned ASVs were excluded. The alpha diversity of each sample was examined with respect to ACE
(abundance-based coverage estimator) and Chao1 estimates, observed ASVs, evenness, Shannon
diversity index, and Simpson diversity index based on rarefied ASV tables using randomly selected 67,151
ASVs per sample. Principal coordinate analysis (PCoA) was conducted based on unweighted and
weighted UniFrac distance matrices to compare the overall dissimilarity of microbiota between the HBCS
and NBCS groups and visualized using Emperor implemented in QIIME2. The number of common and
exclusively identified bacterial taxa between the HBCS and NBCS groups at the phylum, family, and genus
levels were visualized using Venn diagrams. Functional metagenomic features from the 16S rRNA gene
sequences were predicted based on ASVs and the corresponding biological observation matrix (BIOM)
table using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved State
(PICRUSt2) (version 2.3.0-b) with default options [35]. Principal component analysis (PCA) was
performed to compare the overall differences in predicted functional profiles based on the Kyoto
Encyclopedia of Genes and Genomes (KEGG) orthologs between the HBCS and NBCS groups. The PCA
plot was visualized using the ggfortify R package [36]. FastSpar was used to examine the correlations
among the major bacterial genera in the HBCS and NBCS groups [37], and the correlations were
computed based on the sparse correlations for compositional data [38]. Only significant microbial
interactions (P < 0.05) were further analyzed using co-expression differential network analysis
(CoDiNA) [39] to identify specific microbial networks in the HBCS and NBCS groups.

Statistical analysis

Before statistical analysis, all data were tested for normal distribution using the Kolmogorov–Smirnov
test in SAS 9.4 (SAS Institute Inc., NC, USA). Then, normally distributed data, including those on animal
performance (BW, BCS, and BW gain) and alpha diversity indices (ACE, Chao1, observed ASVs, evenness,
Shannon index, and Simpson diversity index) were analyzed using a t-test in SAS 9.4. Data that were not
normally distributed, including major microbial taxa and predicted microbial functional features from
KEGG categories and modules, were analyzed using a non-parametric Wilcoxon rank-sum test.
Permutational multivariate analysis of variance (PERMANOVA) was used to further analyze the PCoA and
PCA results to assess whether the overall microbial community and predicted functional profiles (KEGG
orthologs) differed significantly between the HBCS and NBCS groups. Statistical significance was
declared at P < 0.05, and a statistical trend was speculated at 0.05 ≤ P ≤ 0.10.

Results
                                                 Page 5/25
Relation between BCS and animal performance

Both the initial BW and BCS were higher in the HBCS group than those in the NBCS group (Table 2). After
feeding to meet maintenance energy requirements for individual dogs, the HBCS group had a significantly
higher final BW and BCS than the NBCS group (P < 0.05), resulting in a higher BW gain in the HBCS group
(Table 2).

Quality assessment and sample statistics

A total of 1,826,502 sequences were obtained from 16S rRNA amplicon sequencing analysis, with an
average of 152,209 ± 30,573 sequences per fecal sample (Table S1). After quality filtering using QIIME 2
(Q score > 25), 1,193,079 sequences (65% of the raw reads) were generated, with an average of 99,423 ±
19,814 sequences per sample. Good’s coverage was greater than 99% for all fecal samples.

Relation between BCS and the alpha and beta diversity of fecal microbiota

In fecal microbiota, the difference in BCS did not result in significant changes in any of the alpha diversity
measurements (Table 3). The PCoA results based on both unweighted and weighted UniFrac
distances also did not present any significant difference in the fecal microbiota of the HBCS and NBCS
groups (Figure 1).

Relation between BCS and the fecal microbiota

In this study, only bacterial taxa have been presented because archaeal reads could not be detected. The
bacterial taxa that were detected in over 50% of the samples and had a relative abundance of over 0.05%
in at least one of the BCS groups were defined as major microbial taxa.

The Venn diagram shows common and exclusively identified prokaryotic taxa in the HBCS and NBCS
groups (Figure 2). Eight of the nine phyla were shared between the HBCS and NBCS groups, and only
Desulfobacterota was exclusively detected in the NBCS group (Figure 2A); moreover, Desulfobacterota
was detected in only one of the six samples. At the family level, 50 families (approximately 75.8%) were
shared between the HBCS and NBCS groups and no major families were found exclusively in the BCS
group (Figure 2B). At the genus level, more than 100 genera (approximately 72.7%) were shared between
the HBCS and NBCS groups, and only 14 and 24 genera were exclusively detected in the HBCS and NBCS
groups, respectively (Figure 2C). Succinivibrio and an unclassified genus within Atopobiaceae were the
only major genera (present in at least 50% of the samples in each group) found exclusively in the HBCS
and NBCS groups, respectively; however, they only accounted for a minor proportion (relative abundance
< 0.01%).

The major phyla, families, and genera with relative abundances of over 0.5% in at least one group are
shown in Figure S1. Among them, the major taxa that showed statistical tendency (P < 0.10) are shown in
Figure 3 (relative abundance ≥ 0.05% in at least one group). Firmicutes was the most predominant
phylum (accounting for at least 77% of the total microbiota), followed by Fusobacteriota (Figure S1A). Of
                                                  Page 6/25
the five predominant phyla, the relative abundance of Bacteroidota tended to be higher in the HBCS group
than that in the NBCS group (HBCS, 7.61% vs. NBCS, 2.32%; P = 0.0656), whereas Actinobacteria showed
an opposite trend (HBCS, 2.62% vs. NBCS, 4.32%; P = 0.0927) (Figure 3A). The relative abundance of
Deferribacterota, accounting for a minor proportion, was significantly higher in the HBCS group than that
in the NBCS group (HBCS, 0.051% vs. NBCS, 0.008%; P = 0.0278) (Figure 3A). At the family level,
Lachnospiraceae, Peptostreptococcaceae, and Lactobacillaceae were the three most abundant families
in the HBCS group, whereas Peptostreptococcaceae, Erysipelotrichaceae, and Lactobacillaceae were the
most abundant in the NBCS group (Figure S1B). Among the major classified families, the relative
abundances of Ruminococcaceae, Bacteroidaceae, Acidaminococcaceae, Selenomonadaceae, and
Deferribacteraceae were higher in the HBCS group than those in the NBCS group (Figure 3B). At the genus
level, Peptoclostridium was the most predominant genus in both BCS groups (HBCS, 14.90% vs. NBCS,
23.79%), and more than 84% of the entire fecal microbiota was assigned to 54 major classified genera
(representing > 0.05% of relative abundance in at least one BCS group). At the genus level, the relative
abundances of Bacteroides, Mucispirillum, and the four genera within Firmicutes (Faecalibacterium,
Phascolarctobacterium, Megamonas, and UG_Ruminococcaceae) were significantly higher in the HBCS
group than those in the NBCS group (Figure 3C, P < 0.05). Meanwhile, the relative abundances of
Sutterella, Prevotellaceae Ga6A1 group, and four genera within Firmicutes (UG_Erysipelotrichaceae,
Lachnospiraceae NK4A136 group, Catenibacterium, and Fournierella) tended to be higher in the HBCS
group than those in the NBCS group (Figure 3C, P < 0.10).

Relation between BCS and the predicted functional profiles from fecal microbiota

The overall distribution of the microbial functional profiles predicted based on KEGG orthologs was
visualized using PCA plots (Figure 4A). There was no significant difference in the overall distribution of
KEGG orthologs between the HBCS and NBCS groups based on the PERMANOVA analysis. Among the
top 10 KEGG categories (Level 2), only carbohydrate metabolism was significantly higher in the NBCS
group than that in the HBCS group (Figure 4B, P < 0.05). Among the major KEGG modules, three modules
related to amino acid metabolism (lysine biosynthesis [M00030], serine biosynthesis [M00020], and
tryptophan biosynthesis [M00023]) and five modules related to the metabolism of cofactors and vitamins
(cobalamin biosynthesis [M00122], pyridoxal-P biosynthesis [M00124], phylloquinone biosynthesis
[M00932], and lipoic acid biosynthesis [M00881, M00884]) were markedly enriched in the HBCS group
(Figure 4B, P < 0.10). In contrast, three modules related to carbohydrate metabolism (galactose
degradation [M00632], nucleotide sugar biosynthesis [M00549], and D-galactonate degradation
[M00552]), two modules related to energy metabolism (reductive pentose phosphate cycle [M00165] and
methanogenesis [M00357]), and one related to cysteine biosynthesis [M00609] tended to be higher in the
NBCS group (Figure 4C, P < 0.10).

Differential network analysis based on BCS

The statistics for the differential networks of the major genera in the HBCS and NBCS groups are listed in
Table 4. The total number of nodes for the HBCS and NBCS groups was 31 and 39, respectively, and the

                                                 Page 7/25
total number of edges was 74 and 71, respectively. In the differential network of the HBCS group, positive
edges were detected more frequently than negative edges (Table 4; 51 positive edges vs. 23 negative
edges), whereas the number of positive edges was similar to that of negative edges in the differential
network of the NBCS group (Table 4; 36 positive edges vs. 35 negative edges). In the differential networks
of each treatment, the γ link, which was specific to one of the networks, was assigned to 29% of the HBCS
group and 23.1% of the NBCS group (Table 4). The specific genera (γ link) in the differential network of
the HBCS group were Allobaculum, Lachnospiraceae NK4A136 group, UG_Ruminococcaceae, Collinsella,
UG_Erysipelotrichaceae, Phascolarctobacterium, Muribaculaceae, Megamonas, and Fusobacterium
(Figure 5A), and those of NBCS group were Veillonella, Prevotella, Negativibacillus, Bifidobacterium,
Bacillus, Oribacterium, Turicibacter, Faecalibacterium, and Helicobacter (Figure 5B). Based on the two
centrality measurements (degree and eigenvector centrality) and hub score within each differential
network, UG_Erysipelotrichaceae and Phascolarctobacterium were the keystone genera in HBCS, while
Negativibacillus and Bifidobacterium were those in the NBCS group (Figure 5 and Figure S2).

Discussion
Several human obesity studies have reported that obesity is generally linked to altered gut microbiota,
presenting lower diversity in obese humans than that in lean humans [21, 22]. In addition, a recent
metagenomic study revealed that the dog microbiota was similar to human microbiota in comparison
with pig or mouse microbiota, and the alteration of microbiota in response to dietary changes was also in
agreement with that observed in human studies [40]. Therefore, we expected the HBCS group to have a
lower microbiota diversity than the NBCS group. However, the present study revealed no significant
differences in the alpha and beta diversity measurements between the HBCS and NBCS groups. Similar
to the results of microbial diversity in this study, Handl S, German AJ, Holden SL, Dowd SE, Steiner JM,
Heilmann RM, Grant RW, Swanson KS and Suchodolski JS [17] reported no significant differences in
alpha diversity (Shannon’s index, evenness, and the number of operational taxonomic units, i.e. OTUs)
and beta diversity between lean and obese companion dogs. In addition, a recent study that evaluated
differences in the fecal microbiota among obese, lean, and weight loss dogs showed that there were no
significant differences in alpha diversity (Faith’s phylogenetic diversity, Shannon’s index, and the number
of OTUs) and beta diversity, whereas the evenness index was significantly higher in the obese group than
that in the lean and weight loss groups [20].

Previous studies have reported that the five predominant phyla in the fecal microbiota of dogs are
Firmicutes, Bacteroidota, Proteobacteria, Actinobacteriota, and Fusobacterium [27, 40] as shown in the
present study. The Firmicutes/Bacteroidota ratio is generally considered an indicator of obesity because
several studies in humans and mice have shown that the obese group has a higher
Firmicutes/Bacteroidota ratio than the lean group [23, 24]. In addition, a previous study proposed that
Firmicutes was more effective at harvesting energy from a diet than Bacteroidetes; therefore, a higher
Firmicutes/Bacteroidota ratio may promote weight gain [41]. However, in the present study, the relative
abundance of Bacteroidota tended to be higher in the HBCS group, whereas there was no significant
difference in the relative abundance of Firmicutes, thereby reducing the Firmicutes/Bacteroidota ratio in
                                                 Page 8/25
the HBCS group (data not shown). Meanwhile, some studies on gut microbiota composition in obese
humans or mice have also reported conflicting results, such as a lower Firmicutes/Bacteroidota ratio [42]
or no difference in the Firmicutes/Bacteroidota ratio [43]. In addition, several studies have observed no
significant differences in both Firmicutes and Bacteroidota between obese and lean dogs [15, 17, 44].
Park HJ, Lee SE, Kim HB, Isaacson R, Seo KW and Song KH [26] reported that the relative abundance of
Firmicutes was significantly lower in the obese group than that in the lean group, but no difference was
observed in Bacteroidota. In a recent study, Macedo HT, Rentas MF, Vendramini THA, Macegoza MV,
Amaral AR, Jeremias JT, de Carvalho Balieiro JC, Pfrimer K, Ferriolli E and Pontieri CFF [20] showed a
significantly higher proportion of Bacteroidota and a lower proportion of Firmicutes in obese dogs than
that in lean dogs, which is similar to our findings. Therefore, a high relative abundance of Firmicutes, low
relative abundance of Bacteroidota, and a high Firmicutes/Bacteroidota ratio might not be reliable
biomarkers to identify obese dogs.

Butyrate is a major short-chain fatty acid in the gut environment because it is the main fuel for epithelial
cells in the colon and exerts anti-inflammatory and anti-carcinogenic effects [45]; therefore, butyrate-
producing bacteria are considered important microbiota [46]. In this study, Faecalibacterium and
Lachnospiraceae NK4A136 group were enriched in the HBCS group. Faecalibacterium is the predominant
genus belonging to Clostridial cluster IV and can produce butyrate and D-lactate by degrading a wide
range of carbohydrate substrates [47]. Previous studies have reported conflicting results regarding the
relative abundance of Faecalibacterium between obese and lean humans, such as a lower [48, 49] or
higher [50–52] abundance in obese humans. In the fecal microbiome of dogs, there are few reports on the
relationship between members of Faecalibacterium and obesity. Recently, Macedo HT, Rentas MF,
Vendramini THA, Macegoza MV, Amaral AR, Jeremias JT, de Carvalho Balieiro JC, Pfrimer K, Ferriolli E
and Pontieri CFF [20] reported that the relative abundance of Faecalibacterium was significantly higher in
obese dogs than that in lean dogs, and it decreased significantly after the obese group underwent a
weight loss program. Thus, further studies are needed to confirm the role of Faecalibacterium in dogs
with obesity. Lachnospiraceae NK4A136 group is a potential butyrate producer [53], and Companys J,
Gosalbes MJ, Pla-Pagà L, Calderón-Pérez L, Llauradó E, Pedret A, Valls RM, Jiménez-Hernández N,
Sandoval-Ramirez BA and Del Bas JM [52] previously reported that the relative abundance of the
Lachnospiraceae NK4A136 group was enriched in obese humans, which is supported by our results.
Hence, the presence of butyrate-producing bacteria in the fecal microbiome of dogs may be related to
obesity. The abundance of the family Ruminococcaceae was higher in the HBCS group. Many genera
under Ruminococcaceae belong to cluster IV and possess enzymes that degrade recalcitrant
carbohydrate sources and produce acetate as a major fermentation product that can be used as an
energy source by the host [47]. In addition, four major genera related to acetate production as the major
fermentation product, specifically Fournierella [54], Catenibacterium [55], Megamonas [56], and
unclassified genera within Ruminococcaceae [47], were enriched in the HBCS group. Faecalibacterium
prausnitzii, which is a major species belonging to Faecalibacterium, requires acetate for growth [57].
Therefore, the increase in Faecalibacterium may be linked to the high abundance of acetate-producing
bacteria. Furthermore, several studies have suggested a relationship between Sutterella spp. and obesity

                                                  Page 9/25
in humans [58, 59]. A previous study reported a higher relative abundance of Sutterella in obese dogs [20],
which is supported by our results. Thus, we speculate that the increase in Sutterella may be related to
obesity in dogs. Serena C, Ceperuelo-Mallafré V, Keiran N, Queipo-Ortuño MI, Bernal R, Gomez-Huelgas R,
Urpi-Sarda M, Sabater M, Pérez-Brocal V and Andrés-Lacueva C [60] revealed that obesity in humans is
associated with elevated levels of circulating succinate concomitant with specific changes in succinate-
related microbiota, especially higher succinate producers (Prevotellaceae and Veillonellaceae) and lower
succinate consumers (Odoribacteraceae and Clostridiaceae). In the present study, there were no
significant differences in these families between the HBCS and NBCS groups. Instead, we found that the
relative abundances of Bacteroides [61, 62] and Mucispirillum [63], which are considered to be succinate
producers [63], and that of Phascolarctobacterium, which is a propionate producer through the succinate
pathway, were higher in the HBCS group at the genus level. We thus suggest that obese dogs may have a
higher relative abundance of succinate producers and utilizers than lean dogs.

The alteration in fecal microbiota composition is linked to changes in microbial functional features. In
this study, the difference in BCS did not affect the overall predicted functional features, as expected from
the lack of differences in the overall bacterial diversity (Table 3). An excessive energy supply to the host
by increasing carbohydrate metabolism is an important factor in obese animals [64]. In addition, several
studies have reported an increase in the functional profile of carbohydrate metabolism in obese humans
[65, 66]. Therefore, we expected an increased relative abundance of functional profiles of carbohydrate
metabolism in the HBCS group. However, carbohydrate metabolism among the top 10 KEGG categories
showed a significant increase only in the NBCS group. Some KEGG modules involved in carbohydrate
metabolism (galactose degradation [M00632], nucleotide sugar biosynthesis [M00549], and D-
galactonate degradation [M00552]) were significantly enriched in the NBCS group. The unexpected
increase in carbohydrate metabolism can be explained by two possible reasons: 1) the functional profile
was predicted based on DNA; therefore, the RNA-based results (i.e.. RNA sequencing) in carbohydrate
metabolism may be different from this result, and 2) the lack of a dog microbiome-specific PICRUSt2
database. In this study, some KEGG modules related to B vitamins biosynthesis (cobalamin biosynthesis
[M00122] and pyridoxal-P biosynthesis [M00124]) and amino acid biosynthesis (lysine biosynthesis
[M00030], serine biosynthesis [M00020], and tryptophan biosynthesis [M00023]) were enriched in the
HBCS group. The microbiota produces many types of micronutrients that play an important role in host
homeostasis [67]. B vitamins acts as an enzymatic cofactor or a precursor of cofactors [68] and can be
exchanged and shared between microbiota, thereby stabilizing the gut microbiota by cross-feeding [69]. A
recent study that investigated the effects of co-culturing butyrate producers revealed that
Faecalibacterium prausnitzii, a major species belonging to Faecalibacterium, requires several B vitamins
and tryptophan for growth (auxotrophs) [70]. They also demonstrated that co-culturing Faecalibacterium
prausnitzii with species belonging to Lachnospiraceae, which are prototrophic (capable of de novo
synthesis) for all amino acids and several vitamins, stimulates the growth of Faecalibacterium prausnitzii
through cross-feeding [70]. In this study, the HBCS group showed a higher relative abundance of butyrate
producers (Faecalibacterium and Lachnospiraceae NK4A136 group) and obesity-related genera
(Lachnospiraceae NK4A136 group and Sutterella). Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V,

                                                 Page 10/25
Mardis ER and Gordon JI [25] previously suggested that the microbiota of obese mice may be more
efficient in energy harvesting than that of lean mice. Therefore, we speculate that obesity in dogs may be
influenced by the enhancement of energy-harvesting efficiency through cross-feeding among fecal
microbiota, especially among auxotrophs and prototrophs.

We also compared the overall interactions within the major fecal bacterial genera belonging to each BCS
group. Unclassified genera within Erysipelotrichaceae and Phascolarctobacterium were the keystone
genera within exclusive networks in the HBCS group. Previous studies have reported that several species
belonging to the family Erysipelotrichaceae are associated with obesity in mouse models [71] and
humans [72]. Additionally, an association between Erysipelotrichaceae and host lipid metabolism has
been suggested [73]. In the present study, UG_Erysipelotrichaceae had an exclusive positive interaction
with five major genera, namely Allobaculum, Fusobacterium, Lachnospiraceae NK4A136 group,
Megamonas, and Phascolarctobacterium. Among these genera, Lachnospiraceae NK4A136 group [52]
and Megamonas [74], which were higher in abundance in the HBCS group, were also positively
associated with obesity and weight gain, respectively. Phascolarctobacterium consumes succinate and
produces propionate during carbohydrate fermentation [75]. Succinate is not only a metabolic end-
product of some bacteria but also a primary cross-feeding metabolite in microbial propionate synthesis
[76]. Phascolarctobacterium showed a strong positive correlation with 10 major genera, whereof three
genera, namely Fusobacterium [77], Bacteroides [61, 62], and Mucispirillum [63] are considered to be
succinate producers. Interestingly, Muribaculaceae, which showed no difference in the relative abundance
between the HBCS and NBCS groups, exhibited exclusively negative correlations with four major genera,
namely the Lachnospiraceae NK4A136 group, Megamonas, Phascolarctobacterium, and Prevotellaceae
Ga6A1 group, which had significantly higher abundances in the HBCS group. Although little is known
about the relationship between Muribaculaceae and obesity, a previous study revealed the resistance
effect of Muribaculaceae in lean mice fed with a high-fat diet [78]. Therefore, we speculate that
Muribaculaceae may play an anti-obesity role in dogs because of its negative correlation with obesity-
related genera. In the exclusive network of the NBCS group, Bifidobacterium and Negativibacillus were the
keystone genera. Bifidobacterium is well known as an important probiotic that exerts beneficial health
effects, such as regulating microbial homeostasis, reducing intestinal lipopolysaccharide levels, and
improving the mucosal barrier function [79–82]. In addition, several studies have reported that the relative
abundance of Bifidobacterium is lower in mice fed a high-fat diet [83, 84]. In this study, Bifidobacterium
exhibited a strong negative correlation with Peptoclostridium and Blautia, which have previously been
associated with obesity in dogs [20, 85], and a strong positive correlation with Lactobacillus, Sellimonas,
and Negativibacillus. A previous study suggested that Lactobacillus may serve as a potential probiotic
for dogs [27] and Sellimonas plays an important role in gut recovery homeostasis [86]. Therefore, the co-
occurrence interactions among potential probiotics and mutually exclusive interactions among obesity-
related genera might have maintained the BW of the NBCS group.

The majority of previous studies investigating the fecal microbiome in obese dogs focused on the impact
of dietary changes or physical training to reduce weight rather than on natural differences in the
microbiome and microbial interactions. In this study, we investigated essential differences in the fecal
                                                 Page 11/25
microbiome and their interactions between normal and obese dogs. However, we found several results,
especially regarding microbiota related to succinate metabolism and functional profiles of carbohydrate
metabolism, that were inconsistent with the results of previous studies on obese mice and humans.
Therefore, to better understand obesity in dogs, future studies should be conducted to further investigate
the differences in microbial functional metabolism and metabolites between normal and obese dogs.

Conclusions
Differences in BCS did not affect the overall distribution of the microbiota and predicted functional
features. However, several genera related to short-chain fatty acid production were enriched in the HBCS
group. In the predicted microbial functional profiles, some KEGG modules related to B vitamins and
amino acid biosyntheses were higher in the HBCS group, indicating enhanced energy efficiency through
cross-feeding between auxotrophs and prototrophs. In addition, specific co-occurrence and mutually
exclusive interactions among the different genera were detected among the HBCS and NBCS groups.
Further studies are needed to investigate how the specific networks can be interpreted in the context of
fecal fermentation characteristics and obesity in dogs.

Abbreviations
BW
body weight
BCS
body condition score
KEGG
kyoto encyclopedia of genes and genomes
AAFCO
association of american feed control officials
QIIME2
quantitative insights into microbial ecology 2
ASV
amplicon sequence variant
ACE
abundance-based coverage estimator
PCoA
principal coordinate analysis
BIOM
biological observation matrix
PICRUSt2
phylogenetic investigation of communities by reconstruction of unobserved state
PCA

                                                 Page 12/25
principal component analysis
CoDiNA
co-expression differential network analysis
PERMANOVA
permutational multivariate analysis of variance
OTU
operational taxonomic unit

Declarations
Ethics approval and consent to participate

The animal study was reviewed and approved by the National Institute of Animal Science in Korea. This
experiment was conducted in accordance with the methods approved by the Animal Care and Use
Committee National Institute of Animal (NIAS-2019-370).

Consent for publication

Not applicable.

Availability of data and materials

The raw 16S rRNA sequences have been deposited in the NCBI SRA ubder BioProject PRJNA818097. All
analyzed microbioal datasets in the present study are available from the corresponding author on
reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was carried out with the support of the Cooperative Research Program of Center for
Companion Animal Research (Project No. PJ013984012022), Rural Development Administration, Korea.

This study was supported by the 2022 RDA Fellowship Program of the National Institute of Animal
Science, Rural Development Administration, Korea.

Authors' contributions

The study design was performed by JS and KK. KS, HC, and JC conducted the experiment and analyzed
the animal data. HK and JS did the sequencing based analysis. All the visualization of data and
statistical analyses were performed by HK. HK and JK drafted the manuscript. KK reviewed and edited the
manuscript. All authors read and approved the final manuscript.

                                                  Page 13/25
Acknowledgements

Not applicable

References
  1. Houpt KA, Goodwin D, Uchida Y, Baranyiová E, Fatjó J, Kakuma Y: Proceedings of a workshop to
    identify dog welfare issues in the US, Japan, Czech Republic, Spain and the UK. Appl Anim Behav Sci
    2007, 106:221–233.
  2. Mizelle HL, Edwards TC, Montant J-P: Abnormal cardiovascular responses to exercise during the
     development of obesity in dogs. Am J Hypertens 1994, 7:374–378.
  3. Bach JF, Rozanski EA, Bedenice D, Chan DL, Freeman LM, Lofgren JL, Oura TJ, Hoffman AM:
    Association of expiratory airway dysfunction with marked obesity in healthy adult dogs. Am J Vet
    Res 2007, 68:670–675.
  4. Mattheeuws D, Rottiers R, Kaneko J, Vermeulen A: Diabetes mellitus in dogs: relationship of obesity
     to glucose tolerance and insulin response. Am J Vet Res 1984, 45:98–103.
  5. Kealy RD, Lawler DF, Ballam JM, Mantz SL, Biery DN, Greeley EH, Lust G, Segre M, Smith GK, Stowe
     HD: Effects of diet restriction on life span and age-related changes in dogs. J Am Vet Med Assoc
    2002, 220:1315–1320.
  6. German AJ: The growing problem of obesity in dogs and cats. J Nutr 2006, 136:1940S-1946S.
  7. Edney A, Smith P: Study of obesity in dogs visiting veterinary practices in the United Kingdom. Vet
    Rec 1986, 118:391–396.
  8. Mason E: Obesity to pet dogs. Vet Rec 1970, 86:612–616.
  9. Kienzle E, Bergler R, Mandernach A: A comparison of the feeding behavior and the human–animal
     relationship in owners of normal and obese dogs. J Nutr 1998, 128:2779S-2782S.
10. Burkholder WJ: Precision and practicality of methods assessing body composition of dogs and cats.
    Compend Contin Educ Pract Vet 2001, 23:1–10.
11. Pendergrass PB, Bartley CM, Nagy F, Ream LJ, Stuhlman R: A rapid method for determining normal
    weights of medium-to‐large mongrel dogs. J Small Anim Pract 1983, 24:269–276.
12. Laflamme D: Development and validation of a body condition score system for dogs. Canine Pract
    (Santa Barbara, Calif: 1990)(USA) 1997.
13. Huang Z, Pan Z, Yang R, Bi Y, Xiong X: The canine gastrointestinal microbiota: early studies and
    research frontiers. Gut Microbes 2020, 11:635–654.
14. Li Q, Lauber CL, Czarnecki-Maulden G, Pan Y, Hannah SS: Effects of the dietary protein and
    carbohydrate ratio on gut microbiomes in dogs of different body conditions. mBio 2017, 8:e01703-
    01716.
15. Xu J, Verbrugghe A, Lourenço M, Cools A, Liu DJ, Van de Wiele T, Marzorati M, Eeckhaut V, Van
    Immerseel F, Vanhaecke L: The response of canine faecal microbiota to increased dietary protein is

                                                Page 14/25
influenced by body condition. BMC Vet Res 2017, 13:1–11.
16. Moinard A, Payen C, Ouguerram K, André A, Hernandez J, Drut A, Biourge VC, Suchodolski JS,
   Flanagan J, Nguyen P: Effects of high-fat diet at two energetic levels on fecal microbiota, colonic
   barrier, and metabolic parameters in dogs. Front Vet Sci 2020, 7:699.
17. Handl S, German AJ, Holden SL, Dowd SE, Steiner JM, Heilmann RM, Grant RW, Swanson KS,
    Suchodolski JS: Faecal microbiota in lean and obese dogs. FEMS Microbiol Ecol 2013, 84:332–343.
18. Salas-Mani A, Jeusette I, Castillo I, Manuelian CL, Lionnet C, Iraculis N, Sanchez N, Fernández S,
   Vilaseca L, Torre C: Fecal microbiota composition changes after a BW loss diet in Beagle dogs. J
   Anim Sci 2018, 96:3102–3111.
19. Sanchez SB, Pilla R, Sarawichitr B, Gramenzi A, Marsilio F, Steiner JM, Lidbury JA, Woods GR,
    German AJ, Suchodolski JS: Fecal microbiota in client-owned obese dogs changes after weight loss
   with a high-fiber-high-protein diet. PeerJ 2020, 8:e9706.
20. Macedo HT, Rentas MF, Vendramini THA, Macegoza MV, Amaral AR, Jeremias JT, de Carvalho
   Balieiro JC, Pfrimer K, Ferriolli E, Pontieri CFF: Weight-loss in obese dogs promotes important shifts
   in fecal microbiota profile to the extent of resembling microbiota of lean dogs. Anim Microbiome
   2022, 4:1–13.
21. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe
    BA, Affourtit JP: A core gut microbiome in obese and lean twins. Nat 2009, 457:480–484.
22. Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, Almeida M, Arumugam M, Batto J-M,
    Kennedy S: Richness of human gut microbiome correlates with metabolic markers. Nat 2013,
   500:541–546.
23. Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI: Obesity alters gut microbial
    ecology. Proc Natl Acad Sci 2005, 102:11070–11075.
24. Ley RE, Turnbaugh PJ, Klein S, Gordon JI: Human gut microbes associated with obesity. Nat 2006,
    444:1022–1023.
25. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI: An obesity-associated gut
    microbiome with increased capacity for energy harvest. Nat 2006, 444:1027–1031.
26. Park HJ, Lee SE, Kim HB, Isaacson R, Seo KW, Song KH: Association of obesity with serum leptin,
   adiponectin, and serotonin and gut microflora in beagle dogs. J Vet Intern Med 2015, 29:43–50.
27. Grześkowiak Ł, Endo A, Beasley S, Salminen S: Microbiota and probiotics in canine and feline
   welfare. Anaerobe 2015, 34:14–23.
28. Lee W-J, Hase K: Gut microbiota–generated metabolites in animal health and disease. Nat Chem Biol
   2014, 10:416–424.
29. Simpson J, Martineau B, Jones W, Ballam J, Mackie RI: Characterization of fecal bacterial
   populations in canines: effects of age, breed and dietary fiber. Microb Ecol 2002, 44:186–197.
30. AAFCO: Association of American Feed Control Officials. The Official Publication of the Association of
    American Feed Control Officials 2019.

                                                Page 15/25
31. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glöckner FO: Evaluation of general
    16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity
   studies. Nucleic Acids Res 2013, 41:e1-e1.
32. Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J
   2011, 17:10–12.
33. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ,
   Arumugam M, Asnicar F: Reproducible, interactive, scalable and extensible microbiome data science
   using QIIME 2. Nat Biotechnol 2019, 37:852–857.
34. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO: The SILVA
   ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids
   Res 2012, 41:D590-D596.
35. Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, Huttenhower C, Langille MG:
   PICRUSt2 for prediction of metagenome functions. Nat Biotechnol 2020, 38:685–688.
36. Tang Y, Horikoshi M, Li W: ggfortify: unified interface to visualize statistical results of popular R
   packages. R J 2016, 8:474.
37. Watts SC, Ritchie SC, Inouye M, Holt KE: FastSpar: rapid and scalable correlation estimation for
   compositional data. Bioinform 2019, 35:1064–1066.
38. Friedman J, Alm EJ: Inferring correlation networks from genomic survey data. 2012.
39. Morselli Gysi D, de Miranda Fragoso T, Zebardast F, Bertoli W, Busskamp V, Almaas E, Nowick K:
    Whole transcriptomic network analysis using co-expression differential network analysis (CoDiNA).
   PLoS One 2020, 15:e0240523.
40. Coelho LP, Kultima JR, Costea PI, Fournier C, Pan Y, Czarnecki-Maulden G, Hayward MR, Forslund SK,
    Schmidt TSB, Descombes P: Similarity of the dog and human gut microbiomes in gene content and
   response to diet. Microbiome 2018, 6:1–11.
41. Krajmalnik-Brown R, Ilhan ZE, Kang DW, DiBaise JK: Effects of gut microbes on nutrient absorption
   and energy regulation. Nutr Clin Pract 2012, 27:201–214.
42. Schwiertz A, Taras D, Schäfer K, Beijer S, Bos NA, Donus C, Hardt PD: Microbiota and SCFA in lean
   and overweight healthy subjects. Obe 2010, 18:190–195.
43. Duncan SH, Lobley G, Holtrop G, Ince J, Johnstone A, Louis P, Flint HJ: Human colonic microbiota
    associated with diet, obesity and weight loss. Int J Obe 2008, 32:1720–1724.
44. Chun JL, Ji SY, Lee SD, Lee YK, Kim B, Kim KH: Difference of gut microbiota composition based on
    the body condition scores in dogs. J Anim Sci Technol 2020, 62:239.
45. Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost F, Brummer RJ: The role of butyrate on
    colonic function. Aliment Pharmacol Ther 2008, 27:104–119.
46. Louis P, Flint HJ: Diversity, metabolism and microbial ecology of butyrate-producing bacteria from
   the human large intestine. FEMS Microbiol Lett 2009, 294:1–8.

                                                  Page 16/25
47. Duncan SH, Louis P, Flint HJ: Cultivable bacterial diversity from the human colon. Lett Appl Microbiol
    2007, 44:343–350.
48. Yasir M, Angelakis E, Bibi F, Azhar EI, Bachar D, Lagier JC, Gaborit B, Hassan AM, Jiman-Fatani AA,
    Alshali KZ, et al: Comparison of the gut microbiota of people in France and Saudi Arabia. Nutr
   Diabetes 2015, 5:e153-e153.
49. Andoh A, Nishida A, Takahashi K, Inatomi O, Imaeda H, Bamba S, Kito K, Sugimoto M, Kobayashi T:
    Comparison of the gut microbial community between obese and lean peoples using 16S gene
   sequencing in a Japanese population. J Clin Biochem Nutr 2016, advpub.
50. Chávez-Carbajal A, Nirmalkar K, Pérez-Lizaur A, Hernández-Quiroz F, Ramírez-del-Alto S, García-Mena
   J, Hernández-Guerrero C: Gut Microbiota and Predicted Metabolic Pathways in a Sample of Mexican
   Women Affected by Obesity and Obesity Plus Metabolic Syndrome. Int J Mol Sci 2019, 20:438.
51. Salah M, Azab M, Ramadan A, Hanora A: New Insights on Obesity and Diabetes from Gut
   Microbiome Alterations in Egyptian Adults. Omics 2019, 23:477–485.
52. Companys J, Gosalbes MJ, Pla-Pagà L, Calderón-Pérez L, Llauradó E, Pedret A, Valls RM, Jiménez-
   Hernández N, Sandoval-Ramirez BA, Del Bas JM: Gut microbiota profile and its association with
   clinical variables and dietary intake in overweight/obese and lean subjects: A cross-sectional study.
   Nutrients 2021, 13:2032.
53. Stadlbauer V, Engertsberger L, Komarova I, Feldbacher N, Leber B, Pichler G, Fink N, Scarpatetti M,
   Schippinger W, Schmidt R, Horvath A: Dysbiosis, gut barrier dysfunction and inflammation in
   dementia: a pilot study. BMC Geriatr 2020, 20:248.
54. Togo AH, Durand G, Khelaifia S, Armstrong N, Robert C, Cadoret F, Di Pinto F, Delerce J, Levasseur A,
   Raoult D: Fournierella massiliensis gen. nov., sp. nov., a new human-associated member of the
   family Ruminococcaceae. Int J Syst Evol Microbiol 2017, 67:1393–1399.
55. Kageyama A, Benno Y: Catenibacterium mitsuokai gen. nov., sp. nov., a gram-positive anaerobic
   bacterium isolated from human faeces. Int J Syst Evol Microbiol 2000, 50:1595–1599.
56. Chevrot R, Carlotti A, Sopena V, Marchand P, Rosenfeld E: Megamonas rupellensis sp. nov., an
   anaerobe isolated from the caecum of a duck. Int J Syst Evol Microbiol 2008, 58:2921–2924.
57. Duncan SH, Hold GL, Harmsen HJ, Stewart CS, Flint HJ: Growth requirements and fermentation
   products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii
   gen. nov., comb. nov. Int J Syst Evol Microbiol 2002, 52:2141–2146.
58. Hou Y-P, He Q-Q, Ouyang H-M, Peng H-S, Wang Q, Li J, Lv X-F, Zheng Y-N, Li S-C, Liu H-L: Human gut
   microbiota associated with obesity in Chinese children and adolescents. BioMed Res Int 2017, 2017.
59. Gao X, Zhang M, Xue J, Huang J, Zhuang R, Zhou X, Zhang H, Fu Q, Hao Y: Body mass index
   differences in the gut microbiota are gender specific. Front Microbiol 2018, 9:1250.
60. Serena C, Ceperuelo-Mallafré V, Keiran N, Queipo-Ortuño MI, Bernal R, Gomez-Huelgas R, Urpi-Sarda
   M, Sabater M, Pérez-Brocal V, Andrés-Lacueva C: Elevated circulating levels of succinate in human
   obesity are linked to specific gut microbiota. ISME J 2018, 12:1642–1657.

                                                Page 17/25
61. Miller TL: The pathway of formation of acetate and succinate from pyruvate by Bacteroides
    succinogenes. Arch Microbiol 1978, 117:145–152.
62. Johnson JL, Moore W, Moore LV: Bacteroides caccae sp. nov., Bacteroides merdae sp. nov., and
    Bacteroides stercoris sp. nov. isolated from human feces. Int J Syst Evol Microbiol 1986, 36:499–
   501.
63. Loy A, Pfann C, Steinberger M, Hanson B, Herp S, Brugiroux S, Neto JCG, Boekschoten MV, Schwab C,
    Urich T, et al: Lifestyle and Horizontal Gene Transfer-Mediated Evolution of Mucispirillum schaedleri,
   a Core Member of the Murine Gut Microbiota. mSystems 2017, 2:e00171-00116.
64. Wang T-Y, Zhang X-Q, Chen A-L, Zhang J, Lv B-H, Ma M-H, Lian J, Wu Y-X, Zhou Y-T, Ma C-C: A
   comparative study of microbial community and functions of type 2 diabetes mellitus patients with
   obesity and healthy people. Appl Microbiol Biotechnol 2020, 104:7143–7153.
65. Pedrogo DAM, Jensen MD, Van Dyke CT, Murray JA, Woods JA, Chen J, Kashyap PC, Nehra V: Gut
   microbial carbohydrate metabolism hinders weight loss in overweight adults undergoing lifestyle
   intervention with a volumetric diet. In Mayo Clinic Proceedings. Elsevier; 2018: 1104–1110.
66. Del Chierico F, Abbatini F, Russo A, Quagliariello A, Reddel S, Capoccia D, Caccamo R, Ginanni
    Corradini S, Nobili V, De Peppo F: Gut microbiota markers in obese adolescent and adult patients:
   age-dependent differential patterns. Front Microbiol 2018, 9:1210.
67. Agus A, Clément K, Sokol H: Gut microbiota-derived metabolites as central regulators in metabolic
    disorders. Gut 2021, 70:1174–1182.
68. Madigan MT, Martinko JM, Parker J: Brock biology of microorganisms. Pearson Prentice Hall Upper
    Saddle River, NJ; 2006.
69. Sharma V, Rodionov DA, Leyn SA, Tran D, Iablokov SN, Ding H, Peterson DA, Osterman AL, Peterson
    SN: B-vitamin sharing promotes stability of gut microbial communities. Front Microbiol 2019:1485.
70. Soto-Martin EC, Warnke I, Farquharson FM, Christodoulou M, Horgan G, Derrien M, Faurie J-M, Flint
    HJ, Duncan SH, Louis P: Vitamin biosynthesis by human gut butyrate-producing bacteria and cross-
   feeding in synthetic microbial communities. mBio 2020, 11:e00886-00820.
71. Woting A, Pfeiffer N, Loh G, Klaus S, Blaut M: Clostridium ramosum promotes high-fat diet-induced
    obesity in gnotobiotic mouse models. MBio 2014, 5:e01530-01514.
72. Zhang H, DiBaise JK, Zuccolo A, Kudrna D, Braidotti M, Yu Y, Parameswaran P, Crowell MD, Wing R,
    Rittmann BE: Human gut microbiota in obesity and after gastric bypass. Proc Natl Acad Sci 2009,
   106:2365–2370.
73. Kaakoush NO: Insights into the role of Erysipelotrichaceae in the human host. Front Cell Infect
   Microbiol 2015, 5:84.
74. Kieler IN, Shamzir Kamal S, Vitger AD, Nielsen DS, Lauridsen C, Bjornvad CR: Gut microbiota
    composition may relate to weight loss rate in obese pet dogs. Vet Med Sci 2017, 3:252–262.
75. Del Dot T, Osawa R, Stackebrandt E: Phascolarctobacterium faecium gen. nov, spec. nov., a novel
    taxon of the Sporomusa group of bacteria. Syst Appl Microbiol 1993, 16:380–384.

                                                Page 18/25
76. Fernández-Veledo S, Vendrell J: Gut microbiota-derived succinate: Friend or foe in human metabolic
    diseases? Rev Endocr Metab Disord 2019, 20:439–447.
77. McDonald NC, White RL: Reduction of fumarate to succinate mediated by Fusobacterium varium.
   Appl Biochem Biotechnol 2019, 187:163–175.
78. Cao W, Chin Y, Chen X, Mi Y, Xue C, Wang Y, Tang Q: The role of gut microbiota in the resistance to
   obesity in mice fed a high fat diet. Int J Food Sci Nutr 2020, 71:453–463.
79. Sela D, Price N, Mills D: Bifidobacteria: genomics and molecular aspects. Book Chap. 2010:45–70.
80. Wang Z, Xiao G, Yao Y, Guo S, Lu K, Sheng Z: The role of bifidobacteria in gut barrier function after
    thermal injury in rats. J Trauma Acute Care Surg 2006, 61:650–657.
81. Griffiths EA, Duffy LC, Schanbacher FL, Qiao H, Dryja D, Leavens A, Rossman J, Rich G, Dirienzo D,
    Ogra PL: In vivo effects of bifidobacteria and lactoferrin on gut endotoxin concentration and
   mucosal immunity in Balb/c mice. Dig Dis Sci 2004, 49:579–589.
82. Wang Z-T, Yao Y-M, Xiao G-X, Sheng Z-Y: Risk factors of development of gut-derived bacterial
    translocation in thermally injured rats. World J Gastroenterol 2004, 10:1619.
83. Cani PD, Amar J, Iglesias MA, Poggi M, Knauf C, Bastelica D, Neyrinck AM, Fava F, Tuohy KM, Chabo
    C: Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes 2007, 56:1761–1772.
84. Cani PD, Neyrinck AM, Fava F, Knauf C, Burcelin RG, Tuohy KM, Gibson G, Delzenne NM: Selective
    increases of bifidobacteria in gut microflora improve high-fat-diet-induced diabetes in mice through a
   mechanism associated with endotoxaemia. Diabetologia 2007, 50:2374–2383.
85. Thomson P, Santibáñez R, Rodríguez-Salas C, Flores-Yañez C, Garrido D: Differences in the
    composition and predicted functions of the intestinal microbiome of obese and normal weight adult
   dogs. PeerJ 2022, 10:e12695.
86. Muñoz M, Guerrero-Araya E, Cortés-Tapia C, Plaza-Garrido A, Lawley TD, Paredes-Sabja D:
   Comprehensive genome analyses of Sellimonas intestinalis, a potential biomarker of homeostasis
   gut recovery. Microb Genom 2020, 6.
87. Council NR: Nutrient requirements of dogs and cats. National Academies Press; 2006.
88. Blondel VD, Guillaume J-L, Lambiotte R, Lefebvre E: Fast unfolding of communities in large networks.
    J Stat Mech: Theory and Experiment 2008, 2008:10008.

Tables

                                                Page 19/25
Table 1
                        Analyzed chemical composition of the experimental diet.
Chemical composition (%DM or as stated)

DM, %                                                                       88.0

CP                                                                          37.5

EE                                                                          22.7

CF                                                                          3.41

NFE1                                                                        21.6

Ash                                                                         14. 8

Ca                                                                          1.00

P                                                                           0.80

ME2, Mcal/kg                                                                3.68

DM, dry matter; CP, crude protein; EE, ether extract; CF, crude fiber; NFE, nitrogen-free extract; Ca,
calcium; P, phosphorus; ME, metabolizable energy.

1NFE%    = 100 - (CP + EE + CF + Ash + moisture)

2Estimated    based on National Research Council [87].

                                              Table 2
        Differences in body weight and body condition score between obeses and normal dogs.
                                 Treatment1

Parameters                       HBCS               NBCS              SEM                P-value

Initial BW, kg                   11.88              10.05             0.590              0.0261

Final BW, kg                     13.10              10.07             0.772              0.0045

BW gain, kg                      1.22               0.02              0.375              0.0287

Initial BCS                      7.0                5.2               0.601              0.0262

Final BCS                        8.0                5.0               0.365              < 0.0001

SEM, standard error of the mean.

1HBCS,  high body condition score (BCS) group (BCS range: 7–9, 4 males and 2 females, all 2-year-
old); NBCS, normal BCS group (BCS range: 4–6, 4 males and 2 females, all 2-year-old).

                                                   Page 20/25
Table 3
            Differences in alpha diversity measurements of the fecal microbiota by obesity.
                                      Treatments1

Parameters                            HBCS           NBCS           SEM2              P-value

Observed ASVs                         180            174            19.1              0.8220

Chao1 estimates                       181            175            19.5              0.8268

ACE                                   180            175            19.4              0.8344

Evenness                              0.698          0.636          0.0261            0.1044

Shannon’s index                       5.22           4.72           0.251             0.1830

Simpson’s index                       0.934          0.897          0.0181            0.1361

ASV, amplicon sequencing variant; ACE, abundance-based coverage estimator.

1
 HBCS, high body condition score (BCS) group (BCS range: 7–9, 4 males and 2 females, all 2-year-
old); NBCS, normal BCS group (BCS range: 4–6, 4 males and 2 females, all 2-year-old).

2SEM,   standard error of the mean.

                                               Page 21/25
You can also read