Next Article in Journal
The Defensive Behaviors and Milk Production of Pastured Dairy Cattle in Response to Stable Flies, Horn Flies, and Face Flies
Previous Article in Journal
Pharmacokinetics, Tissue Residues, and Withdrawal Times of Oxytetracycline in Rainbow Trout (Oncorhynchus mykiss) after Single- and Multiple-Dose Oral Administration
Previous Article in Special Issue
A Comparison of the Oral Microbiota in Healthy Dogs and Dogs with Oral Tumors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 13
Issue 24
10.3390/ani13243846
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Longitudinal Analysis of Canine Oral Microbiome Using Whole Genome Sequencing in Aging Companion Dogs

by
Ginger B. Templeton
1,
Gilad Fefer
1,
Beth C. Case
1,
Jeff Roach
2,
M. Andrea Azcarate-Peril
2,
Margaret E. Gruen
1,
Benjamin J. Callahan
3,4 and
Natasha J. Olby
1,*
1
Department of Clinical Sciences, College of Veterinary Medicine, North Carolina State University, Raleigh, NC 27607, USA
2
Department of Medicine, Division of Gastroenterology and Hepatology, and UNC Microbiome Core, Center for Gastrointestinal Biology and Disease, School of Medicine, University of North Carolina, Chapel Hill, NC 27599, USA
3
Department of Population Health and Pathobiology, North Carolina State University, Raleigh, NC 27607, USA
4
Bioinformatics Research Center, North Carolina State University, Raleigh, NC 27695, USA
*
Author to whom correspondence should be addressed.
Animals 2023, 13(24), 3846; https://doi.org/10.3390/ani13243846
Submission received: 29 October 2023 / Revised: 10 December 2023 / Accepted: 11 December 2023 / Published: 14 December 2023
(This article belongs to the Special Issue Investigations of the Oral Microbiota in Health and Disease in Companion Animals)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The prevalence of dental disease and cognitive decline in elderly dogs is extremely high, and, given the known relationship between dental disease and Alzheimer’s Disease in people, this study sought to describe the changes in oral microbiota in aged pet dogs over time. By sequencing oral swabs, we were able to identify bacterial and fungal populations in the dogs’ mouths. The most common bacterial species present, Phorphorymonas spp. is known to produce factors that cause neurodegeneration. Moreover, Leptotrichia, another bacterial species present, correlated to cognition scores in these dogs. We conclude that this small exploratory study shows the importance of defining the oral microbiota in aged dogs with a view to understanding potential therapeutic targets. Larger prospective studies should be undertaken as a priority.

Abstract

Aged companion dogs have a high prevalence of periodontal disease and canine cognitive dysfunction syndrome (CCDS) and the two disorders are correlated. Similarly, periodontal disease and Alzheimer’s Disease are correlated in people. However, little is known about the oral microbiota of aging dogs. The goal of this project was to characterize the longitudinal changes in oral microbiota in aged dogs. Oral swabs were taken from ten senior client-owned dogs on 2–3 occasions spanning 24 months and they underwent whole genome shotgun (WGS) sequencing. Cognitive status was established at each sampling time. A statistically significant increase in alpha diversity for bacterial and fungal species was observed between the first and last study visits. Bacteroidetes and proteobacteria were the most abundant bacterial phyla. Porphyromonas gulae was the most abundant bacterial species (11.6% of total reads). The species Lactobacillus gasseri had a statistically significant increase in relative abundance with age whereas Leptotrichia sp. oral taxon 212 had a statistically significant positive longitudinal association with cognition score. There is an increased fungal and bacterial alpha diversity in aging dogs over time and nearly universal oral dysbiosis. The role of the oral microbiota, particularly Leptotrichia and P. gulae and P. gingivalis, in aging and CCDS warrants further investigation.

1. Introduction

The oral microbiota is a complex ecosystem in both humans and dogs [1,2]. Oral-associated bacteria have been implicated in diseases of inflammation and senescence, including cardiovascular disease, diabetes mellitus, cancers, and Alzheimer’s disease (AD) [3,4,5,6,7]. Companion dogs make an excellent model system for healthy and unhealthy human aging as they share the human environment, have diseases in common with humans, and yet, age at an accelerated rate [8]. Defining the oral microbiota of senior and geriatric companion dogs is critical to our understanding of healthy canine aging and senescence, which, in turn, can be applied to a greater understanding of human aging.
The field of canine microbiome research is in its early stages but it is advancing rapidly. A number of studies have focused on the gastrointestinal bacterial microbiome, with fewer investigating the oral and nasal environments. Early studies of oral bacteria were limited by their reliance on culture-based methods. In contrast, recent studies have taken advantage of next-generation sequencing, primarily 16S rRNA amplicon sequencing analysis, to reveal a much greater diversity of organisms in the canine oral cavity. These studies have primarily focused on populations of young to middle-aged dogs and, in many cases, have been limited to purpose-bred, working, or kenneled dogs of single or similar breeds [9,10,11,12]. One recent study reported a diverse fungal population in the canine mouth using internal transcribed spacer analysis [13]. Still, research into the canine oral mycobiome is even less advanced than the bacterial microbiome.
In humans, the progression of periodontal disease correlates with AD progression [14]. Chronic inflammation has been suggested as an indirect cause of this relationship [15], whereas the keystone periodontal pathogen Porphyromonas gingivalis has been suggested as a direct cause, mediated by virulence factors called gingipains [16].
Canine cognitive dysfunction syndrome (CCDS) bears remarkable similarity at both the cellular and behavioral levels to AD. Dogs with CCDS develop pathognomonic central nervous system features of AD, including cerebral atrophy, ventricular widening, deposition of plaques of misfolded Aβ in the prefrontal cortex and cerebral vasculature as well as neuronal loss affecting the cerebral cortex, hippocampus, and limbic system [17]. CCDS prevalence has been estimated as high as 67% for dogs aged 16–17 years based on owner-reported symptoms in the domains of disorientation, changes in social interactions, house training and the sleep/wake cycle [18]. Further, findings in our laboratory show a significant decline in performance on laboratory tests of cognition in geriatric companion dogs [19]. This decline is accompanied by escalating serum neurofilament light chain concentrations, and elevations followed by reductions in serum concentrations of amyloid beta 42, mirroring changes seen in patients with Alzheimer’s Disease [18]. Based on these studies of histopathology, MRI findings, blood biomarker alterations, behavioral changes and reduced performance on cognitive tests, companion dogs can serve as an excellent real-world model for AD.
The primary objective of this exploratory study was to define the oral microbiota of senior dogs using whole genome sequencing to enhance species-level accuracy over previous studies relying on 16S rRNA amplicon sequencing. A secondary study aim was to evaluate the correlations between the oral microbiome and naturally occurring age-related disease processes, including CCDS.

2. Materials and Methods

2.1. Animals

Dogs in this study were companion animals, and owners signed informed consent prior to enrollment in the study. All procedures were reviewed and approved by the North Carolina State University Institutional Animal Care and Use Committee (protocol number 18-109-O) and all work performed was in accordance with the guidelines and regulations. All methods are reported in accordance with ARRIVE guidelines for the reporting of animal experiments.
Senior companion dogs (small dogs > 10 years, large dogs > 8 years, giant dogs > 6 years of age) were recruited for a longitudinal study of neuro-aging. To be included, dogs were required to have intact hearing, intact vision, and the ability to walk independently at the time of enrollment. Dogs with urinary tract infections and those receiving medications that could impair cognition (e.g., antiepileptic drugs such as phenobarbital, potassium bromide, levetiracetam, zonisamide or behavior modifying drugs) were excluded. Subjects were evaluated at approximately 6-month intervals, each evaluation spanning three non-consecutive hospital day-visits. This study spanned the emergence of COVID-19 and, as a result, some assessment dates were delayed or altogether canceled.
To account for differences in life expectancy due to size and breed, lifespan was predicted with the following formula: L = 13.620 + (0.0702H) − (0.0538W) as described [20] where H = height (cm) and W = weight (lbs.). Fraction of lifespan (FoL) was then calculated: FoL = L/age. Both absolute age and FoL were included as independent covariates in statistical analyses.
Table 1 provides information on the dogs recruited into the study.

