Advertisement
Original Article Immunology/Infection| Volume 136, ISSUE 6, P1182-1190, June 2016

Longitudinal Evaluation of the Skin Microbiome and Association with Microenvironment and Treatment in Canine Atopic Dermatitis

Open ArchivePublished:February 05, 2016DOI:https://doi.org/10.1016/j.jid.2016.01.023
      Host-microbe interactions may play a fundamental role in the pathogenesis of atopic dermatitis, a chronic relapsing inflammatory skin disorder characterized by universal colonization with Staphylococcus species. To examine the relationship between epidermal barrier function and the cutaneous microbiota in atopic dermatitis, this study used a spontaneous model of canine atopic dermatitis. In a cohort of 14 dogs with canine atopic dermatitis, the skin microbiota were longitudinally evaluated with parallel assessment of skin barrier function at disease flare, during antimicrobial therapy, and post-therapy. Sequencing of the bacterial 16S ribosomal RNA gene showed decreased bacterial diversity and increased proportions of Staphylococcus (S. pseudintermedius in particular) and Corynebacterium species compared with a cohort of healthy control dogs (n = 16). Treatment restored bacterial diversity with decreased proportions of Staphylococcus species, concurrent with decreased canine atopic dermatitis severity. Skin barrier function, as measured by corneometry, pH, and transepidermal water loss also normalized with treatment. Bacterial diversity correlated with transepidermal water loss and pH level but not with corneometry results. These findings provide insights into the relationship between the cutaneous microbiome and skin barrier function in atopic dermatitis, show the impact of antimicrobial therapy on the skin microbiome, and highlight the utility of canine atopic dermatitis as a spontaneous nonrodent model of atopic dermatitis.

      Abbreviations:

      AD (atopic dermatitis), ANOSIM (analysis of similarity), cAD (canine atopic dermatitis), DNA (deoxyribonucleic acid), rRNA (ribosomal ribonucleic acid), TEWL (transepidermal water loss)

      Introduction

      Atopic dermatitis (AD) is a chronic inflammatory skin disorder that affects approximately 10% of children (
      • Spergel J.M.
      From atopic dermatitis to asthma: the atopic march.
      ) and is commonly associated with Staphylococcus aureus colonization (
      • Leyden J.J.
      • Marples R.R.
      • Kligman A.M.
      Staphylococcus aureus in the lesions of atopic dermatitis.
      ). Genetic risk conferred by mutations in the gene encoding the epidermal barrier protein filaggrin suggests that barrier dysfunction in part contributes to the disease (
      • O’Regan G.M.
      • Sandilands A.
      • McLean W.H.I.
      • Irvine A.D.
      Filaggrin in atopic dermatitis.
      ). Environmental factors, including staphylococcal colonization and infection, may also contribute to disease etiology and/or severity. Recent studies have highlighted the dysbiotic nature of the AD skin microbiome, including a predominance of S. aureus during active flares, suggesting a role for S. aureus and the skin microbiome in atopic inflammation (
      • Kong H.H.
      • Oh J.
      • Deming C.
      • Conlan S.
      • Grice E.A.
      • Beatson M.A.
      • et al.
      Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis.
      ).
      Mouse models have highlighted the effects of specific genetic changes in AD, but they are limited in their clinical similarity and do not recapitulate the complexity of the human disease (
      • Scharschmidt T.C.
      • Segre J.A.
      Modeling atopic dermatitis with increasingly complex mouse models.
      ,
      • Marsella R.
      • Girolomoni G.
      Canine models of atopic dermatitis: a useful tool with untapped potential.
      ). Canine atopic dermatitis (cAD) occurs spontaneously and exhibits similar immunological and clinical features of human AD, therefore providing a useful intermediate model (
      • Marsella R.
      • Girolomoni G.
      Canine models of atopic dermatitis: a useful tool with untapped potential.
      ). cAD affects approximately 10% of dogs and presents with similar lesion distribution, life stage of onset, IgE-specific immune responses, and predisposition to chronic and recurrent superficial bacterial dermatitis and folliculitis (
      • Hillier A.
      • Griffin C.E.
      The ACVD task force on canine atopic dermatitis (I): incidence and prevalence.
      ,
      • Santoro D.
      • Marsella R.
      • Pucheu-Haston C.M.
      • Eisenschenk M.N.C.
      • Nuttall T.
      • Bizikova P.
      Review: pathogenesis of canine atopic dermatitis: skin barrier and host-micro-organism interaction.
      ). Epidermal barrier function is impaired in cAD, suggesting a route for epicutaneous sensitization (
      • Santoro D.
      • Marsella R.
      • Pucheu-Haston C.M.
      • Eisenschenk M.N.C.
      • Nuttall T.
      • Bizikova P.
      Review: pathogenesis of canine atopic dermatitis: skin barrier and host-micro-organism interaction.
      ). The pathogenesis of cAD is multifactorial, and complex interactions between genetics and environment are hypothesized, as in human AD (
      • Bizikova P.
      • Pucheu-Haston C.M.
      • Eisenschenk M.N.C.
      • Marsella R.
      • Nuttall T.
      • Santoro D.
      Review: role of genetics and the environment in the pathogenesis of canine atopic dermatitis.
      ).
      Flare states in atopic humans and dogs are associated with colonization and/or superficial infection by Staphylococcus species: S. aureus in humans and S. pseudintermedius or S. schleiferi in dogs (
      • Fazakerley J.
      • Nuttall T.
      • Sales D.
      • Schmidt V.
      • Carter S.D.
      • Hart C.A.
      • et al.
      Staphylococcal colonization of mucosal and lesional skin sites in atopic and healthy dogs.
      ,
      • Furiani N.
      • Scarampella F.
      • Martino P.A.
      • Panzini I.
      • Fabbri E.
      • Ordeix L.
      Evaluation of the bacterial microflora of the conjunctival sac of healthy dogs and dogs with atopic dermatitis.
      ,
      • Kong H.H.
      • Oh J.
      • Deming C.
      • Conlan S.
      • Grice E.A.
      • Beatson M.A.
      • et al.
      Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis.
      ,
      • Leyden J.J.
      • Marples R.R.
      • Kligman A.M.
      Staphylococcus aureus in the lesions of atopic dermatitis.
      ,
      • Santoro D.
      • Marsella R.
      • Pucheu-Haston C.M.
      • Eisenschenk M.N.C.
      • Nuttall T.
      • Bizikova P.
      Review: pathogenesis of canine atopic dermatitis: skin barrier and host-micro-organism interaction.
      ). Recent studies in an Adam17-deficient mouse model suggest that S. aureus drives lesion formation (
      • Kobayashi T.
      • Glatz M.
      • Horiuchi K.
      • Kawasaki H.
      • Akiyama H.
      • Kaplan D.H.
      • et al.
      Dysbiosis and Staphylococcus aureus colonization drives inflammation in atopic dermatitis.
      ). Toxins produced by S. aureus are hypothesized to trigger or exacerbate inflammation in AD (
      • Williams M.R.
      • Gallo R.L.
      The role of the skin microbiome in atopic dermatitis.
      ), such as the δ-toxin recently shown to induce mast cell degranulation and promote inflammatory skin disease (
      • Nakamura Y.
      • Oscherwitz J.
      • Cease K.B.
      • Chan S.M.
      • Muñoz-Planillo R.
      • Hasegawa M.
      • et al.
      Staphylococcus δ-toxin induces allergic skin disease by activating mast cells.
      ). Further understanding of host-microbe dynamics during flare, treatment, and resolution is critical for improved therapies to manage atopic inflammation.
      Herein we report an integrated analysis of the canine cutaneous microbiome and the skin barrier in cAD. The skin microbiome was defined using culture-independent sequencing of the 16S ribosomal RNA (rRNA) gene before, during, and after antimicrobial treatment. In parallel, quantitative assessment of the skin barrier was measured using transepidermal water loss (TEWL), epidermal moisture content (corneometry), and pH. The results of this study should inform future studies of the functional relationship between host cutaneous barrier function and the skin microbiome of humans and dogs.

      Results

      Summary of study participants and design

      Thirty-two dogs (n = 15 affected by cAD, n = 17 unaffected) were enrolled from September 3, 2013 to March 7, 2014 at the University of Pennsylvania Matthew J. Ryan Veterinary Hospital (Table 1, see Supplementary Table S1 online). One dog in each cohort was excluded because of unrelated medical problems. All cAD subjects had active lesions of superficial bacterial dermatitis and folliculitis at enrollment (Figure 1a).
      Table 1Signalment data of study cohorts
      CharacteristiccADControl
      n1416
      Median age, mo (range)71.0 (10–120)79.5 (8–144)
      Ratio male:female5:91:1
      Spayed/neutered1315
      Average cumulative lesion score
       Visit 122.000
       Visit 214.920
       Visit 315.930
      Figure 1
      Figure 1Canine atopic dermatitis. (a) Typical clinical findings of cAD include alopecia and erythema of the periocular region and muzzle as well as evidence of chronic dermatitis and folliculitis in regions such as the pinna, axilla, and groin. (b) The anatomic sites sampled for microbiomic analysis included the mouth, axilla, groin, and concave pinna. (c) Lesion scores (x-axis) measuring cumulative site-specific lesion severity at each study visit (y-axis). Each line represents one subject with cAD (n = 14). ID, identification.
      The skin was swabbed to sample microbiota at anatomic sites with a predilection for cAD lesions: the pinna, axilla, and groin (Figure 1b). The mouth was also sampled, because licking of the skin is a manifestation of pruritus in dogs. Assessment and sampling occurred at three study visits: visit 1 at initial presentation with a flare of cAD and concurrent bacterial dermatitis, visit 2 at the conclusion of 4–6 weeks of culture- and susceptibility-directed oral antimicrobial therapy, and visit 3 at 4–6 weeks after the conclusion of antimicrobial therapy. Healthy dogs were assessed and sampled contemporaneously.
      Site-specific (pinna, axilla, and groin) assessment and semiquantitative scoring of lesion severity were performed based on the parameters used for the Canine Atopic Dermatitis Extent and Severity Index (
      • Olivry T.
      • Marsella R.
      • Iwasaki T.
      • Mueller R.
      International Task Force On Canine Atopic Dermatitis. Validation of CADESI-03, a severity scale for clinical trials enrolling dogs with atopic dermatitis.
      ). This scoring was performed to assess the relationship between microbiota, skin barrier, and clinical signs at a given site. Cumulative site-specific lesion scores varied widely in dogs with cAD (range = 5–45) (Figure 1c). The components of the site-specific lesion scores—erythema, lichenification, and alopecia—were strongly positively correlated with each other (see Supplementary Table S2 online). Cumulative lesion scores decreased with antimicrobial treatment from median (± standard deviation) of 19 (±12.1) at visit 1 to a median of 10 (±13.1) at visit 2, although this decrease was not significant (P = 0.56). Subject 16 was identified as a minor outlier, with a cumulative lesion score of 45 and a severe phenotype of chronic cAD in which 4 weeks of therapy was unlikely to alter the degree of chronic dermatitis and scarring. Upon exclusion of this subject from this particular analysis, the median cumulative lesion score significantly dropped with treatment, from 18 (±7.5) at visit 1 to 9.5 (±8.1) at visit 2 (P = 0.035) (Table 1 and Figure 1c, see Supplementary Table S1). This minor outlier was not removed from any of the subsequent analyses. Five of 14 subjects had an elevated cumulative lesion score at visit 3 compared with visit 2, corresponding to recrudescence of bacterial dermatitis in the posttreatment interval, but this change was not significant across the entire cohort (P = 0.26).

      The microbiome of cAD skin differs from that of normal canine skin

      To analyze microbial communities in cAD and control dogs, the V1–V3 region of the 16S rRNA gene was amplified and sequenced from skin swabs. The skin microbiome of cAD and unaffected control dogs differed at visit 1, before treatment. The Shannon Diversity Index score, an alpha diversity metric that takes into account the number of taxa present and their abundance in the community, was significantly lower in the cAD group compared with the control group (P = 0.0001; Figure 2a, see Supplementary Figure S1 online). These differences were significant at the pinna and axilla but not at the groin or mouth (pinna, P = 5.0 × 10-4; axilla, P = 0.006; groin, P = 0.075; mouth, P = 0.363; Figure 2a) and reiterated with numerous other alpha diversity metrics (see Supplementary Figure S1).
      Figure 2
      Figure 2The microbiome of cAD during flares and prior to treatment. (a) Median Shannon Diversity Index score for each anatomic site sampled in cAD dogs (red) compared with unaffected dogs (blue). Line in box represents median, boxes represent interquartile range, whiskers represent lowest and highest values within 1.5 times the interquartile range, and dots represent outliers. (b) Median taxonomic relative abundance (y-axis) of the 12 genera present in greatest abundance. cAD and control dogs are partitioned by anatomical site sampled (x-axis). (c) Principle coordinates analysis of the weighted UniFrac metric comparing cAD (red) and unaffected controls (blue) in the axilla and pinna. Clustering is significant as assessed by the analysis of similarities test (axilla: R = 0.403, P = 0.001; pinna: R = 0.639, P = 0.001). Percent variability explained by each axis is given. ***P < 0.001, **P < 0.01. cAD, canine atopic dermatitis; PC, principle coordinate.
      The predominant bacteria on healthy canine skin were Porphyromonas, Staphylococcus, Streptococcus, Propionibacterium, and Corynebacterium species and genera belonging to the families Neisseriaceae and Moraxellaceae (Figure 2b, see Supplementary Figure S2 online). Although the same taxa were present in cAD and control dogs, the relative abundance of taxa varied dramatically between the two groups. Dogs with cAD flares had significantly increased relative abundance of Staphylococcus species (cAD median = 45 ± 29%, control median = 5 ± 18%, P < 1.0 × 10-4) across all skin sites. There was a decreased relative abundance of Porphyromonas species in the pinna (cAD median = 1 ± 3%, control median = 10 ± 6%, P = 1.0 × 10-4) and axilla (cAD median = 2 ± 5%, control median = 10 ± 7%, P = 1.0 × 10-4). In the groin, there was a significant median increase in the relative abundance of Corynebacterium species in cAD (cAD median = 9 ± 13%, control median = 1 ± 16%, P = 0.003) (Figure 2b, see Supplementary Figure S2). This trend was similar in the axilla and pinna, although it did not reach significance. In the oral cavity the most abundant taxa did not differ significantly between control and cAD dogs and consisted of many anaerobes including Porphyromonas, Conchiformibius, Fusobacterium, unclassified Moraxellaceae, Flavobacterium, and unclassified Prevotellaceae species (Figure 2b, see Supplementary Figure S2). Skin microbial communities of cAD dogs were significantly different from those of control dogs as determined by the weighted UniFrac metric, a distance metric that takes into account shared phylogeny and is weighted for abundance of observed organisms in the community (R = 0.445, P = 0.001; Analysis of Similarity (ANOSIM) test). Similar differences were observed when examining each anatomical site individually and visualizing clustering using principle coordinates analysis (axilla: R = 0.403, P = 0.001; pinna: R = 0.639, P = 0.001; groin: R = 0.297, P = 0.002; ANOSIM test; Figure 2c). Significant clustering was not associated with patient sex, breed, or age or with antimicrobial class used during treatment (see Supplementary Figure S3 online).

      Treatment normalizes the cAD microbiome

      The Shannon Diversity Index score increased across skin sites during treatment of dogs with cAD between visits 1 and 2 (P = 0.004; Figure 3a) and approached mean Shannon Diversity Index values observed in the control group at visit 2 (control = 7.17 ± 1.44;, cAD = 6.53 ± 1.74; Figure 3a). There was no longer a statistically significant difference in Shannon Diversity Index score between control and cAD dogs at visit 2 (P = 0.090), but significant differences re-emerged at visit 3 during the posttreatment interval (P = 0.008; Figure 3a). These findings suggest that antimicrobial therapy restores diversity of the skin microbiome in cAD but that the effects may dissipate once treatment is withdrawn or if bacterial dermatitis recurs.
      Figure 3
      Figure 3Skin microbiome changes during treatment of cAD. (a) Change in Shannon Diversity Index scores longitudinally with treatment. Line in box represents median, boxes represent interquartile range, whiskers represent lowest and highest values within 1.5 times the interquartile range, and dots represent outliers. (b) Median taxonomic relative abundance. See for the taxonomic legend. ***P < 0.001, **P < 0.01.
      To further examine the relationship between treatment, disease severity, and microbial diversity, correlations between microbial diversity and lesion scores were analyzed. The Shannon Diversity Index score was significantly inversely correlated with site-specific erythema (R = –0.467, P < 1.0 × 10-4), alopecia (R = –0.379, P < 1.0 × 10-4), lichenification (R = –0.454, P < 1.0 × 10-4), and total lesion scoring (R = –0.458, P = 1.0 × 10-4) (see Supplementary Table S2), suggesting that decreased microbial diversity is associated with lesion severity at cAD predilection sites.
      At the taxonomic level, concurrent with increased diversity, relative abundance of Staphylococcus species decreased a median of 45% to 6 ± 16% (P < 0.001) across all skin sites of cAD dogs from visits 1 to 2. Corynebacterium species increased in cAD dogs (median = 8 ± 9%) compared with the control group (4 ± 4%) at visit 3 (P = 0.007) relative to control dogs (Figure 3b, see Supplementary Figure S2). Clustering by the weighted UniFrac metric also decreased with treatment when comparing dogs with cAD to control dogs (R = 0.0502, P = 0.002 at visit 3; R = 0.445, P = 0.001 at visit 1; ANOSIM test). Together, these data indicate that antimicrobial therapy normalizes the skin microbiome of cAD.

      Skin barrier dysfunction correlates with cAD severity and skin microbiome

      In parallel with sampling the skin microbiome, barrier function was assessed by TEWL, hydration (corneometry), and pH (see Supplementary Tables S2 and S3 and Supplementary Figure S4 online). Site-specific lesion scores were positively correlated with TEWL (R = 0.365, P = 0.0001), suggesting that more severe cAD is associated with impaired barrier function, and were negatively correlated with pH (R = –0.194, P = 0.03). Similarly, Shannon Diversity Index score was negatively correlated with TEWL (R = –0.249, P = 0.006) but showed a weakly positive correlation with pH (R = 0.18, P = 0.05). Additional alpha diversity metrics followed similar trends and were also statistically significant (P < 0.05; see Supplementary Table S2). Skin moisture (corneometry) did not correlate with site-specific lesion scores or alpha diversity metrics. These results show that concurrent with cAD flares, alpha diversity decreases and is correlated with disease severity, TEWL, and pH.
      During the treatment of cAD, TEWL and corneometry scores decreased and trended toward the distribution seen in the control group but were not statistically different between visits 1 and 2 (see Supplementary Table S3 and Supplementary Figure S4). The proportional change in lesion scores in patients with cAD correlated with the change in TEWL from visit 1 to visit 2 (R = 0.682, P = 0.02). There was no correlation between changes in lesion scores and TEWL between visits 1 and 3, likely because of recrudescence of bacterial dermatitis in some cAD dogs. The proportional change in corneometry did not correlate with the change in lesion scores from visit 1 to visit 2, but there was a significant correlation between these parameters at visits 1 and 3 (R = 0.624, P = 0.03). Skin pH did not differ significantly between dogs with cAD and healthy controls across the three visits.

      Staphylococcus pseudintermedius predominates in cAD

      The relationships between Staphylococcus species, cAD severity, and treatment in the canine skin microbiome were correlated. Relative abundance of Staphylococcus species was significantly inversely correlated with alpha diversity (R = –0.686, P < 1.0 × 10-4; Figure 4a) and decreased in mean relative abundance with treatment (45%, 15%, and 18% at visits 1, 2, and 3, respectively), in contrast to control dogs in which mean relative abundance of Staphylococcus species remained low across skin sites at each visit (10%; Figure 4b). There was a positive correlation with the relative abundance of Staphylococcus species and site-specific lesion scoring (R = 0.609, P < 1.0 × 10-4; Figure 4c), suggesting that the amount of Staphylococcus species increases as disease severity increases.
      Figure 4
      Figure 4Staphylococcus species and cAD. (a) Correlation between relative abundance of Staphylococcus species (y-axis) and Shannon Diversity Index score (x-axis). (b) Change in Staphylococcus species relative abundance (y-axis) longitudinally at each skin site in cAD (red lines) and control (blue lines) dogs. (c) Correlation between relative abundance of Staphylococcus species (y-axis) and lesion severity (x-axis) as measured by site-specific lesion score. Species-level relative abundance (y-axis) of (d) cultured Staphylococcus isolates and (e) 16S rRNA sequences using phylogenetic placement. cAD, canine atopic dermatitis.
      Figure 4
      Figure 4Staphylococcus species and cAD. (a) Correlation between relative abundance of Staphylococcus species (y-axis) and Shannon Diversity Index score (x-axis). (b) Change in Staphylococcus species relative abundance (y-axis) longitudinally at each skin site in cAD (red lines) and control (blue lines) dogs. (c) Correlation between relative abundance of Staphylococcus species (y-axis) and lesion severity (x-axis) as measured by site-specific lesion score. Species-level relative abundance (y-axis) of (d) cultured Staphylococcus isolates and (e) 16S rRNA sequences using phylogenetic placement. cAD, canine atopic dermatitis.
      Because Staphylococcus species in humans can be commensal (S. epidermidis) or potential pathogens (S. aureus), species identification of Staphylococcus was determined. One microbial culture from a lesional site and axilla in control dogs or when lesions had resolved were performed contemporaneously with microbiome swabs at each visit. Most of the isolates were identified as S. pseudintermedius in both dogs with cAD and healthy controls, followed by S. epidermidis (a coagulase-negative Staphylococcus species) and S. schleiferi (a coagulase-variable Staphylococcus species) (Figure 4d). To speciate and compare 16S rRNA sequence data to cultures, we used a most recent common ancestor analysis of phylogenetic placement using the pplacer algorithm (
      • Matsen F.A.
      • Kodner R.B.
      • Armbrust E.V.
      pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree.
      ) and a curated database of Staphylococcus species genomes (
      • Conlan S.
      • Kong H.H.
      • Segre J.A.
      Species-level analysis of DNA sequence data from the NIH Human Microbiome Project.
      ). Most of the sequences identified across all samples were attributed to S. pseudintermedius (59.4%), with lesser contributions from S. aureus, S. epidermidis, S. haemolyticus, S. hominis, S. lugdunensis, and S. saprophyticus (Figure 4e). S. schleiferi was not detected, suggesting that the region of the 16S rRNA gene used in this study may not be able to reliably classify S. schleiferi to the species level.

      Discussion

      Approximately 70 million dogs live in 40 million households across the United States (
      American Veterinary Medical Association. U.S
      Pet Ownership & Demographics Sourcebook.
      ), and 10% are afflicted with cAD (
      • Hillier A.
      • Griffin C.E.
      The ACVD task force on canine atopic dermatitis (I): incidence and prevalence.
      ). In this spontaneous large animal model of AD, alterations in skin microbiome parallel those observed in AD patients, including increased relative abundance of Staphylococcus spp. and decreased microbial diversity compared to healthy controls. The longitudinal dynamics of the skin microbiome of cAD during flare and treatment correlate with changes in cutaneous barrier function. This work unveils the dynamic relationship between cutaneous barrier function and the skin microbiome in mammalian health and disease states.
      The cutaneous microbiome in cAD has been previously reported in a small cross-sectional study of allergic (n = 6) and healthy dogs (n = 12), with a similar decrease in microbial diversity, but there were no differences in the relative abundance of Staphylococcus organisms between the two cohorts (
      • Rodrigues Hoffmann A.
      • Patterson A.P.
      • Diesel A.
      • Lawhon S.D.
      • Ly H.J.
      • Elkins Stephenson C.
      • et al.
      The skin microbiome in healthy and allergic dogs.
      ). The contrast may be due to different methodologies used, geography, and differing selection criteria (with dogs in the present study presenting with evidence of bacterial dermatitis and disease flare). Our results are harmonious with longitudinal studies in human AD, where microbial diversity is decreased and Staphylococcus and Corynebacterium species predominate during flare states (
      • Kong H.H.
      • Oh J.
      • Deming C.
      • Conlan S.
      • Grice E.A.
      • Beatson M.A.
      • et al.
      Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis.
      ). Furthermore, the Adam17fl/flSox9-Cre murine model of AD is characterized by high relative abundance of cutaneous Staphylococcus and Corynebacterium organisms that normalized after antimicrobial therapy (
      • Kobayashi T.
      • Glatz M.
      • Horiuchi K.
      • Kawasaki H.
      • Akiyama H.
      • Kaplan D.H.
      • et al.
      Dysbiosis and Staphylococcus aureus colonization drives inflammation in atopic dermatitis.
      ).
      In cAD there is a predisposition to the development of coagulase-positive Staphylococcus species colonization and dermatitis as in AD. S. aureus is the primary coagulase-positive Staphylococcus species of human skin and mucosal sites. S. pseudintermedius is a skin and mucosal commensal in the dog and, as we substantiate independent of culture, the most frequent pathogen isolated from dogs with skin or ear canal infections (
      • Bannoehr J.
      • Guardabassi L.
      Staphylococcus pseudintermedius in the dog: taxonomy, diagnostics, ecology, epidemiology and pathogenicity.
      ). Human S. pseudintermedius colonization is rare and primarily restricted to those with regular contact with dogs and cats (
      • Talan D.A.
      • Staatz D.
      • Staatz A.
      • Overturf G.D.
      Frequency of Staphylococcus intermedius as human nasopharyngeal flora.
      • Goodacre R.
      • Harvey R.
      • Howell S.A.
      • Greenham L.W.
      • Noble W.C.
      An epidemiological study of Staphylococcus intermedius strains isolated from dogs, their owners and veterinary surgeons.
      • Guardabassi L.
      • Loeber M.E.
      • Jacobson A.
      Transmission of multiple antimicrobial-resistant Staphylococcus intermedius between dogs affected by deep pyoderma and their owners.
      ). By the same token, S. aureus is infrequently isolated from infection and carriage sites of dogs in clinical practice and in epidemiological surveys, and it is considered a comparatively infrequent canine pathogen (
      • Beck K.M.
      • Waisglass S.E.
      • Dick H.L.N.
      • Weese J.S.
      Prevalence of meticillin-resistant Staphylococcus pseudintermedius (MRSP) from skin and carriage sites of dogs after treatment of their meticillin-resistant or meticillin-sensitive staphylococcal pyoderma.
      ,
      • Morris D.O.
      • Rook K.A.
      • Shofer F.S.
      • Rankin S.C.
      Screening of Staphylococcus aureus, Staphylococcus intermedius, and Staphylococcus schleiferi isolates obtained from small companion animals for antimicrobial resistance: a retrospective review of 749 isolates (2003-04).
      ,
      • Morris D.O.
      • Boston R.C.
      • O’Shea K.
      • Rankin S.C.
      The prevalence of carriage of meticillin-resistant staphylococci by veterinary dermatology practice staff and their respective pets.
      ). The dog may act as a potential vector of S. aureus, which raises zoonotic and anthropozoonotic concerns for potential transfer of pathogens, drug resistance, and genetic elements (
      • Boag A.
      • Loeffler A.
      • Lloyd D.H.
      Methicillin-resistant Staphylococcus aureus isolates from companion animals.
      ,
      • Bramble M.
      • Morris D.
      • Tolomeo P.
      • Lautenbach E.
      Potential role of pet animals in household transmission of methicillin-resistant Staphylococcus aureus: a narrative review.
      ,
      • Misic A.M.
      • Davis M.F.
      • Tyldsley A.S.
      • Hodkinson B.P.
      • Tolomeo P.
      • Hu B.
      • et al.
      The shared microbiota of humans and companion animals as evaluated from Staphylococcus carriage sites.
      ,
      • Song S.J.
      • Lauber C.
      • Costello E.K.
      • Lozupone C.A.
      • Humphrey G.
      • Berg-Lyons D.
      • et al.
      (2013). Cohabiting family members share microbiota with one another and with their dogs.
      ).
      Mechanisms of staphylococcal perturbation of the epidermal barrier are still unclear and may be through direct or indirect (immunostimulatory) means. Mitigation of staphylococcal overgrowth clearly ameliorates disease severity, and the epidermal barrier normalizes. However, methicillin and multidrug resistance are now commonplace in both veterinary and human medicine. The findings presented here may inform future efforts to develop alternative (nonantibiotic) approaches to controlling staphylococcal burden in skin disease.

      Materials and Methods

      Animal subjects

      Dogs with cAD (n = 16) and healthy dogs (n = 14) were prospectively enrolled in this pilot study at the University of Pennsylvania Matthew J. Ryan Veterinary Hospital, after examination by a board-certified veterinary dermatologist (CC, EM, DM). The same dermatologist examined each patient across all time points. Dogs with cAD were included in the study by fulfilling standardized criteria (five criteria of
      • Favrot C.
      • Steffan J.
      • Seewald W.
      • Picco F.
      A prospective study on the clinical features of chronic canine atopic dermatitis and its diagnosis.
      ) and by ruling out dermatoses with similar presentations as depicted by
      • Hensel P.
      • Santoro D.
      • Favrot C.
      • Hill P.
      • Griffin C.
      Canine atopic dermatitis: detailed guidelines for diagnosis and allergen identification.
      based on dermatologic examination and clinical histories. Dermatologic examination in all patients involved, but was not limited to, skin and otic cytology, flea combing/direct examination, and skin scraping. Dogs with evidence of other underlying systemic disease, an atopic history that was not clearly documented, or evidence of active ecto-parasitic (including Demodex species, Sarcoptes species, and Ctenocephalides felis) infestations were excluded. All subjects were on a strict flea control regimen. Four dogs with cAD had a history of antibiotic exposure within 45 days before enrollment. Therapy was prescribed as indicated by the attending veterinary dermatologist. Systemic antimicrobial therapy was prescribed based on aerobic culture and sensitivity at enrollment. See Supplementary Table S4a and b for specific therapies prescribed for each patient. Healthy dogs were enrolled in the study during the same time period. A subset of dogs was owned by veterinarians and veterinary technicians. Informed consent was obtained from all owners before enrollment. All experiments were carried out according to approved Institutional Animal Care and Use Committee protocols. Site-specific assessment and semiquantitative scoring of lesion severity were performed and included erythema, lichenification, and self-induced alopecia, each on a 0–5 scale (0 = no lesions to 5 = most severe). The same dermatologist (CC and DM) assessed/scored each patient at each time point. This scoring system is based on Canine Atopic Dermatitis Extent and Severity Index (
      • Olivry T.
      • Marsella R.
      • Iwasaki T.
      • Mueller R.
      International Task Force On Canine Atopic Dermatitis. Validation of CADESI-03, a severity scale for clinical trials enrolling dogs with atopic dermatitis.
      ), because of familiarity of the dermatologists, akin to the Scoring Atopic Dermatitis (SCORAD) metric used in human dermatology. The lesion score was used to assess for changes at a given site that might correlate with microbial shifts or skin barrier function and not as a means for documenting changes in disease flare.

      Skin microenvironmental assessment

      Noninvasive probes for measuring subsurface-epidermal water content by the capacitance method (Corneometer CM 825, Courage+Khazaka, Cologne, Germany), skin pH (Skin-pH-Meter, PH 905, Courage+Khazaka, Cologne, Germany), and TEWL by diffusion (Tewameter TM300, Courage+Khazaka, Cologne, Germany) were used at each visit according to the manufacturer’s instructions. Sampling sites included the axilla, groin, and concave aspect of the pinna. These sites were chosen because they are sparsely haired, allowing for corneometry and tewametry measurements, and are cAD predilection sites. Dogs were sampled in lateral recumbency with minimal physical restraint. The order of measurements at each site was corneometry, TEWL, and pH. TEWL measurements were averaged at 1-second intervals for a 30-second period. Five corneometry measurements were performed over a 3cm2 area and were averaged. Three consecutive pH readings were averaged. Measurements were performed in the dermatology clinic treatment area. During sampling and assessment, the mean indoor ambient temperature was 22.2°C (range = 18.3–23.9°C) and the mean relative humidity was 25.5% (range, = 16–68%).

      Microbiomic sampling

      The oral cavity, axilla, concave pinna, and groin were swabbed (3cm2 region) using Catch-All Sample Collection Swabs (Epicentre Biotechnologies, Madison, WI). Swabs were rubbed vigorously over the skin site or mouth for 10–15 second intervals. Swabs exposed to air in the treatment room and laboratory at extraction were used as negative control samples. Swabs were placed in 300 μl of Yeast Cell Lysis Solution (Epicentre Biotechnologies, Madison, WI), and the tip of the swab was aseptically cut from the handle and stored at –80°C until extraction.
      Contemporaneous swabs were taken from a single lesional site of dogs with dermatitis and submitted for aerobic bacterial culture. Swabs were processed using standard laboratory protocols, and isolates were identified as Staphylococcus species by use of a conventional biochemical identification system (MicroScan Walkaway 40 PC20 Gram-positive combo-panel, Dade Behring, Sacramento, CA) as described by the manufacturer. The results of the culture and antimicrobial sensitivity testing were used to direct 4–6 weeks of systemic (oral) antimicrobial therapy (see Supplementary Table S4a).

      DNA extraction and sequencing

      DNA extractions were performed as previously described (
      • Misic A.M.
      • Davis M.F.
      • Tyldsley A.S.
      • Hodkinson B.P.
      • Tolomeo P.
      • Hu B.
      • et al.
      The shared microbiota of humans and companion animals as evaluated from Staphylococcus carriage sites.
      ). The V1–V3 hypervariable region of the 16S rRNA gene was amplified using a barcoding strategy as described (
      • Fadrosh D.W.
      • Ma B.
      • Gajer P.
      • Sengamalay N.
      • Ott S.
      • Brotman R.M.
      • et al.
      An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform.
      ) and primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 534R (5′-ATTACCGCGGCTGCTGG-3′). Sequencing was performed on the Illumina MiSeq instrument (Illumina, San Diego, CA) using 300–base pair paired-end chemistry at the University of Maryland Institute for Genome Sciences.

      16S rRNA gene analysis

      Paired-end reads were demultiplexed by Flexbar (
      • Dodt M.
      • Roehr J.
      • Ahmed R.
      • Dieterich C.
      FLEXBAR—Flexible barcode and adapter processing for next-generation sequencing platforms.
      ) and assembled using Pear (
      • Zhang J.
      • Kobert K.
      • Flouri T.
      • Stamatakis A.
      PEAR: a fast and accurate Illumina Paired-End reAd mergeR.
      ), resulting in 17,065,344 sequences. Sequences less than 465 base pairs and greater than 535 base pairs and sequences with 10 or more homopolymers were removed, resulting in 13,023,611 sequences. QIIME version 1.8 (
      • Caporaso J.G.
      • Kuczynski J.
      • Stombaugh J.
      • Bittinger K.
      • Bushman F.D.
      • Costello E.K.
      • et al.
      QIIME allows analysis of high-throughput community sequencing data.
      ) was used for further downstream processing and analyses. Sequences were aligned to the Greengenes database, and taxonomy was assigned using the RDP classifier (
      • Wang Q.
      • Garrity G.M.
      • Tiedje J.M.
      • Cole J.R.
      Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.
      ). Operational taxonomic units were picked using 97% sequence similarity with cd-hit (
      • Li W.
      • Godzik A.
      Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences.
      ,
      • Fu L.
      • Niu B.
      • Zhu Z.
      • Wu S.
      • Li W.
      CD-HIT: accelerated for clustering the next-generation sequencing data.
      ), and a representative sequence based on the most abundant sequence was used. Each sample was rarified to 4,000 sequences for alpha- and beta-diversity analyses. The pplacer binary (
      • Matsen F.A.
      • Kodner R.B.
      • Armbrust E.V.
      pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree.
      ) was used as previously described (
      • Gardner S.E.
      • Hillis S.L.
      • Heilmann K.
      • Segre J.A.
      • Grice E.A.
      The neuropathic diabetic foot ulcer microbiome is associated with clinical factors.
      ) for species level assessment of Staphylococcus organisms using a Staphylococcus reference database (
      • Conlan S.
      • Kong H.H.
      • Segre J.A.
      Species-level analysis of DNA sequence data from the NIH Human Microbiome Project.
      ). Water and processed blank samples were sequenced and processed in parallel, and contaminants were identified and removed as previously reported (
      • Misic A.M.
      • Davis M.F.
      • Tyldsley A.S.
      • Hodkinson B.P.
      • Tolomeo P.
      • Hu B.
      • et al.
      The shared microbiota of humans and companion animals as evaluated from Staphylococcus carriage sites.
      ). Published best practices were used as guidelines (
      • Bokulich N.A.
      • Subramanian S.
      • Faith J.J.
      • Gevers D.
      • Gordon J.I.
      • Knight R.
      • et al.
      Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing.
      ). Sequences were deposited in the National Center for Biotechnology Information (NCBI) Short Read Archive under BioProject Accession No. PRJNA302288.

      Statistics

      The R Statistical Package (
      R Core Team
      R: A language and environment for statistical computing. R Foundation for Statistical Computing. R Foundation for Statistical Computing.
      ) was used for all computations. Nonparametric Wilcoxon rank sum tests were used to compare differences between groups. For within-subject comparisons, paired Wilcoxon signed rank tests were used. Pearson product moment correlation coefficients were calculated for correlations and tested using the Student t test.

      Conflict of Interest

      The authors state no conflict of interest.

      Acknowledgments

      The authors thank Dr. Darcie Kunder, Dr. Fiona Lee, Ms. Colleen Walters, and Mr. Joseph Rogosky for exemplary patient care and assessment and members of the Grice laboratory for their underlying contributions. Research reported in this publication was supported by a grant from the University of Pennsylvania, School of Veterinary Medicine, Center for Host-Microbial Interactions and by the National Institute of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of Health under award numbers R00-AR060873 and R01-AR066663 to EAG. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

      Supplementary Material

      References

        • American Veterinary Medical Association. U.S
        Pet Ownership & Demographics Sourcebook.
        American Veterinary Medical Association, Schaumburg, IL2012
        • Bannoehr J.
        • Guardabassi L.
        Staphylococcus pseudintermedius in the dog: taxonomy, diagnostics, ecology, epidemiology and pathogenicity.
        Vet Dermatol. 2012; 23 (253–e52)
        • Beck K.M.
        • Waisglass S.E.
        • Dick H.L.N.
        • Weese J.S.
        Prevalence of meticillin-resistant Staphylococcus pseudintermedius (MRSP) from skin and carriage sites of dogs after treatment of their meticillin-resistant or meticillin-sensitive staphylococcal pyoderma.
        Vet Dermatol. 2012; 23 (e66–7): 369-375
        • Bizikova P.
        • Pucheu-Haston C.M.
        • Eisenschenk M.N.C.
        • Marsella R.
        • Nuttall T.
        • Santoro D.
        Review: role of genetics and the environment in the pathogenesis of canine atopic dermatitis.
        Vet Dermatol. 2015; 26 (95–e26)
        • Boag A.
        • Loeffler A.
        • Lloyd D.H.
        Methicillin-resistant Staphylococcus aureus isolates from companion animals.
        Vet Rec. 2004; 154: 411
        • Bokulich N.A.
        • Subramanian S.
        • Faith J.J.
        • Gevers D.
        • Gordon J.I.
        • Knight R.
        • et al.
        Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing.
        Nat Methods. 2012; 10: 57-59
        • Bramble M.
        • Morris D.
        • Tolomeo P.
        • Lautenbach E.
        Potential role of pet animals in household transmission of methicillin-resistant Staphylococcus aureus: a narrative review.
        Vector-Borne Zoonotic Dis. 2011; 11: 617-620
        • Caporaso J.G.
        • Kuczynski J.
        • Stombaugh J.
        • Bittinger K.
        • Bushman F.D.
        • Costello E.K.
        • et al.
        QIIME allows analysis of high-throughput community sequencing data.
        Nat Methods. 2010; 7: 335-336
        • Conlan S.
        • Kong H.H.
        • Segre J.A.
        Species-level analysis of DNA sequence data from the NIH Human Microbiome Project.
        PloS One. 2012; 7: e47075
        • Dodt M.
        • Roehr J.
        • Ahmed R.
        • Dieterich C.
        FLEXBAR—Flexible barcode and adapter processing for next-generation sequencing platforms.
        Biology. 2012; 1: 895-905
        • Fadrosh D.W.
        • Ma B.
        • Gajer P.
        • Sengamalay N.
        • Ott S.
        • Brotman R.M.
        • et al.
        An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform.
        Microbiome. 2014; 2: 6
        • Favrot C.
        • Steffan J.
        • Seewald W.
        • Picco F.
        A prospective study on the clinical features of chronic canine atopic dermatitis and its diagnosis.
        Vet Dermatol. 2010; 21: 23-31
        • Fazakerley J.
        • Nuttall T.
        • Sales D.
        • Schmidt V.
        • Carter S.D.
        • Hart C.A.
        • et al.
        Staphylococcal colonization of mucosal and lesional skin sites in atopic and healthy dogs.
        Vet Dermatol. 2009; 20: 179-184
        • Fu L.
        • Niu B.
        • Zhu Z.
        • Wu S.
        • Li W.
        CD-HIT: accelerated for clustering the next-generation sequencing data.
        Bioinformatics. 2012; 28: 3150-3152
        • Furiani N.
        • Scarampella F.
        • Martino P.A.
        • 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: 490-496
        • Gardner S.E.
        • Hillis S.L.
        • Heilmann K.
        • Segre J.A.
        • Grice E.A.
        The neuropathic diabetic foot ulcer microbiome is associated with clinical factors.
        Diabetes. 2013; 62: 923-930
        • Goodacre R.
        • Harvey R.
        • Howell S.A.
        • Greenham L.W.
        • Noble W.C.
        An epidemiological study of Staphylococcus intermedius strains isolated from dogs, their owners and veterinary surgeons.
        J Analytical and Applied Pyrolysis. 1997; 44: 49-64
        • Guardabassi L.
        • Loeber M.E.
        • Jacobson A.
        Transmission of multiple antimicrobial-resistant Staphylococcus intermedius between dogs affected by deep pyoderma and their owners.
        Vet Microbiol. 2004; 98: 23-27
        • Hensel P.
        • Santoro D.
        • Favrot C.
        • Hill P.
        • Griffin C.
        Canine atopic dermatitis: detailed guidelines for diagnosis and allergen identification.
        BMC Vet Res. 2015; 11: 196
        • Hillier A.
        • Griffin C.E.
        The ACVD task force on canine atopic dermatitis (I): incidence and prevalence.
        Vet Immunol Immunopathol. 2001; 81: 147-151
        • Kobayashi T.
        • Glatz M.
        • Horiuchi K.
        • Kawasaki H.
        • Akiyama H.
        • Kaplan D.H.
        • et al.
        Dysbiosis and Staphylococcus aureus colonization drives inflammation in atopic dermatitis.
        Immunity. 2015; 42: 756-766
        • Kong H.H.
        • Oh J.
        • Deming C.
        • Conlan S.
        • Grice E.A.
        • Beatson M.A.
        • et al.
        Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis.
        Genome Res. 2012; 22: 850-859
        • Leyden J.J.
        • Marples R.R.
        • Kligman A.M.
        Staphylococcus aureus in the lesions of atopic dermatitis.
        Br J Dermatol. 1974; 90: 525-530
        • Li W.
        • Godzik A.
        Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences.
        Bioinforma Oxf Engl. 2006; 22: 1658-1659
        • Marsella R.
        • Girolomoni G.
        Canine models of atopic dermatitis: a useful tool with untapped potential.
        J Invest Dermatol. 2009; 129: 2351-2357
        • Matsen F.A.
        • Kodner R.B.
        • Armbrust E.V.
        pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree.
        BMC Bioinformatics. 2010; 11: 538
        • Misic A.M.
        • Davis M.F.
        • Tyldsley A.S.
        • Hodkinson B.P.
        • Tolomeo P.
        • Hu B.
        • et al.
        The shared microbiota of humans and companion animals as evaluated from Staphylococcus carriage sites.
        Microbiome. 2015; 3: 2
        • Morris D.O.
        • Rook K.A.
        • Shofer F.S.
        • Rankin S.C.
        Screening of Staphylococcus aureus, Staphylococcus intermedius, and Staphylococcus schleiferi isolates obtained from small companion animals for antimicrobial resistance: a retrospective review of 749 isolates (2003-04).
        Vet Dermatol. 2006; 17: 332-337
        • Morris D.O.
        • Boston R.C.
        • O’Shea K.
        • Rankin S.C.
        The prevalence of carriage of meticillin-resistant staphylococci by veterinary dermatology practice staff and their respective pets.
        Vet Dermatol. 2010; 21: 400-407
        • Nakamura Y.
        • Oscherwitz J.
        • Cease K.B.
        • Chan S.M.
        • Muñoz-Planillo R.
        • Hasegawa M.
        • et al.
        Staphylococcus δ-toxin induces allergic skin disease by activating mast cells.
        Nature. 2013; 503: 397-401
        • Olivry T.
        • Marsella R.
        • Iwasaki T.
        • Mueller R.
        International Task Force On Canine Atopic Dermatitis. Validation of CADESI-03, a severity scale for clinical trials enrolling dogs with atopic dermatitis.
        Vet Dermatol. 2007; 18: 78-86
        • O’Regan G.M.
        • Sandilands A.
        • McLean W.H.I.
        • Irvine A.D.
        Filaggrin in atopic dermatitis.
        J Allergy Clin Immunol. 2009; 124: R2-R6
        • R Core Team
        R: A language and environment for statistical computing. R Foundation for Statistical Computing. R Foundation for Statistical Computing.
        Austria, Vienna2015
        • Rodrigues Hoffmann A.
        • Patterson A.P.
        • Diesel A.
        • Lawhon S.D.
        • Ly H.J.
        • Elkins Stephenson C.
        • et al.
        The skin microbiome in healthy and allergic dogs.
        PloS One. 2014; 9: e83197
        • Santoro D.
        • Marsella R.
        • Pucheu-Haston C.M.
        • Eisenschenk M.N.C.
        • Nuttall T.
        • Bizikova P.
        Review: pathogenesis of canine atopic dermatitis: skin barrier and host-micro-organism interaction.
        Vet Dermatol. 2015; 26 (84–e25)
        • Scharschmidt T.C.
        • Segre J.A.
        Modeling atopic dermatitis with increasingly complex mouse models.
        J Invest Dermatol. 2008; 128: 1061-1064
        • Song S.J.
        • Lauber C.
        • Costello E.K.
        • Lozupone C.A.
        • Humphrey G.
        • Berg-Lyons D.
        • et al.
        (2013). Cohabiting family members share microbiota with one another and with their dogs.
        eLife. 2013; 2: e00458
        • Spergel J.M.
        From atopic dermatitis to asthma: the atopic march.
        Ann Allergy Asthma Immunol. 2010; 105: 99-106
        • Talan D.A.
        • Staatz D.
        • Staatz A.
        • Overturf G.D.
        Frequency of Staphylococcus intermedius as human nasopharyngeal flora.
        J Clin Microbiol. 1989; 27: 2393
        • Wang Q.
        • Garrity G.M.
        • Tiedje J.M.
        • Cole J.R.
        Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.
        Appl Environ Microbiol. 2007; 73: 5261-5267
        • Williams M.R.
        • Gallo R.L.
        The role of the skin microbiome in atopic dermatitis.
        Curr Allergy Asthma Rep. 2015; 15: 1-10
        • Zhang J.
        • Kobert K.
        • Flouri T.
        • Stamatakis A.
        PEAR: a fast and accurate Illumina Paired-End reAd mergeR.
        Bioinformatics. 2014; 30: 614-620

      Linked Article

      • Canine and Human Atopic Dermatitis: Two Faces of the Same Host-Microbe Interaction
        Journal of Investigative DermatologyVol. 136Issue 6
        • Preview
          Host-microbe interaction has been suggested to play a critical role in the pathogenesis of atopic dermatitis. The dog has been shown to be the best model to study both pathogenesis and microbiome modifications in atopic dermatitis. Bradley et al. show a significant correlation between microbiome diversity, clinical signs, and skin barrier function in atopic dogs before, during, and after antimicrobial therapy.
        • Full-Text
        • PDF
        Open Archive