Cited by 159Oral microbiology and periodontitis researchGut microbiota and healthDental Health and Care Utilization

The DOI matches its Crossref record; 62 of the 62 references it cites were found in a registry or an index; 24 of the 53 sources its claims cite were checked for support.

Citizen science charts two major “stomatotypes” in the oral microbiome of adolescents and reveals links with habits and drinking water composition

Open access
Willis, Jesse R., González-Torres, Pedro, Pittis, Alexandros A., Bejarano, Luis A., Cozzuto, Luca, Andreu-Somavilla, Nuria, Alloza-Trabado, Miriam, Valentín, Antonia, Ksiezopolska, Ewa, Company, Carlos, Onywera, Harris, Montfort, Magda, Hermoso, Antonio, Iraola-Guzmán, Susana, Saus, Ester, Labeeuw, Annick, Carolis, Carlo, Hecht, Jochen, Ponomarenko, Julia, Gabaldón, Toni·2018·Microbiome

Abstract

BACKGROUND: The oral cavity comprises a rich and diverse microbiome, which plays important roles in health and disease. Previous studies have mostly focused on adult populations or in very young children, whereas the adolescent oral microbiome remains poorly studied. Here, we used a citizen science approach and 16S profiling to assess the oral microbiome of 1500 adolescents around Spain and its relationships with lifestyle, diet, hygiene, and socioeconomic and environmental parameters. RESULTS: Our results provide a detailed snapshot of the adolescent oral microbiome and how it varies with lifestyle and other factors. In addition to hygiene and dietary habits, we found that the composition of tap water was related to important changes in the abundance of several bacterial genera. This points to an important role of drinking water in shaping the oral microbiota, which has been so far poorly explored. Overall, the microbiome samples of our study can be clustered into two broad compositional patterns (stomatotypes), driven mostly by Neisseria and Prevotella, respectively. These patterns show striking similarities with those found in unrelated populations. CONCLUSIONS: We hypothesize that these stomatotypes represent two possible global optimal equilibria in the oral microbiome that reflect underlying constraints of the human oral niche. As such, they should be found across a variety of geographical regions, lifestyles, and ages.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

A citizen-science 16S survey of 1,319 Spanish adolescents found that oral microbiomes cluster into two broad “stomatotypes” driven mainly by Neisseria and Prevotella, with parallels to unrelated adult populations. Tap water chemistry was the investigated factor most strongly associated with changes in several oral genera, while hygiene, diet, and lifestyle showed more limited links.

Key findings

  • From 1,555 collected samples, 1,319 passed quality filters; 332 genera-level OTUs were identified, with 32 genera present in ≥75% of individuals. The top ten genera made up 84.64% of sequences and appeared in 99.6% of samples; Streptococcus was most abundant in 68% of samples (mean 22.3%), followed by Prevotella (11.9%), Haemophilus (11.4%), Neisseria (10.1%), and Veillonella (9%).
  • Samples clustered into two stomatotypes: stomatotype 1 driven by Neisseria and stomatotype 2 driven by Prevotella, with higher Haemophilus in stomatotype 1 and higher Veillonella and Streptococcus in stomatotype 2. The authors note similarities to two “coinhabiting groups” in a Japanese adult study and propose naming them Neisseria-Haemophilus and Prevotella-Veillonella stomatotypes.
  • Tap water composition was linked to oral microbiome variation: among 17 water parameters, Porphyromonas was positively associated with fluoride and sulfate, while Veillonella, Ralstonia, Rhizobium, Rhodococcus, and Pseudomonas negatively correlated with hardness, alkalinity, conductivity, sulfate, magnesium, sodium, calcium, chloride, and dry matter after boiling.
  • Hygiene and habits had measurable but narrower associations: brushing frequency negatively correlated with Gemella, Streptobacillus, Granulicatella, and Porphyromonas; braces correlated positively with many genera; alcohol consumption among 314 students was associated with higher Mycoplasma, Filifactor, Treponema, and Desulfobulbus, while 108 students who smoked showed no significant microbiome differences.

Methods

  • Citizen-science sampling of adolescents ages 13–15 and their teachers in 40 Spanish schools during Spring 2015, using oral rinse samples, a 54-question questionnaire, and saliva pH measured with test strips.
  • DNA was extracted and the V3–V4 regions of 16S rRNA were sequenced on Illumina MiSeq; reads were processed with Mothur v1.34.4, SILVA v119 alignment, and genus-level OTU grouping.
  • Microbiome composition was analyzed in R with Phyloseq; stomatotypes were defined by PAM clustering across five beta-diversity measures, with JSD used for main analyses, and associations tested using Pearson correlations, Kruskal-Wallis, ANOVA/TukeyHSD, chi-squared tests, and Bonferroni correction.