2.2. Covariate Data Collection

At each 6-month evaluation, diet, medications, and supplements being given were recorded and no restrictions were placed. Owners completed the canine dementia scale (CADES) and the canine brief pain inventory (CBPI) questionnaires. The CADES score is a validated metric of neurobehavioral abnormalities exhibited by dogs developing CCDS. The CADES score ranges from 0–95, and subjects are classified as normal (<8), mild (8–23), moderate (24–44), and severe (45–95) [21]. For this study, subjects were divided into two cognitive groups based on CADES scores: normal to mild impairment (NM, 0–23) and moderate to severe (MdSv, 24–95). The CBPI was designed to measure pain and has been validated in the assessment of pain caused by osteoarthritis (OA) [22] and bone cancer [23] and response to pain medication in OA [24].
On day one of each six-month evaluation, one of two veterinarians performed a physical exam, including body and muscle condition score, orthopedic and neurological examination, and indices of gingivitis (GI), plaque (PI), and calculus (CI) as previously described [25]. Oral examinations were not performed under sedation or anesthesia and so detailed dental assessments were not possible.
Oral swabs were collected on the third day of evaluation. A total of 27 samples from 10 dogs were collected on 2–3 dates spanning 23 months. Swabs were performed by one of two veterinarians wearing sterile gloves and using a sterile oral swab. The supragingival and adjacent buccal mucosa were swabbed for 30–60 s using the same swab. No sedation was used when collecting the samples. Samples were immediately placed in sterile tubes on ice and moved to −80 °C within 30 min of collection.
Behavioral and cognitive testing was performed on days 2 and 3 of the evaluation and included inhibitory control (cylinder task, IC), reversal learning (cylinder detour), and the sustained gaze test (SGT) as previously described [19,26].

2.3. DNA Isolation and Sequencing

