Bacterial communities on classroom surfaces vary with human contact.

BACKGROUND: Humans can spend the majority of their time indoors, but little is known about the interactions between the human and built-environment microbiomes or the forces that drive microbial community assembly in the built environment. We sampled 16S rRNA genes from four different surface types throughout a university classroom to determine whether bacterial assemblages on each surface were best predicted by routine human interactions or by proximity to other surfaces within the classroom. We then analyzed our data with publicly-available datasets representing potential source environments.
RESULTS: Bacterial assemblages from the four surface types, as well as individual taxa, were indicative of different source pools related to the type of human contact each surface routinely encounters. Spatial proximity to other surfaces in the classroom did not predict community composition.
CONCLUSIONS: Our results indicate that human-associated microbial communities can be transferred to indoor surfaces following contact, and that such transmission is possible even when contact is indirect, but that proximity to other surfaces in the classroom does not influence community composition.

Background

In the developed world, humans spend a majority of their lives indoors. While indoors we encounter microorganisms on virtually every surface we touch, and this frequent exposure to indoor microbes carries with it the potential for disease transmission, as well as interactions with our own commensal microbiome [1-3]. Yet we have very little knowledge regarding the ecological processes that drive microbial community assembly indoors, nor do we understand the degree to which humans share microbes with indoor surfaces.

For any given indoor surface, microbial community composition is likely shaped by habitat-specific environmental constraints (such as the type of surface material), and dispersal sources (which include humans, bioaerosols and other surfaces within a space). If dispersal among surfaces is a primary determinant of community structure, then adjacent surfaces should be more similar in community composition than surfaces further apart. If dispersal from humans is a major determinant, then community structure should vary with frequency and nature of human contact.

Because humans harbor distinguishably different microbial assemblages on different parts of their bodies [4-6], it is reasonable to assume that different indoor surface types could harbor different microbial communities due to frequent contact with specific body parts. Several studies suggest that this may be the case. In a recent survey of public restroom surfaces, Flores and colleagues [7] found that microbes on restroom surfaces were similar to those found on specific human body parts, with the strongest association observed between toilet surfaces and gut and vaginal communities. These associations are to be expected given the direct contact common in this environment between surfaces and the human body. Associations between human use and bacterial community composition have also been found on residential kitchen surfaces [8], with bacterial taxa commonly found on human skin predominating on kitchen surfaces, consistent with frequent skin-to-surface contact. It is not known whether more indirect skin-to-surface contact, such as sitting fully clothed on a chair, might result in similar bacterial transmission.

Here we test: (a) if different types of surfaces in a classroom vary in microbial community composition; (b) if their composition varies predictably with the type of human contact; and (c) whether these associations are stronger than the association between community similarity and spatial distance. We sampled desktops, floors, walls and chair seats in a classroom to determine whether these groups of surfaces harbored distinct bacterial assemblages related to their differing types of human contact.

Methods

