Skip to main content

Trait-based analysis of the human skin microbiome

Abstract

Background

The past decade of microbiome research has concentrated on cataloging the diversity of taxa in different environments. The next decade is poised to focus on microbial traits and function. Most existing methods for doing this perform pathway analysis using reference databases. This has both benefits and drawbacks. Function can go undetected if reference databases are coarse-grained or incomplete. Likewise, detection of a pathway does not guarantee expression of the associated function. Finally, function cannot be connected to specific microbial constituents, making it difficult to ascertain the types of organisms exhibiting particular traits—something that is important for understanding microbial success in specific environments. A complementary approach to pathway analysis is to use the wealth of microbial trait information collected over years of lab-based, culture experiments.

Methods

Here, we use journal articles and Bergeys Manual of Systematic Bacteriology to develop a trait-based database for 971 human skin bacterial taxa. We then use this database to examine functional traits that are over/underrepresented among skin taxa. Specifically, we focus on three trait classes—binary, categorical, and quantitative—and compare trait values among skin taxa and microbial taxa more broadly. We compare binary traits using a Chi-square test, categorical traits using randomization trials, and quantitative traits using a nonparametric relative effects test based on global rankings using Tukey contrasts.

Results

We find a number of traits that are over/underrepresented within the human skin microbiome. For example, spore formation, acid phosphatase, alkaline phosphatase, pigment production, catalase, and oxidase are all less common among skin taxa. As well, skin bacteria are less likely to be aerobic, favoring, instead, a facultative strategy. They are also less likely to exhibit gliding motility, less likely to be spirillum or rod-shaped, and less likely to grow in chains. Finally, skin bacteria have more difficulty at high pH, prefer warmer temperatures, and are much less resilient to hypotonic conditions.

Conclusions

Our analysis shows how an approach that relies on information from culture experiments can both support findings from pathway analysis, and also generate new insights into the structuring principles of microbial communities.

Background

The development of rapid, cost-effective sequencing technology has resulted in an explosion of microbiome research over the past decade. Microbial communities are now being sampled in almost every environment imaginable, ranging from the depths of the ocean [1, 2] to outer space [3, 4]. Reflecting the tremendous scope and magnitude of microbiome research are recent initiatives such as the Human Microbiome Project (HMP) [5,6,7,8,9] and the Earth Microbiome Project (EMP) [10,11,12]. The former aims to characterize all microbes on and in the human body, and the latter seeks to describe microbiomes across the entire globe. Already, discoveries from these and other, similar efforts are proving invaluable for understanding human disease [13,14,15,16], developing novel therapeutics [17, 18], and improving agricultural yields [19,20,21].

Existing microbiome research tends to focus on cataloging taxonomic diversity. Microbial function, by contrast, is less well studied [22, 23]. Unfortunately, without an understanding of microbial traits and, in particular, how traits differ among different environments, it is virtually impossible to answer key biological questions, like why certain microbes live where they do [24]. Trait-based analyses, which have a long history in macroscopic ecology [25,26,27], allow researchers to connect ecological traits to environmental associations, helping to explain the mechanisms underlying observed microbial distributions. The sheer diversity of typical microbiomes, however, makes trait-based analysis daunting.

Several strategies have been developed to circumvent challenges associated with trait-based microbial ecology. Shot-gun sequencing studies, for example, have been queried against reference databases, including COG/KOG, KEGG, eggNOG, Pfam, and TIGRFAM, to determine overrepresented genes, proteins, operons, and higher-order cellular processes [28,29,30,31,32,33,34,35] that reflect microbial function. Meanwhile, similar efforts have been extended to amplicon sequencing using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) [36] and Tax4Fun [37]—bioinformatics tools that infer microbial function based on reference databases, along with various assumptions about phylogenetic conservation. Although amplicon and shot-gun sequencing approaches appear comparable [37, 38], neither performs particularly well [38]—likely because of problems with the underlying reference databases, which are coarse-grained [38], represent only a minute fraction of microbial diversity, and are heavily biased toward a few organisms and environments [39]. More recently, machine learning techniques have been applied in an attempt correct for some of these problems and improve accuracy of trait prediction [40, 41].

Despite ongoing improvements in functional reference databases, the gold standard for defining microbial traits remains culture experiments. Decades of lab-based analyses have led to an impressive understanding of the functions of diverse microbial taxa, including many of those prevalent in microbiome studies. This information, however, is largely available through journal articles and Bergeys Manual of Systematic Bacteriology [42,43,44,45], neither of which is methodical in its presentation of data. Recently, there has been an effort to catalog trait information in more manageable and centrally organized databases, including StrainInfo [46], which collects trait data from biological resource centers and the JGI GOLD database, which allows users to enter known information on a handful of traits, including oxygen use, motility, and Gram stain. In addition, a recent text-parsing tool was developed that collects microbial descriptions from six separate sources, and then uses this information to predict microbial traits, including confidence scores [47]. The alternate, more precise but also more work-intensive approach is to link traits determined from lab- and culture-based experiments to output from microbiome sequencing studies directly, by manually curating every organism identified in a particular metagenomics sample. Although the effort involved is immense, if curation is done in a systematic fashion, then the resulting database has added, long-term value.

Here, we introduce such a trait database for human skin microbial communities, and then use it to characterize the bacterial residents of human skin in trait space. Bacterial traits are further compared to characteristics of bacteria more broadly using a similar database generated without any bias toward a particular habitat [48]. Finally, we compare traits across different skin environments to determine whether dry, moist, and sebaceous skin sites have functionally different microbial constituents. Many of the traits that we observe in skin microbiomes are consistent with expectations. For example, skin bacteria prefer warmer habitats and have higher salt requirements, in keeping with abiotic conditions on the skin surface. Several findings, however, suggest novel biological insight. Cocci, for example, are overrepresented on skin. Bacteria that form spores and possess phosphatases, by contrast, are underrepresented. Finally, relative to bacteria as a whole, skin bacteria are more likely to be anaerobic—a feature that is reflected not only in patterns of oxygen use, but also in distributions of oxidase and catalase activity, both of which are primarily beneficial in oxygen-rich environments.

Results

Trait composition of the human skin microbiome

Figure 1a presents binary traits for skin microbes. Spore formation is uncommon, particularly among abundant species, which are five times less likely to sporulate than skin microbes in general. By contrast, over half of skin taxa produce at least one pigment. Enzyme activities are varied. Whereas catalase is present in just under half of skin bacteria, oxidase, urease, alkaline phosphatase, gelatinase, and aesculin hydrolysis are less common, while acid phosphatase, α-galactosidase, arylsulfatase, pyrazinamidase, and tellurite reductase are rare. Catalase is the only enzyme more prevalent in abundant taxa. Gas production by skin bacteria is limited: almost no microbes generate methane, although a small fraction produces hydrogen sulfide and indole. Nitrate reduction is relatively common. This is in keeping with previous findings that skin commensals frequently reduce the nitrate in sweat [49, 50].

Fig. 1
figure 1

Proportion of all taxa (> 0.001% of reads in at least one sample; white) and abundant taxa (> 0.1% of reads in at least one sample; gray) in the human skin microbiome that exhibit (a) a range of different binary traits, (b) different types of oxygen use, (c) different types of motility, (d) different shapes, (e) different Gram stains and (f) different patterns of aggregation