Limitations

  • The paper states that citizen science limits comprehensive clinical evaluation by experts, that saliva pH strips lack pH-meter precision, and that low-abundance genera may reflect sporadic environmental colonization rather than stable oral residents.
  • The paper notes that adolescents may have recently adopted some habits, so effects of smoking or diet may not yet be apparent, and that stomatotypes represent trends along a continuous gradient rather than discrete, well-separated entities.
  • The text makes evident that time since last tooth brushing was not controlled, though the authors argue this is unlikely to drive stomatotype differences because students in one class were sampled together and stomatotypes were not driven by school class.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1186/s40168-018-0592-3

Cites

62 of 62 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Microbiome composition. a Box plots of the relative abundances of the ten most common genera. b Stacked bars of relative abundances of the ten most common genera for all samples, showing the relative proportion of all samples made up of these ten genera. Stacked white bars are meaningless and appear due to lack of image resolution. c Donut chart showing the five most common phyla (inner ring) and the most common genera (outer ring) within each phylum with the average relative abundance per sample
Fig. 1 Microbiome composition. a Box plots of the relative abundances of the ten most common genera. b Stacked bars of relative abundances of the ten most common genera for all samples, showing the relative proportion of all samples made up of these ten genera. Stacked white bars are meaningless and appear due to lack of image resolution. c Donut chart showing the five most common phyla (inner ring) and the most common genera (outer ring) within each phylum with the average relative abundance per sampleCC BY · PubMed Central
Fig. 2 Correlations among genera (all samples). a Heatmap of correlations between relative abundances of genera. Color indicates Pearson correlation coefficient and “+” indicates a statistically significant correlation. While 332 different genera in total were detected, for the sake of visual representation, this figure shows only the 67 genera which were present in at least 1/3 of all samples (436). Correlation coefficient values for significant correlations can be found in Additional file 4. The indexes of genera within Additional file 4 are marked at every fifth position in the figure here so that names can be matched to the figure if so desired. b Co-occurrence network of the 20 most common genera. Edges indicate significant positive (red) or negative (blue) correlations between indicated genera. Edge width is proportional to Pearson correlation coefficient. Only displaying edges for coefficients of 0.25 or greater and − 0.25 or lower. The largest and smallest edge widths are shown with the corresponding absolute value of the correlation coefficient as it appears in the figure
Fig. 2 Correlations among genera (all samples). a Heatmap of correlations between relative abundances of genera. Color indicates Pearson correlation coefficient and “+” indicates a statistically significant correlation. While 332 different genera in total were detected, for the sake of visual representation, this figure shows only the 67 genera which were present in at least 1/3 of all samples (436). Correlation coefficient values for significant correlations can be found in Additional file 4. The indexes of genera within Additional file 4 are marked at every fifth position in the figure here so that names can be matched to the figure if so desired. b Co-occurrence network of the 20 most common genera. Edges indicate significant positive (red) or negative (blue) correlations between indicated genera. Edge width is proportional to Pearson correlation coefficient. Only displaying edges for coefficients of 0.25 or greater and − 0.25 or lower. The largest and smallest edge widths are shown with the corresponding absolute value of the correlation coefficient as it appears in the figureCC BY · PubMed Central
Fig. 3 Stomatotypes. a Principal Coordinates Analysis (PCoA) of samples using a Jensen-Shannon Divergence (JSD). Shows that the samples cluster into 2 groups (stomatotypes). b Boxplots of relative abundances of the five most common genera in samples with stomatotype 1 (red) and stomatotype 2 (blue). Bonferroni-adjusted p values from Wilcoxon tests between samples of stomatotypes 1 and 2 for streptococcus is 1.1e−7, while the values for the other 4 genera here were all less than 2e−16. c Co-occurence networks of 20 most common genera within samples of Stomatotypes 1 and 2 separately. Edges indicate significant positive (red) or negative (blue) correlations between indicated genera. Edge width is proportional to Pearson correlation coefficient. Only displaying edges for coefficients of 0.25 or greater. The largest and smallest edge widths are shown with the corresponding absolute value of the correlation coefficient as it appears in the figure
Fig. 3 Stomatotypes. a Principal Coordinates Analysis (PCoA) of samples using a Jensen-Shannon Divergence (JSD). Shows that the samples cluster into 2 groups (stomatotypes). b Boxplots of relative abundances of the five most common genera in samples with stomatotype 1 (red) and stomatotype 2 (blue). Bonferroni-adjusted p values from Wilcoxon tests between samples of stomatotypes 1 and 2 for streptococcus is 1.1e−7, while the values for the other 4 genera here were all less than 2e−16. c Co-occurence networks of 20 most common genera within samples of Stomatotypes 1 and 2 separately. Edges indicate significant positive (red) or negative (blue) correlations between indicated genera. Edge width is proportional to Pearson correlation coefficient. Only displaying edges for coefficients of 0.25 or greater. The largest and smallest edge widths are shown with the corresponding absolute value of the correlation coefficient as it appears in the figureCC BY · PubMed Central
Fig. 4 Gradients of abundances. Principal Coordinates Analysis (PCoA) of samples using a Jensen-Shannon Divergence (JSD) shows points with the same coordinates as in Fig. 3a. Circles indicate samples of stomatotype 1 and squares indicate samples of stomatotype 2. Colors represent abundance of the indicated genus (or the sum of abundances of indicated genera) for a given sample, where red is higher and blue is lower, with values indicated in the legend to the right. a Abundances for the genus Prevotella. b Sum of abundances of the 19 genera that were found to have significantly higher abundances in samples of stomatotype 1. c Sum of abundances of the 275 genera that did not have significantly higher abundances in samples of either stomatotype over the other. d Points without color or shape in order to display the spread of samples within the PCoA
Fig. 4 Gradients of abundances. Principal Coordinates Analysis (PCoA) of samples using a Jensen-Shannon Divergence (JSD) shows points with the same coordinates as in Fig. 3a. Circles indicate samples of stomatotype 1 and squares indicate samples of stomatotype 2. Colors represent abundance of the indicated genus (or the sum of abundances of indicated genera) for a given sample, where red is higher and blue is lower, with values indicated in the legend to the right. a Abundances for the genus Prevotella. b Sum of abundances of the 19 genera that were found to have significantly higher abundances in samples of stomatotype 1. c Sum of abundances of the 275 genera that did not have significantly higher abundances in samples of either stomatotype over the other. d Points without color or shape in order to display the spread of samples within the PCoACC BY · PubMed Central
Fig. 5 Correlations among stomatotypes. Heatmap of correlations between relative abundances of genera in samples with stomatotype 1 (left) and stomatotype 2 (right). Color indicates Pearson correlation coefficient and “+” indicates a statistically significant correlation. Highlighted boxes indicate genera pairs for which the correlation coefficient in the given stomatotype is at least 0.2 greater (red) or lower (blue) than the correlation coefficient in the other stomatotype. While 332 different genera in total were detected, for the sake of visual representation, this figure, as in Fig. 3, shows only the 67 genera which were present in at least 1/3 of all samples (436). Row and columns are ordered as in Fig. 2 and thus can also be compared with Additional file 4 in the same manner mentioned in Fig. 2
Fig. 5 Correlations among stomatotypes. Heatmap of correlations between relative abundances of genera in samples with stomatotype 1 (left) and stomatotype 2 (right). Color indicates Pearson correlation coefficient and “+” indicates a statistically significant correlation. Highlighted boxes indicate genera pairs for which the correlation coefficient in the given stomatotype is at least 0.2 greater (red) or lower (blue) than the correlation coefficient in the other stomatotype. While 332 different genera in total were detected, for the sake of visual representation, this figure, as in Fig. 3, shows only the 67 genera which were present in at least 1/3 of all samples (436). Row and columns are ordered as in Fig. 2 and thus can also be compared with Additional file 4 in the same manner mentioned in Fig. 2CC BY · PubMed Central
Fig. 6 Geographical patterns. Maps show that most of the measured ion levels follow a similar pattern to the proportion of stomatotype 1 samples. Porphyromonas had a significantly higher abundance in stomatotype 1 samples, while Veillonella had a significantly higher abundance in stomatotype 2 samples. Region names can be seen in Additional file 1. a Percentage of samples from each region that have stomatotype 1. b Mean alkalinity level per sample in each region (an example of one of the tap water measurements compared in this study). c Mean abundance of Porphyromonas per sample in each region. d Mean abundance of Veillonella per sample in each region
Fig. 6 Geographical patterns. Maps show that most of the measured ion levels follow a similar pattern to the proportion of stomatotype 1 samples. Porphyromonas had a significantly higher abundance in stomatotype 1 samples, while Veillonella had a significantly higher abundance in stomatotype 2 samples. Region names can be seen in Additional file 1. a Percentage of samples from each region that have stomatotype 1. b Mean alkalinity level per sample in each region (an example of one of the tap water measurements compared in this study). c Mean abundance of Porphyromonas per sample in each region. d Mean abundance of Veillonella per sample in each regionCC BY · PubMed Central
Fig. 7 Correlations with tap water composition. Heatmap of correlations between relative abundances of genera with measurements of various components of tap water. Samples that primarily drank from bottled water (505 out of 1319) at home were excluded here. Color indicates Pearson correlation coefficients and “+” indicates a statistically significant correlation
Fig. 7 Correlations with tap water composition. Heatmap of correlations between relative abundances of genera with measurements of various components of tap water. Samples that primarily drank from bottled water (505 out of 1319) at home were excluded here. Color indicates Pearson correlation coefficients and “+” indicates a statistically significant correlationCC BY · PubMed Central

