The predominant method by which skin microbiomes have been studied is through sequencing of DNA from swabbed skin areas (12). However, despite the fact that the skin microbiome is often depicted as a well-mixed coating of microbes on the skin surface (8, 13), it remains unclear whether the DNA from the skin surface reflects the underlying biology of the skin microbiome. Additionally, while sequencing and culture-based studies have demonstrated that bacteria extend into deeper portions of the skin, the spatial distribution of the skin microbiome as a whole has not been well characterized (4, 14, 15). To address this gap, we determined the spatial distribution of bacterial cells in the skin using the universal bacterial FISH probe EUB338, which hybridizes to bacterial 16S rRNA (16).

We first used EUB338 FISH on biopsied healthy adult human facial tissues and found that the skin surface contains very few bacteria (Figure 1A). In contrast, clusters of bacteria were found within hair follicles and other cutaneous skin structures like comedos (Figure 1B). There are several possible explanations for the lack of FISH staining on the skin surface. For example, viable surface bacteria may have been eliminated by sterilization of the biopsied area (we address this possibility below). Alternatively, the surface bacteria could exist in a state like stationary phase that has fewer EUB388-hybridizable ribosomes, or by skin microbiome bacterial species being less prone to FISH probe hybridization. However, we confirmed that EUB338 staining works well in stationary-phase cells of multiple skin microbiome bacterial species (Cutibacterium acnes, Staphylococcus epidermidis, Micrococcus luteus, and Corynebacterium striatum) (Figure 1 – figure supplement 1).

To quantify our FISH staining of human skin, we calculated the ratio of the mean fluorescence within an area of interest (a hair follicle, the skin surface, or other cutaneous structures) to the mean fluorescence outside of the area of interest. We refer to this ratio as an enrichment score. The median enrichment score for human follicles was 11.24 (Figure 1D). In contrast, the skin surface (stratum corneum) had an enrichment score of just 0.188. This pattern of observing many bacteria in follicles but few on the skin surface held true when using FISH probes specific for C. acnes, one of the most abundant bacterial species in the human skin microbiome (Figure 1C). For C. acnes, the enrichment score was 3.14 for follicles as compared to 0.67 for the stratum corneum (Figure 1D). We note that while the majority of intact skin-associated bacteria were not associated with the skin surface, there are some visible bacteria on the skin surface, especially in areas near sub-surface structures (Figure 1A-B). This indicates that FISH staining is capable of detecting bacteria on the skin surface. Thus, our results are consistent with previous reports that bacteria can be cultured from skin surface swabs but also extend these findings to demonstrate that the skin surface has fewer intact bacteria than deeper skin structures.

DNA sequencing of samples collected by swabbing the skin surface is the most common method used to evaluate the skin microbiome because it is simple, noninvasive, and has been shown to result in higher consistency than other sampling methods (17). However, traditional metagenomic and 16S rRNA gene sequencing do not differentiate DNA from viable and inviable bacteria. We thus implemented a method that allowed us to quantitatively determine the fraction of DNA in a sample that is specifically from intact bacteria. The ability to molecularly differentiate between viable and inviable cells would enable us to both validate our FISH results using an orthogonal method on untreated, living skin, and address the possibility that the low numbers of intact surface bacteria observed by FISH was due to biopsy sterilization. To this end, we utilized the cell-impermeable small molecule propidium monoazide (PMA), which binds irreversibly to double-stranded DNA upon photoactivation to inhibit PCR amplification (Figure 2A) (18). When PMA photoactivation is performed before the cell lysis step of DNA isolation, the genomic DNA inside viable bacteria is protected from PMA binding because PMA is cell-impermeable, while cell-free DNA or DNA within permeabilized bacteria becomes PMA-bound. To quantify bacterial DNA, we combined the use of PMA with droplet digital PCR (PMA-ddPCR). Calculating the ratio of ddPCR counts between samples without PMA and samples with PMA allowed us to generate a viability score for any given bacterial population. A similar approach has been used to assess the viability of bacteria from environmental waste water samples (19).