Figure 1b–f presents categorical traits for skin microbes. The majority of skin microbes are facultatively anaerobic, although there are sizeable fractions of strictly aerobic and strictly anaerobic organisms as well. Most skin microbes are also non-motile, and this is particularly true of abundant taxa. Still, an unexpectedly large proportion—approximately 40%—have flagella. No other forms of motility are strongly represented. Most skin bacteria are rod-shaped and occur in clumps. Overall, skin microbes are predominantly Gram-negative, although abundant bacteria are split equally between Gram-negative and Gram-positive taxa.

Quantitative microbial traits are given in Table 1. Optimal temperature for growth is between 33.2 and 35.0 °C, which is close to the range of mean skin surface temperature, at 32.5–35.5 °C [51]. Optimal pH is near to neutral, even for abundant bacterial species. This is surprising, because the skin is an acidic environment, with pH values ranging from 4.0 to 7.0, but generally concentrated around pH ~ 5.0 [52,53,54]. In fact, low pH is thought to benefit commensal skin microbes, which adhere better to the skin surface under acidic conditions [54]. Optimal salt concentrations and salt concentration ranges are, likewise, well above salt concentrations measured in sweat [55]. We hypothesize that this may be explained by sweat evaporation at the skin surface, which can concentrate the salt from sweat. Mean GC content is approximately 50%.

Table 1 Mean quantitative trait data for all skin bacteria (>0.001% of reads in at least one sample) and abundant skin bacteria (0.1% of reads in at least one sample)

Figure 2 shows use of carbon substrates by skin bacteria. Here, we include all forms of use, including hydrolysis and fermentation. A wide range of carbon substrates are consumed by multiple skin taxa. This is particularly true of amino acids, with > 50% of the amino acids in our database used by > 70% of abundant skin taxa. Rates of use of monosaccharides and organic acids are lower, but still appreciable, with ~ 40% used by > 70% of abundant skin taxa. Use of alcohols and oligosaccharides/polysaccharides is less widely distributed, with 22% of oligosaccharides and no (0%) alcohols used by > 70% of abundant taxa. Of the carbon compounds considered, the substrates used most often by abundant taxa are glutamate (95%), asparagine (95%), valerate (92%), and glucose (91%).Footnote 1 The substrates used least are gelatin (3%), urea (17%), and xylitol (17%).

Fig. 2
figure 2

Proportion of all taxa (> 0.001% of reads in at least one sample; white) and abundant taxa (> 0.1% of reads in at least one sample; gray) in the human skin microbiome that utilize particular (a) organic acids, (b) amino acids, (c) monosaccharides, (d) oligosaccharides and polysaccharides, (e) alcohols and (f) other compounds

Comparing abundant versus rare skin bacteria, abundant taxa are more likely to use amino and organic acids. Eight amino acids (alanine, asparagine, aspartate, glutamate, glycine, leucine, proline, and serine; see Additional file 1: Supplemental Information II Table S2.3) are used more by abundant microbes than by the skin community as whole. Similarly, nine organic acids (acetate, citrate, formate, gluconate, malate, malonate, pyruvate, succinate, and valerate; see Additional file 1: Supplemental Information II Table S2.3) are used more by abundant microbes. For both amino acids and organic acids, all significant differences indicate that abundant skin taxa use these compounds more than skin taxa as a whole. Differences in consumption of other compounds, including alcohols and saccharides, are less biased toward overuse by abundant species. Indeed, two complex sugars (xylose and cellobiose) are used less by abundant taxa. Glucose, a simple sugar, on the other hand, is used more by abundant taxa (see Additional file 1: Supplemental Information II Table S2.3).

It is well known that certain taxonomic groups, for example Actinobacteria, are overrepresented among skin microbes and, in particular, among abundant skin microbes. While these groups are likely overrepresented because they have traits that make them uniquely adapted to the skin environment, it is possible that the traits that are important for living on skin are not those that we measured. Instead, the skin relevant traits may be other traits and the differences that we observe in the traits that we did measure may merely exist as a result of phylogenetic conservation. For this reason, we performed an additional analysis regressing the probability of a taxon being abundant versus rare against each trait individually, both for a naïve logistic regression and for a regression where phylogenetic relatedness was accounted for using the phylolm package in R [56]. To test the overall significance of a fitted regression, we compared it to a null model using a likelihood ratio test. In general, we found that many of the differences between abundant and rare taxa were preserved when phylogeny was accounted for. For instance, oxygen use, spore formation, Gram stain, type of motility, H2S production, the presence of catalase, aesculin hydrolysis and urease, and use of succinate, acetate, gluconate (organic acids), serine, proline, and glutamate (amino acids) were significantly different among abundant and rare taxa, whether or not phylogeny was considered. A few traits were not significant once phylogeny was included, for example cell shape, the presence of alkaline phosphatase, pyrazinamidase and gelatinase, and use of xylose, glucose, cellobiose (saccharides), malonate, formate, valerate, pyruvate, citrate, aspartate (organic acids), asparagine, alanine, leucine, and glycine (amino acids). Finally, use of 2-ketogluconate (organic acid) and the ability to perform nitrate reduction were only significant when accounting for phylogeny (see Additional file 1: Supplemental Information II, Table S2.1–S2.3).

Trait overrepresentation on human skin

Without comparison to prevalence in the world as a whole, it is impossible to know which traits are generally common versus preferentially selected for in skin environments. Figure 3a presents a comparison of binary traits among abundant skin bacteria versus bacteria more broadly (see “Materials and methods” section; see also Additional file 1: Supplemental Information III Fig. S3.1). Although there is a correlation between prevalence of a trait on skin and in the world as a whole, several traits are underrepresented among abundant skin taxa. Spore formation, for example is 7.5 times less likely among skin taxa as compared to general bacteria. Meanwhile, there is a 4.5-fold reduction in the likelihood of a skin taxon possessing acid phosphatase and a 1.5-fold reduction in the likelihood of a skin taxon possessing alkaline phosphatase as compared to bacteria more broadly. General bacteria are also 23% more likely to produce a pigment, 21% more likely to possess catalase, and 87% more likely to possess oxidase. For categorical traits, we again see significant differences between skin taxa and taxa from the world more broadly. Abundant skin bacteria (see Fig. 3b) are approximately half as likely to be aerobic, favoring, instead, a more flexible, facultative strategy. Likewise, abundant skin bacteria are 8-fold less likely to exhibit gliding motility, and none possess axial filaments, whereas these occur in ~ 0.1% of bacteria overall. Abundant skin taxa are also less likely to be spirillum or rod-shaped, whereas the fraction of cocci and coccibacilli on skin is inflated more than 2-fold. Finally, abundant skin bacteria are half as likely to grow in chains, preferring to aggregate as clumps instead.

Fig. 3
figure 3

Qualitative trait comparison for abundant taxa (> 0.1% of reads in at least one sample; see also Supplementary Information I). a Proportion of taxa with a specific, qualitative trait in skin microbial communities (x-axis) versus the world as a whole (y-axis). Filled symbols represent traits that are significantly different in skin environments; open circles represent traits that are not significantly different; marker size reflects significance. b Plots of trait proportions among skin bacteria (pink) and world bacteria (green). Open red circles denote traits that are overrepresented on skin; filled green circles denote traits that are overrepresented in the world (underrepresented on skin)