All samples were collected in the same classroom at the Lillis Business Complex, University of Oregon, Eugene, OR, USA, during 4–5 August 2011. During the week prior to the sampling period, the classroom was occupied by classes daily, including on the first day of sampling. Desks, floors, chair seats and walls were sampled (n = 18, 18, 18 and 16, respectively) following identical protocols. Samples were distributed throughout the classroom (9.2 m × 15.2 m). Pairwise spatial distance between samples ranged from 0 to 16 m, with sets of chair, desk and floor samples generally taken within 1 m2 of one another. Figure 1 displays the spatial distribution of samples within the classroom. Surfaces were sampled using a nylon flocked swab (copanusa.com; #552C) moistened with sterile buffer solution (0.15 M NaCl, 0.1 Tween20). A 289 cm2 (17 cm × 17 cm) area was swabbed from each surface. Chairs were sampled in the center of their upholstered seat surfaces; desks (plastic laminate surface) and floors (linoleum) were sampled directly above and below chair seats; and walls (latex paint) were sampled around the perimeter of the room (Figure 1). We were not able to ascertain the cleaning schedule prior to sampling, but all surfaces were visibly devoid of standing dust. All samples were frozen (−80°C) until DNA extraction.

Figure 1

Schematic of sampling design. Four different types of surfaces (desks, chairs, floors and walls) were sampled throughout an amphitheater-style classroom.

DNA extraction, amplification and Illumina library preparation followed methods described previously [9,10]. DNA was extracted from swabs using a PowerWater DNA extraction kit (MoBio Laboratories, Inc., Carlsbad, CA, USA) with the following modifications: samples were frozen and thawed for two cycles; bead beating length was extended to 10 minutes; and samples were eluted in 50 μL Solution PW6.

The V4 region of the 16S rRNA gene was amplified using the F515/R806 primer combination (5′-GTGCCAGCMGCCGCGG-3′, 5′-TACNVGGGTATCTAATCC-3′; [11]). Amplification proceeded in two steps using a custom Illumina preparation protocol where PCR1 was performed with forward primers that contained partial unique barcodes and partial Illumina adapters. The remaining ends of the Illumina adapters were attached during PCR2, and barcodes were recombined in silico using paired-end reads. Adapter sequences are detailed in Additional file 1. All extracted samples were amplified in triplicate for PCR1 and triplicates were pooled before PCR2. PCR1 (25 μL total volume per reaction) consisted of the following steps: 5 μL 5× HF buffer (Thermo Fisher Scientific, Waltham, MA, USA), 0.5 μL deoxyribonucleotide triphosphates (10 mM, Invitrogen, Life Technologies, Grand Island, NY, USA), 0.25 μL Phusion Hotstart II polymerase (0.5 units; Thermo Fisher Scientific), 13.25 μL certified nucleic-acid free water, 0.5 μL (10 μM) forward primer, 0.5 μL (10 μM) reverse primer, and 5 μL template DNA. The PCR1 conditions were as follows: initial denaturation for 2 minutes at 98°C; 22 cycles of 20 seconds at 98°C, 30 seconds at 50°C and 20 seconds at 72°C; and 72°C for 2 minutes for final extension. After PCR1, the triplicate reactions were pooled and cleaned with the Qiagen MinElute PCR Purification Kit according to the manufacturer’s protocol (Qiagen, Germantown, MD, USA). Samples were eluted in 11.5 μL Buffer EB (10 mM Tris-Cl, pH 8.5, Qiagen). For PCR2, a single primer pair was used to add the remaining Illumina adapter segments to the ends of the concentrated amplicons of PCR1. The PCR2 (25 μL volume per reaction) consisted of the same combination of reagents that was used in PCR1, along with 5 μL concentrated PCR1 product as a template. The PCR2 conditions were as follows: 2 minutes denaturation at 98°C; 12 cycles of 20 seconds at 98°C, 30 seconds at 66°C and 20 seconds at 72°C; and 2 minutes at 72°C for final extension.

Amplicons were size selected by gel electrophoresis, extracted, quantified on a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), concentrated, combined in equimolar concentrations and sequenced on the Illumina MiSeq platform at the Dana-Farber/Harvard Cancer Center DNA Resource Core (Boston, MA; dnaseq.med.harvard.edu).