To assure that PMA-ddPCR would allow us to reliably gauge the fraction of viable cells in a population, we first validated that it generated the expected results using known ratios of either exponentially-growing or stationary-phase E. coli cultures combined with heat-killed E. coli cells (Figure 2 – figure supplement 1A). To determine if PMA-ddPCR also works for skin-resident bacterial species, we also applied PMA-ddPCR to four of the most common skin microbiome bacteria, S. epidermidis, C. acnes, M. luteus, and C. striatum. We performed serial dilutions of each bacterial culture and in parallel performed both CFU plating and ddPCR with or without PMA. In every case, the amount of DNA in a sample (assessed by PMA-ddPCR) correlated well with the number of culturable bacteria in high and low bacterial abundance scenarios (Figure 2 – figure supplement 1B). Together, these controls confirm that PMA-ddPCR is a reliable method for assessing the amount of DNA in a sample present within intact bacteria.

We next applied PMA-ddPCR to human skin microbiome samples by swabbing the skin of four healthy human volunteers at eight sites (glabella, retroauricular crease, lower back, hair shaft, antecubital fossa, popliteal fossa, nares, and dorsal forearm) (Figure 2B, Figure 2 – figure supplement 2A, B). PMA-ddPCR revealed that the viability scores for these sites ranged between 0.02 and 0.12 (0 represents a fully-nonviable population, 1.0 represents a fully-viable population), indicating that the majority of bacterial DNA found on the skin surface is not associated with viable cells (Figure 2C). To investigate whether this was a skin-specific phenomenon, we tested several non-skin microbiome sites (tongue, saliva, plaque, and feces). We found that in all non-skin microbiome sites, the viability score was significantly higher than for the skin, ranging from 0.4 (saliva) to 0.87 (feces) (Figure 2C).

We next sought to address whether our PMA-ddPCR viability scores accurately represent of the number of viable bacteria on the skin surface. The low number of detectable bacteria on the skin surface (either by FISH or PMA-ddPCR) does not imply that the skin surface is sterile but rather that the majority of the bacterial DNA on the skin surface is not from these viable cells. Indeed, culturing bacteria directly from the skin is common (8). To quantify the number of viable bacteria directly, we plated a small amount of each sample using the standard conditions for culturing skin microbes (5% sheep blood in tryptic soy agarose plates incubated both aerobically and anaerobically). To determine whether PMA-ddPCR or traditional ddPCR better represented the number of viable skin microbiome bacteria, we compared our results to a standard curve generated with known numbers of S. epidermidis. For each sample, the PMA-ddPCR quantification closely matched this standard curve while the samples lacking PMA showed no overlap (Figure 2 – figure supplement 1C, D). The highest-abundance skin microbiome species are readily culturable, such that quantifying the DNA from viable skin surface bacteria should be able to accurately predict colony forming units (CFU) upon plating. Quantifying the bacterial DNA in skin microbiome samples without the use of PMA resulted in DNA quantities that were, on average, 82 times higher than predicted by the standard curve, while the use of PMA brought this value down to just 1.3 (Figure 2 – figure supplement 1E). Using ddPCR counts to predict CFU showed similar results, as ddPCR in the absence of PMA yielded values that predicted CFU counts 58.5 times greater than those measured, while PMA-ddPCR yielded values that predicted CFU counts that were on average only 1.28 times greater than the actual cultured CFU (Figure 2 – figure supplement 1F). Since we expect a ratio near 1, these findings lend further support for our conclusion that PMA-ddPCR reflects the viable microbiome better than approaches lacking PMA.

Our FISH results from biopsied skin suggested that viable bacteria may be protected below the skin surface. To test this hypothesis in untreated skin, we collected facial follicle contents and used PMA-ddPCR to compare the viability of facial follicle contents and skin surface samples (Figure 2D). The average viability for follicle contents was greater than 10-fold higher for follicle contents than for the skin surface, confirming that cells in hair follicles are more viable than those on the skin surface. Together, our results from FISH on biopsies and PMA-ddPCR on skin swabs independently support the conclusion that the skin surface is populated by few viable bacterial cells, indicating that the surface of healthy, non-sterilized human skin is sparsely colonized. While these data indicate that bacterial DNA on the skin surface is predominantly not associated with viable bacterial cells, we note that our results do not suggest that there are no viable cells on the skin surface. Rather, our data indicate that the majority of the bacterial DNA on the skin surface is not within bacteria such that using PMA provides a much more accurate estimation of the viable skin microbiome.