Figure 4 compares quantitative traits among world and skin bacteria (see also Additional file 1: Supplemental Information III, Figure S3.2). Abundant skin bacteria have more difficulty at high pH, tolerating, on average, a pH maximum of 7.97 versus 9.03 for the world in general. Abundant skin taxa also have a smaller range of pH values (2.41 versus 3.38) over which growth occurs. We speculate that this is because skin is a largely acidic environment with a relatively stable pH. Interestingly, however, optimal pH values for skin microbes do not reflect pH ranges measured on skin. Abundant skin bacteria also prefer warmer temperatures, can tolerate warmer temperatures, and have more difficulty at cold temperatures (with all three skin metrics being ~ + 2 °C) as compared to bacteria more broadly. Again, we hypothesize that this is because the skin is, at least relatively speaking, a warmer environment [48]. With respect to salt requirements, abundant skin bacteria are much less resilient to hypotonic conditions, requiring on average 1.1% NaCl, whereas average requirements in the world as a whole are closer to 0.02%. We speculate that this is because the skin is subject to constant excretion of salts through sweating. Finally, skin bacteria have a lower GC content (see also Additional file 1: Supplemental Information I, Figure S2), consistent with previous findings that host-associated organisms are AT-rich [57, 58].

Fig. 4
figure 4

Boxplots comparing quantitative traits among skin bacteria (pink) and bacteria from the world in general (green) for abundant skin microbes (> 0.1% of reads in at least one sample; see also Supplemental Information I). Blue stars are used to denote significant differences between a trait value in the world versus on skin. Box width indicates the relative number of microbes used for the comparison

We do not consider differences in carbon substrate usage between skin and the world because this information was collected differently in the skin database relative to the world database, making comparison impossible (see “Materials and methods” section).

Phylum level differences

As suggested above, one explanation for observed trends in functional traits on human skin is that these result from certain phyla (Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria) being the predominant constituents of the skin microbiome. To address this possibility, we used two separate approaches. First, we determined whether differences in functional traits between skin microbes and microbes more broadly persist when considering each phylum separately (see Tables 2, 3, and 4 and Additional file 1: Supplemental Information IV). For many traits—specifically, spore formation, pigment production, acid phosphatase, catalase (except for Actinobacteria), oxidase (see Table 2, Additional file 1: Table S4.1–S4.3), oxygen requirements, cell aggregation (see Table 3, Additional file 1: Table S4.4–S4.6), GC content, pH, and temperature requirements (see Table 4, Additional file 1: Table S4.7–S4.9)—biases that were apparent at the kingdom level are also apparent across multiple phyla. For other traits—for example, alkaline phosphatase, aeculin hydrolysis, and α-galactosidase (see Table 2, Additional file 1: Table S4.1–S4.3)—biases from the global composition appear driven by a single phylum, usually Proteobacteria, which is the most diverse phylum (see Additional file 1: Table S1.2) and thus most likely to impact overall results. Finally, for a few traits—most notably H2S and indole production (see Table 2, Additional file 1: Table S4.1–S4.3), motility, Gram stain, and cell shape (see Table 3, Additional file 1: Table S4.4–S4.6)—trends vary among phyla. Second, similar to our comparison of abundant versus rare taxa, we regressed the probability of a taxon being on the skin versus in the world more broadly against each trait individually using both a naïve logistic regression and a regression where phylogenetic relatedness was accounted for [56]. We then tested the overall significance of a fitted regression based on a null model using a likelihood ratio test. This analysis showed that all traits significantly over/underrepresented on skin relative to the world remained significant when accounting for phylogeny, while three traits (urease, pyrazinamidase, and nitrate reduction) were only significant under phylogenetic correction (see Additional file 1: Supplemental Information IV, Figure S4.10 and S4.11).

Table 2 Summary of binary trait results across dominant phyla from the human skin microbiome. Black is used for traits that are over-represented in the world; red is used for traits that are over-represented in the human skin microbiome. (See Table S3.1 for more detail)
Table 3 Summary of categorical trait results across dominant phyla from the human skin microbiome. Black is used for traits that are over-represented in the world; red is used for traits that are over-represented in the human skin microbiome. (See Table S3.2 for more detail)
Table 4 Summary of quantitative trait results across dominant phyla from the human skin microbiome. Black is used for traits that take on higher values in the world; red is used for traits that take on higher values in the human skin microbiome. (See Table S3.3 for more detail)

Trait differences among skin sites

Human skin microbiomes generally structure according to skin environment, with three environments—dry, moist, and sebaceous—represented (see Additional file 1: Supplemental Information I, Table S1.1). Because taxonomic composition differs among these three environments, functional diversity may vary as well. To test this hypothesis, we performed pairwise comparisons (dry vs. moist, dry vs. sebaceous, and moist vs. sebaceous) for all traits and substrate utilizations in our database (see Supplemental Information V). Surprisingly, not one difference emerged among skin environments for enzyme activities, gas production, spore formation, pigment production, nitrate reduction, Gram stain, cell aggregation, or pH, temperature, and NaCl requirements (see Additional file 1: Figure S5.1i, iii, S5.2i, iii, S5.3i, iii). Abundant bacteria at sebaceous sites are less likely to be rods as compared to abundant taxa at moist sites (49% versus 68%, see Additional file 1: Figure S5.3iv). As well, anaerobes are slightly underrepresented at dry sites as compared to sebaceous sites (see Additional file 1: Figure S5.2ii), and GC content is slightly lower at dry sites as compared to moist sites (see Additional file 1: Figure S5.5), although these latter two trends only emerge when considering the full skin microbiome, not just abundant taxa. Unfortunately, when accounting for phylogeny, the model for cell shape was degenerate for abundant taxa. However, variation in oxygen use between dry and sebaceous sites was observed even with phylogenetic correction. We did not attempt to control for phylogeny for GC content, since this was a quantitative trait.

Substrate usage (see Additional file 1: Supplementary Information V, Figure S5.6–S5.11) is similarly constant among skin environments, and what few differences do exist only occur between moist and sebaceous sites. Specifically, bacterial use of three organic acids—quinate, malonate, and caprate—as well as glucosamine (a monosaccharide) is overrepresented at sebaceous sites. By contrast, bacterial use of three saccharides—rhamnose, xylose, and cellobiose—as well as glycine (an amino acid) and urea are overrepresented at moist sites.

Our finding of high similarity among skin sites is in keeping with previous studies [6], but contrasts with a KEGG analysis performed in Oh et al. [59]. The discrepancy between our trait database analysis and the KEGG analysis may be because we considered a different set of functions. Alternatively, it may be because of differences in our definition of function prevalence. In particular, Oh et al. [59] quantified commonness of pathways across samples, whereas we quantify commonness of functions across taxa. Defining prevalence across species is not possible using pathway analysis, highlighting a distinction and benefit of our trait-based approach.

Discussion

We have undertaken a comprehensive trait-based analysis of the microbial constituents of human skin. In doing so, we have built an extensive trait-based database that will benefit future endeavors to characterize the functional properties of the skin microbiome. Below, we discuss some of our findings in terms of biological insight and interpretations.

Catalase, oxidase, and oxygen tolerance

