The Human Skin Double-Stranded DNA Virome: Topographical and Temporal Diversity, Genetic Enrichment, and Dynamic Associations with the Host Microbiome

ABSTRACT Viruses make up a major component of the human microbiota but are poorly understood in the skin, our primary barrier to the external environment. Viral communities have the potential to modulate states of cutaneous health and disease. Bacteriophages are known to influence the structure and function of microbial communities through predation and genetic exchange. Human viruses are associated with skin cancers and a multitude of cutaneous manifestations. Despite these important roles, little is known regarding the human skin virome and its interactions with the host microbiome. Here we evaluated the human cutaneous double-stranded DNA virome by metagenomic sequencing of DNA from purified virus-like particles (VLPs). In parallel, we employed metagenomic sequencing of the total skin microbiome to assess covariation and infer interactions with the virome. Samples were collected from 16 subjects at eight body sites over 1 month. In addition to the microenviroment, which is known to partition the bacterial and fungal microbiota, natural skin occlusion was strongly associated with skin virome community composition. Viral contigs were enriched for genes indicative of a temperate phage replication style and also maintained genes encoding potential antibiotic resistance and virulence factors. CRISPR spacers identified in the bacterial DNA sequences provided a record of phage predation and suggest a mechanism to explain spatial partitioning of skin phage communities. Finally, we modeled the structure of bacterial and phage communities together to reveal a complex microbial environment with a Corynebacterium hub. These results reveal the previously underappreciated diversity, encoded functions, and viral-microbial dynamic unique to the human skin virome. IMPORTANCE To date, most cutaneous microbiome studies have focused on bacterial and fungal communities. Skin viral communities and their relationships with their hosts remain poorly understood despite their potential to modulate states of cutaneous health and disease. Previous studies employing whole-metagenome sequencing without purification for virus-like particles (VLPs) have provided some insight into the viral component of the skin microbiome but have not completely characterized these communities or analyzed interactions with the host microbiome. Here we present an optimized virus purification technique and corresponding analysis tools for gaining novel insights into the skin virome, including viral “dark matter,” and its potential interactions with the host microbiome. The work presented here establishes a baseline of the healthy human skin virome and is a necessary foundation for future studies examining viral perturbations in skin health and disease. To date, most cutaneous microbiome studies have focused on bacterial and fungal communities. Skin viral communities and their relationships with their hosts remain poorly understood despite their potential to modulate states of cutaneous health and disease. Previous studies employing whole-metagenome sequencing without purification for virus-like particles (VLPs) have provided some insight into the viral component of the skin microbiome but have not completely characterized these communities or analyzed interactions with the host microbiome. Here we present an optimized virus purification technique and corresponding analysis tools for gaining novel insights into the skin virome, including viral “dark matter,” and its potential interactions with the host microbiome. The work presented here establishes a baseline of the healthy human skin virome and is a necessary foundation for future studies examining viral perturbations in skin health and disease.

[1]  Peter Salamon,et al.  Viral and microbial community dynamics in four aquatic environments , 2010, The ISME Journal.

[2]  Steven J. M. Jones,et al.  Circos: an information aesthetic for comparative genomics. , 2009, Genome research.

[3]  Mya Breitbart,et al.  Marine viruses: truth or dare. , 2012, Annual review of marine science.

[4]  F. Bushman,et al.  The human gut virome: inter-individual variation and dynamic response to diet. , 2011, Genome research.

[5]  S. Sarkar,et al.  The Simes Method for Multiple Hypothesis Testing with Positively Dependent Test Statistics , 1997 .

[6]  Michael G. Kemp,et al.  The histone deacetylase inhibitor trichostatin A alters the pattern of DNA replication origin activity in human cells , 2005, Nucleic acids research.

[7]  Fidel Ramírez,et al.  Computing topological parameters of biological networks , 2008, Bioinform..

[8]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[9]  J. Imhoff,et al.  Phenelfamycins G and H, new elfamycin-type antibiotics produced by Streptomyces albospinus Acta 3619 , 2011, The Journal of Antibiotics.