Oral swabs were submitted to the University of North Carolina (UNC) Microbiome Core and treated as follows: DNA isolation was performed using an optimized version of the QIAamp Fast DNA Stool Mini Kit (Cat No./ID: 51604) protocol supplemented with 60 mg/mL lysozyme (Thermo Fisher Scientific, Grand Island, NY, USA) [27,28,29]. DNA was quantified via Quant-iTTM PicoGreenTM dsDNA quantification reagent.
For the validation of the DNA isolation process, a known bacterial community, ZymoBIOMICS Microbial Community Standard (Cat# D6300), and a blank composed of only DNA isolation reagents was included in the DNA extraction process and again in the library preparation. In addition to the isolation controls, the library preparation also included a library blank composed of library preparation reagents alone.
Swift Whole Genome Shotgun (WGS) sequencing was then performed. For samples with a DNA concentration <2 ng/uL, 5 ng of genomic DNA and for samples with DNA concentration >2 ng/uL, 25 ng of genomic DNA was processed using the Swift 2 S Turbo DNA Library Kit (Swift Biosciences, Ann Arbor, MI, USA; Lot#04030036P). Targeted DNA was treated with the Swift 2 S Turbo Enzyme Mix to induce dsDNA fragmentation, repair, and A-tailing. The thermal profile for fragmentation was 32 °C for 11 and 9 min, respectively, based on input concentration, followed by 60 °C for 30 min. Next, the fragmented DNA was ligated using the Swift 2 S Turbo kit ligation mix, followed by a 20 min incubation at 20 °C. Next, index 1(D7) and index 2(D5) were added along with the Swift 2 S Turbo indexing mix. The thermal profile for the amplification had an initial extension denaturation step at 98 °C for 30 s, followed by 8 and 5 cycles of denaturation, respectively, based on input concentration at 98 °C for 10 s each, annealing at 60 °C for 30 s, and a 60 s extension at 68 °C. The DNA library was then purified using Agencourt® AMPure® XP Reagent (Beckman Coulter, Indianapolis, IN, USA). Each sample was quantified and normalized prior to pooling. The DNA library pool was loaded on the Illumina platform reagent cartridge (Illumina, San Diego, CA, USA) and the Illumina instrument [30].
The prepared pool sequencing suitability was initially verified on the Illumina Nano platform. Following successful preliminary sequencing, prepared pools were submitted to the Illumina NovaSeq SP-XP platform with run configuration PE/150x v1.5. Sequencing output from the Illumina NovaSeq platform was converted to the fastq format and demultiplexed using Illumina Bcl2Fastq 2.20.0 [31]. Quality control of the demultiplexed sequencing reads was verified by FastQC 0.11.9 [32]. Adapters were trimmed using Trim Galore 0.6.7 [33]. The resulting paired-end reads were classified with Kraken 2.1.2 [34], and Bracken 2.5 [35] with the addition of the canine and human reference genomes and canine-relevant bacterial species (Supplementary Table S1), and all reads identified as canine or human were eliminated. Bracken output and covariates were loaded into R 4.1.3 [36]. tidyverse 1.3.1 [37] was used for data organization and phyloseq 1.36.0 [38] and vegan 2.6-2 [39] were used for calculation of alpha and beta diversity. Visualization was performed with ggplot2 3.3.5 [40] and JMP Pro 16 [41].

2.4. Statistical Analysis

Statistical analysis was performed in JMP Pro 16 [41]. Bacterial and fungal communities were analyzed separately from one another. A change in the Shannon Index (H) or Simpson’s Index of Diversity (1-D) within individuals was analyzed using the Wilcoxon Signed-Rank test for first and last samples of bacterial or fungal. The relationship between bacterial or fungal alpha diversity (H or 1-D) and a number of covariates were analyzed using a linear mixed effects model. Specifically, CADES score or CBPI were treated as the dependent variables, subject as random effect, and FoL or age and alpha diversity metrics as fixed effects, Separately, H or 1-D were treated as dependent variables, subject as random effect, and FoL, age, CADES group (NM, MdSv) and dental prophylaxis (dental cleaning) in the preceding six months or FoL, inhibitory control, and reversal learning as fixed effects.
The abundance of individual taxa was longitudinally analyzed using the negative binomial mixed model (NBMM) [42]. Taxa present in fewer than 20% of samples or assigned fewer than 2500 reads across all samples were removed. In total, 3 models were considered with fixed effects: FoL alone, CADES alone, and FoL and CADES concurrently. In all 3 models sequencing pool was included as a fixed effect and subject was included as a random effect. Analogously, 2 additional models were considered replacing CADES with SGT. Adjusted p-values were reported as calculated by NBMM.

3. Results

3.1. Clinical Findings

The first 10 dogs enrolled in a larger longitudinal study of canine neuro-aging to provide 2 to 3 consecutive oral (supragingival/buccal mucosal swabs) samples over 24 months were selected for exploratory analysis of the aging canine oral microbiome. At that time, 24 dogs were enrolled in the longitudinal study but many had missed rechecks due to the COVID pandemic. Physical characteristics at the time of enrollment are summarized (Table 1). Subjects were generally in good health throughout the study, though common age-related comorbidities developed and are summarized in Supplementary Table S2 along with owner-reported history of medications, probiotics, and diet.
All owners denied routine toothbrushing or did not respond to the survey question, which was interpreted as not routinely brushing their teeth. Clinical scores of gingivitis, plaque and calculus and history of recent dental prophylaxis (within four months of sample collection) are summarized (Table 2).
One subject experienced severe cognitive decline over the study, as evidenced by a 45-point CADES score increase. The remaining nine subjects were stable within their CADES group throughout the study, seven were classified as NM, and two were classified as MdSv, Figure 1. CADES and CBPI scores were not obtained at the third (last) study visit for subject G04. CADES was in the MdSv group for this subject at the first and second visits and thus was assumed to be MdSv for the third and treated as the last observation carried forward for both CADES score and CADES group.

3.2. Microbiome Findings

The datasets generated and analysed during the current study are available in the NCBI Sequence Read Archive repository, https://www.ncbi.nlm.nih.gov/sra/PRJNA950991 (accessed on 28 April 2023) Submission ID: PRJNA950991.
Only 1 sample (Subject G02, Sample 1) was discarded due to low reads (34,243). The remaining 26 samples ranged from seven to 16 million (7,056,183–16,211,066) raw reads.
BRACKEN (Bayesian Reestimation of Abundance with KrakEN) estimates of canine host content ranged by sample from 22.2–99.6%, with a median 92.0% and mean 84.2%. Estimated percent human sequence was less than 0.2% in all samples. Of the remaining taxa, bacterial relative abundance was 99.4%, eukaryotic (fungal) relative abundance was 0.5%, and less than 0.001% was archaeal or viral. Due to low abundance, viral and archaeal taxa were removed prior to subsequent analysis. Fungal and bacterial taxa were analyzed and reported as two separate populations.
Bacterial species counts from Kracken2/Bracken analysis are reported by sample (Supplementary Table S3). 5378 bacterial species were identified from 35 phyla. Of these, 10 phyla accounted for greater than 99.7% of relative abundance (Table 3), with the majority of species classified as Bacteroidetes (52.1%) or Proteobacteria (35.3%) (Figure 2).
Actinobacteria (3.7%) and Firmicutes (3.1%) were the third and fourth-most abundant phyla, respectively. These findings were similar to those in a study of young working dogs [12], but in contrast to much of the existing canine oral microbiome literature in young to middle-aged dogs, where the Firmicutes are predominant [9,11] or in much higher abundance [10]. It should be noted that previous reports used 16S rRNA amplicon sequencing analysis and these different methodologies could yield very different results, making direct comparison challenging. Within individual samples, Bacteroidetes or Proteobacteria were the predominant phyla in 25 of 26 samples (Figure 3), and Tenericutes was predominant in the remaining sample.
Shannon Index (H) range for bacterial taxa from oral samples was 1.97–5.57, with a median of 4.15 and mean of 4.18. The Simpson’s Index of Diversity (1-D) range for bacterial taxa was 0.55–0.98, with a median of 0.94 and a mean of 0.90. A statistically significant increase in bacterial alpha diversity was observed between the first and last samples (Wilcoxon Signed-Rank, p = 0.0039 for Shannon Index and p = 0.002 for Simpson’s Index of Diversity) (Figure 4).
In total, 546 bacterial species, each consisting of greater than 0.01% of total counts, represented 89.7% of the total species estimated. The canine oral pathogen, Porphyromonas gulae, was the most prevalent (11.6%), with the closely related human oral pathogen, P. gingivalis at 3.1%. These two organisms are considered keystone species in periodontal disease in their respective hosts [43,44] and will be collectively referred to as Gingipains-producing Porphyromonas species (GPPs). The second-most abundant species was P. cangingivalis (10.5%), which has been reported as the predominant health-associated oral species in adult dogs and is also common in early periodontal disease [45]. In addition to GPPs, the two other “Red Complex” species, Tannerella forsythia (3.4%) and Treponema denticola (1.6%), known for their role in periodontal disease [14], were among the 20 most abundant species in the study (Table 4).

3.3. Relationship between Microbiome, Age, Cognition and Pain

Using a linear mixed effects model, with the subject as a random effect, there was no statistical correlation between bacterial alpha diversity (H or 1-D) and age, FoL, CADES, CADES group, CBPI, SGT (for H), inhibitory control, reversal learning, or recent dental prophylaxis (H and 1-D) (see Supplementary Table S4). There was no significant correlation of bacterial alpha diversity (H or 1-D) with the combined relative abundance of the two known producers of the virulence factor gingipains, P. gulae and P. gingivalis (heretofore collectively referred to as GPPs) or the combined relative abundance of the so-called “red complex bacteria”, a defined group of human periodontal pathogens consisting of P. gingivalis, T. forsythia, and T. denticola, and including P. gulae in this report [15]. There were no apparent trends in beta diversity (Bray–Curtis) with respect to CADES, age, or FoL (see Supplementary Table S5 for the Bray–Curtis matrix). However, the lack of significant relationships could simply represent the small sample size.
Longitudinal analysis of individual taxa using negative binomial mixed models including FoL and sequencing pool as fixed effects and subject as random effect showed the species Lactobacillus gasseri to have a statistically significant (adjusted p = 0.0006) increase in relative abundance with FoL. Models with fixed effects CADES and sequencing pool with subject as a random effect indicated a modest, yet statistically significant (adj. p-value 2.9289 × 10−5) increase in relative abundance of Leptotrichia sp. oral taxon 212 with increased CADES score. Models including both FoL and CADES as well as sequencing pool as fixed effects and subject as a random effect again identified increased Lactobacillus gasseri as significantly associated with increased FoL. However, for Leptotrichia sp. oral taxon 212 neither FoL nor CADES were statistically significant likely due to the correlation between FoL and CADES. Model coefficients and p-values for all models are reported in Table 5. No species were identified to be longitudinally associated with SGT in any of the tested models.

3.4. Fungal Species

The Shannon Index (H) range for fungal taxa was 0.99–3.74, with a median of 1.47 and a mean of 1.71. The Simpson’s Index of Diversity (1-D) range was 0.36–0.86, with a median of 0.49, and a mean of 0.53. A statistically significant increase in fungal alpha diversity was observed between first and last samples (Wilcoxon Signed-Rank, p = 0.07 for Shannon Index and p = 0.01 for Simpson’s Index of Diversity) see Supplementary Table S2. Using a linear mixed effects model, with subject as random effect, there were no statistically significant correlations between fungal alpha diversity (H or 1-D) and age, FoL, dental prophylaxis in prior six months, CADES score, CBPI, inhibitory control, reversal learning or SGT.
In total, 66 fungal species, each of which consisted of at least 0.01% of total fungal counts, representing 100% of total fungal species estimated, were identified. The most abundant species in total was Aspergillus oryzae (67%), which was present in all 26 samples, and the most abundant species in each sample. A. oryzae is considered a non-harmful Aspergillus species used in food production, and, to our knowledge, has not been previously reported in studies of the canine oral mycobiome [46]. The opportunistic pathogen Candida dubliensis, was present in all samples, while the closely related pathogenic yeast, C. albicans, was present in 22 of 26 samples [47]. Eremothecium sinecaudum, a plant pathogen that has been identified as a component of the healthy human lung mycobiome, was present in 24 of 26 of samples [48]. Ten additional fungal species were present in all 26 samples (Table 6), the majority of which have been previously described as plant pathogens [49].

4. Discussion

To our knowledge this study represents the first use of WGS sequencing to define the oral microbiota of aging companion dogs, and the first to investigate the oral microbiota in this population longitudinally. We identified a bacterial community consistent with oral dysbiosis in the majority of our subjects over the span of the study. Two closely-related organisms, P. gulae and P. gingivalis, were predominant. These bacteria are considered keystone species in periodontal disease in dogs [43] and humans [14], respectively. Moreover, P. gingivalis has been implicated in AD in humans, either indirectly through chronic inflammation or directly through virulence factors including small proteases called gingipains [44]. P. gulae produces gingipains and other virulence factors similar to P. gingivalis [50]. A recent study showed a statistically significant correlation between periodontal disease and cognitive decline in dogs [51]. Additionally, in a small study of aged beagles, P. gulae DNA and gingipains antigen was isolated from brain tissue, though cognition was not measured, and a gingipains inhibitor reduced clinical signs associated with periodontal disease in that study [52].
Longitudinal analysis of individual taxa allows changes from baseline to be interrogated with respect to age and cognitive state regardless of initial state. This aspect is essential given the potentially vast variation in canine oral microbiome at later stages of life. We identified one organism with a statistically significant positive correlation with CADES score, Leptotrichia sp. oral taxon 212, and one species with a statistically significant positive correlation with FoL, Lactobacillus gasseri. Leptotrichia species are facultative anaerobic Gram-negative bacilli, and are common inhabitants of the human oral microbiota [3] that have been associated with both oral health and disease states in humans [53] and in gingivitis in dogs [11]. Leptotrichia species have been associated with AD [54,55] and periodontal disease [53]. In humans. Leptotrichia wadei, in particular, has been observed to be significantly enriched in subjects with mild cognitive impairment [54]. However, species of the genus Leptotrichia have been associated with both periodontal health and disease suggesting that individual Leptotrichia species have “distinct pathogenic potential” [53]. The relationship between CCDS and Leptotrichia, P. gulae and P. gingivalis merits further investigation.
Several Lactobacillus species, Lactobacillus gasseri in particular, have been previously identified in dog faeces [56,57] and human breast milk derived Lactobacillus gasseri has been studied as a probiotic treatment for obesity in dogs [58]. The study of association of Lactobacillus gasseri with aging or cognition outside of model organisms [59], however, remains largely limited to the study of probiotic supplementation [60].
The biological significance of the fungal species identified in this study is unclear, but it is interesting that 12 species were found in all 26 samples, including A. oryzaeh, the most abundant species present, and the opportunistic pathogen, C. dubliniensis. This is in contrast to published findings of the canine oral mycobiome in a younger study population where none of the observed species occurred in all samples [13]. We report a statistically significant increase in both bacterial and fungal alpha diversity over the course of the study. In humans, oral alpha diversity has been shown to decrease with age [61]. Our results may represent a difference between the oral ecosystem of aged humans and senior companion dogs, differences in behaviors between humans and dogs, the onset of age-related behaviors (such as an increase or decrease in self-grooming) or study-related factors such as sample size or duration of this study.
This study was an exploratory evaluation of the oral microbiota and the small cohort size limits the conclusions that can be drawn. The fact that only one dog progressed into the MdSv CADES group during the course of the study limited the study’s power to examine the relationship between cognitive change and changes in the oral microbiota. With a larger sample size and a longer timeline, we anticipate a larger number of dogs progressing through CADES groups. This would allow examination of differences in the oral microbiome pre-CCDS and after the diagnosis of CCDS within individuals, potentially elucidating a role for the oral microbiome in CCDS development or progression. Moreover, cognitive aging occurs at different rates in different dogs and comparison of the oral microbiota between successful agers, those that develop mild impairment and those that develop severe impairment could also be evaluated [62].
Because these data were generated as part of a broader study of neuro-aging intended to investigate natural aging in companion dogs, subjects with antibiotic use and dental prophylaxis were not excluded. As our study size and duration continue to increase, so too will the power, potentially enabling further understanding of the relationship between dental care, medications, oral microbiome, and cognition. At present, we found no correlation between recent dental prophylaxis and either alpha diversity, relative abundance of GPPs, or relative abundance of red complex bacteria. However, it is important to note that comprehensive oral examinations under anesthesia were not performed in these dogs. In this longitudinal study of aging dogs’ health, the central tenet is to cause no harm to the dogs during data collection. As such, anesthesia and sedation are avoided at the routine 6-monthly assessments. The use of the calculus, gingival and plaque indices allowed the description of some aspects of periodontal disease but should not be considered a complete evaluation of dental and oral health.
The high prevalence of cognitive decline in aging dogs coupled with the abundance of P. gulae and P. gingivalis in the aging cohort in the present study, leads us to postulate that these species may play a role in CCDS. While it is well-established that periodontal disease is very common in aging dogs [63], the extent to which oral dysbiosis contributes to senescence in dogs has not been fully elucidated. Defining causative organisms in CCDS and elucidating the mechanisms involved will provide potential preventative and therapeutic targets for CCDS.
A limitation and advantage of this study is that the intraoral sites sampled (supragingival and buccal mucosa) is a combination of two distinct niches. This limits our ability to compare findings to prior studies. However, we chose this method for its ability to represent a greater portion of the oral cavity than single-site sampling and thus is likely more representative of the overall health state of the mouth. Moreover, establishing the aging oral microbiome using a method that does not require general anesthesia provides a baseline for future studies of the aging canine oral microbiome that may wish to avoid the risk or expense of anesthesia. Understanding this site is also valuable considering the availability of commercially-available tests of the oral microbiome which rely on owner-administered oral swabs for sample collection.

5. Conclusions

Dysbiosis of oral microbiota is the norm in aging dogs, in spite of dental prophylaxis. The high prevalence of red complex bacterial species points to a potential contributing factor in the development of cognitive decline and we conclude that evaluation of this relationship should take priority, given the potential to intervene therapeutically.
As our knowledge of the aging canine oral microbiota evolves, so too will our understanding of healthy aging and age-related disease. The oral microbiota is a potential diagnostic tool and therapeutic target. Our study is unique in its longitudinal design, use of companion, rather than lab-housed or purpose bred dogs, and, most importantly, the focus on senior dogs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13243846/s1, Table S1: Species added to the Bracken database based on their presence and abundance in Humann2 output and potential relevance to the canine microbiome. Table S2: Demographic, oral health, cognitive and Shannon and Simpson indices for each individual study participant. Table S3: Bacterial species counts for each sample. Table S4: Results of linear mixed models examining the relationship of outcome variables with bacterial diversity (H and 1-D) with subject as the random effect. Table S5: Bray–Curtis Matrix for prokaryotic species.

Author Contributions

Conceptualization: M.A.A.-P., B.J.C., M.E.G. and N.J.O.; Methodology: M.A.A.-P., J.R. and B.J.C.; software: J.R. and G.B.T.; validation: J.R. and G.B.T.; formal analysis: G.B.T., J.R. and N.J.O.; investigation: G.B.T., G.F., B.C.C., M.E.G. and N.J.O.; resources: M.A.A.-P., M.E.G. and N.J.O.; data curation: B.C.C., G.B.T. and G.F.; writing—original draft preparation: G.B.T. and N.J.O.; writing—review and editing: G.B.T., G.F., B.C.C., J.R., M.A.A.-P., B.J.C., M.E.G. and N.J.O.; visualization, G.B.T.; supervision, M.E.G. and N.J.O.; project administration, N.J.O.; funding acquisition, N.J.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding but was funded internally by the Dr. Kady M. Gjessing and Rhanna M. Davidson Distinguished Chair in Gerontology.

Institutional Review Board Statement

All procedures were reviewed and approved by the North Carolina State University Institutional Animal Care and Use Committee (protocol number 18-109-O) and all work performed was in accordance with the guidelines and regulations. “The animal study protocol was approved by the North Carolina State University Institutional Animal Care and Use Committee (protocol number 18-109-O, approved 6 August 2018).

Informed Consent Statement

Informed consent was obtained from all dog owners involved in the study.

Data Availability Statement

The datasets generated and analyzed during the current study are openly available in the NCBI Sequence Read Archive repository, https://www.ncbi.nlm.nih.gov/sra/PRJNA950991 (accessed on 28 April 2023) Submission ID: PRJNA950991. All data used for individual analyses are provided in Supplementary Tables S1–S5.

Acknowledgments

The authors would like to acknowledge the dogs and their owners.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caselli, E.; Fabbri, C.; D’Accolti, M.; Soffritti, I.; Bassi, C.; Mazzacane, S.; Franchi, M. Defining the oral microbiome by whole-genome sequencing and resistome analysis: The complexity of the healthy picture. BMC Microbiol. 2020, 20, 120. [Google Scholar] [CrossRef]
  2. Oh, C.; Lee, K.; Cheong, Y.; Lee, S.W.; Park, S.Y.; Song, C.S.; Choi, I.S.; Lee, J.B. Comparison of the Oral Microbiomes of Canines and Their Owners Using Next-Generation Sequencing. PLoS ONE 2015, 10, e0131468. [Google Scholar] [CrossRef] [PubMed]
  3. Willis, J.R.; Gabaldón, T. The Human Oral Microbiome in Health and Disease: From Sequences to Ecosystems. Microorganisms 2020, 8, 308. [Google Scholar] [CrossRef] [PubMed]
  4. Sureda, A.; Daglia, M.; Argüelles Castilla, S.; Sanadgol, N.; Fazel Nabavi, S.; Khan, H.; Belwal, T.; Jeandet, P.; Marchese, A.; Pistollato, F.; et al. Oral microbiota and Alzheimer’s disease: Do all roads lead to Rome? Pharmacol. Res. 2020, 151, 104582. [Google Scholar] [CrossRef] [PubMed]
  5. Koren, O.; Spor, A.; Felin, J.; Fåk, F.; Stombaugh, J.; Tremaroli, V.; Behre, C.J.; Knight, R.; Fagerberg, B.; Ley, R.E.; et al. Human oral, gut, and plaque microbiota in patients with atherosclerosis. Proc. Natl. Acad. Sci. USA 2010, 108 (Suppl. S1), 4592–4598. [Google Scholar] [CrossRef] [PubMed]
  6. Karpiński, T.M. Role of Oral Microbiota in Cancer Development. Microorganisms 2019, 7, 20. [Google Scholar] [CrossRef]
  7. Mager, D.L.; Haffajee, A.D.; Devlin, P.M.; Norris, C.M.; Posner, M.R.; Goodson, J.M. The salivary microbiota as a diagnostic indicator of oral cancer: A descriptive, non-randomized study of cancer-free and oral squamous cell carcinoma subjects. J. Transl. Med. 2005, 3, 27. [Google Scholar] [CrossRef]
  8. Kaeberlein, M.; Creevy, K.E.; Promislow, D.E. The dog aging project: Translational geroscience in companion animals. Mamm. Genome 2016, 27, 279–288. [Google Scholar] [CrossRef]
  9. Ruparell, A.; Inui, T.; Staunton, R.; Wallis, C.; Deusch, O.; Holcombe, L.J. The canine oral microbiome: Variation in bacterial populations across different niches. BMC Microbiol. 2020, 20, 42. [Google Scholar] [CrossRef]
  10. Wallis, C.; Marshall, M.; Colyer, A.; O’Flynn, C.; Deusch, O.; Harris, S. A longitudinal assessment of changes in bacterial community composition associated with the development of periodontal disease in dogs. Vet. Microbiol. 2015, 181, 271–282. [Google Scholar] [CrossRef]
  11. Davis, I.J.; Wallis, C.; Deusch, O.; Colyer, A.; Milella, L.; Loman, N.; Harris, S. A cross-sectional survey of bacterial species in plaque from client owned dogs with healthy gingiva, gingivitis or mild periodontitis. PLoS ONE 2013, 8, e83158. [Google Scholar] [CrossRef]
  12. Isaiah, A.; Hoffmann, A.R.; Kelley, R.; Mundell, P.; Steiner, J.M.; Suchodolski, J.S. Characterization of the nasal and oral microbiota of detection dogs. PLoS ONE 2017, 12, e0184899. [Google Scholar] [CrossRef]
  13. Niemiec, B.A.; Gawor, J.; Tang, S.; Prem, A.; Krumbeck, J.A. The mycobiome of the oral cavity in healthy dogs and dogs with periodontal disease. Am. J. Vet. Res. 2022, 83, 42–49. [Google Scholar] [CrossRef]
  14. Abbayya, K.; Puthanakar, N.Y.; Naduwinmani, S.; Chidambar, Y.S. Association between Periodontitis and Alzheimer’s Disease. N. Am. J. Med. Sci. 2015, 7, 241–246. [Google Scholar] [CrossRef] [PubMed]
  15. Mohanty, R.; Asopa, S.J.; Joseph, M.D.; Singh, B.; Rajguru, J.P.; Saidath, K.; Sharma, U. Red complex: Polymicrobial conglomerate in oral flora: A review. J. Family Med. Prim. Care 2019, 8, 3480–3486. [Google Scholar] [CrossRef] [PubMed]
  16. Dominy, S.S.; Lynch, C.; Ermini, F.; Benedyk, M.; Marczyk, A.; Konradi, A.; Nguyen, M.; Haditsch, U.; Raha, D.; Griffin, C.; et al. Porphyromonas gingivalis in Alzheimer’s disease brains: Evidence for disease causation and treatment with small-molecule inhibitors. Sci. Adv. 2019, 5, eaau3333. [Google Scholar] [CrossRef] [PubMed]
  17. Mihevc, S.P.; Majdič, G. Canine Cognitive Dysfunction and Alzheimer’s Disease—Two Facets of the Same Disease? Front. Neurosci. 2019, 13, 604. [Google Scholar] [CrossRef] [PubMed]
  18. Neilson, J.C.; Hart, B.L.; Cliff, K.D.; Ruehl, W.W. Prevalence of behavioral changes associated with age-related cognitive impairment in dogs. J. Am. Vet. Med. Assoc. 2001, 218, 1787–1791. [Google Scholar] [CrossRef] [PubMed]
  19. Fefer, G.; Panek, W.K.; Khan, M.Z.; Singer, M.; Westermeyer, H.D.; Mowat, F.M.; Murdoch, D.M.; Case, B.; Olby, N.J.; Gruen, M.E. Use of Cognitive Testing, Questionnaires, and Plasma Biomarkers to Quantify Cognitive Impairment in an Aging Pet Dog Population. J. Alzs. Dis. 2022, 87, 1367–1378. [Google Scholar] [CrossRef]
  20. Greer, K.A.; Canterberry, S.C.; Murphy, K.E. Statistical analysis regarding the effects of height and weight on lifespan of the domestic dog. Res. Vet. Sci. 2007, 82, 208–214. [Google Scholar] [CrossRef]
  21. Madari, A.; Farbakova, J.; Katina, S.; Smolek, T.; Novak, P.; Weissovaa, T.; Novak, M.; Zilka, N. Assessment of severity and progression of canine cognitive dysfunction syndrome using the CAnine DEmentia Scale (CADES). Appl. Anim. Behav. Sci. 2015, 171, 138–145. [Google Scholar] [CrossRef]
  22. Brown, D.C.; Boston, R.C.; Coyne, J.C.; Farrar, J.T. Development and psychometric testing of an instrument designed to measure chronic pain in dogs with osteoarthritis. Am. J. Vet. Res. 2007, 68, 631–637. [Google Scholar] [CrossRef] [PubMed]
  23. Brown, D.C.; Boston, R.; Coyne, J.C.; Farrar, J.T. A novel approach to the use of animals in studies of pain: Validation of the canine brief pain inventory in canine bone cancer. Pain Med. 2009, 10, 133–142. [Google Scholar] [CrossRef] [PubMed]
  24. Brown, D.C.; Boston, R.C.; Coyne, J.C.; Farrar, J.T. Ability of the canine brief pain inventory to detect response to treatment in dogs with osteoarthritis. J. Am. Vet. Med. Assoc. 2008, 233, 1278–1283. [Google Scholar] [CrossRef] [PubMed]
  25. Caiafa, A. Oral Examination/Dental Charting and Diagnostic Tools. In World Small Animal Veterinary Association World Congress Proceedings; North Coast Veterinary Specialists: Sippy Downs, QLD, Australia, 2013. [Google Scholar]
  26. Bray, E.E.; Gruen, M.E.; Gnanadesikan, G.E.; Horschler, D.J.; Levy, K.M.; Kennedy, B.S.; Hare, B.A.; MacLean, E.L. Cognitive characteristics of 8- to 10-week-old assistance dog puppies. Anim. Behav. 2020, 166, 193–206. [Google Scholar] [CrossRef] [PubMed]
  27. Allali, I.; Arnold, J.W.; Roach, J.; Cadenas, M.B.; Butz, N.; Hassan, H.M.; Koci, M.; Ballou, A.; Mendoza, M.; Ali, R.; et al. A Comparison of Sequencing Platforms and Bioinformatics Pipelines for Compositional Analysis of the Gut Microbiome. BMC Microbiol. 2017, 17, 194. [Google Scholar] [CrossRef]
  28. Guadamuro, L.; Azcárate-Peril, M.A.; Tojo, R.; Mayo, B.; Delgado, S. Use of high throughput amplicon sequencing and ethidium monoazide dye to track microbiota changes in an equol-producing menopausal woman receiving a long-term isoflavones treatment. AIMS Microbiol. 2019, 5, 102–116. [Google Scholar] [CrossRef]
  29. Azcarate-Peril, M.A.; Butz, N.; Cadenas, M.B.; Koci, M.; Ballou, A.; Mendoza, M.; Ali, R.; Hassan, H. An Attenuated Salmonella enterica Serovar Typhimurium Strain and Galacto-Oligosaccharides Accelerate Clearance of Salmonella Infections in Poultry through Modifications to the Gut Microbiome. Appl. Environ. Microbiol. 2018, 84, e02526-17. [Google Scholar] [CrossRef]
  30. Seth-Smith, H.M.B.; Bonfiglio, F.; Cuénod, A.; Reist, J.; Egli, A.; Wüthrich, D. Evaluation of Rapid Library Preparation Protocols for Whole Genome Sequencing Based Outbreak Investigation. Front. Public Health 2019, 7, 241. [Google Scholar] [CrossRef]
  31. Bcl2Fastq, version 2.20.0; Illumina, Inc.: San Diego, CA, USA, 2018.
  32. FastQC, version 0.11.9; Babraham Institute: Cambridge, UK, 2014.
  33. Trim Galore, version 0.6.7; Babraham Institute: Cambridge, UK, 2017.
  34. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef]
  35. Lu, J.; Breitwieser, F.P.; Thielen, P.; Salzberg, S.L. Bracken: Estimating species abundance in metagenomics data. PeerJ Comput. Sci. 2017, 3, e104. [Google Scholar] [CrossRef]
  36. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2021. Available online: https://www.R-project.org/ (accessed on 24 March 2023).
  37. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  38. McMurdie, P.J.; Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef]
  39. Oksanen, J.; Simpson, G.; Blanchet, F.; Kindt, R.; Legendre, P.; Minchin, P.; O’Hara, R.; Solymos, P.; Stevens, M.; Szoecs, E.; et al. _vegan: Community Ecology Package_. R package version 2.6-2. 2022. Available online: https://CRAN.R-project.org/package=vegan (accessed on 24 March 2023).
  40. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  41. JMP®, version Pro 15; SAS Institute Inc.: Cary, NC, USA, 1989–2021.
  42. Xinyan, Z.; Yu-Fang, P.; Lei, Z.; Boyi, G.; Pendegraft, A.H.; Wenzhuo, Z.; Nengjun, Y. Negative Binomial Mixed Models for Analyzing Longitudinal Microbiome Data. Front. Microbiol. 2018, 9, 1683. [Google Scholar] [CrossRef]
  43. Senhorinho, G.N.; Nakano, V.; Liu, C.; Song, Y.; Finegold, S.M.; Avila-Campos, M.J. Detection of Porphyromonas gulae from subgingival biofilms of dogs with and without periodontitis. Anaerobe 2011, 17, 257–258. [Google Scholar] [CrossRef] [PubMed]
  44. Fiorillo, L.; Cervino, G.; Laino, L.; D’Amico, C.; Mauceri, R.; Tozum, T.F.; Gaeta, M.; Cicciù, M. Porphyromonas gingivalis, Periodontal and Systemic Implications: A Systematic Review. Dent. J. 2019, 7, 114. [Google Scholar] [CrossRef] [PubMed]
  45. O’Flynn, C.; Deusch, O.; Darling, A.E.; Eisen, J.A.; Wallis, C.; Davis, I.J.; Harris, S.J. Comparative Genomics of the Genus Porphyromonas Identifies Adaptations for Heme Synthesis within the Prevalent Canine Oral Species Porphyromonas cangingivalis. Genome Biol. Evol. 2015, 7, 3397–3413. [Google Scholar] [CrossRef]
  46. Feng-Jie, J.; Shuang, H.; Bao-Teng, W.; Long, J. Advances in Genetic Engineering Technology and Its Application in the Industrial Fungus Aspergillus oryzae. Front. Microbiol. 2021, 12, 644404. [Google Scholar] [CrossRef]
  47. Moran, G.P.; Coleman, D.C.; Sullivan, D.J. Candida albicans versus Candida dubliniensis: Why Is C. albicans More Pathogenic? Int. J. Microbiol. 2012, 2012, 205921. [Google Scholar] [CrossRef]
  48. Nguyen, L.D.; Viscogliosi, E.; Delhaes, L. The lung mycobiome: An emerging field of the human respiratory microbiome. Front. Microbiol. 2015, 6, 89. [Google Scholar] [CrossRef]
  49. Dean, R.; Van Kan, J.A.L.; Pretorius, Z.A.; Hammons-Kosack, K.E.; Di Pietro, A.; Spanu, P.D.; Rudd, J.J.; Dickman, M.; Kahmann, R.; Ellis, J.; et al. The Top 10 fungal pathogens in molecular plant pathology. Molec. Plant Patho. 2012, 13, 414–430. [Google Scholar] [CrossRef]
  50. Lenzo, J.C.; O’Brien-Simpson, N.M.; Orth, R.K.; Mitchell, H.L.; Dashper, S.G.; Reynolds, E.C. Porphyromonas gulae Has Virulence and Immunological Characteristics Similar to Those of the Human Periodontal Pathogen Porphyromonas gingivalis. Infect Immun. 2016, 84, 2575–2585. [Google Scholar] [CrossRef]
  51. Dewey, C.W.; Rishniw, M. Periodontal disease is associated with cognitive dysfunction in aging dogs: A blinded prospective comparison of visual periodontal and cognitive questionnaire scores. Open Vet. J. 2021, 11, 210–216. [Google Scholar] [CrossRef]
  52. Arastu-Kapur, S.; Nguyen, M.; Raha, D.; Ermini, F.; Haditsch, U.; Araujo, J.; De Lannoy, I.; Ryder, M.I.; Dominy, S.S.; Lynch, C.; et al. Treatment of Porphyromonas gulae infection and downstream pathology in the aged dog by lysine-gingipain inhibitor COR388. Pharmacol. Res. Perspect. 2020, 8, e00562. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, C.; Hemme, C.; Beleno, J.; Shi, Z.J.; Ning, D.; Qin, Y.; Tu, Q.; Jorgensen, M.; He, Z.; Wu, L.; et al. Oral microbiota of periodontal health and disease and their changes after nonsurgical periodontal therapy. ISME J. 2018, 12, 1210–1224. [Google Scholar] [CrossRef] [PubMed]
  54. Bathini, P.; Foucras, S.; Dupanloup, I.; Imeri, H.; Perna, A.; Berruex, J.L.; Doucey, M.A.; Annoni, J.M.; Auber Alberi, L. Classifying dementia progression using microbial profiling of saliva. Alzheimers Dement. 2020, 12, e12000. [Google Scholar] [CrossRef] [PubMed]
  55. Liu, X.X.; Jiao, B.; Liao, X.X.; Guo, L.N.; Yuan, Z.H.; Wang, X.; Xiao, X.W.; Zhang, X.Y.; Tang, B.S.; Shen, L. Analysis of Salivary Microbiome in Patients with Alzheimer’s Disease. J. Alzheimers Dis. 2019, 72, 633–640. [Google Scholar] [CrossRef] [PubMed]
  56. Muñana, K.R.; Jacob, M.E.; Callahan, B.J. Evaluation of fecal Lactobacillus populations in dogs with idiopathic epilepsy: A pilot study. Anim. Microbiome. 2020, 2, 19. [Google Scholar] [CrossRef] [PubMed]
  57. Coman, M.M.; Verdenelli, M.C.; Cecchini, C.; Belà, B.; Gramenzi, A.; Orpianesi, C.; Cresci, A.; Silvi, S. Probiotic characterization of Lactobacillus isolates from canine faeces. J. Appl. Microbiol. 2019, 126, 1245–1256. [Google Scholar] [CrossRef]
  58. Lee, H.J.; Cho, J.H.; Cho, W.J.; Gang, S.H.; Park, S.H.; Jung, B.J.; Kim, H.B.; Song, K.H. Effects of Synbiotic Preparation Containing Lactobacillus gasseri BNR17 on Body Fat in Obese Dogs: A Pilot Study. Animals 2022, 12, 642. [Google Scholar] [CrossRef]
  59. Nakagawa, H.; Shiozaki, T.; Kobatake, E.; Hosoya, T.; Moriya, T.; Sakai, F.; Taru, H.; Miyazaki, T. Effects and mechanisms of prolongevity induced by Lactobacillus gasseri SBT2055 in Caenorhabditis elegans. Aging Cell 2016, 15, 227–236. [Google Scholar] [CrossRef] [PubMed]
  60. Yun, S.W.; Kim, J.K.; Lee, K.E.; Oh, Y.J.; Choi, H.J.; Han, M.J.; Kim, D.H. A Probiotic Lactobacillus gasseri Alleviates Escherichia coli-Induced Cognitive Impairment and Depression in Mice by Regulating IL-1β Expression and Gut Microbiota. Nutrients 2020, 12, 3441. [Google Scholar] [CrossRef]
  61. Liu, S.; Wang, Y.; Zhao, L.; Sun, X.; Feng, Q. Microbiome succession with increasing age in three oral sites. Aging 2020, 12, 7874–7907. [Google Scholar] [CrossRef] [PubMed]
  62. EHead, N.W.M.; Cotman, C.W. Neurobiological models of aging in the dog and other vertebrate species. In Functional Neurobiology of Aging; Hof, P.R., Mobbs, C.V., Eds.; Academic Press: San Diego, CA, USA, 2001; pp. 457–468. [Google Scholar]
  63. Wallis, C.; Holcombe, L.J. A review of the frequency and impact of periodontal disease in dogs. J. Sm. Anim. Pract. 2020, 61, 529–540. [Google Scholar] [CrossRef] [PubMed]
