16S Based Microbiome Analysis from Healthy Subjects’ Skin Swabs Stored for Different Storage Periods Reveal Phylum to Genus Level Changes

Microbiome research and improvements in high throughput sequencing technologies revolutionize our current scientific viewpoint. The human associated microbiome is a prominent focus of clinical research. Large cohort studies are often required to investigate the human microbiome composition and its changes in a multitude of human diseases. Reproducible analyses of large cohort samples require standardized protocols in study design, sampling, storage, processing, and data analysis. In particular, the effect of sample storage on actual results is critical for reproducibility. So far, the effect of storage conditions on the results of microbial analysis has been examined for only a few human biological materials (e.g., stool samples). There is a lack of data and information on appropriate storage conditions on other human derived samples, such as skin. Here, we analyzed skin swab samples collected from three different body locations (forearm, V of the chest and back) of eight healthy volunteers. The skin swabs were soaked in sterile buffer and total DNA was isolated after freezing at -80°C for 24 h, 90 or 365 days. Hypervariable regions V1-2 were amplified from total DNA and libraries were sequenced on an Illumina MiSeq desktop sequencer in paired end mode. Data were analyzed using Qiime 1.9.1. Summarizing all body locations per time point, we found no significant differences in alpha diversity and multivariate community analysis among the three time points. Considering body locations separately significant differences in the richness of forearm samples were found between d0 vs. d90 and d90 vs. d365. Significant differences in the relative abundance of major skin genera (Propionibacterium, Streptococcus, Bacteroides, Corynebacterium, and Staphylococcus) were detected in our samples in Bacteroides only among all time points in forearm samples and between d0 vs. d90 and d90 vs. d365 in V of the chest and back samples. Accordingly, significant differences were detected in the ratios of the main phyla Actinobacteria, Firmicutes, and Bacteroidetes: Actinobacteria vs. Bacteroidetes at d0 vs. d90 (p-value = 0.0234), at d0 vs. d365 (p-value = 0.0234) and d90 vs. d365 (p-value = 0.0234) in forearm samples and at d90 vs. d365 in V of the chest (p-value = 0.0234) and back samples (p-value = 0.0234). The ratios of Firmicutes vs. Bacteroidetes showed no significant changes in any of the body locations as well as the ratios of Actinobacteria vs. Firmicutes at any time point. Studies with larger sample sizes are required to verify our results and determine long term storage effects with regard to specific biological questions.