<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Oral. Health</journal-id>
<journal-title>Frontiers in Oral Health</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Oral. Health</abbrev-journal-title>
<issn pub-type="epub">2673-4842</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/froh.2021.641328</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Oral Health</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Biocultural Drivers of Salivary Microbiota in Australian Aboriginal and Torres Strait Islander Children</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Handsley-Davis</surname> <given-names>Matilda</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1088882/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Skelly</surname> <given-names>Emily</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Johnson</surname> <given-names>Newell W.</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/195988/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Kapellas</surname> <given-names>Kostas</given-names></name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Lalloo</surname> <given-names>Ratilal</given-names></name>
<xref ref-type="aff" rid="aff7"><sup>7</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Kroon</surname> <given-names>Jeroen</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1170831/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Weyrich</surname> <given-names>Laura S.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff8"><sup>8</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/830165/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide</institution>, <addr-line>Adelaide, SA</addr-line>, <country>Australia</country></aff>
<aff id="aff2"><sup>2</sup><institution>Australian Research Council Centre of Excellence in Australian Biodiversity and Heritage, University of Wollongong</institution>, <addr-line>Wollongong, NSW</addr-line>, <country>Australia</country></aff>
<aff id="aff3"><sup>3</sup><institution>Menzies Health Institute Queensland, Griffith University</institution>, <addr-line>Gold Coast, QLD</addr-line>, <country>Australia</country></aff>
<aff id="aff4"><sup>4</sup><institution>School of Dentistry and Oral Health, Griffith University</institution>, <addr-line>Gold Coast, QLD</addr-line>, <country>Australia</country></aff>
<aff id="aff5"><sup>5</sup><institution>Faculty of Dentistry, Oral and Craniofacial Sciences, King&#x00027;s College London</institution>, <addr-line>London</addr-line>, <country>United Kingdom</country></aff>
<aff id="aff6"><sup>6</sup><institution>Indigenous Oral Health Unit, Australian Research Centre for Population Oral Health, Adelaide Dental School, University of Adelaide</institution>, <addr-line>Adelaide, SA</addr-line>, <country>Australia</country></aff>
<aff id="aff7"><sup>7</sup><institution>School of Dentistry, University of Queensland</institution>, <addr-line>Brisbane, QLD</addr-line>, <country>Australia</country></aff>
<aff id="aff8"><sup>8</sup><institution>microARCH Laboratory, Department of Anthropology and Huck Institutes of Life Sciences, The Pennsylvania State University</institution>, <addr-line>University Park, PA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Gena D. Tribble, University of Texas Health Science Center at Houston, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Cameron B. Jeter, University of Texas Health Science Center at Houston, United States; Aurea Simon-Soro, University of Seville, Spain</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Laura S. Weyrich <email>lsw132&#x00040;psu.edu</email></corresp>
<corresp id="c002">Newell W. Johnson <email>n.johnson&#x00040;griffith.edu.au</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Oral Infections and Microbes, a section of the journal Frontiers in Oral Health</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>03</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>2</volume>
<elocation-id>641328</elocation-id>
<history>
<date date-type="received">
<day>14</day>
<month>12</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>15</day>
<month>02</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2021 Handsley-Davis, Skelly, Johnson, Kapellas, Lalloo, Kroon and Weyrich.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Handsley-Davis, Skelly, Johnson, Kapellas, Lalloo, Kroon and Weyrich</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license> </permissions>
<abstract><p>Australian Aboriginal and Torres Strait Islander children experience unacceptably high rates of dental caries compared to their non-Indigenous Australian counterparts. Dental caries significantly impacts the quality of life of children and their families, particularly in remote communities. While many socioeconomic and lifestyle factors impact caries risk, the central role of the oral microbiota in mediating dental caries has not been extensively investigated in these communities. Here, we examine factors that shape diversity and composition of the salivary microbiota in Aboriginal and Torres Strait Islander children and adolescents living in the remote Northern Peninsula Area (NPA) of Far North Queensland. We employed 16S ribosomal RNA amplicon sequencing to profile bacteria present in saliva collected from 205 individuals aged 4&#x02013;17 years from the NPA. Higher average microbial diversity was generally linked to increased age and salivary pH, less frequent toothbrushing, and proxies for lower socioeconomic status (SES). Differences in microbial composition were significantly related to age, salivary pH, SES proxies, and active dental caries. Notably, a feature classified as <italic>Streptococcus sobrinus</italic> increased in abundance in children who reported less frequent tooth brushing. A specific <italic>Veillonella</italic> feature was associated with caries presence, while features classified as <italic>Actinobacillus/Haemophilus</italic> and <italic>Leptotrichia</italic> were associated with absence of caries; a <italic>Lactobacillus gasseri</italic> feature increased in abundance in severe caries. Finally, we statistically assessed the interplay between dental caries and caries risk factors in shaping the oral microbiota. These data provide a detailed understanding of biological, behavioral, and socioeconomic factors that shape the oral microbiota and may underpin caries development in this group. This information can be used in the future to improve tailored caries prevention and management options for Australian Aboriginal and Torres Strait Islander children and communities.</p></abstract>
<kwd-group>
<kwd>bacteria</kwd>
<kwd>community dentistry</kwd>
<kwd>dental caries</kwd>
<kwd>ecology</kwd>
<kwd>microbiology</kwd>
</kwd-group>
<contract-num rid="cn001">APP1081320</contract-num>
<contract-num rid="cn002">FT180100407</contract-num>
<contract-num rid="cn002">CE170100015</contract-num>
<contract-sponsor id="cn001">National Health and Medical Research Council<named-content content-type="fundref-id">10.13039/501100000925</named-content></contract-sponsor>
<contract-sponsor id="cn002">Australian Research Council<named-content content-type="fundref-id">10.13039/501100000923</named-content></contract-sponsor>
<counts>
<fig-count count="2"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="62"/>
<page-count count="11"/>
<word-count count="7814"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Dental caries is a highly prevalent oral disease that severely impacts children and families&#x00027; quality of life [<xref ref-type="bibr" rid="B1">1</xref>, <xref ref-type="bibr" rid="B2">2</xref>]. Australian Aboriginal and Torres Strait Islander children have higher rates of dental caries, and of untreated caries, than non-Indigenous Australian children [<xref ref-type="bibr" rid="B3">3</xref>]. Additionally, Aboriginal and Torres Strait Islander children in remote communities experience worse oral health than their urban counterparts [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>]. For example, a 2006 survey of child caries experience in a remote Aboriginal and Torres Strait Islander community in the Northern Peninsula Area (NPA) of Far North Queensland found that caries experience in NPA children was over four times the national average [<xref ref-type="bibr" rid="B6">6</xref>]. Caries incidence in this community remained unacceptably high as of 2015 [<xref ref-type="bibr" rid="B7">7</xref>]. Risk factors for caries include lifestyle (e.g., oral hygiene, diet, fluoride exposure) and underlying host susceptibility (e.g., immune factors, prevalence of caries-promoting oral bacteria) [<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>]. Many of these factors are prevalent for Aboriginal and Torres Strait Islander children [<xref ref-type="bibr" rid="B10">10</xref>].</p>
<p>Recent research describes how caries is mediated by the oral microbiota &#x02013; the community of microorganisms inhabiting the human mouth [<xref ref-type="bibr" rid="B11">11</xref>]. Although the factors shaping Aboriginal and Torres Strait Islander child oral microbiota have not yet been explored, changes in microbiota composition associated with caries initiation and progression have previously been reported in other groups [<xref ref-type="bibr" rid="B12">12</xref>&#x02013;<xref ref-type="bibr" rid="B16">16</xref>]. Many host-related and environmental factors, including vertical transmission, diet, and antibiotic use have been reported to influence the diversity and composition of an individual&#x00027;s oral microbiota [<xref ref-type="bibr" rid="B17">17</xref>&#x02013;<xref ref-type="bibr" rid="B23">23</xref>]. In turn, these factors may contribute to broader-scale microbiota differences among members of human groups with different heritage and lifestyles [<xref ref-type="bibr" rid="B24">24</xref>, <xref ref-type="bibr" rid="B25">25</xref>]. For example, studies from Venezuela, the Philippines and Uganda have reported differences in saliva microbiota diversity and composition between Indigenous groups with hunter-gatherer lifestyles and counterparts living agricultural or industrialized lifestyles [<xref ref-type="bibr" rid="B26">26</xref>&#x02013;<xref ref-type="bibr" rid="B28">28</xref>]. Furthermore, our previous work identified systematic differences between the oral microbiota of Indigenous Australian and non-Indigenous adults, including differences linked to oral health [<xref ref-type="bibr" rid="B29">29</xref>]. These findings raise the possibility that Aboriginal and Torres Strait Islander children may have unique oral microbiota signals linked to caries development.</p>
<p>Here, we investigate the oral microbiota in children aged 4&#x02013;17 from the NPA using stimulated saliva and extensive metadata from dental examinations and participant questionnaires. Our objective was to examine whether any metadata factors were linked to salivary microbiota diversity and composition and thereby to better understand how the microbiota, lifestyle and socioeconomic factors, and oral disease interact in this population. We explored factors influencing the salivary microbiota and identified factors linked to dental caries in the absence of regular professional dental care. Our study provides important baseline data to better understand how the oral microbiota and its relationship to caries is shaped in a remote Aboriginal and Torres Strait Islander community and how the oral microbiota may be harnessed to improve Indigenous oral health. Such data will be crucial to ensure that future microbiota-aware therapies for dental caries and other oral diseases are also relevant to these communities.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and Methods</title>
<sec>
<title>Ethical Approval</title>
<p>In planning the study, extensive consultations were held with Elders and community members. Data reported in this manuscript are from a single baseline survey of oral and general health, conducted in 2015, prior to the implementation of preventive measures against dental caries. Feedback from the overall study was provided to the community in 2018 and 2019. Ethics approval was granted by the Griffith University Human Research Ethics Committee (GU Ref No: DOH/05/15/HREC); the Far North Queensland (FNQ) Human Research Ethics Committee (FNQ HREC/15QCH/39-970); the Department of Education and Training (Queensland Government) to approach participants at the schools; and the Torres and Cape Hospital and Health Service for Site Specific Approval. All surveys were conducted with the full understanding and written consent of parents/guardians of children from the three school campuses in the NPA, and with support of the Principal and teaching staff. We have worked closely with Community Health staff, and regularly consulted with the Mayor and Community over the years.</p>
</sec>
<sec>
<title>Participant Recruitment and Data Collection</title>
<p>Participating children aged 4&#x02013;17 were recruited and sampled at a single timepoint as the baseline for a larger clinical trial examining the impact of an annual caries preventive intervention in the NPA (ANZCTR no. ACTRN12615000693527, registered 3 July 2015) [<xref ref-type="bibr" rid="B4">4</xref>]. Saliva samples were collected by chewing on paraffin wax for 5 min while expectorating into a sterile cup. Total saliva volume produced was recorded and 2 mL of saliva was transferred to an OMNIgene OM-501 collection tube (DNA Genotek) and stored according to manufacturer&#x00027;s instructions. Samples were collected throughout school hours; due to complex field conditions, information on last food intake was not collected. Dental examinations were performed by trained and calibrated examiners as previously described [<xref ref-type="bibr" rid="B30">30</xref>]. For microbiota analyses, caries severity categories were assigned based on International Caries Detection and Assessment System (ICDAS) II codes [<xref ref-type="bibr" rid="B31">31</xref>] as follows: 0 = sound, 1&#x02013;2 = incipient, 3&#x02013;4 = moderate, 5&#x02013;6 = severe; caries status was assigned based on ICDAS II codes as follows: 0 = caries-free, 1&#x02013;6 = caries-active; categories were assigned to saliva samples based on the highest ICDAS II code recorded for that child. As compressed air was not available during examinations, teeth were dried with gauze prior to ICDAS assessment. Questionnaires detailing oral health behaviors, diet, emotional well-being, and oral health impact on quality of life were completed by participants or caregivers (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table 3</xref>). Specifically, the validated CHU-9D [<xref ref-type="bibr" rid="B32">32</xref>] and OHIP-14 [<xref ref-type="bibr" rid="B33">33</xref>] systems were used to collect data on general child quality of life and oral health-related quality of life.</p>
</sec>
<sec>
<title>DNA Extraction and Sequencing</title>
<p>A total of 255 saliva samples were used for DNA extraction and sequencing, from participants whose parents had consented to oral microbiota analysis. Genomic DNA was extracted from saliva samples in a clean facility at the University of Adelaide using the Roche High Pure PCR Template Preparation Kit (Roche Life Sciences). Two extraction blank controls (EBCs, i.e., empty tubes) were included for every 22 saliva samples. The V4 region of the bacterial 16S rRNA gene was amplified using uniquely barcoded reverse primers for each sample, as previously described [<xref ref-type="bibr" rid="B34">34</xref>]. No-template controls (NTCs) were processed alongside each amplification. Amplified, barcoded DNA was quantified using Qubit (ThermoFisher Scientific), pooled at equal relative concentrations, cleaned using Ampure magnetic beads (New England Biolabs), and quantified using TapeStation (Agilent). Paired-end 150 bp sequencing was performed using an Illumina MiSeq.</p>
</sec>
<sec>
<title>Sequence Data Processing</title>
<p>Data were processed using QIIME2 (v. 2018.8) [<xref ref-type="bibr" rid="B35">35</xref>]. Briefly, sequences were demultiplexed, denoised using the <italic>q2-deblur</italic> plugin [<xref ref-type="bibr" rid="B36">36</xref>], and assigned taxonomy using a classifier trained on the SILVA 132 database, selected as the most suitable taxonomic database as it contains sequences from a wide range of sample types. Key taxonomic results were also compared against the Human Oral Microbiome Database (HOMD) v15.1 (<xref ref-type="supplementary-material" rid="SM1">Supplementary Tables 4, 6</xref>). Five samples were removed from further analysis due to insufficient data or withdrawal of consent. Amplicon sequence variants (referred to herein as microbial &#x0201C;features&#x0201D;) observed &#x0003C;3 times were removed from the dataset. Detection and removal of putative contaminant features at the 0.5 threshold (i.e., features that were more prevalent in negative controls than in samples were considered contaminants) was performed using <italic>qiime2R</italic> [<xref ref-type="bibr" rid="B37">37</xref>] and <italic>decontam</italic> [<xref ref-type="bibr" rid="B38">38</xref>].</p>
</sec>
<sec>
<title>Microbiota Analysis</title>
<p>Data were analyzed using QIIME2 (v. 2019.7). Only samples with at least 30,000 sequences per sample were retained for microbiota analysis, leading to the removal of 42 saliva samples due to insufficient sequence depth; three samples from individuals who had contributed two saliva samples each were also removed, leaving a total of 205 samples. Samples with unknown or unrecorded values for a given metadata factor were removed prior to analysis, and categorical metadata factors needed at least 10 samples in each group for significance testing. Alpha diversity (within sample diversity; Faith&#x00027;s phylogenetic diversity [<xref ref-type="bibr" rid="B39">39</xref>]) and beta diversity (between sample diversity; unweighted UniFrac [<xref ref-type="bibr" rid="B40">40</xref>]) were calculated and statistically examined at the feature (amplicon sequence variant) level using the <italic>q2-diversity</italic> plugin, randomly subsampling data to 30,000 sequences per sample. These metrics were chosen because they incorporate phylogenetic as well as non-phylogenetic information about the diversity of the samples. Statistical significance was determined non-parametrically using Spearman correlation (for continuous variables) or Kruskal-Wallis (for categorical variables) [<xref ref-type="bibr" rid="B41">41</xref>] tests for alpha diversity and adonis tests [<xref ref-type="bibr" rid="B42">42</xref>, <xref ref-type="bibr" rid="B43">43</xref>] for beta diversity. Metadata categories that returned a significant or near-significant result in adonis tests were further investigated using PERMANOVA [<xref ref-type="bibr" rid="B42">42</xref>] and permdisp tests [<xref ref-type="bibr" rid="B44">44</xref>]. Features with &#x0003C;10 observations and/or present in &#x0003C;5 samples were removed prior to ANCOM testing using <italic>q2-composition ancom</italic> to identify features that differed significantly in abundance across categorical sample groups [<xref ref-type="bibr" rid="B45">45</xref>]. Figures were constructed in RStudio [<xref ref-type="bibr" rid="B46">46</xref>] using the <italic>qiime2R</italic> [<xref ref-type="bibr" rid="B37">37</xref>] and <italic>ggplot2</italic> packages [<xref ref-type="bibr" rid="B47">47</xref>], or downloaded from QIIME2 View and edited for clarity using Inkscape 2.0.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Sociodemographic Characteristics of Participants</title>
<p>A total of 205 saliva samples with corresponding metadata were used for microbiota analysis (<xref ref-type="table" rid="T1">Table 1</xref>). Although specific ethnicity data was not collected alongside individual saliva samples, within the NPA community 49.5% identified as Aboriginal Australian and Torres Strait Islander, 46.4% identified as Torres Strait Islander only, 1% identified as Aboriginal Australian only, and 3.1% identified as neither. A full list of metadata factors tested in microbiota analysis is given in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 1</xref>; further details on values recorded for continuous metadata variables are given in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 2</xref>.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Sociodemographic characteristics of study participants.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Category</bold></th>
<th valign="top" align="center"><bold>Mean</bold></th>
<th valign="top" align="center"><bold>Standard deviation</bold></th>
<th valign="top" align="center"><bold>Range</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="4"><bold>Continuous variables</bold></td>
</tr>
<tr>
<td valign="top" align="left">Age</td>
<td valign="top" align="center">8.53</td>
<td valign="top" align="center">3.53</td>
<td valign="top" align="center">4&#x02013;17</td>
</tr>
<tr>
<td valign="top" align="left">Saliva flow rate</td>
<td valign="top" align="center">7.05</td>
<td valign="top" align="center">0.49</td>
<td valign="top" align="center">5.4&#x02013;7.8</td>
</tr>
<tr>
<td valign="top" align="left">Saliva pH</td>
<td valign="top" align="center">5.77</td>
<td valign="top" align="center">3.00</td>
<td valign="top" align="center">0.5&#x02013;9.5</td>
</tr>
<tr>
<td valign="top" align="left">Total carious surfaces</td>
<td valign="top" align="center">9.68</td>
<td valign="top" align="center">9.8</td>
<td valign="top" align="center">0&#x02013;62</td>
</tr>
<tr style="border-top: thin solid #000000;border-bottom: thin solid #000000;">
<td valign="top" align="left"><bold>Category</bold></td>
<td valign="top" align="left"><bold>Values</bold></td>
<td valign="top" align="center"><bold>n</bold></td>
<td valign="top" align="center"><bold>%</bold></td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Categorical variables</bold></td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Sex</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Female</td>
<td valign="top" align="center">115</td>
<td valign="top" align="center">56.1%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Male</td>
<td valign="top" align="center">90</td>
<td valign="top" align="center">44.0%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Caries status</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Caries-free</td>
<td valign="top" align="center">17</td>
<td valign="top" align="center">8.3%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Caries-active</td>
<td valign="top" align="center">184</td>
<td valign="top" align="center">90%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">4</td>
<td valign="top" align="center">2.0%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Caries severity</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Caries-free (ICDAS 0)</td>
<td valign="top" align="center">17</td>
<td valign="top" align="center">8.3%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Incipient caries (ICDAS 1&#x02013;2)</td>
<td valign="top" align="center">37</td>
<td valign="top" align="center">18.0%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Moderate caries (ICDAS 3&#x02013;4)</td>
<td valign="top" align="center">47</td>
<td valign="top" align="center">22.9%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Severe caries (ICDAS 5&#x02013;6)</td>
<td valign="top" align="center">100</td>
<td valign="top" align="center">48.8%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">4</td>
<td valign="top" align="center">2.0%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Household size</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center">1&#x02013;5</td>
<td valign="top" align="center">86</td>
<td valign="top" align="center">42.0%</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">6&#x02013;10</td>
<td valign="top" align="center">96</td>
<td valign="top" align="center">46.8%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">More than 10</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">4.9%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">13</td>
<td valign="top" align="center">6.3%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Household employment status</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">No people work</td>
<td valign="top" align="center">18</td>
<td valign="top" align="center">8.8%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">At least one person works</td>
<td valign="top" align="center">172</td>
<td valign="top" align="center">83.9%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">15</td>
<td valign="top" align="center">7.3%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Soft drink consumption</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Yes</td>
<td valign="top" align="center">156</td>
<td valign="top" align="center">76.1%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">No</td>
<td valign="top" align="center">41</td>
<td valign="top" align="center">20.0%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">8</td>
<td valign="top" align="center">3.9%</td>
</tr>
<tr>
<td valign="top" align="left" colspan="4"><bold>Daily toothbrushing</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Less than once</td>
<td valign="top" align="center">16</td>
<td valign="top" align="center">7.8%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Once</td>
<td valign="top" align="center">34</td>
<td valign="top" align="center">16.6%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Twice</td>
<td valign="top" align="center">130</td>
<td valign="top" align="center">63.4%</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">More than twice</td>
<td valign="top" align="center">17</td>
<td valign="top" align="center">8.3%</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td/>
<td valign="top" align="left">Unknown or not recorded</td>
<td valign="top" align="center">8</td>
<td valign="top" align="center">3.9%</td>
</tr>
<tr>
<td valign="top" align="left">Total</td>
<td/>
<td valign="top" align="center">205</td>
<td valign="top" align="center">100%</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>A summary of important sociodemographic characteristics of individuals who donated the 205 saliva samples used for microbiota analysis is presented in this table. For continuous variables, the mean, standard deviation and range values are reported, rounded to two decimal places. A detailed record of data for continuous variables is presented in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 2</xref>. For continuous categories, the number of samples and percentage of the total in each possible value are reported; percentage values are rounded to one decimal place and may not add to exactly 100%</italic>.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec>
<title>Background DNA Had a Limited Effect on Salivary Microbiota</title>
<p>Because environmental and laboratory-based contamination can significantly impact microbiota studies [<xref ref-type="bibr" rid="B48">48</xref>, <xref ref-type="bibr" rid="B49">49</xref>], we used negative controls (EBCs and NTCs) to track contamination in our study. We verified significant differences in diversity (H = 62.67, <italic>p</italic> = 2.44 &#x000D7;10<sup>&#x02212;15</sup>) and composition (<italic>R</italic><sup>2</sup> = 0.21, <italic>p</italic> = 0.001) (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1A</xref>) between saliva and negative controls. We used <italic>decontam</italic> to statistically identify and remove 39 contaminant microbial features [<xref ref-type="bibr" rid="B38">38</xref>] (<xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1B</xref>; <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 3</xref>).</p>
</sec>
<sec>
<title>The Salivary Microbiota of NPA Children Is Dominated by Typical Human Oral Taxa</title>
<p>Following removal of putative contaminant features, we summarized the taxonomic composition of the 205 saliva samples that formed the core of our microbiota analysis (<xref ref-type="fig" rid="F1">Figure 1</xref>). The samples were dominated by the phyla <italic>Proteobacteria</italic> (30%), <italic>Bacteroidetes</italic> (26%), <italic>Firmicutes</italic> (25%), <italic>Actinobacteria</italic> (12%), and <italic>Fusobacteria</italic> (6%), with sequences assigned to <italic>Epsilonbacteraeoata, Spirochaetes, Patescibacteria, Tenericutes, Synergistetes, Cyanobacteria, Chloroflexi</italic>, and unassigned <italic>Bacteria</italic> making up the remaining 1%. At the genus level, the salivary microbiota was dominated by <italic>Prevotella</italic> (18%), <italic>Neisseria</italic> (14%), <italic>Haemophilus</italic> (12%), <italic>Streptococcus</italic> (9%), <italic>Rothia</italic> (8%), <italic>Veillonella</italic> (6%), <italic>Fusobacterium</italic> (4%), <italic>Alloprevotella</italic> (3%), <italic>Porphyromonas</italic> (3%), <italic>Gemella</italic> (2%), <italic>Granulicatella</italic> (2%), <italic>Leptotrichia</italic> (2%), <italic>Actinomyces</italic> (2%), and <italic>Aggregatibacter</italic> (2%), with various genera accounting for 1% or less of total sequences making up the remaining 13%. Taxonomic assignment of sequences using HOMD in place of SILVA yielded nearly identical classifications (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table 4</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Taxonomic summary of saliva samples. Bar chart summarizing the taxonomic composition of each saliva sample at the phylum level. Bars are ordered by participant age; taxonomy was assigned to 16S V4 amplicon sequences using the SILVA 132 database. Overall, samples are dominated by the phyla <italic>Proteobacteria, Bacteroidetes, Firmicutes, Actinobacteria</italic>, and <italic>Fusobacteria</italic>.</p></caption>
<graphic xlink:href="froh-02-641328-g0001.tif"/>
</fig>
</sec>
<sec>
<title>Age Significantly Impacts Salivary Microbiota in NPA Children</title>
<p>We tested all available metadata factors (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table 1</xref>) for associations with changes in microbial diversity (alpha diversity) (<xref ref-type="table" rid="T2">Table 2</xref>; <xref ref-type="fig" rid="F2">Figure 2</xref>), composition (beta diversity) (<xref ref-type="table" rid="T2">Table 2</xref>) and abundance of microbial features (<xref ref-type="table" rid="T3">Table 3</xref>). Participant age (Spearman &#x003C1; = 0.4, <italic>p</italic> = 0.0; <italic>R</italic><sup>2</sup> = 0.048, <italic>p</italic> = 0.001) and examination date (H = 26.7, <italic>p</italic> = 8 &#x000D7;10<sup>&#x02212;4</sup>; <italic>R</italic><sup>2</sup> = 0.101, <italic>p</italic> = 0.001) were significantly related to both diversity and composition (<xref ref-type="table" rid="T2">Table 2</xref>). A feature classified as uncultured <italic>Actinomyces</italic> [identified as <italic>Actinomyces</italic> using HOMD (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table 6</xref>)] was significantly associated with examination date (<xref ref-type="table" rid="T3">Table 3</xref>); this feature was observed on 16 of the 26 different examination dates and varied in abundance. Additionally, salivary pH and flow rate were significantly associated with microbiota diversity (Spearman &#x003C1; = 0.18, <italic>p</italic> = 0.011 and &#x003C1; = 0.28, <italic>p</italic> = 1 &#x000D7;10<sup>&#x02212;4</sup>, respectively) and composition (<italic>R</italic><sup>2</sup> = 0.013, <italic>p</italic> = 0.007 and <italic>R</italic><sup>2</sup> = 0.028, <italic>p</italic> = 0.001) (<xref ref-type="table" rid="T2">Table 2</xref>). However, multi-factor adonis tests indicated that both examination date (<italic>R</italic><sup>2</sup> = 0.063, <italic>p</italic> = 0.025) and salivary flow rate (<italic>R</italic><sup>2</sup> = 0.013, <italic>p</italic> = 0.005) were confounded with age although examination date retained some independent explanatory power (R<sup>2</sup> = 0.075, <italic>p</italic> = 0.001). Salivary pH remained independently significant in the multi-factor test (<italic>R</italic><sup>2</sup> = 0.011, <italic>p</italic> = 0.014) (<xref ref-type="table" rid="T2">Table 2</xref>). We accounted for age in subsequent beta diversity (i.e., composition) analyses using multi-factor adonis tests (<xref ref-type="table" rid="T2">Table 2</xref>), and report only age-adjusted results (where applicable) in the remainder of the text.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Metadata factors linked to significant differences in alpha and beta diversity across sample groups.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left" colspan="7" style="border-bottom: thin solid #000000;"><bold>Alpha diversity (Faith&#x00027;s phylogenetic diversity)</bold></th>
</tr>
<tr>
<th valign="top" align="left" colspan="2"><bold>Category</bold></th>
<th valign="top" align="left" colspan="3" style="border-bottom: thin solid #000000;"><bold>Spearman correlation test results</bold></th>
<th/>
<th/>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="2">Age</td>
<td valign="top" align="center">0.39</td>
<td valign="top" align="center" colspan="2">0.0</td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Saliva pH</td>
<td valign="top" align="center">0.18</td>
<td valign="top" align="center" colspan="2">0.011</td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Saliva flow rate</td>
<td valign="top" align="center">0.28</td>
<td valign="top" align="center" colspan="2">1 &#x000D7; 10<sup>&#x02212;4</sup></td>
<td/>
<td/>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Total carious surfaces</td>
<td valign="top" align="center">0.2</td>
<td valign="top" align="center" colspan="2">0.005</td>
<td/>
<td/>
</tr>
<tr>
<td colspan="2"/>
<td valign="top" align="left" colspan="3" style="border-bottom: thin solid #000000;"><bold>Kruskal-Wallis test results</bold></td>
<td/>
<td/>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td colspan="2"/>
<td valign="top" align="left"><bold>H</bold></td>
<td valign="top" align="left" colspan="2"><bold><italic>P</italic></bold></td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Daily toothbrushing</td>
<td valign="top" align="center">8.68</td>
<td valign="top" align="center" colspan="2">0.034</td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Examination date</td>
<td valign="top" align="center">26.7</td>
<td valign="top" align="center" colspan="2">8 &#x000D7; 10<sup>&#x02212;4</sup></td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left" colspan="2">Household employment status</td>
<td valign="top" align="center">4.79</td>
<td valign="top" align="center" colspan="2">0.029</td>
<td/>
<td/>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="2">Household size</td>
<td valign="top" align="center">6.58</td>
<td valign="top" align="center" colspan="2">0.037</td>
<td/>
<td/>
</tr>
<tr>
<td valign="top" align="left"><bold>Household employment status</bold></td>
<td valign="top" align="center" colspan="2" style="border-bottom: thin solid #000000;"><bold>Household employment status</bold></td>
<td valign="top" align="center" colspan="2" style="border-bottom: thin solid #000000;"><bold>Household employment status</bold></td>
<td valign="top" align="center" colspan="2" style="border-bottom: thin solid #000000;"><bold>Household employment status</bold></td>
</tr>
<tr>
<td valign="top" align="left">Age</td>
<td valign="top" align="center">0.046</td>
<td valign="top" align="center">0.001</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">0.043</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Examination date</td>
<td valign="top" align="center">0.101</td>
<td valign="top" align="center">0.001</td>
<td valign="top" align="center">0.075</td>
<td valign="top" align="center">0.001</td>
<td valign="top" align="center">0.093</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Saliva pH</td>
<td valign="top" align="center">0.013</td>
<td valign="top" align="center">0.007</td>
<td valign="top" align="center">0.011</td>
<td valign="top" align="center">0.014</td>
<td valign="top" align="center">0.015</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="left">Saliva flow rate</td>
<td valign="top" align="center">0.028</td>
<td valign="top" align="center">0.001</td>
<td valign="top" align="center">0.006</td>
<td valign="top" align="center">0.23</td>
<td valign="top" align="center">0.025</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Soft drink consumption</td>
<td valign="top" align="center">0.008</td>
<td valign="top" align="center">0.05</td>
<td valign="top" align="center">0.007</td>
<td valign="top" align="center">0.067</td>
<td valign="top" align="center">0.008</td>
<td valign="top" align="center">0.074</td>
</tr>
<tr>
<td valign="top" align="left">Household size</td>
<td valign="top" align="center">0.019</td>
<td valign="top" align="center">0.011</td>
<td valign="top" align="center">0.015</td>
<td valign="top" align="center">0.054</td>
<td valign="top" align="center">0.016</td>
<td valign="top" align="center">0.043</td>
</tr>
<tr>
<td valign="top" align="left">Household employment status</td>
<td valign="top" align="center">0.012</td>
<td valign="top" align="center">0.013</td>
<td valign="top" align="center">0.011</td>
<td valign="top" align="center">0.009</td>
<td valign="top" align="center">0.012</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="left">Caries status</td>
<td valign="top" align="center">0.009</td>
<td valign="top" align="center">0.039</td>
<td valign="top" align="center">0.008</td>
<td valign="top" align="center">0.052</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
</tr>
<tr>
<td valign="top" align="left">Caries severity</td>
<td valign="top" align="center">0.025</td>
<td valign="top" align="center">0.007</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.034</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
</tr>
<tr>
<td valign="top" align="left">Total carious surfaces</td>
<td valign="top" align="center">0.025</td>
<td valign="top" align="center">0.001</td>
<td valign="top" align="center">0.014</td>
<td valign="top" align="center">0.002</td>
<td valign="top" align="center">0.023</td>
<td valign="top" align="center">0.001</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>To test for significant differences in alpha diversity, a Spearman correlation test was used for numerical metadata categories; Kruskal-Wallis tests were used for categorical metadata categories. To test for significant differences in beta diversity, adonis tests were used for both numerical and categorical metadata categories. The results of each significant test for both alpha and beta diversity are displayed; categories that did not return a significant result for a given test are excluded fromthe relevant section of the table</italic>.</p>
</table-wrap-foot>
</table-wrap>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Factors linked to significant differences in alpha diversity. Scatter <bold>(A)</bold> and box <bold>(B-D)</bold> plots illustrating the range of Faith&#x00027;s phylogenetic diversity values for selected metadata categories identified as significantly linked to alpha diversity (<bold>A</bold>: age, <bold>B</bold>: daily toothbrushing behavior, <bold>C</bold>: household employment status, <bold>D</bold>: household size). Results of overall Spearman correlation (continuous variables) and Kruskal-Wallis (categorical variables) tests are also displayed, along with any statistically significant (FDR <italic>p</italic> &#x0003C;0.05) pairwise Kruskal-Wallis test comparisons between groups. All other pairwise Kruskal-Wallis test comparisons were not statistically significant.</p></caption>
<graphic xlink:href="froh-02-641328-g0002.tif"/>
</fig>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Microbial features differing significantly in abundance across groups identified by ANCOM.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Category</bold></th>
<th valign="top" align="left"><bold>Feature Taxonomy</bold></th>
<th valign="top" align="center"><bold>ANCOM W-value</bold></th>
<th valign="top" align="left"><bold>Prevalence (no. of samples detected in)</bold></th>
<th valign="top" align="left"><bold>Abundance (no. of sequences across all samples)</bold></th>
<th valign="top" align="left"><bold>Group association</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Daily toothbrushing</td>
<td valign="top" align="left"><italic>Streptococcus sobrinus</italic></td>
<td valign="top" align="center">279</td>
<td valign="top" align="left">51</td>
<td valign="top" align="left">433</td>
<td valign="top" align="left">Brush teeth less than once per day</td>
</tr>
<tr>
<td valign="top" align="left">Caries status</td>
<td valign="top" align="left"><italic>Actinobacillus porcinus</italic></td>
<td valign="top" align="center">81</td>
<td valign="top" align="left">5</td>
<td valign="top" align="left">8,649</td>
<td valign="top" align="left">Caries-free</td>
</tr>
<tr>
<td/>
<td valign="top" align="left"><italic>Leptotrichia</italic></td>
<td valign="top" align="center">26</td>
<td valign="top" align="left">11</td>
<td valign="top" align="left">121</td>
<td valign="top" align="left">Caries-free</td>
</tr>
<tr>
<td/>
<td valign="top" align="left"><italic>Veillonella</italic> uncultured organism</td>
<td valign="top" align="center">19</td>
<td valign="top" align="left">195</td>
<td valign="top" align="left">24,662</td>
<td valign="top" align="left">Caries-active</td>
</tr>
<tr>
<td valign="top" align="left">Caries severity</td>
<td valign="top" align="left"><italic>Lactobacillus gasseri</italic></td>
<td valign="top" align="center">208</td>
<td valign="top" align="left">44</td>
<td valign="top" align="left">1,660</td>
<td valign="top" align="left">Severe caries</td>
</tr>
<tr>
<td valign="top" align="left">Examination date</td>
<td valign="top" align="left"><italic>Actinomyces</italic> uncultured bacterium</td>
<td valign="top" align="center">622</td>
<td valign="top" align="left">59</td>
<td valign="top" align="left">3,978</td>
<td valign="top" align="left">Unclear</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Specific microbial features significantly associated with daily toothbrushing behavior, caries status, caries severity, and examination date are listed below. The metadata category, brief feature taxonomy as assigned using the SILVA 132 database, the reported W-value from ANCOM testing, overall prevalence (i.e., number of samples the feature was detected in), overall abundance (i.e., total number of sequences associated with the feature), and group association (the sample group in which the feature was most abundant) are displayed for each significant feature. Samples with unknown or unrecorded values for a given category were removed prior to ANCOM testing. Individual IDs and full taxonomy strings for significant features are given in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 5</xref>; comparison of taxonomic classification of the significant features using the SILVA and HOMD databases is given in <xref ref-type="supplementary-material" rid="SM1">Supplementary Table 6</xref></italic>.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec>
<title>Behavioral Factors Impact Salivary Microbiota Diversity and Composition</title>
<p>We investigated associations between the salivary microbiota and behavioral factors known to impact caries risk. The number of times per day a child reported brushing their teeth was significantly associated with alpha diversity (H = 8.68, <italic>p</italic> = 0.034) (<xref ref-type="table" rid="T2">Table 2</xref>). Children who reported brushing their teeth less than once per day harbored higher diversity and exhibited less inter-individual variation in alpha diversity than children who reported more frequent tooth brushing (<xref ref-type="fig" rid="F2">Figure 2</xref>). A feature classified as <italic>Streptococcus sobrinus</italic> was significantly more abundant in the saliva of children who reported brushing their teeth less than once per day (<xref ref-type="table" rid="T3">Table 3</xref>). However, tooth brushing was not associated with significant change in overall composition. Of dietary variables collected, only self-reported soft (i.e., carbonated) drink consumption approached a significant association with microbiota composition (<italic>R</italic><sup>2</sup> = 0.007, <italic>p</italic> = 0.067), an effect that may have been confounded by other factors such as age and caries severity (<xref ref-type="table" rid="T2">Table 2</xref>).</p>
</sec>
<sec>
<title>Socioeconomic Factors Are Linked to Salivary Microbiota</title>
<p>We used self-reported questionnaire data on household size and employment status in the child&#x00027;s household as proxies for socioeconomic status (SES). Employment status was significantly related to alpha diversity (H = 4.79, <italic>p</italic> = 0.029), with children who reported no people in their household working having higher average diversity and lower inter-individual variability than those who reported at least one person working (<xref ref-type="table" rid="T2">Table 2</xref>; <xref ref-type="fig" rid="F2">Figure 2</xref>). Household size was also generally associated with higher average diversity (H = 6.58, <italic>p</italic> = 0.037) (<xref ref-type="table" rid="T2">Table 2</xref>; <xref ref-type="fig" rid="F2">Figure 2</xref>); however, no pairwise significant differences in diversity were observed between groups. Both employment status (<italic>R</italic><sup>2</sup> = 0.011, <italic>p</italic> = 0.009; permdisp F = 8.68, <italic>p</italic> = 0.008) and household size (<italic>R</italic><sup>2</sup> = 0.015, <italic>p</italic> = 0.054) were significantly related to microbiota composition (<xref ref-type="table" rid="T2">Table 2</xref>). Specifically, significant differences in composition were found between children who reported living with 1&#x02013;5 people and those who reported living with 6&#x02013;10 (pairwise PERMANOVA FDR-corrected <italic>p</italic> = 0.039) or more than 10 people (pairwise PERMANOVA FDR-corrected <italic>p</italic> = 0.041). However, ANCOM testing did not identify any microbial features that changed significantly in abundance across these groups.</p>
</sec>
<sec>
<title>Dental Caries Is Associated With Salivary Microbiota Composition</title>
<p>Finally, we explored associations between caries and the salivary microbiota. Total number of carious surfaces was significantly associated with alpha diversity (Spearman &#x003C1; = 0.2, <italic>p</italic> = 0.005) (<xref ref-type="table" rid="T2">Table 2</xref>). Variations in composition were associated with caries status (i.e., caries-free or caries-active) (<italic>R</italic><sup>2</sup> = 0.008, <italic>p</italic> = 0.052), caries severity (<italic>R</italic><sup>2</sup> = 0.02, <italic>p</italic> = 0.034), and number of carious surfaces (<italic>R</italic><sup>2</sup> = 0.014, <italic>p</italic> = 0.002) (<xref ref-type="table" rid="T2">Table 2</xref>). Pairwise tests comparing the caries-free or incipient caries groups with the moderate or severe caries groups approached significance (pairwise PERMANOVA FDR-corrected <italic>p</italic> = 0.072). We used ANCOM differential abundance testing to identify microbial features driving these compositional changes (<xref ref-type="table" rid="T3">Table 3</xref>). The <italic>Actinobacillus porcinus</italic> feature previously associated with examination date was significantly associated with caries status and was more abundant in the caries-free group (W = 81) (<xref ref-type="table" rid="T3">Table 3</xref>). Members of the genera <italic>Leptotrichia</italic> (W = 26) and <italic>Veillonella</italic> (W = 19) were significantly more abundant in the caries-free and caries-active groups, respectively (<xref ref-type="table" rid="T3">Table 3</xref>). A microbial feature identified as <italic>Lactobacillus gasseri</italic> was significantly associated with caries severity (W = 208), specifically with the severe caries group (<xref ref-type="table" rid="T3">Table 3</xref>).</p>
<p>We further sought to understand whether caries severity interacted with other factors significantly or near-significantly related to salivary microbiota composition using multi-factor adonis tests (<xref ref-type="table" rid="T2">Table 2</xref>). Caries severity consistently explained 2.1&#x02013;3.2% of variation in the dataset and was significantly associated with composition (<italic>p</italic> &#x0003C;0.05). Self-reported soft drink consumption co-varied with caries severity (<italic>R</italic><sup>2</sup> = 0.008, <italic>p</italic> = 0.074) (<xref ref-type="table" rid="T2">Table 2</xref>). However, saliva pH (<italic>R</italic><sup>2</sup> = 0.015, <italic>p</italic> = 0.004), saliva flow rate (<italic>R</italic><sup>2</sup> = 0.025, <italic>p</italic> = 0.001), household size (<italic>R</italic><sup>2</sup> = 0.016, <italic>p</italic> = 0.043) household employment status (<italic>R</italic><sup>2</sup> = 0.012, <italic>p</italic> = 0.004), and total carious surfaces (<italic>R</italic><sup>2</sup> = 0.023, <italic>p</italic> = 0.001) were significantly related to composition, independent of caries severity (<xref ref-type="table" rid="T2">Table 2</xref>).</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>This is the first study to examine the biological, behavioral, and socioeconomic factors driving overall salivary microbiota diversity and composition in Aboriginal and Torres Strait Islander children and adolescents. Although Aboriginal and Torres Strait Islander children, especially those living in remote communities, are at high risk of caries [<xref ref-type="bibr" rid="B3">3</xref>&#x02013;<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B10">10</xref>], little is known about the role of the oral microbiota and how it may mediate disease in this group [<xref ref-type="bibr" rid="B50">50</xref>]. While several recent studies have investigated oral microbiota in Indigenous individuals around the world, many have focused primarily on characterizing differences between human groups living different lifestyles [<xref ref-type="bibr" rid="B26">26</xref>&#x02013;<xref ref-type="bibr" rid="B28">28</xref>]. Only a handful of studies have examined links between oral microbiota and oral health in Indigenous populations [<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B50">50</xref>]. In our study, we found that salivary microbiota diversity in Aboriginal and Torres Strait Islander children and adolescents from the NPA is linked to age, salivary characteristics, number of carious surfaces, toothbrushing behaviors, and SES (<xref ref-type="table" rid="T2">Table 2</xref>; <xref ref-type="fig" rid="F2">Figure 2</xref>). The composition of salivary microbiota is related to age, salivary characteristics, self-reported soft drink consumption, SES, and dental caries (<xref ref-type="table" rid="T2">Table 2</xref>). We identified microbial features that significantly varied in abundance according to examination date, toothbrushing behavior, caries status, and caries severity (<xref ref-type="table" rid="T3">Table 3</xref>). We acknowledge that the number of caries-active individuals in this study clearly outweighs the caries-free group; however, this was a population-based study, not focused on the recruitment of caries-active children. Accordingly, our results should be interpreted to identify these factors within a population where caries is highly prevalent. While saliva does not represent a singular, structured oral microbial community, it is easy and non-invasive to collect and provides a broad overview of the microbes present in the oral cavity [<xref ref-type="bibr" rid="B51">51</xref>]. Our findings contribute to a novel understanding of the mechanistic associations between biological, behavioral and socioeconomic factors, the oral microbiota, and dental caries in this population of children, with implications for the microbiota-aware treatment and prevention of oral diseases.</p>
<p>Behavioral factors, such as toothbrushing frequency and soft drink consumption, were linked to changes in microbiota diversity and composition in our dataset. The ability of microbes to grow at a given site depends on environmental factors, such as pH and oxygen or nutrient availability, which may be altered by host behavior. Less frequent brushing was linked to higher microbial alpha diversity in our study (<xref ref-type="fig" rid="F2">Figure 2</xref>). Regular tooth brushing interrupts the accumulation of microbial species on the surfaces of the oral cavity and thereby lowers the overall diversity of microbes found in saliva. Further, a feature classified as <italic>Streptococcus sobrinus</italic> was significantly more abundant in the salivary microbiota of children who reported brushing their teeth less than once per day (<xref ref-type="table" rid="T3">Table 3</xref>). <italic>S. sobrinus</italic> has long been associated with dental caries in the literature and is thought to aggravate caries when found in association with other cariogenic species such as <italic>Streptococcus mutans</italic> [<xref ref-type="bibr" rid="B52">52</xref>]. Of interest, a recent study of supragingival plaque microbiota in non-Indigenous Australian children also reported that <italic>Streptococcus</italic> abundance decreased according to tooth brushing frequency [<xref ref-type="bibr" rid="B53">53</xref>]. However, <italic>S. mutans</italic> itself was not significantly associated with dental caries or other metadata factors in our microbiota analysis. Soft drink consumption, which was widespread within our study population (<xref ref-type="table" rid="T1">Table 1</xref>), decreases environmental pH in the mouth through the introduction of free sugars that microbes ferment to acid [<xref ref-type="bibr" rid="B11">11</xref>]. The established impact of environmental pH on oral microbial communities is further supported by our result that salivary pH was significantly associated with salivary microbiota composition and diversity (<xref ref-type="table" rid="T2">Table 2</xref>). Overall, our findings suggest that daily behavioral activities linked to caries development can impact child salivary microbiota in this population. Understanding which microbes respond to specific factors such as tooth brushing could be translated into personalized medical approaches in the future, providing tailored recommendations relevant to the individual&#x00027;s microbiota composition.</p>
<p>We also found that socioeconomic factors, such as household size and employment status of household members, were associated with salivary microbiota diversity and composition in NPA children (<xref ref-type="table" rid="T2">Table 2</xref>). It is well-established that low SES increases caries risk in children [<xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B10">10</xref>]. However, whether microbial mechanisms mediate this risk is less well-understood. Johansson et al. demonstrated differences in dental plaque microbiota richness and composition between a low-SES, high-caries population in Romania and a high-SES, low-caries population in Sweden [<xref ref-type="bibr" rid="B8">8</xref>]. This comparison is confounded by large cultural, geographic, and historical differences, so it is of interest that we, for the first time, demonstrate that a similar pattern in salivary microbiota in a comparatively homogeneous population. Here, children who reported no one in their household working had less variability in microbial composition compared to the group who reported at least one person working. However, given that household size co-varied with caries severity in our dataset and that significant differences in dispersion linked to household employment status were identified by permdisp testing, these findings require further investigation to determine whether these factors are mechanistically linked to caries.</p>
<p>In our study population, multiple measures of caries were significantly associated with salivary microbiota composition (<xref ref-type="table" rid="T2">Table 2</xref>), supporting previous findings that salivary microbiota correlates with caries status [<xref ref-type="bibr" rid="B14">14</xref>]. In particular, pairwise comparisons of the caries-free or incipient caries groups to the moderate or severe caries groups approached significance, suggesting that progression of caries to more advanced stages may involve a distinct shift in microbiota composition. Studies of dental plaque have identified shifts in microbial community composition related to caries progression [<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B54">54</xref>]; such patterns may be less obvious, but still present, when sampling saliva rather than tooth surfaces [<xref ref-type="bibr" rid="B55">55</xref>, <xref ref-type="bibr" rid="B56">56</xref>]. A feature classified as <italic>Lactobacillus gasseri</italic> was significantly associated with caries severity in our dataset and was most abundant in the severe caries group (<xref ref-type="table" rid="T3">Table 3</xref>). <italic>L. gasseri</italic> has previously been identified in the human mouth [<xref ref-type="bibr" rid="B57">57</xref>]; as a lactic acid bacterium, it likely participates in sugar fermentation and hence, caries promotion. A <italic>Veillonella</italic> feature was significantly more abundant in caries-active samples, while features classified as <italic>Leptotrichia</italic> and <italic>Actinobacillus porcinus</italic> were more abundant in the caries-free group (<xref ref-type="table" rid="T3">Table 3</xref>). Numerous studies have reported that <italic>Veillonella</italic> or <italic>Veillonellaceae</italic> species are associated with dental caries [<xref ref-type="bibr" rid="B12">12</xref>&#x02013;<xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B16">16</xref>, <xref ref-type="bibr" rid="B54">54</xref>]. <italic>Veillonella</italic> species use lactic acid as their primary energy source and therefore are closely associated with caries-promoting species that produce lactate [<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B58">58</xref>]. Because of this association, <italic>Veillonella</italic> levels in plaque and saliva have been suggested as a biomarker for future caries risk even at apparently healthy sites [<xref ref-type="bibr" rid="B16">16</xref>, <xref ref-type="bibr" rid="B54">54</xref>, <xref ref-type="bibr" rid="B58">58</xref>]. Of interest, a <italic>Veillonella</italic> species was previously identified as significantly more abundant in the dental plaque microbiota of Canadian First Nations children with early childhood caries compared to caries-free counterparts [<xref ref-type="bibr" rid="B13">13</xref>]. The association of <italic>Leptotrichia</italic> species with oral health and disease is less clear; some species in this genus may be disease-associated and others health-associated [<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B54">54</xref>]. The importance of the <italic>Actinobacillus porcinus</italic> feature in oral health is also difficult to interpret, as this species is not typically found in the human oral microbiota. However, the 16S amplicon sequence associated with this feature was classified as <italic>Haemophilus</italic> using the Human Oral Microbiome Database (<xref ref-type="supplementary-material" rid="SM1">Supplementary Table 6</xref>). Other recent publications have reported that <italic>Haemophilus</italic> is found in higher relative abundance in the saliva of caries-free children and adults compared to those with caries [<xref ref-type="bibr" rid="B56">56</xref>, <xref ref-type="bibr" rid="B59">59</xref>], although this association is not universal [<xref ref-type="bibr" rid="B13">13</xref>]. Better characterization of the oral health relevance of <italic>Leptotrichia</italic> and <italic>Haemophilus</italic> strains present in the NPA child population could be useful in understanding microbial oral health and informing new therapeutic strategies.</p>
<p>While this study is the first of its kind, there are several limitations. While our study used saliva samples to profile the oral microbiota, samples of plaque biofilm from specific tooth sites might reveal closer associations between the microbiota and disease [<xref ref-type="bibr" rid="B51">51</xref>, <xref ref-type="bibr" rid="B60">60</xref>]. However, plaque collection was not practicable under the field conditions of our study. Saliva samples, especially after a period of chewing wax, which would dislodge much adherent biofilm [<xref ref-type="bibr" rid="B61">61</xref>], give an overview of oral microbiota [<xref ref-type="bibr" rid="B51">51</xref>, <xref ref-type="bibr" rid="B55">55</xref>, <xref ref-type="bibr" rid="B60">60</xref>, <xref ref-type="bibr" rid="B62">62</xref>] and can be related to overall caries experience and activity of the individual child. In addition, some metadata information was collected using self-administered questionnaires, data on last meal before sample collection was not collected, and our dataset lacked sampling controls to detect contamination at the sampling site. In relation to this last point, investigation of sampling controls by our group in later years of the clinical trial demonstrated minimal overlap between saliva and sampling controls, indicating negligible contamination of saliva during sampling [<xref ref-type="bibr" rid="B50">50</xref>]. We also saw that some significant factors co-varied, making it difficult to distinguish the specific mechanisms underpinning each factor, i.e., whether these factors work together or independently to influence microbiota composition, or are both correlated with some other, unmeasured factor. For example, examination date was significant in explaining variation in both alpha and beta diversity (<xref ref-type="table" rid="T2">Table 2</xref>) and was associated with a significant change in abundance of a specific <italic>Actinomyces</italic> feature (<xref ref-type="table" rid="T3">Table 3</xref>). However, further investigation suggested that examination date was partially confounded with age, meaning that these associations may not be directly linked to examination date but an artifact of collecting samples from children of different ages at different schools on different days. In another example, we found that soft drink consumption co-varied with caries severity (<xref ref-type="table" rid="T2">Table 2</xref>), making it difficult to discern whether either of these factors acts independently. Further research into each of these factors is needed to better understand their contributions to oral microbiota diversity and composition.</p>
<p>For the first time, we describe the salivary microbiota of Aboriginal and Torres Strait Islander children living in a remote location with limited access to dental care. We identify relationships between the salivary microbiota, dental caries, and known caries risk factors such as behavioral activities and SES. Given the importance of the oral microbiota for oral health, refining our understanding of oral microbial communities and how they mediate oral health and disease could be key to informing better treatment and prevention strategies, particularly in populations at high risk of oral disease. This understanding may be especially important for the oral health of Australian Aboriginal and Torres Strait Islander peoples, as early evidence suggests a distinctive relationship between these peoples and their associated oral microbes [<xref ref-type="bibr" rid="B29">29</xref>]. Datasets such as ours form a baseline for longitudinal studies of caries prevention and will be key in ensuring that new microbiome-based or microbiome-aware therapies are also applicable to Indigenous communities and do not damage or disrupt Indigenous microbiota. Future research toward this goal could include the investigation of different sample types, such as dental plaque, that allow for a more structured view of oral microbial communities; inclusion of more Aboriginal and Torres Strait Islander communities across Australia; employment of more precise sequencing techniques, such as shotgun metagenomics, to obtain species-level identification of oral microbes; and collection of more detailed metadata to support a finer-scale understanding of the relationship between oral microbiota and, for example, diet.</p>
</sec>
<sec sec-type="data-availability-statement" id="s5">
<title>Data Availability Statement</title>
<p>The datasets presented in this article are not readily available because of ethical and privacy constraints and respect for principles of Indigenous data sovereignty. Data may be made available for further analyses on a case-by-case basis subject to HREC review and community approval. Requests to access the datasets should be directed to Newell W. Johnson, <email>n.johnson&#x00040;griffith.edu.au</email>.</p>
</sec>
<sec id="s6">
<title>Ethics Statement</title>
<p>The studies involving human participants were reviewed and approved by Griffith University Human Research Ethics Committee, Far North Queensland (FNQ) Human Research Ethics Committee, Department of Education and Training (Queensland Government), and Torres and Cape Hospital and Health Service. Written informed consent to participate in this study was provided by the participants&#x00027; legal guardian/next of kin.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>MH-D performed data analysis and interpretation and wrote the manuscript. ES performed lab work and data processing and critically revised the manuscript. NJ contributed to study conception and design, data acquisition and interpretation, and critically revised the manuscript. KK and RL contributed to data acquisition and critically revised the manuscript. JK contributed to study conception and design, ethical approvals, data acquisition, and critically revised the manuscript. LW contributed to study conception and design, contributed to data interpretation, and critically revised the manuscript. All authors gave their final approval and agree to be accountable for all aspects of the work.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<ack><p>We pay our respects to the traditional owners of the lands on which this research took place and on which we work: the Anggamuthi, Atambay, Wuthathi, Yadhaykenu and Gudang peoples of the Northern Peninsula Area (NPA); the Kaurna people of the Adelaide Plains; the Yugambeh/Kombumerri peoples of the Gold Coast; and the Turrbal and Jagara people of the Brisbane (Meanjin) area. We are grateful to the schoolchildren, parents and guardians of the Northern Peninsula Area of Far North Queensland, to the teaching staff of the NPA State College and to the Community Health staff for their contributions and participation in the project. We also gratefully acknowledge the practical contributions of Associate Professor Robyn Boase, Elizabeth Cobbledick, Joel Rogers, and other members of the epidemiological and treatment teams. We thank Dr. Santosh Kumar Tadakamadla for assistance with data management and processing, Dr. Raymond Tobler with statistical interpretation, and Ms. Caitlin Selway, Dr. Raphael Eisenhofer, Dr. Justine Debelius, Dr. Nicholas Bokulich, Mr. Colin Brislawn, and other users and moderators of the QIIME2 forum (<ext-link ext-link-type="uri" xlink:href="https://forum.qiime2.org">https://forum.qiime2.org</ext-link>) with QIIME2 troubleshooting.</p></ack>
<sec sec-type="supplementary-material" id="s8">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/froh.2021.641328/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/froh.2021.641328/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gilchrist</surname> <given-names>F</given-names></name> <name><surname>Marshman</surname> <given-names>Z</given-names></name> <name><surname>Deery</surname> <given-names>C</given-names></name> <name><surname>Rodd</surname> <given-names>HD</given-names></name></person-group>. <article-title>The impact of dental caries on children and young people: what they have to say?</article-title> <source>Int J Paediatr Dent.</source> (<year>2015</year>) <volume>25</volume>:<fpage>327</fpage>&#x02013;<lpage>38</lpage>. <pub-id pub-id-type="doi">10.1111/ipd.12186</pub-id><pub-id pub-id-type="pmid">26153526</pub-id></citation></ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Butten</surname> <given-names>K</given-names></name> <name><surname>Johnson</surname> <given-names>NW</given-names></name> <name><surname>Hall</surname> <given-names>KK</given-names></name> <name><surname>Toombs</surname> <given-names>M</given-names></name> <name><surname>King</surname> <given-names>N</given-names></name> <name><surname>O&#x00027;Grady</surname> <given-names>K-AF</given-names></name></person-group>. <article-title>Impact of oral health on Australian urban Aboriginal and Torres Strait Islander families: a qualitative study</article-title>. <source>Int J Equity Health.</source> (<year>2019</year>) <volume>18</volume>:<fpage>34</fpage>. <pub-id pub-id-type="doi">10.1186/s12939-019-0937-y</pub-id><pub-id pub-id-type="pmid">30777079</pub-id></citation></ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Ha</surname> <given-names>DH</given-names></name> <name><surname>Roberts-Thomson</surname> <given-names>KF</given-names></name> <name><surname>Arrow</surname> <given-names>P</given-names></name> <name><surname>Peres</surname> <given-names>KG</given-names></name> <name><surname>Do</surname> <given-names>LG</given-names></name></person-group>. <article-title>Children&#x00027;s oral health status in Australia, 2012-2014. In: Do LG, Spencer AJ, editors</article-title>. <source>National Survey of Child Oral Health 2012-2014</source>. Adelaide: University of Adelaide Press (<year>2016</year>). p. <fpage>86</fpage>&#x02013;<lpage>152</lpage>.</citation></ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lalloo</surname> <given-names>R</given-names></name> <name><surname>Kroon</surname> <given-names>J</given-names></name> <name><surname>Tut</surname> <given-names>O</given-names></name> <name><surname>Kularatna</surname> <given-names>S</given-names></name> <name><surname>Jamieson</surname> <given-names>LM</given-names></name> <name><surname>Wallace</surname> <given-names>V</given-names></name> <etal/></person-group>. <article-title>Effectiveness, cost-effectiveness and cost-benefit of a single annual professional intervention for the prevention of childhood dental caries in a remote rural Indigenous community</article-title>. <source>BMC Oral Health.</source> (<year>2015</year>) <volume>15</volume>:<fpage>99</fpage>. <pub-id pub-id-type="doi">10.1186/s12903-015-0076-9</pub-id></citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Roberts-Thomson</surname> <given-names>KF</given-names></name> <name><surname>Kapellas</surname> <given-names>K</given-names></name> <name><surname>Ha</surname> <given-names>D</given-names></name> <name><surname>Jamieson</surname> <given-names>LM</given-names></name> <name><surname>Arrow</surname> <given-names>P</given-names></name> <name><surname>Do</surname> <given-names>LG</given-names></name></person-group>. <article-title>Oral health status and behaviours of Indigenous Australian children. In: Do LG, Spencer AJ, editors</article-title>. <source>Oral health of Australian children: The National Child Oral Health Study 2012&#x02013;14</source>. Adelaide: University of Adelaide Press (<year>2016</year>). p. <fpage>264</fpage>&#x02013;<lpage>87</lpage>.</citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hopcraft</surname> <given-names>M</given-names></name> <name><surname>Chow</surname> <given-names>W</given-names></name></person-group>. <article-title>Dental caries experience in aboriginal and torres strait islanders in the Northern Peninsula Area, Queensland</article-title>. <source>Aust Dent J.</source> (<year>2007</year>) <volume>52</volume>:<fpage>300</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1111/j.1834-7819.2007.tb00506.x</pub-id></citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Kroon</surname> <given-names>J</given-names></name> <name><surname>Lalloo</surname> <given-names>R</given-names></name> <name><surname>Tadakamadla</surname> <given-names>SK</given-names></name> <name><surname>Johnson</surname> <given-names>NW</given-names></name></person-group>. <article-title>Dental caries experience in children of a remote Australian Indigenous community following passive and active preventive interventions</article-title>. <source>Community Dent Oral Epidemiol.</source> (<year>2019</year>) 47: 1&#x02212;7. <pub-id pub-id-type="doi">10.1111/cdoe.12486</pub-id></citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johansson</surname> <given-names>I</given-names></name> <name><surname>Witkowska</surname> <given-names>E</given-names></name> <name><surname>Kaveh</surname> <given-names>B</given-names></name> <name><surname>Lif</surname> <given-names>Holgerson P</given-names></name> <name><surname>Tanner</surname> <given-names>ACR</given-names></name></person-group>. <article-title>The microbiome in populations with a low and high prevalence of caries</article-title>. <source>J Dent Res.</source> (<year>2016</year>) <volume>95</volume>:<fpage>80</fpage>&#x02013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1177/0022034515609554</pub-id><pub-id pub-id-type="pmid">26442950</pub-id></citation></ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fernando</surname> <given-names>S</given-names></name> <name><surname>Kumar</surname> <given-names>S</given-names></name> <name><surname>Bakr</surname> <given-names>M</given-names></name> <name><surname>Speicher</surname> <given-names>D</given-names></name> <name><surname>Lea</surname> <given-names>R</given-names></name> <name><surname>Scuffham</surname> <given-names>PA</given-names></name> <etal/></person-group>. <article-title>Children&#x00027;s untreated decay is positively associated with past caries experience and with current salivary loads of mutans Streptococci; negatively with self-reported maternal iron supplements during pregnancy: a multifactorial analysis</article-title>. <source>J Public Health Dent.</source> (<year>2019</year>) <volume>79</volume>:<fpage>109</fpage>&#x02013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1111/jphd.12301</pub-id><pub-id pub-id-type="pmid">30551255</pub-id></citation></ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Butten</surname> <given-names>K</given-names></name> <name><surname>Johnson</surname> <given-names>N</given-names></name> <name><surname>Anderson</surname> <given-names>J</given-names></name> <name><surname>Toombs</surname> <given-names>M</given-names></name> <name><surname>King</surname> <given-names>N</given-names></name> <name><surname>O&#x00027;Grady</surname> <given-names>K</given-names></name> <etal/></person-group>. <article-title>Risk factors for oral health in young, urban, Aboriginal and Torres Strait Islander children</article-title>. <source>Aust Dent J.</source> (<year>2019</year>) <volume>64</volume>:<fpage>72</fpage>&#x02013;<lpage>81</lpage>. <pub-id pub-id-type="doi">10.1111/adj.12662</pub-id><pub-id pub-id-type="pmid">30375649</pub-id></citation></ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tanner</surname> <given-names>ACR</given-names></name> <name><surname>Kressirer</surname> <given-names>CA</given-names></name> <name><surname>Rothmiller</surname> <given-names>S</given-names></name> <name><surname>Johansson</surname> <given-names>I</given-names></name> <name><surname>Chalmers</surname> <given-names>NI</given-names></name></person-group>. <article-title>The caries microbiome: implications for reversing dysbiosis</article-title>. <source>Adv Dent Res.</source> (<year>2018</year>) <volume>29</volume>:<fpage>78</fpage>&#x02013;<lpage>85</lpage>. <pub-id pub-id-type="doi">10.1177/0022034517736496</pub-id><pub-id pub-id-type="pmid">29355414</pub-id></citation></ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gross</surname> <given-names>EL</given-names></name> <name><surname>Beall</surname> <given-names>CJ</given-names></name> <name><surname>Kutsch</surname> <given-names>SR</given-names></name> <name><surname>Firestone</surname> <given-names>ND</given-names></name> <name><surname>Leys</surname> <given-names>EJ</given-names></name> <name><surname>Griffen</surname> <given-names>AL</given-names></name></person-group>. <article-title>Beyond Streptococcus mutans: dental caries onset linked to multiple species by 16S rRNA community analysis</article-title>. <source>PLoS ONE.</source> (<year>2012</year>) <volume>7</volume>:<fpage>e47722</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0047722</pub-id><pub-id pub-id-type="pmid">23091642</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Agnello</surname> <given-names>M</given-names></name> <name><surname>Marques</surname> <given-names>J</given-names></name> <name><surname>Cen</surname> <given-names>L</given-names></name> <name><surname>Mittermuller</surname> <given-names>B</given-names></name> <name><surname>Huang</surname> <given-names>A</given-names></name> <name><surname>Chaichanasakul</surname> <given-names>Tran N</given-names></name> <etal/></person-group>. <article-title>Microbiome associated with severe caries in Canadian first nations children</article-title>. <source>J Dent Res.</source> (<year>2017</year>) <volume>96</volume>:<fpage>1378</fpage>&#x02013;<lpage>85</lpage>. <pub-id pub-id-type="doi">10.1177/0022034517718819</pub-id><pub-id pub-id-type="pmid">28709393</pub-id></citation></ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eriksson</surname> <given-names>L</given-names></name> <name><surname>Holgerson</surname> <given-names>PL</given-names></name> <name><surname>Johansson</surname> <given-names>I</given-names></name></person-group>. <article-title>Saliva and tooth biofilm bacterial microbiota in adolescents in a low caries community</article-title>. <source>Sci Rep.</source> (<year>2017</year>) <volume>7</volume>:<fpage>5861</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-06221-z</pub-id><pub-id pub-id-type="pmid">28724921</pub-id></citation></ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Espinoza</surname> <given-names>JL</given-names></name> <name><surname>Harkins</surname> <given-names>DM</given-names></name> <name><surname>Torralba</surname> <given-names>M</given-names></name> <name><surname>Gomez</surname> <given-names>A</given-names></name> <name><surname>Highlander</surname> <given-names>SK</given-names></name> <name><surname>Jones</surname> <given-names>MB</given-names></name> <etal/></person-group>. <article-title>Supragingival plaque microbiome ecology and functional potential in the context of health and disease</article-title>. <source>mBio.</source> (<year>2018</year>) <volume>9</volume>:<fpage>e01631</fpage>&#x02013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1128/mBio.01631-18</pub-id><pub-id pub-id-type="pmid">30482830</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>L</given-names></name> <name><surname>Chen</surname> <given-names>X</given-names></name> <name><surname>Wang</surname> <given-names>Y</given-names></name> <name><surname>Jiang</surname> <given-names>W</given-names></name> <name><surname>Wang</surname> <given-names>S</given-names></name> <name><surname>Ling</surname> <given-names>Z</given-names></name> <etal/></person-group>. <article-title>Dynamic alterations in salivary microbiota related to dental caries and age in preschool children with deciduous dentition: a 2-year follow-up study</article-title>. <source>Front Physiol.</source> (<year>2018</year>) <volume>9</volume>:<fpage>342</fpage>. <pub-id pub-id-type="doi">10.3389/fphys.2018.00342</pub-id><pub-id pub-id-type="pmid">29670544</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blekhman</surname> <given-names>R</given-names></name> <name><surname>Goodrich</surname> <given-names>JK</given-names></name> <name><surname>Huang</surname> <given-names>K</given-names></name> <name><surname>Sun</surname> <given-names>Q</given-names></name> <name><surname>Bukowski</surname> <given-names>R</given-names></name> <name><surname>Bell</surname> <given-names>JT</given-names></name> <etal/></person-group>. <article-title>Host genetic variation impacts microbiome composition across human body sites</article-title>. <source>Genome Biol.</source> (<year>2015</year>) <volume>16</volume>:<fpage>191</fpage>. <pub-id pub-id-type="doi">10.1186/s13059-015-0759-1</pub-id><pub-id pub-id-type="pmid">26374288</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Corby</surname> <given-names>PM</given-names></name> <name><surname>Bretz</surname> <given-names>WA</given-names></name> <name><surname>Hart</surname> <given-names>TC</given-names></name> <name><surname>Schork</surname> <given-names>NJ</given-names></name> <name><surname>Wessel</surname> <given-names>J</given-names></name> <name><surname>Lyons-Weiler</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Heritability of oral microbial species in caries-active and caries-free twins</article-title>. <source>Twin Res Hum Genet.</source> (<year>2007</year>) <volume>10</volume>:<fpage>821</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1375/twin.10.6.821</pub-id><pub-id pub-id-type="pmid">18179393</pub-id></citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Demmitt</surname> <given-names>BA</given-names></name> <name><surname>Corley</surname> <given-names>RP</given-names></name> <name><surname>Huibregtse</surname> <given-names>BM</given-names></name> <name><surname>Keller</surname> <given-names>MC</given-names></name> <name><surname>Hewitt</surname> <given-names>JK</given-names></name> <name><surname>McQueen</surname> <given-names>MB</given-names></name> <etal/></person-group>. <article-title>Genetic influences on the human oral microbiome</article-title>. <source>BMC Genomics.</source> (<year>2017</year>) <volume>18</volume>:<fpage>659</fpage>. <pub-id pub-id-type="doi">10.1186/s12864-017-4008-8</pub-id><pub-id pub-id-type="pmid">28836939</pub-id></citation></ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gomez</surname> <given-names>A</given-names></name> <name><surname>Espinoza</surname> <given-names>JL</given-names></name> <name><surname>Harkins</surname> <given-names>DM</given-names></name> <name><surname>Leong</surname> <given-names>P</given-names></name> <name><surname>Saffery</surname> <given-names>R</given-names></name> <name><surname>Bockmann</surname> <given-names>M</given-names></name> <etal/></person-group>. <article-title>Host genetic control of the oral microbiome in health and disease</article-title>. <source>Cell Host Microbe.</source> (<year>2017</year>) <volume>22</volume>:<fpage>269</fpage>&#x02013;<lpage>78</lpage>.e3. <pub-id pub-id-type="doi">10.1016/j.chom.2017.08.013</pub-id><pub-id pub-id-type="pmid">28910633</pub-id></citation></ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hansen</surname> <given-names>TH</given-names></name> <name><surname>Kern</surname> <given-names>T</given-names></name> <name><surname>Bak</surname> <given-names>EG</given-names></name> <name><surname>Kashani</surname> <given-names>A</given-names></name> <name><surname>Allin</surname> <given-names>KH</given-names></name> <name><surname>Nielsen</surname> <given-names>T</given-names></name> <etal/></person-group>. <article-title>Impact of a vegan diet on the human salivary microbiota</article-title>. <source>Sci Rep.</source> (<year>2018</year>) <volume>8</volume>:<fpage>5847</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-24207-3</pub-id><pub-id pub-id-type="pmid">29643500</pub-id></citation></ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kato</surname> <given-names>I</given-names></name> <name><surname>Vasquez</surname> <given-names>A</given-names></name> <name><surname>Moyerbrailean</surname> <given-names>G</given-names></name> <name><surname>Land</surname> <given-names>S</given-names></name> <name><surname>Djuric</surname> <given-names>Z</given-names></name> <name><surname>Sun</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Nutritional correlates of human oral microbiome</article-title>. <source>J Am Coll Nutr.</source> (<year>2017</year>) <volume>36</volume>:<fpage>88</fpage>&#x02013;<lpage>98</lpage>. <pub-id pub-id-type="doi">10.1080/07315724.2016.1185386</pub-id><pub-id pub-id-type="pmid">27797671</pub-id></citation></ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marsh</surname> <given-names>PD</given-names></name></person-group>. <article-title>Microbiology of dental plaque biofilms and their role in oral health and caries</article-title>. <source>Dent Clin North Am.</source> (<year>2010</year>) <volume>54</volume>:<fpage>441</fpage>&#x02013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.1016/j.cden.2010.03.002</pub-id><pub-id pub-id-type="pmid">20630188</pub-id></citation></ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Weyrich</surname> <given-names>LS</given-names></name> <name><surname>Duchene</surname> <given-names>S</given-names></name> <name><surname>Soubrier</surname> <given-names>J</given-names></name> <name><surname>Arriola</surname> <given-names>L</given-names></name> <name><surname>Llamas</surname> <given-names>B</given-names></name> <name><surname>Breen</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus</article-title>. <source>Nature.</source> (<year>2017</year>) <volume>544</volume>:<fpage>357</fpage>&#x02013;<lpage>61</lpage>. <pub-id pub-id-type="doi">10.1038/nature21674</pub-id><pub-id pub-id-type="pmid">28273061</pub-id></citation></ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Handsley-Davis</surname> <given-names>M</given-names></name> <name><surname>Jamieson</surname> <given-names>L</given-names></name> <name><surname>Kapellas</surname> <given-names>K</given-names></name> <name><surname>Hedges</surname> <given-names>J</given-names></name> <name><surname>Weyrich</surname> <given-names>LS</given-names></name></person-group>. <article-title>The role of the oral microbiota in chronic non-communicable disease and its relevance to the Indigenous health gap in Australia</article-title>. <source>BMC Oral Health.</source> (<year>2020</year>) <volume>20</volume>:<fpage>327</fpage>. <pub-id pub-id-type="doi">10.1186/s12903-020-01308-y</pub-id><pub-id pub-id-type="pmid">33198712</pub-id></citation></ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clemente</surname> <given-names>JC</given-names></name> <name><surname>Pehrsson</surname> <given-names>EC</given-names></name> <name><surname>Blaser</surname> <given-names>MJ</given-names></name> <name><surname>Sandhu</surname> <given-names>K</given-names></name> <name><surname>Gao</surname> <given-names>Z</given-names></name> <name><surname>Wang</surname> <given-names>B</given-names></name> <etal/></person-group>. <article-title>The microbiome of uncontacted Amerindians</article-title>. <source>Sci Adv.</source> (<year>2015</year>) <volume>1</volume>:<fpage>e1500183</fpage>. <pub-id pub-id-type="doi">10.1126/sciadv.1500183</pub-id><pub-id pub-id-type="pmid">26229982</pub-id></citation></ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lassalle</surname> <given-names>F</given-names></name> <name><surname>Spagnoletti</surname> <given-names>M</given-names></name> <name><surname>Fumagalli</surname> <given-names>M</given-names></name> <name><surname>Shaw</surname> <given-names>L</given-names></name> <name><surname>Dyble</surname> <given-names>M</given-names></name> <name><surname>Walker</surname> <given-names>C</given-names></name> <etal/></person-group>. <article-title>Oral microbiomes from hunter-gatherers and traditional farmers reveal shifts in commensal balance and pathogen load linked to diet</article-title>. <source>Mol Ecol.</source> (<year>2018</year>) <volume>27</volume>:<fpage>182</fpage>&#x02013;<lpage>95</lpage>. <pub-id pub-id-type="doi">10.1111/mec.14435</pub-id><pub-id pub-id-type="pmid">29165844</pub-id></citation></ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nasidze</surname> <given-names>I</given-names></name> <name><surname>Li</surname> <given-names>J</given-names></name> <name><surname>Schroeder</surname> <given-names>R</given-names></name> <name><surname>Creasey</surname> <given-names>JL</given-names></name> <name><surname>Li</surname> <given-names>M</given-names></name> <name><surname>Stoneking</surname> <given-names>M</given-names></name></person-group>. <article-title>High diversity of the saliva microbiome in batwa pygmies</article-title>. <source>PLoS ONE.</source> (<year>2011</year>) <volume>6</volume>:<fpage>e23352</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0023352</pub-id><pub-id pub-id-type="pmid">21858083</pub-id></citation></ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Handsley-Davis</surname> <given-names>M</given-names></name></person-group>. <article-title><italic>Investigating Oral Microbial Communities in Aboriginal Australians</italic> (BSc</article-title>. Honours thesis). Australia: University of Adelaide (<year>2016</year>).</citation></ref>
<ref id="B30">
<label>30.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lalloo</surname> <given-names>R</given-names></name> <name><surname>Tadakamadla</surname> <given-names>SK</given-names></name> <name><surname>Kroon</surname> <given-names>J</given-names></name> <name><surname>Jamieson</surname> <given-names>LM</given-names></name> <name><surname>Ware</surname> <given-names>RS</given-names></name> <name><surname>Johnson</surname> <given-names>NW</given-names></name></person-group>. <article-title>Carious lesions in permanent dentitions are reduced in remote Indigenous Australian children taking part in a non-randomised preventive trial</article-title>. <source>PLoS ONE</source>. (<year>2021</year>) <volume>16</volume>:<fpage>e0244927</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0244927</pub-id><pub-id pub-id-type="pmid">33507984</pub-id></citation></ref>
<ref id="B31">
<label>31.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gugnani</surname> <given-names>N</given-names></name> <name><surname>Pandit</surname> <given-names>I</given-names></name> <name><surname>Srivastava</surname> <given-names>N</given-names></name> <name><surname>Gupta</surname> <given-names>M</given-names></name> <name><surname>Sharma</surname> <given-names>M</given-names></name></person-group>. <article-title>International Caries Detection and Assessment System (ICDAS): a new concept</article-title>. <source>Int J Clin Pediatr Dent.</source> (<year>2011</year>) <volume>4</volume>:<fpage>93</fpage>&#x02013;<lpage>100</lpage>. <pub-id pub-id-type="doi">10.5005/jp-journals-10005-1089</pub-id><pub-id pub-id-type="pmid">27672245</pub-id></citation></ref>
<ref id="B32">
<label>32.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stevens</surname> <given-names>K</given-names></name></person-group>. <article-title>Developing a descriptive system for a new preference-based measure of health-related quality of life for children</article-title>. <source>Qual Life Res.</source> (<year>2009</year>) <volume>18</volume>:<fpage>1105</fpage>&#x02013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1007/s11136-009-9524-9</pub-id><pub-id pub-id-type="pmid">19693703</pub-id></citation></ref>
<ref id="B33">
<label>33.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Slade</surname> <given-names>GD</given-names></name> <name><surname>Spencer</surname> <given-names>AJ</given-names></name></person-group>. <article-title>Development and evaluation of the oral health impact profile</article-title>. <source>Community Dent Health.</source> (<year>1994</year>) <volume>11</volume>:<fpage>3</fpage>&#x02013;<lpage>11</lpage>.<pub-id pub-id-type="pmid">8193981</pub-id></citation></ref>
<ref id="B34">
<label>34.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname> <given-names>JG</given-names></name> <name><surname>Lauber</surname> <given-names>CL</given-names></name> <name><surname>Walters</surname> <given-names>WA</given-names></name> <name><surname>Berg-Lyons</surname> <given-names>D</given-names></name> <name><surname>Huntley</surname> <given-names>J</given-names></name> <name><surname>Fierer</surname> <given-names>N</given-names></name> <etal/></person-group>. <article-title>Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms</article-title>. <source>ISME J.</source> (<year>2012</year>) <volume>6</volume>:<fpage>1621</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2012.8</pub-id><pub-id pub-id-type="pmid">22402401</pub-id></citation></ref>
<ref id="B35">
<label>35.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bolyen</surname> <given-names>E</given-names></name> <name><surname>Rideout</surname> <given-names>JR</given-names></name> <name><surname>Dillon</surname> <given-names>MR</given-names></name> <name><surname>Bokulich</surname> <given-names>NA</given-names></name> <name><surname>Abnet</surname> <given-names>CC</given-names></name> <name><surname>Al-Ghalith</surname> <given-names>GA</given-names></name> <etal/></person-group>. <article-title>Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2</article-title>. <source>Nat Biotechnol.</source> (<year>2019</year>) <volume>37</volume>:<fpage>852</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-019-0209-9</pub-id><pub-id pub-id-type="pmid">31399723</pub-id></citation></ref>
<ref id="B36">
<label>36.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amir</surname> <given-names>A</given-names></name> <name><surname>McDonald</surname> <given-names>D</given-names></name> <name><surname>Navas-Molina</surname> <given-names>JA</given-names></name> <name><surname>Kopylova</surname> <given-names>E</given-names></name> <name><surname>Morton</surname> <given-names>JT</given-names></name> <name><surname>Xu</surname> <given-names>ZZ</given-names></name> <etal/></person-group>. <article-title>Deblur rapidly resolves single-nucleotide community sequence patterns</article-title>. <source>mSystems.</source> (<year>2017</year>) <volume>2</volume>:<fpage>e00191</fpage>&#x02013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1128/mSystems.00191-16</pub-id><pub-id pub-id-type="pmid">28289731</pub-id></citation></ref>
<ref id="B37">
<label>37.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bisanz</surname> <given-names>JE</given-names></name></person-group>.<source>qiime2R: Importing QIIME2 Artifacts and Associated Data into R Sessions R</source> <publisher-name>GitHub</publisher-name> (<year>2018</year>).</citation></ref>
<ref id="B38">
<label>38.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davis</surname> <given-names>NM</given-names></name> <name><surname>Proctor</surname> <given-names>DM</given-names></name> <name><surname>Holmes</surname> <given-names>SP</given-names></name> <name><surname>Relman</surname> <given-names>DA</given-names></name> <name><surname>Callahan</surname> <given-names>BJ</given-names></name></person-group>. <article-title>Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data</article-title>. <source>Microbiome.</source> (<year>2018</year>) <volume>6</volume>:<fpage>226</fpage>. <pub-id pub-id-type="doi">10.1186/s40168-018-0605-2</pub-id><pub-id pub-id-type="pmid">30558668</pub-id></citation></ref>
<ref id="B39">
<label>39.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Faith</surname> <given-names>DP</given-names></name> <name><surname>Baker</surname> <given-names>AM</given-names></name></person-group>. <article-title>Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges</article-title>. <source>Evol Bioinforma Online.</source> (<year>2007</year>) <volume>2</volume>:<fpage>121</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1177/117693430600200007</pub-id><pub-id pub-id-type="pmid">19455206</pub-id></citation></ref>
<ref id="B40">
<label>40.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lozupone</surname> <given-names>C</given-names></name> <name><surname>Knight</surname> <given-names>R</given-names></name></person-group>. <article-title>UniFrac: a new phylogenetic method for comparing microbial communities</article-title>. <source>Appl Environ Microbiol.</source> (<year>2005</year>) <volume>71</volume>:<fpage>8228</fpage>&#x02013;<lpage>35</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.71.12.8228-8235.2005</pub-id><pub-id pub-id-type="pmid">16332807</pub-id></citation></ref>
<ref id="B41">
<label>41.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kruskal</surname> <given-names>WH</given-names></name> <name><surname>Wallis</surname> <given-names>WA</given-names></name></person-group>. <article-title>Use of ranks in one-criterion variance analysis</article-title>. <source>J Am Stat Assoc.</source> (<year>1952</year>) <volume>47</volume>:<fpage>583</fpage>&#x02013;<lpage>621</lpage>. <pub-id pub-id-type="doi">10.1080/01621459.1952.10483441</pub-id></citation></ref>
<ref id="B42">
<label>42.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>MJ</given-names></name></person-group>. <article-title>A new method for non-parametric multivariate analysis of variance</article-title>. <source>Austral Ecol.</source> (<year>2001</year>) <volume>26</volume>:<fpage>32</fpage>&#x02013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.1046/j.1442-9993.2001.01070.x</pub-id></citation></ref>
<ref id="B43">
<label>43.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Oksanen</surname> <given-names>J</given-names></name> <name><surname>Blanchet</surname> <given-names>FG</given-names></name> <name><surname>Friendly</surname> <given-names>M</given-names></name> <name><surname>Kindt</surname> <given-names>R</given-names></name> <name><surname>Legendre</surname> <given-names>P</given-names></name> <name><surname>McGlinn</surname> <given-names>D</given-names></name> <etal/></person-group>. <source>vegan: Community Ecology Package</source> <publisher-name>CRAN</publisher-name> (<year>2018</year>).</citation></ref>
<ref id="B44">
<label>44.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>MJ</given-names></name></person-group>. <article-title>Distance-based tests for homogeneity of multivariate dispersions</article-title>. <source>Biometrics.</source> (<year>2006</year>) <volume>62</volume>:<fpage>245</fpage>&#x02013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.1111/j.1541-0420.2005.00440.x</pub-id><pub-id pub-id-type="pmid">16542252</pub-id></citation></ref>
<ref id="B45">
<label>45.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mandal</surname> <given-names>S</given-names></name> <name><surname>Van</surname> <given-names>Treuren W</given-names></name> <name><surname>White</surname> <given-names>RA</given-names></name> <name><surname>Eggesb&#x000F8;</surname> <given-names>M</given-names></name> <name><surname>Knight</surname> <given-names>R</given-names></name> <name><surname>Peddada</surname> <given-names>SD</given-names></name></person-group>. <article-title>Analysis of composition of microbiomes: a novel method for studying microbial composition</article-title>. <source>Microb Ecol Health Dis.</source> (<year>2015</year>) <volume>26</volume>:<fpage>27663</fpage>. <pub-id pub-id-type="doi">10.3402/mehd.v26.27663</pub-id><pub-id pub-id-type="pmid">26028277</pub-id></citation></ref>
<ref id="B46">
<label>46.</label>
<citation citation-type="book"><person-group person-group-type="author"><collab>RStudio Team</collab></person-group>. <source>RStudio: Integrated development for R</source>. <publisher-loc>Boston, MA</publisher-loc>: <publisher-name>RStudio, Inc</publisher-name>. (<year>2015</year>).</citation></ref>
<ref id="B47">
<label>47.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Wickham</surname> <given-names>H</given-names></name></person-group>. <source>ggplot2: Elegant Graphics for Data Analysis</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer-Verlag</publisher-name> (<year>2009</year>).</citation></ref>
<ref id="B48">
<label>48.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eisenhofer</surname> <given-names>R</given-names></name> <name><surname>Minich</surname> <given-names>JJ</given-names></name> <name><surname>Marotz</surname> <given-names>C</given-names></name> <name><surname>Cooper</surname> <given-names>A</given-names></name> <name><surname>Knight</surname> <given-names>R</given-names></name> <name><surname>Weyrich</surname> <given-names>LS</given-names></name></person-group>. <article-title>Contamination in low microbial biomass microbiome studies: issues and recommendations</article-title>. <source>Trends Microbiol.</source> (<year>2019</year>) <volume>27</volume>:<fpage>105</fpage>&#x02013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1016/j.tim.2018.11.003</pub-id><pub-id pub-id-type="pmid">30497919</pub-id></citation></ref>
<ref id="B49">
<label>49.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Weyrich</surname> <given-names>LS</given-names></name> <name><surname>Farrer</surname> <given-names>AG</given-names></name> <name><surname>Eisenhofer</surname> <given-names>R</given-names></name> <name><surname>Arriola</surname> <given-names>LA</given-names></name> <name><surname>Young</surname> <given-names>J</given-names></name> <name><surname>Selway</surname> <given-names>CA</given-names></name> <etal/></person-group>. <article-title>Laboratory contamination over time during low-biomass sample analysis</article-title>. <source>Mol Ecol Resour.</source> (<year>2019</year>) <volume>19</volume>:<fpage>982</fpage>&#x02013;<lpage>96</lpage>. <pub-id pub-id-type="doi">10.1111/1755-0998.13011</pub-id><pub-id pub-id-type="pmid">30887686</pub-id></citation></ref>
<ref id="B50">
<label>50.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Skelly</surname> <given-names>E</given-names></name> <name><surname>Johnson</surname> <given-names>NW</given-names></name> <name><surname>Kapellas</surname> <given-names>K</given-names></name> <name><surname>Kroon</surname> <given-names>J</given-names></name> <name><surname>Lalloo</surname> <given-names>R</given-names></name> <name><surname>Weyrich</surname> <given-names>L</given-names></name></person-group>. <article-title>Response of salivary microbiota to caries preventive treatment in aboriginal and torres strait islander children</article-title>. <source>J Oral Microbiol.</source> (<year>2020</year>) <volume>12</volume>:<fpage>1830623</fpage>. <pub-id pub-id-type="doi">10.1080/20002297.2020.1830623</pub-id><pub-id pub-id-type="pmid">33149844</pub-id></citation></ref>
<ref id="B51">
<label>51.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eren</surname> <given-names>AM</given-names></name> <name><surname>Borisy</surname> <given-names>GG</given-names></name> <name><surname>Huse</surname> <given-names>SM</given-names></name> <name><surname>Welch</surname> <given-names>JLM</given-names></name></person-group>. <article-title>Oligotyping analysis of the human oral microbiome</article-title>. <source>Proc Natl Acad Sci USA.</source> (<year>2014</year>) <volume>111</volume>:<fpage>E2875</fpage>&#x02013;<lpage>84</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1409644111</pub-id><pub-id pub-id-type="pmid">24965363</pub-id></citation></ref>
<ref id="B52">
<label>52.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Okada</surname> <given-names>M</given-names></name> <name><surname>Kawamura</surname> <given-names>M</given-names></name> <name><surname>Oda</surname> <given-names>Y</given-names></name> <name><surname>Yasuda</surname> <given-names>R</given-names></name> <name><surname>Kojima</surname> <given-names>T</given-names></name> <name><surname>Kurihara</surname> <given-names>H</given-names></name></person-group>. <article-title>Caries prevalence associated with <italic>Streptococcus mutans</italic> and <italic>Streptococcus sobrinus</italic> in Japanese schoolchildren</article-title>. <source>Int J Paediatr Dent.</source> (<year>2012</year>) <volume>22</volume>:<fpage>342</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-263X.2011.01203.x</pub-id><pub-id pub-id-type="pmid">22225789</pub-id></citation></ref>
<ref id="B53">
<label>53.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Freire</surname> <given-names>M</given-names></name> <name><surname>Moustafa</surname> <given-names>A</given-names></name> <name><surname>Harkins</surname> <given-names>DM</given-names></name> <name><surname>Torralba</surname> <given-names>MG</given-names></name> <name><surname>Zhang</surname> <given-names>Y</given-names></name> <name><surname>Leong</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>Longitudinal study of oral microbiome variation in twins</article-title>. 1. <source>Sci Rep</source>. (<year>2020</year>) <volume>10</volume>:<fpage>7954</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-020-64747-1</pub-id><pub-id pub-id-type="pmid">32409670</pub-id></citation></ref>
<ref id="B54">
<label>54.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Richards</surname> <given-names>VP</given-names></name> <name><surname>Alvarez</surname> <given-names>AJ</given-names></name> <name><surname>Luce</surname> <given-names>AR</given-names></name> <name><surname>Bedenbaugh</surname> <given-names>M</given-names></name> <name><surname>Mitchell</surname> <given-names>ML</given-names></name> <name><surname>Burne</surname> <given-names>RA</given-names></name> <etal/></person-group>. <article-title>Microbiomes of site-specific dental plaques from children with different caries status</article-title>. <source>Infect Immun.</source> (<year>2017</year>) <volume>85</volume>:<fpage>e00106</fpage>&#x02013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1128/IAI.00106-17</pub-id><pub-id pub-id-type="pmid">28507066</pub-id></citation></ref>
<ref id="B55">
<label>55.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Teng</surname> <given-names>F</given-names></name> <name><surname>Yang</surname> <given-names>F</given-names></name> <name><surname>Huang</surname> <given-names>S</given-names></name> <name><surname>Bo</surname> <given-names>C</given-names></name> <name><surname>Xu</surname> <given-names>ZZ</given-names></name> <name><surname>Amir</surname> <given-names>A</given-names></name> <etal/></person-group>. <article-title>Prediction of early childhood caries via spatial-temporal variations of oral microbiota</article-title>. <source>Cell Host Microbe.</source> (<year>2015</year>) <volume>18</volume>:<fpage>296</fpage>&#x02013;<lpage>306</lpage>. <pub-id pub-id-type="doi">10.1016/j.chom.2015.08.005</pub-id><pub-id pub-id-type="pmid">26355216</pub-id></citation></ref>
<ref id="B56">
<label>56.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Belstr&#x000F8;m</surname> <given-names>D</given-names></name> <name><surname>Holmstrup</surname> <given-names>P</given-names></name> <name><surname>Fiehn</surname> <given-names>N-E</given-names></name> <name><surname>Kirkby</surname> <given-names>N</given-names></name> <name><surname>Kokaras</surname> <given-names>A</given-names></name> <name><surname>Paster</surname> <given-names>BJ</given-names></name> <etal/></person-group>. <article-title>Salivary microbiota in individuals with different levels of caries experience</article-title>. <source>J Oral Microbiol.</source> (<year>2017</year>) <volume>9</volume>:<fpage>1270614</fpage>. <pub-id pub-id-type="doi">10.1080/20002297.2016.1270614</pub-id><pub-id pub-id-type="pmid">28326153</pub-id></citation></ref>
<ref id="B57">
<label>57.</label>
<citation citation-type="journal"><person-group person-group-type="author"><collab>Human Oral Microbiome Database (HOMD)</collab></person-group>. <source>Lactobacillus gasseri Human Oral Microbiome Taxon Description</source>.</citation></ref>
<ref id="B58">
<label>58.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mashima</surname> <given-names>I</given-names></name> <name><surname>Nakazawa</surname> <given-names>F</given-names></name></person-group>. <article-title>Interaction between Streptococcus spp. and <italic>Veillonella tobetsuensis</italic> in the early stages of oral biofilm formation</article-title>. <source>J Bacteriol.</source> (<year>2015</year>) <volume>197</volume>:<fpage>2104</fpage>&#x02013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1128/JB.02512-14</pub-id><pub-id pub-id-type="pmid">25917902</pub-id></citation></ref>
<ref id="B59">
<label>59.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>W</given-names></name> <name><surname>Jiang</surname> <given-names>Q</given-names></name> <name><surname>Yan</surname> <given-names>G</given-names></name> <name><surname>Yang</surname> <given-names>D</given-names></name></person-group>. <article-title>The oral microbiome and salivary proteins influence caries in children aged 6 to 8 years</article-title>. <source>BMC Oral Health.</source> (<year>2020</year>) <volume>20</volume>:<fpage>295</fpage>. <pub-id pub-id-type="doi">10.1186/s12903-020-01262-9</pub-id><pub-id pub-id-type="pmid">33115458</pub-id></citation></ref>
<ref id="B60">
<label>60.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Segata</surname> <given-names>N</given-names></name> <name><surname>Haake</surname> <given-names>SK</given-names></name> <name><surname>Mannon</surname> <given-names>P</given-names></name> <name><surname>Lemon</surname> <given-names>KP</given-names></name> <name><surname>Waldron</surname> <given-names>L</given-names></name> <name><surname>Gevers</surname> <given-names>D</given-names></name> <etal/></person-group>. <article-title>Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples</article-title>. <source>Genome Biol.</source> (<year>2012</year>) <volume>13</volume>:<fpage>R42</fpage>. <pub-id pub-id-type="doi">10.1186/gb-2012-13-6-r42</pub-id><pub-id pub-id-type="pmid">22698087</pub-id></citation></ref>
<ref id="B61">
<label>61.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dawes</surname> <given-names>C</given-names></name> <name><surname>Tsang</surname> <given-names>RWL</given-names></name> <name><surname>Suelzle</surname> <given-names>T</given-names></name></person-group>. <article-title>The effects of gum chewing, four oral hygiene procedures, and two saliva collection techniques, on the output of bacteria into human whole saliva</article-title>. <source>Arch Oral Biol.</source> (<year>2001</year>) <volume>46</volume>:<fpage>625</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1016/S0003-9969(01)00017-6</pub-id><pub-id pub-id-type="pmid">11369317</pub-id></citation></ref>
<ref id="B62">
<label>62.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Baca</surname> <given-names>P</given-names></name> <name><surname>Castillo</surname> <given-names>AM</given-names></name> <name><surname>Baca</surname> <given-names>AP</given-names></name> <name><surname>Li&#x000E9;bana</surname> <given-names>MJ</given-names></name> <name><surname>Junco</surname> <given-names>P</given-names></name> <name><surname>Li&#x000E9;bana</surname> <given-names>J</given-names></name></person-group>. <article-title>Genotypes of <italic>Streptococcus mutans</italic> in saliva versus dental plaque</article-title>. <source>Arch Oral Biol.</source> (<year>2008</year>) <volume>53</volume>:<fpage>751</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1016/j.archoralbio.2008.02.011</pub-id><pub-id pub-id-type="pmid">18374899</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> Field work for this project was conducted under grant APP from the Australian National Health and Medical Research Council APP1081320 (Principal Investigators NJ, JK, and RL). Additional funding support for laboratory and sequencing costs was provided by the Australian Research Council (ARC) through an ARC Future Fellowship (FT180100407) and an ARC Centre of Excellence (CE170100015) grant (LW). MH-D was supported by a Westpac Future Leaders Scholarship and an Australian Government Research Training Program scholarship.</p></fn>
</fn-group>
</back>
</article>