Sequence processing was performed using the FastX Toolkit (http://hannonlab.cshl.edu/fastx_toolkit) and QIIME [12]. Quality filtering settings were as follows: a minimum 30 quality score over at least 75% of the sequence read, no ambiguous bases, and 1 primer mismatch allowed. After quality control and barcode assignment, the remaining high-quality sequences (mean 13,786 sequences per sample ± 4,735 SD) were binned into operational taxonomic units (OTUs) at a 97% sequence similarity cutoff using uclust [13]. High-quality sequences from each OTU cluster were taxonomically identified using reference sequences from the Greengenes database (2011 release; [14,15]), and plant-chloroplast OTUs were removed by name (“Streptophyta”) based on Greengenes taxonomic classifications (mean = 8% of sequences in each sample ± 10% SD). All samples retaining more than 4,000 sequences (n = 15, 14, 15 and 14 for desks, walls, floors, and chairs, respectively) were rarefied to that level for even sampling depth. Sequence files and metadata for all samples used in this study have been deposited in Figshare (http://dx.doi.org/10.6084/m9.figshare.687155). Metadata, the unrarefied OTU table, and corresponding taxonomic classifications have all been included as Additional files 2, 3 and 4, respectively.

All statistical analyses were performed in R [16], primarily utilizing multivariate community ecology procedures in the labdsv and vegan packages [17,18]. A full record of all statistical analysis is included as Additional file 5, and was created using the knitr package in R [19]. Community distance between samples was calculated using the Canberra distance, implemented in the vegan package. Given that the four surface types we sampled represent potentially similar microbial habitats, we did not expect to see major phylogenetic differences among the communities. Thus, we used the Canberra taxonomic metric as a way to emphasize community differences driven by relatively rare taxa with a decreased emphasis on total abundance and phylogenetic community distance. Distance-based redundancy analysis (DB-RDA) was performed with the capscale function in vegan. Individual OTUs exert varying degrees of influence, or weight, on DB-RDA clustering patterns, and we selectively tested the strength of association between the strongest weighted individual bacterial OTUs and relative surface types (that is, indicator value, as defined in [20]) with the indval function in labdsv. Mantel tests, using the mantel function in vegan, were used to test for a community dissimilarity correlation with spatial distance within the classroom for each of the four surface types, as well as all sites regardless of surface type. To assess community similarity to various potential source habitats (soil, aquatic, phyllosphere, human skin and human gut), we compared all samples at the class level to publicly-available datasets from representative potential source environments [4,21-24] using principal components analysis. Positions of each point in relation to source environments along the first principal component are used to compare their similarity to either human skin or phyllosphere bacterial communities.

Results and discussion

After rarefaction to 4,000 sequences per sample, 58 surface samples were represented by a total of 3,745 bacterial OTUs. Approximately half (51.4%) of OTUs were found only once or twice, and these were included in analysis. The most commonly detected OTU (Sphingomonas sp.) accounted for 1.9% of all sequences. The four surface types harbored significantly different communities (P = 0.001; F = 2.34; from permutational multivariate analysis of variation on Canberra distances; Figure 2). Within each surface type, bacterial community composition was not significantly predicted by spatial distance among samples (all four P values >0.1; from Mantel tests). Even when all surfaces were considered together, spatial distance was not a significant predictor of bacterial community composition (P = 0.4). These results suggest that site-specific factors (for example, habitat selection or dispersal from specific body sites) are more important than dispersal among sites for bacterial community assembly in the built environment.

Figure 2

Surfaces harbored significantly different bacterial communities and were linked to differential human contact. (a) Bacterial communities were constrained by four different surface types using distance-based redundancy analysis (DB-RDA; constrained inertia = 11.4%) and were significantly different among types based on Canberra taxonomic distances (P = 0.001 from permutational multivariate analysis of variation). (b) Bacterial operational taxonomic units (OTUs) from DB-RDA are shown weighting communities in the same four primary directions. The first and second axes from DB-RDA are used in both ordinations (CAP 1 and CAP 2) . The strongest ten weighting OTUs for each surface type are highlighted if they were also significant indicator OTUs (all P values <0.05). (c) All samples were compared to potential source environments using principal components analysis (PCA), and the first principal component (37.8% of variance explained) was used as a surrogate for community similarity to either phyllosphere or human skin bacterial communities. Boxplots delineate (from bottom) minimum value, Q1, median (Q2), Q3, maximum value; notches approximate 95% confidence around median value, and outliers fall outside of the quartile range. Letters above each box indicate significant groupings after Tukey’s hones significant difference (HSD) test (adjusted P value <0.05).

We found that surface communities were significantly associated with taxa specific to distinct microbial sources (Table 1). Bacterial OTUs labeled in Figure 2b are the most strongly weighted OTUs for each surface type that are also significant (P < 0.05) indicator OTUs. Three of the strongest indicators for chair seats were Lactobacillus OTUs closely related to taxa commonly found in the human gut and vagina, in addition to other human-associated Corynebacterium and Staphylococcus species. Desk surfaces were significantly associated with several indicator taxa from human body habitats; two of the strongest desk surface indicator OTUs were Streptococcus species commonly found in human skin and oral samples, as well as another Streptococcus originally isolated from purulent infections. Floors in the classroom also harbored skin-associated OTUs, but were more strongly indicated by a cyanobacterial OTU, presumably from a non-human environmental source. Walls are likely the surfaces in this study that have the least contact with humans; this was consistent with our observation that these surfaces were associated with indicator taxa related to Sphingomonas and Alicyclobacillus species that are commonly abundant in airborne bacterial assemblages [10,25]. Both walls and floors held larger relative proportions of Cyanobacteria than seats and desktops, likely reflecting soil and bioaerosols as microbial sources. Few other major phylogenetic differences are evident at a broad taxonomic level (Figure 3).

Table 1

Closest known isolates related to indicator operational taxonomic units

Greengenes genus	P value	Surface type	Closest 16S NCBI isolate and accession	Isolate source environment	Sequence similarity to isolate (%)
Lactobacillus	0.001*	Chairs	Lactobacillus johnsonii NR_075064.1	Human gut	99
Corynebacterium	0.001*	Chairs	Corynebacterium resistens NR_040999.1	Human infection	99
Corynebacterium	0.001*	Chairs	Corynebacterium confusum NR_026449.1	Human clinical specimens	99
Staphylococcus	0.011*	Chairs	Staphylococcus epidermidis NR_074995.1	Human skin	99
Corynebacterium	0.001*	Chairs	Corynebacterium riegelii NR_026434.1	Human urinary tract	99
Staphylococcus	0.019*	Chairs	Staphylococcus saprophyticus NR_074999.1	Human urinary tract	99
Lactobacillus	0.001*	Chairs	Lactobacillus crispatus NR_074986.1	Human vagina	99
Lactobacillus	0.003*	Chairs	Lactobacillus acidophilus NR_075049.1	Human gut	99
Streptococcus	0.001*	Desks	Streptococcus oralis NR_102809.1	Human oral	99
Streptococcus	0.001*	Desks	Streptococcus salivarius NR_102816.1	Human oral	99
Brevundimonas	0.002*	Desks	Brevundimonas variabilis NR_037106.1	Pond water	99
Streptococcus	0.001*	Desks	Streptococcus intermedius NR_102797.1	Human purulent infection	99
CandidatusPhytoplasma	0.001*	Desks	None**	-	-
Alicyclobacillus	0.001*	Walls	Tumebacillus permanentifrigoris NR_043849.1	Soil	99
Chroococcidiopsis	0.028*	Walls	Halospirulina tapeticola NR_026510.1	Saline aquatic	96
Alicyclobacillus	0.001*	Walls	Tumebacillus permanentifrigoris NR_043849.1	Soil	98
Rhodopseudomonas	0.001*	Walls	Methylobacterium adhaesivum NR_042409	Drinking water	98
Salmonella	0.001*	Floors	Pantoea ananatis NR_103927.1	Phyllosphere	99
Roseomonas	0.001*	Floors	Roseomonas gilardii NR_029061.1	Human blood	99
Roseomonas	0.001*	Floors	Roseomonas frigidaquae NR_044455.1	Water-cooling system	99
Salmonella	0.001*	Floors	Pantoea ananatis NR_103927.1	Phyllosphere	99

All extant operational taxonomic units labeled in Figure 2 (and thus influential in distance-based redundancy analysis, as well as significant indicator taxa for their respective surface type) were related to their closest known bacterial isolate using 16S rRNA sequences in the NCBI Bacteria & Archaea Isolate Database. Source environments are from each isolate’s respective published source environment. *Unadjusted P value < 0.05. **Closest known isolate 89% similar. NCBI: National Center for Biotechnology Information.

Figure 3

Taxonomic composition of all 58 samples used in this study. Samples are grouped by surface type. All taxonomic groups representing <1% (Phylum and Class) and <4% (Order) of sequences were grouped into ‘other’.

The community variation we observed among surface types could be the result of either surface-specific environmental filtering or dispersal. Evidence to date suggests that microbial communities on indoor surfaces, with a few exceptions, are likely primarily dispersal-driven given the well-documented inputs from dust, settled airborne particles, outdoor air and surface contact [7-9,25-27]. The relative importance of in situ community dynamics and habitat filtering remain to be fully understood in the built environment.

Dispersal between human skin and contacted surfaces is well documented, primarily in a medical context using known pathogenic strains [1,28,29]. However, the vast majority of bacteria on and in the human body are not pathogens, but rather appear to be commensal components of our own microbiome [30]. Recent evidence has suggested that the unique microbial assemblages detected on different human body parts can be transferred to indoor surfaces following contact [26,31], and that these bacterial traces of human contact are evident in places such as restrooms [7]. Our results suggest that this transmission is possible even when that contact is indirect, and that such contact has a greater impact on microbial community structure indoors than dispersal among surfaces.

Conclusion

We are constantly surrounded by diverse microbial communities indoors, and we are just beginning to understand how our interactions in the built environment shape those communities and our own human microbiome. Our results indicate that human-associated microbial communities can be transferred to indoor surfaces following contact, and that such transmission is possible even when contact is indirect, but that proximity to other surfaces in the classroom does not influence community composition.

Availability of supporting data

Sequence files and metadata for all samples used in this study have been deposited in Figshare (http://dx.doi.org/10.6084/m9.figshare.687155). A full record of all statistical analysis is included as Additional file 5, and was created using the knitr package in R [19]. Original R scripts are available in GitHub (https://github.com/jfmeadow/Meadow_etal_Surfaces). Metadata, the unrarefied OTU table, and corresponding taxonomic classifications have all been included as Additional files 2, 3 and 4, respectively.

Abbreviations

DB-RBA: distance-based redundancy analysis; OTU: operational taxonomic unit; PCR: polymerase chain reaction.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JFM analyzed the data and wrote the manuscript. AEA performed laboratory assays. SWK conceived the study and collected samples. MM conceived the study, prepared figures and collected samples. TKO conceived the study, performed laboratory assays and collected samples. AMW conceived the study and collected samples. GZB contributed resources. JLG conceived the study, contributed resources and collected samples. BJMB conceived the study, contributed resources and collected samples. All authors edited the manuscript and approved the final draft.

Additional file 1

Explanation of partial adapter sequences used in Illumina amplicon library preparation.

Click here for file

Additional file 2

Metadata associated with all samples used in this study.

Click here for file

Additional file 3

Unrarefied operational taxonomic unit table used for this study.

Click here for file

Additional file 4

Taxonomic assignments for all operational taxonomic units in Additional file.

Click here for file

Additional file 5

Full account of statistical analysis performed in R.

Click here for file