Like many microbiomes, the existing knowledge of the skin microbiome is heavily based upon bacterial 16S rRNA gene amplicon sequencing, which was developed to assess bacterial populations while avoiding biases introduced by culturing methods. However, our findings suggest that using 16S rRNA gene amplicon sequencing to study the skin microbiome is not entirely unbiased, as most of the DNA in these samples is not from viable bacteria and traditional 16S rRNA gene amplicon sequencing does not differentiate between DNA originating from live or dead cells. The inability of 16S rRNA gene amplicon sequencing to differentiate between these two types of bacterial populations has been mentioned as a potential downfall of the method (8). To evaluate how accurately traditional 16S rRNA gene amplicon sequencing captures the living skin microbiome composition, we utilized PMA followed by 16S rRNA gene amplicon sequencing (an approach we refer to as PMA-seq) (20). By sequencing pairs of matched samples with PMA treatment (PMA-seq) and without PMA treatment (traditional sequencing), we were able to explore how closely the microbiome compositions obtained from traditional sequencing methods resembled the viable microbiome composition obtained by PMA-seq (Figure 3A). These experiments established that at each skin site sampled, as compared to traditional sequencing, the PMA-treated samples were less rich (richness, R, is a measure of the number of identifiable bacterial taxa) and less diverse (diversity, H, is measured by the Shannon diversity index) (Figure 3B, C). Furthermore, samples that had greater richness in traditional sequencing (R_trad) showed proportionally larger decreases in richness and Shannon diversity with PMA-seq (R_PMA and H_PMA) (Figure 3B, C). These results suggest that, although it appears by traditional sequencing that there is a wide range of richness values at different skin sites (1-30 different taxa), in reality the richness across the skin microbiome at different body sites is relatively similar and low (1-10 different taxa). Thus, though there appears to be a wide range of diversity in the skin microbiome by traditional sequencing, PMA-seq indicates that this is generally an overestimation at any given skin site. In fact, our results suggest that the viable skin microbiome tends to be dominated by a relatively small number of taxa at most sites. An independent analysis of glabella samples from 5 healthy volunteers including contamination controls revealed the same trends that traditional sequencing overrepresented both diversity and richness with or without contaminant removal (Figure 6 – figure supplement 3).

The metrics of richness and diversity offer important information regarding how the composition of the skin microbiome changes at different body sites between traditional sequencing and PMA-seq. We further wanted to understand the changes in the relative abundance of specific bacterial taxa. To quantify taxon-level PMA-dependent changes, we developed a PMA-index (I_PMA) for each bacterial taxon, which is calculated as follows: , where A_PMA is the relative abundance by PMA-seq and A_trad is the relative abundance by traditional sequencing. A low PMA-index (below 0.5) indicates that the taxon in question is overrepresented by traditional sequencing, while a high PMA index (above 0.5) indicates that the taxon in question is underrepresented by traditional sequencing (Figure 4A). We note that these values represent enrichment relative to the rest of the sequences, a measurement that is distinct from the viability score. Calculating PMA-indices revealed that the abundances of most bacterial taxa at any given body site are overestimated by traditional sequencing, as most taxa had PMA-index values close to 0 (Figure 4A). 16S rRNA copy numbers could affect these values but the 16S copy numbers of the most abundant species were all within a roughly two-fold range, such that this effect is minor compared to the trends observed (21).

The four most abundant bacterial families (Propionibacteriaceae, Staphylococcaceae, Corynebacteriaceae, and Micrococcaceae) made up 93% of total sequencing reads (96% of PMA-seq reads and 91% of traditional sequencing reads) and demonstrated interesting family-level PMA-index patterns. The family Propionibacteriaceae includes a major component of the skin microbiome, C. acnes, which has been shown by traditional sequencing to comprise upwards of 50% of the skin microbiome irrespective of site type (8). PMA-seq revealed that traditional sequencing accurately represents Propionibacteriaceae abundance in sebaceous sites (demonstrated by a PMA-index close to 0.5), but overrepresents Propionibacteriaceae in moist and dry sites (PMA-indices of 0.2-0.3). Furthermore, Propionibacteriaceae dominated sebaceous sites (accounting for >75% of all viable bacteria in most sebaceous samples), but did not dominate moist or dry sites (their viable abundance did not exceed 50% of all viable bacteria in any of those samples) (Figure 4 A, B). Similar to the Propionibacteriaceae, bacteria in the family Staphylococcaceae appeared to be well-represented across all body sites (Figure 4A, B and Figure 4 – figure supplement 2).