Catalase is the most broadly distributed enzyme across the entire skin microbiome, and the only enzyme present in a significantly higher fraction of abundant skin taxa as compared to skin taxa as a whole. This suggests that catalase may be particularly beneficial for survival on skin, which should not be surprising. The majority of human skin is exposed to oxygen, while the role of catalase is to protect cells against hydrogen peroxide (H2O2)—an oxidant primarily generated as a result of reaction between oxygen and growth substrates [60]. Interestingly, however, catalase is still less common in skin bacteria as compared to bacteria as a whole. We speculate that this is because of the existence of one or more diverse, low-oxygen niches on human skin. Further evidence for such niches comes from the markedly lower prevalence of oxidase and the increased fraction of facultative and strict anaerobes and microaerophiles found on skin (see Additional file 1: Figure S3.1). One potential low oxygen niche is sebaceous follicles. These house the classic skin anaerobe, Propionibacterium acnes [61], and have been previously shown to be dominated by anaerobic taxa [62]. Sequencing studies, however, have pointed to low microbial diversity within follicles [63], which is not consistent with our finding that ~ 1/3 of culturable bacterial diversity on skin is either anaerobic or microaerophilic. Thus, we hypothesize that there are additional, low-oxygen environments hosting anaerobic taxa. One potential candidate is mixed-species biofilms [64]. Another is lower dermal layers, which may have been collected through scraping of the skin [59].

Several previous studies have considered the anaerobic portion of the skin microbiome, which is of interest because of its role in wound infections [65, 66]. These studies have found that counts of aerobes outnumber counts of anaerobes [67]. Although this may seem to contradict our conclusions, our analysis is based on diversity, rather than absolute counts. Based on our work, we theorize that, though anaerobes and microaerophiles may be less abundant, they must still be quite diverse. Consistent with previous findings, we observe evidence of increased anaerobicity among microbes at sebaceous sites (see Additional file 1: Figure S5.2) [67]. Similarly, our conclusion that anaerobes are less common at dry sites (see Additional file 1: Figure S5.2) accords with the KEGG analysis performed in [59], which found that dry sites harbored an abundance of citrate cycle modules.

Acid and alkaline phosphatases

Phosphatases enable bacteria to utilize certain components of soluble organic phosphorus [68], and thus are prevalent in environments where inorganic phosphorus is limiting. Almost 50% of microorganisms in soil and plant roots possess phosphatases [69,70,71]. By contrast, we find acid phosphatase in 7–8% of skin bacteria, and alkaline phosphatase in 12–13%; thus, we conjecture that phosphorus limitation is not significant in skin environments. This is surprising, because an experiment designed to measure loss of inorganic elements through healthy skin did not detect any phosphorus [72], nor is phosphorus abundant in human sweat [73, 74]. One explanation could be that skin bacteria rely on host-produced phosphatases [75, 76] to meet their needs. This would circumvent the metabolic cost of producing phosphatases, highlighting potentially unique aspects of microbial strategies in human-associated environments.

Spore formation

In a recent review article, Lennon and Jones [77] outlined factors promoting bacterial dormancy, with spore formation being an extreme case. Unlike the human gut, where few microbial genomes (~ 15%) show evidence of sporulation [77], human skin satisfies many of the conditions for dormancy. Skin, for example, is a highly inhospitable, exposed environment, lacking in resource availability [78]. By contrast, the gut is well-fed and generally protected. Furthermore, residence times on skin are long as compared to in the gut. Despite these differences, we find that the prevalence of sporulation is similar on skin and in the gut, both of which are significantly lower than rates among bacteria more broadly (see Fig. 3). Only ~ 20% of skin taxa produce spores, and this number is drastically lower (3%) when considering abundant taxa. Clearly, then, human microbiomes favor species without sporulation. We surmise that this is a result of the constant environment provided by host homeostasis.

Cell shape and aggregation

Relative to the broader world, skin microbiomes are enriched for cocci and coccobacilli (see Fig. 3). There are several hypotheses for why this might occur. First, rods allow for increased surface-to-volume ratios, improving nutrient uptake by passive diffusion [79] or when nutrients are directly acquired from a surface [80]. The fact that relatively fewer skin bacteria are elongated may thus indicate that nutrients on skin are readily available or, at the very least, are not acquired by passive diffusion (but see [81]). Second, although rods and filamentous cells are predicted to perform better under shear stress [82], cocci may be better able to fit into small pockets and pores of the stratum corneum. This is an alternate strategy for protection [82] that may be particularly advantageous on skin. Third, rod-shaped cells are more hydrodynamic, and thus can propel through liquid more efficiently [83]. This, however, may be of minimal importance in skin environments (although it is worth noting that rods appear to be enriched in moist regions). By contrast, cocci move much faster under conditions of Brownian motion [84]. Because skin bacteria frequently spread from one person to another through airborne release [85], a coccoid shape could facilitate interpersonal dispersal. Interestingly, coccoid cells can acquire some of the advantages of a rod shape (e.g., increased surface attachment) by growing in chains [82]. Despite this, chains, like rods, are underrepresented on human skin, further supporting our conclusion that skin selects for a spherical, rather than elongated shape.

Substrate utilization

Although many different substrates are consumed by skin bacteria, several stand out as being particularly important for success. Bacterial use of organic and amino acids, for example, shows enrichment in abundant skin bacteria. Interestingly, all eight of the amino acids that we find used significantly more by successful skin species have been positively identified in fingerprint samples [86]. This is consistent with our conclusion that these are important skin nutrients. Similar to amino acids, many of the organic acids that are used by a greater fraction of abundant skin taxa also appear commonly on human skin. This includes lactate, pyruvate [73], formate [87], caprate, and valerate [88]. In other cases, nutrients whose use is overrepresented among abundant taxa may not be produced by human skin, but rather, by dominant skin constituents. Succinate, for example, is a skin fermentation product of Staphylococcus epidermidis, meaning that it is likely widely available on the skin surface [89]. Further analysis of the chemical composition of skin secretions, not only by the human host but also by the entire skin microbiome, will help elucidate our findings regarding preferential substrate use.

Substrates that are less used by abundant skin taxa tend to be plant sugars, for example cellobiose [90], rhamnose [91], and xylose [92]. It is not difficult to understand why the ability to consume plant compounds provides little advantage on skin. Surprisingly, however, consumption of these sugars seems to be preferentially concentrated at moist sites, at least relative to sebaceous sites (see Additional file 1: Supplemental Information V, Figure S5.8 and S5.9). It is not obvious why there would be any benefit of plant sugar consumption in these regions. Urea use is also more common at moist sites (see Additional file 1: Supplemental Information IV, Figure S5.11), again for reasons that are unclear. In fact, urea use in general is surprising. Despite being prevalent on human skin [93], urea is one of the least commonly used substrates in our study (see Figs. 1 and 2). Why urea is not used by more skin bacteria, and why it seems to be used most at moist sites, highlights how trait-based analyses can uncover new, and unexpected trends, opening novel lines of inquiry that will ultimately help to elucidate factors governing skin microbiome composition.

Comparison to ProTrait

