http://orcid.org/0000-0002-0684-7283
Figures
Abstract
The ocular surface microbiome of veterinary species has not been thoroughly characterized using molecular-based techniques, such as next generation sequencing (NGS), as the vast majority of studies have utilized traditional culture-based techniques. To date, there is one pilot study evaluating the ocular surface of healthy dogs using NGS. Furthermore, alterations in the ocular surface microbiome over time and after topical antibiotic treatment are unknown. The objectives of this study were to describe the bacterial composition of the ocular surface microbiome in clinically normal dogs, and to determine if microbial community changes occur over time or following topical antibiotic therapy. Topical neomycin-polymyxin-bacitracin ophthalmic ointment was applied to one eye each of 13 adult dogs three times daily for seven days, while contralateral eyes served as untreated controls. The inferior conjunctival fornix of both eyes was sampled via swabbing at baseline prior to antibiotic therapy (day 0), after 1 week of treatment (day 7), and 4 weeks after discontinuing treatment (day 35). Genomic DNA was extracted from the conjunctival swabs and primers targeting the V4 region of bacterial 16S rRNA genes were used to generate amplicon libraries, which were then sequenced on an Illumina platform. Data were analyzed using Quantitative Insights Into Molecular Ecology (QIIME 2.0). At baseline, the most relatively abundant phyla sequenced were Proteobacteria (49.7%), Actinobacteria (25.5%), Firmicutes (12%), Bacteroidetes (7.5%), and Fusobacteria (1.4%). The most common families detected were Pseudomonadaceae (13.2%), Micrococcaceae (12%), Pasteurellaceae (6.9%), Microbacteriaceae (5.2%), Enterobacteriaceae (3.9%), Neisseriaceae (3.5%), and Corynebacteriaceae (3.3%). Alpha and beta diversity measurements did not differ in both control and treatment eyes over time. This report examines the temporal stability of the canine ocular surface microbiome. The major bacterial taxa on the canine ocular surface remained consistent over time and following topical antibiotic therapy.
Citation: Rogers CM, Scott EM, Sarawichitr B, Arnold C, Suchodolski JS (2020) Evaluation of the bacterial ocular surface microbiome in ophthalmologically normal dogs prior to and following treatment with topical neomycin-polymyxin-bacitracin. PLoS ONE 15(6): e0234313. https://doi.org/10.1371/journal.pone.0234313
Editor: Simon Clegg, University of Lincoln, UNITED KINGDOM
Received: October 23, 2019; Accepted: May 24, 2020; Published: June 9, 2020
Copyright: © 2020 Rogers et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession number SRP161472.
Funding: The authors received no external funding for this work. Internal funding from Texas A&M in the form of start-up for EMS were received.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The ocular surface microbiota consists of microorganisms residing on the corneoconjunctival surface and within the tear film. The canine eye is susceptible to ocular surface diseases such as keratoconjunctivitis sicca and infectious ulcerative keratitis, which are often treated with topical broad-spectrum antibiotics such as neomycin-polymyxin-bacitracin [1–7]. Changes in the composition of the ocular surface microbiota may be associated with ocular surface disease, as evidence suggests these commensal microorganisms protect the eye against the proliferation of opportunistic and pathogenic species [8,9]. In addition, antibiotic usage may negatively impact the structure and stability of microbial communities [10–13].
The vast majority of previous studies describing the bacterial ocular surface communities of healthy dogs used traditional culture-based methods [3,5,7,14–16]. The percent of total positive cultures from healthy canine eyes was low, ranging from 29–45%, with Gram-positive bacteria such as Staphylococcus, Streptococcus, Bacillus, and Micrococcus spp. most commonly isolated regardless of geography or climate [3,5,7,14–16]. Gram-negative bacteria were less frequently cultured, with Pseudomonas, Moraxella, Acinetobacter, Neisseria, E. coli, Klebsiella, and Enterococcus spp. comprising less than 4% of isolates [3,5,7,14–16].
Of the traditional culture-based studies involving the canine eye, few evaluated the effect of seasonality or topical antibiotic use on the frequency or type of bacteria sampled from the ocular surface [4,7,12,13]. The percent of positive bacterial cultures from the canine ocular surface were increased during spring and summer compared to autumn and winter, with July having the highest percent of total positive cultures (60%) [7]. Use of the topical fluoroquinolone antibiotic, ofloxacin, for three weeks following cataract surgery in 16 dogs altered the composition of the bacterial community cultured from their conjunctival surfaces, producing >10 colonies per streak compared to <5 colonies per streak prior to treatment [12]. Additionally, increased colony formation was mirrored by increased resistance to topical fluoroquinolones, with the highest percentage of resistant organisms identified 3 weeks post-operatively [12].
There are limitations to culture-dependent studies as cultivable bacteria only represent a fraction of the organisms present in a community. Molecular-based methods, such as 16S rRNA gene sequencing, allow for a more accurate overview of the diversity of bacterial populations colonizing the ocular surface of humans [17–23], horses [24], cats [25–27] and dogs [28].
In the single previous study describing the ocular surface microbiome in dogs using culture-independent molecular-based methods, samples from both eyes were taken at a solitary time point from 10 research-bred dogs [28]. The most common phyla identified were Firmicutes, Actinobacteria, Proteobacteria, and Bacteroidetes. Additionally, the most commonly cultured bacterial genera from canine eyes, Staphylococcus, Streptococcus, and Bacillus spp., comprised only 2.63% of the bacterial community sequenced from the ocular surface [28]. While this preliminary report provides a detailed description of the organisms present, sampling from a single time point limits the usefulness of data to a sole set of environmental and laboratory conditions that may not be repeatable over time. This is compounded by the usage of subjects from a research colony that share many of the same living and dietary conditions, as well as physical characteristics. Samples collected from the previous study, while relevant to a baseline description of organisms present in canines, may not be clinically applicable to dogs living and interacting in diverse environments such as with privately-owned pets.
This study aimed to describe the bacterial compostition of the ocular surface microbiome in healthy privately-owned dogs, and to determine if microbial community changes occur over time or following topical antibiotic therapy. Continued investigations on the effects of antibiotics and other influences on the ocular suface microbiome are necessary to improve our understanding of ocular diseases in both veterinary and physician ophthalmology.
Materials and methods
Participants
The study was approved, with owner consent, by the Texas A&M University Institutional Animal Care and Use Committee (Animal Use Protocol #2018–0079). Thirteen dogs, free of ocular disease, were admitted to the Texas A&M Veterinary Medical Teaching Hospital on an out-patient basis for ocular examination and conjunctival swabbing (Table 1). All dogs in the study were recruited by a university-wide email and owned by faculty, students, and staff of Texas A&M College of Veterinary Medicine & Biomedical Sciences. Dogs 8 and 9 as well as dogs 10 and 11 were respectively housed together. The study was performed in May and June in east-central Texas. Sample size was determined from previously published reports indicating 10–12 animals per control and experimental groups would be adequate [29–31].
Sample collection
All dogs had a complete ophthalmic examination performed by a board-certified veterinary ophthalmologist (EMS), as previously described [24,27]. A routine minimal ophthalmic database was performed, consisting of Schirmer tear test measurements (Intervet Inc., Summit, NJ), fluorescein staining (Amcon Laboratories Inc., St. Louis, MO), and tonometry (TonoVet, Jorgensen Laboratories, Loveland, CO). Any dog with evidence of ocular disease was excluded from the study.
Baseline conjunctival swab samples were collected after the Schirmer tear test and before tonometry and fluorescein staining in order to prevent sample dilution or contamination, as previously described [24,27]. One drop of 0.5% proparacaine (Bausch & Lomb Inc., Bridgewater, NJ) was placed on the ocular surface of each eye to provide topical anesthesia. The inferior conjunctival fornix of both eyes was sampled with Isohelix buccal swabs (Boca Scientific, Inc. Westwood, MA), as previously described [24,27]. Two swabs were used per eye, and each side of the swab was rubbed in the conjunctival fornix 10 times. Swabs were then collected in DNeasy Powerbead tubes with 750 μl buffer containing guanidine thiocyanate (QIAGEN, Inc., Germantown, MD). One drop of 0.5% proparacaine was placed on an unused swab and inserted into an empty PowerBead tube at each of the three time points to serve as three negative controls and confirm lack of environmental contamination. All samples were immediately stored at 4 degrees C and extractions were performed within 24 hours.
Following the collection of baseline samples, one eye of each dog was randomly selected for treatment with topical neomycin-polymyxin B-bacitracin (Dechra Veterinary Products, Overland Park, KS). Online software (https://www.randomizer.org) was employed to randomize left or right eyes into treatment and control groups for each dog. Owners were instructed to instill 1/4-inch strip of the triple antibiotic ointment directly to the ocular surface of the selected eye of their dog three times daily for 7 days. Treatment sheets with detailed instructions and a checklist for treatment application times were distributed to the owners and returned to the investigators at the end of the study to confirm medication compliance (S1 Checklist). Owners returned to the Veterinary Medical Teaching Hospital with their dogs for repeat conjunctival swabs that occurred at the completion of topical antimicrobial therapy (Day 7), and 4 weeks after discontinuing therapy (Day 35). On day 7, sampling occurred approximately 4 hours following the last instillation of medication.
DNA extraction and sequencing
Genomic DNA was extracted from conjunctival swabs and three negative controls using a single 100 tube DNeasy Powersoil DNA isolation kit (QIAGEN, Inc., Germantown, MD) according to the manufacturer’s instructions. Sequencing of the 16S rRNA gene V4 variable region was performed at MR DNA Laboratory (www.mrdnalab.com, Shallowater, TX, USA) on an Illumina MiSeq platform (Illumina Inc., San Diego, CA) to construct 2x300 paired-end reads using 515F (5’ -GTGYCAGCMGCCGCGGTAA- 3’) and 806R (5´-GGACTACNVGGGTWTCTAAT- 3´) primers [24,27]. The 16S rRNA gene V4 variable region PCR primers 515/806 with barcode on the forward primer were used in a 35 cycle PCR using the HotStarTaq Plus Master Mix Kit (QIAGEN, Inc., Germantown, MD) under the following conditions: 94°C for 3 minutes, followed by 30–35 cycles of 94°C for 30 seconds, 53°C for 40 seconds and 72°C for 1 minute, after which a final elongation step at 72°C for 5 minutes was performed. After amplification, PCR products were checked in 2% agarose gel to determine the success of amplification and the relative intensity of bands. Multiple samples were pooled together (e.g., 100 samples) in equal proportions based on their molecular weight and DNA concentrations. Pooled samples were purified using calibrated Ampure XP beads. Then the pooled and purified PCR product was used to prepare Illumina Truseq nano DNA library.
Data analysis
Statistical analysis was performed as previously described [24,27]. Sequences were processed and analyzed using Quantitative Insights Into Microbial Ecology (QIIME2 2018.6) [32]. Raw sequences were demultiplexed and the amplicon sequence variant (ASV) table was created using DADA2 [33]. Prior to downstream analysis, sequences assigned as chloroplast, mitochondria, and low abundance ASVs, containing less than 0.01% of the total reads in the dataset were removed. Data were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession number SRP161472.
Alpha diversity (observed ASVs, Shannon, Chao1) was calculated in QIIME2 and analyzed to compare species richness and evenness between control and treatment eyes at baseline and among control and treatment eyes over time. Data were tested for normality using a Shapiro-Wilk test (JMP Pro 14, SAS, Marlow, Buckinghamshire), and followed a non-normal distribution. Therefore, a non-parametric Wilcoxon matched-pairs signed-ranks test was applied for comparison between treatment and control eyes at baseline (PRISM 7, GraphPad Software Inc., San Diego, CA). A non-parametric Friedman test, followed by a Dunn’s multiple comparison post-test were utilized to assess differences in treatment and control eyes over three time points [34].
Beta diversity was calculated in QIIME2 to compare bacterial community composition between samples, and evaluated with the weighted and unweighted phylogeny-based UniFrac distance metric [35] and visualized using Principal Coordinate Analysis (PCoA) plots. An Analysis of Similarity test (ANOSIM) was performed within PRIMER 6 software (PRIMER-E Ltd. Luton, UK) to assess differences in bacterial community composition between samples, where an R statistic near 1 indicates a difference in composition while a value near 0 indicates no difference in composition. Unweighted Unifrac distances over time in both treatment and control eyes were compared using a two-way repeated measures ANOVA (PRISM 7, GraphPad Software Inc., San Diego, CA).
Differences in bacterial taxa relative abundance between eyes at baseline, and among control and treatment eyes over time, were explored. Data did not meet the assumptions for normality using a Shapiro-Wilk test (JMP Pro 14, SAS, Marlow, Buckinghamshire). Therefore, a non-parametric Mann-Whitney U test was elected for statistical comparison between treatment and control eyes at baseline; and a non-parametric Friedman test was used to assess differences in treatment and control eyes over three time points (PRISM 7, GraphPad Software Inc., San Diego, CA). A Dunn’s multiple comparison post-test was then used to determine which time points were significantly different. P-values were adjusted for multiple comparisons and corrected for false discovery rate [36]. P- and q-values <0.05 were considered statistically significant. Lastly, the association of bacterial taxa abundance with each time point in both treatment and control eyes was analyzed by linear discriminant analysis effect size (LEfSe) using Calypso software [37, 38]. This additional step helps to further detect bacterial organisms and functional characteristics differentially abundant between two or more microbial environments [37].
Results
Sequence analysis
A total of three negative controls, collected at each time point, failed to show amplification on PCR, indicating the sampling and DNA extraction processes were not contaminated. Sample controls (unused swabs that were processed after the extraction protocol) were also included in sequencing and did not generate any data. All 78 samples (conjunctival swabs collected from 26 eyes at three time points) were positive for PCR amplification and yielded sufficient quality sequences. A total of 2,574,531 sequences were amplified (Min: 13,338.00; Max: 61,820.00; Median: 31,079.50; Mean: 33,006.81; SD: 11,570.95), and rarified to an even sequencing depth, based on the lowest read depth of samples, to 13,338 sequences per sample (S1 Fig). The relative abundance of bacteria was defined for each individual sample.
Healthy canine eyes at baseline
Species richness and diversity at baseline.
Species richness (observed ASVs and Chao1) and abundance and evenness (Shannon) were analyzed to examine taxonomic diversity within a sample. There was no significant difference between control eyes and treatment eyes at baseline (prior to antibiotic treatment) for all three alpha diversity metrics (S1 Table). Hence, eyes at baseline exhibited similar results with regard to species richness, evenness, or abundance (Fig 1A).
Scatter plots comparing alpha diversity results from (A) treatment and control eyes at baseline (day 0), (B) control eyes of healthy dogs (n = 13) at 3 time points: day 0, day 7, day 35, and (C) treatment eyes of healthy dogs (n = 13) at 3 time points: baseline (day 0), following one week of topical antibiotic therapy (day 7), four weeks after discontinued topical antibiotic therapy (day 35). Each dot corresponds to one eye from 13 healthy dogs. There is no difference in alpha diversity between control and treatment eyes at baseline (Wilcoxon match-pairs signed-ranks test). Alpha diversity did not differ in control eyes or treatment eyes over time (Friedman test and Dunn’s post-test).
Microbial community structure at baseline.
Beta diversity measures were analyzed to examine taxonomic diversity between samples. Weighted and unweighted UniFrac distance matrices showed no difference in community structure between treatment and control eyes at baseline (R = -0.06, R = -0.05, respectively, p > 0.05), with a lack of clustering observed on principal coordinate analysis plots (PCoA) (Fig 2). Weighted UniFrac PCoA plots gave us similar results (data not shown). Clustering was also not apparent by individual dog or shared housing (S2 Fig and S3 Fig).
Principal coordinate analysis (PCoA) plot based on the unweighted UniFrac distance metric illustrating no difference in beta diversity by a lack of clustering between eyes at baseline. Each dot corresponds to the microbial composition of one eye.
Microbial community composition at baseline.
The relative abundance of bacterial taxa did not differ between treatment and control eyes at baseline (Mann-Whitney U test). Results from all sampled eyes were averaged to describe the bacterial composition of the ocular surface from 13 healthy dogs. A total of 10 bacterial phyla were identified, with 4 phyla consistently present in all 26 eyes (S2 Table). The most prevalent phyla detected were Proteobacteria (49.7%), followed by Actinobacteria (25.5%), Firmicutes (12.0%), and Bacteroidetes (7.5%) (S2 Table, Fig 3).
Average relative taxa abundances from the ocular surface of healthy dogs, annotated to the level of bacterial phylum, at baseline (day 0). Each bar chart represents the left (OS) or right (OD) eyes of 13 dogs.
A total of 46 bacterial families were identified at ≥1% relative abundance, and 14 families were present in at least 20/26 eyes. Only two families, Microbacteriaceae and Pseudomonadaceae, were present in all 26 eyes (S2 Table). The most prevalent bacterial families identified were Pseudomonadaceae (13.2%), Micrococcaceae (12.0%), Pasteurellaceae (6.9%), Microbacteriaceae (5.2%), Enterobacteriaceae (3.9%), Neisseriaceae (3.5%), and Corynebacteriaceae (3.3%) (S2 Table, Fig 4). Frequently cultured bacterial families from the ocular surface of healthy dogs, such as Micrococcaceae, Corynebacteriaceae, Staphylococcaceae, Bacillaceae, and Streptococcaceae, represented 12.0%, 3.3%, 2.3%, 1.5% and 1.2% of the bacterial families sequenced, respectively. Individual variation in relative abundances of bacterial taxa was observed both between eyes and between dogs (Figs 3 and 4). Throughout all samples, an average of 107 different ASVs were sequenced.
Average relative taxa abundances from the ocular surface of healthy dogs, annotated to the level of bacterial family, at baseline (day 0). Each bar represents the left (OS) or right (OD) eyes of 13 dogs.
Temporal variability in control eyes
In order to investigate temporal stability of the ocular surface microbiome in clinically normal dogs, additional conjunctival swabs were collected from control eyes one week (day 7) and five weeks (day 35) after the initial samples (day 0, baseline).
Species richness and diversity over time.
Alpha diversity metrics were unchanged in control eyes between the three sampled time points (S3 Table and Fig 1B).
Microbial community structure over time.
There was no significant change in beta diversity associated with time, evidenced by the lack of clustering in PCoA plots (Fig 5A). The global community structure of control eyes did not vary over time based on ANOSIM (weighted UniFrac, R = 0.01, R = 0.06, R = -0.03 for day 0 vs. 7, day 0 vs. 35, and day 7 vs. 35, respectively, p > 0.05); (unweighted UniFrac, R = -0.03, R = -0.003, R = 0.01 for day 0 vs. 7, day 0 vs. 35, and day 7 vs. 35, respectively, p > 0.05). Weighted UniFrac PCoA plots gave us similar results (data not shown). There was no difference in unweighted Unifrac distances over time in control eyes, when measured from baseline (day 0) (p > 0.05) (Fig 6).
Principal coordinate analysis comparing (A) control eyes over time, and (B) treatment eyes over time. Principal coordinate analysis (PCoA) plot based on the unweighted UniFrac distance metric illustrating no difference in beta diversity by a lack of clustering among (A) 13 control eyes and (B) 13 treatment eyes at three time points: day 0, day 7, day 35. Each dot corresponds to the microbial composition of one eye.
There is no significant difference in Unifrac distances measured from baseline (day 0), suggesting bacterial communities remain stable over time (Two-way ANOVA; p > 0.05).
Microbial community composition over time.
Fig 7A and 7B displays the mean relative abundance of bacteria in control eyes at each of the time points sampled. Variance in bacterial abundance was observed between individual dogs (S4 Fig). No significant changes in abundance were detected at the phylum, family, or genus level (q > 0.05) (Table 2).
Bacterial phylum and family composition of (A, B) control eyes and (C, D) treatment eyes over time. Bars represent mean percentage of taxa present at ≥ 1% mean relative abundance and totaling 100% at each time point.
Median percentages and ranges of relative bacterial abundance, annotated to the level of family and genus with p-values <0.05 are shown. Phylum data are included for comparison.
Similar to the trends reported in Table 2, LEfSe analysis indicated that Burkholderia and unclassified Actinomycetales were increased on day 35 among control eyes in comparison to day 0 (baseline) and day 7 (Table 3). LEfSe demonstrated additional changes in relative taxa abundance of few bacterial families and genera over time (Table 3).
LDA scores >3.0 are included.
Temporal variability in eyes treated with neomycin-polymyxin-bacitracin
In order to evaluate the temporal stability of the ocular surface microbiome in clinically normal dogs following the application of topical antibiotics, additional conjunctival swabs were collected from treatment eyes subsequent to the baseline samples (day 0). Swabbing was repeated after eyes were treated three times daily for one week (day 7), and four weeks after the completion of antibiotic therapy (day 35).
Species richness and diversity over time.
Alpha diversity metrics were unchanged in treatment eyes between the three sampled time points (S3 Table and Fig 1C).
Microbial community structure over time.
There was no significant change in beta diversity associated with time, evidenced by the lack of clustering in PCoA plots (Fig 5B). The global community structure of treatment eyes did not vary over time based on ANOSIM (weighted UniFrac, R = -0.004, R = 0.11, R = 0.02 for day 0 vs. 7, day 0 vs. 35, and day 7 vs. 35, respectively, p > 0.05); (unweighted UniFrac, R = -0.01, R = -0.02, R = -0.005 for day 0 vs. 7, day 0 vs. 35, and day 7 vs. 35, respectively, p > 0.05). Weighted UniFrac PCoA plots gave us similar results (data not shown). There was no difference in unweighted Unifrac distances over time in treatment eyes, when measured from baseline (day 0) (p > 0.05) (Fig 6). Therefore, bacterial communities on day 7 compared to day 0 were no more different than communities at day 35 compared to day 0. Additionally, Unifrac distances between control eyes and treatment eyes did not vary, suggesting bacterial communities remained similar over time regardless of treatment (p > 0.05) (Fig 6).
Microbial community composition over time.
Fig 7B and 7C illustrates the mean relative abundance of bacteria in treatment eyes at each of the time points sampled. Variance in bacterial abundance was observed between individual dogs (S5 Fig). Friedman and Dunn’s multiple comparison tests revealed the abundance of some bacterial families and genera changed on the ocular surface of treatment eyes over time (q < 0.05) (Table 4). No significant changes were detected at the phylum level (q > 0.05). At the family level, treatment eyes had significantly more Microbacteriaceae on day 0 compared to day 7 and day 35 (p = 0.001, q = 0.029). At the genus level, treatment eyes had significantly more Salinibacterium on day 0 compared to day 7 and day 35 (p = 0.001, q = 0.032). Unclassified Enterobacteriaceae were significantly enriched in treatment eyes on day 35, compared to day 0 and day 7 (p = 0.001, q = 0.040).
Median percentages and ranges of relative bacterial abundance, annotated to the level of family and genus with p-values < 0.05 are shown. Phylum data are included for comparison.
LEfSe analysis exhibited similar temporal changes in relative taxa abundance of some bacterial families and genera (Table 5). Concurrent with Table 4, Microbacteriaceae and Salinibacterium were comparatively increased among treatment eyes on day 0 (baseline) compared to day 7 and day 35, while unclassified Enterobacteriaceae were enriched on day 35 compared to day 0 and day 7 (Table 5 and Fig 7C and 7D).
LDA scores >3.0 are shown.
Discussion
The bacterial community associated with the canine ocular surface is more complex and diverse than previously reported with culture-dependent studies. All 26 eyes sampled contained bacteria from at least 4 phyla and 2 families at >1% relative abundance (S2 Table). Furthermore, the ocular surface microbiome in our population of privately-owned dogs is distinct to the eye when compared to previous reports on the canine skin, mouth, and nasal mucosa [39–41].
The most common bacterial phyla and their relative proportions detected on the canine ocular surface in the present study, Proteobacteria (49.7%), Actinobacteria (25.5%), Firmicutes (12.0%), and Bacteroidetes (7.5%), are analogous to descriptions of the human [17,22,23] and equine [24] ocular surface microbiome. Preliminary studies examining the ocular surface microbiome of cats and dogs using culture-independent methods identified the same bacterial phyla at different proportions, with a relatively increased abundance of Firmicutes across all feline (30–43%) [25,27] and canine samples (34.9%) [28]. It is important to note; however, that limitations exist when comparing microbiome studies as numerous variations in methodologies for DNA extraction, sequencing, analysis, contaminant filtering, and clustering strategies may influence the results.
The most common bacterial families and their relative proportions sequenced from the majority of canine eyes sampled were Pseudomonadaceae (13.2%), Micrococcaceae (12.0%), Pasteurellaceae (6.9%), Microbacteriaceae (5.2%), Enterobacteriaceae (3.9%), Neisseriaceae (3.5%), and Corynebacteriaceae (3.3%) (Fig 4). Gram-positive bacteria regularly cultured from the canine ocular surface, such as Micrococcaceae, Corynebacteriaceae, Staphylococcaceae, Bacillaceae, and Streptococcaceae, represented 12.0%, 3.3%, 2.3%, 1.5% and 1.2% of the bacterial families sequenced, respectively, and were detected in 65–96% of the eyes sampled (S2 Table). A larger proportion of the most abundantly sequenced microorganisms were Gram-negative, contradicting the previous culture-based evidence of a primarly Gram-positive ocular microbiota. This finding may be explained by inherent bias of culture-based detection toward fast growing bacteria that can be easily cultivated in standard laboratory settings [17]; however, additional NGS studies with appropriate sequencing of negative controls and rigourous contaminant filtering are necessary to validate our findings.
This study identified numerous taxa that were never before linked to the canine ocular surface, possibly due to lack of cultivability. Previously unrecognized organisms that may inhabit the canine ocular surface include families from 4 phyla: Proteobacteria (Sphingomonadaceae, Oxalobacteraceae, Unclassified Oxalobacteraceae, Rhodobacteraceae, Methylobacteriaceae, Xanthomonadeceae, Unclassified Xanthomonadeceae, Bradyrhizobiaceae, Unclassified Bradyrhizobiaceae, Burkholderiaceae, Comamonadaceae), Actinobacteria (Unclassified Micrococcaceae, Intrasporangiaceae), Firmicutes (Lachnospiraceae) and Bacteroidetes (Cytophagaceae, Sphingobacteriaceae, Weeksellaceae) (S2 Table). Although NGS has allowed us to detect the presence of these taxa on the canine ocular surface, their role in the microbiome and impact on the health of the canine eye is currently unknown.
In control eyes, the relative abundances of bacterial taxa, along with alpha and beta diversity metrics, did not vary when compared over three separate time points: day 0, day 7, and day 35. This finding supports the notion that species richness, community structure, and community composition of the canine ocular surface microbiome are not significantly altered over time.
The ocular surface, similar to the skin, mouth, and nasal cavity, is an open and exposed system in constant contact with its environment and external microbes. However, the ocular surface is thought to have a relatively low microbial abundance compared to other open systems in the body [18,23], primarily due to innate defense mechanisms that protect against infection such as blinking, tearing and the presence of antimicrobial secretions [42]. One may argue that such mechanisms prevent the development of a stable core bacterial population, only allowing for a haphazard collection of transient organisms to temporarily reside on the ocular surface. Individual variation in the relative abundances of taxa both between eyes and between dogs at baseline were observed at the family level (S4 and S5 Figs). This finding; however, is not unique to the ocular surface, as a high degree of interindividual variability exists within human and animal microbiomes, and is likely attributed to environmental factors and host genetics [23]. A shift in bacterial communities; however, was not measured over time or between treatment and control eyes in our study, suggesting the general bacterial composition of the ocular surface remains relatively unchanged (and is not a haphazard collection) despite this variability (Fig 6). This study, therefore, proposes the presence of both a core and transient microbial population on the canine ocular surface. Similar to the equine [24] and feline [27] eye, several bacteria taxa were present in the majority of canine eyes at each time point sampled, consistent with the presence of a stable core community of microbes. Many of the bacterial families listed in Fig 7B likely constitute the core microbiome of the canine ocular surface; however, additional longitudinal investigations are necessary with more rigorous contaminant filtering to support our findings and define its composition.
In treatment eyes, alpha and beta diversity metrics did not vary when compared between three separate time points: prior to antibiotic therapy (day 0), following one week of treatment with neomycin-polymyxin-bacitracin (day 7), and four weeks after stopping antibiotic therapy (day 35). As noted in control eyes, there were no significant differences in the relative abundance of bacterial phyla in treatment eyes over time. Although the current study detected minor fluctuations in the relative abundance of few bacterial families and genera in treatment eyes over time (Table 4 and Table 5), comparable trends were noted in control eyes. We suspect these minor variations were more likely caused by a combination of host and environmental factors rather than antibiotic treatment. Similar studies of the equine and feline ocular surface also recognized a stable bacterial community that remained consistent following short-term topical broad-spectrum antibiotic use [24,27]. This demonstrates the bacterial microbiome of the ocular surface is not altered by the short term use of broad spectrum topical antibiotics, as species richeness, community structure, and global community composition remained stable.
Neomycin-polymyxin-bacitracin is a broad-spectrum triple antibiotic ointment commonly prescribed in veterinary medicine to prevent or treat bacterial infections of the ocular surface. Neomycin disrupts bacterial protein synthesis, while polymyxin B increases the permeability of bacterial cell membranes [43]. Both antibiotics provide a Gram-negative spectrum of activity [4,43]. Bacitracin inhibits bacterial cell wall synthesis and is primarily active against Gram-positive bacteria [4,43]. This combination of antibiotics is most effective for surface disease, as transcorneal and ocular penetration of the drugs are poor [44]. Additionally, the bioavailability of topically applied medication is generally low [45]. Therefore, systemic absorption of neomycin-polymycin-bacitracin is not expected to occur nor reach therapeutic concentrations on the ocular surface of the contralateral (control) eye. However, in the absence of data showing antibiotic levels in both eyes or longitudinal data in completely treatment-naïve dogs, we cannot completely rule out the possibility of systemic absorption contributing to the similar results in both treatment and control eyes.
There are several reasons neomycin-polymyxin-bacitracin did not shift the ocular surface microbiome as we often see when profiling the gastrointestinal microbiome following parenteral antibiotics. For one, topically applied ophthalmic drugs have a limited retention time on the ocular surface due to lacrimation and drainage. After application to the eye, ointment is retained in the conjunctival fornix until the petrolatum base melts, exposing and delivering the drug to the tear film [46]. Although ointments have increased retention times compared to aqueous solutions and suspensions, the continual turnover of tears will eventually eliminate the medication in 2–3 hours [47]. Therefore, with three times daily application, a common dosage frequency for ophthalmic ointments, the regional microbiome of healthy eyes may be able to rebound and repopulate between dosages. This broad-spectrum combination of antibiotics is frequently used in veterinary medicine as prophylactic therapy to prevent infection by opportunistic bacteria following corneal ulceration, and to treat conjunctival bacterial overgrowth in patients with ocular surface disease such as dry eye [1–7]. While effective in these scenarios, more intensive antibiotic therapy is required with the treatment of bacterial keratitis. Dogs, especially brachycephalic breeds, are susceptible to severe ocular surface diseases that require prolonged antibiotic therapies for several weeks to months [1–7]. Chronic and frequent use of ophthalmic antibiotics of several week duration may have more notable effects on the ocular surface microbiota of diseased eyes, and enable the development of antibiotic resistence [10–13]. Therefore, additional studies are warranted to further evaluate chronic topical antimicrobial use and its impact on the canine ocular surface microbiome.
This study evaluated a relatively small and heterogeneous canine population, representing a diverse sampling of privately-owned dogs with different housing environments. Although homogeneity is recommended when targeting a study population to eliminate compounding variables, our findings are more likely to represent what is found in the general canine population compared to a research colony of dogs. Although two pairs of dogs shared the same household in our study, their samples did not display apparent clustering on principal coordinate analysis plots (S3 Fig), nor did we find clustering by eyes of individual dogs (S2 Fig). Currently, the effect of population variations such as age, sex, season, geography, and other environmental factors on the composition of the ocular surface microbiome of veterinary species is unknown. Future investigations incorporating larger and more geographically diverse canine populations are warranted to limit bias and more accurately describe the variability present on the ocular surface microbiome.
In addition to a relatively small sample size, owner compliance was a potential limitation to this study. All dogs were owned by personnel of the Veterinary Medical Teaching Hospital whose training and medical skills allowed for the appropriate administration of ophthalmic ointment to the eye. While treatment sheets with detailed instructions and a checklist for treatment application times were distributed to the owners and returned to the investigators at the conclusion of the study, compliance was not verified further.
Additional limitations are intrinsic to the interpretation of microbiome studies, including the assessment of relative abundance, which does not take into account the absolute bacterial quantities present in a sample. The microbiome datasets generated from high-throughput sequencing are compositional, yet many microbiome analyses, including those in this study, use non-compositional models where an increase in proportion does not always translate into an increase in absolute abundance of bacteria [48]. Quantitative PCR of known, previously identified organisms is required to detect absolute quantities and should be considered, along with the compositional analyses of microbiome datasets, with future NGS studies. The viability of organisms detected via NGS cannot be determined, possibly leading to an overrepresentation of nonviable organisms controlled by the host’s immune response [18].
Relatively low biomass environments, such the ocular surface, are more susceptible to having contaminating DNA from laboratory reagents negatively affect the results, and some bacteria noted in this study have been reported as common contaminants [49]. Although DNA extraction and sequencing of our sample collection blanks (unused swabs) was performed and did not yield any data, a limitation to this study is the failure to sequence DNA extraction blanks (reagents only with no swab included). The sequencing laboratory in our study reported the likelihood of background noise interfering with the sample data to be remote at best. However, a lack of subtractive filtering to remove contaminants, even when negative sample controls fail to show amplification on PCR or data on sequencing, may still allow potential reagent contaminants to confound the interpretation of microbiome data [50]. This limitation may impact the reliability of our results and other studies that failed to sequence extraction and PCR reagents blanks. Future studies will include sequencing of our negative sample and reagents controls, regardless of PCR amplification, to eliminate any concerns for contamination.
The usage of a topical anesthetic, while necessary to apply greater pressure during swabbing to obtain a more accurate sample of the ocular surface [17,21], may also limit microbial diversity by diluting the ocular surface [51]. Despite these inherent limitations, there remains a wealth of valuable information to improve our insight into the ocular surface microbiome and elucidate its influence on the eye.
Conclusion
This report investigates the temporal stability of the canine ocular surface microbiome in clinically normal eyes with and without topical broad-spectrum antibiotic therapy. In contrast to previous culture-dependent studies, all canine eyes demonstrated the presence of bacteria, several of which were Gram-negative and formerly extraneous to the canine eye. Similar to equine and feline studies using molecular-based techniques, a diverse and stable bacterial community was identified to inhabit the canine ocular surface, and recognized to remain consistent over time and following short-term topical broad-spectrum antibiotic use. Investigations into the canine ocular surface microbiome in diseased eyes are ongoing to reveal if changes to bacterial microbial communites are associated with ocular disease.
Supporting information
S1 Checklist. Example of treatment sheet with instructions and checklist for medication application times distributed to owners for increased compliance.
https://doi.org/10.1371/journal.pone.0234313.s001
(PDF)
S1 Table. Alpha diversity averages for eyes at baseline (day 0) measured at 13,338 sequences per sample.
https://doi.org/10.1371/journal.pone.0234313.s002
(DOCX)
S2 Table. Taxa present at ≥1% mean relative abundance in healthy dogs at baseline (day 0).
Mean percentages and standard deviation of relatively abundant bacteria, annotated to the level of phylum, family, and genus, are represented.
https://doi.org/10.1371/journal.pone.0234313.s003
(DOCX)
S3 Table. Alpha diversity averages for control eyes and treatment eyes over time measured at 13,338 sequences per sample.
https://doi.org/10.1371/journal.pone.0234313.s004
(DOCX)
S1 Fig. Rarefaction analysis of 16S rRNA gene sequences from healthy dogs, comparing treatment and control eyes.
Lines represent the mean of each group for all three time points sampled.
https://doi.org/10.1371/journal.pone.0234313.s005
(TIF)
S2 Fig. Principal coordinate analysis plot (PCoA) of unweighted UniFrac distance matrices of left (OS) and right (OD) eyes from 13 healthy dogs at baseline.
Clustering was not observed between left and right eyes of individual dogs. Based on ANOSIM, pairwise comparisons between dogs were not signifant (Unweighted UniFrac, 999 permutations: mean R = -0.008, median R = -0.008, SD = 0.111, p > 0.05 (ranging from 0.333–1).
https://doi.org/10.1371/journal.pone.0234313.s006
(TIF)
S3 Fig. Principal coordinate analysis plot (PCoA) of unweighted UniFrac distance matrices of both eyes (n = 26) from healthy dogs (n = 13) at baseline.
Clustering was not observed between dogs with shared households. Based on ANOSIM, pairwise comparison between households was not significant (Unweighted Unifrac, R = 0.313, p > 0.05).
https://doi.org/10.1371/journal.pone.0234313.s007
(TIF)
S4 Fig. Bacterial family composition of individual control eyes (n = 13) over time.
Bars represent mean relative abundance of all taxa present in ≥ 6/13 eyes at each time point.
https://doi.org/10.1371/journal.pone.0234313.s008
(TIF)
S5 Fig. Bacterial family composition of individual treatment eyes (n = 13) over time.
Bars represent mean relative abundance of all taxa present in ≥ 6/13 eyes at each time point.
https://doi.org/10.1371/journal.pone.0234313.s009
(TIF)
References
- 1. Murphy JM, Lavach JD, Severin GA. Survey of conjunctival flora in dogs with clinical signs of external eye disease. J Am Vet Med Assoc. 1978; 172(1): 66–8. pmid:624664
- 2. Gerding JP, McLaughlin SA, Troop MW. Pathogenic bacteria and fungi associated with external ocular diseases in dogs: 131 cases (1981–1986). J Am Vet Med Assoc. 1988; 193(2): 242–4. pmid:3403357
- 3. Whitley RD. Canine and feline primary ocular bacterial infections. Vet Clin Small Anim Pract. 2000; 30(5): 1151–67.
- 4. Tolar EL, Hendrix DV, Rohrbach BW, Plummer CE, Brooks DE, Gelatt KN. Evaluation of clinical characteristics and bacterial isolates in dogs with bacterial keratitis: 97 cases (1993–2003). J Am Vet Med Assoc. 2006; 228(1): 80–5. pmid:16426172
- 5. Prado MR, Rocha MF, Brito EH, Girao MD, Monteiro AJ, Teixeira MF, et al. Survey of bacterial microorganisms in the conjunctival sac of clinically normal dogs and dogs with ulcerative keratitis in Fortaleza, Ceará, Brazil. Vet Ophthalmol. 2005; 8(1):33–37. pmid:15644098
- 6. Lin CT, Petersen‐Jones SM. Antibiotic susceptibility of bacterial isolates from corneal ulcers of dogs in Taiwan. J Small Anim Pract 2007; 48(5): 271–274. pmid:17425695
- 7. Wang L, Pan Q, Zhang L, Xue Q, Cui J, Qi C. Investigation of bacterial microorganisms in the conjunctival sac of clinically normal dogs and dogs with ulcerative keratitis in Beijing, China. Vet Ophthalmol. 2008; 11(3):145–9. pmid:18435654
- 8. Gilger BC. Immunology of the ocular surface. Vet Clin North Am Small Anim Pract. 2008; 38(2): 223–231. pmid:18299004
- 9. Kugadas A, Gadjeva M. Impact of microbiome on ocular health. Ocul Surf. 2016; 14(3): 342–349. pmid:27189865
- 10. Suchodolski JS, Dowd SE, Westermarck E, Steiner JM, Wolcott RD, Spillmann T, et al. The effect of the macrolide antibiotic tylosin on microbial diversity in the canine small intestine as demonstrated by massive parallel 16S rRNA gene sequencing. BMC Micro. 2009; 9(1): 210.
- 11. Langdon A, Crook N, Dantas G. The effects of antibiotics on the microbiome throughout development and alternative approaches for therapeutic modulation. Genome Med. 2016; 8(1):39 pmid:27074706
- 12. Sandmeyer LS, Bauer BS, Mohaghegh Poor SM, Feng CX, Chirino-Trejo M. Alterations in conjunctival bacteria and antimicrobial susceptibility during topical administration of ofloxacin after cataract surgery in dogs. Am J Vet Res. 2017; 78(2): 207–14. pmid:28140645
- 13. Varges R, Penna B, Martins G, Martins R, Lelenbaum W. Antimicrobial susceptibility of Staphylococci isolated from naturally occurring canine external ocular diseases. Vet Ophthalmol. 2009; 12(4):216–220. pmid:19604336
- 14. McDonald PJ, Watson AD. Microbial flora of normal canine conjunctivae. J Small Anim Pract. 1976; 17(12): 809–12. pmid:1011806
- 15. Gerding PA, Kakoma I. Microbiology of the canine and feline eye. Vet Clin North Am Small Anim Pract. 1990; 20(3): 615–25. pmid:2194350
- 16. Furiani N, Scarampella F, Anna Martino P, Panzini I, Fabbri E, Ordeix L. Evaluation of the bacterial microflora of the conjunctival sac of healthy dogs and dogs with atopic dermatitis. Vet Dermatol. 2011; 22(6): 490–6. pmid:21535255
- 17. Dong Q, Brulc JM, Iovieno A, Bates B, Garoutte A, Miller D, et al. Diversity of bacteria at healthy human conjunctiva. Invest Ophthalmol Vis Sci. 2011;52(8): 5408–13. pmid:21571682
- 18. Willcox MDP. Characterization of the normal microbiota of the ocular surface. Exp Eye Res. 2013; 117:99–105. pmid:23797046
- 19. Zegans ME, Van Gelder RN. Considerations in understanding the ocular surface microbiome. Am J Ophthalmol. 2014; 158(3):420–2. pmid:25132249
- 20. Doan T, Akileswaran L, Andersen D, Johnson B, Ko N, Shrestha A, et al. Paucibacterial microbiome and resident DNA Virome of the healthy conjunctiva. Invest Ophthalmol Vis Sci. 2016; 57(13):5116–5126. pmid:27699405
- 21. Lu LJ, Liu J. Human microbiota and ophthalmic disease. Yale J Bio Med. 2016; 89(3):325–330.
- 22. Huang Y, Yang B, Li W. Defining the normal core microbiome of conjunctival microbial communities. Clinical Microbiology and Infection. 2016 Jul 1;22(7):643–e7.
- 23. Ozkan J, Nielsen S, Diez-Vives C, Coroneo M, Thomas T, Willcox M. Temporal stability and composition of the ocular surface microbiome. Nature Sci Rep. 2017; pmid:28852195 7(1):9880.
- 24. Scott EM, Arnold C, Dowell S, Suchodolski JS. Evalutation of the bacterial ocular surface microbiome in clinically normal horses before and after treatment with topical neomycin-polymyxin-bacitracin. PLoS ONE. 2019; 14(4): e0214877. pmid:30943258
- 25. Weese SJ, Nichols J, Jalali M, Litster A. The oral and conjunctival microbiotas in cats with and without feline immunodeficiency virus infection. Vet Res. 2015; pmid:25879465
- 26. Weese JS. The canine and feline skin microbiome in health and disease. Vet Dermatol. 2013; 24(1):137–145. pmid:23331690
- 27. Darden JE, Scott EM, Arnold C, Scallan EM, Simon BT, Suchodolski JS. Evaluation of the bacterial ocular surface microbiome in clinically normal cats before and after treatment with topical erythromycin. PLoS ONE. 2019; 14(10): e0223859. pmid:31603921
- 28. Leis ML, Costa MO. Initial description of the core ocular surface microbiome in dogs: bacterial community diversity and composition in a defined canine population. Vet Ophthalmol. 2018; pmid:30095241
- 29. Meason-Smith C, Diesel A, Patterson AP, Older CE, Mansell JM, Suchodolski JS, et al. What is living on your dog’s skin? Characterization of the canine cutaneous mycobiota and fungal dysbiosis in canine allergic dermatitis. FEMS Microbiol Ecol. 2015; 91: 1–12.
- 30. Older CE, Diesel A, Patterson AP, Meason-Smith C, Johnson TJ, Mansell J, et al. The feline skin microbiota: the bacteria inhabiting the skin of healthy and allergic cats. PLoS ONE. 2017; 12(6): e0178555. pmid:28575016
- 31. Meason-Smith C, Diesel A, Patterson AP, Older CE, Johnson TJ, Mansell JM, et al. Characterization of the cutaneous mycobiota in healthy and allergic cats using next generation sequencing. Vet Dermatol. 2017; 28: 71–e17. pmid:27553477
- 32. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, et al. QIIME2: Reproducible, interactive, scalable, extensible microbiome data science. Peer J Preprints. 2018;
- 33. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016; 13(7): 581–583. pmid:27214047
- 34. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006; 72(7): 5069–5072. pmid:16820507
- 35. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. App Environ Microbiol. 2005; 71(12): 8228–8235.
- 36. Benjamini Y, Hochberg Y. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995; 57(1): 289–300.
- 37. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011; 12(6) pmid:21702898
- 38. Zakrzewski M, Proietti C, Ellis JJ, Hasan S, Brion MJ, Berger B, et al. Calypso: a user-friendly webserver for mining and visualizing microbiome-environment interactions. Bioinfo. 2017; 33(5): 782–783.
- 39. Hoffmann R, Patterson AP, Diesel A, Lawhon SD, Ly HJ, Elkins Stephenson C, et al. The skin microbiome in healthy and allergic dogs. PLoS ONE. 2014; 9(1):e83197. pmid:24421875
- 40. Isaiah A, Hoffmann AR, Kelley R, Mundell P, Steiner JM, Suchodolski JS. Characterization of the nasal and oral microbiota of detection dogs. PloS one. 2017 Sep 21;12(9):e0184899. pmid:28934260
- 41. Tress B, Dorn ES, Suchodolski JS, Nisar T, Ravindran P, Weber K, et al. Bacterial microbiome of the nose of healthy dogs and dogs with nasal disease. PloS one. 2017 May 1;12(5):e0176736. pmid:28459886
- 42. McDermott AM. Antimicrobial compounds in tears. Exp Eye Res. 2013; 117:53–61. pmid:23880529
- 43.
Clode A. Antibacterial agents, antifungal agents, and antiviral agents. In: Gelatt KN, Gilger BC, Kern TJ, editors. Veterinary Ophthalmology. 5th ed. Ames: John Wiley & Sons; 2013.
- 44.
Mauger T. Antimicrobials. In: Mauger T, Craig E, editors. Havener’s Ocular Pharmacology. 6th ed. St. Louis: Mosby; 1994.
- 45. Dey S, Gunda S, Mitra AK. Pharmacokinetics of erythryomycin in rabbit corneas after single-dose infusion: role of p-glycoprotein as a barrier to in vivo ocular drug absorption. J Pharmacol Exp Ther. 2004; 311(1): 246–255. pmid:15175422
- 46. Frangie JP. Clinical pharmacokinetics of various topical ophthalmic delivery systems. Clin Pharmacokinet. 1995; 29: 130–138. pmid:7586901
- 47. Sasaki H, Yamamura K, Nishida K, Nakamura J, Ichikawa M. Delivery of drugs to the eye by topical application. Prog Retin Eye Res. 1996; 15: 583–620.
- 48. Gloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ. Microbiome datasets are compostional: and this is not optional. Front Microbiol 2017; 8: https://doi.org/10.3389/fmicb.2017.02224
- 49. Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biology. 2014; 12(87): http://www.biomedcentral.com/1741-7007/12/87
- 50. Eisenhofer R, Minich JJ, Marotz C, Cooper A, Knight R, Weyrich LS. Contamination in low microbioal biomass microbiome studies: issues and recommendations. Trends Microbiol. 2019; 27(2): 105–117. pmid:30497919
- 51. Shin H, Price K, Albert L, Dodick J, Park L, Dominguez-Bello MG. Changes in the eye microbiota associated with contact lens wearing. mBio. 2016; 7(2):e00198–16. pmid:27006462