Bacteria in the family Corynebacteriaceae are also considered main constituents of the skin microbiome, but our results showed that traditional sequencing overestimates the abundance of Corynebacteriaceae at every skin site except for the nares. For example, traditional sequencing identified a high abundance of Corynebacteriaceae in the popliteal fossa, but PMA-seq showed that these reads were largely of inviable origin (Figure 4B). Previous studies have demonstrated that Corynebacteria are readily cultured from nasal isolates, which supports our PMA-seq finding that viable members of this taxon are abundant in the nares but not at most other skin sites (22).

Interestingly, Micrococcaceae were overrepresented by traditional sequencing at every site except for the hair shaft. In the hair shaft, Micrococcaceae were abundant by PMA-seq but almost undetectable by traditional sequencing. As shown in Figure 4B, the increase in the relative proportion of viable Micrococcaceae detected by PMA-seq corresponds to a decrease in viable Propionibacteriaceae, suggesting that Micrococcaceae may not be detected by traditional sequencing because of the high abundance of DNA from inviable Propionibacteriaceae. These results suggest that most skin sites are colonized by a relatively small number of bacterial families, that different families distinctly colonize different skin sites, and that the majority of the additional bacterial DNA on the skin surface comes from inviable bacteria.

While our human skin microbiome FISH and PMA studies strongly agree, we are unable to directly compare the two approaches on the same samples. For a system in which we could perform both FISH and PMA-ddPCR on the same samples we turned to mouse skin, which also enabled us to ask if our findings are human-specific. We assessed the spatial distribution of bacterial cells in mouse skin tissue using the universal bacterial EUB338 FISH probe with tissue from K14-H2B-GFP mice. We observed the same bacterial distributions as seen in the human tissues: a high abundance of bacteria in hair follicles (enrichment score of 15.26) with relatively few bacteria on the skin surface (enrichment score of 0.21) (Figure 5A, D). To test how fur impacts the presence of bacteria on the skin surface, we performed FISH staining on skin from nude mice (SKH1-Hrhr Elite) and found similar bacterial distributions (follicle-associated enrichment score of 10.79 compared to 1.13 for the skin surface) (Figure 5B, D). The significant numbers of bacteria observed in hair follicles suggests that the absence of surface bacteria is not merely due to the relatively clean conditions in which laboratory mice are housed. Because the mouse tissue was untreated with any sterilizing agents prior to dissection, these results also suggest that the low numbers of skin-surface-associated bacteria found in the human tissue samples are not simply a result of pre-biopsy sterilization.

Since we did not detect many bacteria on the surface of any of the skin samples tested, we sought a positive control to confirm that our FISH staining can visualize bacteria on the skin surface if they are present. For example, it is possible that the dry, acidic nature of the skin affects FISH efficacy, leading to little staining at the skin surface. To this end, we applied E. coli cells to dorsal mouse skin tissue after removing it from the animal. This tissue was then processed in the same way as the human and other mouse tissue. FISH staining revealed many bacteria on the surface of these samples, confirming that this technique can be used to reliably visualize bacteria on the skin surface (Figure 5C). As a negative control, we also confirmed that a probe encoding the reverse complement of the EUB338 FISH probe (NONEUB338) did not significantly hybridize to the skin surface or follicles (Figure 5C). These results support the conclusions from our FISH experiments on human tissues by demonstrating EUB338 staining works well in the biological context of the skin surface.

We next assessed the viability of bacteria in the mouse skin microbiome using PMA-ddPCR. As further evidence that our previous findings are not human-specific, the PMA-ddPCR-based viability score for mouse skin microbiome sites was similar to the average viability score for human skin sites (0.066 and 0.045 respectively) and was much lower than the viability score for the mouse or human fecal microbiome (0.98 and 0.66 respectively). (Figure 5E). These results indicate that, despite having distinct skin biology, both humans and mice have an abundance of bacterial DNA on the skin surface that is not associated with viable cells.