[10]  Frederic D Bushman,et al.  Rapid evolution of the human gut virome , 2013, Proceedings of the National Academy of Sciences.

[11]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[12]  C. Deming,et al.  Topographical and Temporal Diversity of the Human Skin Microbiome , 2009, Science.

[13]  Robert C. Edgar,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

[14]  P. Legendre,et al.  vegan : Community Ecology Package. R package version 1.8-5 , 2007 .

[15]  Timothy L. Tickle,et al.  Computational meta'omics for microbial community studies , 2013, Molecular systems biology.

[16]  William A. Walters,et al.  QIIME allows analysis of high-throughput community sequencing data , 2010, Nature Methods.

[17]  J. Tukey,et al.  Variations of Box Plots , 1978 .

[18]  J. Collins,et al.  Antibiotic Treatment Expands the Resistance Reservoir and Ecological Network of the Phage Metagenome , 2013, Nature.

[19]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[20]  Curtis Huttenhower,et al.  Microbial Co-occurrence Relationships in the Human Microbiome , 2012, PLoS Comput. Biol..

[21]  Shoshana J. Wodak,et al.  ACLAME: A CLAssification of Mobile genetic Elements , 2004, Nucleic Acids Res..

[22]  Jo Handelsman,et al.  A statistical toolbox for metagenomics: assessing functional diversity in microbial communities , 2008, BMC Bioinformatics.

[23]  Three proposed new bacteriophage genera of staphylococcal phages: “3alikevirus”, “77likevirus” and “Phietalikevirus” , 2014, Archives of Virology.

[24]  Forest Rohwer,et al.  Laboratory procedures to generate viral metagenomes , 2009, Nature Protocols.

[25]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[26]  Florent E. Angly,et al.  Viral diversity and dynamics in an infant gut. , 2008, Research in microbiology.

[27]  D. Caron,et al.  Marine bacterial, archaeal and protistan association networks reveal ecological linkages , 2011, The ISME Journal.

[28]  Steven Salzberg,et al.  Identifying bacterial genes and endosymbiont DNA with Glimmer , 2007, Bioinform..

[29]  G. Weinstock,et al.  Metagenomic analysis of double-stranded DNA viruses in healthy adults , 2014, BMC Biology.

[30]  C. Huttenhower,et al.  Metagenomic microbial community profiling using unique clade-specific marker genes , 2012, Nature Methods.

[31]  P. Salamon,et al.  Metagenomic Analyses of an Uncultured Viral Community from Human Feces , 2003, Journal of bacteriology.

[32]  François Enault,et al.  Assessment of viral community functional potential from viral metagenomes may be hampered by contamination with cellular sequences , 2013, Open Biology.

[33]  O. Dereure,et al.  Human Skin Microbiota: High Diversity of DNA Viruses Identified on the Human Skin by High Throughput Sequencing , 2012, PloS one.

[34]  Andrew C. Pawlowski,et al.  The Comprehensive Antibiotic Resistance Database , 2013, Antimicrobial Agents and Chemotherapy.

[35]  D. Relman,et al.  Evidence of a robust resident bacteriophage population revealed through analysis of the human salivary virome , 2011, The ISME Journal.

[36]  Shane S. Sturrock,et al.  Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data , 2012, Bioinform..

[37]  N. Moran,et al.  Bacteriophages Encode Factors Required for Protection in a Symbiotic Mutualism , 2009, Science.

[38]  Bonnie L Hurwitz,et al.  Depth-stratified functional and taxonomic niche specialization in the ‘core’ and ‘flexible’ Pacific Ocean Virome , 2014, The ISME Journal.

[39]  Jesse R. Zaneveld,et al.  Phage-bacteria network analysis and its implication for the understanding of coral disease. , 2015, Environmental microbiology.

