BACKGROUND
dental caries and the periodontal diseases are the commonest bacterial diseases of man and can result in the loss of the teeth and their associated function, and are also significant risk factors for disease at other body sites, particularly cardiovascular disease  <cit> .
treatment of oral diseases is expensive and efforts have therefore focused on prevention, particularly through the use of mouthrinses and toothpastes containing active agents that control the bacteria found in dental plaque.
the evaluation of novel active agents and formulations in humans is time-consuming and difficult and toxicological data may not be available for new compounds.
consequently, attempts have been made to develop in-vitro assays that accurately predict efficacy in vivo  <cit> .
these assays enable researchers to perform preliminary screening of active agents in order to identify candidates for subsequent testing in clinical trials, wherein the therapeutic effects, as well as issues such as substantivity in the oral cavity  <cit> , can be determined.
many existing in-vitro assays cultivate oral bacteria for testing in planktonic suspension but it has been shown that bacteria naturally form biofilms  <cit> .
bacteria in biofilms exhibit greater resistance to antimicrobial agents than cells growing in planktonic culture  <cit> .
in addition, the composition of oral biofilms is typically highly complex and includes a substantial number of species which have yet to be cultivated  <cit> .
recent deep sequencing studies have, for example, detected hundreds of species in dental plaque samples from individual subjects  <cit> .
although some in-vitro oral biofilm assays have been developed, these have typically used relatively simple defined inocula  <cit> .
given the high richness and diversity of oral biofilms, it would be preferable to use natural inocula in order to more accurately represent the in-vivo ecosystem.
one in-vitro system that has been previously developed and used to grow bacteria as biofilms is the calgary biofilm device   <cit> .
in this system, biofilms are grown on pegs protruding from the lid of a 96-well plate.
the pegs are immersed in a growth medium that can easily be replaced by transferring the lid to a new baseplate, thereby enabling the long-term growth of biofilms.
the cbd was originally developed to determine the susceptibility of bacterial biofilms to antibiotics  <cit>  for applications such as medical device-related infections, and is therefore commercially available as the ‘minimum biofilm eradication concentration’  assay .
previous work has demonstrated that uniform biofilms with reproducible total viable counts can be obtained when using simple defined bacterial inocula  <cit> .
the cbd has also been used to examine the interactions among five common oral species when grown anaerobically for up to 36 hours  <cit> .
using quantitative-pcr, the authors showed that porphyromonas gingivalis cell counts increased when grown together with a veillonella sp., fusobacterium nucleatum, or aggregatibacter actinomycetemcomitans, suggesting mutualism between these species.
the overall aim of this study was to develop an oral microbial biofilm model derived from a natural inoculum: the saliva of healthy individuals, using the cbd.
the biofilms thus generated will be suitable for assessing the impact of different oral care product components on oral biofilm composition.
the specific aims of this study were to determine the reproducibility, stability and variability of the model by using pyrosequencing of partial 16s rrna genes to compare the composition of the biofilms over time derived from saliva from  the same individuals at different time points and  from different panels of volunteers.
methods
participants
ethical approval for this study was granted by queen mary university of london ethics of research committee .
informed consent was obtained from all of the individuals who participated.
all of the participants were between  <dig> and 65 years of age and were medically healthy volunteers who were staff or postgraduate students at queen mary university of london.
any subjects with systemic conditions that may have affected their immune or inflammatory status were excluded from the study.
a total of  <dig> subjects participated in the study.
sample collection
un-stimulated saliva samples were obtained from the participants by expectoration into sterile universal tubes.
saliva was collected between approximately 14: <dig> and 15: <dig> on the days that the biofilms were to be inoculated.
participants were grouped into panels of six and their saliva samples pooled together in equal volumes: 1 ml was used from each individual to produce a 6 ml pooled sample.
pooled saliva was placed on ice and processed within an hour.
one panel was sampled at three different time points, a week apart, and two panels were sampled at one time point.
inoculation of the calgary biofilm device
the saliva was vortexed for 15 s and 200 μl was pipetted per well of a 96-well microplate, up to the required number of wells.
wells around the outside of the microplate were not used.
the lid of the cbd was fitted onto the microplate so that the hydroxyapatite-coated pegs were bathed in the saliva.
the cbd plate was then incubated at 37°c under anaerobic conditions  for 18 hours, after which the lid was transferred to a new baseplate containing 200 μl of pre-reduced brain heart infusion broth  growth medium supplemented with hog gastric mucin , haemin , and vitamin k .
the growth medium was changed after every  <dig>  days of anaerobic incubation.
biofilms were harvested from half of the pegs after seven days and from the remaining half after 14 days.
removal of pegs and propidium monoazide treatment of samples
pegs with biofilms were snapped off the lid with sterile pliers and washed by dipping into sterile phosphate-buffered saline  three times.
all of the visible biofilm material was then removed using a sterile curette and suspended into 500 μl of sterile pbs.
the material from three pegs was pooled to produce one sample for analysis, and three samples were processed for each incubation time.
each sample, and half of the original saliva sample, was subjected to propidium monoazide  treatment to prevent subsequent pcr amplification of extracellular dna and dna from dead or damaged cells  <cit> :  <dig>  μl of pma was added  to the cells suspended in pbs and incubated in the dark with occasional shaking for  <dig> mins at room temperature.
the samples were then exposed to light from a 500 w halogen lamp for  <dig> mins at a distance of 20 cm in order to form a covalent linkage between the pma and the dna.
during the exposure time the samples were placed on ice to avoid excessive heating and subjected to occasional shaking.
the samples were used for dna extractions immediately after the pma treatment.
dna extraction
dna was extracted from the saliva and biofilm samples using the genelute bacterial dna extraction kit .
extractions were performed on both pma-treated and -untreated aliquots  of pooled saliva samples.
dna extraction was performed following the manufacturer’s instructions with an additional cell lysis step to increase the recovery of dna from gram-positive cells, in which samples were incubated in a 45 mg/ml lysozyme solution at 37°c for  <dig> mins.
molecular microbiological analysis
the bacterial composition of the biofilms and saliva was determined using  <dig> pyrosequencing of partial 16s rrna genes as described previously  <cit> , with some minor modifications.
pcr amplification of a fragment of the 16s rrna gene, approximately 500 bp in length covering the v1-v <dig> hypervariable regions, was performed for each dna sample using composite fusion primers.
the fusion primers comprised the broad-range 16s rrna gene primers  <dig> fym  <cit>  and  <dig> r  <cit>  along with roche gs-flx titanium series adapter sequences  for 454-pyrosequencing using the lib-l emulsion-pcr method.
the forward primers included previously described 12-base error-correcting golay barcodes.
pcr reactions were performed using extensor hi-fidelity pcr mastermix  along with the appropriate barcoded forward primer and the reverse primer.
the pcr conditions were as follows:  <dig> mins initial denaturation at 95°c, followed by 25 cycles of 95°c for 45 s, 53°c for 45 s and 72°c for 45 s and a final extension of 72°c for  <dig> mins.
pcr amplicons were then purified using the qiaquick pcr purification kit  according to the manufacturer’s instructions.
the size and purity of the amplicons was checked using the agilent dna  <dig> kit and the agilent  <dig> bioanalyzer.
quantitation of the amplicons was performed by means of a fluorometric assay using the quant-it picogreen fluorescent nucleic acid stain .
the amplicons were then pooled together at equimolar concentrations .
emulsion-pcr and unidirectional sequencing of the samples was performed using the lib-l kit and the roche  <dig> gs-flx + titanium series sequencer by the department of biochemistry, cambridge university, cambridge, uk.
the raw sequence data were deposited with the ncbi sra database as accession srp <dig>
sequence analysis
sequence analysis was performed using the ‘mothur’ software suite version  <dig>   <cit> , following the  <dig> standard operating procedure  <cit>  on mothur.org.
the sequences were denoised using the ampliconnoise algorithm  <cit> , as implemented by mothur.
sequences that were less than  <dig> bases in length and/or had one of the following: > <dig> mismatches to the primer, > <dig> mismatch to the barcode regions, and homopolymers of > <dig> bases in length, were discarded.
the remaining sequences were trimmed to remove primers and barcodes and aligned to the silva 16s rrna reference alignment  <cit> .
the uchime algorithm  <cit>  was used to identify chimeric sequences, which were then removed from the dataset.
sequences were clustered into operational taxonomic units  at a genetic distance of  <dig>  using the average neighbour algorithm and identified using a naïve bayesian classifier  <cit>  with the human oral microbiome database  reference set .
for those otus that could not be identified using the bayesian classifier, representative sequences were obtained in mothur  and identified using blast against the homd reference set .
the possible alternatives for the species identification were then provided.
analysis of alpha and beta diversity
the sequences for each sample were randomly sub-sampled to the same number  for the alpha and beta diversity analyses.
the extent of sampling of the communities was assessed using good’s non-parametric coverage estimator  <cit> .
the diversity of the communities was calculated using simpson’s inverse diversity index  <cit> .
the beta-diversity of the samples was analysed using distance matrices generated using the jaccard index and the thetayc calculator  <cit> .
the distance matrices were visualised using principal coordinates analysis  plots generated in r .
statistical analysis
analysis of molecular variance   <cit> , as implemented in mothur, was used to determine if there were statistically significant differences in jaccard index and thetayc distances between saliva and biofilm samples and between biofilms from different time points, subject panels, and incubation times.
the mean relative abundances of phyla, genera, and species-level otus in biofilms was determined as follows: the proportions of sequences assigned to a particular taxon were calculated for biofilms derived from different panels at a single time point only  after 7- and 14-day incubations.
the pma-treated saliva samples were used for statistical comparisons of the taxonomic composition of saliva to biofilms.
to determine if there were significant differences in otu richness and diversity between the incubation times, paired t-tests were performed in r.
RESULTS
pyrosequencing
a total of  <dig>  sequences with a mean length of  <dig> bases were obtained for analysis after quality filtering, screening of the sequence alignment, and removal of chimeras.
a mean of  <dig> sequences  were obtained per sample.
one replicate biofilm sample derived from panel  <dig> after seven days of incubation  was not included in the sequencing run due to poor pcr amplification.
the number of otus  detected in individual biofilm and pooled saliva samples ranged between  <dig> and  <dig>  the mean number of otus detected was  <dig> in the biofilms and  <dig> in the saliva samples.
a table summarising the alpha diversity of the biofilms and saliva is shown in additional file  <dig>  the number of otus detected in the seven-day biofilms  was significantly higher  than in the 14-day biofilms .
however, there was no significant difference in diversity  between the seven- and 14-day biofilms.
reproducibility of biofilms and shifts in biofilm otu composition over time
comparison of the community membership and structure of biofilms using principal coordinates analysis  plots indicated that replicate biofilm samples, derived from the same saliva pool after the same incubation time but harvested from different pegs, were highly similar and clustered together .
in addition, biofilms derived from the same panel  at different times were similar  and amova tests found there to be no statistically significant difference between the time points.
however, there were significant differences in both the membership and structure of biofilms derived from the three different subject panels.
the most significant differences by amova were between biofilms from panel  <dig> and panel  <dig> .
interestingly, pcoa indicated that the dissimilarity in community structure between panels was greater after 14 days than after seven days .figure  <dig>
principal coordinates analysis of biofilms derived from different panels after different incubation times.
plots are based on community membership using the jaccard index  and community structure using the thetayc calculator .
blue - panel 1; red - panel 2; black - panel  <dig>  labels indicate the incubation time.
a: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
b: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
principal coordinates analysis of biofilm replicates from panel  <dig> at different time points.
plots are based on community membership using the jaccard index  and community structure using the thetayc calculator .
blue - time 1; red - time 2; black - time  <dig>  a: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
b: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
there was a directional shift along the axes in the pcoa plots between seven and 14 days of incubation for biofilms derived from all three panels .
amova confirmed that there was a significant difference between the 7-day and 14-day incubations both in terms of membership and structure .
analysis using lefse identified a total of  <dig> otus that were significantly differentially abundant between incubation times.
the identities of the otus with lda effect size scores of > <dig>  are shown in figure  <dig>  an otu identified as veillonella parvula was most strongly associated with the 7-day incubations, whilst parvimonas micra was most strongly associated with the 14-day incubations.figure  <dig>
linear discriminant analysis effect size  analysis showing those otus that were significantly differentially abundant between seven- and 14-day incubation times, ranked by effect size .
hierarchical cluster analysis showed that the biofilm samples clustered by panel in a dendrogram based on community membership  and predominantly by incubation time in a dendrogram based on structure .
otu-based comparisons of biofilms to saliva
comparison of saliva samples with the biofilms using pcoa plots revealed a differing community membership and structure .
amova tests confirmed that there were significant differences in both the membership and structure of saliva compared to seven and 14-day biofilms .
analysis using lefse identified  <dig> otus that were significantly differentially abundant between saliva and biofilms that had been grown for seven days.
a list of differentially abundant otus with lda effect size scores of > <dig>  is shown in figure  <dig>  an otu identified as neisseria flavescens/subflava was most strongly associated with saliva, whilst veillonella parvula was most strongly associated with the biofilms.figure  <dig>
principal coordinates analysis of saliva and seven-day biofilms.
plots are based on community membership using the jaccard index  and community structure using the thetayc calculator .
blue - biofilms; red - saliva.
labels indicate the panel number.
a: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
b: pc1 =  <dig> % of variance, pc2 =  <dig> % of variance.
linear discriminant analysis effect size  analysis showing those otus that were significantly differentially abundant between saliva and seven-day biofilms, ranked by effect size .
comparison of pma-treated and untreated saliva samples
comparisons using pcoa showed that pma-treated and untreated saliva were similar and samples clustered by the panel from which they were obtained, rather than by treatment .
there were no significant differences in community membership or structure between pma-treated and untreated saliva using amova tests.
taxonomic composition of the biofilms
the predominant phyla detected in all of the biofilms in order of mean relative abundance were: firmicutes, bacteroidetes, synergistetes, fusobacteria, proteobacteria and actinobacteria.
other phyla that were detected in minor relative proportions, and not in every sample, included: sr <dig>  spirochaetes, tm <dig> and tenericutes.
a total of  <dig> genera were detected in the biofilms, the most abundant of which were: veillonella, streptococcus and prevotella in seven-day biofilms, and streptococcus, prevotella and parvimonas in 14-day biofilms.
figure  <dig> shows the relative abundances of the predominant genera detected in the biofilms after seven and 14 days incubation.
the relative abundances of the predominant genera detected in the saliva samples, from which biofilms were derived, are shown in additional file  <dig>  a table detailing all of the taxa identified down to the species level, and their relative proportions in individual biofilm and saliva samples, can be found in additional file  <dig> figure  <dig>
predominant genera detected in the biofilms.
the graph shows the mean relative abundances of genera that were detected in seven and 14-day biofilms derived from three panels.
genera shown are those with mean relative abundances of > 1%.
error bars show the standard error of the mean .
discussion
this study has demonstrated that complex oral biofilms derived from a natural saliva inoculum can be successfully grown and maintained using the cbd.
the results showed that the biofilms had a richness and diversity close to that of the pooled saliva inocula, with a mean of  <dig> species-level otus detected per biofilm sample compared to a mean of  <dig> in saliva.
there was a significant difference in community membership and structure between the saliva and the biofilms, with some otus detected in saliva not present in the biofilms.
this is not surprising because the bacterial composition of saliva is known to differ to that of dental plaque  <cit> .
saliva, however, does include representatives of the various surfaces found in the mouth and is thus a useful inoculum for biofilms.
hydroxyapatite-coated pegs were used for the purpose of mimicking the teeth in order to obtain biofilms with a similar composition to plaque.
another possible reason for differences between the inocula and the biofilms is that specific nutrients or growth factors required by certain species could have been absent in the growth medium used.
a brain-heart infusion  based medium was chosen in this study because it has been successfully used to cultivate a broad range of fastidious and non-fastidious oral bacteria  <cit> .
the bhi was supplemented with mucin, vitamin k, and haemin, as some oral species grow poorly or not at all in the absence of one or more of these substances.
for example, a number of black-pigmented species of prevotella and porphyromonas require haemin and vitamin k for growth  <cit> .
in addition, hog gastric mucin, a high molecular weight glycoprotein, has been shown to support the growth of mixed communities of oral bacteria when used as the principal source of carbon and energy  <cit> .
future work could investigate the use of different media, such as an artificial saliva-based medium, with this model.
another reason that certain salivary species may have been lost is that the anaerobic atmosphere in which the cbd was incubated would have selected against the growth of aerobic species.
the absence of host immune cells and molecules in the second phase of growth may also have had an impact on the community composition.
nevertheless, a highly diverse community of oral bacteria was maintained which included the genera known to be predominant in plaque and also a variety of fastidious and uncultivated taxa, such as un-named bacteroidetes, lachnospiraceae, clostridiales and peptostreptococcaceae species.
both saliva and biofilm samples were treated with propidium monoazide prior to dna extraction to avoid detection of bacterial cells that were non-viable; this method has been shown to prevent pcr amplification of dna from dead or damaged cells  <cit> .
interestingly, there was no significant difference between pma-treated and untreated saliva in terms of community membership and structure.
this suggests that the vast majority of taxa detected in the saliva samples were viable.
this could be due to the rapid processing of the samples performed in order to avoid loss of cell viability.
in addition, human saliva has been shown to contain dnase i produced by the parotid glands  <cit> , which could rapidly break down extracellular bacterial dna from dead cells.
the ability to grow biofilms with a richness and diversity close to that detected in oral habitats, and that include previously uncultivated taxa, is a major strength of this model over others.
for instance, the well-established zürich biofilm model  <cit>  has been developed with defined inocula consisting of five, or more recently,  <dig> cultivable species  <cit> .
defined biofilms consisting of a low number of selected cultivable species are less representative of the in-vivo ecosystem and may, therefore, be less accurate in predicting the efficacy of an active agent.
another study recently reported the development of an in-vitro biofilm model in which a natural saliva inoculum was used and the composition of the samples determined by pyrosequencing of 16s rrna genes  <cit> .
whilst the authors also reported a high microbial diversity, the incubation times used were relatively short  and this may explain why certain slow-growing oral taxa, including fretibacterium spp.
and tannerella spp.
were not detected, whilst streptococci were dominant with s. vestibularis constituting approximately 40% of the communities.
the taxonomic composition of the biofilms grown in this study was similar to that of dental plaque.
streptococcus, veillonella and prevotella were the predominant genera detected in the biofilms, all of which have been shown to be major constituents of plaque  <cit> .
the otu detected with the highest mean relative abundance in the biofilms was veillonella parvula, which was also detected with the highest rank abundance in an extensive cloning and sanger sequencing study of the human oral microbiome  <cit> .
in addition, periodontitis-associated species, including the ‘red complex’: porphyromonas gingivalis, treponema denticola, tannerella forsythia, and the gram-positive anaerobes filifactor alocis and parvimonas micra, were all detected in the biofilms.
these organisms have been strongly associated with deep periodontal pockets in individuals with severe chronic periodontitis  <cit> .
in the cbd model, p. micra was the organism most strongly associated with 14-day biofilms.
species that had a significantly lower relative abundance in 14-day biofilms than 7-day biofilms included the streptococcal species s. cristatus and s. salivarius / vestibularis, which have previously been associated with health  <cit> .
species among the genera neisseria and rothia were detected at only very low proportions in the biofilms, despite being abundant in saliva.
this is likely explained by the anaerobic incubation of the cbd, as these organisms grow optimally under aerobic conditions  <cit> .
it has been shown that the redox potential  of dental plaque rapidly falls as the biofilm develops in vivo  <cit> .
in addition, experimental gingivitis studies have shown that plaque accumulating in the absence of oral hygiene supports the growth of increasing numbers of anaerobic species, many of which are gingivitis-associated  <cit> .
this study aimed to grow biofilms that were similar in composition to biofilms that would develop naturally in vivo without oral hygiene intervention, and anaerobic incubation was chosen in order to reproducibly obtain a biofilm typical of mature plaque.
however, future work could examine the composition of cbd oral biofilms grown under aerobic conditions and compare them to those grown anaerobically.
if using the biofilm model to screen antimicrobial agents or oral care product components, it would be useful to grow the biofilms under both aerobic and anaerobic conditions in order to determine the effect of a given substance on as diverse a range of oral taxa as possible.
future studies could also compare the similarity of the in-vitro biofilms to dental plaque biofilms that form naturally in vivo in the same individuals abstaining from oral hygiene.
this would further confirm that this model generates oral biofilms that are representative of those formed in vivo.
the results of this study showed that the bacterial composition of the biofilms was highly reproducible for sample replicates from different pegs derived from the same saliva pool and incubated for the same length of time.
moreover, the biofilms were similar in composition when derived from the same panel at different time points.
after 14 days of incubation, when dissimilarity in the biofilms might have been expected to increase, the biofilms from the same panel clustered closely in the pcoa plots.
the differences between biofilms derived from different panels was not surprising given the high inter-individual variation in bacterial diversity found in the normal human oral microbiome  <cit> , although, an attempt was made to reduce this variability by pooling saliva from six individuals for use as the inoculum.
hierarchical clustering of the biofilm samples in dendrograms indicated that the panel was the primary determinant of community membership, but that incubation time had a stronger influence on community structure.
this is likely to be because the relative abundances of otus would be expected to change over time as the biofilms mature.
due to the differences in both membership and structure of biofilms derived from different panels, the same panel should be used to provide the saliva inoculum in future studies that aim to compare biofilms grown under different conditions, or after exposure to different challenges e.g.
antimicrobial agents.
the reproducible biofilms that can be obtained using this model will enable relatively small changes in bacterial composition to be detected.
this will be particularly useful for assessing the impact of oral care products that aim to manipulate or alter, rather than eradicate, plaque.
this includes active agents or bacteriocin-producing probiotics that target particular taxa, or prebiotics that could promote the growth of health-associated bacteria.
in the case of probiotics, the model could also be useful in helping to predict whether or not a particular strain is likely to colonise and persist within oral biofilms.
in addition to determining changes in the community composition of the biofilms, future work could also investigate changes in community function using metagenomics, metabolomics and metatranscriptomics, in response to different active agents or changes in key environmental parameters.
CONCLUSIONS
this study has successfully developed an oral biofilm model using the cbd seeded with a natural saliva inoculum.
the biofilms generated were highly complex and comprised of microbial taxa that are commonly found in dental plaque.
in addition, their composition was shown to be reproducible when derived from the pooled saliva of the same panel of individuals.
this model will therefore be useful for screening novel antimicrobial agents and also pre- or probiotics that aim to modify plaque composition.
additional files
additional file 1:
table showing the alpha diversity of pooled saliva and biofilm samples.
all samples were sub-sampled to  <dig> sequences.
p - panel; t - time point; s – saliva; sp – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate.
additional file 2:
dendrogram showing the similarity of biofilm and saliva samples based on community membership .
p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate.
additional file 3:
dendrogram showing the similarity of biofilm and saliva samples based on community structure .
p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate.
additional file 4:
bar chart of the predominant genera detected in pma-treated saliva samples.
the chart shows the mean relative abundances of genera that were detected in all three of the pooled saliva samples from different panels.
genera shown are those with mean relative abundances of > 1%.
error bars show the standard error of the mean .
additional file 5:
table showing the classification of sequences to the species level in each sample.
the table shows the relative abundances of the phylotypes in the different samples.
p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
jok contributed to the conception and design of the study, analysis of the data, and preparation of the manuscript.
mp contributed to the conception and design of the study and critically revised the manuscript.
wgw contributed to the conception and design of the study, analysis of the data, and preparation of the manuscript.
all authors read and approved the final manuscript.

 dental caries and the periodontal diseases are the commonest bacterial diseases of man and can result in the loss of the teeth and their associated function, and are also significant risk factors for disease at other body sites, particularly cardiovascular disease. treatment of oral diseases is expensive and efforts have therefore focused on prevention, particularly through the use of mouthrinses and toothpastes containing active agents that control the bacteria found in dental plaque. the evaluation of novel active agents and formulations in humans is time-consuming and difficult and toxicological data may not be available for new compounds. consequently, attempts have been made to develop in-vitro assays that accurately predict efficacy in vivo. these assays enable researchers to perform preliminary screening of active agents in order to identify candidates for subsequent testing in clinical trials, wherein the therapeutic effects, as well as issues such as substantivity in the oral cavity, can be determined. many existing in-vitro assays cultivate oral bacteria for testing in planktonic suspension but it has been shown that bacteria naturally form biofilms. bacteria in biofilms exhibit greater resistance to antimicrobial agents than cells growing in planktonic culture. in addition, the composition of oral biofilms is typically highly complex and includes a substantial number of species which have yet to be cultivated. recent deep sequencing studies have, for example, detected hundreds of species in dental plaque samples from individual subjects. although some in-vitro oral biofilm assays have been developed, these have typically used relatively simple defined inocula. given the high richness and diversity of oral biofilms, it would be preferable to use natural inocula in order to more accurately represent the in-vivo ecosystem. one in-vitro system that has been previously developed and used to grow bacteria as biofilms is the calgary biofilm device. in this system, biofilms are grown on pegs protruding from the lid of a 96-well plate. the pegs are immersed in a growth medium that can easily be replaced by transferring the lid to a new baseplate, thereby enabling the long-term growth of biofilms. the cbd was originally developed to determine the susceptibility of bacterial biofilms to antibiotics for applications such as medical device-related infections, and is therefore commercially available as the ‘minimum biofilm eradication concentration’ assay. previous work has demonstrated that uniform biofilms with reproducible total viable counts can be obtained when using simple defined bacterial inocula. the cbd has also been used to examine the interactions among five common oral species when grown anaerobically for up to 36 hours. using quantitative-pcr, the authors showed that porphyromonas gingivalis cell counts increased when grown together with a veillonella sp., fusobacterium nucleatum, or aggregatibacter actinomycetemcomitans, suggesting mutualism between these species. the overall aim of this study was to develop an oral microbial biofilm model derived from a natural inoculum: the saliva of healthy individuals, using the cbd. the biofilms thus generated will be suitable for assessing the impact of different oral care product components on oral biofilm composition. the specific aims of this study were to determine the reproducibility, stability and variability of the model by using pyrosequencing of partial 16s rrna genes to compare the composition of the biofilms over time derived from saliva from the same individuals at different time points and from different panels of volunteers. methods participants ethical approval for this study was granted by queen mary university of london ethics of research committee. informed consent was obtained from all of the individuals who participated. all of the participants were between and 65 years of age and were medically healthy volunteers who were staff or postgraduate students at queen mary university of london. any subjects with systemic conditions that may have affected their immune or inflammatory status were excluded from the study. a total of subjects participated in the study. sample collection un-stimulated saliva samples were obtained from the participants by expectoration into sterile universal tubes. saliva was collected between approximately 14: and 15: on the days that the biofilms were to be inoculated. participants were grouped into panels of six and their saliva samples pooled together in equal volumes: 1 ml was used from each individual to produce a 6 ml pooled sample. pooled saliva was placed on ice and processed within an hour. one panel was sampled at three different time points, a week apart, and two panels were sampled at one time point. inoculation of the calgary biofilm device the saliva was vortexed for 15 s and 200 μl was pipetted per well of a 96-well microplate, up to the required number of wells. wells around the outside of the microplate were not used. the lid of the cbd was fitted onto the microplate so that the hydroxyapatite-coated pegs were bathed in the saliva. the cbd plate was then incubated at 37°c under anaerobic conditions for 18 hours, after which the lid was transferred to a new baseplate containing 200 μl of pre-reduced brain heart infusion broth growth medium supplemented with hog gastric mucin, haemin, and vitamin k. the growth medium was changed after every  days of anaerobic incubation. biofilms were harvested from half of the pegs after seven days and from the remaining half after 14 days. removal of pegs and propidium monoazide treatment of samples pegs with biofilms were snapped off the lid with sterile pliers and washed by dipping into sterile phosphate-buffered saline three times. all of the visible biofilm material was then removed using a sterile curette and suspended into 500 μl of sterile pbs. the material from three pegs was pooled to produce one sample for analysis, and three samples were processed for each incubation time. each sample, and half of the original saliva sample, was subjected to propidium monoazide treatment to prevent subsequent pcr amplification of extracellular dna and dna from dead or damaged cells:  μl of pma was added to the cells suspended in pbs and incubated in the dark with occasional shaking for mins at room temperature. the samples were then exposed to light from a 500 w halogen lamp for mins at a distance of 20 cm in order to form a covalent linkage between the pma and the dna. during the exposure time the samples were placed on ice to avoid excessive heating and subjected to occasional shaking. the samples were used for dna extractions immediately after the pma treatment. dna extraction dna was extracted from the saliva and biofilm samples using the genelute bacterial dna extraction kit. extractions were performed on both pma-treated and -untreated aliquots of pooled saliva samples. dna extraction was performed following the manufacturer’s instructions with an additional cell lysis step to increase the recovery of dna from gram-positive cells, in which samples were incubated in a 45 mg/ml lysozyme solution at 37°c for mins. molecular microbiological analysis the bacterial composition of the biofilms and saliva was determined using pyrosequencing of partial 16s rrna genes as described previously, with some minor modifications. pcr amplification of a fragment of the 16s rrna gene, approximately 500 bp in length covering the v1-v hypervariable regions, was performed for each dna sample using composite fusion primers. the fusion primers comprised the broad-range 16s rrna gene primers fym and r along with roche gs-flx titanium series adapter sequences for 454-pyrosequencing using the lib-l emulsion-pcr method. the forward primers included previously described 12-base error-correcting golay barcodes. pcr reactions were performed using extensor hi-fidelity pcr mastermix along with the appropriate barcoded forward primer and the reverse primer. the pcr conditions were as follows: mins initial denaturation at 95°c, followed by 25 cycles of 95°c for 45 s, 53°c for 45 s and 72°c for 45 s and a final extension of 72°c for mins. pcr amplicons were then purified using the qiaquick pcr purification kit according to the manufacturer’s instructions. the size and purity of the amplicons was checked using the agilent dna kit and the agilent bioanalyzer. quantitation of the amplicons was performed by means of a fluorometric assay using the quant-it picogreen fluorescent nucleic acid stain. the amplicons were then pooled together at equimolar concentrations. emulsion-pcr and unidirectional sequencing of the samples was performed using the lib-l kit and the roche gs-flx + titanium series sequencer by the department of biochemistry, cambridge university, cambridge, uk. the raw sequence data were deposited with the ncbi sra database as accession srp sequence analysis sequence analysis was performed using the ‘mothur’ software suite version, following the standard operating procedure on mothur.org. the sequences were denoised using the ampliconnoise algorithm, as implemented by mothur. sequences that were less than bases in length and/or had one of the following: > mismatches to the primer, > mismatch to the barcode regions, and homopolymers of > bases in length, were discarded. the remaining sequences were trimmed to remove primers and barcodes and aligned to the silva 16s rrna reference alignment. the uchime algorithm was used to identify chimeric sequences, which were then removed from the dataset. sequences were clustered into operational taxonomic units at a genetic distance of using the average neighbour algorithm and identified using a naïve bayesian classifier with the human oral microbiome database reference set. for those otus that could not be identified using the bayesian classifier, representative sequences were obtained in mothur and identified using blast against the homd reference set. the possible alternatives for the species identification were then provided. analysis of alpha and beta diversity the sequences for each sample were randomly sub-sampled to the same number for the alpha and beta diversity analyses. the extent of sampling of the communities was assessed using good’s non-parametric coverage estimator. the diversity of the communities was calculated using simpson’s inverse diversity index. the beta-diversity of the samples was analysed using distance matrices generated using the jaccard index and the thetayc calculator. the distance matrices were visualised using principal coordinates analysis plots generated in r. statistical analysis analysis of molecular variance, as implemented in mothur, was used to determine if there were statistically significant differences in jaccard index and thetayc distances between saliva and biofilm samples and between biofilms from different time points, subject panels, and incubation times. the mean relative abundances of phyla, genera, and species-level otus in biofilms was determined as follows: the proportions of sequences assigned to a particular taxon were calculated for biofilms derived from different panels at a single time point only after 7- and 14-day incubations. the pma-treated saliva samples were used for statistical comparisons of the taxonomic composition of saliva to biofilms. to determine if there were significant differences in otu richness and diversity between the incubation times, paired t-tests were performed in r.  pyrosequencing a total of sequences with a mean length of bases were obtained for analysis after quality filtering, screening of the sequence alignment, and removal of chimeras. a mean of sequences were obtained per sample. one replicate biofilm sample derived from panel after seven days of incubation was not included in the sequencing run due to poor pcr amplification. the number of otus detected in individual biofilm and pooled saliva samples ranged between and the mean number of otus detected was in the biofilms and in the saliva samples. a table summarising the alpha diversity of the biofilms and saliva is shown in additional file the number of otus detected in the seven-day biofilms was significantly higher than in the 14-day biofilms. however, there was no significant difference in diversity between the seven- and 14-day biofilms. reproducibility of biofilms and shifts in biofilm otu composition over time comparison of the community membership and structure of biofilms using principal coordinates analysis plots indicated that replicate biofilm samples, derived from the same saliva pool after the same incubation time but harvested from different pegs, were highly similar and clustered together. in addition, biofilms derived from the same panel at different times were similar and amova tests found there to be no statistically significant difference between the time points. however, there were significant differences in both the membership and structure of biofilms derived from the three different subject panels. the most significant differences by amova were between biofilms from panel and panel. interestingly, pcoa indicated that the dissimilarity in community structure between panels was greater after 14 days than after seven days.figure principal coordinates analysis of biofilms derived from different panels after different incubation times. plots are based on community membership using the jaccard index and community structure using the thetayc calculator. blue - panel 1; red - panel 2; black - panel labels indicate the incubation time. a: pc1 =  % of variance, pc2 =  % of variance. b: pc1 =  % of variance, pc2 =  % of variance. principal coordinates analysis of biofilm replicates from panel at different time points. plots are based on community membership using the jaccard index and community structure using the thetayc calculator. blue - time 1; red - time 2; black - time a: pc1 =  % of variance, pc2 =  % of variance. b: pc1 =  % of variance, pc2 =  % of variance. there was a directional shift along the axes in the pcoa plots between seven and 14 days of incubation for biofilms derived from all three panels. amova confirmed that there was a significant difference between the 7-day and 14-day incubations both in terms of membership and structure. analysis using lefse identified a total of otus that were significantly differentially abundant between incubation times. the identities of the otus with lda effect size scores of > are shown in figure  an otu identified as veillonella parvula was most strongly associated with the 7-day incubations, whilst parvimonas micra was most strongly associated with the 14-day incubations.figure linear discriminant analysis effect size analysis showing those otus that were significantly differentially abundant between seven- and 14-day incubation times, ranked by effect size. hierarchical cluster analysis showed that the biofilm samples clustered by panel in a dendrogram based on community membership and predominantly by incubation time in a dendrogram based on structure. otu-based comparisons of biofilms to saliva comparison of saliva samples with the biofilms using pcoa plots revealed a differing community membership and structure. amova tests confirmed that there were significant differences in both the membership and structure of saliva compared to seven and 14-day biofilms. analysis using lefse identified otus that were significantly differentially abundant between saliva and biofilms that had been grown for seven days. a list of differentially abundant otus with lda effect size scores of > is shown in figure  an otu identified as neisseria flavescens/subflava was most strongly associated with saliva, whilst veillonella parvula was most strongly associated with the biofilms.figure principal coordinates analysis of saliva and seven-day biofilms. plots are based on community membership using the jaccard index and community structure using the thetayc calculator. blue - biofilms; red - saliva. labels indicate the panel number. a: pc1 =  % of variance, pc2 =  % of variance. b: pc1 =  % of variance, pc2 =  % of variance. linear discriminant analysis effect size analysis showing those otus that were significantly differentially abundant between saliva and seven-day biofilms, ranked by effect size. comparison of pma-treated and untreated saliva samples comparisons using pcoa showed that pma-treated and untreated saliva were similar and samples clustered by the panel from which they were obtained, rather than by treatment. there were no significant differences in community membership or structure between pma-treated and untreated saliva using amova tests. taxonomic composition of the biofilms the predominant phyla detected in all of the biofilms in order of mean relative abundance were: firmicutes, bacteroidetes, synergistetes, fusobacteria, proteobacteria and actinobacteria. other phyla that were detected in minor relative proportions, and not in every sample, included: sr spirochaetes, tm and tenericutes. a total of genera were detected in the biofilms, the most abundant of which were: veillonella, streptococcus and prevotella in seven-day biofilms, and streptococcus, prevotella and parvimonas in 14-day biofilms. figure  shows the relative abundances of the predominant genera detected in the biofilms after seven and 14 days incubation. the relative abundances of the predominant genera detected in the saliva samples, from which biofilms were derived, are shown in additional file a table detailing all of the taxa identified down to the species level, and their relative proportions in individual biofilm and saliva samples, can be found in additional file figure predominant genera detected in the biofilms. the graph shows the mean relative abundances of genera that were detected in seven and 14-day biofilms derived from three panels. genera shown are those with mean relative abundances of > 1%. error bars show the standard error of the mean. discussion this study has demonstrated that complex oral biofilms derived from a natural saliva inoculum can be successfully grown and maintained using the cbd. the results showed that the biofilms had a richness and diversity close to that of the pooled saliva inocula, with a mean of species-level otus detected per biofilm sample compared to a mean of in saliva. there was a significant difference in community membership and structure between the saliva and the biofilms, with some otus detected in saliva not present in the biofilms. this is not surprising because the bacterial composition of saliva is known to differ to that of dental plaque. saliva, however, does include representatives of the various surfaces found in the mouth and is thus a useful inoculum for biofilms. hydroxyapatite-coated pegs were used for the purpose of mimicking the teeth in order to obtain biofilms with a similar composition to plaque. another possible reason for differences between the inocula and the biofilms is that specific nutrients or growth factors required by certain species could have been absent in the growth medium used. a brain-heart infusion based medium was chosen in this study because it has been successfully used to cultivate a broad range of fastidious and non-fastidious oral bacteria. the bhi was supplemented with mucin, vitamin k, and haemin, as some oral species grow poorly or not at all in the absence of one or more of these substances. for example, a number of black-pigmented species of prevotella and porphyromonas require haemin and vitamin k for growth. in addition, hog gastric mucin, a high molecular weight glycoprotein, has been shown to support the growth of mixed communities of oral bacteria when used as the principal source of carbon and energy. future work could investigate the use of different media, such as an artificial saliva-based medium, with this model. another reason that certain salivary species may have been lost is that the anaerobic atmosphere in which the cbd was incubated would have selected against the growth of aerobic species. the absence of host immune cells and molecules in the second phase of growth may also have had an impact on the community composition. nevertheless, a highly diverse community of oral bacteria was maintained which included the genera known to be predominant in plaque and also a variety of fastidious and uncultivated taxa, such as un-named bacteroidetes, lachnospiraceae, clostridiales and peptostreptococcaceae species. both saliva and biofilm samples were treated with propidium monoazide prior to dna extraction to avoid detection of bacterial cells that were non-viable; this method has been shown to prevent pcr amplification of dna from dead or damaged cells. interestingly, there was no significant difference between pma-treated and untreated saliva in terms of community membership and structure. this suggests that the vast majority of taxa detected in the saliva samples were viable. this could be due to the rapid processing of the samples performed in order to avoid loss of cell viability. in addition, human saliva has been shown to contain dnase i produced by the parotid glands, which could rapidly break down extracellular bacterial dna from dead cells. the ability to grow biofilms with a richness and diversity close to that detected in oral habitats, and that include previously uncultivated taxa, is a major strength of this model over others. for instance, the well-established zürich biofilm model has been developed with defined inocula consisting of five, or more recently, cultivable species. defined biofilms consisting of a low number of selected cultivable species are less representative of the in-vivo ecosystem and may, therefore, be less accurate in predicting the efficacy of an active agent. another study recently reported the development of an in-vitro biofilm model in which a natural saliva inoculum was used and the composition of the samples determined by pyrosequencing of 16s rrna genes. whilst the authors also reported a high microbial diversity, the incubation times used were relatively short and this may explain why certain slow-growing oral taxa, including fretibacterium spp. and tannerella spp. were not detected, whilst streptococci were dominant with s. vestibularis constituting approximately 40% of the communities. the taxonomic composition of the biofilms grown in this study was similar to that of dental plaque. streptococcus, veillonella and prevotella were the predominant genera detected in the biofilms, all of which have been shown to be major constituents of plaque. the otu detected with the highest mean relative abundance in the biofilms was veillonella parvula, which was also detected with the highest rank abundance in an extensive cloning and sanger sequencing study of the human oral microbiome. in addition, periodontitis-associated species, including the ‘red complex’: porphyromonas gingivalis, treponema denticola, tannerella forsythia, and the gram-positive anaerobes filifactor alocis and parvimonas micra, were all detected in the biofilms. these organisms have been strongly associated with deep periodontal pockets in individuals with severe chronic periodontitis. in the cbd model, p. micra was the organism most strongly associated with 14-day biofilms. species that had a significantly lower relative abundance in 14-day biofilms than 7-day biofilms included the streptococcal species s. cristatus and s. salivarius / vestibularis, which have previously been associated with health. species among the genera neisseria and rothia were detected at only very low proportions in the biofilms, despite being abundant in saliva. this is likely explained by the anaerobic incubation of the cbd, as these organisms grow optimally under aerobic conditions. it has been shown that the redox potential of dental plaque rapidly falls as the biofilm develops in vivo. in addition, experimental gingivitis studies have shown that plaque accumulating in the absence of oral hygiene supports the growth of increasing numbers of anaerobic species, many of which are gingivitis-associated. this study aimed to grow biofilms that were similar in composition to biofilms that would develop naturally in vivo without oral hygiene intervention, and anaerobic incubation was chosen in order to reproducibly obtain a biofilm typical of mature plaque. however, future work could examine the composition of cbd oral biofilms grown under aerobic conditions and compare them to those grown anaerobically. if using the biofilm model to screen antimicrobial agents or oral care product components, it would be useful to grow the biofilms under both aerobic and anaerobic conditions in order to determine the effect of a given substance on as diverse a range of oral taxa as possible. future studies could also compare the similarity of the in-vitro biofilms to dental plaque biofilms that form naturally in vivo in the same individuals abstaining from oral hygiene. this would further confirm that this model generates oral biofilms that are representative of those formed in vivo. the results of this study showed that the bacterial composition of the biofilms was highly reproducible for sample replicates from different pegs derived from the same saliva pool and incubated for the same length of time. moreover, the biofilms were similar in composition when derived from the same panel at different time points. after 14 days of incubation, when dissimilarity in the biofilms might have been expected to increase, the biofilms from the same panel clustered closely in the pcoa plots. the differences between biofilms derived from different panels was not surprising given the high inter-individual variation in bacterial diversity found in the normal human oral microbiome, although, an attempt was made to reduce this variability by pooling saliva from six individuals for use as the inoculum. hierarchical clustering of the biofilm samples in dendrograms indicated that the panel was the primary determinant of community membership, but that incubation time had a stronger influence on community structure. this is likely to be because the relative abundances of otus would be expected to change over time as the biofilms mature. due to the differences in both membership and structure of biofilms derived from different panels, the same panel should be used to provide the saliva inoculum in future studies that aim to compare biofilms grown under different conditions, or after exposure to different challenges e.g. antimicrobial agents. the reproducible biofilms that can be obtained using this model will enable relatively small changes in bacterial composition to be detected. this will be particularly useful for assessing the impact of oral care products that aim to manipulate or alter, rather than eradicate, plaque. this includes active agents or bacteriocin-producing probiotics that target particular taxa, or prebiotics that could promote the growth of health-associated bacteria. in the case of probiotics, the model could also be useful in helping to predict whether or not a particular strain is likely to colonise and persist within oral biofilms. in addition to determining changes in the community composition of the biofilms, future work could also investigate changes in community function using metagenomics, metabolomics and metatranscriptomics, in response to different active agents or changes in key environmental parameters.  this study has successfully developed an oral biofilm model using the cbd seeded with a natural saliva inoculum. the biofilms generated were highly complex and comprised of microbial taxa that are commonly found in dental plaque. in addition, their composition was shown to be reproducible when derived from the pooled saliva of the same panel of individuals. this model will therefore be useful for screening novel antimicrobial agents and also pre- or probiotics that aim to modify plaque composition. additional files additional file 1: table showing the alpha diversity of pooled saliva and biofilm samples. all samples were sub-sampled to sequences. p - panel; t - time point; s – saliva; sp – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate. additional file 2: dendrogram showing the similarity of biofilm and saliva samples based on community membership. p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate. additional file 3: dendrogram showing the similarity of biofilm and saliva samples based on community structure. p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate. additional file 4: bar chart of the predominant genera detected in pma-treated saliva samples. the chart shows the mean relative abundances of genera that were detected in all three of the pooled saliva samples from different panels. genera shown are those with mean relative abundances of > 1%. error bars show the standard error of the mean. additional file 5: table showing the classification of sequences to the species level in each sample. the table shows the relative abundances of the phylotypes in the different samples. p - panel; t - time point; s – saliva; s_p – pma-treated saliva; 7d - 7 days incubation; 14d - 14 days incubation; a, b, c - replicate. competing interests the authors declare that they have no competing interests. authors’ contributions jok contributed to the conception and design of the study, analysis of the data, and preparation of the manuscript. mp contributed to the conception and design of the study and critically revised the manuscript. wgw contributed to the conception and design of the study, analysis of the data, and preparation of the manuscript. all authors read and approved the final manuscript. 
