- Open Access
Alteration of the cutaneous microbiome in psoriasis and potential role in Th17 polarization
Microbiomevolume 6, Article number: 154 (2018)
Psoriasis impacts 1–3% of the world’s population and is characterized by hyper-proliferation of keratinocytes and increased inflammation. At the molecular level, psoriasis is commonly driven by a Th17 response, which serves as a major therapeutic target. Microbiome perturbations have been associated with several immune-mediated diseases such as atopic dermatitis, asthma, and multiple sclerosis. Although a few studies have investigated the association between the skin microbiome and psoriasis, conflicting results have been reported plausibly due to the lack of standardized sampling and profiling protocols, or to inherent microbial variability across human subjects and underpowered studies. To better understand the link between the cutaneous microbiota and psoriasis, we conducted an analysis of skin bacterial communities of 28 psoriasis patients and 26 healthy subjects, sampled at six body sites using a standardized protocol and higher sequencing depth compared to previous studies. Mouse studies were employed to examine dermal microbial-immune interactions of bacterial species identified from our study.
Skin microbiome profiling based on sequencing the 16S rRNA V1–V3 variable region revealed significant differences between the psoriasis-associated and healthy skin microbiota. Comparing the overall community structures, psoriasis-associated microbiota displayed higher diversity and more heterogeneity compared to healthy skin bacterial communities. Specific microbial signatures were associated with psoriatic lesional, psoriatic non-lesional, and healthy skin. Specifically, relative enrichment of Staphylococcus aureus was strongly associated with both lesional and non-lesional psoriatic skin. In contrast, Staphylococcus epidermidis and Propionibacterium acnes were underrepresented in psoriatic lesions compared to healthy skin, especially on the arm, gluteal fold, and trunk. Employing a mouse model to further study the impact of cutaneous Staphylcoccus species on the skin T cell differentiation, we found that newborn mice colonized with Staphylococcus aureus demonstrated strong Th17 polarization, whereas mice colonized with Staphylococcus epidermidis or un-colonized controls showed no such response.
Our results suggest that microbial communities on psoriatic skin is substantially different from those on healthy skin. The psoriatic skin microbiome has increased diversity and reduced stability compared to the healthy skin microbiome. The loss of community stability and decrease in immunoregulatory bacteria such as Staphylococcus epidermidis and Propionibacterium acnes may lead to higher colonization with pathogens such as Staphylococcus aureus, which could exacerbate cutaneous inflammation along the Th17 axis.
Psoriasis is an immune-mediated inflammatory skin disease that impacts 1–3% of the world’s population. The pathogenesis of psoriasis is multifactorial with notable contributions from patient genetics and environmental factors such as lifestyle, diet, and health history [1, 2]. Psoriasis can be mediated by an overactive Th17 response leading to skin inflammation and hyper-proliferation of keratinocytes . In the clinic, blocking components of the Th17 pathway effectively dampens the aberrant immune response in psoriasis patients and controls symptoms, but these treatments do are not curative and disease management effectiveness varies across patients. This highlights the need to further understand the pathogenesis of psoriasis and the factors associated with disease initiation and progression.
The skin is the human body’s largest organ which serves not only as a physical protective barrier against environmental insults, but also as a dynamic interface for host dermal-microbial interactions. The microbial community that inhabits the human skin is highly complex and consists of highly diverse microorganisms including bacteria, fungi, viruses, and archaea . Bacteria have been shown to be essential for skin health by restricting pathogen colonization and fine-tuning resident T cell function [5, 6]. As a result, perturbations to the skin microbial community have the potential to contribute to altered skin immune function. Indeed, dysbiosis of the skin microbiome has been linked to several inflammatory and autoimmune diseases including atopic dermatitis and vitiligo [7, 8], suggesting the importance of the cutaneous microbiome in the health of the skin.
Interestingly, throat and nasal Streptococcal infection have been shown to trigger initiation and exacerbation of psoriasis [9, 10], suggesting a microbial contribution to disease. Moreover, keratinocytes, the most prominent cell type in the epidermis, can trigger innate and adaptive immune responses in psoriasis through interactions with skin bacteria . To date, several studies have sought to characterize the psoriasis-associated skin microbiome and identify bacterial species that might contribute to the pathogenesis of psoriasis [12,13,14,15,16]. However, these studies revealed a lack of consensus on psoriasis-associated microbial signatures plausibly due to the inherent heterogeneity of microbial species that promote immune dysfunction in psoriatic patients and or to different study designs. For example, collecting samples using skin swabs [12, 14] or skin biopsies  introduces significant variability since these methods sample different cutaneous anatomical compartments with likely different associated bacteria . Moreover, these studies used different 16S rRNA primers amplifying different variable regions of the 16S rRNA gene, which may contribute to variance across studies, making cross study comparisons difficult. Therefore, application of a standardized protocol to allow for a better understanding in the relationship between microbiome and disease is critical [17, 18].
In this study, we surveyed the skin microbiome from 28 psoriasis patients and 26 healthy subjects using the standardized protocol recommended by the NIH Human Microbiome Project [19,20,21]. In contrast to some previous studies targeting the V4 region of the 16S rRNA gene , we profiled the skin microbial community using primers targeting the V1–V3 region, which results in more accurate bacterial identities of the skin microbiome at the genus and species levels compared to the traditional V4 approach [20, 22]. We also used higher sequencing depth to ensure high-quality data. Our data revealed significant alterations in the psoriasis skin microbiome and identified Staphylococcus aureus as a potential contributor to psoriasis pathogenesis.
Cohort of patients and skin sampling
The cohort in this study consisted of 28 patients with plaque psoriasis and 26 healthy individuals. To avoid any confounding demographic effects, gender and age composition were matched between the two groups (Table 1). All psoriasis patients were clinically diagnosed with psoriasis at the UCSF Psoriasis and Skin Treatment Center. The psoriasis patients in this study had a mean Psoriasis Area and Severity Index (PASI) of 11.1 representing moderate-to-severe disease. To avoid the variabilities introduced by treatments, we excluded subjects with recent antibiotic treatment and/or other biologic and systemic therapy. In addition, all subjects required to undergo a 2-week wash-out period for topical therapy. Different anatomic sites in the human skin can be categorized into three major groups: dry, moist, and sebaceous. The biogeographical differences across different skin sites provide different environments that support distinct microbial communities [23,24,25,26]. In order to gain a comprehensive view of the psoriasis-associated skin microbiome, we sampled the microbiome across six different skin sites: scalp, trunk, arm, leg, axilla, and gluteal fold, which covers all three skin groups (Table 2). Three different “disease states” were sampled for each skin site: healthy skin from healthy subjects (Healthy), unaffected or non-lesional skin from psoriasis patients (PSO_N), and lesional skin from psoriasis patients (PSO_L). We sampled all six sites for both healthy (Healthy) and unaffected skin (PSO_N). Only sites with psoriasis lesions present were sampled for psoriatic lesional samples (PSO_L). The psoriasis subjects in our cohort most frequently had psoriatic plaques on the arms, legs, and scalp, whereas there was lowest frequency in the axilla (armpit). Intermediate frequency of plaques was found on the trunk and gluteal fold (Table 2). The sampling of these six skin sites from psoriatic lesional skin, psoriatic non-lesional skin, and healthy control skin allowed for an examination of how the psoriatic microbiome differs at different sites as well as how it changes with disease progression (lesional vs non-lesional).
Alteration in psoriatic skin microbiome diversity is site specific and exhibits an increasing trend in alpha diversity and greater heterogeneity compared with healthy skin
The diversity of the microbial community in a given human body site reflects the structure and composition of the community. Alterations in human microbiome diversity has been linked to disease states. For example, reduced bacterial alpha diversity in the gut microbiome has been linked to obesity and inflammatory bowel disease (IBD) [27,28,29] while increased diversity in the vaginal microbiome is associated with bacterial vaginosis [30, 31]. To understand if diversity of the skin microbial community is altered in psoriasis patients, we first examined alpha diversity of psoriatic lesional skin, psoriatic unaffected skin and healthy control skin using four different metrics to measure community richness (chao1 and observed OTUs), evenness (Simpson diversity index) and overall diversity (Shannon index) (Fig. 1a–d, and Table 3a). Overall, we observed increasing diversity in all four measures going from healthy skin to non-lesional skin to lesional skin, with a statistically significant trend for the Simpson (Fig. 1c, p-value = 0.005) and Shannon indices (Fig. 1d, p-value = 0.005). This unidirectional trend in microbiome diversity suggests that the skin microbiome community diversifies as psoriatic disease progresses. To evaluate alpha diversity at different skin sites, we further examined the four metrics at each skin site. Interestingly, we found significantly increased community richness (chao1) in scalp psoriatic lesions compared to healthy scalp and increased community evenness (Simpson and Shannon indices) in arm psoriatic lesional and non-lesional skin compared to arm healthy control skin, with a significant trend test in the arm for the Simpson and Shannon indices (Table 4). When we grouped samples by skin type (Table 5), we observed higher alpha diversity in all four indices at dry psoriatic skin sites (arm, leg, trunk combined) relative to healthy skin, but no difference in alpha diversity for moist sites (axilla, gluteal fold combined). Overall, these results indicate that increased alpha diversity in psoriasis is mostly observed at dry skin sites, with a trend at the sebaceous (scalp) site, and no increase at moist sites. Our data demonstrates that the association between skin microbiome and psoriasis is complex and sometimes site and/or skin type specific. This highlights the need for comprehensive sampling at various skin sites and skin types to study the skin microbiome in association with cutaneous disease.
We further explored the relationship among bacterial communities isolated from psoriatic and healthy skin by calculating beta diversity using weighted Unifrac distance . There was no distinct difference between bacterial communities isolated from the healthy skin and psoriatic skin as there was not a distinctive separation between bacterial communities isolated according to skin status (Fig. 1c) and the first PC is not significantly different in both psoriasis disease states and healthy skin (PSO_L vs Healthy: p value = 0.109, PSO_N vs. Healthy: p value = 0.128). Although we did not observe distinct clusters associated with disease states, the bacterial communities isolated from psoriatic skin were more dispersed in the principal coordinate analysis than those from healthy skin (Fig. 1e). Indeed, we assess the community dispersion of each disease status by calculating the mean weighted Unifrac distance between bacterial microbiota found in either healthy, psoriatic non-lesional or psoriatic lesional skin and noted that psoriatic non-lesional skin or psoriatic lesional skin exhibited significantly higher mean distances compared with healthy skin (Fig. 1f), indicating greater heterogeneity in the composition of skin microbiota of psoriatic patients irrespective of lesions. We observed a similar trend of increasing heterogeneity by disease state in skin bacterial communities isolated from the arm, trunk, and leg (Additional file 1: Figure S1A, S1B, S1C) as well as in the dry skin group (Additional file 1: Figure S2A). In the moist skin group, bacterial communities of psoriasis lesional skin also exhibited higher heterogeneity compared to healthy and non-lesional skin (Additional file 1: Figure S2B). The heterogeneity differences in moist skin group were largely driven by samples from the gluteal fold (Additional file 1: Figure S1E), as there was little difference in heterogeneity for the axilla (Additional file 1: Figure S1D). Interestingly, the scalp skin microbiome displayed no differences in heterogeneity among different disease states. These results indicate increasing beta heterogeneity for all dry skin sites in psoriasis and for the gluteal fold in psoriasis. Taken with the previous results for alpha diversity, there appears to be an overall loss of stability in the skin microbial community as psoriatic disease progresses, particularly for dry skin sites.
Psoriasis skin microbiota is enriched for Staphylococcus aureus and Staphylococcus pettenkoferi
We next examined bacterial composition of skin microbial communities from psoriatic and healthy skin at various taxonomic levels. The skin microbiome of all disease states (healthy, psoriatic lesional, and psoriatic non-lesional) consisted of four dominant phyla: Actinobacteria (53.8–66.5%) Firmicutes (23.9–28.3%), Proteobacteria (5.8–12.0%), and Bacteroidetes (2.1–2.9%) (Fig. 2a), consistent with previous descriptions of skin microbiota composition . At the genus level, skin microbiome is dominated by Propionibacterium (22.8–38.1%), Corynebacterium (21.4–23.9%), Staphylococcus (5.3–9.2%) in all disease states (Fig. 2b). Although the dominant taxa are similar in different disease states, we observed a gradual shift of taxonomic composition from healthy skin to psoriatic non-lesional skin to psoriatic lesional skin at both phylum and genus levels, suggesting that these microbial community shifts may precede the appearance of lesions in patients and have potential roles in disease progression. To further associate the taxonomic shift to different disease states, we identified bacterial taxa that discriminate each disease group using Lefse . At the phylum level, Actinobacteria and Proteobacteria served as strong discriminants for the skin microbiome from healthy and psoriatic lesions respectively (Fig. 2c). Lefse identified three genera, Propionibacterium, Ethanoligenens, and Macrococcus, as additional discriminative signatures for healthy skin (Fig. 2d). Lefse also identified 18 microbial genera that are discriminatively associated with psoriatic lesional skin including the genus Pseudomonas, which includes many opportunistic pathogens (Fig. 2d, Table 6). Four genera, Conchiformibius, Lactococcus, Moraxella, and Acetobacter, were associated with psoriatic unaffected skin (Fig. 2d). The combination of these genera can serve as potential markers for distinguishing skin from different disease states.
Our 16S rRNA sequencing also provided species-level resolution for some but not all of the sequencing reads. Lefse analysis identified several species-level bacterial signatures specific for different disease states (Fig. 2e). Consistent with what we observed in the genus level, the healthy skin microbiome was more enriched in both Propionibacterium acnes (P. acnes) and Propionibacterium granulosum (P. granulosum) compared to the psoriasis-associated skin microbiome (Fig. 3a, b). Staphylococcus sciuri was enriched in psoriatic non-lesional skin (Fig. 3c). Interestingly, two Staphylococcus species, S. aureus and S. pettenkoferi were significantly enriched in the psoriatic lesions while the genera Staphylococcus as a whole was not significantly enriched with any skin condition in our analysis (Fig. 2d, e).
We further explored the relative abundance of the Staphylococcus species across all samples with different disease states (Table 7). Strikingly, Staphylococcus aureus was more abundant in both lesional and non-lesional psoriatic skin compared to healthy skin (Fig. 3d). Although a low level of S. aureus was detected in 102 out of 147 healthy control samples and was detected in at least one skin swab of every healthy control subject, increased S. aureus abundance was exclusively observed in psoriasis samples (Fig. 3d). A similar trend was observed for Staphylococcus pettenkoferi, although to a lesser degree (Fig. 3e). In contrast, Staphylococcus epidermidis was more abundant in healthy skin compared to psoriatic skin (Fig. 3f) which is consistent with the previously reported competitive relationship between the Staphylococcus epidermidis and Staphylococcus aureus . The dynamic inter-microbe relationship between different Staphylococcus species might contribute to the distinct microbial communities associated with healthy and psoriatic skin.
Anatomic skin site is one of the major determinants of skin microbiome composition [24, 35]. Therefore, we further used Lefse to identify bacterial species at each skin site associated with healthy, non-lesional psoriatic, and lesional psoriatic skin (Table 8). We found that a reduced abundance of P. acnes is associated with psoriasis lesional skin at the arm, trunk, and gluteal fold (Additional file 1: Figure S3A), with a similar trend for the scalp and axilla. We did not observe a decrease for P. acnes in leg psoriasis samples, which is possibly due to the low abundance of P. acnes in healthy leg skin (Additional file 1: Figure S3A). Together, our data suggest that P. acnes may play a crucial role to maintain skin health at most skin sites besides the leg. Surprisingly, we did not observe a statistically significant increase in S. aureus abundance in psoriasis compared to healthy skin at any individual skin site (Table 8), whereas when anatomic sites were combined, S. aureus was highly associated with psoriasis lesional skin (Figs. 2e and 3d). Therefore, we defined a group of psoriasis samples with S. aureus abundance above the highest level of S. aureus colonization in healthy skin (baseline level = 0.0068) as “S. aureus high samples” (Table 9). We found that S. aureus high samples were observed exclusively in psoriasis patients and were seen at all skin sites (Additional file 1: Figure S3C), but that the number of S. aureus high samples at each skin site is modest, between 2 and 8 (Additional file 1: Figure S3D). This indicates that the association of S. aureus with psoriasis is not driven by any single anatomic site and that the presence of abundant S. aureus in only a subset of psoriasis patients, at different anatomic locations, results in an underpowered sample size for detection of S. aureus at any single body site.
Correlations between different bacterial species
Like any ecosystem, the composition of skin microbiome is modulated by both environmental factors (i.e., nutrient availability and host immune response) and interactions between different bacterial species. Inter-microbial interactions can be a major driver of microbial community composition, and understanding this interaction can yield important insights regarding the establishment and maintenance of psoriasis-associated microbial communities. We further investigated this microbe-microbe interaction by correlating microbial abundances with each other. At the genus level, we identified three clusters of bacterial communities, each constituting a group of bacteria significantly correlated in abundance (Fig. 4a). Cluster A was the largest cluster and consisted of Corynebacterium, Porphyromonas, Prevotella, Peptoniphilus, Finegoldia, and Anaercoccus. Cluster B was composed of Kocuria, Paracoccus, Micrococcus, and Janibacter. Lastly, Cluster C consisted of strongly correlated Streptococcus and Rothia. Given the previous reports of the potential role of Streptococcus in driving psoriasis [9, 10], it would be interesting to further investigate the role of Rothia spp. in psoriasis since it is highly co-abundant with Streptococcus. At the species level, P. acnes, which was more abundant in healthy skin, was negatively correlated with S. sciuri and S. pettenkoferi, both of which were enriched in the skin microbiota of psoriasis patients (Fig. 4b). Consistent with this observation, we also found P. acnes and S. epidermidis to be significantly enriched in S. aureus low psoriasis samples and S. pettenkoferi was enriched in S. aureus high psoriasis samples (Additional file 1: Figure S3E), suggesting that the antagonistic interaction among these bacteria may contribute to pathogenesis. Interestingly, Pseudoclavibacter bifida was negatively correlated with P. acnes and positively correlated to S. sciuri (Fig. 4b). The abundance of Pseudoclavibacter bifida was also enriched in S. aureus high psoriasis samples (Additional file 1: Figure S3E). Moreover, P. acnes and P. granulosum serve as two predominant Propionibacterium species and our data shows that they are positively correlated with each other (Fig. 4b). The strong co-correlation of P. acnes and P. granulosum and their association to healthy skin suggests that these Propionibacterium spp. may have a role in maintaining skin health.
Psoriatic lesions are characterized by thick and highly inflamed skin plaques, so the psoriatic lesions, psoriasis non-lesional skin, and healthy skin represent very distinct microbial habitats that may affect the quality of interactions between different microbes. Consistent with this, we observed distinct species correlations in these disease states, supporting the hypothesis that different microbe-microbe interactions occur in each disease state. We found the most numerous and strongest microbe-microbe correlations in healthy skin samples (Fig. 4c). Surprisingly, species correlations in the microbial community associated with psoriatic lesions (Fig. 4e) were more similar to those in healthy skin than psoriatic non-lesional skin (Fig. 4d). In healthy skin, P. acnes was negatively correlated with several bacterial species (Fig. 4c), suggesting it may inhibit the growth of these bacteria. Fewer microbes were negatively correlated with P. acnes in psoriatic non-lesional skin (Fig. 4d) and only Pseudoclavibacter bifida was anti-correlated with P. acnes in lesional skin (Fig. 4e). Overall, our data suggests the possibility that P. acnes may have a role in influencing the skin microbial community by keeping the growth of some microbes under control and that perturbation of this balance in psoriatic skin could serve as a potential disease driver.
S. aureus triggers Th17 immune response in a murine model
The increased prevalence of S. aureus in both lesional and non-lesional skin of psoriasis patients suggested the possibility that S. aureus might play a role in early stages of psoriasis pathogenesis. Despite its undesirable role in the context of psoriasis, the Th17 response serves as one of the major arms of host defense against bacterial infection through promotion of B cell activation and attraction of neutrophils [36, 37]. IL-17 is crucial in clearance of S. aureus at nasal, skin, and soft tissue sites . Previous studies have shown that S. aureus proteins promote Th17 differentiation in vitro , suggesting that colonization by S. aureus can lead to increased Th17 activation and IL-17 secretion. To assess the effect of S. aureus colonization on Th17 response in the skin, we performed skin colonization of newborn-specific pathogen-free (SPF) mice with S. aureus strain USA300 and assessed the cutaneous effector CD4+ T (Teff) cell response using RNAseq in comparison with SPF mice colonized with the commensal S. epidermidis, or un-colonized SPF controls (Table 10). We found significantly stronger Th17 transcriptomic signals in Teff cells sorted from the skin S. aureus-colonized mice. Teff cells isolated from mice exposed to S. aureus expressed significantly higher levels of IL-17A and IL-17F cytokine transcripts (Fig. 5a, b). IL-17A has been well characterized as one of the major drivers for psoriasis pathogenesis whereas IL-17F shares some redundant functions to IL-17A but its role in psoriasis is less defined . Besides IL-17, other components of Th17 responses including IL23R and IL22 were also increased upon S. aureus exposure (Fig. 5c, d). While S. aureus exposure during early life triggers a strong Th17 response in mice, the same treatment did not elicit consistent activation of a Th1 response (Additional file 1: Figure S4). Although S. aureus colonization has been strongly associated with atopic dermatitis, which is driven by Th2 responses [8, 40, 41], most components of the Th2 response such as IL-4, IL-5, and IL-13 were not induced by early colonization of S. aureus (Additional file 1: Figure S5). However, we did observe a strong induction in expression of the Th2-promoting transcription factor, GATA3 (Additional file 1: Figure S5E). Together, our data suggests that S. aureus colonization can specifically trigger activation of Th17 response, which might contribute to IL-17-driven inflammation in psoriasis.
In this study, we profiled the skin microbiota of psoriasis patients and healthy controls using the NIH standardized protocol and with higher sequencing depth to gain a more comprehensive understanding in psoriasis-associated microbiome. Our data demonstrate that the psoriasis skin microbiome is more heterogeneous compared to that of healthy skin. The compositional variance in the psoriatic skin community could be attributable to local environmental changes that accompany or immediately precede psoriatic disease. Proliferating keratinocytes in psoriasis patients are a rich source of antimicrobial peptides such as LL37, β-defensin, and psoriasin . The constant presence of these antimicrobial peptides could undermine equilibrium of the skin microbiome community and select for microbial species resistant to these antimicrobials. Based on our data, we speculate that a healthy skin microbial community consists of key stabilizer species, which may prevent growth of other species in the local microenvironment. In psoriatic skin, these stabilizer species may be outcompeted by invading pathogenic species and/or inhibited by chronic exposure to antimicrobial peptides, enabling colonization by pathogenic bacteria normally excluded from this niche. This could explain the higher heterogeneity that we observed in psoriatic skin. In contrast to our result, Alekseyenko et al. and Gao et al. observed decreased bacterial diversity in psoriatic skin compared with healthy skin [12, 14] while Fahlen et al. found no difference . Consistent with all the previous studies, we observed a decrease in relative abundance of Actinobacter associated with psoriasis skin [12,13,14]. Similar to Fahlen et al., we observed overrepresentation of Proteobacteria in psoriasis skin while both Gao et al. and Alekseyenko et al. showed a reduced abundance of Proteobacteria in psoriasis skin. These discrepancies might be due to the inherent heterogeneity in microbiota composition observed on the skin of psoriatic patients or to different experimental designs. This highlights the need to use standardized protocols among different studies to enhance reproducibility and to allow for meta-analysis of study cohorts. It is important to note that all studies mentioned above including ours profile the skin microbial community using an OTU (operational taxonomic unit) approach which groups reads from part of the 16S rRNA gene in order to account for artifact variance introduced by sequencing error. The major limitation of this approach is that by grouping different sequence variants, subtle inter-species variance can be sacrificed, which can reduce the resolution of taxonomical assignment at the species level . Despite the inherent limitation of OTU-based profiling, we were still able to gain species insights from our dataset.
Increased colonization of S. aureus in psoriatic skin has been reported previously in several small studies, but only a few of these examined unaffected skin from psoriasis patients . Our data revealed a significant increase of colonization in both psoriatic lesional and non-lesional sites compared to the baseline levels of S. aureus colonization found in healthy skin. This suggests that the increase in S. aureus is less likely a consequence of structural change in the skin from psoriasis but rather might be an important factor in initiating disease. Indeed, the potential role of bacterial infection in initiating and exacerbating psoriasis has been shown in Streptococcus infection [9, 10]. Our data from a murine model of Staphylococcal skin colonization suggests that cutaneous exposure to S. aureus triggers a strong Th17-type response in skin Teff cells as evident by induction of IL-17A, IL-17F, and IL22 cytokines. IL22 not only triggers a pro-inflammatory response, it also inhibits terminal differentiation of keratinocytes which is one of the characteristics of psoriasis. This suggests a potential capability of S. aureus to initiate psoriasis through upregulating a Th17 response. It is important to note that S. aureus consists of many strains and some strains are more virulent than others depending on their expression on a variety of different toxins and other molecules. Our murine experiments utilized S. aureus strain SF8300 from the USA300 lineage, which contains Panton-Valentine leukocidin (PVL) and phenol-soluble modulins (PSMs) contributing to its virulence in skin and soft tissue infections. Colonization of S. aureus has long been implicated in the pathogenesis of atopic dermatitis [8, 45]. Consistent with a recent study demonstrating that the capacity to induce Th2-type inflammation is limited to specific S. aureus strains isolated from severe atopic dermatitis patients, we observed little induction in components of Th2 response except for GATA3, which is the transcription factor that is required for Th2 polarization . Compared to the baseline level of S. aureus colonization in healthy controls, we found that S. aureus levels in some of the psoriasis samples were increased up to 90-fold. The increased S. aureus colonization in psoriasis could be due to the expansion of S. aureus strains seen in healthy controls or colonization with new S. aureus strains. However, due to lack of strain resolution of 16S rRNA-based profiling, we are not able to distinguish these two possibilities. In addition to strain-specific immunomodulating effects, our observation might suggest a temporal relationship between Th17 and Th2 polarization in response to S. aureus colonization and/or during neonatal development, as has been reported in murine models of atopic dermatitis . Future studies examining the S. aureus strain diversity in psoriasis skin and skin immune response to S. aureus strains specifically isolated from patients with psoriasis would be of significant interest in further dissecting the role of this bacterium in psoriasis disease pathogenesis.
Propionibacterium is one of the most dominant skin commensal bacteria . P. acnes has long been linked to acne vulgaris [47, 48], but recent studies suggest that P. acnes is also highly abundant in healthy skin and specific pathogenic P. acnes strains is one of the key determinants for acne vulgaris [49,50,51]. In this study, we found Propionibacterium to be abundant in healthy, lesional, and non-lesional skin but with highest abundance in healthy skin. Moreover, P. acnes and P. granulosum were two of the strongest microbial species associated with healthy skin. A possible explanation for the reduced abundance in Propionibacterium species in psoriatic lesions might be that reduced sebaceous content in psoriatic plaques contributes to a less favorable environment for Propionibacterium growth. The potential consequence of this reduction in Propionibacterium species in psoriasis is less clear. On one hand, P. acnes is known to produce propionate, a short chain fatty acid which can promote regulatory T cell in the colon , as well as RoxP (radical oxygenase of Propionibacterium acnes), which can potentially reduce oxidative stress and prevent skin inflammation . In contrast, certain strains of P. acnes isolated from acne patients have the potential to induce higher IL-17 production compared to strains isolated from healthy subjects based on unknown mechanisms . Although our current study does not have the resolution for P. acnes strain identification, identifying psoriasis-specific P. acnes isolates and compare them to those in healthy subject and assessing their differential genomic content and ability to modulate host T cell responses will be crucial in understanding whether the abundance or type of P. acnes in psoriasis patients contributes to their propensity for disease.
The association between Staphylococcus sciuri and psoriasis non-lesional skin is rather surprising since S. sciuri is better known as an animal-associated bacteria [55, 56]. S. sciuri has also been found in the human skin in both healthy and hospitalized individuals . The clinical relevance of S. sciuri has become important since several studies have isolated S. sciuri from hospitalized patients and methicillin-resistant strains of S. sciuri can be a health hazard for hospitalized patients [58,59,60,61,62]. Our study provides the first observation of S. sciuri in the context of psoriasis. It is possible that the S. sciuri carriers of our cohort obtained the bacteria from a previous hospital visit since S. sciuri has been found to be persistently present in the hospital environment . While the possible role of S. sciuri in psoriasis is unclear, we have observed an interesting pattern of S. sciuri in our cohort. Our data show that an increase in S. sciuri abundance is exclusively associated with psoriasis skin, particularly in non-lesional skin (Fig. 3c). Moreover, S. sciuri abundance is negatively correlated with P. acnes, which is highly enriched in healthy skin (Fig. 4b). Together, our results suggest S. sciuri may have a potential role in psoriasis pathogenesis.
In order to understand the roles of microbiome in human health, it is important to consider the microbial community as whole. Bacterial interactions are as important as the host environment in shaping the skin microbial community. The microbe-microbe relationship can be competitive or symbiotic. We performed a correlation analysis on the most abundant microbial genera and species to elucidate possible microbe-microbe interactions in healthy and psoriasis skin. We found an anti-correlation between S. sciuri and P. acnes, consistent to their respective disease state associations. To our surprise, although S. aureus is enriched in psoriatic skin and is known from other studies to have a competitive relationship with S. epidermidis and P. acnes, our data did not corroborate these negative associations. Possible explanations for this might include strain-specific interactions as well as the impact of the skin environment on inter-species interactions, the latter being supported by our finding that inter-bacterial correlation clusters differed by different disease state. It is important to note that microbe-microbe relationships suggested in our study are only correlative, further experimentation on isolated microbes will be needed to validate these relationships. Nonetheless, our work predicts strong candidates for the microbe-microbe relationships that may be crucial for psoriasis pathogenesis. Taken together, our correlation analysis demonstrates the highly complex relationship among skin bacteria by showing that these inter-microbial relationships are altered in psoriasis, possibly due to changes in the biochemical changes in skin environment and/or ecological pressure imposed by an altered host immune response.
In this study, we adhered to a stringent sampling protocol and measured skin microbiome profiles associated with psoriasis skin at six different skin sites. Our data revealed higher diversity and heterogeneity in psoriatic skin relative to healthy skin. Taxonomic analyses revealed specific microbial signatures associated with each disease state at the genus and species levels. Intriguingly, we found Staphylococcus aureus to be more abundant in both psoriatic non-lesional and lesional skin while Staphylococcus epidermidis, Propionibacterium acnes, and Propionibacterium granulosum were more abundant in healthy skin. We further tested the impact of Staphylococcus aureus colonization on host response in murine skin and validated its capacity for Th17 polarization. Finally, we demonstrated that disease state can alter microbe-microbe interactions and co-associations possibly due to differences in the physical and chemical environment of the skin. Our study confirms that psoriasis is accompanied by a shift in the skin-resident microbial community and raises intriguing possibilities worthy of further exploration for how this might directly impact the host immune response and psoriasis pathogenesis.
Twenty-eight adult psoriasis patients and 26 healthy volunteers recruited from the San Francisco Bay area were enrolled in the study after providing informed consent. Individuals with abnormal coagulation studies, positive HIV screening test, or a known history of bleeding disorders, abdominal surgery, gastrointestinal cancer, inflammatory bowel disease, AIDS, or other immunosuppressive condition, or concurrent inflammatory skin condition were also excluded. All psoriasis patients had a diagnosis of psoriasis from a physician for at least 6 months prior to study enrollment, which was verified by study staff. To assess the psoriatic microbiome in an untreated state, subjects were excluded if they had received systemic biologic therapy in the last 6 months, non-biologic systemic medications (methotrexate, cyclosporine, corticosteroids, cyclophosphamide, retinoids, photochemotherapy) or antibiotics in the last month, or phototherapy or topical therapy in the last 2 weeks prior to skin swabbing. Healthy volunteers had no personal or family history of psoriasis.
Skin swabs and stool samples were collected according to the protocol outlined in the Manual of Procedures for the NIH Human Microbiome Project . Study participants were asked to refrain from showering and using any substances on their skin (lotion, perfume, make-up, etc.) for at least 24 h prior to skin swabbing. Samples of the skin microbiome were collected using individually packed, sterile cotton swabs (Epicentre Catch-All Swabs). For each subject, the skin was swabbed at six standardized sites: scalp, trunk, axilla, arm, leg, and gluteal fold for healthy samples and psoriasis unaffected samples. Psoriatic lesional samples were only taken when psoriasis plaques were present at one of the six sites. Negative controls were obtained by exposing swabs to room air for 10 s. All samples were stored in − 80 °C while they awaited further processing.
DNA was extracted from the skin swab samples using the MasterPure Yeast DNA Purification Kit (Epicentre) with bead beating method to lyse the bacterial cells. To prepare skin microbiome library for sequencing, 16S rRNA were amplified at V1 to V3 hypervariable region using a universal forward primer (V1_27F primer): 5′-AGAGTTTGATCCTGGCTCAG-3′ attached to 5′ Illumina adapter and indexed reverse primer (V3_534R primer): 5′-ATTACCGCGGCTGCTGG-3′ attached with 3′ Illumina adapter sequence. Amplicon PCR reactions were completed as follows (per reaction): 2 μl of gDNA, 1× final concentration of 10× LA PCR Buffer ll (Mg2 + free) (Takara Bio USA), 0.4 mM dNTPs, 0.4 uM forward and reverse primers, 1.25 U of TaKaRa LA Taq polymerase high fidelity (Takara Bio USA), and nuclease-free water to bring the final volume to 25 μl. PCR cycling protocol consisted of an initial denaturation of 5 min at 95 °C, 30 cycles of 30 s at 95 °C, and 30 s at 56 °C followed by 5 min at 72 °C. PCR reactions were subsequently cleaned up using Agencourt AMPure XP kit (Beckman Coulter), and the purified amplicons were quantified using Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen). Samples with amplicon concentrations less than three times above the average air control amplicon concentrations were excluded from the subsequent sequencing. In average, 60–70 samples were pooled in equal molar quantities and the pooled library was purified using minElute PCR purification kit (Qiagen) for the final purification. The pooled library was sequenced on a Miseq sequencer (Illumina) as described in commercially provided protocol with 25% phiX DNA added as spike-in. Miseq reagent kit V3 (Illumina) was used to generate paired-end 300-bp reads. To avoid confounding from batch effects or possible external contamination that could affect our comparison between psoriasis and healthy skin microbiome, each Miseq sequencing run was balanced with samples from three or four healthy subjects and four psoriasis samples.
Data process and OTU picking
After quality check, the high-quality pair-end reads were assembled into ~ 550-bp fragments using FLASH  and we performed the subsequent analysis using Qiime scripts . The pair-end-assembled sequence fragments were first aligned against the Qiime supplied reference database for picking close OTUs. Subsequently, the unaligned sequences clustered into 97% identity operational taxonomical units (OTUs) using UCLUST . A representative sequence from each OTU cluster was aligned against the GreenGenes core set alignment template using PyNAST . The chimera due to the PCR errors were identified by ChimeraSlayer and excluded from the subsequent analyses . The remaining chimera-free OTUs were then used to approximate the phylogenetic tree using FastTree . We removed samples with less than 10,000 sequences to ensure adequate sample depth after rarefaction. The chimera-free sequences were rarefied into 11,766 per samples using the custom script analyses (https://github.com/alifar76/MicroNorm) to account for library size differences between samples.
Subsequent analysis for community diversities, microbial signatures, and inter-microbial correlations
The rarefied OTU table was used for all the subsequent analyses using Qiime 1.8.0  and R . Four alpha diversity metrics—chao1 index, observed OTUs, Shannon’s diversity index, and Simpson’s diversity index—were calculated for samples in each disease state, and the significance between different disease states were evaluated by a linear mixed effect model to account for multiple sampling from the same patient. The linear mixed effect model was performed by using R package lmerTest(v.3.0.1) . The R package Kendall (v.2.2)  was used to perform the Mann-Kendall trend test to detect significant trends of alpha diversity across different disease states.
Weighted UniFrac dissimilarity matrices were calculated to determine beta diversity using Qiime script. To account for multiple sampling of different body sites from the same subject, we applied linear mixed effects model to the first principal component coordinate from weighted UniFrac (scripts available in Additional file 2). The heterogeneity of microbial communities within each disease state were determined by the average weighted UniFrac distances between site-matched samples within each disease state, and the significance was determined by one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparisons test done by GraphPad Prism 7 (statistics output in Additional file 3). The relative abundance of each microbial community was summarized by Qiime script at the phylum and genus level. The microbial signatures associated to each disease state at both genus and species level were identified by Lefse, which combines non-parametric Kruskal-Wallis test and linear discriminant analysis (LDA) . Taxonomical features with p value < 0.05 and LDA effect size > 2.0 were considered significant microbial signatures. To evaluate the skin site- and skin type-specific diversity changes and microbial signatures, the OTU tables were further subsetted into different skin sites or different skin types. The diversity and marker analyses were performed as described above. We further investigated the difference between psoriasis samples with high S. aureus abundance and those with low S. aureus abundance seen in healthy skin. We defined the highest S. aureus abundance in healthy skin as baseline. Psoriasis samples (both lesional and non-lesional) with higher S. aureus abundance than the baseline level were considered “S. aureus high” and psoriasis samples with lower or equal S. aureus abundance as the baseline level were considered “S. aureus low”. Lefse analysis was performed to identify bacterial signatures associated with S. aureus high and S. aureus low samples. Lastly, the rarefied OTU table was summarized to genus (L6) and species (L7) by Qiime script and Spearman correlation was calculated using a R package, Hmisc (v. 4.1.1) , for the top 25 most abundant genera and top 30 most abundant species (scripts available in Additional file 4). The species-species correlations were further investigated within different disease states (scripts available in Additional file 4). Summarized taxonomic tables used in correlation analyses were available in Additional files 5 and 6.
Mouse skin bacterial colonization
C57BL/6 mice were purchased from Jackson Laboratories (Bar Harbor, ME) then bred and maintained in the UCSF-specific pathogen-free facility. All animal experiments were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals and the guidelines of the Laboratory Animal Resource Center and Institutional Animal Care and Use Committee of the University of California, San Francisco.
Bacterial strains Staphylococcus epidermidis Tü3298 and Staphylococcus aureus SF8300 were grown in tryptic soy broth for 24–48 h; pelleted and cellular mass from 2.5 ml of saturated culture was re-suspended in 100 μl of PBS for colonization of each animal. Newborn C57BL6 mice were colonized starting on day 3 of life and every other day thereafter until post-natal day 19 by pipetting 100 μl of bacterial suspension onto their skin and distributing evenly using a sterile PBS-soaked cotton-tipped swab.
Isolation and RNA sequencing of CD4+ effector T (Teff) and regulatory T (Treg) cells from skin of bacterially associated mice
To isolate skin T cells for staining and FACS sorting, mice were sacrificed at 21 days of age and the entire trunk skin harvested and lightly defatted. The skin was then minced with scissors and re-suspended in a 50 ml conical with 1–2 ml of digestion media comprised of 2 mg/ml collagenase XI, 0.5 mg/ml hyaluronidase, and 0.1 mg/ml DNase in RPMI with 1% HEPES, 1% penicillin-streptomycin, and 10% fetal calf serum. This mixture was incubated in a shaking incubated at 37 °C at 250 rpm for 45 min. An additional 15 ml of RPMI/HEPES/P-S/FCS media was then added, and the 50 ml conical was shaken vigorously by hand for 30 s. Another 15 ml of media was added, and then, the entire suspension was filtered through a sterile 100-mm cell strainer followed by a 40-mm cell strainer into a new 50 ml conical. The suspension was then pelleted, and the cell pellet was re-suspended in sort buffer (RPMI, 2 mM EDTA, 25 mM HEPES, 2% FBS) with U/ml RiboLock RNase inhibitor (Thermo Scientific) for staining for 30 min at 4 °C with fluorophore-conjugated antibodies specific for CD3, CD4, CD8, CD25, CD45, ICOS, TCRβ, and Tonbo Live-dead Ghost Dye. Teff (Live, CD45+, CD3int, CD4+, CD8neg, TCRβ+ CD25neg, ICOSneg) cells were then isolated via cell sorting on a MoFlo XDP (Beckman Coulter) in the UCSF Flow Cytometry Core. Cells were pelleted and flash frozen. RNA isolation was performed by Expression Analysis Q2 Solutions using QIAGEN RNeasy Spin columns and was quantified via Nanodrop ND-8000 spectrophotometer. RNA quality was checked by Agilent Bioanalyzer Pico Chip. The SMARTer Ultra Low input kit was used to generate cDNA libraries which were then sequenced to a 25-M read depth with Illumina RNASeq.
Reads were aligned to UCSC GRCm38/mm10 reference genome with STAR software (2.4.2a) . SAM files were generated with SAMtools from alignment results . Read counts were obtained with htseq-count (0.6.1p1) with the union option . Differential expression was determined using the R/Bioconducter package DESeq2 . RPKM table was generated by Cuffdiff Cuffdiff (2.2.1)  (Additional file 7).
Roberson EDO, Bowcock AM. Psoriasis genetics: breaking the barrier. Trends Genet. 2010;26:415–23. Elsevier Ltd
Afifi L, Danesh MJ, Lee KM, Beroukhim K, Farahnik B, Ahn RS, et al. Dietary behaviors in psoriasis: patient-reported outcomes from a U.S. national survey. Dermatol Ther (Heidelb). 2017;7:227–42. Springer
Cai Y, Fleming C, Yan J. New insights of T cells in the pathogenesis of psoriasis. Cell Mol Immunol. 2012;9:302–9.
Grice EA, Segre JA. The skin microbiome. Nat Rev Microbiol. 2011;9:244–53.
Belkaid Y, Tamoutounour S. The influence of skin microorganisms on cutaneous immunity. Nat Rev Immunol. 2016;16:353–66. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved
Naik S, Bouladoux N, Wilhelm C, Molloy MJ, Salcedo R, Kastenmuller W, et al. Compartmentalized control of skin immunity by resident commensals. Science. 2013;337:1115–9.
Ganju P, Nagpal S, Mohammed M, Nishal Kumar P, Pandey R, Natarajan VT, et al. Microbial community profiling shows dysbiosis in the lesional skin of vitiligo subjects. Sci Rep. 2016;6:18761.
Kong HH, Oh J, Deming C, Conlan S, Grice EA, Beatson MA, et al. Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis. Genome Res. 2012;22:850–9. Cold Spring Harbor Laboratory Press
Gudjonsson JE, Thorarinsson AM, Sigurgeirsson B, Kristinsson KG, Valdimarsson H. Streptococcal throat infections and exacerbation of chronic plaque psoriasis: a prospective study. Br J Dermatol. 2003;149:530–4. Blackwell Science Ltd
Telfer NR, Chalmers RJG, Whale K, Colman G. The role of streptococcal infection in the initiation of guttate psoriasis. Arch Dermatol. 1992;128:39. American Medical Association
Lowes MA, Suárez-Fariñas M, Krueger JG. Immunology of psoriasis. Annu Rev Immunol. 2014;32:227–55. NIH Public Access
Alekseyenko AV, Perez-Perez GI, De Souza A, Strober B, Gao Z, Bihan M, et al. Community differentiation of the cutaneous microbiota in psoriasis. Microbiome. 2013;1:31.
Fahlén A, Engstrand L, Baker BS, Powles A, Fry L. Comparison of bacterial microbiota in skin biopsies from normal and psoriatic skin. Arch Dermatol Res. 2012;304:15–22.
Gao Z, Tseng C, Strober BE, Pei Z, Blaser MJ. Substantial alterations of the cutaneous bacterial biota in psoriatic lesions. PLoS One. 2008;3:e2719.
Tett A, Pasolli E, Farina S, Truong DT, Asnicar F, Zolfo M, et al. Unexplored diversity and strain-level structure of the skin microbiome associated with psoriasis. NPJ Biofilms Microbiomes. 2017;3:14.
Loesche MA, Farahi K, Capone K, Fakharzadeh S, Blauvelt A, Duffin KC, et al. Longitudinal study of the psoriasis-associated skin microbiome during therapy with ustekinumab in a randomized phase 3b clinical trial. J Invest Dermatol. 2018; Elsevier. https://doi.org/10.1016/j.jid.2018.03.1501
Kong HH, Andersson B, Clavel T, Common JE, Jackson SA, Olson ND, et al. Performing skin microbiome research: a method to the madness. J Invest Dermatol. 2017;137:561–8. NIH Public Access
Jo J-H, Kennedy EA, Kong HH. Research techniques made simple: bacterial 16S ribosomal RNA gene sequencing in cutaneous research. J Invest Dermatol. 2016;136:e23–7. NIH Public Access
Cutting, Mary A, Mclnnes P. Manual of procedures for Human Microbiome Project core microbiome sampling Protocol A HMP Protocol # 07–001 Version 12.0. NC.fB Inf. 2010;
Kong HH. Details matter: designing skin microbiome studies. J Invest Dermatol. 2016;136:900–2. NIH Public Access
Oh J. Human, Bacterial and Fungal Amplicon Collection and Processing for Sequencing. Bio-protocol. 2015;5(10): e1477. https://doi.org/10.21769/BioProtoc.1477.
Meisel JS, Hannigan GD, Tyldsley AS, SanMiguel AJ, Hodkinson BP, Zheng Q, et al. Skin microbiome surveys are strongly influenced by experimental design. J Invest Dermatol. 2016;136:947–56. Elsevier Ltd
Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science. 2009;326:1694–7. NIH Public Access
Grice EA, Kong HH, Conlan S, Deming CB, Davis J, Young AC, et al. Topographical and temporal diversity of the human skin microbiome. Science. 2009;324:1190–2. NIH Public Access
Human Microbiome Project Consortium THMP. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14. NIH Public Access
Oh J, Byrd AL, Park M, Kong HH, Segre JA. Temporal stability of the human skin microbiome. Cell. 2016;165:854–66.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4. Nature Publishing Group
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65. Nature Publishing Group
Simonis N, Rual J-FJ, Carvunis A-R, Tasan M, Lemmens I, Hirozane-kishikawa T, et al. Empirically controlled mapping of the Caenorhabditis elegans protein-protein interactome network. Nat Methods. 2008;6:47–54.
Fredricks DN, Fiedler TL, Marrazzo JM. Molecular identification of bacteria associated with bacterial vaginosis. N Engl J Med. 2005;353:1899–911. Massachusetts Medical Society
Srinivasan S, Hoffman NG, Morgan MT, Matsen FA, Fiedler TL, Hall RW, et al. Bacterial communities in women with bacterial vaginosis: high resolution phylogenetic analyses reveal relationships of microbiota to clinical criteria. PLoS One. 2012;7:e37818. Ratner AJ, editor. Public Libr Sci
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35. American Society for Microbiology (ASM)
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60. BioMed Central
Otto M, Echner H, Voelter W, Götz F. Pheromone cross-inhibition between Staphylococcus aureus and Staphylococcus epidermidis. Infect Immun. 2001;69:1957–60. American Society for Microbiology
Perez Perez GI, Gao Z, Jourdain R, Ramirez J, Gany F, Clavaud C, et al. Body site is a more determinant factor than human population diversity in the healthy skin microbiome. PLoS One. 2016;11:e0151990. Public Libr Sci
McHeyzer-Williams M, Okitsu S, Wang N, McHeyzer-Williams L. Molecular programming of B cell memory. Nat Rev Immunol. 2011;12:24–34. NIH Public Access
Bröker BM, Mrochen D, Péton V. The T cell response to Staphylococcus aureus. Pathog (Basel, Switzerland). 2016;5 Multidisciplinary Digital Publishing Institute (MDPI)
Kolata JB, Kühbandner I, Link C, Normann N, Vu CH, Steil L, et al. The fall of a dogma? Unexpected high T-cell memory response to Staphylococcus aureus in humans. J Infect Dis. 2015;212:830–8. Oxford University Press
Ishigame H, Kakuta S, Nagai T, Kadoki M, Nambu A, Komiyama Y, et al. Differential roles of interleukin-17A and -17F in host defense against mucoepithelial bacterial infection and allergic responses. Immunity. 2009;30:108–19.
Byrd AL, Deming C, Cassidy SKB, Harrison OJ, Ng W-I, Conlan S, et al. Staphylococcus aureus and Staphylococcus epidermidis strain diversity underlying pediatric atopic dermatitis. Sci Transl Med. 2017;9:eaal4651. American Association for the Advancement of Science
Brandt EB, Sivaprasad U. Th2 cytokines and atopic dermatitis. J Clin Cell Immunol. 2011;2:110. NIH Public Access
Nestle FO, Kaplan DH, Barker J. Psoriasis. N Engl J Med. 2009;361:496–509.
Knight R, Vrbanac A, Taylor BC, Aksenov A, Callewaert C, Debelius J, et al. Best practices for analysing microbiomes. Nat Rev Microbiol. 2018;1 Nature Publishing Group
Totté JEE, van der Feltz WT, Bode LGM, van Belkum A, van Zuuren EJ, Pasmans SGMA. A systematic review and meta-analysis on Staphylococcus aureus carriage in psoriasis, acne and rosacea. Eur J Clin Microbiol Infect Dis. 2016;35:1069–77. Springer
Leyden JJ, Marples RR, Kligman AM. Staphylococcus aureus in the lesions of atopic dermatitis. Br J Dermatol. 1974;90:525. Blackwell Publishing Ltd
Oyoshi MK, Murphy GF, Geha RS. Filaggrin-deficient mice exhibit TH17-dominated skin inflammation and permissiveness to epicutaneous sensitization with protein antigen. J Allergy Clin Immunol. 2009;124:485–93, 493.e1. NIH Public Access
Bojar RA, Holland KT. Acne and propionibacterium acnes. Clin Dermatol. 2004;22:375–9.
Leeming JP, Holland KT, Cuncliffe WJ. The microbial colonization of inflamed acne vulgaris lesions. Br J Dermatol. 1988;118:203–8. Blackwell Publishing Ltd
Tomida S, Nguyen L, Chiu B-H, Liu J, Sodergren E, Weinstock GM, et al. Pan-genome and comparative genome analyses of propionibacterium acnes reveal its genomic diversity in the healthy and diseased human skin microbiome. MBio. 2013;4:e00003–13. American Society for Microbiology (ASM)
Fitz-Gibbon S, Tomida S, Chiu B-H, Nguyen L, Du C, Liu M, et al. Propionibacterium acnes strain populations in the human skin microbiome associated with acne. J Invest Dermatol. 2013;133:2152–60. NIH Public Access
Barnard E, Shi B, Kang D, Craft N, Li H. The balance of metagenomic elements shapes the skin microbiome in acne and health. Sci Rep. 2016;6:39491. Nat Publ Group
Smith PM, Howitt MR, Panikov N, Michaud M, Gallini CA, Bohlooly-Y M, et al. The microbial metabolites, short-chain fatty acids, regulate colonic Treg cell homeostasis. Science. 2013;341:569–73. NIH Public Access
Allhorn M, Arve S, Brüggemann H, Lood R. A novel enzyme with antioxidant capacity produced by the ubiquitous skin colonizer Propionibacterium acnes. Sci Rep. 2016;6:36412. Nat Publ Group
Agak GW, Kao S, Ouyang K, Qin M, Moon D, Butt A, et al. Phenotype and Antimicrobial Activity of Th17 Cells Induced by Propionibacterium acnes Strains Associated with Healthy and Acne Skin. J Invest Dermatol. Elsevier; 2018;138:316–24.
Kloos WE. Natural populations of the genus Staphylococcus. Annu Rev Microbiol. 1980;34:559–92. Annual Reviews 4139 El Camino Way, P.O. Box 10139, Palo Alto, CA 94303–0139, USA
Hauschild T, Schwarz S. Differentiation of Staphylococcus sciuri strains isolated from free-living rodents and insectivores. J Vet Med Ser B. 2003;50:241–6. Wiley/Blackwell (10.1111)
Couto I, Sanches IS, Sá-Leão R, de Lencastre H. Molecular characterization of Staphylococcus sciuri strains isolated from humans. J Clin Microbiol. 2000;38:1136–43. American Society for Microbiology (ASM)
Bansal S, Jain A, Agarwal J. Prevalence of methicillin-resistant, coagulase-negative staphylococci in neonatal intensive care units: findings from a tertiary care hospital in India. J Med Microbiol. 2004;53:941–4.
Lang S, Livesley MA, Lambert PA, Elliott J, Elliott TTJ. The genomic diversity of coagulase-negative Staphylococci associated with nosocomial infections. J Hosp Infect. 1999;43:187–93. WB Saunders
Petinaki E, Kontos F, Miriagou V, Maniati M, Hatzi F, Maniatis A. Survey of methicillin-resistant coagulase-negative staphylococci in the hospitals of central Greece. Int J Antimicrob Agents. 2001;18:563–6. Elsevier
Ahoyo TA, Yehouenou Pazou E, Baba-Moussa L, Attolou Gbohou A, Boco M, Dramane KL, et al. Staphylococcus sciuri outbreak at tertiary hospital in Benin. J Med Microbiol Diagnosis. 2013;2 OMICS International
de Sousa MA, Santos Sanches I, Ferro ML, de Lencastre H. Epidemiological study of Staphylococcal colonization and cross-infection in two West African Hospitals. Microb Drug Resist. 2000;6:133–41. Mary Ann Liebert, Inc
Dakić I, Morrison D, Vuković D, Savić B, Shittu A, Jezek P, et al. Isolation and molecular characterization of Staphylococcus sciuri in the hospital environment. J Clin Microbiol. 2005;43:2782–5. American Society for Microbiology (ASM)
Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63. Oxford University Press
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1. Oxford University Press
Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26:266–7. Oxford University Press
Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21:494–504. Cold Spring Harbor Laboratory Press
Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490. Poon AFY, editor. Public Libr Sci
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21. Oxford University Press
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9. Oxford University Press
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9. Oxford University Press
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. BioMed Central
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78. Nature Publishing Group
We thank the patients who generously provided samples for this study. This study was supported in part by a National Psoriasis Foundation Translation Research Award and NIH grants to WL (R01AR065174, U01AI119125) and National Psoriasis Foundation Fellowship Award to DY. TCS is supported by K08AR068409 and Burroughs Wellcome Fund CAMS-1015631. Flow cytometry was performed in the UCSF Parnassus Flow Cytometry Core which is supported by the Diabetes Research Center (DRC) grant, NIH P30 Dk063720. NGS sequencing was performed in Gladstone Genomics Core. We thank Mr. Kevin Lai for his computational assistance on this project.
Availability of data and materials
16S rRNA sequence data is deposited in the European Nucleotide Archive (http://www.ebi.ac.uk/ena) under study accession number PRJEB25915. RNAseq sequence data is deposited in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) under study accession number GSE114398.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Box plot shows the average weighted UniFrac distances among samples within each disease state at (a) arm, (b)trunk, (c)leg, (d)axilla, (e)gluteal fold and (f)scalp (*: p-value < 0.05, **: p-value < 0.01***: p-value < 0.001****: p-value < 0.0001). Figure S2. Box plot shows the average weighted UniFrac distances among samples within each disease state in (a) dry and (b)moist skin group (*: p-value < 0.05, **: p-value < 0.01***: p-value < 0.001****: p-value < 0.0001). Figure S3. Relative abundance of P. acnes and S. aureus in each disease state at different body sites. Relative abundance of P. acnes in healthy (red), psoriasis lesional (green), psoriasis non-lesional (blue) skin at (a) different skin sites and (b) different skin types (*: p-val < 0.05, **: p-val < 0.01, ****: p-val < 0.0001). (c) Box plot showing S. aureus abundance in S. aureus high samples. S. aureus high samples were defined as samples with higher S. aureus abundance than the highest S. aureus abundance among the healthy samples (baseline = 0.0068). (d) Bar graph depicts the prevalence of S. aureus high samples at each skin site in psoriasis lesional (blue bars) and psoriasis non-lesional (orange bars) skin. (e) Bacterial species associated with S. aureus high samples (red bars) and S. aureus low samples (green bars). Figure S4. Expression of Th1 components in effector T cells in response to Staphylococcus aureus colonization. The expression of Th1 components (a) T-bet (b) IFNγ and (c)IL-2 are comparable in all experimental groups. Figure S5. Expression of Th2 components in effector T cells in response to Staphylococcus aureus colonization. The expression of Th2 cytokines (a) IL-4 (b) IL-5, (c) IL-13 (d) IL-9 are comparable in all experimental groups. The expression of Th2 promoting transcription factor (e) GATA3 is induced by early colonization of S. aureus (adj. p-value = 1.49e-16); whereas, another Th2 transcription factor (f) STAT6 is not significantly induced. (PDF 405 kb)
R notebook for applying linear mixed effects model to test the significance of difference in microbial diversity among psoriatic lesional, psoriatic unaffected, and healthy skin. (HTML 772 kb)
Output of GraphPad Prism for statistics of heterogeneity of microbial communities isolated from skin associated with different disease state. (XML 18 kb)
R notebook for microbe-microbe correlation analyses at genus and species level. (HTML 1518 kb)
OTU table summarized to genus level used for correlation analysis. (CSV 1031 kb)
OTU table summarized to species level used for correlation analysis. (CSV 295 kb)
Mann-Kendall Test to detect trend in alpha diversity. (HTML 785 kb)