[40]  Geoffrey D. Hannigan,et al.  Culture‐independent pilot study of microbiota colonizing open fractures and association with severity, mechanism, location, and complication from presentation to early outpatient follow‐up , 2014, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[41]  Wolf-Dietrich Hardt,et al.  Phages and the Evolution of Bacterial Pathogens: from Genomic Rearrangements to Lysogenic Conversion , 2004, Microbiology and Molecular Biology Reviews.

[42]  F. Raymond,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Ray Meta: scalable de novo metagenome assembly and profiling , 2012 .

[43]  Peter Salamon,et al.  PHACCS, an online tool for estimating the structure and diversity of uncultured viral communities using metagenomic information , 2005, BMC Bioinformatics.

[44]  Jian Yang,et al.  VFDB 2012 update: toward the genetic diversity and molecular evolution of bacterial virulence factors , 2011, Nucleic Acids Res..

[45]  Qi Zheng,et al.  GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis , 2008, Nucleic Acids Res..

[46]  S. Schuster,et al.  Integrative analysis of environmental sequences using MEGAN4. , 2011, Genome research.

[47]  R. Edwards,et al.  Fast Identification and Removal of Sequence Contamination from Genomic and Metagenomic Datasets , 2011, PloS one.

[48]  J. Banfield,et al.  Rapidly evolving CRISPRs implicated in acquired resistance of microorganisms to viruses. , 2007, Environmental microbiology.

[49]  María Martín,et al.  Activities at the Universal Protein Resource (UniProt) , 2013, Nucleic Acids Res..

[50]  Florent E. Angly,et al.  The Marine Viromes of Four Oceanic Regions , 2006, PLoS biology.

[51]  R. Knight,et al.  Bacterial Community Variation in Human Body Habitats Across Space and Time , 2009, Science.

[52]  Albert,et al.  Emergence of scaling in random networks , 1999, Science.

[53]  Richard E. Lenski,et al.  Effect of resource enrichment on a chemostat community of bacteria and bacteriophage , 1997 .

[54]  C. Alteri,et al.  Prophage Induction Is Enhanced and Required for Renal Disease and Lethality in an EHEC Mouse Model , 2013, PLoS pathogens.

[55]  Allyson L. Byrd,et al.  Biogeography and individuality shape function in the human skin metagenome , 2014, Nature.

[56]  Katherine H. Huang,et al.  Structure, Function and Diversity of the Healthy Human Microbiome , 2012, Nature.

[57]  Robert C. Edgar,et al.  PILER-CR: Fast and accurate identification of CRISPR repeats , 2007, BMC Bioinformatics.

[58]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[59]  Ye Deng,et al.  Functional Molecular Ecological Networks , 2010, mBio.

[60]  W. Jacobs,et al.  Origins of Highly Mosaic Mycobacteriophage Genomes , 2003, Cell.

[61]  Jun Yu,et al.  VFDB: a reference database for bacterial virulence factors , 2004, Nucleic Acids Res..

[62]  Paul Turner,et al.  Reagent and laboratory contamination can critically impact sequence-based microbiome analyses , 2014, BMC Biology.

[63]  Bernard Henrissat,et al.  Metabolic Reconstruction for Metagenomic Data and Its Application to the Human Microbiome , 2012, PLoS Comput. Biol..

[64]  Forest Rohwer,et al.  The GAAS Metagenomic Tool and Its Estimations of Viral and Microbial Average Genome Size in Four Major Biomes , 2009, PLoS Comput. Biol..

[65]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[66]  Julia Oh,et al.  Topographic diversity of fungal and bacterial communities in human skin , 2013, Nature.

[67]  Jaysheel D. Bhavsar,et al.  Metagenomics: Read Length Matters , 2008, Applied and Environmental Microbiology.

[68]  Geoffrey D. Hannigan,et al.  Microbial ecology of the skin in the era of metagenomics and molecular microbiology. , 2013, Cold Spring Harbor perspectives in medicine.

[69]  Charles A. Bowman,et al.  Propionibacterium acnes Bacteriophages Display Limited Genetic Diversity and Broad Killing Activity against Bacterial Skin Isolates , 2012, mBio.