Both our database and the ProTrait database [47] draw from a vast literature of culture-based experiments. Whereas we manually curate our data, the ProTrait database uses a text-mining algorithm. Not surprisingly, our database contains information on fewer bacterial species (971 vs. 3046, with 25 unique to our database). Coverage of traits, however, is similar. We include several enzymes and carbon sources (for example arylsulfatase, pyrazinamidase, tellurite reductase, caprate, itaconate, suberate, succinate, urocanate, valerate, 3-hydroxybutyric acid, 3-hydroxybenzoate, asparagine, ornithine, phenylalanine, proline, threonine, tryptophan, glucosamine, methyl-B-d-glucoside, butanol, xylitol, 2,3-butanediol, carnitine, phenethylamine, putrescine, thymidine, uridine, and 2-aminethanol) that are not in ProTrait; however, the ProTrait database contains other enzymes and substrates that are not in our database. Interestingly, there do not appear to be significant differences in error rates between the two databases, at least for traits whose values are specified. The databases do, however, substantially differ in trait coverage. In particular, our database specifies the values of traits for a greater number of organisms, whereas the ProTrait database is more likely to report traits as unknown, at least using a precision of ≥ 0.9 (see Supplemental Information VI for several example comparisons).

Potential limitations

Our curated trait-based approach has many benefits, but also some draw-backs. First, we only consider well-defined taxa, ignoring detected taxa that have not been fully characterized, as well as all “dark matter” [59]. This could bias some of our predictions. While functional database methods are not as restricted in this way, they still rely on detection of orthologous genes. Consequently, both approaches are likely to miss at least some traits, particularly when these arise from poorly characterized taxonomic groups. Another complication of our approach is that it relies on conservation of functional traits within a species. Though our assumptions are likely less severe than tools like PICRUSt, functional traits are not always conserved. In compiling our database, we recorded evidence of strain variation, which suggested that interstrain differences in carbon source utilization are most common (14% of taxa), followed by differences in enzyme activity (11% of taxa). Although such variability complicates our analyses, it is more likely to obscure patterns than create them. Thus, when a pattern is detected, it likely reflects true biology.

Conclusions

Many opportunities exist for increased trait-based analysis of microbiome communities. Future studies considering additional human and non-human environments will help elucidate the structuring principles and biological mechanisms driving patterns in worldwide microbial distributions. Meanwhile, extended analyses of skin microbiomes will further highlight the principles governing community assembly. Analyses that quantitatively account for microbial abundance, for example, could clarify differences among dry, moist, and sebaceous sites, while further gradation by body location is also possible. Another extension would be to consider functional trait differences between different people—something that would be particularly informative when comparing individuals with skin disease to healthy controls.

Trait-based analyses and functional comparisons are the next step in microbiome research. Although most studies attempting to do this have taken a functional database/pathway analysis approach, culture and lab-based studies afford unique benefits. Our analysis of the skin microbiome has elucidated some of these benefits, detecting different patterns than were observed using KEGG [59]. This, in turn, has opened up a range of questions about why specific microbes exist in certain skin environments, and what they are doing to survive.

Materials and methods

Species list for the human skin microbiome

We defined a list of skin bacterial species using a recent study [59] that employed shotgun sequencing (see Additional file 1: Supplemental Information I, Table S1.1). Specifically, whole genome shotgun data from the NCBI Sequence Read Archive (SRA) project SRP002480 was obtained from the SRA FTP site and converted to paired-end FASTQ format using the splitsra script in our Git repository hosted at: https://bitbucket.org/skinmicrobiome/metagenomics-scripts. FASTQ data originating from the same BioSample were consolidated into the same file using a custom shell script and the SRA RunInfo table found here: http://www.ncbi.nlm.nih.gov/Traces/study/?acc=SRP002480.