Animals 13 03846 g001
Figure 1. Progression of CADES score over study by subject. Subjects with CADES above 23 (black line) are considered to have moderate to severe (MdSv) cognitive dysfunction. Scores below 23 are considered normal cognition to mild (NM) dysfunction.
Figure 1. Progression of CADES score over study by subject. Subjects with CADES above 23 (black line) are considered to have moderate to severe (MdSv) cognitive dysfunction. Scores below 23 are considered normal cognition to mild (NM) dysfunction.
Animals 13 03846 g001
Animals 13 03846 g002
Figure 2. Relative abundance of bacteria by phyla in combined samples. The dominance of Bacteroidetes and Proteobacteria is clearly visible in this pie chart. Actinobacteria and Firmicutes account for nearly 7% of the bacterial phyla present.
Figure 2. Relative abundance of bacteria by phyla in combined samples. The dominance of Bacteroidetes and Proteobacteria is clearly visible in this pie chart. Actinobacteria and Firmicutes account for nearly 7% of the bacterial phyla present.
Animals 13 03846 g002
Animals 13 03846 g003
Figure 3. Relative abundance of Bacteroidetes and Proteobacteria in each sample from each dog over time. It is evident that in all but dog G06, S1 these 2 bacterial phyla account for the majority of bacteria present. * Indicates that dental prophylaxis was undertaken within the 6-month period prior to oral sampling. Blue bars represent Bacteroidetes and red represent Proteobacteria.
Figure 3. Relative abundance of Bacteroidetes and Proteobacteria in each sample from each dog over time. It is evident that in all but dog G06, S1 these 2 bacterial phyla account for the majority of bacteria present. * Indicates that dental prophylaxis was undertaken within the 6-month period prior to oral sampling. Blue bars represent Bacteroidetes and red represent Proteobacteria.
Animals 13 03846 g003
Animals 13 03846 g004
Figure 4. Simpson’s Index of Diversity (1-D) for bacteria increases in all subjects from first to last sample. This increase in bacterial alpha diversity over time was significant (p = 0.002).
Figure 4. Simpson’s Index of Diversity (1-D) for bacteria increases in all subjects from first to last sample. This increase in bacterial alpha diversity over time was significant (p = 0.002).
Animals 13 03846 g004
Table 1. Physical characteristics of study subjects at time of enrollment in trial. M/N: male neutered; F/S: female spayed.
Table 1. Physical characteristics of study subjects at time of enrollment in trial. M/N: male neutered; F/S: female spayed.
Study IDAge (years)Weight (kg)SexBreed
G0211.622.2M/NBorder Collie
G039.533.7F/SLabrador Retriever
G0414.89.9F/SCorgi, Pembroke
G0510.829.7M/NGerman Shepherd Dog
G0611.110.3M/NJack Russell Terrier
G0711.425.4F/SGolden Retriever
G1110.127.6F/SMixed Breed (Labrador Mix)
G1311.011.1M/NMixed Breed (Boxer/Pitbull)
G1711.811.1M/NBeagle
G249.98M/NDachshund
Table 2. Dental health at time of sampling. Calculus index (CI), gingival index (GI), and plaque index (PI) were recorded at time of microbiome sampling. (nr, not recorded). For the 4 dogs that underwent dental prophylaxis within 6 months of sample collection, the number of months prior is provided.
Table 2. Dental health at time of sampling. Calculus index (CI), gingival index (GI), and plaque index (PI) were recorded at time of microbiome sampling. (nr, not recorded). For the 4 dogs that underwent dental prophylaxis within 6 months of sample collection, the number of months prior is provided.
Study IDSampleCalculus Index (CI)Gingival Index (GI)Plaque Index (PI)Prior Dental Prophylaxis (months)
G021nrnrnr
2223
3212
G0311113.7
2112
3111
G0410210.7
2nrnrnr
3112
G051212
2111
3111
G061212
2212
3231
G071313
2222
3223
G1111012.4
2111
31221.6
G131111
2120
G1712220.5
2121
G241nrnrnr
2112
Table 3. Relative abundance of bacterial phyla, listed in order of greatest to least total abundance, by sample. Bacteroidetes and Proteobacteria account for more than 70% of phyla present in each sample.
Table 3. Relative abundance of bacterial phyla, listed in order of greatest to least total abundance, by sample. Bacteroidetes and Proteobacteria account for more than 70% of phyla present in each sample.
IDSample NumberBacteroidetesProteo-bacteriaFirm-icutesActino-bacteriaSpiro-chaetesTener-icutesFuso-bacteriaCandi-datusChloro-flexiCyano-bacteriaOther Phyla
G02S20.39910.36790.09860.02440.06740.00040.03050.00170.00190.00100.0071
S30.39660.44510.04130.04280.03810.00100.02850.00150.00100.00150.0026
G03S10.40310.48720.01960.01530.01130.00040.06160.00040.00000.00060.0005
S20.46140.41300.04560.02820.01940.00140.01840.00220.00260.00230.0054
S30.36620.50520.03860.03180.01760.00070.02750.00280.00280.00170.0049
G04S10.27930.60990.02390.05690.00650.00130.02130.00030.00000.00010.0003
S20.38310.37230.07600.03770.08730.00280.01880.00310.01010.00230.0066
S30.44600.36810.04390.02420.08000.00310.02120.00210.00560.00140.0044
G05S10.35650.60390.01890.00290.00540.00060.00180.00970.00000.00010.0002
S20.59510.18020.09340.02710.03680.00570.04850.00750.00140.00140.0028
S30.44340.30160.03290.11970.06550.00120.01460.00120.01040.00170.0079
G06S10.07190.19770.01940.00850.02760.66930.00290.00060.00000.00070.0014
S20.38720.47050.02560.04900.04100.00100.01850.00300.00170.00060.0020
S30.58160.28030.05330.01680.05030.00130.00990.00250.00180.00070.0015
G07S10.12120.76740.04780.04170.01090.00190.00680.00000.00000.00000.0023
S20.34380.48820.03230.03210.09040.00080.00890.00140.00020.00090.0009
S30.55970.26160.03500.05880.04900.00180.02120.00360.00480.00130.0033
G11S10.43220.49600.04400.00940.00210.00000.00050.01550.00000.00000.0002
S20.68250.21080.03280.02950.01140.00040.02470.00280.00310.00040.0015
S30.44230.41290.03530.07550.00690.00080.01030.00410.00560.00180.0045
G13S10.66890.26300.02160.01370.00960.00250.01790.00180.00010.00040.0005
S20.44900.40460.05600.06940.00180.00090.01470.00060.00070.00020.0021
G17S10.59340.35030.00940.01290.00400.00110.02780.00000.00010.00020.0007
S20.65730.28350.01390.01290.00040.00020.03030.00010.00020.00030.0009
G24S10.73040.12620.07600.00800.01860.00240.03310.00330.00050.00030.0011
S20.49560.32850.02300.10020.02140.00090.02390.00130.00170.00060.0029
Table 4. Relative abundance of bacterial species. Twenty most abundant oral bacterial species in the combined samples. Porphorymonas species are the most prevalent bacteria in these samples.
Table 4. Relative abundance of bacterial species. Twenty most abundant oral bacterial species in the combined samples. Porphorymonas species are the most prevalent bacteria in these samples.
SpeciesRelative Abundance
Porphyromonas gulae0.116
Porphyromonas cangingivalis0.105
Bergeyella zoohelcum0.076
Capnocytophaga canimorsus0.054
Pasteurella multocida0.047
Tannerella forsythia0.034
Capnocytophaga cynodegmi0.033
Porphyromonas gingivalis0.030
Neisseria weaveri0.028
Desulfomicrobium orale0.024
Frederiksenia canicola0.021
Treponema denticola0.016
Conchiformibius steedae0.016
Fusobacterium russii0.015
Neisseria shayeganii0.013
Capnocytophaga sp. H29310.013
Capnocytophaga sp. H43580.012
Campylobacter sp. CCUG 573100.010
Porphyromonas crevioricanis0.007
Pseudomonas aeruginosa0.006
Table 5. Longitudinal analysis of bacterial species, fraction of lifespan and CADES. Bacterial species identified as statistically significant using the Negative Binomial Mixed Model. FoL: Fraction of Lifespan. CADES; Canine Dementia Scale. Three different models are included, one each for FoL and CADES and a third in which both FoL and CADES were included. Lactobacillus gasseri is significantly associated with FoL when considered alone or with CADES in the model, and Leptotrichia sp. are significant associated with CADES when modeled without FoL. Statistically significant adjusted p values are in bold.
Table 5. Longitudinal analysis of bacterial species, fraction of lifespan and CADES. Bacterial species identified as statistically significant using the Negative Binomial Mixed Model. FoL: Fraction of Lifespan. CADES; Canine Dementia Scale. Three different models are included, one each for FoL and CADES and a third in which both FoL and CADES were included. Lactobacillus gasseri is significantly associated with FoL when considered alone or with CADES in the model, and Leptotrichia sp. are significant associated with CADES when modeled without FoL. Statistically significant adjusted p values are in bold.
InterceptFoLCADESPool
ModelTaxaCoefficientadj. p-ValueCoefficientadj. p-ValueCoefficientadj. p-ValueCoefficientadj. p-Value
FoL + PoolLactobacillus gasseri14.64502.4732 × 10−73.85860.0006NANA1.92960.5512
CADES + PoolLeptotrichia sp. oral taxon 212−9.28082.6547 × 10−16NANA2.25532.9289 × 10−5−0.58530.3086
FoL + CADES + PoolLeptotrichia sp. oral taxon 212−10.16953.8961 × 10−14−0.24270.92830.92700.3838−0.64150.5627
Lactobacillus gasseri−15.01342.0346 × 10−63.38700.00011.42960.29992.23780.5627
Table 6. Oral fungal species. Twenty most abundant fungal species in total. * Represents species found in all 26 samples.
Table 6. Oral fungal species. Twenty most abundant fungal species in total. * Represents species found in all 26 samples.
SpeciesRelative Abundance
Aspergillus oryzae *0.6699
Colletotrichum higginsianum *0.0964
Ustilago maydis *0.0182
Botrytis cinerea *0.0168
Eremothecium sinecaudum0.0102
Pochonia chlamydosporia *0.0101
Talaromyces rugulosus *0.0081
Candida dubliniensis *0.0080
Kluyveromyces marxianus *0.0078
Thermothielavioides terrestris *0.0076
Zymoseptoria tritici0.0074
Fusarium verticillioides *0.0068
Neurospora crassa0.0067
Ustilaginoidea virens0.0066
Aspergillus luchuensis *0.0065
Thermothelomyces thermophilus *0.0061
Fusarium pseudograminearum0.0057
Pyricularia oryzae0.0056
Naumovozyma dairenensis0.0052
Drechmeria coniospora0.0050
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Templeton, G.B.; Fefer, G.; Case, B.C.; Roach, J.; Azcarate-Peril, M.A.; Gruen, M.E.; Callahan, B.J.; Olby, N.J. Longitudinal Analysis of Canine Oral Microbiome Using Whole Genome Sequencing in Aging Companion Dogs. Animals 2023, 13, 3846. https://doi.org/10.3390/ani13243846

AMA Style

Templeton GB, Fefer G, Case BC, Roach J, Azcarate-Peril MA, Gruen ME, Callahan BJ, Olby NJ. Longitudinal Analysis of Canine Oral Microbiome Using Whole Genome Sequencing in Aging Companion Dogs. Animals. 2023; 13(24):3846. https://doi.org/10.3390/ani13243846

Chicago/Turabian Style

Templeton, Ginger B., Gilad Fefer, Beth C. Case, Jeff Roach, M. Andrea Azcarate-Peril, Margaret E. Gruen, Benjamin J. Callahan, and Natasha J. Olby. 2023. "Longitudinal Analysis of Canine Oral Microbiome Using Whole Genome Sequencing in Aging Companion Dogs" Animals 13, no. 24: 3846. https://doi.org/10.3390/ani13243846

APA Style

Templeton, G. B., Fefer, G., Case, B. C., Roach, J., Azcarate-Peril, M. A., Gruen, M. E., Callahan, B. J., & Olby, N. J. (2023). Longitudinal Analysis of Canine Oral Microbiome Using Whole Genome Sequencing in Aging Companion Dogs. Animals, 13(24), 3846. https://doi.org/10.3390/ani13243846

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop