Preliminary characterization of the oral microbiota of Chinese adults with and without gingivitis
- Shi Huang†1,
- Fang Yang†4, 5,
- Xiaowei Zeng1,
- Jie Chen1,
- Rui Li2Email author,
- Ting Wen2,
- Chun Li2,
- Wei Wei2,
- Jiquan Liu2,
- Lan Chen2,
- Catherine Davis3 and
- Jian Xu1Email author
© Huang et al; licensee BioMed Central Ltd. 2011
Received: 6 May 2011
Accepted: 12 December 2011
Published: 12 December 2011
Microbial communities inhabiting human mouth are associated with oral health and disease. Previous studies have indicated the general prevalence of adult gingivitis in China to be high. The aim of this study was to characterize in depth the oral microbiota of Chinese adults with or without gingivitis, by defining the microbial phylogenetic diversity and community-structure using highly paralleled pyrosequencing.
Six non-smoking Chinese, three with and three without gingivitis (age range 21-39 years, 4 females and 2 males) were enrolled in the present cross-sectional study. Gingival parameters of inflammation and bleeding on probing were characterized by a clinician using the Mazza Gingival Index (MGI). Plaque (sampled separately from four different oral sites) and salivary samples were obtained from each subject. Sequences and relative abundance of the bacterial 16 S rDNA PCR-amplicons were determined via pyrosequencing that produced 400 bp-long reads. The sequence data were analyzed via a computational pipeline customized for human oral microbiome analyses. Furthermore, the relative abundances of selected microbial groups were validated using quantitative PCR.
The oral microbiomes from gingivitis and healthy subjects could be distinguished based on the distinct community structures of plaque microbiomes, but not the salivary microbiomes. Contributions of community members to community structure divergence were statistically accessed at the phylum, genus and species-like levels. Eight predominant taxa were found associated with gingivitis: TM7, Leptotrichia, Selenomonas, Streptococcus, Veillonella, Prevotella, Lautropia, and Haemophilus. Furthermore, 98 species-level OTUs were identified to be gingivitis-associated, which provided microbial features of gingivitis at a species resolution. Finally, for the two selected genera Streptococcus and Fusobacterium, Real-Time PCR based quantification of relative bacterial abundance validated the pyrosequencing-based results.
This methods study suggests that oral samples from this patient population of gingivitis can be characterized via plaque microbiome by pyrosequencing the 16 S rDNA genes. Further studies that characterize serial samples from subjects (longitudinal study design) with a larger population size may provide insight into the temporal and ecological features of oral microbial communities in clinically-defined states of gingivitis.
Keywordsoral microbiota gingivitis saliva plaque pyrosequencing
Metagenomic techniques have recently revolutionized our understanding of the plethora of microbes that co-inhabit the human body, collectively known as the human microbiome. Various body sites (e.g. the skin, the gastrointestinal and vaginal tracts and the oral cavity) harbor distinct communities of microbes that vary among host individuals as well as among the ecological niches within each body site . Interactions among resident microbiota and between the microbiota and the human host underlie human health and disease. Within the oral cavity, the tongue, soft and hard palates, buccal mucosa, supragingival and subgingival surfaces of the teeth and saliva may represent different ecological niches or habitats . The composition and diversity of microbiota in these habitats may contribute to oral health [3–6] and oral diseases such as dental caries, periodontitis, and gingivitis [7, 8].
Gingivitis is inflammation of the soft tissues of the gum surrounding the teeth. It is believed to result from the build-up of plaque  and the ensuing interactions between the plaque microbiota and host tissues [10, 11]. These tissues become erythematous and bleed upon probing, but no apical migration of the junctional epithelium occurs. Previous studies of gingival plaque showed that as gingivitis develops, the microbial constituents of subgingival plaque shift from a population dominated by Gram-positive streptococci to one with elevated levels of Gram-negative anaerobes such as Actinobacillus, Capnocytophaga, Campylobacter, Eikenella, Fusobactrium and Prevotella [2, 12, 13]. However, these studies have been based on culture-based and molecular methods that target only a limited and partial number of culturable microbes, a bias that can be overcome by metagenomic approaches. During the last decade, high-throughput sequencing approaches based on 16 S rDNA amplicons have been used to survey the diversity of human oral microbiota in health and disease. Notably, these techniques revealed the microbial diversity within the healthy subgingival crevice exceeds far beyond that was characterized previously. Kroes et al. noted that less than a quarter (24%) of phylotypes identified with metagenomic techniques could be recovered by cultivation and that almost half of the subgingival phylotypes identified with a combination of molecular and culture-based techniques had not been characterized previously . In fact, another study estimated that in human oral cavity approximately 68% of all bacterial taxa were still uncultivated .
Although molecular techniques have been used to compare subgingival plaques in healthy hosts and those with oral diseases such as periodontitis , few studies have investigated in depth the oral microbiota associated with gingivitis. There are several reasons. First, the depth and breadth of sampling for oral microbiota have been insufficient in general, and the optimal parameters not determined for that of gingivitis patients in particular. Second, regarding the selection of gingival sites for plaque sampling, it was still not clear whether or which of the different sites are clinically relevant (e.g., anterior teeth or posterior teeth? supragingival plaque or subgingival plaque?). Such ambiguity severely limits meaningful data analysis and comparisons across studies, and delays the translation of the findings clinically. Third, most oral microbial surveys that enumerated 16 S sequences of PCR-amplicons have ignored potential PCR artifacts [16–19]; as a result, a comprehensive and accurate organismal landscape of most oral microbiomes, particularly those related to diseases, remained largely elusive. All these factors have confounded the assessment of microbial factors associated with gingivitis.
Employing pyrosequencing of 16 S rDNA amplicons, this article elucidated the diversity and population structure of the oral microbiota, sampled respectively from five oral ecological niches from each of the three Chinese adults with gingivitis and three without the disease. Microbiota of supragingival plaque, subgingival plaque, and saliva were characterized to test whether and how the microbiomes from the various oral ecological niches distinguished healthy hosts and those with gingivitis. Our study pinpointed a number of organisms as potential biomarkers of gingivitis, and provided important insights for the sampling and analysis strategies for unraveling gingivitis-associated microbial risk factors in human populations.
Metadata for the six subjects sampled in this study.
Estimates of species diversity for the samples.
OTUs at 3% difference
Richness and diversity analysis based on Operational Taxonomic Units
Comparisons of community structures
Taxonomy-based characterization of oral microbiota
A total of 70 genera were identified in the oral microbiota from the five sampled sites (Additional file 1). The most frequently detected taxa at the genus level (the 12 most abundant genera that each represents at least 2% of oral microbiome) were Streptococcus, Neisseria, Leptotrichia, Actinomyces, Prevotella , Veillonella, Rothia, Fusobacterium, Lautropia, Selenomonas, Haemophilus, Granulicatella.
Bacterial genera differentially distributed between H and U groups based on the oral "CORE" database.
mean abundance (%)
mean abundance (%)
OTU-level comparisons of oral microbiota
As the above comparisons on the relative abundance of microbial taxa were derived from only those reads with confidence value above 0.8 (Methods), essentially all those reads without any reliable phylogenetic assignments (0.33~27.86% of the total reads in each sample) had been masked. Therefore, we further compared, between the two host-groups, the relative abundance of all species-level OTUs in plaque microbiota, independently of their phylogenetic identity assignments. In total, 98 OTUs (accounted for 5.38% of all OTUs) were found differentially distributed at the significance cutoff of 0.05 (both p-value and q-value of Metastats) in not only each of the four plaque sites but in all the plaque sites. Interestingly, among those 98 OTUs, 36 of them were gingivitis-depleted and the remaining 62 were gingivitis-enriched (Additional file 2). Consensus taxonomy of the OTUs was interrogated by MOTHUR based on oral "CORE" 16 S rDNA Gene Database (Methods). Twenty-four out of the 36 gingivitis-deplete OTUs while 57 out of the 62 gingivitis-enriched OTUs were supported by the taxonomy-assignment based results. Thus results from the two methodologies are largely consistent. These gingivitis-depleted or gingivitis-enriched OTUs represent a novel set of potential organismal markers for evaluation and prognosis of gingivitis, although their validity and significance remain to be further tested in larger human populations.
Correlation of microbial quantifications between pyrosequencing and quantitative PCR (qPCR)
The abundances, in gene copies per ng of DNA, of two genera (Streptococcus and Fusobacterium) in all of the plaque and saliva samples were determined with qPCR. The relative abundance of these two genera as determined by pyrosequencing showed a positive correlation to the qPCR-based gene copies per ng of DNA (Streptococcus: r = 0.554; p < 0.002; Fusobacterium: r = 0.813; p < 0.001) (Additional file 3).
This study employed highly paralleled pyrosequencing of 16 S rDNA to assess and compare the diversity and population structure of microbiota associated with gingivitis in Chinese adults. The microbial diversity in plaque and saliva estimated in our study, 464~737 OTUs (97% identity cutoff) in each sample, was similar to that reported by Zaura et al (saliva; ). The Zaura study employed a stringent and conservative read-trimming strategy, where only those reads present at least five times in one sample were taken into analysis. In our analysis, stringent quality-based read-trimming suggested by MOTHUR was performed, requiring average quality score of over 35 in a 50 bp moving window along the whole read (Methods). This conservative selection criterion http://www.mothur.org/wiki/ significantly reduced the number of OTUs from the estimates based on alternative read-trimming criteria such as requiring average base quality score > 25 (data not shown). Thus potential sequencing artifacts might inflate the observed bacterial diversity. Furthermore, the estimated Good's coverage showed that most of the bacterial phylotypes (> 97%) in the saliva and plaque of these healthy and gingivitis hosts were already identified in this study. The richness estimator of ACE and Chao1 also suggested that the majority of phylotypes (> 97%) were already represented by the sequences in our study.
Our study firstly aimed to assess whether communities from healthy and gingivitis-associated host populations differ in any specific site(s) of oral cavity. Both FastUnifrac-based and thetaYC-based analysis showed that saliva and plaque samples represented distinct microbiomes in the oral cavity. Regardless of disease status, salivary microbiota clustered distinctly from plaque microbiota, in each of the two distance matrixes tested. This likely reflected the different environmental conditions characterizing the two habitats. Plaque microbiota reside in biofilms on the tooth enamel surface and are affected by dietary composition, oral hygiene practices , microbial interactions within the biofilm  and interactions between microbes and host epithelial cells [10, 11, 23]. In contrast, the salivary habitat was shaped by food intake flux, transient microbiota, mucins, serous exudate, sloughed epithelial cells, etc [3, 4, 24, 25]. Interestingly, a survey of global diversity of the human salivary microbiome in ten individuals from each of twelve geographic locations worldwide (including China) reported a high diversity within and between host individuals but little geographic structure in the saliva microbiomes .
Secondly, members of the bacterial communities were identified. Furthermore, their contributions to the structural segregation of plaque microbiota between the two host populations were evaluated. When plaque microbiota were considered at the level of phylum, Fusobacteria and TM7, two of the predominant phyla, were more abundant in microbiota associated with gingivitis, while Actinobacteria and Bacteroidetes were less abundant in gingivitis-associated microbiomes. At the level of genus, several genera such as Leptotrichia and Selenomonas were more abundant in gingivitis plaque (21 such genera in total; Table 3), whereas only five genera, Streptococcus, Veillonella, Prevotella, Lautropia and Haemophilus, were less abundant. At species level, phylogeny-assignment independent comparison of relative abundances of OTUs between the healthy and gingivitis hosts was performed for not only each of the four plaque sites but also all of the plaque sites. Consistent with the above findings, 98 gingivitis-associated (both gingivitis-enriched and gingivitis-depleted) OTUs were pinpointed and found distributed in all sampled sites of plaque. Moreover, 58 OTUs affiliated to the genera of Leptotrichia (16), Selenomonas (12), Streptococcus (7), Veillonella (6), Prevotella (6), Lautropia (2), Haemophilus (3) and the candidate division TM7 (6) were found to be associated with gingivitis.
Notably, several members of these gingivitis-associated taxa were known to play a role in both oral health and disease. The gingivitis-enriched genus Leptotrichia, of the Fusobacteria phylum and Fusobacteriaceae family, were Gram-negative non-sporing-forming, anaerobic, saccharolytic rods. They were among the normal microbiota in the healthy oral cavity  and intestine . Leptotrichia buccalis was found in high prevalence in a study of the gingival crevice of Chinese patients with gingivitis and necrotizing ulcerative gingivitis . In a model of experimentally induced gingivitis, children harbored three-fold greater proportions of Leptotrichia species and 2.3-fold greater proportions of Selenomonas species in subgingival plaque than adults treated in the same way . Similarly, Selenomonas species are Gram negative anaerobes normally found in the buccal flora and associated with gingivitis [31, 32] and periodontitis [33, 34]. TM7 is a prominent bacterial phylum of over 200 phylotypes without cultivated representatives [35–37] and found in diverse environmental habitats (such as soil, freshwater, deep sea and hydrothermal vents). Members of the TM7 candidate division have been recently detected in various human body sites [6, 38–40], and associated with the inflammatory pathogenesis of several diseases (periodontitis , vaginosis  and inflammatory bowel diseases ). The subgroup I025 was found in subgingival plaque primarily at diseased sites in periodontitis hosts, suggesting their potential role in the multifactorial process leading to periodontitis [41, 44].
On the other hand, only five gingivitis-depleted genera were detected in the current study. Streptococcus is one of the most predominant genera in the human oral cavity. However, the "oral streptococci" are a highly heterogeneous group genetically . Although most are opportunistic pathogens and have been linked with a variety of oral diseases [46–48], they are also considered commensals. Similar to our results, Streptococcus sanguinis, as well as Lautropia mirabilis and Haemophilus parainfluenzae, were recently associated with oral health [34, 47]. The genus Veillonella represents a group of small, usually non-fermentative, strict anaerobic, Gram-negative cocci. They are found in the human oral cavity, the upper respiratory tract, small intestines and vagina. In a survey of subgingival plaque from 22 subjects, the majority of the subgingival Veillonella isolates were identified as Veillonella parvula . Prevotella species are part of the normal human oral microbiota and are frequently isolated from oral infections such as periodontitis, dental caries and abscesses [15, 50, 51]. Black-pigmenting members of Prevotella were associated with oral diseases. Consistently, in this study, most Prevotella OTUs detected in healthy hosts belonged to non-pigmenting species except Prevotella tannerae. Once validated in larger surveys, these gingivitis-associated genera, including both gingivitis-enriched and depleted ones, may represent valuable biomarkers for gingivitis.
Pyrosequencing techniques, such as the one employed in this study, revealed vast phylogenetic diversity and variability of bacterial communities in the human oral ecosystem [20, 52]. Characterization and quantification of community components allowed distinctions in community structure between healthy and diseased states to be explored for disease biomarkers and specific-microbe-targeted therapy. To our knowledge, this is the first organismal survey of gingivitis-associated microbiota using deep-sequencing techniques. Our preliminary findings formulate the basis for further studies that feature a longitudinal design and include a larger number of subjects. Ongoing technical improvements on phylogenetic-marker amplification (such as those targeting DNA-extraction bias, sequence chimerism caused by PCR, bias of PCR amplification, sequencing errors, unequal amplification of community members and the typically unknown variations in the rDNA-gene copy numbers among different residents [16–19, 53]) and the increasing coverage of oral 16 S rDNA reference databases  should allow the dissection of gingivitis-associated microbial factors at even higher sensitivity and resolution.
This study revealed that, first, microbiota from the four sampling sites for plaque (supragingival plaque and subgingival plaque from anterior teeth; supragingival plaque and subgingival plaque from posterior teeth; Methods) were similar to each other, yet were distinguishable from salivary microbiota. Second, community structures of plaque microbiota, but not saliva microbiota, can be used to distinguish gingivitis. Thus plaque should serve as the sampling site of choice in providing a microbial perspective for the disease. Third, a number of organisms were identified as gingivitis-associated (with several low-abundance gingivitis-specific genera detected; Table 3), which can serve as potential biomarkers. These results have important implications in the sampling and analysis strategies for surveying gingivitis-associated microbial risk factors. Our findings now enable further studies that examine the temporal development and epidemiology of microbial risk factors behind gingivitis. Furthermore, based on the gingivitis-associated microbiota identified in this study, an integrated organism- and gene-based survey of oral microbiomes at various clinically defined states should unravel the nature of microbial contribution to the development of gingivitis.
All oral samples were collected at the Hai Tai He Chang Clinical Research Center in Beijing with approval from P&G Beijing Technical Center (China) Institutional Review Board and in accordance with the World Medical Association Declaration of Helsinki (1996 amendment). ICH Guidelines for Good Clinical Practice (GCPs) were followed. Healthy subjects aged 18 years or older who had a minimum of 18 natural teeth were recruited from the Beijing area. Voluntary informed consent was provided. Individuals meeting the following criteria were excluded: current participation in another clinical study; use of antibiotic, anti-inflammatory or anticoagulant therapy within 30 days prior to examination; self-reported pregnancy or lactation; diabetes; a history of hepatitis or blood disorders such as hemophilia or leukocythemia; the presence of orthodontic appliances or removable partial dentures; significant oral pathology, such as advanced periodontal disease, hard or tissue tumors, or other conditions considered significant by the study director. Gingivitis was assessed using Mazza Gingival Index (MGI) as defined by Mazza in 1981 . Specifically, probing was performed by a dentist on the mesiobuccal and the distolingual sites of each tooth, for a maximum of 56 sites. BOP (Bleeding on probing) frequency and mean MGI were recorded for each subject. The MGI is similar to the Loe and Silness Gingival index; both are validated indices for describing gingivitis . The merit, however, of using MGI is that it combines measurements that address both the signs of inflammation as well as the degree of the severity of bleeding. Scores range from 0-5, with 0 assigned for normal appearing and healthy gingiva up to a score of 5 for spontaneous bleeding (without provocation). Five individuals with healthy gums and another five with extensive gingivitis were enrolled. Subjects were assigned to the healthy group (H) if there were less than ≤5 bleeding sites and to the unhealthy (gingivitis) group (U) when the frequency of bleeding sites was ≥ 20. No randomization among groups was performed. Two subjects from each group did not return for follow-up examinations and were excluded from further analyses. In the end, a total of six subjects (three in each group) completed the full study.
Samples of dental plaque and saliva were collected in the morning, 12 hours after evening tooth brushing. No oral hygiene or intakes of food and drink were allowed in the morning before sampling. Five samples were collected from each subject: supragingival dental plaque from anterior teeth (3~4 upper incisors), denoted A-sup; subgingival plaque [2 mm below gingival margin] from the same teeth, denoted A-sub; supragingival plaque from posterior teeth (2~3 upper molars), denoted P-sup; subgingival plaque from the same teeth, denoted P-sub; and saliva, denoted S. In the healthy group, plaque samples were collected from non-bleeding sites; in the unhealthy (gingivitis) group, incisor plaques were collected from non-bleeding sites and molar plaque were collected from bleeding sites. For unhealthy subjects, there are both non-bleeding incisor sites and bleeding molar sites for collecting plaque. We have considered the possibility that samples from the bleeding sites might not represent a complete picture of the microbiome of unhealthy gum. Therefore, plaque samples from both non-bleeding and bleeding sites were collected in our study.
Dental plaque samples were collected with sterile Gracey curettes and then removed from the curettes with a cotton-tipped swab. The tip of the swab was then placed into 0.4 mL lysis buffer (20 mM pH 8.0 Tris, 2 mM EDTA, 1.2% Triton X-100) and vortexed for 30 s. To collect salivary samples, subjects rinsed the mouth with 10 mL 0.9% saline buffer for 1 min and expectorated into a 50 ml tube. All samples were stored under -70°C before total genomic DNA extraction.
DNA extraction and PCR amplification
Bacterial pellets collected from dental plaque and saliva were suspended in lysis buffer with lysozyme (20 mg/ml) and incubated with proteinase K. Bacterial DNA was extracted using QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) following the manufacturer's instructions. PCR amplicon libraries of the small subunit ribosomal (16 S) RNA gene V1-V3 hypervariable region (Escherichia coli positions 5-534) were generated for each individual sample. PCR were performed using the forward primer (NNNNNNN-TGGAGAGTTTGATCCTGGCTCAG) and reverse primer (NNNNNNN-TACCGCGGCTGCTGGCAC). Unique heptad-nucleotide sequences (seven bases) were synthesized at 5' end of each pair of primers as barcodes, which helped to assign sequences to different samples.
The amplification mix contained 12.5 ul of Gotaq Hotstart polymerase 2 × mix (Promega, USA), a 1 ul of each primer (5 pM), 1 ul genomic DNA (0.1-10 ng μl-1) and 9.5 ul H2O in a total volume of 25 μl. Cycling conditions were an initial denaturation at 95°C for 2 min, 25 cycles at 94°C for 30 s, at 56°C for 25 s, and at 72°C for 25 s, followed by a final 5 minute extension at 72°C. Samples were processed via separate PCR reactions (ABI StepOnePlus™ Real-Time PCR Systems) and then pooled. Each sample was amplified using one specific barcoded primer. To assess quality, the PCR product for each sample was subjected to electrophoresis (1.2% agarose, 5 V cm-1, for 40 min). Gels were stained with a buffer containing SYBR Gold Nucleic Acid Gel Stain (Invitrogen, USA); DNA fragments of approximately 500 bp were excised from the gel and further purified using Qiagen MiniElute kit. Concentrations of DNA in purified PCR products were further analyzed with PicoGreen (Invitrogen, USA). The amplicons were pooled into a single tube in equimolar ratios. Pyrosequencing of the 16 S PCR-amplicons was carried out on Genome Sequencer FLX Titanium (Roche, USA) where, on average, 400 bp-long reads were produced.
The sequences generated from pyrosequencing were mainly analyzed with MOTHUR  for preprocessing, identification of operational taxonomic units (OTU), taxonomic assignment and community-structure comparisons. To minimize the effects of random sequencing errors and avoid overestimates of the phylogenetic diversity , relatively stringent quality-based trimming of the reads was performed. First, the 454-reads were removed if they were < 150 bp, had an average quality score < 35 in each 50-bp window rolling along the whole read, had an ambiguous base call (N), had any homopolymers of more than eight bases or did not contain the primer sequence; reads were then sorted by the tag sequences. To reduce sequencing noise from pyrosequencing data, we performed the pre-clustering step  with the "pre.cluster" script in MOTHUR . We also removed chimeric sequences detected by UCHIME .
Operational Taxonomical Units (OTU) assignment and taxonomic classification
The trimmed reads were assigned to clusters using UCLUST http://www.drive5.com/uclust/. An in-house perl script was used to convert UCLUST output into a format recognized by MOTHUR  http://www.mothur.org/ for further analysis. Reads were assigned to OTUs (species-level). Calculation of coverage percentage (Good), species richness estimators (ACE and Chao1) and rarefaction analysis were performed using MOTHUR . The relative abundance of OTUs with 97%-identity between pair-wise samples or between groups of samples were compared.
For taxonomic assignments, we used the "classify.seqs" script in MOTHUR  to classify all trimmed reads based on Naive Bayesian method with oral "CORE"  taxonomy sequences as the reference database. The confidence score threshold was set to 0.8, such that those with bootstrap value below 0.8 were assigned as unclassified. Relevant abundances of the bacterial taxa at the phylum and genus level were calculated and compared.
The OTUs defined by a 3% distance level were phylogenetically classified using the "classify.otu" script in MOTHUR  with oral "CORE" database  and a taxonomy file describing the complete taxonomic information of each sequence in the database from domain to species (using a 51% confidence threshold). The consensus taxonomy for each OTU was obtained in this step.
Comparisons of microbiota community structures
FastUnifrac-based community structure comparisons were performed . In each sample, representative sequences from each OTU were chosen by selecting the longest sequence based on UCLUST. Each sequence was assigned to its closest relative in a phylogeny of the Greengenes core set  using BLAST's megablast protocol. The resulting sample ID mapping file and category mapping file were used as inputs to the unweighted and weighted FastUniFrac . FastUniFrac allows pairwise comparisons of distances between two microbial communities in terms of the fraction of evolutionary history that separates the organisms. A distance (a measurement of the similarity in community structure between two microbiota) was computed for each pair of samples, both within a single population and across the two populations, to create a matrix of pairwise distances among all samples. These distances were then clustered to reduce dimensionality using PCoA . PCoA is a multivariate statistical technique for finding the most important axes along which the samples vary. The principal coordinates (PC), in descending order, describe of the degree of variation that each of the axes in the new space explains. ThetaYC-based community structure comparisons were performed in parallel with MOTHUR . ThetaYC () measures the dissimilarity between the structures of two communities , where S T is the total number of OTUs in communities A and B, a i is the relative abundance of OTU i in community A, b i is the relative abundance of OTU i in community B. A matrix of pairwise thetaYC-based distances among all samples was created for clustering and PCoA analysis.
Validation of 16 S rDNA pyrosequencing data by qPCR
Quantitative PCR assays on selected species were performed to test the degree of correlation with 16 S rDNA pyrosequencing data. Two genera, Streptococcus and Fusobacterium, were frequently identified based on our taxonomy assignments of the reads. Therefore, we chose two pairs of primers and probes targeting these two genera to perform the quantitative assays for comparisons to the pyrosequencing data.
Genus-specific primers and TaqMan probes were used, as listed in Additional file 4. The oligonucleotide probes were labeled with the fluorescent dyes 6-carboxyfluorescein (FAM) at the 5' end and 6-carboxytetramethylrhodamine (TAMRA) at the 3' end. The specificities of the probe and primer sets for their target DNA were tested in duplicate with the TaqMan Universal PCR Master Mix. The optimized concentrations of the forward primer, the reverse primer, and the fluorogenic probe in the 20-μl reaction volume were selected to be 300 nM, 300 nM, and 200 nM, respectively. Amplification and detection by quantitative PCR were performed with the StepOnePlus™ Real-Time PCR Systems (Applied Biosystems, Foster City, CA, USA). For each quantitative PCR, 20 μl reaction mixtures containing 2-μl sample DNA, forward primer, reverse primer and TaqMan probe at the optimized concentrations (as described above) were placed in each well of a 96-well plate. Following the fast TaqMan thermocycling protocol, reaction conditions were set at 95°C for 20 seconds, followed by 40 cycles of 95°C for 1 second and 58°C for 20 seconds. Standard curves for each organism were plotted in duplicate for each primer-probe set using the Ct (the cycle number at which the threshold fluorescence was reached) values, which were obtained by amplifying successive 10-fold dilutions of a known concentration of bacterial DNA (Streptococcus mutans UA159 and Fusobacterium nucleatum subsp. nucleatum ATCC25586). Copy-numbers of the target genes (tuf-elongation factor Tu and 16 S rDNA) in standard samples were calculated by the genome sizes (S. mutans 2.0 Mb and F. nucleatum 2.2 Mb) and the copy-number per genome (one copy of tuf gene per cell of S. mutans and five copies of 16 S rDNA gene per cell of F. nucleatum [46, 63]. One ng of S. mutans genomic DNA contains 4.63 × 105 copies of tuf gene while 1 ng F. nucleatum genome DNA contains 2.10 × 106 copies of 16 S rDNA gene. Based on these assumptions, the absolute copy number of a target gene was determined by referring Ct value to a standard cure measured on the same plate. The relative abundance of these bacteria in the 30 different oral specimens was normalized by the absolute quantity of DNA in the clinical samples.
AMOVA (Analysis of Molecular Variance) were used to test whether two communities from H and U populations have the same centroid [64, 65]. HOMOVA (Homogeneity of Molecular Variance) was employed to test whether the genetic diversity are similar between the communities from the H and U populations [65, 66].
Relative abundance of OTUs and phylotypes were reported as mean ± SEM. Due to the small sample sizes of these oral-site-specific datasets, features that are differentially distributed (i.e. abundant) between populations were statistically detected using Metastats  via a web interface http://metastats.cbcb.umd.edu/detection.html. Frequency data of OTUs and phylotypes were converted into a Feature Frequency Matrix as the input to this analysis tool. To exclude the extremely sparsely-sampled features (OTUs/phylotypes), tests were applied only if the total number of observations of a feature (OTU/phylotype) in either population is greater than the total number of subjects in the population (i.e. the average number of observations across subjects for a given feature is greater than one). Metastats was performed using 1000 permutations to compute p-values in statistical tests. We set p-value threshold of significance as 0.05. To control the FDR (False Discovery Rate) within the entire set of tests, we only took those features whose q-values and p-values were both below 0.05 into considerations. Levels of confidence were denoted as: *: 0.01 < p < 0.05; **: p < 0.01.
In validating the pyrosequencing results, the relative abundance of selected genera (Streptococcus and Fusobacterium) as measured via 16 S-amplicon pyrosequencing was compared to the corresponding gene copy number as determined by qPCR. The Shapiro-Wilk statistics of the variables were statistically significant. The degrees of correlation between the two measured parameters were determined from the Spearman's nonparametric correlation coefficient, r. Statistical analyses were performed with R (version 2.13.1). All reported p values were two-sided, at a 95% confidence level.
We thank Jay Tiesman, Kenton Juhlin, Jun Xu, Bruce Tepper and Tao He for their review and helpful criticism of the manuscript, Lily Sun for coordinating clinical study, Deborah Hutchins for technical writing and formatting and Xinping Cui for advice on statistical analysis. This work was partially funded by National High-Tech R&D Program of China (No. 2009 AA02 Z310).
- Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R: Bacterial community variation in human body habitats across space and time. Science. 2009, 326 (5960): 1694-1697. 10.1126/science.1177486.View ArticlePubMedPubMed CentralGoogle 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 (11): 5721-5732. 10.1128/JCM.43.11.5721-5732.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Keijser BJ, Zaura E, Huse SM, van der Vossen JM, Schuren FH, Montijn RC, ten Cate JM, Crielaard W: Pyrosequencing analysis of the oral microflora of healthy adults. J Dent Res. 2008, 87 (11): 1016-1020. 10.1177/154405910808701104.View ArticlePubMedGoogle Scholar
- Lazarevic V, Whiteson K, Huse S, Hernandez D, Farinelli L, Osteras M, Schrenzel J, Francois P: Metagenomic study of the oral microbiota by Illumina high-throughput sequencing. J Microbiol Methods. 2009, 79 (3): 266-271. 10.1016/j.mimet.2009.09.012.View ArticlePubMedPubMed CentralGoogle Scholar
- Zaura E, Keijser BJ, Huse SM, Crielaard W: Defining the healthy "core microbiome" of oral microbial communities. BMC Microbiol. 2009, 9: 259-10.1186/1471-2180-9-259.View ArticlePubMedPubMed CentralGoogle Scholar
- Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, Lakshmanan A, Wade WG: The human oral microbiome. J Bacteriol. 2010, 192 (19): 5002-5017. 10.1128/JB.00542-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Marsh PD: Microbial ecology of dental plaque and its significance in health and disease. Adv Dent Res. 1994, 8 (2): 263-271.PubMedGoogle Scholar
- Kuramitsu HK, He XS, Lux R, Anderson MH, Shi WY: Interspecies interactions within oral microbial communities. Microbiol Mol Biol Rev. 2007, 71 (4): 653-670. 10.1128/MMBR.00024-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Moore LVH, Moore WEC, Cato EP, Smibert RM, Burmeister JA, Best AM, Ranney RR: Bacteriology of Human Gingivitis. J Dent Res. 1987, 66 (5): 989-995. 10.1177/00220345870660052401.View ArticlePubMedGoogle Scholar
- Handfield M, Baker HV, Lamont RJ: Beyond good and evil in the oral cavity: Insights into host-microbe relationships derived from transcriptional profiling of gingival cells. J Dent Res. 2008, 87 (3): 203-223. 10.1177/154405910808700302.View ArticlePubMedPubMed CentralGoogle Scholar
- Offenbacher S, Barros SP, Paquette DW, Winston JL, Biesbrock AR, Thomason RG, Gibb RD, Fulmer AW, Tiesman JP, Juhlin KD, et al: Gingival transcriptome patterns during induction and resolution of experimental gingivitis in humans. J Periodontol. 2009, 80 (12): 1963-1982. 10.1902/jop.2009.080645.View ArticlePubMedGoogle Scholar
- Slots J, Moenbo D, Langebaek J, Frandsen A: Microbiota of gingivitis in man. Scand J Dent Res. 1978, 86 (3): 174-181.PubMedGoogle Scholar
- Savitt ED, Socransky SS: Distribution of certain subgingival microbial species in selected periodontal conditions. J Periodontal Res. 1984, 19 (2): 111-123. 10.1111/j.1600-0765.1984.tb00800.x.View ArticlePubMedGoogle Scholar
- Kroes I, Lepp PW, Relman DA: Bacterial diversity within the human subgingival crevice. Proc Natl Acad Sci USA. 1999, 96 (25): 14547-14552. 10.1073/pnas.96.25.14547.View ArticlePubMedPubMed CentralGoogle Scholar
- Paster BJ, Boches SK, Galvin JL, Ericson RE, Lau CN, Levanos VA, Sahasrabudhe A, Dewhirst FE: Bacterial diversity in human subgingival plaque. J Bacteriol. 2001, 183 (12): 3770-3783. 10.1128/JB.183.12.3770-3783.2001.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashelford KE, Chuzhanova NA, Fry JC, Jones AJ, Weightman AJ: At least 1 in 20 16 S rRNA sequence records currently held in public repositories is estimated to contain substantial anomalies. Appl Environ Microbiol. 2005, 71 (12): 7724-7736. 10.1128/AEM.71.12.7724-7736.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu Z, DeSantis TZ, Andersen GL, Knight R: Accurate taxonomy assignments from 16 S rRNA sequences produced by highly parallel pyrosequencers. Nucleic Acids Res. 2008, 36 (18): e120-10.1093/nar/gkn491.View ArticlePubMedPubMed CentralGoogle Scholar
- Petrosino JF, Highlander S, Luna RA, Gibbs RA, Versalovic J: Metagenomic pyrosequencing and microbial identification. Clin Chem. 2009, 55 (5): 856-866. 10.1373/clinchem.2008.107565.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang X, Xie L, Li YX, Wei CC: More than 9,000,000 Unique Genes in Human Gut Bacterial Community: Estimating Gene Numbers Inside a Human Body. PloS One. 2009, 4 (6):
- Hamady M, Lozupone C, Knight R: Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 2010, 4 (1): 17-27. 10.1038/ismej.2009.97.View ArticlePubMedGoogle Scholar
- Griffen AL, Beall CJ, Firestone ND, Gross EL, Difranco JM, Hardman JH, Vriesendorp B, Faust RA, Janies DA, Leys EJ: CORE: a phylogenetically-curated 16 S rDNA database of the core oral microbiome. PloS One. 2011, 6 (4): e19051-10.1371/journal.pone.0019051.View ArticlePubMedPubMed CentralGoogle Scholar
- Teles RP, Bogren A, Patel M, Wennstrom JL, Socransky SS, Haffajee AD: A three-year prospective study of adult subjects with gingivitis II: microbiological parameters. J Clin Periodontol. 2007, 34 (1): 7-17. 10.1111/j.1600-051X.2006.01015.x.View ArticlePubMedGoogle Scholar
- Mans JJ, von Lackum K, Dorsey C, Willis S, Wallet SM, Baker HV, Lamont RJ, Handfield M: The degree of microbiome complexity influences the epithelial response to infection. BMC Genomics. 2009, 10: 380-10.1186/1471-2164-10-380.View ArticlePubMedPubMed CentralGoogle Scholar
- Lazarevic V, Whiteson K, Hernandez D, Francois P, Schrenzel J: Study of inter- and intra-individual variations in the salivary microbiota. BMC Genomics. 2010, 11: 523-10.1186/1471-2164-11-523.View ArticlePubMedPubMed CentralGoogle Scholar
- Crielaard W, Zaura E, Schuller AA, Huse SM, Montijn RC, Keijser BJ: Exploring the oral microbiota of children at various developmental stages of their dentition in the relation to their oral health. BMC medical genomics. 2011, 4: 22-10.1186/1755-8794-4-22.View ArticlePubMedPubMed CentralGoogle Scholar
- Nasidze I, Li J, Quinque D, Tang K, Stoneking M: Global diversity in the human salivary microbiome. Genome Res. 2009, 19 (4): 636-643. 10.1101/gr.084616.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Kononen E, Asikainen S, Saarela M, Karjalainen J, Jousimies-Somer H: The oral gram-negative anaerobic microflora in young children: longitudinal changes from edentulous to dentate mouth. Oral Microbiol Immunol. 1994, 9 (3): 136-141. 10.1111/j.1399-302X.1994.tb00049.x.View ArticlePubMedGoogle Scholar
- Eribe ER, Olsen I: Leptotrichia species in human infections. Anaerobe. 2008, 14 (3): 131-137. 10.1016/j.anaerobe.2008.04.004.View ArticlePubMedGoogle Scholar
- Gmur R, Wyss C, Xue Y, Thurnheer T, Guggenheim B: Gingival crevice microbiota from Chinese patients with gingivitis or necrotizing ulcerative gingivitis. Eur J Oral Sci. 2004, 112 (1): 33-41. 10.1111/j.0909-8836.2004.00103.x.View ArticlePubMedGoogle Scholar
- Moore WE, Holdeman LV, Smibert RM, Cato EP, Burmeister JA, Palcanis KG, Ranney RR: Bacteriology of experimental gingivitis in children. Infect Immun. 1984, 46 (1): 1-6.PubMedPubMed CentralGoogle Scholar
- van Palenstein Helderman WH: Total viable count and differential count of vibrio (campylobacter) sputorum, fusobacterium nucleatum, selenomonas sputigena, bacteroides ochraceus and veillonella in the inflamed and non inflamed human gingival crevice. J Periodontal Res. 1975, 10 (5): 294-305. 10.1111/j.1600-0765.1975.tb00037.x.View ArticlePubMedGoogle Scholar
- Bisiaux-Salauze B, Perez C, Sebald M, Petit JC: Bacteremias caused by Selenomonas artemidis and Selenomonas infelix. J Clin Microbiol. 1990, 28 (1): 140-142.PubMedPubMed CentralGoogle Scholar
- Drescher J, Schlafer S, Schaudinn C, Riep B, Neumann K, Friedmann A, Petrich A, Gobel UB, Moter A: Molecular epidemiology and spatial distribution of Selenomonas spp. in subgingival biofilms. Eur J Oral Sci. 2010, 118 (5): 466-474. 10.1111/j.1600-0722.2010.00765.x.View ArticlePubMedGoogle Scholar
- Colombo AP, Boches SK, Cotton SL, Goodson JM, Kent R, Haffajee AD, Socransky SS, Hasturk H, Van Dyke TE, Dewhirst F, et al: Comparisons of subgingival microbial profiles of refractory periodontitis, severe periodontitis, and periodontal health using the human oral microbe identification microarray. J Periodontol. 2009, 80 (9): 1421-1432. 10.1902/jop.2009.090185.View ArticlePubMedPubMed CentralGoogle Scholar
- Rappe MS, Giovannoni SJ: The uncultured microbial majority. Annu Rev Microbiol. 2003, 57: 369-394. 10.1146/annurev.micro.57.030502.090759.View ArticlePubMedGoogle Scholar
- Marcy Y, Ouverney C, Bik EM, Losekann T, Ivanova N, Martin HG, Szeto E, Platt D, Hugenholtz P, Relman DA, et al: Dissecting biological "dark matter" with single-cell genetic analysis of rare and uncultivated TM7 microbes from the human mouth. Proc Natl Acad Sci USA. 2007, 104 (29): 11889-11894. 10.1073/pnas.0704662104.View ArticlePubMedPubMed CentralGoogle Scholar
- Dinis JM, Barton DE, Ghadiri J, Surendar D, Reddy K, Velasquez F, Chaffee CL, Lee MC, Gavrilova H, Ozuna H, et al: In search of an uncultured human-associated TM7 bacterium in the environment. Plos One. 2011, 6 (6): e21280-10.1371/journal.pone.0021280.View ArticlePubMedPubMed CentralGoogle Scholar
- Gao Z, Tseng CH, Pei Z, Blaser MJ: Molecular analysis of human forearm superficial skin bacterial biota. Proc Natl Acad Sci USA. 2007, 104 (8): 2927-2932. 10.1073/pnas.0607077104.View ArticlePubMedPubMed CentralGoogle Scholar
- Bik EM, Long CD, Armitage GC, Loomer P, Emerson J, Mongodin EF, Nelson KE, Gill SR, Fraser-Liggett CM, Relman DA: Bacterial diversity in the oral cavity of 10 healthy individuals. ISME J. 2010, 4 (8): 962-974. 10.1038/ismej.2010.30.View ArticlePubMedPubMed CentralGoogle Scholar
- Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308 (5728): 1635-1638. 10.1126/science.1110591.View ArticlePubMedPubMed CentralGoogle Scholar
- Brinig MM, Lepp PW, Ouverney CC, Armitage GC, Relman DA: Prevalence of bacteria of division TM7 in human subgingival plaque and their association with disease. Appl Environ Microbiol. 2003, 69 (3): 1687-1694. 10.1128/AEM.69.3.1687-1694.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Fredricks DN, Fiedler TL, Marrazzo JM: Molecular identification of bacteria associated with bacterial vaginosis. N Engl J Med. 2005, 353 (18): 1899-1911. 10.1056/NEJMoa043802.View ArticlePubMedGoogle Scholar
- Kuehbacher T, Rehman A, Lepage P, Hellmig S, Folsch UR, Schreiber S, Ott SJ: Intestinal TM7 bacterial phylogenies in active inflammatory bowel disease. J Med Microbiol. 2008, 57 (Pt 12): 1569-1576.View ArticlePubMedGoogle Scholar
- Ouverney CC, Armitage GC, Relman DA: Single-cell enumeration of an uncultivated TM7 subgroup in the human subgingival crevice. Appl Environ Microbiol. 2003, 69 (10): 6294-6298. 10.1128/AEM.69.10.6294-6298.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Itzek A, Gillen CM, Fulde M, Friedrichs C, Rodloff AC, Chhatwal GS, Nitsche-Schmitz DP: Contribution of plasminogen activation towards the pathogenic potential of oral streptococci. PloS One. 2010, 5 (11): e13826-10.1371/journal.pone.0013826.View ArticlePubMedPubMed CentralGoogle Scholar
- Picard FJ, Ke D, Boudreau DK, Boissinot M, Huletsky A, Richard D, Ouellette M, Roy PH, Bergeron MG: Use of tuf sequences for genus-specific PCR detection and phylogenetic analysis of 28 streptococcal species. J Clin Microbiol. 2004, 42 (8): 3686-3695. 10.1128/JCM.42.8.3686-3695.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang W, Jiang Y, Li C, Liang J: Investigation of Supragingival Plaque Microbiota in Different Caries Status of Chinese Preschool Children by Denaturing Gradient Gel Electrophoresis. Microb Ecol. 2011, 61 (2): 342-352. 10.1007/s00248-010-9753-z.View ArticlePubMedGoogle Scholar
- Ling Z, Kong J, Jia P, Wei C, Wang Y, Pan Z, Huang W, Li L, Chen H, Xiang C: Analysis of oral microbiota in children with dental caries by PCR-DGGE and barcoded pyrosequencing. Microb Ecol. 2010, 60 (3): 677-690. 10.1007/s00248-010-9712-8.View ArticlePubMedGoogle Scholar
- Leuckfeld I, Paster BJ, Kristoffersen AK, Olsen I: Diversity of Veillonella spp. from subgingival plaque by polyphasic approach. APMIS. 2010, 118 (3): 230-242. 10.1111/j.1600-0463.2009.02584.x.View ArticlePubMedGoogle Scholar
- Dymock D, Weightman AJ, Scully C, Wade WG: Molecular analysis of microflora associated with dentoalveolar abscesses. J Clin Microbiol. 1996, 34 (3): 537-542.PubMedPubMed CentralGoogle Scholar
- Nadkarni MA, Caldon CE, Chhour KL, Fisher IP, Martin FE, Jacques NA, Hunter N: Carious dentine provides a habitat for a complex array of novel Prevotella-like bacteria. J Clin Microbiol. 2004, 42 (11): 5238-5244. 10.1128/JCM.42.11.5238-5244.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Hamady M, Knight R: Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009, 19 (7): 1141-1152. 10.1101/gr.085464.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Roh C, Villatte F, Kim BG, Schmid RD: Comparative study of methods for extraction and purification of environmental DNA from soil and sludge samples. Appl Biochem Biotechnol. 2006, 134 (2): 97-112. 10.1385/ABAB:134:2:97.View ArticlePubMedGoogle Scholar
- Sundquist A, Bigdeli S, Jalili R, Druzin ML, Waller S, Pullen KM, El-Sayed YY, Taslimi MM, Batzoglou S, Ronaghi M: Bacterial flora-typing with targeted, chip-based Pyrosequencing. BMC Microbiol. 2007, 7: 108-10.1186/1471-2180-7-108.View ArticlePubMedPubMed CentralGoogle Scholar
- Mazza JE, Newman MG, Sims TN: Clinical and antimicrobial effect of stannous fluoride on periodontitis. J Clin Periodontol. 1981, 8 (3): 203-212. 10.1111/j.1600-051X.1981.tb02031.x.View ArticlePubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75 (23): 7537-7541. 10.1128/AEM.01541-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Kunin V, Engelbrektson A, Ochman H, Hugenholtz P: Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. Environ Microbiol. 2010, 12 (1): 118-123. 10.1111/j.1462-2920.2009.02051.x.View ArticlePubMedGoogle Scholar
- Huse SM, Welch DM, Morrison HG, Sogin ML: Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol. 2010, 12 (7): 1889-1898. 10.1111/j.1462-2920.2010.02193.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27 (16): 2194-2200. 10.1093/bioinformatics/btr381.View ArticlePubMedPubMed CentralGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16 S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72 (7): 5069-5072. 10.1128/AEM.03006-05.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71 (12): 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Yue JC, Clayton MK, Lin FC: A nonparametric estimator of species overlap. Biometrics. 2001, 57 (3): 743-749. 10.1111/j.0006-341X.2001.00743.x.View ArticlePubMedGoogle Scholar
- Martin FE, Nadkarni MA, Jacques NA, Hunter N: Quantitative microbiological study of human carious dentine by culture and real-time PCR: association of anaerobes with histopathological changes in chronic pulpitis. J Clin Microbiol. 2002, 40 (5): 1698-1704. 10.1128/JCM.40.5.1698-1704.2002.View ArticlePubMedPubMed CentralGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMedPubMed CentralGoogle Scholar
- Schloss PD: Evaluating different approaches that test whether microbial communities have the same structure. ISME J. 2008, 2 (3): 265-275. 10.1038/ismej.2008.5.View ArticlePubMedGoogle Scholar
- Stewart CN, Excoffier L: Assessing population genetic structure and variability with RAPD data: Application to Vaccinium macrocarpon (American Cranberry). J Evol Biol. 1996, 9 (2): 153-171. 10.1046/j.1420-9101.1996.9020153.x.View ArticleGoogle Scholar
- White JR, Nagarajan N, Pop M: Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009, 5 (4): e1000352-10.1371/journal.pcbi.1000352.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6831/11/33/prepub