- 1Department of Research and Development, Affinity Petcare, L'Hospitalet de Llobregat, Spain
- 2Gastrointestinal Laboratory, Department of Small Animal Clinical Sciences, Texas A&M University, College Station, TX, United States
- 3School of Veterinary Science, Universitat Autònoma de Barcelona, Bellaterra, Spain
The gut microbiome is critical for maintaining host health. In healthy humans, the aging process is one of the main factors modulating the changes in the intestinal microbiota. However, little is known about the relationship between gut health, microbiota, and the aging process in dogs. The present study aims to explore the differences in the intestinal microbiota and intestinal health based on fecal biomarkers in a population of dogs of different ages. The study involved 106 dogs of different breeds aged between 0.2 and 15 years categorized as senior (>7 years; n = 40), adult (2–7 years; n = 50), and junior (<2 years; n = 16). Fecal samples were collected during the same period at the same facilities. The analysis included the following gut health indicators: 16S rRNA gene sequencing to investigate the differences in the fecal microbiota; qPCR to determine the dysbiosis index; fecal short-chain fatty acid concentrations; fecal calprotectin; and immunoglobulin A. Beta diversity analysis revealed a significant difference with a small effect size (p = 0.003; R = 0.087) among age categories based on the unweighted UniFrac metric, but no significance was observed based on the weighted UniFrac metric or Bray–Curtis distances. There were no significant differences in the alpha diversity measures or the fecal dysbiosis index among age categories. Senior dogs had significantly higher relative abundance proportions in phyla Bacteroidota and Pseudomonadota and the genus Faecalibacterium, but not on qPCR analysis. At the family level, Ruminococcaceae, Uncl. Clostridiales.1, Veillonellaceae, Prevotellaceae, Succinivibrionaceae, and Bacteroidaceae abundances were higher in the senior category than in the adult and/or junior categories. Relative proportions, but not concentrations of fecal acetate, were higher in the senior category, while butyrate, isovaleric acid, and valeric acid were lower. The valeric acid concentration was significantly lower in the senior category than in the adult category. Calprotectin and immunoglobulin A levels did not differ significantly across groups. In conclusion, this study observed multiple minor changes in the fecal microbiota composition and the relative amount of short-chain fatty acids in dogs among different age groups, but studies in larger populations representative of all ages are warranted to refine the present results.
1. Introduction
The gut microbiome is defined as the entire collection of microorganisms living in the gastrointestinal tract, dominated by bacteria and complemented by commensal populations of fungi, viruses, archaea, and protists (1). It has the largest microbiota population in the body when compared to other colonized organs. The intestinal microbiota is estimated to be composed of 1013-1014 cells (2), and its genome is 150 times larger than the human genome.
The gut microbiome plays an important role in the physiological and pathological state of its host. It participates in multiple metabolic functions, protects against pathogens, and plays a crucial role in the immune response. These processes are critical to maintaining the health status of the host (3). The composition of the microbiota is crucial to maintaining balanced gut functionality. An imbalanced or disrupted microbiota is related to several intestinal and extraintestinal diseases (4).
There is an increasing number of studies evaluating the different elements modulating the intestinal microbiota (5–8). The main factors can be classified as external and internal, including, for example, diet and pharmacological treatments and age and genetics, respectively. In addition, pathological disorders (e.g., inflammation and type 2 diabetes) may induce deviations in the gut microbiota, causing an imbalance in the gut microbiota, defined as dysbiosis. Among all these factors, age is one of the most important variables that must be considered when studying the evolution and changes of the gut microbiota.
Aging and its relationship with gut health and gut microbiota are currently being explored in mammals, and they are already well-established in humans. Badal et al. (9) summarized the current knowledge of the human gut microbiota, considering composition, function, and metabolic products from the microbiota in relation to aging and lifespan. Changes in the human microbiota related to age have been described and are apparently associated with the host's health scenario. For example, bifidobacteria belonging to members of the genus Bifidobacterium, which are considered a beneficial bacterial group, decrease during the shift from middle age to old age. Contrarily, Clostridium perfringens, lactobacilli, enterococci, and Enterobacteriaceae increase during the aging process (10). In contrast to the human microbiota, bifidobacteria may not play an important role in dog gut health, according to the results reported by Masuoka et al. (10), in which bifidobacteria were only found in half of the youngest dogs and none of the adult dogs sampled, although differences in methodologies need to be considered when interpreting these results.
The existence of a core microbiota in healthy dogs has already been described (11). Moreover, Ziese and Suchodolski (12) associated shifts in the canine fecal microbiota with certain pathologies. Garrigues et al. (13) recently reviewed the development of the gut microbiota during the early stages of canine life, supporting observed changes in bacterial communities from day 2 of age to up to 52 weeks. On day 2 after birth, gut microbiota richness increases, and from day 2 to 21, Bacillota predominance is substituted by a codominance of Bacillota, Bacteroidota, and Fusobacteriota (14). Puppies in the first few weeks of life have immature microbiota, characterized by an increased dysbiosis index (DI), Clostridium difficile abundance, and decreased Clostridium hiranonis compared to adult dogs (15). Approximately after 4–6 months of age, the microbiota resembles that observed in adult dogs and remains largely stable in adulthood (13). However, the relationship between the intestinal microbiota, gut health conditions, and the aging process in older dogs is not well established. Studies in this field are needed first to understand this relationship and then to explore effective strategies to improve the quality of the aging process in dogs from this perspective.
The present study aimed to investigate the differences in the intestinal microbiota and intestinal health based on fecal biomarkers in a population of dogs of different ages, with a special focus on the senior category.
2. Materials and methods
All samples collected from dogs were in accordance with the 2010/63/UE directive, and any additional procedures to the dogs' daily routine were undertaken. Before the study, the protocol was shared with the Affinity Ethical Committee to ensure good practices and the animals' welfare. All dogs were housed at the Affinity Nutrition Center (Barcelona, Spain).
2.1. Animals
The study involved 106 multi-breed dogs (n = 53 male dogs; n = 53 female dogs) aged between 0.2 and 15.0 years (mean age = 6.3 years). The dog breeds that participated in the study were as follows: Beagle, Bichon Maltese, Jack Russell Terrier, Boxer, German Shorthaired Pointer, Poodle, Chihuahua, Cocker Spaniel, Dalmatian, German Shepherd, Spanish Water Dog, Miniature Pinscher, Andalusian Podenco, Pointer, Pomeranian, Miniature, Schnauzer, Shih Tzu, Brittany Spaniel, and Yorkshire Terrier [detailed information is in the Supplementary Table 1]. Animals were grouped into junior (J), adult (A), and senior (S) categories, considering the respective age ranges: up to 2 years, from 2 to 7 years, and over 7 years (Table 1).
Animals lived in kennels in groups of between two and seven individuals, depending on the size of the dogs. Dogs living in pairs shared an area of 15 m2, whereas those living in larger groups shared an area of 57–60 m2. All the animals had free access to the outdoor and indoor parts of their kennels, and they spent between 2 and 4 h per day in a bigger outdoor park where they could socialize with a larger group of dogs.
The water supply was ad libitum. Non-specific diets were manufactured and offered because of the study. Dogs were fed commercially available, dry-complete, and balanced diets adequate to their needs. Furthermore, 54% of dogs were fed a constant diet, whereas the diet was changed for 46% of dogs before and during the collection period. All diets met all requirements for protein, fat, vitamins, and minerals recommended by FEDIAF (16). Most of the dogs in the study (89%) were fed diets formulated with the following composition range: 7–8% moisture, 22–30% protein, 15–21% fat, and 2–5% fiber. The rest of the dogs (11%) were out of range because of the protein, fat, and/or fiber composition [detailed information is in the Supplementary Table 2]. Mineral and vitamin levels were similar in all the diets that were present throughout the study.
The inclusion criteria for enrollment were generally healthy animals and comprised sterilized male and female dogs of all pure breeds and ages. Based on the literature, the exclusion criteria included the presence of clinical disease or any pharmacological treatment potentially interfering with the gut microbiota and/or intestinal health. Specifically, no antibiotics, NSAIDs, PPIs, prebiotics, probiotics, or deworming treatments were administered during or 4 weeks before the sampling. None of the animals participating in this study were fed diets containing probiotics. Dogs with watery, soft, or unformed stools in three consecutive samples were excluded from the study. No clinical signs were compatible with gastrointestinal disorders (such as vomiting, regurgitation, and hyporexia) were observed in any of the animals participating in the study.
2.2. Fecal collection
Samples were collected during 3 months (March, April, and May) without interfering with the dogs' daily routines. Feces were collected through direct observation while the animals were kept in their kennels or outdoor recreational areas. The majority of samples were obtained during the first month, but due to the volume of samples and the physiologically atypical routines of some animals, the collection extended up to 3 months. In all cases, individualizing the animals and/or using cages to collect the samples were avoided.
One stool sample was collected per dog, and it was first scored according to an adapted five-point scale of fecal consistency (0 = watery stool; 25 = soft unformed stool; 50 = soft formed and moist stool retaining shape when being collected; 75 = hard formed stool remaining shaped but soft; and 100 = hard dry stool). The scoring system was adapted from the five-point scale previously described by Strickling et al. (17).
If fecal consistency was equal to or lower than 50, the sample was not collected, and the researchers waited for a higher fecal score. The entire stool or enough quantity to fill the sterile fecal collection tube with 250 ml was taken. External contaminants (stones, grass, sand, etc.) were discarded. All samples were processed within 4 h after the deposition. All the materials used, from sample collection to separation into aliquots, were sterile and meant for single use.
Feces were divided into two different aliquots according to the analysis to be conducted. For the microbiota analysis, between 0.5 and 1 g of feces were prepared in a 1.5-ml sterile, RNAse-free tube. For short-chain fatty acid (SCFA) analysis, calprotectin (cCP) and immunoglobulin A (IgA) were filled in one sterile stool tube with 5 g of feces. Both tubes were stored at −80°C before being sent and analyzed at the Small Animal Clinical Sciences Department at Texas A&M University (College Station, TX).
2.3. Fecal biomarkers
2.3.1. Microbiota analysis
Illumina sequencing of the bacterial 16S rRNA genes was performed using primers 515F (5′-GTGYCAGCMGCCGCGGTAA) (18) to 806RB (5′-GGACTACNVGGGTWTCTAAT) (19) at the MR DNA laboratory (Shallowater, TX, USA).
Sequences were processed and analyzed using the Quantitative Insights Into Microbial Ecology 2 (QIIME 2) (20) v 2021.8 pipeline. The raw sequences were uploaded to the NCBI Sequence Read Archive under the BioProject identification PRJNA901473. In brief, the sequences were demultiplexed, and the ASV table was created using DADA2 (21). Before downstream analysis, sequences assigned as chloroplast, mitochondria, and low abundance ASVs, containing < 0.01% of the total reads in the dataset, were removed. All samples were rarefied to even sequencing depth, based on the lowest read depth of samples, to 21,025 sequences per sample.
Alpha diversity was measured with the Chao1 (richness) and Shannon diversity metrics within QIIME2. Beta diversity was evaluated with the unweighted and weighted phylogeny-based UniFrac (22) distance metric (measures that consider phylogenetic information) and the Bray-Curtis distance metric and visualized using Principal Coordinate Analysis (PCoA) plots, generated within QIIME2.
The dysbiosis index (logDNA) is based on a validated algorithm that considers a panel of eight bacterial groups identified by a Quantitative PCR assay for total bacteria, Faecalibacterium spp., Turicibacter spp., Escherichia coli, Streptococcus spp., Blautia spp., Fusobacterium spp., and Clostridium hiranonis (23). A DI of <0 was defined as normal, a DI between 0 and 2 was defined as mild to moderate microbiota shift, and a DI of >2 was considered significant dysbiosis.
2.3.2. Intestinal health indicators
Concentrations of SCFAs (acetate, propionate, butyrate, isobutyric acid, isovaleric acid, and valeric acid) in feces were measured using a stable isotope dilution gas chromatography-mass spectrometry (GC-MS) assay with some modifications (24). To consider the difference in water content between fecal samples, the final concentrations of fecal SCFAs were adjusted by fecal dry matter (DM) and expressed as μmol/g of fecal DM (25).
The fecal IgA (mg/g) was quantified by using the commercial ELISA kit for canine IgA determination (Bethyl Laboratories, Montgomery, TX, USA). The quantification of the fecal cCP (ng/g) was also analyzed by using the ELISA kit previously validated in dogs (26, 27).
2.4. Statistical analysis
For the statistical analysis of microbiota data, an analysis of similarity (ANOSIM) test within the PRIMER 7 software package (PRIMER-E Ltd., Luton, UK) was performed to analyze significant differences in microbial communities and the size effect (R-values between 0 and 1; a higher R-value indicates a larger size effect) among age groups (J, A, and S). All datasets were tested for normality by performing the Shapiro–Wilk test (JMP Pro 11, SAS Software Inc.). The Kruskal–Wallis test was performed (Prism v. 9, GraphPad Software Inc.), followed by a post-hoc Dunn's multiple comparison tests, to determine the age group differences in bacterial taxa (including phylum, class, order, family, genus, and species). All p-values were adjusted for multiple comparisons using Benjamini and Hochberg's False Discovery Rate (28) at each taxonomic level, and an adjusted p-value of <0.05 (q-value) was considered statistically significant.
For the general statistical analysis, the model included demographic information as an explanatory variable and fecal indicators as response variables. The primary explanatory variable was age, which was grouped into the three different categories described. The potential explanatory variables included the following: body weight, sex, “feeding routine,” “housing,” and “breed.” Two categories were established for feeding routines: rotation (R), if animals were changing the type of diet in periods shorter than 2 weeks, and stable (S), if animals were eating the same diet for at least 4 weeks. No animals had a change of diet in the 2–4 week range. Each animal was only assigned to one category during the whole study. Housing was the variable defined to identify the three different buildings where animals were allocated. Finally, as for the breed, the dogs were categorized into two groups: beagles and non-beagle dogs. The different breeds were not statistically analyzed due to the high number of different breeds but due to the low number of dogs within each breed; however, the dog size was introduced in the statistical analysis to reduce the heterogeneity of dogs participating in the study. The response variables considered in the general statistical analysis were as follows: fecal consistency; SCFA concentration (acetate, propionate, butyrate, isobutyric acid, isovaleric acid, valeric acid, and total SCFA); SCFA relative amounts (acetate, propionate, butyrate, isobutyric acid, and isovaleric acid); calprotectin; immunoglobulin A; alpha diversity (Chao 1 and Shannon diversity); and finally, from qPCR, total bacteria, Faecalibacterium, Turicibacter, Streptococcus, E. coli, Blautia, Fusobacterium, Clostridium hiranonis, Bifidobacterium, Bacteroides, Lactobacillus, and dysbiosis index.
For the general statistical analysis, a first summary statistic was performed, in which quantitative variables were analyzed using the mean and standard deviation. Qualitative variables were tested using relative and absolute frequencies. The existence of differences among age groups was tested by performing the appropriate tests (ANOVA and Kruskal–Wallis test), considering the equality among groups' null hypothesis. The compliance of the application criteria was assessed by performing the Shapiro–Wilk normality test. The relationship between quantitative variables was analyzed using Spearman's correlation.
To analyze the relationship between the age categories and fecal markers, the appropriate linear model was considered, including potential explanatory variables such as body weight, sex, feeding routines, and housing breed. Estimated means (emmeans) for age categories were calculated using the adjusted model. Pairwise post-hoc comparisons were also performed, and the model validation was analyzed by performing a graphical residual analysis. For the general statistical analyses, differences were considered significant with a p-value of <0.05.
3. Results
3.1. Quantitative real-time PCR analysis and dysbiosis index
The results from the qPCR analysis are summarized in Table 2. Bifidobacterium abundance was different between J and A, being higher in A. However, no difference was found when considering S vs. J and A. The calculated dysbiosis index was not different among the age categories (Figure 1). In J, 0 of 16 (0%) dogs had a DI of >0, while in A and S, 8 of 50 (16%) and 5 of 40 (13%) dogs had a DI of >0.
Table 2. Results of the quantitative real-time PCR (abundance expressed as LogDNA) and values of the calculated dysbiosis index (ratio).
3.2. 16s rRNA sequencing
A total of 3,856,065 quality bacterial 16S rRNA sequences were obtained from the 106 fecal samples analyzed. The range count per sample was between 21,030 and 60,827 (median: 33,688 and mean: 36,378). The beta diversity was analyzed by the weighted and unweighted UniFrac distances and the Bray–Curtis distance. The unweighted UniFrac analysis of similarities, which only considers the presence or absence of individual taxa, revealed significant differences in the microbial communities among age categories (PCoA plot shown in Figure 2). However, no significance was observed based on weighted UniFrac metric or Bray–Curtis distances (Table 3). Alpha diversity indices (Chao1 and Shannon) were not significantly different among the age groups (Figures 3, 4).
Figure 2. 3D beta diversity patterns of canine fecal microbiota comparing Adult (red), Junior (blue) and Senior (green) categories based on unweighted UniFrac distance.
Table 3. Results of beta diversity analysis including unweighted UniFrac, weighted UniFrac, and Bray–Curtis distances.
3.2.1. Phylum taxonomic level
When the bacterial relative abundance was studied at the tax level, the major phylum identified was Bacillota, followed by Actinomycetota, Fusobacteriota, Bacteroidota, and Pseudomonadota. The same distribution, considering relative abundance, was described in the three age categories. Significant differences were found for the phylum Bacteroidota and Pseudomonadota, in which S had a higher abundance than A and J (Table 4). However, the corrected q-value did not reach significance.
Table 4. Fecal microbiota composition (% of relative abundance from rarefied data) at the phylum level split by age category.
3.2.2. Family taxonomic level
Falling into the main phylum (Bacillota), Ruminococcaceae, Uncl. Clostridiales.1, Veillonellaceae, and Lachnospiraceae families were significantly different among the age categories (Figures 5–8). Specifically, the abundance of Ruminococcaceae and Uncl. Clostridiales.1 was higher in the S category than in the A category, and Veillonellaceae showed a higher abundance in the S category than in the J category. Lachnospiraceae were significantly higher in the J category than in the A and S categories and did not differ significantly. Within the Bacteroidota phylum, the relative abundance of Bacteroidaceae and Prevotellaceae was different among the age categories; Bacteroidaceae was higher in the S category than in the A and J categories, and Prevotellaceae was higher in the S category than in the A. For Pseudomonadota, the relative abundance of Succinivibrionaceae was higher in the S category than in the A category. Finally, the Bifidobacteriaceae family (Actinomycetota) was significantly higher in the A category than in the J category; however, the adjusted q-value was not significant (Table 5).
Figure 5. A box plot representation of the relative abundance from rarefied data of the family Ruminococcaceae family by age category.
Figure 6. A box plot representation of the relative abundance from rarefied data of the Uncl. Clostridiales.1 family by age category.
Figure 7. A box plot representation of the relative abundance from rarefied data of the Veillonellaceae family by age category.
Figure 8. A box plot representation of the relative abundance from rarefied data of the Lachnospiraceae family by age category.
Table 5. Fecal microbiota composition (% of relative abundance from rarefied data) at the family level split by age group.
3.2.3. Genus taxonomic level
The two predominant genera identified in the dogs' population belonged to the phylum Bacillota, Uncl. Clostridiaceae.1, and Allobaculum. Within Bacillota, Faecalibacterium, and Uncl. Clostridiales.3 showed the highest values in S, being significantly different when compared with A (p = 0.006 and 0.008, respectively). The relative abundance observed in the genus Faecalibacterium was 18.6 times higher in the S category than in the A category. However, the difference in Faecalibacterium was not confirmed by a targeted qPCR. Within this genus, the relative abundance of F. prausnitzii was 18.4 times higher in the S category (median 1.3% of sequences) than in the A category (median 0.07%; p = 0.006 and q = 0.030). Detailed information at the species level is in the Supplementary Table 3. In relation to Phascolarctobacterium abundance, S was higher than the A and J categories (p = 0.011). Although the S category also showed higher values for Ruminococcus when compared with A (p = 0.035), the adjusted q-value did not reach significance. Contrarily, the abundance of Eubacterium, Uncl. Erysipelotrichaceae, and Uncl. Lachnospiraceae was not significantly different when the S category vs. the J and A categories were considered. However, J showed significantly higher values for the three genera when compared with the A category (p = 0.005, p = 0.001, and p = 0.002). Specifically, J showed 8.6 and 7.2 times more relative abundance in Eubacterium and Uncl. Erysipelotrichaceae genera, respectively, compared with A. Considering the Uncl. Peptostreptococcaceae genus, the relative abundance was higher in J than in S and A categories (p = 0.009).
The relative abundance of Bacteroides (phylum Bacteroidota) was 24.3 and 10.4 times higher in the S category than in the J and A categories, respectively (p < 0.001). Prevotella (Prevotellaceae family) was higher in S than in the A category (p = 0.016); however, the adjusted q-value was not significant. The relative abundance of Fusobacterium (phylum Fusobacteriota) was higher in the S category than in the A category when considering only the significance of the p-value (p = 0.045), but this was not confirmed by targeted qPCR. Detailed information at the genus level is in the Supplementary Table 4.
3.3. Fecal SCFAs, IgA, and cCP
Valeric acid was significantly lower in the S category than in the A category. However, the concentrations for total SCFA were numerically higher (but not significant) in the S category and lower in the A and J categories (Table 6). Acetate, butyrate, isovaleric acid, and valeric acid reached significant results when they were analyzed as relative percentages (Table 6). In the post-hoc analysis, the S category showed higher values than A for acetate (p < 0.001) and lower values for butyrate, isovaleric acid, and valeric acid (p = 0.001, p = 0.045, and p = 0.026, respectively). For isovaleric acid, the S category had lower values than J. Finally, no significant differences among age categories were found for the fecal IgA and cCP concentrations (Table 7).
Table 6. Fecal SCFA values expressed as concentration (μmol/g DM) and relative percentages (%) reported by age category.
3.4. Correlations among qPCR, IgA, cCP, and SCFA
Considering significant correlations (rho) equal to or >0.3, Faecalibacterium abundance correlated positively with propionate and valeric acid concentrations, while Fusobacterium correlated positively with isobutyric concentration. Clostridium hiranonis correlated positively with isovaleric acid percentage, and Bifidobacterium correlated positively with valeric acid concentration. This genus correlated positively with the acetate relative percentage, while the correlation with the butyrate relative percentage was negative. Bacteroides genus correlated positively with valeric acid concentration and its relative percentage. IgA concentrations correlated negatively with isobutyric, isovaleric, and valeric acid concentrations. This indicator also correlated negatively with isobutyric and valeric absolute percentages. Correlation indexes are detailed in the Supplementary Table 5.
4. Discussion
Our study supports previously published results on the relationship between the aging process and gut health and gut microbiota in dogs. Notably, we included a wider age range (between 0.2 and 15 years old) than previous studies on puppies and adult dogs. In addition, we studied our dog population from different perspectives based on the analysis of fecal microbiota together with other intestinal health biomarkers to broadly understand the overall changes. Our main focus for the discussion of the results is the senior life stage, the period when the quality of life is threatened by the aging process itself. Exploring gut health and microbiota based on fecal indicators may help find specific diet interventions to improve the quality of life of aging dogs.
To compare the selected indicators along the different life stages, this study included dogs with specific diet needs because of their age, which may be considered a confounder to understanding the results among the age categories. However, the diet, together with body weight, sex, feeding routine, housing, and breed, were considered potential explanatory variables for the general statistical analysis.
4.1. Microbiota results
Bacillota, Actinomycetota, Fusobacteriota, Bacteroidota, and Pseudomonadota were the five main phyla found in the fecal samples analyzed in our studied population of dogs using 16S rRNA gene sequencing. Previous studies described that most bacterial sequences identified in the gut belong to these five different phyla (29–31). In this regard, previous literature may have used the old phylum terminology for Bacillota, Actinomycetota, Fusobacteriota, Bacteroidota, and Pseudomonadota, which were previously called Firmicutes, Actinobacteria, Fusobacteria, Bacterioidetes, and Proteobacteria, respectively (32). Although the names have recently changed, these phyla have been previously identified as the main contributors to gut microbial composition in dogs. More specifically, Bacillota, Bacteroidota, and Fusobacteriota were described as the three predominant phyla in the healthy canine fecal microbiome in previous studies (33, 34). However, in our study, within the three main bacterial groups, the phylum Actinomycetota was found instead of Bacteroidota. Interestingly, the distribution of relative abundance was the same across the three categories. The low abundance in Bacteroidota phyla could be because of obesity and gastrointestinal disease (35–37). Although Actinomycetota is present in the small intestine of healthy animals, the relative abundance found in feces is lower, and its increase may be associated with pathological conditions (29). Considering the health status of the population studied, obesity and gastrointestinal pathologies do not explain the low numbers obtained, and the differences may be explained by the analysis itself. The influence of methodological aspects in 16S rRNA gene sequencing on the results is important to consider, so different studies are difficult to compare. Previous studies have reported significant differences in the abundance of bacterial groups based on variations in labs and methodologies (38–40).
The Bifidobacteriaceae family is an important bacterial group for humans. Considering the results for Bifidobacterium in the qPCR and the sequencing analysis, this genus was present in all three age categories (senior dogs, adult dogs, and puppies), being more representative in adults and less representative in puppies. Although intestinal microbiota composition has been studied in different animal species, including dogs, the transition of the intestinal microbiota with age has not been thoroughly investigated. Masuoka et al. (10) studied a different age population, from pre-weanling puppies to senior dogs, to evaluate age-dependent differences and changes in the intestinal microbiota by employing a culture-based method. Bifidobacteria were detected in puppies but were not present when analyzing dogs older than 3 years old (adult dogs and senior dogs). Based on the literature on human beings, the presence of this genus in the three categories should be understood as positive. Bifidobacteria is the most prevalent bacteria in infants and adults, and it is presumed to promote health benefits in the host (41). However, the extrapolation of the results among animal species should be interpreted carefully. Bifidobacterium abundance in the qPCR correlated positively with valeric acid concentration and acetate relative percentage. Positive correlations between Bifidobacterium and acetate have been previously described (42, 43), although further studies investigating links with valeric acid may be warranted.
Senior dogs had a higher relative abundance of the families Ruminococcaceae, Veillonellaceae, Bacteroidaceae, and Lachnospiraceae (phylum Bacillota). Previous studies have found decreases at this level in dogs with inflammatory bowel disease (IBD) (25, 44), but there is no clear relationship with aging. Senior dogs also had higher values in the relative percentage of acetate, one of the main products of bacterial fermentation. These two findings are in the same direction as previously published data; there is a well-known positive correlation between the bacterial groups classified in Bacillota phyla and the production of SCFA (25). However, our results supported the fact that relative percentages of butyrate and valeric acid were significantly lower in the population of senior dogs. The relative percentages are logically influenced by the high numbers in the acetate concentration. Focusing on the valeric result, the relative percentage was supported by the concentration result, which was also lower in the senior group. These findings may be interpreted as negative considering gut health since SCFA production provides an appropriate pH environment to maintain healthy microbiota, inhibits the growth of pH-sensitive pathogenic bacteria, and preserves gut integrity (45). Moreover, valeric acid has been positively associated with the gut-brain axis and cognitive function in mice (46). Therefore, senior dogs may not benefit from the mechanisms surrounding the production of SCFA that promote intestinal health and positively contribute to cognitive function.
Although Faecalibacterium genera were not significantly different in the qPCR analysis, Faecalibacterium prausnitzii was found to be higher in older dogs, which may be considered a minor positive change for this aged population. This discrepancy may be caused by the fact that F. prausnitzii falls into the Faecalibacterium genera together with other species. Another possible reason could be that Facecalibacterium is frequently undetected in dogs during sequencing due to its low abundance, but it could be detected using qPCR. In addition, Faecalibacterium qPCR was significantly positively correlated with propionate and valeric acid concentrations, which is positive for gut health. The absence of this bacterium has been associated with gastrointestinal disorders in dogs (36) and humans (47). Focusing on the species level, F. prausnitzii has been clearly associated with increased fiber utilization and different fiber types in the diets of dogs (8, 48). However, the diet's fiber content was not considered a cause of the higher abundance described in the S category for F. prausintzii. In our study, only 5% of the dogs were out of the defined fiber composition range (2–5%), and they were similarly distributed in the three different age categories (J = 2/106; A = 2/106; and S = 1/106). From this information, considering the importance of F. prausnitzii within the Faecalibacterium genera, we could state that the difference found in senior dogs is a positive change in the microbiota population related to aging.
At the genus level, the two predominant groups identified belong to the phylum Bacillota (Uncl. Clostridiaceae.1 and Allobaculum). These findings support previously published data on dogs, in which major taxa fell within the Bacillota group and the bacterial class Clostridia was considered the most abundant taxon (29). In our study, the Fusobacterium genus (within the Fusobacteriota phylum) was present in all the age populations and was specifically higher in dogs over 7 years old. This finding can be considered positive for the microbiome of senior dogs since this genus has previously been related to a healthy intestinal microbiome in dogs (13).
4.2. Beta diversity, alpha diversity, and dysbiosis index
For beta diversity, unweighted UniFrac was significantly different among age categories when ANOSIM analysis was performed. However, the size effect as estimated by the R-value was extremely small, and other distance metrics were not significant. The changes found in our population, and especially the categories defined, may not be enough to change the weighted UniFrac and Bray–Curtis beta diversity. Another possible explanation for the observed differences in the unweighted UniFrac but not in the weighted UniFrac and Bray–Curtis beta diversity could be attributed to the fact that these two metrics consider relative abundances instead of only the presence/absence of different taxa. Furthermore, alpha diversity analysis measured by the Chao1 and Shannon indexes was not significantly different among the age groups. Our findings are consistent with previous studies. You and Kim (31) studied different individual traits in healthy dogs and also reported no differences at the age level when studying microbial diversity patterns. Interestingly, none of the dogs belonging to the junior category had a moderate or significant change in the index (DI > 0). Blake et al. (15) showed that the DI increased in younger puppies, but many dogs at 9 weeks and most dogs after 6 months of age had normal DI. In our study, the junior age ranged from 2 to 19 months; therefore, our data are consistent with previous studies, as most dogs over 6 months are expected to have a DI of <0. From those results, we could state that, in our study, small changes in the microbiota may have occurred due to aging. To conclude with the microbiota data, the performance of the Kruskal–Wallis test for assessing differential abundance between groups could be a limitation. This statistical test considers the overall abundance between groups but not the frequency of organism identification in samples within a group. Further statistical tests for assessing differential abundance between groups (e.g., ANCOM and DESeq) could complement the information analyzed by the Kruskal–Wallis test.
4.3. Fecal SCFAs, IgA, and cCP
Fecal cCP and IgA concentrations in feces have been studied as non-invasive biomarkers of gut health in adult dogs. Fecal cCP is an indicator of intestinal inflammation in dogs (49, 50), whereas IgA is associated more with the gut immune response (51). In the present study, the fecal concentration of cCP and IgA was not related to the aging process since no differences were observed among the three age categories. Based on the literature, aged mice and humans did not show a clear relationship between age and the concentration of mucosal IgA (52). Zaine et al. (53) reported no significant differences in IgA fecal values in senior dogs compared with adults and puppies. However, 5-month-old puppies had lower IgA concentrations than adult dogs. Fecal cCP concentrations were recently correlated with age in healthy humans in Korea (54), but to our knowledge, the results on cCP concentrations in healthy dogs of different ages are not consistent. The lack of clear reference values for fecal cCP and IgA considering the different age stages is a clear limitation to thoroughly understanding our results in relation to age and the SCFA correlation found in our study. However, in our study, age did not change the intestinal inflammation status or gut immune response based on the two indicators analyzed.
4.4. Main features
Considering the overall results of the present study, minor changes in the intestinal health biomarkers analyzed were found among the junior, adult, and senior categories. The microbiota community and SCFA production showed age-dependent differences. However, the IgA and cCP were not significantly different across the three categories.
Our age category definition may be a limitation to thoroughly understanding the continuous changes of the microbiome in relation to the aging process. Significantly, no clear consensus exists about this classification in animal research or veterinary medicine (55). The definition of age categories may reduce the sensitivity of understanding whether different ages equally affect gut health and microbiota. Thus, a larger number of animals representative of all age categories may help define more categories to increase the sensitivity of these first results.
In addition, the extrapolation of this trial in a field setting study may help better explore effective strategies, such as dietary interventions, in the dog population.
5. Conclusion
In conclusion, in this study, the canine gut microbiota and certain SCFAs showed minor variations among age groups in dogs. Other intestinal health biomarkers (IgA and cCP) were similar among age groups. Considering all the parameters, we can state that, once the microbiota becomes stable in healthy dogs, minor changes occur during the aging process.
The modulation of these minor changes, considering the intestinal microbiota and short-chain fatty acid production, could improve the overall gut health during aging in dogs.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA901473.
Ethics statement
The animal study was reviewed and approved by Affinity Petcare Ethical Committee.
Author contributions
AF-P, AS-M, and CT designed and performed the study. AF-P and RP analyzed the results and prepared the tables and graphs for inclusion in the manuscript. AF-P wrote the manuscript. Overall, all authors contributed to the paper, discussed its results, reviewed the drafts, and approved the final version submitted for publication.
Funding
This study received funding from Affinity Petcare. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Conflict of interest
AF-P, AS-M, and CT were employed by the pet food company Affinity Petcare, which funded the study. JS and RP were employed by the Gastrointestinal Laboratory at Texas A&M University, which provides microbiome assessment on a fee-for-service basis.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1213287/full#supplementary-material
Abbreviations
DI, Dysbiosis index; J, Junior; A, Adult; S, Senior; SCFA, Short Chain Fatty Acid; cCP, Calprotectin; IgA, Immunoglobulin A; PCoA, Principal Coordinate Analysis; DM, Dry Matter; GC-MS, Gas Chromatography-Mass Spectrometry; RIA, Radio Radioimmunoassay; ANOSIM, Analysis Of Similarity; R, Rotation; S, Stable; Emmeans, Estimated Means; FDR, False Discovery Rate; ANCOM, Analysis of Composition of Microbiomes; IBD, Inflammatory Bowel Disease.
References
1. Matijašić M, Meštrović T, Paljetak HC, Perić M, Barešić A, Verbanac D. Gut microbiota beyond bacteria—mycobiome, virome, archaeome, and eukaryotic parasites in IBD. Int J Mol Sci. (2020) 21:2668. doi: 10.3390/ijms21082668
2. Gill SR, Pop M, DeBoy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, et al. Metagenomic analysis of the human distal gut microbiome. Science. (2006) 312:1355–9. doi: 10.1126/science.1124234
3. Jandhyala SM, Talukdar R, Subramanyam C, Vuyyuru H, Sasikala M, Reddy DN. Role of the normal gut microbiota. World J Gastroentero. (2015) 21:8787–803. doi: 10.3748/wjg.v21.i29.8787
4. de Vos WM, Tilg H, Hul MV, Cani PD. Gut microbiome and health: mechanistic insights. Gut. (2022) 71:1020–32. doi: 10.1136/gutjnl-2021-326789
5. Bajinka O, Darboe A, Tan Y, Abdelhalim KA, Cham LB. Gut microbiota and the human gut physiological changes. Ann Microbiol. (2020) 70:65. doi: 10.1186/s13213-020-01608-2
6. Wen L, Duffy A. Factors Influencing the gut microbiota, inflammation, and type 2 diabetes. J Nutrition. (2017) 147:1468S–75S. doi: 10.3945/jn.116.240754
7. Hasan N, Yang H. Factors affecting the composition of the gut microbiota, and its modulation. PeerJ. (2019) 7:e7502. doi: 10.7717/peerj.7502
8. Pilla R, Suchodolski JS. The gut microbiome of dogs and cats, and the influence of diet. Vet Clin North Am Small Animal Pract. (2021) 51:605–21. doi: 10.1016/j.cvsm.2021.01.002
9. Badal VD, Vaccariello ED, Murray ER, Yu KE, Knight R, Jeste DV, et al. The gut microbiome, aging, and longevity: a systematic review. Nutrients. (2020) 12:3759. doi: 10.3390/nu12123759
10. Masuoka H, Shimada K, Kiyosue-Yasuda T, Kiyosue M, Oishi Y, Kimura S, et al. Transition of the intestinal microbiota of dogs with age. Biosci Microbiota Food Heal. (2016) 36:27–31. doi: 10.12938/bmfh.BMFH-2016-021
11. Suchodolski JS. Analysis of the gut microbiome in dogs and cats. Vet Clin Path. (2022) 50:6–17. doi: 10.1111/vcp.13031
12. Ziese A-L, Suchodolski JS. Impact of changes in gastrointestinal microbiota in canine and feline digestive diseases. Vet Clin North Am Small Animal Pract. (2021) 51:155–69. doi: 10.1016/j.cvsm.2020.09.004
13. Garrigues Q, Apper E, Chastant S, Mila H. Gut microbiota development in the growing dog: a dynamic process influenced by maternal, environmental and host factors. Front Vet Sci. (2022) 9:964649. doi: 10.3389/fvets.2022.964649
14. Guard BC, Mila H, Steiner JM, Mariani C, Suchodolski JS, Chastant-Maillard S. Characterization of the fecal microbiome during neonatal and early pediatric development in puppies. PLoS ONE. (2017) 12:e0175718. doi: 10.1371/journal.pone.0175718
15. Blake AB, Cigarroa A, Klein HL, Khattab MR, Keating T, Coevering PVD, et al. Developmental stages in microbiota, bile acids, and clostridial species in healthy puppies. J Vet Intern Med. (2020) 34:2345–56. doi: 10.1111/jvim.15928
16. FEDIAF. The European Pet Food Industry. Nutritional Guidelines for Complete and Complementary Pet Food for Cats and Dogs. (2021). Available online at: https://europeanpetfood.org/wp-content/uploads/2022/03/Updated-Nutritional-Guidelines.pdf
17. Strickling JA, Harmon DL, Dawson KA, Gross KL. Evaluation of oligosaccharide addition to dog diets: influences on nutrient digestion and microbial populations. Anim Feed Sci Technol. (2000) 86:205–19. doi: 10.1016/S0377-8401(00)00175-9
18. Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. (2015) 18:1403–14. doi: 10.1111/1462-2920.13023
19. Apprill A, McNally S, Parsons R, Weber L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol. (2015) 75:129–37. doi: 10.3354/ame01753
20. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. (2019) 37:852–7. doi: 10.1038/s41587-019-0209-9
21. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high resolution sample inference from Illumina amplicon data. Nat Methods. (2016) 13:581–3. doi: 10.1038/nmeth.3869
22. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microb. (2005) 71:8228–35. doi: 10.1128/AEM.71.12.8228-8235.2005
23. AlShawaqfeh MK, Wajid B, Minamoto Y, Markel M, Lidbury JA, Steiner JM, et al. A dysbiosis index to assess microbial changes in fecal samples of dogs with chronic inflammatory enteropathy. FEMS Microbiol Ecol. (2017) 93:fix136. doi: 10.1093/femsec/fix136
24. Moreau NM, Goupry SM, Antignac JP, Monteau FJ, Bizec BJL, Champ MM, et al. Simultaneous measurement of plasma concentrations and 13C-enrichment of short-chain fatty acids, lactic acid and ketone bodies by gas chromatography coupled to mass spectrometry. J Chromatogr B. (2003) 784:395–403. doi: 10.1016/S1570-0232(02)00827-9
25. Minamoto Y, Minamoto T, Isaiah A, Sattasathuchana P, Buono A, Rangachari VR, et al. Fecal short-chain fatty acid concentrations and dysbiosis in dogs with chronic enteropathy. J Vet Intern Med. (2019) 33:1608–18. doi: 10.1111/jvim.15520
26. Heilmann RM, Suchodolski JS, Steiner JM. Development and analytic validation of a radioimmunoassay for the quantification of canine calprotectin in serum and feces from dogs. Am J Vet Res. (2008) 69:845–53. doi: 10.2460/ajvr.69.7.845
27. Heilmann RM, Ruaux CG, Grützner N, Cranford SM, Bridges CS, Steiner JM. Biological variation of serum canine calprotectin concentrations as measured by ELISA in healthy dogs. Vet J. (2019) 247:61–4. doi: 10.1016/j.tvjl.2019.03.002
28. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
29. Pilla R, Suchodolski JS. The role of the canine gut microbiome and metabolome in health and gastrointestinal disease. Front Vet Sci. (2020) 6:498. doi: 10.3389/fvets.2019.00498
30. Deschamps C, Humbert D, Zentek J, Denis S, Priymenko N, Apper E, et al. From Chihuahua to Saint-Bernard: how did digestion and microbiota evolve with dog sizes. Int J Biol Sci. (2022) 18:5086–102. doi: 10.7150/ijbs.72770
31. You I, Kim MJ. Comparison of gut microbiota of 96 healthy dogs by individual traits: breed, age, and body condition score. Anim Open Access J. (2021) 11:2432. doi: 10.3390/ani11082432
32. Oren A, Garrity GM. Valid publication of the names of forty-two phyla of prokaryotes. Int J Syst Evol Microbiol. (2021) 71:005056. doi: 10.1099/ijsem.0.005056
33. Hand D, Wallis C, Colyer A, Penn CW. Pyrosequencing the canine faecal microbiota: breadth and depth of biodiversity. PLoS ONE. (2013) 8:e53115. doi: 10.1371/journal.pone.0053115
34. Middelbos IS, Boler BMV, Qu A, White BA, Swanson KS, Fahey GC. Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing. PLoS ONE. (2010) 5:e9768. doi: 10.1371/journal.pone.0009768
35. Salas-Mani A, Jeusette I, Castillo I, Manuelian CL, Lionnet C, Iraculis N, et al. Fecal microbiota composition changes after a BW loss diet in Beagle dogs. J Anim Sci. (2018) 96:3102–11. doi: 10.1093/jas/sky193
36. Minamoto Y, Otoni CC, Steelman SM, Büyükleblebici O, Steiner JM, Jergens AE, et al. Alteration of the fecal microbiota and serum metabolite profiles in dogs with idiopathic inflammatory bowel disease. Gut Microbes. (2015) 6:33–47. doi: 10.1080/19490976.2014.997612
37. Söder J, Wernersson S, Höglund K, Hagman R, Lindåse S, Dicksved J. Composition and short-term stability of gut microbiota in lean and spontaneously overweight healthy Labrador retriever dogs. Acta Vet Scand. (2022) 64:8. doi: 10.1186/s13028-022-00628-z
38. Handl S, Dowd SE, Garcia-Mazcorro JF, Steiner JM, Suchodolski JS. Massive parallel 16S rRNA gene pyrosequencing reveals highly diverse fecal bacterial and fungal communities in healthy dogs and cats. Fems Microbiol Ecol. (2011) 76:301–10. doi: 10.1111/j.1574-6941.2011.01058.x
39. Garcia-Mazcorro JF, Lanerie DJ, Dowd SE, Paddock CG, Grützner N, Steiner JM, et al. Effect of a multi-species synbiotic formulation on fecal bacterial microbiota of healthy cats and dogs as evaluated by pyrosequencing. Fems Microbiol Ecol. (2011) 78:542–54. doi: 10.1111/j.1574-6941.2011.01185.x
40. Guard BC, Barr JW, Reddivari L, Klemashevich C, Jayaraman A, Steiner JM, et al. Characterization of microbial dysbiosis and metabolomic changes in dogs with acute diarrhea. PLoS ONE. (2015) 10:e0127259. doi: 10.1371/journal.pone.0127259
41. Mitsuoka T. History and evolution of probiotics. Jpn J Lact Acid Bact. (2011) 22:26–37. doi: 10.4109/jslab.22.26
42. Oliphant K, Allen-Vercoe E. Macronutrient metabolism by the human gut microbiome: major fermentation by-products and their impact on host health. Microbiome. (2019) 7:91. doi: 10.1186/s40168-019-0704-8
43. Koh A, De Vadder F, Kovatcheva-Datchary P, Bäckhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. (2016) 165:1332–45. doi: 10.1016/j.cell.2016.05.041
44. Suchodolski JS, Dowd SE, Wilke V, Steiner JM, Jergens AE. 16S rRNA gene pyrosequencing reveals bacterial dysbiosis in the duodenum of dogs with idiopathic inflammatory bowel disease. PLoS ONE. (2012) 7:e39333. doi: 10.1371/journal.pone.0039333
45. Mackei M, Talabér R, Müller L, Sterczer Á, Fébel H, Neogrády Z, et al. Altered intestinal production of volatile fatty acids in dogs triggered by lactulose and psyllium treatment. Vet Sci. (2022) 9:206. doi: 10.3390/vetsci9050206
46. Lai Z, Shan W, Li J, Min J, Zeng X, Zuo Z. Appropriate exercise level attenuates gut dysbiosis and valeric acid increase to improve neuroplasticity and cognitive function after surgery in mice. Mol Psychiatr. (2021) 26:7167–87. doi: 10.1038/s41380-021-01291-y
47. Varela E, Manichanh C, Gallart M, Torrejón A, Borruel N, Casellas F, et al. Colonisation by faecalibacterium prausnitzii and maintenance of clinical remission in patients with ulcerative colitis. Aliment Pharm Therap. (2013) 38:151–61. doi: 10.1111/apt.12365
48. Mori A, Goto A, Kibe R, Oda H, Kataoka Y, Sako T. Comparison of the effects of four commercially available prescription diet regimens on the fecal microbiome in healthy dogs. J Vet Med Sci. (2019) 81:1783–90. doi: 10.1292/jvms.19-0055
49. Grellet A, Heilmann RM, Lecoindre P, Feugier A, Day MJ, Peeters D, et al. Fecal calprotectin concentrations in adult dogs with chronic diarrhea. Am J Vet Res. (2013) 74:706–11. doi: 10.2460/ajvr.74.5.706
50. Otoni CC, Heilmann RM, García-Sancho M, Sainz A, Ackermann MR, Suchodolski JS, et al. Serologic and fecal markers to predict response to induction therapy in dogs with idiopathic inflammatory bowel disease. J Vet Intern Med. (2018) 32:999–1008. doi: 10.1111/jvim.15123
51. Grellet A, Mila H, Heilmann RM, Feugier A, Gruetzner N, Suchodolski JS, et al. Effect of age, gestation and lactation on faecal IgA and calprotectin concentrations in dogs*. J Nutritional Sci. (2014) 3:e41. doi: 10.1017/jns.2014.44
52. Fujihashi K, McGhee JR. Mucosal immunity and tolerance in the elderly. Mech Ageing Dev. (2004) 125:889–98. doi: 10.1016/j.mad.2004.05.009
53. Zaine L, Ferreira C, de OS Gomes M, Monti M, Tortola L, Vasconcellos RS, et al. Faecal IgA concentration is influenced by age in dogs. Brit J Nutr. (2011) 106:S183–6. doi: 10.1017/S0007114511000559
54. Park SY. Age-related fecal calprotectin concentrations in healthy adults. Korean J Clin Laboratory Sci. (2020) 52:181–7. doi: 10.15324/kjcls.2020.52.3.181
Keywords: fecal microbiota, canine, aging, intestinal health, nutrition
Citation: Fernández-Pinteño A, Pilla R, Manteca X, Suchodolski J, Torre C and Salas-Mani A (2023) Age-associated changes in intestinal health biomarkers in dogs. Front. Vet. Sci. 10:1213287. doi: 10.3389/fvets.2023.1213287
Received: 27 April 2023; Accepted: 13 July 2023;
Published: 22 August 2023.
Edited by:
Shengfa F. Liao, Mississippi State University, United StatesReviewed by:
Zhongyue Yang, University of California, Davis, United StatesEmily Coffey, University of Minnesota Twin Cities, United States
Copyright © 2023 Fernández-Pinteño, Pilla, Manteca, Suchodolski, Torre and Salas-Mani. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Anna Fernández-Pinteño, YW5mZXJuYW5kZXomI3gwMDA0MDthZmZpbml0eS1wZXRjYXJlLmNvbQ==