Finally, we sought to determine the functional significance of our findings in the context of skin microbiome surface perturbation and repopulation. Specifically, we performed a perturbation-recovery experiment in which the skin microbiome of the forehead from healthy volunteers was sampled at time T=0 to establish a baseline microbial community, sterilized with benzalkonium chloride, and then sampled 3, 24, and 48 hours later to monitor recovery (Figure 6A). At each sampling timepoint, half of the sample was treated with PMA prior to DNA isolation and subsequent 16S rRNA gene sequencing. Focusing on the PMA-untreated samples, we found that 4 of the 5 volunteers demonstrated significant reduction in the absolute amount of DNA upon sterilization (Figure 6B), but that the overall bacterial composition of the skin microbiome remained relatively stable (Figure 6A, all bacteria listed in Figure 6 – figure supplement 2.). It should be noted that, although the skin microbiome composition found in volunteer 4 is markedly different than the communities identified in the other volunteers, the species identified are still common constituents of a healthy skin microbiome. Additionally, volunteer 4 had a lower abundance of bacterial DNA overall, which might lead to the lower relative abundance of Propionibacteriaceae.

Our findings recapitulate the previously-reported paradox that the skin microbiome is both unstable at short timescales and stable at long timescales (11, 23). Importantly, analysis of the PMA-treated samples enabled us to understand the recovery dynamics. Specifically, we used Bray-Curtis dissimilarity to monitor the extent of population recovery (Figure 6C). We found that comparing samples with and without the use of PMA to the PMA-treated baseline community resulted in a consistent pattern across all individuals. Over the course of the experiment, the communities converged back to the PMA-treated baseline samples but showed no consistent patterns of similarity relative to the PMA-untreated baseline samples (Figure 6 – figure supplement 2A). This result supports the hypothesis that bacterial DNA on the skin surface can be readily removed and is continuously replenished by viable bacterial populations in protected reservoirs below the skin surface.

Examining repopulation dynamics more finely at the ASV level revealed that ASVs that were lost upon surface sterilization often re-appeared over time, consistent with the model that surface repopulation can be driven by subsurface populations. Specifically, we found that after surface sterilization, the skin is repopulated by the same ASVs as were present before sterilization (Figure 6 – figure supplement 2B). In Figure 6 – figure supplement 2B we highlight specific examples in which a given ASV appears in the live-cell population (+PMA) prior to surface sterilization (T0), disappears from the live-cell population following sterilization (T3), and later re-appears in the live-cell population either 24 or 48 hours later. Together, these results suggest that bacteria thrive in protected areas below the skin surface and that the DNA of their dead remains accumulates on the skin surface. In this way, the bacterial DNA on the skin surface acts as a fingerprint of the communities below.

Human tissue was obtained through the Skin Translational Research Core within the Skin Biology and Diseases Resource-based Center (SBDRC) at the University of Pennsylvania. Normal adult human skin was obtained during skin cancer surgery procedures. After the cancerous tissue is removed, normal (non-cancerous) tissue must also be removed to ensure a cosmetic repair. These tissue portions are normally discarded, but were collected for imaging. Surgical scrub solution was used prior to the procedure (alcohol/iodine or chlorhexidine). Tissues were embedded in room temperature OCT immediately upon excision and stored at - 80°C. Tissues were obtained from the forehead of a 32-year-old female and the cheek of an 81- year-old male. Tissues were obtained through the Core B-Skin Translational Research Core at the SBDRC as part of the “Fresh normal and diseased human skin” service.

All mice were housed and maintained in a certified animal facility and all experiments were conducted according to USA Public Health Service Policy of Humane Care and Use of Laboratory Animals. All protocols were approved by the Institutional Animal Care and Use Committee, protocol #1867-17 (Princeton University). Dorsal skin from K14-H2B-GFP, C57BL/6, and SKH1-Elite nude mice was used for fluorescence experiments. The fur from K14- H2B-GFP and C57BL/6 mice was shaved using clippers prior to dissection. Dorsal tissue was removed, cut into thin (∼5mm) strips, and embedded in room temperature OCT. OCT blocks were placed on dry ice to freeze and subsequently stored at −80°C. Both sexes were used. All animals were between 3 and 6 weeks in age and were fed standard, non-sterilized rodent chow. No pathogens were identified in the regular health screenings conducted by facility staff. Mice were housed with 1-5 animals per cage with bedding and enrichment using Ventilated Thoren racks. Mice were euthanized using CO₂ inhalation (as per the AVMA Guidelines for the Euthanasia of Animals: 2020 Edition), followed by cervical dislocation. Dissections were performed with ethanol-sterilized tools.