this model could be used to screen both antimicrobial-containing oral care products and also novel approaches aiming to modify plaque composition, such as pre- or probiotics. the principal method of prevention is the mechanical removal of dental plaque augmented by active agents incorporated into toothpastes and mouthrinses. in-vitro assays that include complex oral bacterial biofilms are required to accurately predict the efficacy of novel active agents in vivo. a mean of and otus were detected in pooled saliva and biofilm samples, respectively. principal coordinates analysis plots based on community membership and structure showed that replicate biofilm samples were highly similar and clustered together. the specific objectives were to determine the reproducibility and stability of the model by comparing the composition of the biofilms over time derived from the same volunteers at different time points, and different panels of volunteers. there were significant differences between biofilms from different panels. in addition, there were no significant differences between biofilms derived from the same panel at different times using analysis of molecular variance. veillonella parvula, veillonella atypica/dispar/parvula and peptostreptococcus stomatis were the predominant otus detected in seven-day biofilms, whilst prevotella oralis, v. parvula and streptococcus constellatus were predominant in 14-day biofilms. the aim of this study was to develop an oral biofilm model using the calgary biofilm device seeded with a natural saliva inoculum and analysed by next generation sequencing. dental caries and periodontal disease are the commonest bacterial diseases of man and can result in tooth loss.
 dental caries and periodontal disease are the commonest bacterial diseases of man and can result in tooth loss. the principal method of prevention is the mechanical removal of dental plaque augmented by active agents incorporated into toothpastes and mouthrinses. in-vitro assays that include complex oral bacterial biofilms are required to accurately predict the efficacy of novel active agents in vivo. the aim of this study was to develop an oral biofilm model using the calgary biofilm device seeded with a natural saliva inoculum and analysed by next generation sequencing. the specific objectives were to determine the reproducibility and stability of the model by comparing the composition of the biofilms over time derived from the same volunteers at different time points, and different panels of volunteers.  pyrosequencing yielded sequences with a mean length of bases after filtering. a mean of and otus were detected in pooled saliva and biofilm samples, respectively. principal coordinates analysis plots based on community membership and structure showed that replicate biofilm samples were highly similar and clustered together. in addition, there were no significant differences between biofilms derived from the same panel at different times using analysis of molecular variance. there were significant differences between biofilms from different panels. pcoa revealed that there was a shift in biofilm composition between seven and 14 days. veillonella parvula, veillonella atypica/dispar/parvula and peptostreptococcus stomatis were the predominant otus detected in seven-day biofilms, whilst prevotella oralis, v. parvula and streptococcus constellatus were predominant in 14-day biofilms.  diverse oral biofilms were successfully grown and maintained using the cbd. biofilms derived from the same panel of volunteers were highly reproducible. this model could be used to screen both antimicrobial-containing oral care products and also novel approaches aiming to modify plaque composition, such as pre- or probiotics. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords 16s rrnabacteriasalivaplaquemicrobiomemicrobiotaissue-copyright-statement© the author 2015 