Transcriptome analysis of bacteriophage communities in periodontal health and disease
© Santiago-Rodriguez et al. 2015
Received: 1 July 2014
Accepted: 15 July 2015
Published: 28 July 2015
The role of viruses as members of the human microbiome has gained broader attention with the discovery that human body surfaces are inhabited by sizeable viral communities. The majority of the viruses identified in these communities have been bacteriophages that predate upon cellular microbiota rather than the human host. Phages have the capacity to lyse their hosts or provide them with selective advantages through lysogenic conversion, which could help determine the structure of co-existing bacterial communities. Because conditions such as periodontitis are associated with altered bacterial biota, phage mediated perturbations of bacterial communities have been hypothesized to play a role in promoting periodontal disease. Oral phage communities also differ significantly between periodontal health and disease, but the gene expression of oral phage communities has not been previously examined.
Here, we provide the first report of gene expression profiles from the oral bacteriophage community using RNA sequencing, and find that oral phages are more highly expressed in subjects with relative periodontal health. While lysins were highly expressed, the high proportion of integrases expressed suggests that prophages may account for a considerable proportion of oral phage gene expression. Many of the transcriptome reads matched phages found in the oral cavities of the subjects studied, indicating that phages may account for a substantial proportion of oral gene expression. Reads homologous to siphoviruses that infect Firmicutes were amongst the most prevalent transcriptome reads identified in both periodontal health and disease. Some genes from the phage lytic module were significantly more highly expressed in subjects with periodontal disease, suggesting that periodontitis may favor the expression of some lytic phages.
As we explore the contributions of viruses to the human microbiome, the data presented here suggest varying expression of bacteriophage communities in oral health and disease.
Periodontal disease is widespread and affects approximately half of the United States adult population . It has been associated with systemic diseases such as atherosclerosis [2, 3], obesity , anemia of chronic disease , and diabetes [6–8], among other inflammatory conditions. The etiology of periodontal disease is likely multifactorial, with the oral microbiota playing a major role in its development. Many oral microbes have been linked to periodontitis, including Porphyromonas gingivalis , Tannerella forsythia , Aggregatibacter actinomycetemcomitans [11, 12], Streptococcus mutans , and Treponema denticola . However, there are cases of periodontal disease in which none of these potential pathogens have been isolated from affected oral sites. The realization that there is a significant community of indigenous microbes inhabiting the human body, including many uncultivable bacteria [15–18], archaea , fungi , and viruses [21–24] has led some to hypothesize that periodontal disease may develop as a result of the oral microbial communities rather than the presence of any single pathogen [25–28].
Not only are bacterial communities in the human oral biofilm altered in periodontal disease [29, 30], but so too are oral viral communities . Most studies of individual viruses associated with periodontal disease have been limited to cultivable eukaryotic viruses, such as Human Herpesvirus, Cytomegalovirus, and Epstein-Barr virus . While studies have shown associations between herpesviruses and severe periodontal disease , a causal link has yet to be definitively demonstrated [34, 35]. We now know that human oral viruses are highly diverse, and most that have been identified are bacteriophages that predate upon bacteria and archaea instead of the human host [31, 36]. Phage communities are highly personalized features of individual human subjects and have distinct membership based on the oral biogeographic site from which they were derived [36, 37]. These communities of phages likely play an active role in shaping the ecology of oral cellular microbiota through lysis of their hosts and lysogenic conversions . However, little is known about the transcriptional activity of these viral communities, or whether their patterns of expression reflect oral health status.
Despite the lack of studies characterizing oral phage gene expression, there have been efforts to characterize their bacterial counterparts [29, 38] and to associate patterns of transcription with disease severity . Patterns of cellular microbial gene expression are altered when subjects with relative periodontal health are compared with those with significant periodontal disease. Among those genes upregulated in Fusobacterium nucleatum, P. gingivalis, and A. actinomycetemcomitans are those coding for chaperones, transposases and ABC transporters . The abundance of some bacterial species is not associated with a similar increase in transcription from that species in the oral transcriptome, as was demonstrated for Acinetobacter baumannii; however, for other species such as P. gingivalis, T. forsythia, and T. denticola, relative abundance is associated with concomitant increases in transcription. These data suggest that there are factors other than oral health status that affect the transcription of the cellular microbiota, but it is unclear whether such a phenomenon exists for oral phage communities.
Although viruses make up a significant portion of the genes in any community, there are few studies focusing on the contribution viruses make to community gene expression. We sought to provide the first characterization of the contributions of viral communities to the oral transcriptome by comparing meta-transcriptomes with the virome sequences isolated from the same group of subjects. Our goals were to: 1) characterize viral transcription in the oral transcriptome, 2) identify specific viruses whose genes are transcribed in a cohort of human subjects, and 3) discern whether there are differences in viral gene expression patterns that correspond to oral health status.
Human subjects and RNA enrichment
Taxonomic distributions found in oral phage RNA
Percentage of reads assigned to phage families
Expression of individual viruses in the oral cavity
Gene expression in lytic and lysogenic cycles
The study of human viral communities is in its relative infancy compared to human cellular microbiota. The majority of studies characterize bacterial communities by 16S rRNA amplicon sequencing, while far fewer utilize metagenome or transcriptome sequencing. In contrast to microarrays, RNA sequencing does not require pre-existing genome sequence data from each microbe present, but rather can characterize gene expression based on sequence homology to existing databases. The data presented here provide a first glimpse into the in vivo behavior of the diverse communities of phages that inhabit the human oral cavity. However, despite the benefits, its main limitation is in a relative inability to discern phage sequences from the abundant bacterial and archaeal sequences present. Our data illustrate this point, where only a minority of the sequence reads have homology to known viruses, including phages in the SEED database compared to bacteria (Additional file 1: Table S2), while a far greater proportion of the reads match viruses assembled from the viromes of each subject (Additional file 2: Figure S3). We identified only 0.16 % of the reads homologous to known viruses in this study, which is greater than a recent oral microbiome study where only 0.003 % of their reads were homologous to known viruses . That study, however, did not compare their metagenome reads to viromes from the same subjects, which likely would have demonstrated that a far greater proportion of their metagenomes were attributable to viruses. Viral communities are known in a number of different environments to harbor gene functions also present in their bacterial hosts [21, 42], which could result in an overestimation of the gene functions we ascribed to phages in this study. Conversely, when these genes are not properly annotated as phages in their host bacterial genomes, they could result in an underestimation of the proportion of the transcriptome sequences attributable to phages.
There are significant differences between viral communities in periodontal health and disease , but prior studies have had limited insight into the gene expression of viruses in the oral cavity. The data produced in this study suggest that some lytic modules are active in subjects with significant periodontal disease (Fig. 5). Lysins were among the most highly expressed phage genes, and they were differentially expressed in periodontal health and disease, although the results were not statistically significant. We hypothesize from these data that the inflammatory environment in periodontal disease may support the lytic cycle of some temperate phages, which may be related to bacterial stress. Shifts in pH, the presence of heavy metals and toxic compounds are known to cause bacterial stress and consequently the production of virus-like particles in certain bacteria . Periodontal disease is also known to promote bacterial stress, resulting in bacterial resistance to acidic or alkaline conditions; yet, how this potentially triggers the lytic cycle of temperate phages remains to be addressed . While a link between oral health status and viral lytic/lysogenic functions needs to be demonstrated experimentally, the data from this study provides an intriguing basis to support further study. A similar study with the same group of subjects showed a high proportion of putative myoviruses in the subgingival crevices of subjects with periodontal disease , suggesting that the lytic cycle may be preferred in subjects with significant periodontal disease. Perhaps the significant proportion of lysins identified in our cohort suggests a “kill the winner” dynamic in subjects with significant periodontal disease, but a longitudinal study with a larger group of subjects would be necessary to support that hypothesis.
While this study provides evidence for differential gene phage expression in subjects with periodontal health and disease, the are some confounders in the subject population that could affect the results. The age of the subjects studied was significantly different, where subjects with periodontal health had a mean age of 35 years, and subjects with periodontal disease had a mean age of 60. Very little is known about the effects of age on the human oral virome; however, in a prior study of subjects with periodontal health, we observed no impacts of age on the oral virome . In fact, we found that there were significant numbers of shared phage between household members regardless of their ages. Similarly, we found no specific effect of ethnicity on virome contents, but none of our virome studies [21, 31, 36, 37, 45, 46] have been powered adequately to characterize any potential effects of ethnicity on the human oral virome. In this study, most of the subjects with periodontal disease were Caucasian compared to a number of different ethnicities amongst those subjects with periodontal health; these subject characteristics largely reflect the patient population at our study site. A much larger study focused on the effects of age and ethnicity is needed to characterize any effects of potential confounders on the study of the human oral virome.
This study offers some benefits not found in other RNA-sequencing studies, including that we sequenced the RNA from a diverse group of different subjects and demonstrated trends across subject groups, and we enriched the RNA to remove many ribosomal and human RNAs to characterize a more robust dataset. Because the phages characterized across each subject or subject group were highly diverse, traditional techniques such as RT-PCR to verify trends seen in sequence data were not possible. For example, all subjects have integrases, lysins, repressors, and tail fibers, but most of these gene sequences are unique to each subject and phage; therefore, gene expression could only be reproduced from a single phage using RT-PCR, which would not capture trends seen throughout that community. Thus, while the expression of any individual phage may vary, the overall expression in the community may show significant trends.
There were some gene expression differences across the phage communities in periodontal health and disease (Fig. 5), despite there being roughly similar proportions of phage families in each state (Table 1). Genes associated with penetration of the cell membrane (holins), promotion of lytic gene expression (antirepressors), and lysis (lysins) were expressed at higher levels in periodontal disease than in periodontal health (Fig. 5). This may support a preference of certain phages for the lytic cycle in periodontal disease, although experimental demonstration would be necessary to prove such a phenomenon. Most genes associated with lysogeny, including integrases and repressors, were not significantly different based on oral health status. Notably, there was no concerted expression of phage holins and lysins in periodontal health or disease. Some phages lyse their hosts using a holin-endolysin system; however, in others, there is no interaction between the holin and the endolysin genes in host lysis . Some of these endolysins are secreted to the extracellular environment before completion of the viral lysis cycle independent of holin activity. Additionally, some holin genes are not located downstream of lysins , suggesting that holin and lysin gene expression in these phages is independent.
We believed that by analyzing the viromes in each of these subjects, we would find that the most abundant viruses in each subject also would have the greatest gene expression. Indeed, a relative minority of the phages present in each subject had detectable transcripts. Many of the transcriptome reads from different subjects with periodontal health and disease matched the same viruses, suggesting that there is some conservation of phage genes across all subjects. This phenomenon has previously been demonstrated in separate cohorts in the saliva of healthy human subjects and in the dental plaque of subjects with periodontal health and disease [31, 37]. We believe that part of the explanation why only a minority of the phages in each cohort had identifiable transcripts is related to the relative longevity of viruses in the oral cavity. We recently demonstrated that human oral viruses are remarkably persistent over short and long time periods [37, 49], but transcription from these viruses was not necessarily detected. The longevity of these viruses may not be directly related to their ongoing transcription and production of virions, but rather their ability to attach and persist on mucosal surfaces . Thus, the persistence of oral viruses may be secondary to their ability to attach to surfaces, which would reduce their need for ongoing replication and transcription to produce virions.
We know very little about the role of viruses on human body surfaces, but the oral virome now is coming into greater context. The human oral cavity is populated by numerous bacteriophages that carry substantial gene functions putatively involved in the pathogenesis of their cellular hosts [21, 22]. Many phages are highly persistent members of the oral microbiome, with transience the exception rather than the rule . Oral phages are also highly specific to individual, biogeographic site, host-sex, and oral health status [31, 37]. The RNA-sequencing data here indicates that there is substantial transcriptional activity from oral phage communities, and that oral health status may play some role in the activity of viruses. While there does not appear to be great differences in the expression of the lysogenic modules of oral phages, lytic module expression is somewhat increased in subjects with significant periodontal disease. The development of a model oral environment to investigate the interactions between phages and their cellular hosts would be ideal to substantiate the role and activity of the oral phage community in periodontal health and disease.
Subject recruitment and enrollment was approved by the University of California, San Diego and the Western University Administrative Panels on Human Subjects in Medical Research. All subjects signed informed consent indicating their willingness to participate in this study.
Study criteria and sample collection
Each subject underwent a baseline periodontal examination including measurements of probing depths, clinical attachment loss, Gingival Index, Plaque Index, and gingival irritation , and their oral health status was recorded (Additional file 1: Table S1). We used the 1999 International Workshop for Classification of Periodontal Diseases and Conditions, where periodontitis including juvenile forms of periodontitis is defined by loss of attachment. For a diagnosis of healthy, all sites had to have an attachment level of 0 mm and an absence of bleeding on probing. We excluded attachment levels from sites that were located next to 3rd molars, edentulous areas and sites where attachment loss was clearly caused by factors other than periodontal disease such as chronic toothbrush trauma. For those subjects in this study diagnosed as having moderate gingivitis and mild periodontitis, they met criteria for periodontal health according to this classification scheme. Approximately 3 ml of saliva was collected from each subject without stimulation, immediately frozen on dry ice and stored at −20 °C until use in this study. All subjects completed a questionnaire detailing their dietary habits and other co-morbidities. Exclusion criteria included preexisting medical conditions that could result in substantial immunosuppression.
RNA extraction, enrichment and sequencing
Total RNA was extracted from each subject using the Mirvana kit (Life Technologies, Grand Island, NY), with the inclusion of a bead-beating step for 20 min with Lysing-Matrix B (MP Bio, Santa Ana, CA). Total RNA then was enriched for microbial RNA using MicrobEnrich (Life Technologies), and further enriched for mRNA using MicrobExpress (Life Technologies) and MegaClear (Life Technologies), which are designed to remove ribosomal RNAs. Enriched RNA then was prepared for sequencing through the construction of cDNA libraries using the Ion Total RNA-Seq kit (Life Technologies), and subjected to successive rounds of Ampure bead purification (Beckman-Coulter, Brea, CA) to remove small cDNAs. Libraries were quantified using an Agilent Bioanalyzer HS DNA Kit (Agilent, Santa Clara, CA) and then were sequenced on 316 chips using an Ion Torrent Personal Genome Machine , producing an average of 1,500,304 reads per subject of mean length 96 nucleotides. All sequences are available for download in the MG-RAST database (metagenomics.anl.gov/) under the project #5691 or under the name ‘Health_And_Disease_Transcriptome_Project.’
Processing of RNA sequences
Sequencing reads were trimmed according to modified Phred quality scores of 0.5 using CLC Genomics Workbench 6.51 (CLC bio USA, Cambridge, MA). The remaining reads were further processed for quality control by removing reads with substantial length variation (reads <50 nucleotides or >200 nucleotides), or reads where ≥25 % of the length was due to homopolymers tracts. Each transcriptome was screened for contaminating human nucleic acids using BLASTN analysis (E-score <10−5) against the human reference database available at ftp://ftp.ncbi.nlm.nih.gov/genomes/H_sapiens/. Any reads homologous to human sequences were removed prior to further analysis.
Analysis of transcriptomes
Analysis of gene expression was determined by BLASTX analysis of the SEED database using MG-RAST (E-score <10−5) . Statistical significance was determined by comparing the means for all subjects for all subsystems between subject groups by two-tailed t-tests using Microsoft Excel 2007 (Microsoft Corp., Redman, WA). Analysis of the expression differences between subject groups for individual genes also was determined by comparisons of means. The data for each individual gene were compared between MG-RAST and the normalized RPKM values obtained from CLC Genomics Workbench 6.51 to verify that they produced similar results. The relative abundances of virus families were determined by BLASTX analysis of the SEED database using MG-RAST . Phages were classified into Caudovirus families Myoviridae, Siphoviridae, and Podoviridae based on homology to known viruses from these families as previously described . Analysis of transcription of individual viruses was determined by mapping transcriptome reads against a database of virome contigs from each subject . Read mappings were performed at 98 % identity over the entire length of each sequenced read, which were defined as matches. Heatmaps were generated based on the normalized proportions of reads mapping to pooled viromes from each subject, and were visualized using JAVA Treeview . Principal coordinates analysis was performed using binary Sorensen distances using Qiime 1.1 . Statistical comparisons of means were performed by two-tail t-tests and equality of variances for data that was not normally distributed was performed by F-tests. Each was determined using Microsoft Excel 2007 (Microsoft Corp., Redman, WA).
Supported by the Robert Wood Johnson Foundation, the Burroughs Wellcome Fund, and NIH 1K08AI085028 to DTP.
- Eke PI, Thornton-Evans G, Dye B, Genco R. Advances in surveillance of periodontitis: the Centers for Disease Control and Prevention periodontal disease surveillance project. J Periodontol. 2012;83:1337–42.PubMedView ArticleGoogle Scholar
- Gocke C, Holtfreter B, Meisel P, Grotevendt A, Jablonowski L, Nauck M, et al. Abdominal obesity modifies long-term associations between periodontitis and markers of systemic inflammation. Atherosclerosis. 2014;235:351–7.PubMedView ArticleGoogle Scholar
- Bokhari SA, Khan AA, Butt AK, Hanif M, Izhar M, Tatakis DN, et al. Periodontitis in coronary heart disease patients: strong association between bleeding on probing and systemic biomarkers. J Clin Periodontol. 2014;41(11):1048–54.PubMedView ArticleGoogle Scholar
- Suvan J, D'Aiuto F, Moles DR, Petrie A, Donos N. Association between overweight/obesity and periodontitis in adults. A systematic review. Obes Rev. 2011;12:e381–404.PubMedView ArticleGoogle Scholar
- Patel MD, Shakir QJ, Shetty A. Interrelationship between chronic periodontitis and anemia: a 6-month follow-up study. J Indian Soc Periodontol. 2014;18:19–25.PubMed CentralPubMedView ArticleGoogle Scholar
- Perayil J, Suresh N, Fenol A, Vyloppillil R, Bhaskar A, Menon S. Comparison of Hba1c levels in non-diabetic healthy subjects and subjects with periodontitis before and after non-surgical periodontal therapy. J Periodontol. 2014;1–17.Google Scholar
- Raman RP, Taiyeb-Ali TB, Chan SP, Chinna K, Vaithilingam RD. Effect of nonsurgical periodontal therapy verses oral hygiene instructions on Type 2 diabetes subjects with chronic periodontitis: a randomised clinical trial. BMC Oral Health. 2014;14:79.PubMed CentralPubMedView ArticleGoogle Scholar
- Corbella S, Francetti L, Taschieri S, De Siena F, Fabbro MD. Effect of periodontal treatment on glycemic control of patients with diabetes: a systematic review and meta-analysis. J Diabetes Investig. 2013;4:502–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Darveau RP, Hajishengallis G, Curtis MA. Porphyromonas gingivalis as a potential community activist for disease. J Dent Res. 2012;91:816–20.PubMed CentralPubMedView ArticleGoogle Scholar
- Gersdorf H, Meissner A, Pelz K, Krekeler G, Gobel UB. Identification of Bacteroides forsythus in subgingival plaque from patients with advanced periodontitis. J Clin Microbiol. 1993;31:941–6.PubMed CentralPubMedGoogle Scholar
- Ximenez-Fyvie LA, Haffajee AD, Socransky SS. Microbial composition of supra- and subgingival plaque in subjects with adult periodontitis. J Clin Periodontol. 2000;27:722–32.PubMedView ArticleGoogle Scholar
- Fine DH, Markowitz K, Fairlie K, Tischio-Bereski D, Ferrendiz J, Furgang D, et al. A consortium of Aggregatibacter actinomycetemcomitans, Streptococcus parasanguinis, and Filifactor alocis is present in sites prior to bone loss in a longitudinal study of localized aggressive periodontitis. J Clin Microbiol. 2013;51:2850–61.PubMed CentralPubMedView ArticleGoogle Scholar
- Loesche WJ. Role of Streptococcus mutans in human dental decay. Microbiol Rev. 1986;50:353–80.PubMed CentralPubMedGoogle Scholar
- Sela MN. Role of Treponema denticola in periodontal diseases. Crit Rev Oral Biol Med. 2001;12:399–413.PubMedView ArticleGoogle Scholar
- Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI. The human microbiome project. Nature. 2007;449:804–10.PubMed CentralPubMedView ArticleGoogle Scholar
- 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:1355–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Huse SM, Ye Y, Zhou Y, Fodor AA. A core human microbiome as viewed through 16S rRNA sequence clusters. PLoS One. 2012;7, e34242.PubMed CentralPubMedView ArticleGoogle Scholar
- Aas JA, Paster BJ, Stokes LN, Olsen I, Dewhirst FE. Defining the normal bacterial flora of the oral cavity. J Clin Microbiol. 2005;43:5721–32.PubMed CentralPubMedView ArticleGoogle Scholar
- Eckburg PB, Lepp PW, Relman DA. Archaea and their potential role in human disease. Infect Immun. 2003;71:591–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Ding T, Schloss PD. Dynamics and associations of microbial community types across the human body. Nature. 2014;509:357–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Pride DT, Salzman J, Haynes M, Rohwer F, Davis-Long C, White 3rd RA, et al. Evidence of a robust resident bacteriophage population revealed through analysis of the human salivary virome. ISME J. 2012;6:915–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Willner D, Furlan M, Schmieder R, Grasis JA, Pride DT, Relman DA, et al. Metagenomic detection of phage-encoded platelet-binding factors in the human oral cavity. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4547–53.PubMed CentralPubMedView ArticleGoogle Scholar
- Breitbart M, Hewson I, Felts B, Mahaffy JM, Nulton J, Salamon P, et al. Metagenomic analyses of an uncultured viral community from human feces. J Bacteriol. 2003;185:6220–3.PubMed CentralPubMedView ArticleGoogle Scholar
- Breitbart M, Haynes M, Kelley S, Angly F, Edwards RA, Felts B, et al. Viral diversity and dynamics in an infant gut. Res Microbiol. 2008;159:367–73.PubMedView ArticleGoogle Scholar
- Belda-Ferre P, Alcaraz LD, Cabrera-Rubio R, Romero H, Simon-Soro A, Pignatelli M, et al. The oral metagenome in health and disease. ISME J. 2012;6:46–56.PubMed CentralPubMedView ArticleGoogle Scholar
- Lepp PW, Brinig MM, Ouverney CC, Palm K, Armitage GC, Relman DA. Methanogenic Archaea and human periodontal disease. Proc Natl Acad Sci U S A. 2004;101:6176–81.PubMed CentralPubMedView ArticleGoogle Scholar
- Griffen AL, Beall CJ, Campbell JH, Firestone ND, Kumar PS, Yang ZK, et al. Distinct and complex bacterial profiles in human periodontitis and health revealed by 16S pyrosequencing. ISME J. 2012;6:1176–85.PubMed CentralPubMedView ArticleGoogle Scholar
- Abusleme L, Dupuy AK, Dutzan N, Silva N, Burleson JA, Strausbaugh LD, et al. The subgingival microbiome in health and periodontitis and its relationship with community biomass and inflammation. ISME J. 2013;7:1016–25.PubMed CentralPubMedView ArticleGoogle Scholar
- Li Y, He J, He Z, Zhou Y, Yuan M, Xu X, et al. Phylogenetic and functional gene structure shifts of the oral microbiomes in periodontitis patients. ISME J. 2014;8(9):1879–91.PubMedView ArticleGoogle Scholar
- Belstrom D, Fiehn NE, Nielsen CH, Kirkby N, Twetman S, Klepac-Ceraj V, et al. Differences in bacterial saliva profile between periodontitis patients and a control cohort. J Clin Periodontol. 2014;41:104–12.PubMedView ArticleGoogle Scholar
- Ly M, Abeles SR, Boehm TK, Robles-Sikisaka R, Naidu M, Santiago-Rodriguez T, et al. Altered oral viral ecology in association with periodontal disease. MBio. 2014;5:e01133–01114.PubMed CentralPubMedView ArticleGoogle Scholar
- Cappuyns I, Gugerli P, Mombelli A. Viruses in periodontal disease - a review. Oral Dis. 2005;11:219–29.PubMedView ArticleGoogle Scholar
- Imbronito AV, Okuda OS, Maria de Freitas N, Moreira Lotufo RF, Nunes FD. Detection of herpesviruses and periodontal pathogens in subgingival plaque of patients with chronic periodontitis, generalized aggressive periodontitis, or gingivitis. J Periodontol. 2008;79:2313–21.PubMedView ArticleGoogle Scholar
- Chalabi M, Rezaie F, Moghim S, Mogharehabed A, Rezaei M, Mehraban B. Periodontopathic bacteria and herpesviruses in chronic periodontitis. Mol Oral Microbiol. 2010;25:236–40.PubMedView ArticleGoogle Scholar
- Stein JM, Said Yekta S, Kleines M, Ok D, Kasaj A, Reichert S, et al. Failure to detect an association between aggressive periodontitis and the prevalence of herpesviruses. J Clin Periodontol. 2013;40:1–7.PubMedView ArticleGoogle Scholar
- Naidu M, Robles-Sikisaka R, Abeles SR, Boehm TK, Pride DT. Characterization of bacteriophage communities and CRISPR profiles from dental plaque. BMC Microbiol. 2014;14:175.PubMed CentralPubMedView ArticleGoogle Scholar
- Abeles SR, Robles-Sikisaka R, Ly M, Lum AG, Salzman J, Boehm TK, et al. Human oral viruses are personal, persistent and gender-consistent. ISME J. 2014;8(9):1753–67.PubMedView ArticleGoogle Scholar
- Jorth P, Turner KH, Gumus P, Nizam N, Buduneli N, Whiteley M. Metatranscriptomics of the human oral microbiome during health and disease. MBio. 2014;5:e01012–01014.PubMed CentralPubMedView ArticleGoogle Scholar
- Kebschull M, Demmer RT, Grun B, Guarnieri P, Pavlidis P, Papapanou PN. Gingival tissue transcriptomes identify distinct periodontitis phenotypes. J Dent Res. 2014;93:459–68.PubMed CentralPubMedView ArticleGoogle Scholar
- Lucchini S, Desiere F, Brussow H. Similarly organized lysogeny modules in temperate Siphoviridae from low GC content gram-positive bacteria. Virology. 1999;263:427–35.PubMedView ArticleGoogle Scholar
- Sandaa RA. Environmental Viral Pool. In: Schmidt TM, Schaechter M, editors. Topics in Environmental and Ecological Microbiology. Elsevier; 2012. p.101–3.Google Scholar
- Angly FE, Felts B, Breitbart M, Salamon P, Edwards RA, Carlson C, et al. The marine viromes of four oceanic regions. PLoS Biol. 2006;4, e368.PubMed CentralPubMedView ArticleGoogle Scholar
- Choi J, Kotay SM, Goel R. Various physico-chemical stress factors cause prophage induction in Nitrosospira multiformis 25196-an ammonia oxidizing bacteria. Water Res. 2010;44:4550–8.PubMedView ArticleGoogle Scholar
- Chavez de Paz LE, Bergenholtz G, Dahlen G, Svensater G. Response to alkaline stress by root canal bacteria in biofilms. Int Endod J. 2007;40:344–55.PubMedView ArticleGoogle Scholar
- Robles-Sikisaka R, Ly M, Boehm T, Naidu M, Salzman J, Pride DT. Association between living environment and human oral viral ecology. ISME J. 2013;7:1710–24.PubMed CentralPubMedView ArticleGoogle Scholar
- Robles-Sikisaka R, Naidu M, Ly M, Salzman J, Abeles SR, Boehm TK, et al. Conservation of streptococcal CRISPRs on human skin and saliva. BMC Microbiol. 2014;14:146.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang IN, Smith DL, Young R. Holins: the protein clocks of bacteriophage infections. Annu Rev Microbiol. 2000;54:799–825.PubMedView ArticleGoogle Scholar
- Catalão MJ, Gil F, Moniz-Pereira J, Sao-Jose C, Pimentel M. Diversity in bacterial lysis systems: bacteriophages show the way. FEMS Microbiol Rev. 2013;37:554–71.PubMedView ArticleGoogle Scholar
- Barr JJ, Auro R, Furlan M, Whiteson KL, Erb ML, Pogliano J, et al. Bacteriophage adhering to mucus provide a non-host-derived immunity. Proc Natl Acad Sci U S A. 2013;110:10771–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Loe H. The gingival index, the plaque index and the retention index systems. J Periodontol. 1967;38:Suppl:610–6.View ArticleGoogle Scholar
- Rothberg JM, Hinz W, Rearick TM, Schultz J, Mileski W, Davey M, et al. An integrated semiconductor device enabling non-optical genome sequencing. Nature. 2011;475:348–52.PubMedView ArticleGoogle Scholar
- Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386.PubMed CentralPubMedView ArticleGoogle Scholar
- Rohwer F, Edwards R. The Phage Proteomic Tree: a genome-based taxonomy for phage. J Bacteriol. 2002;184:4529–35.PubMed CentralPubMedView ArticleGoogle Scholar
- Saldanha AJ. Java Treeview--extensible visualization of microarray data. Bioinformatics. 2004;20:3246–8.PubMedView ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.