Verification

Claims the paper makes

56 Own result · 28 Own interpretation · 15 Not checked yet · 15 Own method · 6 Supported · 5 Partly supported · 5 No full text of the sources · 3 Supported so far · 3 No basis given · accuracy not yet measured

Verify all 20 references — free account

Show all 56 claims

  • It may state two separate claims as one. 57
  • Probably not a substantial claim. 40
80 of 80
  • The paper relies on this without support.

    Questionnaire responses accurately reflect participants' habits and health.

    “Donors were asked to answer a questionnaire, including 54 questions [see Additional file 2], some of which proposed by citizens, about their health, and their dietary and hygiene habits.”
  • The paper relies on this without support.

    The chemical composition of tap water in the nearest town represents the tap water consumed by participants.

    “Hence, we decided to investigate this factor in more detail and linked our samples to the chemical composition of tap water of the nearest town, as reported in recent studies [40–42].”
  • The paper relies on this without support.

    The ionic composition of tap water from donors' hometowns, as obtained from national studies, reflects the tap water actually consumed by the donors.

    “Data on tap water hardness was obtained from several national ionic composition studies [40–42].”
  • The paper relies on this without support.

    The questionnaire responses accurately reflect participants' diet and hygiene habits.

    “Participants were asked to answer one questionnaire inquiring about aspects relevant to their hygiene and dietary habits. These questions were adapted from questionnaires available at the PhenX toolkit (consensus measures for Phenotypes and eXposures), which provides recommended standard data collection protocols for conducting biomedical research [60] and which has been recommended by the microbi”
