<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<?covid-19-tdm?>
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Immunol.</journal-id>
<journal-title>Frontiers in Immunology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Immunol.</abbrev-journal-title>
<issn pub-type="epub">1664-3224</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fimmu.2022.1079995</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Immunology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Persistence of salivary antibody responses after COVID-19 vaccination is associated with oral microbiome variation in both healthy and people living with HIV</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Ghorbani</surname>
<given-names>Mahin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2074804"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Al-Manei</surname>
<given-names>Khaled</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="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Naud</surname>
<given-names>Sabrina</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Healy</surname>
<given-names>Katie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1665018"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gabarrini</surname>
<given-names>Giorgio</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1056282"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sobkowiak</surname>
<given-names>Michal Jacek</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/584234"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Puran</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ray</surname>
<given-names>Shilpa</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1319547"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Akber</surname>
<given-names>Mira</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Muschiol</surname>
<given-names>Sandra</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bogdanovic</surname>
<given-names>Gordana</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bergman</surname>
<given-names>Peter</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff6">
<sup>6</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/81668"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ljungman</surname>
<given-names>Per</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff7">
<sup>7</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1180575"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Buggert</surname>
<given-names>Marcus</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/188979"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ljunggren</surname>
<given-names>Hans-Gustaf</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/56029"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pin</surname>
<given-names>Elisa</given-names>
</name>
<xref ref-type="aff" rid="aff8">
<sup>8</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Nowak</surname>
<given-names>Piotr</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff9">
<sup>9</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1861983"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Aleman</surname>
<given-names>Soo</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff9">
<sup>9</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>S&#xe4;llberg Chen</surname>
<given-names>Margaret</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">*</xref>
<uri xlink:href="https://loop.frontiersin.org/people/526795"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Dental Medicine, Karolinska Institutet</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Restorative Dental Science, College of Dentistry, King Saud University</institution>, <addr-line>Riyadh</addr-line>, <country>Saudi Arabia</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Department of Laboratory Medicine, Karolinska Institutet</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Department of Medicine, Huddinge, Karolinska Institutet</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Department of Clinical Microbiology, Karolinska University Hospital</institution>, <addr-line>Stockholm</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff6">
<sup>6</sup>
<institution>Department of Clinical Immunology and Transfusion Medicine, Karolinska University Hospital</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff7">
<sup>7</sup>
<institution>Department of Cellular Therapy and Allogeneic Stem Cell Transplantation, Karolinska Comprehensive Cancer Center, Karolinska University Hospital</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff8">
<sup>8</sup>
<institution>Department of Protein Science, SciLifeLab, KTH Royal Institute of Technology</institution>, <addr-line>Stockholm</addr-line>, <country>Sweden</country>
</aff>
<aff id="aff9">
<sup>9</sup>
<institution>Department of Infectious Diseases, Karolinska University Hospital</institution>, <addr-line>Huddinge</addr-line>, <country>Sweden</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Mahbuba Rahman, McMaster University, Canada</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Rupsa Basu, Humane Genomics, United States; Afsheen Raza, Hamad Medical Corporation, Qatar</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Margaret S&#xe4;llberg Chen, <email xlink:href="mailto:Margaret.Chen@ki.se">Margaret.Chen@ki.se</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Vaccines and Molecular Therapeutics, a section of the journal Frontiers in Immunology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>10</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>13</volume>
<elocation-id>1079995</elocation-id>
<history>
<date date-type="received">
<day>25</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>28</day>
<month>11</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Ghorbani, Al-Manei, Naud, Healy, Gabarrini, Sobkowiak, Chen, Ray, Akber, Muschiol, Bogdanovic, Bergman, Ljungman, Buggert, Ljunggren, Pin, Nowak, Aleman and S&#xe4;llberg Chen</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Ghorbani, Al-Manei, Naud, Healy, Gabarrini, Sobkowiak, Chen, Ray, Akber, Muschiol, Bogdanovic, Bergman, Ljungman, Buggert, Ljunggren, Pin, Nowak, Aleman and S&#xe4;llberg Chen</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>Coevolution of microbiome and immunity at mucosal sites is essential for our health. Whether the oral microbiome, the second largest community after the gut, contributes to the immunogenicity of COVID-19 vaccines is not known. We investigated the baseline oral microbiome in individuals in the COVAXID clinical trial receiving the BNT162b2 mRNA vaccine. Participants (n=115) included healthy controls (HC; n=57) and people living with HIV (PLHIV; n=58) who met the study selection criteria. Vaccine-induced Spike antibodies in saliva and serum from 0 to 6 months were assessed and comparative analyses were performed against the individual salivary 16S ASV microbiome diversity. High- versus low vaccine responders were assessed on general, immunological, and oral microbiome features. Our analyses identified oral microbiome features enriched in high- <italic>vs</italic>. low-responders among healthy and PLHIV participants. In low-responders, an enrichment of Gram-negative, anaerobic species with proteolytic activity were found including <italic>Campylobacter, Butyrivibrio, Selenomonas, Lachnoanaerobaculum, Leptotrichia, Megasphaera, Prevotella</italic> and <italic>Stomatobaculum</italic>. In high-responders, enriched species were mainly Gram-positive and saccharolytic facultative anaerobes: <italic>Abiotrophia, Corynebacterium, Gemella, Granulicatella, Rothia</italic>, and <italic>Haemophilus</italic>. Combining identified microbial features in a classifier using the area under the receiver operating characteristic curve (ROC AUC) yielded scores of 0.879 (healthy controls) to 0.82 (PLHIV), supporting the oral microbiome contribution in the long-term vaccination outcome. The present study is the first to suggest that the oral microbiome has an impact on the durability of mucosal immunity after Covid-19 vaccination. Microbiome-targeted interventions to enhance long-term duration of mucosal vaccine immunity may be exploited.</p>
</abstract>
<kwd-group>
<kwd>Oral microbiome</kwd>
<kwd>mRNA vaccination</kwd>
<kwd>saliva</kwd>
<kwd>SARS-CoV-2</kwd>
<kwd>campylobacter</kwd>
<kwd>leptotrichia</kwd>
<kwd>gemella</kwd>
<kwd>granulicatella</kwd>
</kwd-group>
<contract-sponsor id="cn001">Knut och Alice Wallenbergs Stiftelse<named-content content-type="fundref-id">10.13039/501100004063</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Center for Innovative Medicine<named-content content-type="fundref-id">10.13039/501100018713</named-content>
</contract-sponsor>
<counts>
<fig-count count="4"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="30"/>
<page-count count="11"/>
<word-count count="5039"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>SARS-CoV-2 infection and replication take place in the oral cavity and saliva (<xref ref-type="bibr" rid="B1">1</xref>). The presence of local immune control is important in limiting the viral infection and transmission. A durable mucosal immunity at this site is therefore a highly desired outcome in a COVID-vaccination. One of the most immunogenic COVID-19 vaccines - the BNT162b2 (Comirnaty<sup>&#xae;</sup>) mRNA vaccine induces high titres of systemic SARS-CoV-2 spike-specific antibodies (<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B3">3</xref>). In the oral cavity and saliva, it conveys virus-binding antibodies highly correlative to neutralizing capacity. The magnitude of this acquired local immunity however varies among both healthy and immunocompromised vaccinees (<xref ref-type="bibr" rid="B4">4</xref>), and the durability of those translocated specific antibodies at the oral mucosa is presently unknown. Recent developments in techniques for microbiome sequencing have enabled a comprehensive analysis of commensal microbiota, with over 700 bacterial species detected in the oral cavity (<xref ref-type="bibr" rid="B5">5</xref>). The variants of microorganisms form unique communities, and the oral microbiome composition is believed to be resilient and remain rather stable within individuals longitudinally (<xref ref-type="bibr" rid="B6">6</xref>). The normal oral microbiota is continuously in contact with the oral mucosa and plays a role in modulating immune surveillance mechanisms (<xref ref-type="bibr" rid="B7">7</xref>). Dysbiotic oral microbiota on the other hand is associated with oral diseases, medical conditions, dietary habits, and lifestyles (<xref ref-type="bibr" rid="B8">8</xref>&#x2013;<xref ref-type="bibr" rid="B10">10</xref>). Recently, a negative link between vaccine efficacy and gut microbiota variation linked to BMI and systemic inflammation was reported (<xref ref-type="bibr" rid="B11">11</xref>), other studies have reported that decreased seroconversion in BNT162b2 mRNA recipients is associated with antibiotic use (<xref ref-type="bibr" rid="B12">12</xref>). Although the gut microbiota is recognized as a key factor in supporting gut homeostasis and health (<xref ref-type="bibr" rid="B13">13</xref>), the contribution of oral microbiota to COVID vaccination is still unknown.</p>
<p>In the present study, we hypothesize that oral microbiota composition has a role in the maintenance of the salivary immunity induced by mRNA vaccination. We investigated the salivary microbiota using 16S rRNA sequencing analysis in samples obtained from healthy participants (HC) and people living with HIV (PLHIV). The aim was to characterize salivary microbiota signatures associated with a durable vaccine response capable of persisting at the oral mucosal site.</p>
</sec>
<sec id="s2" sec-type="results">
<title>Results</title>
<sec id="s2_1">
<title>Study design and participants</title>
<p>Eligible participants were healthy controls (HC; n=57) and people living with HIV (PLHIV; n=58) who had received two vaccine doses and fulfilled the study inclusion criteria including baseline and 6-month screenings for negativity for SARS-CoV-2 exposures. As described earlier, all participants seroconverted during the one-month follow-up (day 35) i.e.,14 days after dose 2 in both groups (<xref ref-type="bibr" rid="B4">4</xref>). Vaccine-induced spike IgG in serum and saliva among the participants correlated significantly during the 6-month follow-up, although a proportion showed a larger reduction in salivary spike IgG in than others (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S1</bold>
</xref>). Therefore, we sub-grouped the participants as High- or Low-responders, respectively, based on the expected convalescence spike IgG level at this time-point (<xref ref-type="bibr" rid="B14">14</xref>). The subgroup characteristics and vaccine response data are summarized in <xref ref-type="table" rid="T1">
<bold>Tables&#xa0;1</bold>
</xref>, <xref ref-type="table" rid="T2">
<bold>2</bold>
</xref>. As shown, Low-responders in both the HC and the PLHIV cohorts demonstrated a reduced ability to maintain salivary IgG responses to the full-length trimeric spike (S.f) as well as the S1 spike antigen during this 6-month follow-up. In PLHIV Low-responders, reduced magnitude and duration of spike-IgG was also found in both saliva and serum through the entire follow-up period (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). However, High- and Low responders showed no significant differences in any general or medical variable such as age, gender, diet, BMI, ethnicity, anti-inflammatory medications, or clinical immunological variables (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). The peak and long-term spike IgG levels in saliva were both about 6-fold lower in Low-responders (n=29) than in the High-responders (n=28). Similarly, in PLHIV Low-responders (n=38), a significant reduction of anti-spike IgG level both in serum and saliva was noted compared to High-responders (n=20). A detailed assessment on the anti-capsid serology and PCR-reports further confirmed that no break-through infection had occurred during the follow-up that could have influenced these results (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>).</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Demographics and medical report of study participants.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Variables</th>
<th valign="top" colspan="2" align="center">HC (n = 57)</th>
<th valign="top" align="center">P-value</th>
<th valign="top" colspan="2" align="center">PLHIV (n = 58)</th>
<th valign="top" align="center">P-value</th>
</tr>
<tr>
<th valign="top" align="left"/>
<th valign="top" align="center">High responder(n = 28)</th>
<th valign="top" align="center">Low responder(n = 29)</th>
<th valign="top" align="center"/>
<th valign="top" align="center">High responder(n = 20)</th>
<th valign="top" align="center">Low responder(n = 38)</th>
<th valign="top" align="center"/>
</tr>
</thead>
<tbody>
<tr>
<th valign="top" colspan="7" align="left">Demographic variables</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>Age <sup>A</sup> (years)</bold>
</td>
<td valign="top" align="center">47 (32.3)</td>
<td valign="top" align="center">58 (30)</td>
<td valign="top" align="center">0.247</td>
<td valign="top" align="center">54 (22.3)</td>
<td valign="top" align="center">54 (18.3)</td>
<td valign="top" align="center">0.518</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Gender</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.881</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.258</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Male <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">12 (42.9)</td>
<td valign="top" align="center">13 (44.8)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">9 (45)</td>
<td valign="top" align="center">23 (60.5)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Female <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">16 (57.1)</td>
<td valign="top" align="center">16 (55.2)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">11 (55)</td>
<td valign="top" align="center">15 (39.5)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Age-group</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.190</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.327</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>&#x2265; 60 <sup>B</sup> (years)</bold>
</td>
<td valign="top" align="center">7 (25)</td>
<td valign="top" align="center">12 (41.4)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">8 (40)</td>
<td valign="top" align="center">14 (36.8)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>40-59 <sup>B</sup> (years)</bold>
</td>
<td valign="top" align="center">8 (28.6)</td>
<td valign="top" align="center">10 (34.5)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">8 (40)</td>
<td valign="top" align="center">21 (55.3)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>18-39 <sup>B</sup> (years)</bold>
</td>
<td valign="top" align="center">13 (46.4)</td>
<td valign="top" align="center">7 (24.1)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">4 (20)</td>
<td valign="top" align="center">3 (7.9)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>BMI <sup>A</sup> (kg/m<sup>2</sup>)</bold>
</td>
<td valign="top" align="center">25 (7)</td>
<td valign="top" align="center">25.5 (5)</td>
<td valign="top" align="center">0.600</td>
<td valign="top" align="center">25.1 (4.3)</td>
<td valign="top" align="center">25.2 (5.7)</td>
<td valign="top" align="center">0.390</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Race</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.721</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Caucasian <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">9 (45)</td>
<td valign="top" align="center">23 (60.5)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Black/African American <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">5 (25)</td>
<td valign="top" align="center">7 (18.4)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Latin <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">5 (25)</td>
<td valign="top" align="center">7 (18.4)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Unknown <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">1 (5)</td>
<td valign="top" align="center">1 (2.7)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Diet</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Normal <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38 (100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Vegetarian <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<th valign="top" colspan="7" align="left">Medical variables</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>Non-steroidal anti-inflammatory drugs</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Yes <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>No <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38 (100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>IgG baseline <sup>A</sup>(g/L)</bold>
</td>
<td valign="top" align="center">11 (3.3)</td>
<td valign="top" align="center">11 (2)</td>
<td valign="top" align="center">0.293</td>
<td valign="top" align="center">13.5 (3)</td>
<td valign="top" align="center">12.5 (5)</td>
<td valign="top" align="center">0.310</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Lymphocytes baseline <sup>A</sup> (x10<sup>9</sup>/L)</bold>
</td>
<td valign="top" align="center">1.8 (1)</td>
<td valign="top" align="center">1.7 (1)</td>
<td valign="top" align="center">0.712</td>
<td valign="top" align="center">1.8 (1.1)</td>
<td valign="top" align="center">1.6 (2.8)</td>
<td valign="top" align="center">0.969</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Creatinine baseline <sup>A</sup> (&#xb5;mol/L)</bold>
</td>
<td valign="top" align="center">69 (2.7)</td>
<td valign="top" align="center">71 (22.5)</td>
<td valign="top" align="center">0.527</td>
<td valign="top" align="center">77 (35.5)</td>
<td valign="top" align="center">83 (32.5)</td>
<td valign="top" align="center">0.423</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Mann&#x2013;Whitney U test was applied to compare the continuous variables and X<sup>2</sup> test for analyzing categorical variables using SPSS version 24.0 (IBM Crop, Armonk, NY, USA) software. <bold>
<sup>A</sup>
</bold> variable is illustrated as the median (IQR), whereas <bold>
<sup>B</sup>
</bold> variable is illustrated as the number (%).  IQR, interquartile range.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Virological data, serum and saliva antibody status of study participants.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Variables</th>
<th valign="top" colspan="2" align="center">HC (n = 57)</th>
<th valign="top" align="center">P-value</th>
<th valign="top" colspan="2" align="center">PLHIV (n = 58)</th>
<th valign="top" align="center">P-value</th>
</tr>
<tr>
<th valign="top" align="left"/>
<th valign="top" align="center">High responder(n = 28)</th>
<th valign="top" align="center">Low responder(n = 29)</th>
<th valign="top" align="center"/>
<th valign="top" align="center">High responder(n = 20)</th>
<th valign="top" align="center">Low responder(n = 38)</th>
<th valign="top" align="center"/>
</tr>
</thead>
<tbody>
<tr>
<th valign="top" colspan="7" align="left">Virological data</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>PCR-positive at baseline</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">Yes <sup>B</sup>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">No <sup>B</sup>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38 (100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Antibody-positive at baseline</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">Yes <sup>B</sup>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">No <sup>B</sup>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38 (100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>CD4 count</bold>
</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.765</td>
</tr>
<tr>
<td valign="top" align="left">&#x2264; 300 <sup>B</sup> (cells/mm<sup>3</sup>)</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center"/>
<td valign="top" align="center">6 (30)</td>
<td valign="top" align="center">10 (26.3)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&gt; 300 <sup>B</sup> (cells/mm<sup>3</sup>)</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center">NA</td>
<td valign="top" align="center"/>
<td valign="top" align="center">14 (70)</td>
<td valign="top" align="center">28 (73.7)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<th valign="top" colspan="7" align="left">Serum antibody level (AU)</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S d0 <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">0.4 (0)</td>
<td valign="top" align="center">0.4 (0)</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.4 (0)</td>
<td valign="top" align="center">0.4 (0)</td>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S d10 <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">0.8 (1.4)</td>
<td valign="top" align="center">0.4 (0.9)</td>
<td valign="top" align="center">0.090</td>
<td valign="top" align="center">0.4 (1.6)</td>
<td valign="top" align="center">0.9 (1.7)</td>
<td valign="top" align="center">0.395</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S d21 <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">57.4 (119.9)</td>
<td valign="top" align="center">62.4 (149)</td>
<td valign="top" align="center">0.786</td>
<td valign="top" align="center">43.4 (81.6)</td>
<td valign="top" align="center">19.6 (39.7)</td>
<td valign="top" align="center">
<bold>0.015*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S d35 <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">2368.5 (1189)</td>
<td valign="top" align="center">1625 (2878.5)</td>
<td valign="top" align="center">0.231</td>
<td valign="top" align="center">1972 (1526)</td>
<td valign="top" align="center">1049 (1338)</td>
<td valign="top" align="center">
<bold>0.005*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S mo.6 <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">739.5 (691.3)</td>
<td valign="top" align="center">559 (737.5)</td>
<td valign="top" align="center">0.078</td>
<td valign="top" align="center">782 (967)</td>
<td valign="top" align="center">249 (412.5)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Positive seroconversion</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.999</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>d.35</bold> Yes <sup>B</sup>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">29 (100)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38 (100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>mo 6.</bold>Yes <sup>B</sup>
</td>
<td valign="top" align="center">28 (100)</td>
<td valign="top" align="center">28 (96.6)</td>
<td valign="top" align="center"/>
<td valign="top" align="center">20 (100)</td>
<td valign="top" align="center">38(100)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<th valign="top" colspan="7" align="left">Saliva antibody level (MFI)</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.0</bold>
<break/>
<bold>Over cutoff (55 MFI)</bold> <sup>B</sup>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.10 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">78 (306)</td>
<td valign="top" align="center">30 (32)</td>
<td valign="top" align="center">
<bold>0.003*</bold>
</td>
<td valign="top" align="center">78.5 (34.8)</td>
<td valign="top" align="center">41 (65.3)</td>
<td valign="top" align="center">0.098</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.21 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">540 (654)</td>
<td valign="top" align="center">203 (415)</td>
<td valign="top" align="center">
<bold>0.028*</bold>
</td>
<td valign="top" align="center">340 (438.8)</td>
<td valign="top" align="center">159 (315.5)</td>
<td valign="top" align="center">
<bold>0.012*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.35 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">3248 (3287)</td>
<td valign="top" align="center">571 (1994.5)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
<td valign="top" align="center">3875 (6557)</td>
<td valign="top" align="center">1251 (2166.5)</td>
<td valign="top" align="center">
<bold>0.001*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f mo.6 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">438.00 (402.50)</td>
<td valign="top" align="center">68. (86.5)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
<td valign="top" align="center">567 (362.5)</td>
<td valign="top" align="center">116.5 (120)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S1 d.0</bold>
<break/>
<bold>Over cutoff (98 MFI) <sup>B</sup>
</bold>
</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">0 (0)</td>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S1 d.10 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">52 (34)</td>
<td valign="top" align="center">32 (18)</td>
<td valign="top" align="center">
<bold>0.003*</bold>
</td>
<td valign="top" align="center">53 (66)</td>
<td valign="top" align="center">38 (49.3)</td>
<td valign="top" align="center">0.105</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S1 d.21 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">332 (426)</td>
<td valign="top" align="center">161 (273)</td>
<td valign="top" align="center">0.076</td>
<td valign="top" align="center">244.5 (308)</td>
<td valign="top" align="center">107 (163.5)</td>
<td valign="top" align="center">
<bold>0.020*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S1 d.35 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">2512 (2141)</td>
<td valign="top" align="center">400 (1618)</td>
<td valign="top" align="center">
<bold>0.001*</bold>
</td>
<td valign="top" align="center">2362 (4715.8)</td>
<td valign="top" align="center">941 (1859.5)</td>
<td valign="top" align="center">
<bold>0.004*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S1 mo.6 <sup>A</sup> (MFI)</bold>
</td>
<td valign="top" align="center">273 (216)</td>
<td valign="top" align="center">56 (59.5)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
<td valign="top" align="center">321.5 (200)</td>
<td valign="top" align="center">85 (61.3)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
</tr>
<tr>
<th valign="top" colspan="7" align="left">Fold change relative to baseline</th>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.10-fold change <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">2.3 (10.6)</td>
<td valign="top" align="center">1.1 (0.9)</td>
<td valign="top" align="center">
<bold>0.010*</bold>
</td>
<td valign="top" align="center">2.3 (4.6)</td>
<td valign="top" align="center">1.3 (2.2)</td>
<td valign="top" align="center">0.351</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.21-fold change <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">18 (22.5)</td>
<td valign="top" align="center">7.8 (13.4)</td>
<td valign="top" align="center">
<bold>0.028*</bold>
</td>
<td valign="top" align="center">12.4 (13.7)</td>
<td valign="top" align="center">5.4 (11.4)</td>
<td valign="top" align="center">
<bold>0.033*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f d.35-fold change <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">102 (116.8)</td>
<td valign="top" align="center">23.8 (67.3)</td>
<td valign="top" align="center">
<bold>0.002*</bold>
</td>
<td valign="top" align="center">133.5 (274.8)</td>
<td valign="top" align="center">37.2 (66.5)</td>
<td valign="top" align="center">
<bold>0.002*</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>Anti-S-f mo.6-fold change <sup>A</sup>
</bold>
</td>
<td valign="top" align="center">15.2 (15.4)</td>
<td valign="top" align="center">2.5 (3)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
<td valign="top" align="center">18 (15.4)</td>
<td valign="top" align="center">4.2 (4)</td>
<td valign="top" align="center">
<bold>&lt;0.001*</bold>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Mann&#x2013;Whitney U test was applied to compare the continuous variables and X<sup>2</sup> test for analyzing categorical variables using SPSS version 24.0 (IBM Crop, Armonk, NY, USA) software. <bold>
<sup>A</sup>
</bold> variable is illustrated as the median (IQR), whereas <bold>
<sup>B</sup>
</bold> variable is illustrated as the number (%). <bold>*</bold> Bolded values denotes statistical significance when P-value is &lt; 0.05. NA, not applicable; IQR, interquartile range; d, day; mo, months; S1, spike antigen; S-f, full-length trimeric spike; MFI, mass fluorescence intensity.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s2_2">
<title>Microbiome richness and diversity in saliva of participants</title>
<p>We next investigated if the oral microbiome composition correlated with the vaccine responses noted in these participants. Baseline saliva samples were sequenced to address this question by subjecting salivary DNA to Illumina 16S rRNA gene sequencing (V3-V4 region). Output data in ASV-classified format were used for all downstream analysis. As shown, the salivary microbiome composition of healthy and PLHIV at genus level showed an overall dominance of <italic>Prevotella</italic>, <italic>Veillonella</italic>, and <italic>Neisseria</italic>, with the top identified 15 genera present in all participants (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>). There were no significant differences in richness and evenness in the microbial communities between High- <italic>vs</italic>. Low-responders in either cohort (alpha diversity: Observed and Chao1 indices or Shannon and Simpson&#x2019;s indices) (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>). Consistently with previously reports (<xref ref-type="bibr" rid="B15">15</xref>), higher microbial diversity (both alpha and beta) was found in HC as compared to the PLHIV participants (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S2</bold>
</xref>). The beta diversity analysis further indicated that there were interpersonal variations between the high- and low-responders in PLHIV participants (Bray Curtis and PERMANOVA: p=0.019; Jaccard index and PERMANOVA: p=0.022) as shown by the non-metric multidimensional scaling (NMDS)-based Bray-Curtis dissimilarity distances and Jaccard index <bold>(</bold>
<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>, left), while no significant difference was noted between low and high responder in HC participants (Bray Curtis and PERMANOVA: p=0.371; Jaccard index PERMANOVA: p=0.252) (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>, right).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>The oral microbiome richness and diversity and interpersonal variations in High- <italic>vs</italic>. Low responders. <bold>(A)</bold> Microbiome composition profiles at the genus level in low- respectively high-responders in the healthy- and PLHIV participants. It illustrates the frequently occurring genera, the &#x201c;top 15&#x201d; are present in all subjects in both subgroups of healthy and PLHIV, and the overall dominance of <italic>Prevotella</italic> (green), <italic>Veillonella</italic> (yellow), and <italic>Neisseria</italic> (blue). Vertical bars represent individual samples. <bold>(B)</bold> Scatterplots of alpha diversity of Observed and Chao1 indices of the ASV abundance, Shannon and Simpson&#x2019;s indices of the diversity of ASV among the participants. Lines and error bars indicate geographic means and standard deviation. <bold>(C)</bold> Non-metric multidimensional scaling (NMDS) plots visualising the beta diversity represented with Bray-Curtis dissimilarity distances and Jaccard index, validated by PERMANOVA test. ns, not significant.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1079995-g001.tif"/>
</fig>
</sec>
<sec id="s2_3">
<title>Taxonomic differences of oral microbiome in Low- and High-responders</title>
<p>To find the bacteria taxa that were differentially abundant between the High- and Low responders, Linear discriminant analysis Effect Size (LEfSe) analysis was applied for the HC and PLHIV data separately. Shown in <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>, the analysis of the taxonomic cladograms identified several significant differentially abundant taxa between the groups. The assigned LDA scores further showed that Healthy Low-responders (Healthy_Lo) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>) had an increased abundance of <italic>Campylobacter</italic>, <italic>Selenomonas and Butyrivibrio 2</italic>, and reduced abundance of <italic>Abiotrophia</italic> and <italic>Corynebacterium</italic>, as compared to the High responders (Healthy_Hi). As for the PLHIV participants (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>), Low-responders (PLHIV_Lo) had an increased abundance of <italic>Prevotella</italic>, <italic>Stomatobaculum, Lachnoanaerobaculum, Megasphaera</italic> and <italic>Leptotrichia</italic>, while <italic>Rothia, Granulicatella, Haemophilus</italic>, and <italic>Gemella</italic> were lower, as compared to High-responders (PLHIV_Hi).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Oral microbial signatures associated with long-term antibody responses in HC and PLWH. <bold>(A)</bold> LEfSe&#x2019;s cladogram shows the taxonomic levels, with the outer circle representing the phyla and the inner circle the genera. Each circle represents a taxa member within that taxonomic level. Green label indicates the high responders and red label the low responders. <bold>(B)</bold> Linear discriminant analysis (LDA) effect size analysis (LEfSe) identified the most differentially abundant genera between high and low responders in the healthy and PLHIV, respectively (P &lt;.05; LDA score &gt; 2). High responder-associated genera are indicated with negative LDA scores (green), and low responder-associated genera indicated with positive LDA scores (red). <bold>(C)</bold> The Area Under the Receiver Operating Characteristic (AUROC) scores show the predictive values of the individual microbial feature, or the combined significant features to predict the type of response e.g. Healthy_Hi (n=28) <italic>vs</italic>. Healthy_Lo (n=29), or PLHIV_Hi (n=20) <italic>vs</italic> PLHIV_Lo (n=38) among all vaccinees. Orange dots indicate AUC &gt; 0.70 and gray dots indicate AUC &lt; 0.70.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1079995-g002.tif"/>
</fig>
<p>We further determined the predictive values of the identified microbial features as a validation and to address how well they could distinguish the outcome of vaccination responses in the participants. The results obtained from area under the receiver operating characteristic curves (AUROC) indicated that these microbial features individually (<italic>Campylobacter</italic> and <italic>Butyrivibrio 2</italic>) provided predictive value of 0.810 and 0.703 respectively (AUROC; 95% CI, p&lt;0.0001 and p&lt;0.005) for salivary spike-IgG responses among HC participants. When combining all five significant bacteria features, the AUROC score increased to 0.879 (p&lt;0.0001) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>, right panel). Consistent with these findings, we found in PLHIV that <italic>Stamotabaculum</italic>, <italic>Granulicatella, Prevotella</italic> individually yielded predictive values of 0.710, 0705 to 0.703, respectively (p=0.008, 0.0112, 0.0107, respectively). When combining all nine significant bacteria from PLHIV, the AUC score increased to 0.82 (p&lt;0.05) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>, right panel).</p>
<p>We considered that there could be functional resemblances beyond the identified bacteria taxa. Striking, we observed functional associations indicating that, among Low-responders of HC as well as of PLHIV, the enriched taxa were mainly of anaerobic, gram-negative (lipopolysaccharide LPS+) bacteria species with known proteolytic activities (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). In a subsequent KEGG-pathway analysis (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>
<bold>)</bold>, they also showed significant positive associations with processes of amino acid metabolism, nucleotide metabolism and biosynthesis of other secondary metabolism (p &lt; 0.05, FDR &lt;0.05). On the contrary, in High-responders of both HC and PLHIV, the enriched species were instead mainly gram-positive bacteria of facultative genera with rather limited proteolytic activities. These bacteria were positively associated with carbohydrate metabolism, metabolism of other amino acids, vitamin and cofactor metabolism, and xenobiotics biodegradation and metabolism (p &lt; 0.05, FDR &lt;0.05) (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). Altogether, these data suggest that oral microbiome signatures in Low-responders among both HC and PLHIV cohorts resemble those described for a dysbiotic salivary community (<xref ref-type="bibr" rid="B9">9</xref>), and are distinguishable from the High-responders.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Overview of microbial features, enriched in high- and low-responders to COVID-19 vaccine, per study group (HC and PLHIV).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1079995-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Oral microbial signature and KEGG metabolism associated with High <italic>vs</italic>. Low responders in HC <bold>(A)</bold> or PLHIV <bold>(B)</bold>. Heatmap representation of KEGG metabolism pathways result on the microbial signatures of High (Hi) and Low responder (Low) that were selected by LEfSe v1.1.01 (linear discriminant analysis (LDA) at genus level (score &gt; 2, P &lt; 0.05) identified in Healthy and PLHIV participants. Red = positive correlations and blue = negative correlations by Spearman correlation (rho) at significance level of p &lt; 0.05 and FRD &lt;0.05.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1079995-g004.tif"/>
</fig>
</sec>
</sec>
<sec id="s3" sec-type="discussion">
<title>Discussion</title>
<p>The microbiota is fundamental for health and the evolution of the immune system with the microbiota is interconnected (<xref ref-type="bibr" rid="B13">13</xref>). We here describe the baseline oral microbiome composition with immunogenicity follow-up in a longitudinal cohort of healthy or people living with HIV vaccinated with BNT162b2 mRNA vaccine. We found that numerous features in the oral microbiome, including the diversity, specific bacteria taxa, and functions, correlated significantly with the immunogenicity of mRNA COVID-19 vaccination - in a primary anatomical site that needs immune protection against the SARS-CoV-2 infection (<xref ref-type="bibr" rid="B1">1</xref>). Our finding is in line with recent studies indicating that the gut microbiome or antibiotics usage could influence the systemic immunogenicity of the same mRNA vaccine (<xref ref-type="bibr" rid="B11">11</xref>) studied in the present study. We found that the durability of the salivary as well as the serum immunity through the mRNA-vaccination (two doses) lasted for at least six months. But it was negatively affected by the presence of a dysbiotic oral microbiome signature in both healthy and HIV participants who participated in our clinical trial (<xref ref-type="bibr" rid="B16">16</xref>). Here, we identified that the oral microbiota signature of high responders in HC and PLHIV individuals harboured significantly higher abundances of gram-positive facultative anaerobes saccharolytic bacteria, particularly <italic>Corynebacterium</italic> spp., <italic>Abiotrophia</italic> spp., <italic>Rothia</italic> spp., <italic>Gemella</italic> spp., and <italic>Granulicatella</italic> spp. In contrast, the low responders showed significantly higher abundances of gram-negative rod-shaped anaerobic proteolytic bacteria, including <italic>Campylobacter</italic> spp., <italic>Selenomonas</italic> spp., <italic>Butyrivibrio</italic> spp., <italic>Leptotrichia</italic> spp., <italic>Megasphaera</italic> spp., and <italic>Prevotella</italic> spp. Our results further indicate that a combined oral bacterial panel has the highest ability to predict the antibody magnitude and duration in saliva following the mRNA vaccination, which is consistent with the recent gut microbiome study on a one-month follow-up of COVID-19 vaccinees (<xref ref-type="bibr" rid="B11">11</xref>). Besides that, our findings also suggest that the durability of vaccine-induced immunity in the oral cavity could be influenced by the baseline oral microbiome of the vaccinees up to six months.</p>
<p>That the oral microbiome community might regulate local vaccine-induced mucosal immunity is intriguing and to our knowledge, similar results have not been reported before for other COVID-vaccines. The findings are of interest for vaccine strategies that aim at improving the mucosal immune memory. Although the precise mechanistic role of oral microbiota in vaccine response is unknown, several potential mechanisms could explain the link between deferentially enriched oral taxa and the persistence of salivary spike-IgG in the vaccinees. While most salivary IgG antibodies are derived from the bloodstream by passive leakage of the periodontal epithelium, there can also be some local production by salivary gland plasma cells (<xref ref-type="bibr" rid="B17">17</xref>). Data from molecular studies suggest that extracellular receptors such as toll-like receptor 2 (TLR2) and cytoplasmic receptors such as nucleotide-binding oligomerization domain protein 2 (NOD2) can recognize Gram-positive peptide fragments or metabolites released by bacteria that trigger the NF&#x3ba;B pathway (<xref ref-type="bibr" rid="B18">18</xref>). This pathway induces further phosphorylation of tight junction proteins, which promotes the expression of polymeric Ig receptor (pIgR), controlling the rate of production, translocation, and secretion of salivary immunoglobulins. Therefore, the high abundance of Gram-positive facultative bacteria in high responders may explain the durability of salivary spike IgG levels. Previous clinical studies have also shown an inverse relationship between oral diseases such as tooth decay, oral mucositis, and salivary immunoglobulin antibodies (<xref ref-type="bibr" rid="B19">19</xref>, <xref ref-type="bibr" rid="B20">20</xref>). Further, pathobionts &#x201c;formerly known as periodontitis-associated bacteria&#x201d; and abundant also in patients with edentulism, play unique and synergistic roles in dysbiosis of oral microbiota (<xref ref-type="bibr" rid="B21">21</xref>). Intriguingly, such pathobionts, especially the strictly anaerobe LPS-producing rods (<italic>Campylobacter</italic> spp., <italic>Selenomonas</italic> spp., <italic>Butyrivibrio</italic> spp., <italic>Leptotrichia</italic> spp., <italic>Megasphaera</italic> spp., and <italic>Prevotella</italic> spp.) are enriched in the saliva of the low responders. These pathobionts are known for potent proteolytic activity that can break down immunoglobulins, complements, and other innate defence proteins (<xref ref-type="bibr" rid="B22">22</xref>), and an inflammation-type dysbiosis has similarly been associates with long COVID (<xref ref-type="bibr" rid="B23">23</xref>). Whether oral health-related interventions or inhibition of specific microbial adhesion, pro-inflammatory mechanisms can further improve the duration of protective antibodies in the oral cavity during systemic or even mucosal vaccination, therefore deserves further investigation.</p>
<p>Our study design had considered age, sex, and BMI, which are all associated with vaccine efficacy against SARS-CoV-2 (<xref ref-type="bibr" rid="B24">24</xref>), as well as antibiotic use that could reduce COVID-19 vaccine immunogenicity (<xref ref-type="bibr" rid="B12">12</xref>). Those with exposure to SARS-CoV-2 infection before and during the study were omitted to reduce the potential bias of natural infection affecting vaccine-related antibody levels. All test platforms and clinical trial related procedures were also highly standardised and monitored (<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B16">16</xref>). The present study is not without limitations, only baseline microbiome timepoint was characterised only, the sample size is relatively small, oral health data are lacking, the microbiome analysis was a targeted approach, and only binding IgG antibodies up to six-months after vaccination were measured.</p>
<p>In conclusion, our study is the first to demonstrate that oral microbiome may have a role in maintaining the long-term antibody persistence in saliva and in blood after COVID-19 mRNA vaccination. Including oral health- or microbiome-targeted interventions to improve the long-term mucosal memory after vaccination should be further explored.</p>
</sec>
<sec id="s4" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s4_1">
<title>Study participants</title>
<p>Among 177 participants, including healthy controls (HC) or people living with HIV (PLHIV) who received two doses of mRNA vaccine Comirnaty<sup>&#xae;</sup> in COVAXID Clinical Trial (NCT04780659), 115 participants (HC, n=57; PLHIV, n=58) met study criteria for the current microbiome study. They were tested on day 0 and were confirmed negative on nasopharynx SARS-CoV-2 by PCR, and were seronegative following serological analysis with the Elecsys<sup>&#xae;</sup> anti-SARS CoV 2 S assay. None had recent or ongoing antibiotic treatment. All participants received two vaccine doses on study day 1 and 21, completed all baseline samplings, and tested negative for SARS-CoV-2 nucleocapsid antibodies at 6 months. Approval has been obtained by the Swedish Medical Product Agency (ID 5.1-2021-5881) for conducting the COVAXID clinical trial, and ethical permit was granted by the Swedish Ethical Review Authority (ID 2021-00451 and 2020-06381). All participants provided written informed consents.</p>
</sec>
<sec id="s4_2">
<title>Sample collection</title>
<sec id="s4_2_1">
<title>Sample collection and SARS-CoV-2 antibody detection in saliva</title>
<p>All saliva samples were processed by a standardized protocol in the same laboratory. Briefly, unstimulated whole saliva was self-collected by fasted study participants as described earlier using standardized picture instructions (<xref ref-type="bibr" rid="B4">4</xref>). After five minutes of passive drooling, the saliva was aliquoted in tubes using sterile transfer pipettes. All samples were immediately placed at 4&#xb0;C upon the same day and stored at -80&#xb0;C. Prior to antibody analysis, saliva samples were thawed at 4&#xb0;C and centrifuged at 400 xg for one min at 4&#xb0;C to separate any debris. Antibody analysis was performed using inactivated saliva (56&#xb0;C for 30&#xa0;min) as described earlier (<xref ref-type="bibr" rid="B14">14</xref>). Briefly, antibodies binding to the full-length spike glycoprotein in trimeric form (S-f) as well as the S1 subunit were measured by means of a multiplex bead-based assay in the 384-well plate format. The antigens were immobilized on the surface of uniquely color-coded bead identities (IDs) (MagPlex-C, Luminex corp.), and the IDs pooled to generate the bead-array. Saliva samples were diluted 1:5 in assay buffer and incubated with the array. After cross-linking the antibody-antigen complexes, a R-phycoerythrine-conjugated anti-human IgG antibody (H10104, Invitrogen) was applied for detection of IgG bound to spike. The assay readout was performed using a FlexMap3D instrument and the Luminex xPONENT software (Luminex Corp.). Each assay run included the same set of 12 negative and 4 positive saliva controls. Positive controls were samples from convalescent individuals with mild COVID-19 showing clear reactivity to spike. Negative controls were pre- pandemic saliva samples that were used to calculate the assay specific cutoffs and inter-assay variability. The inter-assay variability, evaluated as the % CV of the 16 control samples included in each assay run, was 10.8% for Spike-f and 12% for Spike S1 on average.</p>
</sec>
<sec id="s4_2_2">
<title>Sample collection and SARS-CoV-2 antibody detection in serum</title>
<p>Serum samples were analyzed for detection of antibodies to SARS-CoV-2 spike using the quantitative Elecsys<sup>&#xae;</sup> Anti-SARS-CoV-2 S test (Roche Diagnostics) (<xref ref-type="bibr" rid="B25">25</xref>) on the Cobas 8000 e801pro. The measuring range was between 0.40 to 250 U/mL, and the cut-off value for positive results is &#x2265; 0.80 U/mL Positive samples with antibody titers of &gt;250 U/mL were re-tested following 1/10 dilution, and in some cases 1/100 dilution with the upper level of measuring range 25,000 U/mL.</p>
</sec>
</sec>
</sec>
<sec id="s5">
<title>Antibody quantification and data analysis</title>
<p>The salivary antibody data were acquired as median fluorescence intensities (MFI) for each sample and antigen. The antigen and assay specific cutoff for positivity was calculated as the mean plus 6x standard deviation (SD) of the intensity signals of the 12 selected negative controls. The inter-assay variability was estimated for Spike-f and S1 as the average percent CV of the 16 control samples included in all 6 assay runs required to test the samples included in the current study. Statistical analyses (except microbiome analysis) were performed with Prism software v.9 (GraphPad) and SPSS version 24.0 (IBM Corp, Armonk, NY, USA) Statistics. Datasets initially underwent a data normality distribution test. Differences between groups of samples were analyzed by Mann-Whitney U test for univariate analysis. Correlations were determined using Spearman rank correlation. Two-sided p values &lt;0.05 were considered significant.</p>
<sec id="s5_1">
<title>DNA extraction and sequencing</title>
<p>DNA from saliva samples has been extracted using the ZymoBIOMICS&#x2122; DNA Mini Kit (Zymo Research, Irvine, California, United States) according to the manufacturer&#x2019;s instructions. Elution was performed in 60 mL of RNAse-free water and resulting DNA was stored at &#x2212;20&#xb0;C prior to preparation for sequencing. All samples were subsequently normalized to a standard concentration of 1 ng/mL and a volume of 15 mL (for a total of 15 ng of DNA) prior to Illumina 16S sequencing (V3-V4) at KI SciLife Laboratory on MiSeq. The 16S rRNA analysis was performed using the nf-core/ ampliseq analysis pipeline (<xref ref-type="bibr" rid="B26">26</xref>, <xref ref-type="bibr" rid="B27">27</xref>).</p>
<sec id="s5_1_1">
<title>Bioinformatics and statistics analysis</title>
<p>Raw RNA-Seq data were quality checked using FastQC v0.11.8 and then pre-processed using Cutadapt v2.8 to remove adapter sequences and poor-quality bases. The pre-processed sequencing reads were processed using QIIME2 v2019.10.0 and were denoised using DADA2 and converted to amplicon sequence variants (ASVs). Non-bacteria taxa e.g., eukaryota, chloroplast and mitochondria in the sequence output by DADA2 were removed. Taxonomic assignments were made using the Silva v132 annotation. DADA2 detected 2456 ASVs of which 278 bacterial taxa were recovered. A total of 2653288 feature counts were detected from all samples with an average of 23072 count per sample (range: 2101 - 318474 counts). A filter was applied to remove counts &lt; 2 and minimum prevalence of 10% per sample with low variance (5%) (<xref ref-type="bibr" rid="B28">28</xref>). The total sum of scales (TSS) normalisation and subsequent rarefication for sample depth normalization were applied to create the final ASV-based feature count table for downstream analysis. Tax4Fun R package (<xref ref-type="bibr" rid="B29">29</xref>) was used for predicting the functional profiles. Comparative analysis of the abundance at the individual level of subgroups was done using TSS normalized, ASV-based feature count table and visualized by the miaViz 4.2 R package. Alpha diversity metrics such as Observed, Chao1, Shannon and Simpson were analysed in QIIME2 v2019.10.0. Differential microbial communities (in - between or beta diversity) were assessed by Microbiome Analyst, using Bray&#x2013;Curtis- and Jaccard index-based non-metric multidimensional scaling (NMDS) and significance test of PERMANOVA. Differentially abundant genera analysis was done using Linear Discriminant Analysis Effect Size (LEfSe v1.1.01) in in the Galaxy web application (<uri xlink:href="http://huttenhower.sph.harvard.edu/galaxy">http://huttenhower.sph.harvard.edu/galaxy</uri>) (<xref ref-type="bibr" rid="B30">30</xref>). Correlation analysis was performed with Spearman correlation test and the receiver operating characteristic curve (ROC) was applied. The CombiROC (<uri xlink:href="https://ingmbioinfo.github.io/combiroc/">https://ingmbioinfo.github.io/combiroc/</uri>) was used for selection of combination of biomarker features and ROC was visualized in MATLAB. Subgroup analysis of KEGG (Kyoto Encyclopedia of Genes and Genomes) metabolisms was performed with LEfSe v1.1.01 (linear discriminant analysis (LDA) selected features at genus level (score &gt; 2, P &lt; 0.05), followed by Spearman correlation test (P &lt; 0.05, FDR &lt; 0.05).</p>
</sec>
</sec>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: <uri xlink:href="https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA900274">https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA900274</uri>.</p>
</sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The studies involving human participants were reviewed and approved by Swedish Ethical Review Authority. The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>Conception, planning, funds acquisition: PN, H-GL, PL. SA, MSC. Planning and organised clinical work: PB, PL, PC, PN, SA. Performed laboratory analyses or statistics: MG, KA-M, SN, KH, GG, MS, SR, MA, SM, GB, EP. Conclusion and interpretation of result: MG, KA-M, EP, PN, SA, MSC. Wrote first draft: MG, KA-M, MSC. Reviewed and edited final version: All co-authors. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>Swedish Research Council: 2021-04779. Knut and Alice Wallenberg Foundation: VC-2021-0018. Center for Innovative Medicine (CIMED): FoUI-963079.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We thank all the participants enrolled in this study, and all the research and clinical staff at Karolinska University Hospital. The authors acknowledge support from the Autoimmunity and Serology Profiling Unit at Science for Life Laboratory, as well as the National Genomics Infrastructure in Stockholm, and SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science for assistance with massively parallel sequencing and access to the UPPMAX computational infrastructure.</p>
</ack>
<sec id="s10" sec-type="COI-statement">
<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>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12" sec-type="supplementary-material">
<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/fimmu.2022.1079995/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fimmu.2022.1079995/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>N</given-names>
</name>
<name>
<surname>P&#xe9;rez</surname> <given-names>P</given-names>
</name>
<name>
<surname>Kato</surname> <given-names>T</given-names>
</name>
<name>
<surname>Mikami</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Okuda</surname> <given-names>K</given-names>
</name>
<name>
<surname>Gilmore</surname> <given-names>R. C</given-names>
</name>
<etal/>
</person-group>. <article-title>SARS-CoV-2 infection of the oral cavity and saliva</article-title>. <source>Nat Med</source> (<year>2021</year>) <volume>27</volume>:<fpage>892</fpage>&#x2013;<lpage>903</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41591-021-01296-8</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pilishvili</surname> <given-names>T</given-names>
</name>
<name>
<surname>Gierke</surname> <given-names>R</given-names>
</name>
<name>
<surname>Fleming-Dutra</surname> <given-names>KE</given-names>
</name>
<name>
<surname>Farrar</surname> <given-names>JL</given-names>
</name>
<name>
<surname>Mohr</surname> <given-names>NM</given-names>
</name>
<name>
<surname>Talan</surname> <given-names>DA</given-names>
</name>
<etal/>
</person-group>. <article-title>Effectiveness of mRNA covid-19 vaccine among U.S. health care personnel</article-title>. <source>N Engl J Med</source> (<year>2021</year>) <volume>385</volume>:<elocation-id>e90</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1056/NEJMoa2106599</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dagan</surname> <given-names>N</given-names>
</name>
<name>
<surname>Barda</surname> <given-names>N</given-names>
</name>
<name>
<surname>Kepten</surname> <given-names>E</given-names>
</name>
<name>
<surname>Miron</surname> <given-names>O</given-names>
</name>
<name>
<surname>Perchik</surname> <given-names>S</given-names>
</name>
<name>
<surname>Katz</surname> <given-names>MA</given-names>
</name>
<etal/>
</person-group>. <article-title>BNT162b2 mRNA covid-19 vaccine in a nationwide mass vaccination setting</article-title>. <source>N Engl J Med</source> (<year>2021</year>) <volume>384</volume>:<page-range>1412&#x2013;23</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1056/NEJMoa2101765</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Healy</surname> <given-names>K</given-names>
</name>
<name>
<surname>Pin</surname> <given-names>E</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>P</given-names>
</name>
<name>
<surname>S&#xf6;derdahl</surname> <given-names>G</given-names>
</name>
<name>
<surname>Nowak</surname> <given-names>P</given-names>
</name>
<name>
<surname>Mielke</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Salivary IgG to SARS-CoV-2 indicates seroconversion and correlates to serum neutralization in mRNA-vaccinated immunocompromised individuals</article-title>. <source>Med (N Y)</source> (<year>2022</year>) <volume>3</volume>:<fpage>137</fpage>&#x2013;<lpage>153.e133</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.medj.2022.01.001</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaan</surname> <given-names>AMM</given-names>
</name>
<name>
<surname>Kahharova</surname> <given-names>D</given-names>
</name>
<name>
<surname>Zaura</surname> <given-names>E</given-names>
</name>
</person-group>. <article-title>Acquisition and establishment of the oral microbiota</article-title>. <source>Periodontol 2000</source> (<year>2021</year>) <volume>86</volume>:<page-range>123&#x2013;41</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/prd.12366</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wade</surname> <given-names>WG</given-names>
</name>
</person-group>. <article-title>Resilience of the oral microbiome</article-title>. <source>Periodontol 2000</source> (<year>2021</year>) <volume>86</volume>:<page-range>113&#x2013;22</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/prd.12365</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>Crosstalk between the oral microbiota, mucosal immunity, and the epithelial barrier regulates oral mucosal disease pathogenesis</article-title>. <source>Mucosal Immunol</source> (<year>2021</year>) <volume>14</volume>:<page-range>1247&#x2013;58</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41385-021-00413-7</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Acharya</surname> <given-names>C</given-names>
</name>
<name>
<surname>Sahingur</surname> <given-names>SE</given-names>
</name>
<name>
<surname>Bajaj</surname> <given-names>JS</given-names>
</name>
</person-group>. <article-title>Microbiota, cirrhosis, and the emerging oral-gut-liver axis</article-title>. <source>JCI Insight</source> (<year>2017</year>) <volume>2</volume>:<fpage>94416</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1172/jci.insight.94416</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>J</given-names>
</name>
<name>
<surname>Domingue</surname> <given-names>JC</given-names>
</name>
<name>
<surname>Sears</surname> <given-names>CL</given-names>
</name>
</person-group>. <article-title>Microbiota dysbiosis in select human cancers: Evidence of association and causality</article-title>. <source>Semin Immunol</source> (<year>2017</year>) <volume>32</volume>:<fpage>25</fpage>&#x2013;<lpage>34</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.smim.2017.08.001</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kinane</surname> <given-names>DF</given-names>
</name>
<name>
<surname>Stathopoulou</surname> <given-names>PG</given-names>
</name>
<name>
<surname>Papapanou</surname> <given-names>PN</given-names>
</name>
</person-group>. <article-title>Periodontal diseases</article-title>. <source>Nat Rev Dis Primers</source> (<year>2017</year>) <volume>3</volume>:<fpage>17038</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nrdp.2017.38</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ng</surname> <given-names>SC</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L</given-names>
</name>
<name>
<surname>Mok</surname> <given-names>CK</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>S</given-names>
</name>
<name>
<surname>Li</surname> <given-names>A</given-names>
</name>
<etal/>
</person-group>. <article-title>Gut microbiota composition is associated with SARS-CoV-2 vaccine immunogenicity and adverse events</article-title>. <source>Gut</source> (<year>2022</year>) <volume>71</volume>:<page-range>1106&#x2013;16</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1136/gutjnl-2021-326563</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheung</surname> <given-names>KS</given-names>
</name>
<name>
<surname>Lam</surname> <given-names>LK</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>R</given-names>
</name>
<name>
<surname>Ooi</surname> <given-names>PH</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>JT</given-names>
</name>
<name>
<surname>To</surname> <given-names>WP</given-names>
</name>
<etal/>
</person-group>. <article-title>Association between recent usage of antibiotics and immunogenicity within six months after&#xa0;COVID-19 vaccination</article-title>. <source>Vaccines (Basel)</source> (<year>2022</year>) <volume>10</volume>:<fpage>1122</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/vaccines10071122</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hooper</surname> <given-names>LV</given-names>
</name>
<name>
<surname>Littman</surname> <given-names>DR</given-names>
</name>
<name>
<surname>Macpherson</surname> <given-names>AJ</given-names>
</name>
</person-group>. <article-title>Interactions between the microbiota and the immune system</article-title>. <source>Science</source> (<year>2012</year>) <volume>336</volume>:<page-range>1268&#x2013;73</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.1223490</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alkharaan</surname> <given-names>H</given-names>
</name>
<name>
<surname>Bayati</surname> <given-names>S</given-names>
</name>
<name>
<surname>Hellstr&#xf6;m</surname> <given-names>C</given-names>
</name>
<name>
<surname>Aleman</surname> <given-names>S</given-names>
</name>
<name>
<surname>Olsson</surname> <given-names>A</given-names>
</name>
<name>
<surname>Lindahl</surname> <given-names>K</given-names>
</name>
<etal/>
</person-group>. <article-title>Persisting salivary IgG against SARS-CoV-2 at 9 months after mild COVID-19: A complementary approach to population surveys</article-title>. <source>J Infect Dis</source> (<year>2021</year>) <volume>224</volume>(<issue>3</issue>):<page-range>407&#x2013;14</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/infdis/jiab256</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>S</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Su</surname> <given-names>B</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>H</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>F</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>H</given-names>
</name>
<etal/>
</person-group>. <article-title>Alteration in oral microbiome among men who have sex with men with acute and chronic HIV infection on antiretroviral therapy</article-title>. <source>Front Cell Infect Microbiol</source> (<year>2021</year>) <volume>11</volume>:<elocation-id>695515</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fcimb.2021.695515</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bergman</surname> <given-names>P</given-names>
</name>
<name>
<surname>Blennow</surname> <given-names>O</given-names>
</name>
<name>
<surname>Hansson</surname> <given-names>L</given-names>
</name>
<name>
<surname>Mielke</surname> <given-names>S</given-names>
</name>
<name>
<surname>Nowak</surname> <given-names>P</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>P</given-names>
</name>
<etal/>
</person-group>. <article-title>Safety and efficacy of the mRNA BNT162b2 vaccine against SARS-CoV-2 in five groups of immunocompromised patients and healthy controls in a prospective open-label clinical trial</article-title>. <source>EBioMedicine</source> (<year>2021</year>) <volume>74</volume>:<elocation-id>103705</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ebiom.2021.103705</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brandtzaeg</surname> <given-names>P</given-names>
</name>
</person-group>. <article-title>Do salivary antibodies reliably reflect both mucosal and systemic immunity</article-title>? <source>Ann N Y Acad Sci</source> (<year>2007</year>) <volume>1098</volume>:<fpage>288</fpage>&#x2013;<lpage>311</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1196/annals.1384.012</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brandtzaeg</surname> <given-names>P</given-names>
</name>
</person-group>. <article-title>Secretory immunity with special reference to the oral cavity</article-title>. <source>J Oral Microbiol</source> (<year>2013</year>) <volume>5</volume>:<fpage>20401</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3402/jom.v5i0.20401</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nogueira</surname> <given-names>RD</given-names>
</name>
<name>
<surname>Alves</surname> <given-names>AC</given-names>
</name>
<name>
<surname>Napimoga</surname> <given-names>MH</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>DJ</given-names>
</name>
<name>
<surname>Mattos-Graner</surname> <given-names>RO</given-names>
</name>
</person-group>. <article-title>Characterization of salivary immunoglobulin a responses in children heavily exposed to the oral bacterium streptococcus mutans: influence of specific antigen recognition in infection</article-title>. <source>Infect Immun</source> (<year>2005</year>) <volume>73</volume>:<page-range>5675&#x2013;84</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/IAI.73.9.5675-5684.2005</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pels</surname> <given-names>EJ</given-names>
</name>
</person-group>. <article-title>Oral mucositis and saliva IgA, IgG and IgM concentration during anti-tumor treatment in children suffering from acute lymphoblastic leukemia</article-title>. <source>Adv Clin Exp Med</source> (<year>2017</year>) <volume>26</volume>:<page-range>1351&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.17219/acem/64940</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hunter</surname> <given-names>MC</given-names>
</name>
<name>
<surname>Pozhitkov</surname> <given-names>AE</given-names>
</name>
<name>
<surname>Noble</surname> <given-names>PA</given-names>
</name>
</person-group>. <article-title>Microbial signatures of oral dysbiosis, periodontitis and edentulism revealed by gene meter methodology</article-title>. <source>J Microbiol Methods</source> (<year>2016</year>) <volume>131</volume>:<fpage>85</fpage>&#x2013;<lpage>101</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.mimet.2016.09.019</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiao</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Hasegawa</surname> <given-names>M</given-names>
</name>
<name>
<surname>Inohara</surname> <given-names>N</given-names>
</name>
</person-group>. <article-title>The role of oral pathobionts in dysbiosis during periodontitis development</article-title>. <source>J Dent Res</source> (<year>2014</year>) <volume>93</volume>:<page-range>539&#x2013;46</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1177/0022034514528212</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haran</surname> <given-names>JP</given-names>
</name>
<name>
<surname>Bradley</surname> <given-names>E</given-names>
</name>
<name>
<surname>Zeamer</surname> <given-names>AL</given-names>
</name>
<name>
<surname>Cincotta</surname> <given-names>L</given-names>
</name>
<name>
<surname>Salive</surname> <given-names>MC</given-names>
</name>
<name>
<surname>Dutta</surname> <given-names>P</given-names>
</name>
<etal/>
</person-group>. <article-title>Inflammation-type dysbiosis of the oral microbiome associates with the duration of COVID-19 symptoms and long COVID</article-title>. <source>JCI Insight</source> (<year>2021</year>) <volume>6</volume>:<elocation-id>e152346</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1172/jci.insight.152346</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levin</surname> <given-names>EG</given-names>
</name>
<name>
<surname>Lustig</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Cohen</surname> <given-names>C</given-names>
</name>
<name>
<surname>Fluss</surname> <given-names>R</given-names>
</name>
<name>
<surname>Indenbaum</surname> <given-names>V</given-names>
</name>
<name>
<surname>Amit</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Waning immune humoral response to BNT162b2 covid-19 vaccine over 6 months</article-title>. <source>N Engl J Med</source> (<year>2021</year>) <volume>385</volume>:<elocation-id>e84</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1056/NEJMoa2114583</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Higgins</surname> <given-names>V</given-names>
</name>
<name>
<surname>Fabros</surname> <given-names>A</given-names>
</name>
<name>
<surname>Kulasingam</surname> <given-names>V</given-names>
</name>
</person-group>. <article-title>Quantitative measurement of anti-SARS-CoV-2 antibodies: Analytical and clinical evaluation</article-title>. <source>J Clin Microbiol</source> (<year>2021</year>) <volume>59</volume>:<elocation-id>e03149&#x2013;20</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JCM.03149-20</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ewels</surname> <given-names>PA</given-names>
</name>
<name>
<surname>Peltzer</surname> <given-names>A</given-names>
</name>
<name>
<surname>Fillinger</surname> <given-names>S</given-names>
</name>
<name>
<surname>Patel</surname> <given-names>H</given-names>
</name>
<name>
<surname>Alneberg</surname> <given-names>J</given-names>
</name>
<name>
<surname>Wilm</surname> <given-names>A</given-names>
</name>
<etal/>
</person-group>. <article-title>The nf-core framework for community-curated bioinformatics pipelines</article-title>. <source>Nat Biotechnol</source> (<year>2020</year>) <volume>38</volume>:<page-range>276&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41587-020-0439-x</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Straub</surname> <given-names>D</given-names>
</name>
<name>
<surname>Blackwell</surname> <given-names>N</given-names>
</name>
<name>
<surname>Langarica-Fuentes</surname> <given-names>A</given-names>
</name>
<name>
<surname>Peltzer</surname> <given-names>A</given-names>
</name>
<name>
<surname>Nahnsen</surname> <given-names>S</given-names>
</name>
<name>
<surname>Kleindienst</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Interpretations of environmental microbial community studies are biased by the selected 16S rRNA (Gene) amplicon sequencing pipeline</article-title>. <source>Front Microbiol</source> (<year>2020</year>) <volume>11</volume>:<elocation-id>550420</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmicb.2020.550420</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nearing</surname> <given-names>JT</given-names>
</name>
<name>
<surname>Douglas</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Hayes</surname> <given-names>MG</given-names>
</name>
<name>
<surname>MacDonald</surname> <given-names>J</given-names>
</name>
<name>
<surname>Desai</surname> <given-names>DK</given-names>
</name>
<name>
<surname>Allward</surname> <given-names>N</given-names>
</name>
<etal/>
</person-group>. <article-title>Microbiome differential abundance methods produce different results across 38 datasets</article-title>. <source>Nat Commun</source> (<year>2022</year>) <volume>13</volume>:<fpage>342</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-022-28034-z</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>A&#xdf;hauer</surname> <given-names>KP</given-names>
</name>
<name>
<surname>Wemheuer</surname> <given-names>B</given-names>
</name>
<name>
<surname>Daniel</surname> <given-names>R</given-names>
</name>
<name>
<surname>Meinicke</surname> <given-names>P</given-names>
</name>
</person-group>. <article-title>Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data</article-title>. <source>Bioinformatics</source> (<year>2015</year>) <volume>31</volume>:<page-range>2882&#x2013;4</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btv287</pub-id>
</citation>
</ref>
<ref id="B30">
<label>30</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Segata</surname> <given-names>N</given-names>
</name>
<name>
<surname>Izard</surname> <given-names>J</given-names>
</name>
<name>
<surname>Waldron</surname> <given-names>L</given-names>
</name>
<name>
<surname>Gevers</surname> <given-names>D</given-names>
</name>
<name>
<surname>Miropolsky</surname> <given-names>L</given-names>
</name>
<name>
<surname>Garrett</surname> <given-names>WS</given-names>
</name>
<etal/>
</person-group>. <article-title>Metagenomic biomarker discovery and explanation</article-title>. <source>Genome Biol</source> (<year>2011</year>) <volume>12</volume>:<fpage>R60</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>