A reference database was constructed for the Kraken classifier [94] using the complete genomes in RefSeq for the bacterial (2199 taxonomic IDs), archaeal (165 taxonomic IDs), and viral (4011 taxonomic IDs) domains, as well as eight representative fungal taxonomic IDs, the Plasmodium falciparum 3D7 genome, the human genome, and the UniVec Core database (ftp://ftp.ncbi.nlm.nih.gov/pub/UniVec). Low complexity regions of the microbial reference sequences were masked using the dustmasker program with a DUST level of 20 [http://www.ncbi.nlm.nih.gov/pubmed/16796549]. After masking, every 31-mer nucleotide sequence present in the collection of reference FASTA sequences was stored at the taxonomic ID of the lowest common ancestor among the leaf nodes that share that 31-mer (see [94] for details). The total size of the database plus index was 110 GB.

Each input read from SRA project SRP002480 was assigned a taxonomic ID using Kraken by finding exact matches between every 31-mer nucleotide sequence present in that read and the database of 31-mers constructed above. Because of the hierarchical storage of k-mers in the database, reads can be classified at more general taxonomic levels than the specific strain sequences that were used to build the database. Output from the Kraken classification was summarized by taxonomic ID along with the number of unique k-mers detected in the data using the kraken-report-modif script (present in the metagenomics-scripts repository linked above). The total number of unique k-mers for each taxonomic ID in the database was obtained using the count_kmers.pl script, and full taxonomic strings were generated using the taxid2taxstring script, both included in the metagenomics-scripts git repository linked above.

Two separate lists were constructed from the above output (see Additional file 1: Supplemental Information I, Table S3.1). The first list, representing all human skin taxa, was determined by recording any species that occurred in at least one sample with a relative abundance > 0.001% of reads. We set a lower bound on the percentage of reads because taxa with only a handful of reads may be spurious and/or may represent incorrect taxonomic assignments. The second list, representing abundant skin taxa, was determined by recording any species that occurred in at least one sample with a relative abundance of 0.1% of reads. We chose to consider abundance classes (all taxa vs. abundant taxa), rather than specifically accounting for abundance because abundance estimation from shotgun sequencing data is notoriously difficult.

Skin database compilation

Using the lists of taxa generated above, we compiled a database of microbial traits. For this, we relied on Bergeys Manual of Systematic Bacteriology [42,43,44,45] and the initial journal articles describing each species. We only considered validly described species and did not include Candidatus taxa, since little information was available for these. Our database contains information for 971 species.

World database compilation

We used a database compiled from species descriptions in the International Journal of Systematic and Evolutionary Microbiology. A full description of this database, including its availability, can be found at [48] (see also, Additional file 1: Supplemental Information I, Table S1.2).

Statistical analyses

Depending on the variable, we performed three types of comparisons: binary, categorical, and quantitative, across two sets of contrasts: skin vs. world and within skin bacteria, among the three skin environments: dry, moist, and sebaceous. These comparisons were conducted across all Bacteria and the four major phyla, separately considering abundant (> 0.1% of reads) and all taxa (> 0.001% of reads) respectively.

Binary comparisons were performed on variables that had two outcomes (e.g., positive and negative). When making two-way binary comparisons, we estimated proportion of occurrence with standard errors using a standard binomial model. For an overall test of difference in proportion, we used a Chi-square test. Pairwise comparisons were made using the standard errors of the binomial proportion. We visualized the comparisons with scatter plots of point estimates and error bars, using the 45° equality line as a guide for relative prevalence of the variables.

Categorical comparisons were performed on variables with multiple discrete, unordered outcomes (e.g., chain, clump, or singly). We compared the relative frequencies of the different outcomes in skin vs. world (or pairwise across skin environments) using a randomization test in which we resampled the data 105 times and computed a p value for the null hypothesis of equality of proportions by computing the number of randomized samples that were less extreme than the observed proportion.

Quantitative outcomes (e.g., volume, pH tolerance) were compared using a nonparametric relative effects test based on global rankings using Tukey contrasts [95]. We chose this test because it is robust to highly non-normal distributions and non-uniform variances and controls appropriately for multiple comparisons. We used box-and-whisker plots of each variable for visualization of the medians and deviations in the data.

Finally, to explore the role of phylogenetic conservation as an explanation for observed trends, for all binary and qualitative traits, we regressed the probability of a taxon being abundant versus rare or being from skin versus the world against each trait individually, both for a naïve logistic regression and for a regression where phylogenetic relatedness was accounted for. For the latter, we used the phylolm package in R [56] and the phylogenetic tree from Yarza et al. [96]. A handful of taxa were missing from the tree, and these were ignored in subsequent analysis. To test the overall significance of a fitted regression, we compared the regression to a null model using a likelihood ratio test. We then compared p-values for the naïve logistic regression and the regression with phylogenetic correction.

All statistical analysis was performed using the R programming language (R Code Team 2016), with the quantitative analysis performed using the nparcomp package [95].

Notes

  1. Note that we have ignored several compounds (e.g., carnitine, phenylethylamine, methyl-pyruvate) where results were only reported for a handful (< 20) of species.

References

  1. Huber JA, Welch DBM, Morrison HG, Huse SM, Neal PR, Butterfield DA, et al. Microbial population structures in the deep marine biosphere. Science. 2007;318(5847):97–100.

    Article  CAS  PubMed  Google Scholar 

  2. Xie W, Wang F, Guo L, Chen Z, Sievert SM, Meng J, et al. Comparative metagenomics of microbial communities inhabiting deep-sea hydrothermal vent chimneys with contrasting chemistries. The ISME journal. 2011;5(3):414–26.

    Article  PubMed  Google Scholar 

  3. Venkateswaran K, Vaishampayan P, Cisneros J, Pierson DL, Rogers SO, Perry J. International Space Station environmental microbiome—microbial inventories of ISS filter debris. Appl Microbiol Biotechnol. 2014;98(14):6453–66.

    Article  CAS  PubMed  Google Scholar 

  4. Ichijo T, Hieda H, Ishihara R, Yamaguchi N, Nasu M. Bacterial monitoring with adhesive sheet in the international space station-“Kibo”, the Japanese experiment module. Microbes Environ. 2013;28(2):264–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett C, Knight R, Gordon JI. The human microbiome project: exploring the microbial part of ourselves in a changing world. Nature. 2007;449(7164):804.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Consortium HMP. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.

    Article  CAS  Google Scholar 

  7. Peterson J, Garges S, Giovanni M, McInnes P, Wang L, Schloss JA, et al. The NIH human microbiome project. Genome Res. 2009;19(12):2317–23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Consortium HMP. A framework for human microbiome research. Nature. 2012;486(7402):215–21.

    Article  CAS  Google Scholar 

  9. Gevers D, Knight R, Petrosino JF, Huang K, McGuire AL, Birren BW, et al. The human microbiome project: a community resource for the healthy human microbiome. PLoS Biol. 2012;10(8):e1001377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Gilbert JA, Meyer F, Jansson J, Gordon J, Pace N, Tiedje J, et al. The earth microbiome project: meeting report of the “1 st EMP meeting on sample selection and acquisition” at Argonne National Laboratory October 6 th 2010. Stand Genomic Sci. 2010;3(3):249.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Gilbert JA, Meyer F, Antonopoulos D, Balaji P, Brown CT, Brown CT, et al. Meeting report: the terabase metagenomics workshop and the vision of an Earth microbiome project. Stand Genomic Sci. 2010;3(3):243.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Gilbert JA, Jansson JK, Knight R. The Earth microbiome project: successes and aspirations. BMC Biol. 2014;12(1):1.

    Article  CAS  Google Scholar 

  13. Cho I, Blaser MJ. The human microbiome: at the interface of health and disease. Nat Rev Genet. 2012;13(4):260–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13(9):1.

    Article  CAS  Google Scholar 

  15. 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(5):850–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Greenblum S, Turnbaugh PJ, Borenstein E. Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc Natl Acad Sci. 2012;109(2):594–9.

    Article  CAS  PubMed  Google Scholar 

  17. Bakken JS, Borody T, Brandt LJ, Brill JV, Demarco DC, Franzos MA, et al. Treating Clostridium difficile infection with fecal microbiota transplantation. Clin Gastroenterol Hepatol. 2011;9(12):1044–9.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Cani PD, Delzenne NM. The gut microbiome as therapeutic target. Pharmacol Ther. 2011;130(2):202–12.

    Article  CAS  PubMed  Google Scholar 

  19. Chaparro JM, Sheflin AM, Manter DK, Vivanco JM. Manipulating the soil microbiome to increase soil health and plant fertility. Biol Fertil Soils. 2012;48(5):489–99.

    Article  Google Scholar 

  20. Bhardwaj D, Ansari MW, Sahoo RK, Tuteja N. Biofertilizers function as key player in sustainable agriculture by improving soil fertility, plant tolerance and crop productivity. Microb Cell Factories. 2014;13(1):1.

    Article  Google Scholar 

  21. Barea J. Future challenges and perspectives for applying microbial biotechnology in sustainable agriculture based on a better understanding of plant-microbiome interactions. J Soil Sci Plant Nutr. 2015;15(2):261–82.

    Google Scholar 

  22. Stulberg E, Fravel D, Proctor LM, Murray DM, LoTempio J, Chrisey L, et al. An assessment of US microbiome research. Nat Microbiol. 2016;1:15015.

    Article  CAS  PubMed  Google Scholar 

  23. Waldor MK, Tyson G, Borenstein E, Ochman H, Moeller A, Finlay BB, et al. Where next for microbiome research? PLoS Biol. 2015;13(1):e1002050.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Green JL, Bohannan BJ, Whitaker RJ. Microbial biogeography: from taxonomy to traits. Science. 2008;320(5879):1039–43.

    Article  CAS  PubMed  Google Scholar 

  25. Schimper AFW. Plant Geography Upon a Physiological Basis, eds Groom P, Balfour IB, Fischer WR. Oxford: Clarendon Press; 1903.

  26. McGill BJ, Enquist BJ, Weiher E, Westoby M. Rebuilding community ecology from functional traits. Trends Ecol Evol. 2006;21(4):178–85.

    Article  PubMed  Google Scholar 

  27. Westoby M, Wright IJ. Land-plant ecology on the basis of functional traits. Trends Ecol Evol. 2006;21(5):261–8.

    Article  PubMed  Google Scholar 

  28. Tringe SG, Von Mering C, Kobayashi A, Salamov AA, Chen K, Chang HW, et al. Comparative metagenomics of microbial communities. Science. 2005;308(5721):554–7.

    Article  CAS  PubMed  Google Scholar 

  29. Gill SR, Pop M, DeBoy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, et al. Metagenomic analysis of the human distal gut microbiome. Science. 2006;312(5778):1355–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kurokawa K, Itoh T, Kuwahara T, Oshima K, Toh H, Toyoda A, et al. Comparative metagenomics revealed commonly enriched gene sets in human gut microbiomes. DNA Res. 2007;14(4):169–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Shabat SKB, Sasson G, Doron-Faigenboim A, Durman T, Yaacoby S, Miller MEB, et al. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. The ISME journal. 2016;10.12:2958.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2011:gkr988.

  33. Muller J, Szklarczyk D, Julien P, Letunic I, Roth A, Kuhn M, et al. eggNOG v2. 0: extending the evolutionary genealogy of genes with enhanced non-supervised orthologous groups, species and functional annotations. Nucleic Acids Res. 2010;38(suppl(1):D190–D5. https://doi.org/10.1093/nar/gkp951.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Finn RD, Mistry J, Schuster-Böckler B, Griffiths-Jones S, Hollich V, Lassmann T, et al. Pfam: clans, web tools and services. Nucleic Acids Res. 2006;34(suppl(1):D247–D51.

    Article  CAS  PubMed  Google Scholar 

  35. Selengut JD, Haft DH, Davidsen T, Ganapathy A, Gwinn-Giglio M, Nelson WC, et al. TIGRFAMs and genome properties: tools for the assignment of molecular function and biological process in prokaryotic genomes. Nucleic Acids Res. 2007;35(suppl(1):D260–D4.

    Article  CAS  PubMed  Google Scholar 

  36. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Aßhauer KP, Wemheuer B, Daniel R, Meinicke P. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics. 2015;31(17):2882–4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Xu Z, Malmer D, Langille MG, Way SF, Knight R. Which is more important for classifying microbial communities: who’s there or what they can do. ISME J. 2014;8(12):2357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Blaser MJ, Cardon ZG, Cho MK, Dangl JL, Donohue TJ, Green JL, et al. Toward a predictive understanding of Earth’s microbiomes to address 21st century challenges. mBio. 2016;7(3):e00714–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Weimann A, Mooren K, Frank J, Pope PB, Bremges A, McHardy AC. From genomes to phenotypes: Traitar, the microbial trait analyzer. MSystems. 2016;1(6):e00101–16.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Davis JJ, Gerdes S, Olsen GJ, Olson R, Pusch GD, Shukla M, et al. PATtyFams: protein families for the microbial genomes in the PATRIC database. Front Microbiol. 2016;7:118.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Garrity G, Staley JT, Boone DR, De Vos P, Goodfellow M, Rainey FA, et al. Bergey's manual® of systematic bacteriology: volume two: the Proteobacteria: Springer New York; 2006.

  43. Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, et al. Bergey's manual of systematic bacteriology: Volume 3: the Firmicutes: Springer New York; 2011.

  44. Bergey DHB, Garrity GM, Boone DR, Brenner DJ, Castenholz RW, Goodfellow M, et al. Bergey's manual of systematic bacteriology: the Bacteroidetes, Spirochaetes, Tenericutes (mollicutes), Acidobacteria, Fibrobacteres, Fusobacteria, Dictyoglomi, Gemmatimonadetes, Lentisphaerae, Verrucomicrobia, Chlamydiae, and Planctomycetes: Springer New York; 2011;4.

  45. Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Suzuki K-i, Ludwig W, et al. Bergey’s manual® of systematic bacteriology: volume five the Actinobacteria, part a: Springer New York; 2012.

  46. Verslyppe B, De Smet W, De Baets B, De Vos P, Dawyndt P. StrainInfo introduces electronic passports for microorganisms. Syst Appl Microbiol. 2014;37(1):42–50.

    Article  PubMed  Google Scholar 

  47. Brbić M, Piškorec M, Vidulin V, Kriško A, Šmuc T, Supek F. The landscape of microbial phenotypic traits and associated genes. Nucleic Acids Res. 2016:gkw964.

  48. Barberán A, Velazquez HC, Jones S, Fierer N. Hiding in plain sight: mining bacterial species records for phenotypic trait information. MSphere. 2017;2(4):e00237–17.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Lundberg JO, Weitzberg E, Cole JA, Benjamin N. Nitrate, bacteria and human health. Nat Rev Microbiol. 2004;2(7):593–602.

    Article  CAS  PubMed  Google Scholar 

  50. Weller R, Pattullo S, Smith L, Golden M, Ormerod A, Benjamin N. Nitric oxide is generated on the skin surface by reduction of sweat nitrate. J Investig Dermatol. 1996;107(3):327–31.

    Article  CAS  PubMed  Google Scholar 

  51. Ramanathan N. A new weighting system for mean surface temperature of the human body. J Appl Physiol. 1964;19(3):531–3.

    Article  CAS  PubMed  Google Scholar 

  52. Yosipovitch G, Xiong GL, Haus E, Sackett-Lundeen L, Ashkenazi I, Maibach HI. Time-dependent variations of the skin barrier function in humans: transepidermal water loss, stratum corneum hydration, skin surface pH, and skin temperature. J Investig Dermatol. 1998;110(1):20–3.

    Article  CAS  PubMed  Google Scholar 

  53. Zlotogorski A. Distribution of skin surface pH on the forehead and cheek of adults. Arch Dermatol Res. 1987;279(6):398–401.

    Article  CAS  PubMed  Google Scholar 

  54. Lambers H, Piessens S, Bloem A, Pronk H, Finkel P. Natural skin surface pH is on average below 5, which is beneficial for its resident flora. Int J Cosmet Sci. 2006;28(5):359–70.

    Article  CAS  PubMed  Google Scholar 

  55. Bates GP, Miller VS. Sweat rate and sodium loss during work in the heat. J Occup Med Toxicol. 2008;3(1):1.

    Article  Google Scholar 

  56. Tung Ho L, Ané C. A linear-time algorithm for Gaussian and non-Gaussian trait evolution models. Syst Biol. 2014;63(3):397–408.

    Article  Google Scholar 

  57. McCutcheon JP, Moran NA. Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol. 2012;10(1):13–26.

    Article  CAS  Google Scholar 

  58. Bentley SD, Parkhill J. Comparative genomic structure of prokaryotes. Annu Rev Genet. 2004;38:771–91.

    Article  CAS  PubMed  Google Scholar 

  59. Oh J, Byrd AL, Deming C, Conlan S, Kong HH, Segre JA, et al. Biogeography and individuality shape function in the human skin metagenome. Nature. 2014;514(7520):59–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Imlay JA. How oxygen damages microbes: oxygen tolerance and obligate anaerobiosis. Adv Microb Physiol. 2002;46:111–53.

    Article  CAS  PubMed  Google Scholar 

  61. Brüggemann H, Henne A, Hoster F, Liesegang H, Wiezer A, Strittmatter A, et al. The complete genome sequence of Propionibacterium acnes, a commensal of human skin. Science. 2004;305(5684):671–3.

    Article  PubMed  CAS  Google Scholar 

  62. Puhvel SM, Reisner RM, Amirian DA. Quantification of bacteria in isolated pilosebaceous follicles in normal skin. J Investig Dermatol. 1975;65(6):525–31.

    Article  CAS  PubMed  Google Scholar 

  63. Bek-Thomsen M, Lomholt HB, Kilian M. Acne is not associated with yet-uncultured bacteria. J Clin Microbiol. 2008;46(10):3355–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Sun Y, Smith E, Wolcott R, Dowd S. Propagation of anaerobic bacteria within an aerobic multi-species chronic wound biofilm model. J Wound Care. 2009;18(10).

    Article  CAS  PubMed  Google Scholar 

  65. Wolcott R, Gontcharova V, Sun Y, Zischakau A, Dowd S. Bacterial diversity in surgical site infections: not just aerobic cocci any more. J Wound Care. 2009;18(8):317–23.

    Article  CAS  PubMed  Google Scholar 

  66. Martin JM, Zenilman JM, Lazarus GS. Molecular microbiology: new dimensions for cutaneous biology and wound healing. J Investig Dermatol. 2010;130(1):38–48.

    Article  CAS  PubMed  Google Scholar 

  67. Selwyn S, Ellis H. Skin bacteria and skin disinfection reconsidered. Br Med J. 1972;1(5793):136–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Karl D, Craven D. Effects of alkaline phosphatase activity on nucleotide measurements in aquatic microbial communities. Appl Environ Microbiol. 1980;40(3):549–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Cosgrove D. Metabolism of organic phosphates in soil. Soil Biochem. 1967;1:216–28.

    Google Scholar 

  70. Tarafdar J, Claassen N. Organic phosphorus compounds as a phosphorus source for higher plants through the activity of phosphatases produced by plant roots and microorganisms. Biol Fertil Soils. 1988;5(4):308–12.

    Article  CAS  Google Scholar 

  71. Khan AA, Jilani G, Akhtar MS, Naqvi SMS, Rasheed M. Phosphorus solubilizing bacteria: occurrence, mechanisms and their role in crop production. J Agric Biol Sci. 2009;1(1):48–58.

    Google Scholar 

  72. Freyberg R, Grant R. Loss of minerals through the skin of normal humans when sweating is avoided. J Clin Investig. 1937;16(5):729.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Robinson S, Robinson AH. Chemical composition of sweat. Physiol Rev. 1954;34(2):202–20.

    Article  CAS  PubMed  Google Scholar 

  74. Mitchell H, Hamilton T, Haines W. The dermal excretion under controlled environmental conditions of nitrogen and minerals in human subjects, with particular reference to calcium and iron. J Biol Chem. 1949;178:345–61.

    Article  CAS  PubMed  Google Scholar 

  75. Loewenthal L, Politzer W. Alkaline phosphatase in human eccrine sweat; 1962.

    Book  Google Scholar 

  76. Loewenthal L, Hins S. Phosphomonoesterases in eccrine sweat and their localization after prolonged sweating. Dermatology. 1963;126(5):311–8.

    Article  CAS  Google Scholar 

  77. Lennon JT, Jones SE. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol. 2011;9(2):119–30.

    Article  CAS  PubMed  Google Scholar 

  78. Fyhrquist N, Salava A, Auvinen P, Lauerma A. Skin biomes. Curr Allergy Asthma Rep. 2016;16(5):1–7.

    Article  Google Scholar 

  79. Schulz HN, Jørgensen BB. Big bacteria. Annus Rev Microbiol. 2001;55(1):105–37.

    Article  CAS  Google Scholar 

  80. Steinberger R, Allen A, Hansma H. Holden pm. Elongation correlates with nutrient deprivation in Pseudomonas aeruginosa unsaturated biofilms. Microb Ecol. 2002;43(4):416–23.

    Article  CAS  PubMed  Google Scholar 

  81. Young KD. Bacterial morphology: why have different shapes? Curr Opin Microbiol. 2007;10(6):596–600.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Young KD. The selective value of bacterial shape. Microbiol Mol Biol Rev. 2006;70(3):660–703.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Cooper S, Denny MW. A conjecture on the relationship of bacterial shape to motility in rod-shaped bacteria. FEMS Microbiol Lett. 1997;148(2):227–31.

    Article  CAS  Google Scholar 

  84. Dusenbery DB. Fitness landscapes for effects of shape on chemotaxis and other behaviors of bacteria. J Bacteriol. 1998;180(22):5978–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Meadow JF, Altrichter AE, Bateman AC, Stenson J, Brown G, Green JL, et al. Humans differ in their personal microbial cloud. PeerJ. 2015;3:e1258.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  86. Croxton RS, Baron MG, Butler D, Kent T, Sears VG. Variation in amino acid and lipid composition of latent fingerprints. Forensic Sci Int. 2010;199(1):93–102.

    Article  CAS  PubMed  Google Scholar 

  87. Poulos A, Sharp P, Singh H, Johnson D, Carey W, Easton C. Formic acid is a product of the α-oxidation of fatty acids by human skin fibroblasts: deficiency of formic acid production in peroxisome-deficient fibroblasts. Biochem J. 1993;292(2):457–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Ara K, Hama M, Akiba S, Koike K, Okisaka K, Hagura T, et al. Foot odor due to microbial metabolism and its control. Can J Microbiol. 2006;52(4):357–64.

    Article  CAS  PubMed  Google Scholar 

  89. Wang Y, Kuo S, Shu M, Yu J, Huang S, Dai A, et al. Staphylococcus epidermidis in the human skin microbiome mediates fermentation to inhibit the growth of Propionibacterium acnes: implications of probiotics in acne vulgaris. Appl Microbiol Biotechnol. 2014;98(1):411–24.

    Article  CAS  PubMed  Google Scholar 

  90. Sensuła BM, Derrick PJ, Bickerton JC, Pazdur A. Mass spectrometric study of glucose and cellobiose produced during enzymatic hydrolysis of α-cellulose extracted from oak late-wood annual rings. Rapid Commun Mass Spectrom. 2009;23(13):2070–4.

    Article  PubMed  CAS  Google Scholar 

  91. Linhardt RJ, Bakhit R, Daniels L, Mayerl F, Pickenhagen W. Microbially produced rhamnolipid as a source of rhamnose. Biotechnol Bioeng. 1989;33(3):365–8.

    Article  CAS  PubMed  Google Scholar 

  92. Jeffries TW, Grigoriev IV, Grimwood J, Laplaza JM, Aerts A, Salamov A, et al. Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis. Nat Biotechnol. 2007;25(3):319–26.

    Article  CAS  PubMed  Google Scholar 

  93. Hamilton PB. Amino-acids on hands; 1965.

    Book  Google Scholar 

  94. Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15(3):R46.

    Article  PubMed  PubMed Central  Google Scholar 

  95. Konietschke F, Placzek M, Schaarschmidt F, Hothorn LA. Nparcomp: an R software package for nonparametric multiple comparisons and simultaneous confidence intervals. J Stat Softw. 2015;64(9):1–17.

  96. Yarza P, Richter M, Peplies J, Euzeby J, Amann R, Schleifer K-H, et al. The all-species living tree project: a 16S rRNA-based phylogenetic tree of all sequenced type strains. Syst Appl Microbiol. 2008;31(4):241–50.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Stuart Jones and Noah Fierer for guidance and for providing their database of world microbes.

Funding

This work was support by the U.S. Army Research Laboratory and the U.S. Army Research Office under contract/grant number #W911NF-14-1-0490 (DK).

Availability of data and materials

https://bewicklab.weebly.com/databases--packages.html.

Author information

Authors and Affiliations

Authors

Contributions

SB, DK, and WF conceived of the paper, and SB oversaw trait database generation. EG and JW performed all statistical analyses. Each of JB, CD, and RF entered data on over 200 organisms. TM, PT, and FB performed all analyses related to generation of skin taxon lists. SB, DK, WF, FB, and EG wrote the paper.

Corresponding author

Correspondence to Sharon Bewick.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Supplementary Information I through VI, containing information on the dataset and database, phylogenetically corrected analyses, additional analyses of full microbiomes, including rare members, body-site comparisons, and a comparison to ProTrait. (DOCX 43176 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bewick, S., Gurarie, E., Weissman, J. et al. Trait-based analysis of the human skin microbiome. Microbiome 7, 101 (2019). https://doi.org/10.1186/s40168-019-0698-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-019-0698-2

Keywords