- Research article
- Open Access
Microbiota diversity and gene expression dynamics in human oral biofilms
- Alfonso Benítez-Páez†1, 2Email author,
- Pedro Belda-Ferre†1,
- Aurea Simón-Soro1 and
- Alex Mira1Email author
© Benítez-Páez et al.; licensee BioMed Central Ltd. 2014
- Received: 22 September 2013
- Accepted: 10 April 2014
- Published: 27 April 2014
Micro-organisms inhabiting teeth surfaces grow on biofilms where a specific and complex succession of bacteria has been described by co-aggregation tests and DNA-based studies. Although the composition of oral biofilms is well established, the active portion of the bacterial community and the patterns of gene expression in vivo have not been studied.
Using RNA-sequencing technologies, we present the first metatranscriptomic study of human dental plaque, performed by two different approaches: (1) A short-reads, high-coverage approach by Illumina sequencing to characterize the gene activity repertoire of the microbial community during biofilm development; (2) A long-reads, lower-coverage approach by pyrosequencing to determine the taxonomic identity of the active microbiome before and after a meal ingestion. The high-coverage approach allowed us to analyze over 398 million reads, revealing that microbial communities are individual-specific and no bacterial species was detected as key player at any time during biofilm formation. We could identify some gene expression patterns characteristic for early and mature oral biofilms. The transcriptomic profile of several adhesion genes was confirmed through qPCR by measuring expression of fimbriae-associated genes. In addition to the specific set of gene functions overexpressed in early and mature oral biofilms, as detected through the short-reads dataset, the long-reads approach detected specific changes when comparing the metatranscriptome of the same individual before and after a meal, which can narrow down the list of organisms responsible for acid production and therefore potentially involved in dental caries.
The bacteria changing activity during biofilm formation and after meal ingestion were person-specific. Interestingly, some individuals showed extreme homeostasis with virtually no changes in the active bacterial population after food ingestion, suggesting the presence of a microbial community which could be associated to dental health.
- Dental plaque
- Biofilm formation
- Human microbiome
The study of microbial communities from environment- and human-derived samples through Next Generation Sequencing (NGS) methods has revealed a vast complexity in those ecological niches where hundreds or thousands of microbial species co-inhabit and functionally interact. One of these complex communities is that found in the human oral dental plaque (hereinafter, human oral biofilm). Although some studies, using NGS methods and 16S rRNA-based analysis, estimate that microbial diversity of the oral cavity is composed by thousands of species , more recent data have limited these estimates to a few hundreds [2–4]. Contrary to Koch's postulates, dental caries is not considered etiologically the outcome of a single-agent but is associated to an unbalance of microbial species that synergistically cause enamel demineralization by their acidogenic activity [5, 6]. Thus, characterizing the composition of whole bacterial communities that actively engage in biofilm formation and sugar fermentation after the ingestion of food is vital for understanding community dynamics under health and disease conditions .
Although the set of species present in the human oral biofilm is almost fully depicted, new efforts have to be conducted to establish microbial agonistic or antagonistic associations, to distinguish actively-growing bacteria from inactive or transient species, as well as to outline the role of individual species during biofilm formation on tooth surfaces. The co-aggregation detected to occur between streptococci and Actinomyces species has been proposed to be a major promoter of human oral biofilm formation . Like most biofilms, the dental plaque is built in a continued process characterized by succession of different bacterial species, each one with relevant roles in every step of biofilm construction . Formation of the oral biofilm could be dissected in three major stages, namely: i) attachment; ii) colonization; and iii) biofilm development . However, species participating of the entire process are traditionally characterized as “early” and “late” colonizers, where early colonizers would be responsible of the two first stages . Among early colonizers the viridans streptococci group is considered as a cornerstone of the oral biofilm puzzle given its ability to bind saliva proteins through Antigens I and II. In this manner, streptococci species become the first colonizers able to bind tooth surfaces and promoting arrival of secondary colonizers by intergeneric coaggregation (reviewed in ). Actinomyces naeslundii is one of the secondary colonizers and a well known coaggregation partner of streptococci [8, 11]. Fusobacterium nucleatum is considered a key player given its capability to coaggregate both with early and late colonizers of the oral biofilm , the latter group characterized by species belonging to Bacteroidetes and Spirochaetes [6, 9]. It is noteworthy to highlight that intergeneric coaggregation not only contributes to bacterial growth and colonization , but it is thought to facilitate the genetic and metabolic exchange among species, and even to create the adequate environment for arrival of some obligate anaerobic bacteria . Therefore, any disruption in the development of the oral biofilm caused by impairing of early colonizers tooth attachment or inability to recruit other key players during biofilm formation, would affect the entire process avoiding presence of pathogens responsible for periodontal disease or caries .
Although few attempts to link specific gene expression profiles in oral bacteria with the establishment and maturation of oral biofilm have been done , further studies are needed to understand global gene dynamics and intracellular signalling which are the basis for cell-to-cell communication among oral bacteria and to promote biofilm formation on tooth surfaces. There are important limitations to study gene expression from in vivo oral samples, including RNA instability and amounts of sampling material, but a sequencing approach of total cDNA from an in vitro oral biofilm model has recently been performed .
Because gene transcripts typically occupy a small fraction of total bacterial RNA, even after mRNA-enriching protocols, a massive coverage is normally required to quantify gene expression by RNAseq technologies. On the other hand, high-coverage sequencing technologies are normally coupled to short read lengths, which jeopardize accurate taxonomic assignment of the sequences. The latter can be achieved through the use of longer reads, at the expense of a lower coverage of mRNA transcripts. In the present manuscript we present the first metatranscriptome analysis of in vivo human oral biofilm samples through two approaches: A short read-length, high coverage Illumina® approach to study oral biofilm formation through time, and a long read-length, lower coverage pyrosequencing strategy to study changes in community composition before and after a meal. For the first approach, a total of 16 samples of supragingival plaque from 4 healthy individuals were collected at four different time points (6, 12, 24 and 48 hours after a professional ultrasound cleaning) to disclose the microbiota and gene expression dynamics during oral biofilm formation. For the second experiment, the metatranscriptome of dental plaque from five individuals was studied 30 minutes before and after a controlled meal, in order to characterize the potential shifts in the active bacterial community when dietary nutrients are available for growth.
Metagenome vs metatranscriptome
Low-coverage, long-reads approach
Active communities before and after a meal
Actinomyces was the only genus found at a proportion over 10% in all samples and was found to be significantly more abundant in healthy individuals (Figure 2B) according to a high-dimensional class comparison test . Actinomyces is an early colonizer of the oral biofilm, usually present in the first layers of dental plaque in contact with enamel . It is able to increase local pH by producing ammonia through the degradation of arginine, lysilarginine and urea . As a consequence, individual Actinomyces species could represent biomarkers of healthy biofilms with a protective role against acidogenic bacteria . On the other hand, late colonizers being strictly anaerobes like Porphyromonas, Fusobacterium, Capnocytophaga, Tannerella and Leptotrichia were found significantly more abundant in oral biofilm from caries-bearing individuals (Figure 2B). This should be further studied given the observed over-representation of species belonging to the red complex of periodontal disease .
Although the number of assigned mRNA transcripts detected by this approach is quite small, we could identify some genes being expressed in all samples such as those encoding ribosomal proteins and basic housekeeping machinery like elongation factors and cell division proteins. Interestingly, we found expression of multiple sugar transport systems and central metabolism genes such as glyceraldehyde 3-phosphate dehydrogenase, L-lactate dehydrogenase, citrate synthase, enolase, and malate dehydrogenase.
High-coverage, short-reads approach
Microbial activity during oral biofilm formation
Although the high-coverage approach (~25 million Illumina reads per sample) was mainly aimed at determining gene expression patterns during biofilm formation, an attempt was first made to utilize the reads mapping to rRNA genes to characterize microbiota composition at different times of dental plaque formation. Potential errors in taxonomic assignment were minimized by 1) assigning reads at the family and genus taxonomic level only; 2) selecting matches against 16S rRNA database of 100% sequence identity; and 3) eliminating hits to conserved regions of the 16S rRNA gene, keeping only hypervariable, informative regions.
Microbial interactions during oral biofilm formation
In vivo gene expression and functional analysis of the oral biofilm
For the functional analysis of gene expression patterns during oral biofilm formation the KEGG functional classification  was used. We detected expression of 19,519 genes (ORFs with aligned reads) on average per sample and the set of KO categories represented was 2,266, indicating that ~ 12% of expressed genes were functionally annotated. A distinguishable clustering according to biofilm stage was not fully depicted, although a pairwise clustering among early and mature biofilm samples was observed (data not shown). Some molecular pathways were differentially expressed between early and mature biofilm (i.e. Ribosome; Purine Metabolism; and Glycolysis). A comparison was performed between early (6 h and 12 h) and mature (24 h and 48 h) biofilm samples, determining the False Discovery Rate (FDR) with q-value ≤ 0.05. We detected a set of 271 KO categories differentially expressed (35 over-expressed in early biofilm and 236 in mature biofilm).
Over-expression in early biofilm
Over-expression of KO categories in early biofilm was predominantly grouped in genes involved in the metabolism of Carbohydrates, Energy, Amino Acids, Cofactor/Vitamins, and Xenobiotic Degradation. Translation functions were also overrepresented because several ribosome proteins showed higher expression during early biofilm as well as the Elongation Factors Tu (EF-Tu, K02358) and G (EF-G, K02355). Central role of Translation in early steps of oral biofilm formation was also evidenced by over-expression of K00566 category corresponding to the MnmA tRNA-modifying protein. MnmA is an evolutionarily conserved enzyme and incorporates the posttranscriptional modification s2U at the wobble position of several tRNAs  which reads A/G ending codons during translation. The tRNA modifications are largely associated to control the fine-tuning of protein synthesis, thus improving ribosome accuracy . As a consequence, they appear to be involved in controlling a wide range of bacterial phenotypic traits including biofilm formation by multi-drug resistant human pathogens  and could be important for oral biofilm formation as well. Our results are in agreement with previous reports where amino acid metabolism is critical for growth of early colonizers such as Streptococcus gordonii, which needs coaggregation to stabilize expression of genes involved in amino acid synthesis and membrane transporters .
Over-expression in mature biofilm
On the other hand, genes over-expressed at the late biofilm stage had a more variable functional profile. New functional categories over-expressed in late oral biofilms included the ABC transporters, the Cell motility represented by the orthology groups K03407 associated to bacterial chemotaxis, K02676 and K02390 involved in pilus and flagella assembly, respectively, and finally genes involved in Base Excision Repair, Mismatch Repair, and Homologous Recombination systems; and some putative transposases. The group of tRNA-modifying enzymes was also present in late oral biofilm with a larger set of genes over-expressed such as mnmC, yfiC, cmoA, tadA, queE, trmK, and the hydrouridine synthase gene dusC. Consequently, this molecular pathway seems to be involved in controlling the expression of proteins along the full oral biofilm process. Strikingly, a set of genes involved in competition between bacterial species also showed over-expression during this late stage of oral biofilm formation. The comFC, comFA, comGB, comGC, and comGA orthology groups are members of the Type II secretion system of which other members appear to be annotated as Competence-related DNA transformation transporters. These genes have been reported to be involved in quorum sensing response to produce mutacins; these are non-lantibiotic bacteriocins able to induce lysis and consequently DNA liberation in related species possibly supporting DNA horizontal transfer . Likewise, tfoX orthologues were also found to be over-expressed in late oral biofilm, thus strengthening the idea of natural genetic transformation [34, 35] occurring among close species in the mature oral biofilm. Globally, over-expression of these competence-related genes, permitting DNA transformations in vivo, could support the specific low ratio between functional diversity of genes and operational taxonomic units detected in supragingival plaque, thus indicating high functional redundancy and microbial population homogenization . Other important functional categories over-expressed in late oral biofilm included those involved in Environmental Information Processing and membrane transporters, such as those belonging to the Phosphotransferase System (PTS) as well as MFS membrane receptors specialized in the importing/exporting of small molecules. Both major families of membrane transporters were found to be preferentially expressed from Actinomyces species indicating a high level of metabolic exchange between this genus and its environment. However, subfamilies such a salX-like ABC transporters associated to bacteriocin export and defense were detected to be predominantly active in streptococci species. Over-expression of some KEGG orthology categories belonging to Two-Component family of proteins indicate an active role of cells in perceiving external signals of nutrient availability in the environment. An over-expression was found of the PTS-Ntr-EIIA enzyme and the GlnB protein, both involved in nitrogen regulation, and the sigma factor 54 of the RNA polymerase involved in expression of genes for nitrogen metabolism. Therefore, processes related to nitrogen uptake/metabolism appear to be very relevant in the mature stage of oral biofilm probably indicating that nitrogen is a limiting factor for oral biofilm progression. In recent studies of cDNA massive sequencing from an in vitro five-species oral biofilm microbial community, similar results were obtained in terms of overrepresented functions in mature biofilms . Finally, luxS homologue in Neisseria spp. was found to be significantly over-expressed in early biofilm. The luxS genes are responsible of Autoinducer-2 (AI-2) synthesis, a molecule considered as a major interspecies signal for cell-cell communication [9, 10]. Evidence for AI-2 role to control biofilm formation was previously observed when a luxS null strain of S. gordonii was unable to form a mixed-species biofilm with P. gingivalis.
Gene expression by qPCR
Our study shows for the first time the microbial diversity and gene expression dynamics in the complex oral microbial community in vivo. We could follow oral biofilm formation and determine proportions of active microbiota through time, including before and after a carbohydrate-rich meal, when the process of acid production, responsible of enamel demineralization, takes place. We present a large set of correlations among bacterial groups and genera being in agreement with biological and classical interactions reported to be central for biofilm installation and development [8, 48, 49]. In the functional exploration of genes expressed during human oral biofilm formation, we present a quantitative analysis, further supported by results obtained by qPCR, demonstrating several functional categories of prevalence at different oral biofilm stages. Among them we showed that translation machinery is predominantly expressed in early biofilm stages whereas more specialized genes are required in mature biofilm. Some genes involved in competence, and reported to be involved in quorum sensing response and functionally related to mutacin production and DNA uptake, were over-expressed in late biofilm supporting the intricate level of cell-to-cell interactions in mature biofilm and suggesting strong competition for colonization. More than 70% of the genetic information compiled from this oral metatranscriptome has no functional assignment; therefore, further efforts must be conducted for classification and characterization of genes and their involvement in biofilm development and/or cell-to-cell communication. From an applied point of view, the identification of active bacterial species after food uptake can be considered a first step to narrow down the list of potential etiological agents of dental caries from the large set of micro-organisms found in the metagenome of dental plaque and cavities . The striking homeostasis found in one of the individuals who had never suffered from dental caries, and where virtually no changes were found in the active microbiota before and after a meal, could indicate that the microbiota of some individuals is not affected by food ingestion, potentially reducing the risk of acidic pH and promoting dental health.
Sample collection and RNA processing
The sampling procedure was approved by the Ethical Committee for Clinical Research from the DGSP-CSISP (Valencian Health Authority, Spain) and all donors signed an informed consent. The oral health status of each individual was evaluated before sampling and following recommendations and nomenclature from the WHO. Donors were 20-30 years of age, had all 28 teeth present (excluding third molars) had not suffered from any systemic disease and had not taken systemic antimicrobials in the previous 6 months. Dental plaque samples were taken with autoclaved spoon excavators from vestibular and lingual surfaces of teeth excluding a 1 mm region on the edges.
For the biofilm formation experiments, 16 supragingival dental plaque samples were obtained from 4 caries-free volunteers (DMFT = 0 [decayed, Missing, Filling Teeth], OHI = 1 [oral Hygiene Index]. GI = 1 [gingival Index]). The volunteers were subjected to professional teeth ultrasound cleaning. Oral biofilm (supragingival plaque) from all teeth surfaces was pooled and collected from each volunteer at 6, 12, 24, and 48 h of biofilm formation. After every sampling a professional brushing was performed to reset biofilm formation for next sampling. Total RNA was extracted using the MasterPureTM RNA Purification Kit (Epicentre®). Samples were collected and processed for elimination of 5S rRNA and tRNAs through ion exchange chromatography with KCl gradient in Nucleobond AX 20 columns (Macherey-Nagel). Pre- and post-processed RNA were loaded in RNA chip (Agilent Technologies) and analyzed for integrity using Agilent Bioanalyzer 2100 (Agilent Technologies). The first strand of cDNA from processed RNA was synthesized using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). For this aim, two cDNA reactions were prepared for each RNA sample and modifying some manufacturer’s instructions to obtain best performance during synthesis. Specifically, each cDNA reaction had 100U of Multiscribe Reverse transcriptase and synthesis was completed at 48°C during 210 min using 5-10 ug of RNA as template. Doubled stranded cDNA (ds-cDNA) was achieved in 100 uL of reaction containing 35U E.coli DNA Polymerase I (New England Biolabs), 5U E. coli DNA Ligase (New England Biolabs), 5U RNase H (Epicentre®), 300 uM dNTP's, and two reactions of first strand cDNA synthesis. The ds-cDNA synthesis was initiated by incubation during 150 min at 16°C and completed by adding 4.5U of T4 DNA Polymerase (New England Biolabs), and 1X BSA (New England Biolabs) followed by incubation during additional 30 min at 16°C. Purified ds-cDNAs were obtained using High Pure PCR Product Purification Kit (Roche®) and sent to GATC Biotech AG (Konstanz, Germany) for parallel single-end sequencing using HiSeq2000 system (Illumina®).
The five donors for the before/after meal transcriptome were asked not to brush their teeth for 16 hours. Three of them had active caries at the moment of sampling (Decayed Teeth = 3, OHI = 1, GI = 1) and the other two had no history of dental caries (DMFT = 0, OHI = 1, GI = 1). None of the donors had periodontal disease. All donors ingested the same meal, whose nutritional characteristics are indicated in Additional file 4: Table S1. Supragingival dental plaque was obtained from the right maxillary and left mandibular quadrant free teeth surfaces for the sample 30 minutes before eating and from the left maxillary and right mandibular quadrant 30 minutes after food intake, without touching caries lesions, if they were present. Opposite quadrants were sampled because, when analysing PCR-amplified 16S rRNA from cDNA, an equivalence in terms of taxonomic composition was found when sampling opposite mandibular and maxillary quadrants (Additional file 5: Figure S3). The obtained total ds-cDNA (as described above) was purified and enriched in fragments longer than 400 bp using AMPure beads (Agencourt). Those long cDNA fragments were sequenced using 454 GS-FLX technology with titanium chemistry (Roche).
Taxonomic assignment and correlations
Filtering and trimming of original data set was assisted by Galaxy Web Server [16–18], filtering by quality using the sliding window method (window size 25 with a minimum quality in the window of 20), and sequences shorter than 200 bp were removed. For the high-coverage biofilm samples, microbial diversity was established by taxonomic assignment using reads matching 16S rRNA sequences. For this aim we constructed a RDP-based (Release 10, Update 29) database containing almost 10,000 reference sequences of 16S rDNA annotated according to NCBI taxonomy . This reference database was processed to filter out the conserved regions of 16S rDNA genes using Hidden Markov Models . Then, using MegaBlast v2.2.21 algorithm  and selecting alignments for 48 nt in length and 100% identity we could assign taxonomy at genus level using only hypervariable regions of 16S rDNA sequences, thus determining predominant microbiota. Heat maps of taxonomic composition were generated using the gplots library of R , frequencies were log2 transformed and clustered with Euclidean distance. In the case of samples before/after a meal, microbial diversity was established using the 16S and 23S rRNA gene. 16S and 23S sequences were binned using META-RNA 1.0 . 16S sequences were assigned using the online RDP assigner . 23S sequences were assigned using the SILVA database and SINA assigner [55, 56]. All statistical analyses were conducted on R v2.15. Non-parametric Spearman rank correlation was calculated among top 40 most frequent genera to associate frequency fluctuations during biofilm formation between genera. Then, Spearman's (ρ) coefficient and t-test significance was calculated for pairs of genera from all patients, using a Bonferroni correction for multiple comparisons.
Based on predominant microbiota present at all states of biofilm formation, available complete and WGS genomes were retrieved from the RefSeq and the Human Oral Microbiome databases [57, 58]. More than 80 genomes of oral related microorganisms were downloaded and used to build a local database with almost 300,000 coding sequences. More than 800 small predicted ORFs (100-400 nt) were removed, being 98-100% identical to different regions of 16S or 23S rDNAs . The remaining set of ORFs were then submitted to the KEGG Automatic Annotation Server  for KEGG Orthology (KO) assignment. Using MegaBlast v2.2.21 algorithm  with e-value cutoff 1e-08 and selecting alignments longer than 60% of read with >80% of identity, we assigned KO numbers and PATH categories to the BRITE functional hierarchy . Negative binomial distribution contained in DESeq  bioconductor v2.10 package (default parameters) was employed for differential expression analysis. KO over-representation was determined by comparison between early (6-12 h) and late (24-48 h) biofilm samples with q values ≤ 0.05. Counting of reads per gene and genome were normalized against genus frequency and size dataset and then transformed in log2 for comparison with qPCR expression data.
Primers for qPCR were designed submitting the respective ORF sequences from S. gordonii and A. naeslundii to the Primer3Plus webserver  (Additional file 6: Table S2). Gene amplification was performed using LightCycler® 480 System (Roche), SYBR Green I Master (Roche), and a small aliquot from the respective sample sequenced by Illumina. The Cp values were calculated from three replicates using the LightCycler® 480 SW software v1.5 (Roche). Expression was normalized against 16S rRNA expression from S. gordonii and A. naeslundii, respectively, and referred to expression seen for every gene at 6h for all patients in average using the ΔΔCt method.
All sequence data derived from 454 pyrosequencing of cDNA from samples after/before meal experiments, microbial diversity associated to dental quadrants, and Illumina HiSeq2000 sequencing of cDNA from oral biofilm are stored in the MG-RAST server to be publicly available by accessing to the “Oral Metatranscriptome” project, id 935 (http://metagenomics.anl.gov/linkin.cgi?project=935). Sequence data is also available at the European Nucleotide Archive (ENA-EBML) with provisional accession number ERP003984.
Authors thank nutritionist Mercedes Mora-Ruiz for advice in designing the carbohydrate-rich diet for the before/after meal experiment. Authors also thank Anny Camelo for assistance in qPCR procedures and Bernadent Dental Clinic (Valencia, Spain) for facilities during sample collection. This study was funded by projects MICROGEN CSD2009-00006 and 2012-40007 from Spanish MECO.
- 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.PubMedView ArticleGoogle Scholar
- Ahn J, Yang L, Paster BJ, Ganly I, Morris L, Pei Z, Hayes RB: Oral microbiome profiles: 16S rRNA pyrosequencing and microarray assay comparison. PLoS One. 2011, 6 (7): e22788-10.1371/journal.pone.0022788.PubMed CentralPubMedView ArticleGoogle Scholar
- Belda-Ferre P, Alcaraz LD, Cabrera-Rubio R, Romero H, Simon-Soro A, Pignatelli M, Mira A: The oral metagenome in health and disease. ISME J. 2012, 6 (1): 46-56. 10.1038/ismej.2011.85.PubMed CentralPubMedView ArticleGoogle 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Marsh PD: Dental plaque as a biofilm and a microbial community - implications for health and disease. BMC Oral Health. 2006, 6 (Suppl 1): S14-10.1186/1472-6831-6-S1-S14.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilson M: The oral cavity and its indigenous microbiota. Microbial inhabitants of humans. Edited by: Wilson M. 2005, New York: Cambridge University Press, 318-374.Google Scholar
- Jenkinson HF, Lamont RJ: Oral microbial communities in sickness and in health. Trends Microbiol. 2005, 13 (12): 589-595. 10.1016/j.tim.2005.09.006.PubMedView ArticleGoogle Scholar
- Palmer RJ, Gordon SM, Cisar JO, Kolenbrander PE: Coaggregation-mediated interactions of streptococci and actinomyces detected in initial human dental plaque. J Bacteriol. 2003, 185 (11): 3400-3409. 10.1128/JB.185.11.3400-3409.2003.PubMedView ArticleGoogle Scholar
- Kolenbrander PE, Andersen RN, Blehert DS, Egland PG, Foster JS, Palmer RJ: Communication among oral bacteria. Microbiol Mol Biol Rev. 2002, 66 (3): 486-505. 10.1128/MMBR.66.3.486-505.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Hojo K, Nagaoka S, Ohshima T, Maeda N: Bacterial interactions in dental biofilm development. J Dent Res. 2009, 88 (11): 982-990. 10.1177/0022034509346811.PubMedView ArticleGoogle Scholar
- Dige I, Raarup MK, Nyengaard JR, Kilian M, Nyvad B: Actinomyces naeslundii in initial dental biofilm formation. Microbiology. 2009, 155 (Pt 7): 2116-2126.PubMedView ArticleGoogle Scholar
- Kolenbrander PE, Andersen RN, Moore LV: Coaggregation of Fusobacterium nucleatum, Selenomonas flueggei, Selenomonas infelix, Selenomonas noxia, and Selenomonas sputigena with strains from 11 genera of oral bacteria. Infect Immun. 1989, 57 (10): 3194-3203.PubMed CentralPubMedGoogle Scholar
- Periasamy S, Kolenbrander PE: Aggregatibacter actinomycetemcomitans builds mutualistic biofilm communities with Fusobacterium nucleatum and Veillonella species in saliva. Infect Immun. 2009, 77 (9): 3542-3551. 10.1128/IAI.00345-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Jakubovics NS, Gill SR, Iobst SE, Vickerman MM, Kolenbrander PE: Regulation of gene expression in a mixed-genus community: stabilized arginine biosynthesis in Streptococcus gordonii by coaggregation with Actinomyces naeslundii. J Bacteriol. 2008, 190 (10): 3646-3657. 10.1128/JB.00088-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Frias-Lopez J, Duran-Pinedo A: Effect of periodontal pathogens on the metatranscriptome of a healthy multispecies biofilm model. J Bacteriol. 2012, 194 (8): 2082-2095. 10.1128/JB.06328-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, Nekrutenko A, Taylor J: Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol. 2010, 89: 19.10.1-19.10.21.Google Scholar
- Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, Zhang Y, Blankenberg D, Albert I, Taylor J, Miller W, Kent WJ, Nekrutenko A: Galaxy: a platform for interactive large-scale genome analysis. Genome Res. 2005, 15 (10): 1451-1455. 10.1101/gr.4086505.PubMed CentralPubMedView ArticleGoogle Scholar
- Goecks J, Nekrutenko A, Taylor J: Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010, 11 (8): R86-10.1186/gb-2010-11-8-r86.PubMed CentralPubMedView ArticleGoogle Scholar
- Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12 (6): R60-10.1186/gb-2011-12-6-r60.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Y, Hu T, Zhang J, Zhou X: Characterization of the Actinomyces naeslundii ureolysis and its role in bacterial aciduricity and capacity to modulate pH homeostasis. Microbiol Res. 2006, 161 (4): 304-310. 10.1016/j.micres.2005.11.002.PubMedView ArticleGoogle Scholar
- Takahashi N, Nyvad B: The role of bacteria in the caries process: ecological perspectives. J Dent Res. 2011, 90 (3): 294-303. 10.1177/0022034510379602.PubMedView ArticleGoogle Scholar
- Socransky SS, Haffajee AD, Cugini MA, Smith C, Kent RL: Microbial complexes in subgingival plaque. J Clin Periodontol. 1998, 25 (2): 134-144. 10.1111/j.1600-051X.1998.tb02419.x.PubMedView ArticleGoogle Scholar
- Chalmers NI, Palmer RJ, Kolenbrander PE, Cisar JO: Characterization of a Streptococcus sp.-Veillonella sp. community micromanipulated from dental plaque. J Bacteriol. 2008, 190 (24): 8145-8154. 10.1128/JB.00983-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Chung H, Pamp SJ, Hill JA, Surana NK, Edelman SM, Troy EB, Reading NC, Villablanca EJ, Wang S, Mora JR, Umesaki Y, Mathis D, Benoist C, Relman DA, Kasper DL: Gut immune maturation depends on colonization with a host-specific microbiota. Cell. 2012, 149 (7): 1578-1593. 10.1016/j.cell.2012.04.037.PubMed CentralPubMedView ArticleGoogle Scholar
- Dethlefsen L, McFall-Ngai M, Relman DA: An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature. 2007, 449 (7164): 811-818. 10.1038/nature06245.PubMedView ArticleGoogle Scholar
- Filoche S, Wong L, Sissons CH: Oral biofilms: emerging concepts in microbial ecology. J Dent Res. 2010, 89 (1): 8-18. 10.1177/0022034509351812.PubMedView ArticleGoogle Scholar
- Rajilic-Stojanovic M, Heilig HG, Tims S, Zoetendal EG, de Vos WM: Long-term monitoring of the human intestinal microbiota composition. Environ Microbiol. 2012, 15: 1146-1159.View ArticleGoogle Scholar
- Abiko Y, Sato T, Mayanagi G, Takahashi N: Profiling of subgingival plaque biofilm microflora from periodontally healthy subjects and from subjects with periodontitis using quantitative real-time PCR. J Periodontal Res. 2010, 45 (3): 389-395. 10.1111/j.1600-0765.2009.01250.x.PubMedView ArticleGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35 (Web Server issue): W182-W185.PubMed CentralPubMedView ArticleGoogle Scholar
- Björk GR, Hagervall TG: Transfer RNA modification. EcoSal—Escherichia coli and Salmonella: cellular and molecular biology. Edited by: Böck RCI JB, Neidhardt FC, Nyström T, Rudd KE, Squires CL. 2005, Washington, D.C: ASM PressGoogle Scholar
- Grosjean H: Fine tuning of RNA functions by modification and editing. Topics in Current Genetics. Edited by: Hohmann S. 2005, New York: Springer VerlagGoogle Scholar
- Shin JH, Lee HW, Kim SM, Kim J: Proteomic analysis of Acinetobacter baumannii in biofilm and planktonic growth mode. J Microbiol. 2009, 47 (6): 728-735. 10.1007/s12275-009-0158-y.PubMedView ArticleGoogle Scholar
- Kreth J, Merritt J, Shi W, Qi F: Co-ordinated bacteriocin production and competence development: a possible mechanism for taking up DNA from neighbouring species. Mol Microbiol. 2005, 57 (2): 392-404. 10.1111/j.1365-2958.2005.04695.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhattacharjee MK, Fine DH, Figurski DH: tfoX (sxy)-dependent transformation of Aggregatibacter (Actinobacillus) actinomycetemcomitans. Gene. 2007, 399 (1): 53-64. 10.1016/j.gene.2007.04.026.PubMed CentralPubMedView ArticleGoogle Scholar
- Pollack-Berti A, Wollenberg MS, Ruby EG: Natural transformation of Vibrio fischeri requires tfoX and tfoY. Environ Microbiol. 2010, 12 (8): 2302-2311.PubMed CentralPubMedGoogle Scholar
- Human Microbiome Project Consortium: A framework for human microbiome research. Nature. 2012, 486 (7402): 215-221. 10.1038/nature11209.View ArticleGoogle Scholar
- McNab R, Ford SK, El-Sabaeny A, Barbieri B, Cook GS, Lamont RJ: LuxS-based signaling in Streptococcus gordonii: autoinducer 2 controls carbohydrate metabolism and biofilm formation with Porphyromonas gingivalis. J Bacteriol. 2003, 185 (1): 274-284. 10.1128/JB.185.1.274-284.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Nobbs AH, Vajna RM, Johnson JR, Zhang Y, Erlandsen SL, Oli MW, Kreth J, Brady LJ, Herzberg MC: Consequences of a sortase A mutation in Streptococcus gordonii. Microbiology. 2007, 153 (Pt 12): 4088-4097.PubMedView ArticleGoogle Scholar
- Ton-That H, Marraffini LA, Schneewind O: Protein sorting to the cell wall envelope of Gram-positive bacteria. Biochim Biophys Acta. 2004, 1694 (1–3): 269-278.PubMedView ArticleGoogle Scholar
- Ton-That H, Mazmanian SK, Faull KF, Schneewind O: Anchoring of surface proteins to the cell wall of Staphylococcus aureus. Sortase catalyzed in vitro transpeptidation reaction using LPXTG peptide and NH(2)-Gly(3) substrates. J Biol Chem. 2000, 275 (13): 9876-9881. 10.1074/jbc.275.13.9876.PubMedView ArticleGoogle Scholar
- Mishra A, Wu C, Yang J, Cisar JO, Das A, Ton-That H: The Actinomyces oris type 2 fimbrial shaft FimA mediates co-aggregation with oral streptococci, adherence to red blood cells and biofilm development. Mol Microbiol. 2010, 77: 841-854.PubMed CentralPubMedGoogle Scholar
- Jakubovics NS, Kerrigan SW, Nobbs AH, Stromberg N, van Dolleweerd CJ, Cox DM, Kelly CG, Jenkinson HF: Functions of cell surface-anchored antigen I/II family and Hsa polypeptides in interactions of Streptococcus gordonii with host receptors. Infect Immun. 2005, 73 (10): 6629-6638. 10.1128/IAI.73.10.6629-6638.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Jakubovics NS, Stromberg N, van Dolleweerd CJ, Kelly CG, Jenkinson HF: Differential binding specificities of oral streptococcal antigen I/II family adhesins for human or bacterial ligands. Mol Microbiol. 2005, 55 (5): 1591-1605. 10.1111/j.1365-2958.2005.04495.x.PubMedView ArticleGoogle Scholar
- Kolenbrander PE, Palmer RJ, Periasamy S, Jakubovics NS: Oral multispecies biofilm development and the key role of cell-cell distance. Nat Rev Microbiol. 2010, 8 (7): 471-480. 10.1038/nrmicro2381.PubMedView ArticleGoogle Scholar
- McNab R, Holmes AR, Clarke JM, Tannock GW, Jenkinson HF: Cell surface polypeptide CshA mediates binding of Streptococcus gordonii to other oral bacteria and to immobilized fibronectin. Infect Immun. 1996, 64 (10): 4204-4210.PubMed CentralPubMedGoogle Scholar
- McNab R, Jenkinson HF, Loach DM, Tannock GW: Cell-surface-associated polypeptides CshA and CshB of high molecular mass are colonization determinants in the oral bacterium Streptococcus gordonii. Mol Microbiol. 1994, 14 (4): 743-754. 10.1111/j.1365-2958.1994.tb01311.x.PubMedView ArticleGoogle Scholar
- Stratul S, Didilescu A, Hanganu C, Greabu M, Totan A, Spinu T, Onisei D, Rusu D, Jentsch H, Sculean A: On the molecular basis of biofilm formation. Oral biofilms and systemic infections. TMJ. 2008, 58: 118-123.Google Scholar
- Loozen G, Ozcelik O, Boon N, De Mol A, Schoen C, Quirynen M, Teughels W: Inter-bacterial correlations in subgingival biofilms: a large-scale survey. J Clin Periodontol. 2014, 41 (1): 1-10. 10.1111/jcpe.12167.PubMedView ArticleGoogle Scholar
- Ammann TW, Belibasakis GN, Thurnheer T: Impact of early colonizers on in vitro subgingival biofilm formation. PLoS One. 2013, 8 (12): e83090-10.1371/journal.pone.0083090.PubMed CentralPubMedView ArticleGoogle Scholar
- Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37 (Database issue): D141-D145.PubMed CentralPubMedView ArticleGoogle Scholar
- Hartmann M, Howes CG, Abarenkov K, Mohn WW, Nilsson RH: V-Xtractor: an open-source, high-throughput software tool to identify and extract hypervariable regions of small subunit (16S/18S) ribosomal RNA gene sequences. J Microbiol Methods. 2010, 83 (2): 250-253. 10.1016/j.mimet.2010.08.008.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1016/S0022-2836(05)80360-2.PubMedView ArticleGoogle Scholar
- Warnes G, Bolker B, Bonebakker L, Gentleman R, Liaw WHA, Lumley T, Maechler M, Magnusson A, Moeller S, Schwartz M, Venables B: gplots: Various R programming tools for plotting data. The Comprehensive R Archive Network. 2009Google Scholar
- Huang Y, Gilna P, Li W: Identification of ribosomal RNA genes in metagenomic fragments. Bioinformatics. 2009, 25 (10): 1338-1340. 10.1093/bioinformatics/btp161.PubMed CentralPubMedView ArticleGoogle Scholar
- Pruesse E, Peplies J, Glockner FO: SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012, 28 (14): 1823-1829. 10.1093/bioinformatics/bts252.PubMed CentralPubMedView ArticleGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35 (21): 7188-7196. 10.1093/nar/gkm864.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen T, Yu WH, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE: The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database (Oxford). 2010, 2010: baq013-View ArticleGoogle Scholar
- Pruitt KD, Tatusova T, Brown GR, Maglott DR: NCBI Reference Sequences (RefSeq): current status, new features and genome annotation policy. Nucleic Acids Res. 2012, 40 (Database issue): D130-D135.PubMed CentralPubMedView ArticleGoogle Scholar
- Tripp HJ, Hewson I, Boyarsky S, Stuart JM, Zehr JP: Misannotations of rRNA can now generate 90% false positive protein matches in metatranscriptomic studies. Nucleic Acids Res. 2011, 39 (20): 8792-8802. 10.1093/nar/gkr576.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JA: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35 (Web Server issue): W71-W74.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.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.