Human and mouse tissues were processed identically. 30 µm tissue sections were sectioned and mounted on slides using a Leica CM3050S cryostat. Tissues on slides were fixed in 4% formaldehyde in 1X PBS for 20 mins. Following fixation, tissues were washed for 5 mins in 1X PBS and then incubated in hybridization buffer (0.9 M NaCl, 20 mM Tris-HCl, 10% formamide) containing 10 mg/mL lysozyme and FISH probes at 1 µg/µL for 2-3 hours at 47°C in the dark. Nucleotide sequences of FISH probes can be found in Table S1. After hybridization, slides were incubated in wash buffer (0.9 M NaCl, 20 mM Tris-HCl) for 1 hr at 47°C in the dark. Slides were then washed in 1X PBS for 10 mins. To visualize tissue nuclei, tissues were then stained with 1 µg/mL DAPI for 10 mins at room temperature. Slides were then washed 3X in 1X PBS for 10 mins each. Following the final wash step, tissues were mounted with glycerol-based anti-fade non-curing mounting media. Coverslips were sealed using a 1:1:1 mixture of petroleum jelly, lanolin, and paraffin. Images were acquired on a Nikon A1R-Si HD confocal microscope controlled by NIS Elements software. ImageJ and MATLAB (Mathworks, Natwick, MA) were used for image processing. To calculate fluorescence enrichment scores, maximum projections of Z-stacks were used. A rectangular region of interest (ROI) of width W was drawn around a follicle using only the DAPI channel. A second ROI was drawn near the follicle opening with dimensions 0.5W x 2W in order to capture follicle-associated fluorescence near the follicle opening. In an orthogonal processing step, the pixel intensities for the entire image in the FISH fluorescent channel were fit using a 3 component Gaussian mixture model (GMM) corresponding roughly to non-tissue background, non-probe autofluorescence, and probe-based signal. Using this GMM as the background subtraction value for each image, the mean intensity for each ROI was calculated. This threshold value was used to calculate the mean fluorescence value inside the ROI and outside the ROI. The ratio of the mean fluorescence value inside the ROI to the mean fluorescence value outside of the ROI was used for quantifying follicle-and stratum corneum-associated fluorescence (enrichment score). Significance was calculated using the Mann-Whitney test.

Stationary phase cultures were grown according to the following conditions:

Microbiome samples were collected under Princeton University IRB #13003 at the Princeton University Department of Molecular Biology. Healthy volunteers were recruited via informational flyers. Volunteers gave informed consent prior to sample collection. Participants were healthy volunteers, male and female, white, aged 26-35, with no history of chronic skin conditions or autoimmune disease and were not using antibiotics at the time of sampling or one month prior to sample collection.