4 of 4

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

2 of 2

Citations

Cites

Research works printed but not linked: 1

  1. Maraver F, Vitoria Miñana I, Ferreira-Pêgo C, Armijo F, Salas-Salvadó J. Magnesio en el agua de consume público y aguas minerales naturals en España y su contribución en cubrir las necesidades nutricionales. Nutr Hosp. 2015;31:2297–312.

Cited by

Sorted by citations

See all 159

Metadata

Authors
Willis, J. R., González-Torres, P., Pittis, A. A., Bejarano, L. A., Cozzuto, L., Andreu-Somavilla, N., Alloza-Trabado, M., Valentín, A., Ksiezopolska, E., Company, C., Onywera, H., Montfort, M., Hermoso, A., Iraola-Guzmán, S., Saus, E., Labeeuw, A., Carolis, C., Hecht, J., Ponomarenko, J., & Gabaldón, T.
Affiliations
  • Barcelona Institute of Science and Technology (ES)
  • University of Cape Town (ZA)
  • Institució Catalana de Recerca i Estudis Avançats (ES)
Corresponding author
Gabaldón, Toni
Venue
Microbiome
Published
2018-12
Volume
6
Issue
1
Pages
218-218
ISSN
2049-2618 (ISSN-L)
OpenAlex
W2901982083
Cited by
158
References
61
Citations by year
2026: 12 · 2025: 18 · 2024: 19 · 2023: 20 · 2022: 15 · 2021: 37 · 2020: 24 · 2019: 13
FWCI
13.63
Type
journal article
Language
en
Open access
Locations
  • Microbiome (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RECERCAT (Consorci de Serveis Universitaris de Catalunya) (publishedVersion, other-oa)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Open University of Cape Town (University of Cape Town) (submittedVersion)
  • Repositori digital de la UPF (Universitat Pompeu Fabra) (submittedVersion, cc-by)
Funders
  • Secretaría de Estado de Investigación, Desarrollo e Innovación (BFU2015-67107)
  • Ministerio de Economía, Industria y Competitividad, Gobierno de España (SEV-2012-02-08)
  • “la Caixa” Foundation (Saca la Lengua)
Topics
  • Health Sciences › Dentistry › Periodontics › Oral microbiology and periodontitis research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gut microbiota and health
  • Health Sciences › Dentistry › General Dentistry › Dental Health and Care Utilization
Keywords
oral microbiota, adolescents, tap water, Neisseria, Prevotella, diet and lifestyle, oral hygiene, socioeconomic factors, 16S rRNA gene profiling, citizen science, microbial community composition
MeSH terms
School Teachers, Adolescent, Adult, Bacteria, DNA, Bacterial, DNA, Ribosomal, Feeding Behavior, Humans, Hygiene, Life Style, Mouth, Neisseria, Phylogeny, RNA, Ribosomal, 16S, Spain, Sequence Analysis, DNA, Prevotella, Metagenomics, Drinking Water

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-11-10, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-11-10 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review