Skin microbiomes from healthy volunteers were collected using sterile foam-tipped collection swabs pre-moistened with sterile 1X DPBS. Though often included, we chose not to use detergent in the swabbing buffer in order to avoid negatively affecting bacterial cell membranes and altering viability scores. Figure 2 – figure supplement 1G shows how swabbing with and without 0.1% Triton X-100 (Sigma) affects viability scores, ddPCR results, and CFU. Areas of interest were sampled for 60 seconds before being re-suspended in sterile 1X DPBS. Tongue microbiome samples were collected using sterile foam-tipped collection swabs. Hair shaft samples were collected by plucking hairs and using only the bulb portion. Follicle contents were collected using Bioré Deep Cleansing Pore Strips (Kao USA Inc., Cincinnati, OH). Saliva was collected in sterile 50 mL conicals from healthy volunteers. Plaque was collected by scraping the teeth of healthy volunteers using sterile toothpicks and re-suspending the collection in sterile 1X DPBS. Murine fecal samples from C57BL/6 mice were collected during dissection. Human skin microbiome samples shown in Figure 2 – figure supplement 1 were plated for CFU calculations prior to the addition of PMA. Samples were plated on blood agar plates (5% sheep blood in tryptic soy agar, VWR International) and grown for 24-48 hours aerobically or anaerobically. For perturbation recovery experiments, volunteers were asked to refrain from using any type of products on their face for the duration of the experiment (48 hours) as well as 24 hours prior to baseline sampling. Volunteers were allowed to bathe/shower, but were asked to refrain from washing their faces. Skin surface sterilization was done by scrubbing the sample area for 60 seconds with an antiseptic wipe containing benzalkonium chloride 0.13% w/v (North by Honeywell, #150910). To ensure that the process of sampling the skin surface would not affect the bacterial repopulation dynamics, forehead of each volunteer into four sections of equal areas: T0 (baseline, before sterilization), T1 (3hrs post sterilization), T2 (24 hrs post-sterilization), and T3 (48 hrs post-sterilization). Each section was sampled for 60 seconds using a sterile foam-tipped collection swab pre-moistened with sterile 1X DPBS. As described above, each swab was then resuspended in sterile 1X DPBS.

The human fecal sample was collected by the Donia lab under Princeton University IRB #11606 and was gifted to the Gitai lab. This sample and the methods of collection were described in Javdan, Bahar et al. “Personalized Mapping of Drug Metabolism by the Human Gut Microbiome.” Cell vol. 181,7 (2020): 1661-1679.e22. doi:10.1016/j.cell.2020.05.001 (27).

In order to demonstrate the efficacy of PMA, known ratios of live and heat-killed E. coli cultures were mixed and subjected to PMA treatment. First, an overnight culture of E. coli NCM3722 was back-diluted into fresh LB media at a ratio of 1:1000 and grown at 37°C for 4 hours to reach mid-log. Stationary-phase experiments were done with overnight, 18h cultures. Cultures were then washed 3X and suspended in sterile PBS. Half of the culture was heat-killed by incubating at 70°C for 20 minutes while the other half remained at room temperature. The heat-killed E. coli cultures were allowed to cool to room temperature before combining with non-heat-killed E. coli cultures to achieve a 50% (by volume) heat-killed mixture. For the 0% heat-killed, no heat-killed bacteria were added. Likewise, for the 100% heat-killed, only heat-killed bacteria were used.

Each condition was mixed well and then split evenly between two sterile 1.5 mL microcentrifuge tubes.

After collection, samples were split evenly between two sterile 1.5 mL microcentrifuge tubes. PMA (Biotium Inc.) was added to one of the two tubes to a final concentration of 50 µM. All tubes were incubated in the dark at room temperature for 10 mins before being exposed to light to cross-link PMA molecules using the PMA-Lite^TM LED Photolysis Device (Biotium Inc.). DNA was then isolated from all samples using either the DNeasy PowerSoil Kit (Qiagen) or according to the protocol outlined in Meisel et al. (28). The Meisel et al. protocol was modified such that sonication (Branson Sonifier 250, VWR Scientific) was used instead of bead-beating for the purposes of cell-lysis. Both DNA isolation methods work similarly (Figure 2 – figure supplement 1H). If lysostaphin (Sigma-Aldrich) was used, it was added following PMA activation and before DNA isolation to a final concentration of 0.1 mg/mL and incubated at room temperature for 30 mins.

The Bio Rad QX200 AutoDG Droplet Digital PCR System was used to quantify extracted DNA from microbiome samples and from pure bacterial cultures. Reaction mixtures contained 2x QX200 ddPCR EvaGreen Supermix and universal 16S qPCR primers at 10 nM concentrations in a total volume of 25 µL. Primer sequences can be found in Table S1. Reaction mixtures were transferred to sterile ddPCR 96-well plates (BioRad #12001925) which were loaded into the QX200 Automated Droplet Generator. After droplet generation, the plate was heat-sealed using the PX1 PCR Plate Sealer (BioRad #1814000) and PCR was performed with a pre-step of 95 C for 5 minutes followed by 40 rounds of amplification with 60C, 1 minute extensions and a final hold temperature of 12 C using a C1000 Touch Thermal Cycler (BioRad #1851197). Samples were subsequently loaded into the QX200 Droplet Reader for quantification. Automatic thresholding was performed using the Quantasoft software and subsequently exported to Microsoft Excel for analysis. Significance was calculated using a Student’s T-Test. To calculate the viability score for a given pair of “-PMA” and “+PMA” matched samples, the following calculation was done:

Cultures of S. epidermidis EGM 2-06 were grown in tryptic soy broth (TSB) overnight and diluted 1:1000 the following morning in TSB and grown for 4 hours until a final OD of 0.4. Tenfold dilutions of S. epidermidis culture were then prepared, plated for CFU on 5% sheep blood in tryptic soy agar (VWR International) and divided between two 1.5 mL microcentrifuge tubes. Stationary-phase cultures of C. acnes, S. epidermidis, M. luteus, and C. striatum were grown according to the table above. PMA was added to one tube for a final concentration of 50 µM and the other tube was left untreated. PMA activation and DNA isolation was then done according to the methods outlined above.

DNA was isolated from microbiome samples (with and without PMA) using the DNeasy PowerSoil Kit (Qiagen) (Figures 2-5) or according to the protocol outlined in Meisel et al (Figure 6). The V1-V3 region of the 16S gene was amplified using the primers 27F (5’- AGAGTTTGATCCTGGCTCAG) and 534R (5’- ATTACCGCGGCTGCTGG). Illumina sequencing libraries were prepared using previously published primers (29). Libraries were then pooled at equimolar ratios and sequenced on an Illumina MiSeq Micro 500 nt (Figure 3, Figure 4) or MiSeq V3 600 nt (Figure 6) as paired-end reads. Reads were 2X250 bp or 2X300 bp with an average depth of ∼33,616 reads for 500 nt and 245,488 reads for 600 nt (Tables S2, S3). Also included were 8 bp Index reads, following the manufacturer’s protocol (Illumina, USA). Raw sequencing reads were filtered by Illumina HiSeq Control Software to generate Pass-Filter reads for further analysis. Index reads were used for sample de-multiplexing. Amplicon sequencing variants (ASVs) were then inferred from the sequences using the DADA2 plugin within QIIME2 version 2018.6 (30, 31). Reverse reads were trimmed to 245 bp for the 2X250 bp data, and the reverse reads were trimmed to 275 bp for the 2X300 bp data. For the 2X250 bp data, two different ASV inference methods were utilized: one using only the forward reads (to capture the longer-than-average Staphylococcus amplicons) and one using both paired-end reads. For the 2X300 bp data, ASV inference was performed using both paired-end reads. ASV inference using forward reads was performed using DADA2’s denoise-single function, while ASV inference using both reads was performed using the denoise-paired function.. Taxonomy was assigned to the resulting ASVs with a naive Bayes classifier trained on the GreenGenes database version using only the regions of the 16S rRNA gene spanned by the ASVs (32, 33). All downstream analyses were performed using family-level taxonomy assignments unless specified to be at the ASV-level. PBS control samples were included at the library preparation step for 500 nt sequencing data and air swab DNA isolation controls were included for 600 nt sequencing data. Contamination removal was performed on 600 nt sequencing data using the Decontam package for R (Figure 6 and Figure 6 – figure supplement 2 show data without contamination removal). (34). Figure 6 – figure supplement 3 shows how Shannon diversity and richness change between traditional sequencing and PMA-seq (similar to what is shown in Figure 3B and Figure 3C).

Removal of contaminants (using either 0.1 or 0.2 threshold) has little effect on the data, suggesting that, although contamination removal cannot be performed on the 500 nt sequencing data due to the absence of air swab controls, this should not change the overall findings shown in Figure 3 and Figure 4. Relative abundance, richness, Shannon diversity, Bray-Curtis dissimilarity, and PMA-index were assessed using the Vegan package for R or Microsoft Excel and plotted using R, Prism, and MATLAB (35). The PMA-index was calculated using relative abundance and was not calculated for any bacterial taxa that was present in fewer than four samples.

Sequencing datasets generated and analyzed during this study are available at NCBI Sequence Read Archive (SRA) (BioProject numbers (PRINA918671 and PRINA918959). Oligonucleotide sequences used are included in table S1.