Next Article in Journal
Impact of Hydrogen Sulfide on Mitochondrial and Bacterial Bioenergetics
Next Article in Special Issue
Gut Microbiota Extracellular Vesicles as Signaling Molecules Mediating Host-Microbiota Communications
Previous Article in Journal
Pathophysiological Potentials of NRF3-Regulated Transcriptional Axes in Protein and Lipid Homeostasis
Previous Article in Special Issue
A Journey on the Skin Microbiome: Pitfalls and Opportunities
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 23
10.3390/ijms222312682
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Early-Life Adversity Leaves Its Imprint on the Oral Microbiome for More Than 20 Years and Is Associated with Long-Term Immune Changes

by
Eleftheria G. Charalambous
1,2,
Sophie B. Mériaux
1,
Pauline Guebels
1,
Claude P. Muller
1,
Fleur A. D. Leenen
1,
Martha M. C. Elwenspoek
1,
Ines Thiele
3,4,5,6,
Johannes Hertel
3,7 and
Jonathan D. Turner
1,*
1
Immune Endocrine and Epigenetics Research Group, Department of Infection and Immunity, Luxembourg Institute of Health (LIH), 29 Rue Henri Koch, L-4354 Esch-sur-Alzette, Luxembourg
2
Faculty of Science, Technology and Medicine, University of Luxembourg, 2 Avenue de Université, L-4365 Esch-sur-Alzette, Luxembourg
3
School of Medicine, National University of Ireland, H91 YR71 Galway, Ireland
4
Ryan Institute, National University of Galway, H91 TK33 Galway, Ireland
5
Division of Microbiology, National University of Galway, H91 TK33 Galway, Ireland
6
APC Microbiome Ireland, T12 HW58 Cork, Ireland
7
Department of Psychiatry and Psychotherapy, University Medicine Greifswald, 17489 Greifswald, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(23), 12682; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222312682
Submission received: 27 October 2021 / Revised: 19 November 2021 / Accepted: 22 November 2021 / Published: 24 November 2021
(This article belongs to the Special Issue The Interplay of Microbiome and Immune Response in Health and Diseases-2nd Edition)
Browse Figures
Versions Notes

Abstract

:
The early-life microbiome (ELM) interacts with the psychosocial environment, in particular during early-life adversity (ELA), defining life-long health trajectories. The ELM also plays a significant role in the maturation of the immune system. We hypothesised that, in this context, the resilience of the oral microbiomes, despite being composed of diverse and distinct communities, allows them to retain an imprint of the early environment. Using 16S amplicon sequencing on the EpiPath cohort, we demonstrate that ELA leaves an imprint on both the salivary and buccal oral microbiome 24 years after exposure to adversity. Furthermore, the changes in both communities were associated with increased activation, maturation, and senescence of both innate and adaptive immune cells, although the interaction was partly dependent on prior herpesviridae exposure and current smoking. Our data suggest the presence of multiple links between ELA, Immunosenescence, and cytotoxicity that occur through long-term changes in the microbiome.

1. Introduction

Early-life adversity (ELA) is defined by a poor environment and conditions in early life that induce intense psychophysiological stress [1]. It is mostly observed together with low socioeconomic status and is pathophysiologically correlated with a lifelong imbalance of health and disease [2]. The first 1000 days from conception to 2 years is the most vulnerable life period [3]. At birth, the body is almost fully formed; however, many biological systems continue to mature over the following years. Research on the lifelong health and disease balance has shown the significance of the environment during this period on multiple disease phenotypes [4], including cardiovascular, allergic, and autoimmune disorders, as well as mental disorders [5,6,7,8,9,10,11,12,13,14,15,16,17]. There has been a focus on the molecular mechanisms and the cellular phenotype behind the effect of stress and adversity on immune and endocrine systems as well as epigenetic modifications [5,6,18,19].
ELA has been reported to influence health trajectories via the immune system [18,20,21], with a clear ELA-associated immunophenotype centred around the activation and functional status of T lymphocytes. In the institutionalisation model of early-life stress, strong T-cell immunosenescence has been reported [18,21,22]. Immunosenescence is a form of accelerated immune ageing. The CD57 T- and NK- cell immunosenescence marker is absent in early life and increases with age, with high numbers of such cells in the elderly population. Immunosenescence is driven by chronic inflammation or recurrent viral infections such as CMV [23]. NK functionality is also highly impacted by recurrent reactivation of CMV inducing NK cell exhaustion, increased cytotoxicity, and senescence [24]. Additionally, such viral infections potentially program the immune system [21,22]. Latent CMV infection of haematopoietic progenitor cells reduces GR transcription and translation, impacting immune cell maturation, which can be dependent on CMV reactivation [19,25,26].
The environment is the most critical factor shaping ELA. In the immediate postnatal period, the immune system starts maturing and the first body-area-specific microbial communities are established. Once established, the microbiome modulates the host [27], a mechanism to protect symbiotic microbial communities, where cases of microbial dysbiosis can be fatal [28,29]. The ELM plays important roles in an infant’s subsequent development [5,6] and a long-term health trajectory [5,6,7,8,9,10,11]. Nevertheless, the development of the ELM is critically impacted by the environment, mode of birth, early-life nutrition, and environmental exposure, which leave a clear lifelong trace [16,30]. For example, maternal milk is rich in essential nutrients, protective antibodies, and components essential for the developing microbiome, such as human milk oligosaccharides (HMOs) and short-chain fatty acids (SCFAs) [31,32,33,34,35,36].
The oral microbiome (OM) is composed of various distinct, smaller communities within the oral cavity [37,38,39,40,41] that are robust, stable, and particularly resilient [42,43], particularly to antibiotic therapy [43,44,45,46,47]. Moreover, saliva contains actively secreted components such as cortisol, glucose, lactate, urea, and proteins, such as polypeptides, glycoproteins (cystatins, mucins, and immunoglobulins) and antimicrobial peptides (histatins, defensins, and immunoglobulins-IgA). Many of these are energy sources for the OM, and salivary glycoproteins are the principal nutrient source. These substrates are crucial for the development of multispecies communities and their preservation [41,46,48], and enhance the resistance of the community to environmental stressors [41,49,50]. The long-term stability of the OM leads to the hypothesis that, once established in early life, it remains stable, robust, and resilient, retaining an imprint of the early environment [42,43].
We previously reported higher virally mediated activation and senescence of the immune system in association with ELA in the EpiPath cohort [22,51]. This cohort consists of young adults exposed to ELA by either institutionalisation or separation from their parents at birth, and were subsequently adopted into Luxembourg, while controls were raised by their biological parents. With the growing evidence of a microbiome–immune–system interaction, we attempted to identify if institutionalisation left a mark on the OM of the adoptees. In this study, we sequenced the 16s-rRNA from the buccal and salivary bacterial communities from our cohort. Integrating this with the full immunophenotype, we identified associations with various taxa and analysed how the microbiome interacts with the immune system.

2. Results

The V4 region of the bacterial 16S gene was successfully amplified from both buccal swabs and salimetrics oral swabs for the 115 members of the EpiPath cohort, and a total of 45 Gbp sequencing data were obtained. All samples were successfully processed using mothur and a total of 288 and 371 genera from 24 phyla were identified for buccal and salivary samples respectively. The saliva and buccal swabs from the EpiPath cohort were examined independently as they are two distinct oral communities from the same participants.

2.1. Microbial Diversity and Overall Microbial Composition

2.1.1. Salivary Microbiome (α- and β-Diversity)

We identified sequences from all of the 24 principal bacterial phyla in the overall salivary microbial community. Within these 24 phyla, the most abundant in both control and ELA groups were Actinobacteria, Proteobacteria, Firmicutes, and Fusobacteriota (Figure 1A). The most abundant genera of the salivary composition were Acinetobacter, Micrococcaceae, Actinobacillus, Rothia, Corynebacterium, Micrococcus, Actinomyces, Alloprevotela, Porphymonas, Fusobacterium, Weeksellaceae, Flavobacterium, Bradyrhizobium, Porphyromonas, _ Comamonas, Olsenella, Fluviicola, Fusobacterium, Absconditabacteriales _(SR1)_ge, Streptobacillus, Fretibacterium, JGI_0000069-P22_ge, Capnocytophaga, Pseudarcicella, Tannerella, Prevotella, and Campylobacter (Figure 1B). There was no difference in α-diversity between the controls and the adoptees in terms of diversity (controls: mean = 13.87872, SD = 7.127721, adoptees: mean = 14.14869, SD = 6.601899, Wilcoxon rank sum test p = 0.7579) and evenness (controls: mean = 0.5619667, SD = 0.06864806, adoptees: mean = 0.5646776, SD = 0.06172601, Wilcoxon rank sum test p = 0.9461). Plotting the Shannon evenness index against the inverse Simpson diversity index confirmed that there analogous diversity and evenness between the controls and adoptees (Supplementary Figure S1). Principal coordinate analysis could not detect systematic differences either (Supplementary Figure S1).

2.1.2. Buccal Microbiome (α- and β-Diversity)

As for the salivary microbiome, we identified sequences from all the 24 principal bacterial phyla in the buccal microbial community. The most abundant phyla were the same as in the salivary microbiome: Actinobacteria, Proteobacteria, Firmicutes, and Fusobacteriota (Figure 2A). The most abundant genera were Actinomyces, Corynebacterium, Micrococcaceae, Rothia, Alloprevotela, Porphymonas, Fusobacterium, Actinobacillus, Bradyrhizobium, Haemophilus, Methylobacterium-Methylorubrum, Oxalobacteracea, Actinomyces, Neisseria, Paucibacter, Lautropia Cardiobacterium, Brucella, Alysiella, and Campylobacter, therefore revealing a substantial overlap in the detected genera between the saliva and the buccal microbiome (Figure 2B).
As for the salivary microbiome, we observed a very similar diversity and evenness between controls and adoptees as measured by the inverse Simpson index (controls: mean = 12.85828, SD = 10.117292, adoptees: mean = 14.72315, SD = 8.991065, Wilcoxon rank sum test p = 0.08803). The Shannon evenness index was similar between adoptees and controls, and, in both, it was higher than in the salivary microbiome (controls: mean = 0.5137724, SD = 0.10141122, adoptees: mean = 0.465861, SD = 0.08482336, Wilcoxon rank sum test p = 0.09024). α-diversity was again similar in the controls and the adoptees in terms of evenness (Wilcoxon rank sum test p = 0.08803) and diversity (Wilcoxon rank sum test p = 0.09024). Plotting the Shannon evenness index against the inverse Simpson diversity index revealed no systematic differences in diversity or evenness between the controls and adoptees (Supplementary Figure S1). Similarly, principal coordinate analysis indicated no differences between adoptees and control (Supplementary Figure S1).

2.1.3. Salivary and Buccal Microbiomes Are Two Separate Entities

To ensure that sample collection was performed correctly and that we had two distinct communities, we compared the diversity and evenness of the salivary and buccal samples. We found a low correlation between the salivary and buccal communities in both the inverse Simpson diversity index (Figure 3A, p = 0.47, rho = −0.07372058, Spearman’s rank correlation test) and Shannon evenness index (Figure 3B, p = 0.8759, rho = 0.01595802, Spearman’s rank correlation test), giving evidence in favour of the hypothesis that, despite their close physical proximity, they can be seen as distinct communities. Comparing the salivary and buccal microbiomes by group, the diversity ratios of the salivary against buccal communities were similar between the controls and adoptees (Figure 3A, controls: p = 0.3311, rho = −0.1222465; adoptees: p = 0.7898, rho = 0.04812834; Spearman’s rank correlation test, Figure 3B; controls: p = 0.9578, rho = 0.009692513; adoptees: p = 0.7898, rho = −0.04812834; Spearman’s rank correlation test). This suggests that the overall composition between the controls and adoptees may be similar, but differences would be seen at the phyla level.

2.2. ELA Induces Differences in Specific Taxa in Both Salivary and Buccal Communities

Investigating the abundance levels of phyla and genera highlighted differences in the community composition across the ELA group and healthy controls (Table 1). While there were no differences in the phyla level in the buccal data, Proteobacteria and Verrucomicrobiota were significantly lower (FDR < 0.05 for both) in the adoptees in comparison to the controls (Supplementary Table S1A, Figure 4A,B) in the saliva microbiome as detected in fractional regression analyses. Analyses of deeper taxonomy revealed two of the most abundant genera of the Proteobacteria phylum, Comamonas and Acinetobacter, to be significantly lower in the saliva of adoptees compared to controls alongside Aquabacterium and unclassified Comamonadaceae (Table S1B, Figure 4C–F). In conclusion, while we could not detect systematic differences in the buccal microbiome between the ELA group and the controls, the saliva microbiome was structurally different in its composition, with a prominent role for Proteobacteria genera (Figure 4A).

2.3. Environmental Covariates

Next, we investigated the impact of environmental factors on the OM to potentially explain the effects of ELA described above (Table 1).

2.3.1. Smoking

As lifestyle has a pivotal role in the development of the microbiome, we assessed the effect of smoking on the OM by including smoking status (binary: smokers vs. non-smokers) into the regression modelling (Table 1). No significant genera were detected for the salivary community, whereas from the buccal community, we exposed five genera negatively associated with smoking: three from the Proteobacteria phylum, Neisseria, Neisseriaceae_unclassified, and Pasteurellaceae_unclassified; 1 from the Bacteroidetes phylum and Capnocytophaga genus; as well as one from the Firmicutes phylum and Planococcaceae_unclassified genus (FDR < 0.05, Figure 5A–E). In sensitivity analyses, we removed smoking as a covariate from the regression equations for the FDR-corrected significant genera to explore potential effect mediation through smoking, but the results remained virtually unchanged. The full results for the buccal and saliva microbiomes can be found in Supplementary Table S2.

2.3.2. Prior Viral Infections

We previously reported that viral infections may mediate the early-life immunophenotype [22]. Consequently, we tested whether prior viral infection, measured as anti- HSV, EBV, and CMV seropositivity, affected the oral bacterial communities. We achieved this via fractional regressions with the antibody titre (binary: positive vs. negative) as the predictor of interest and the genus abundance as the response variable (Table 1). Due to low case numbers of positive titre results for CMV, we could not adjust for basic covariates without inducing numerical instability in the fitting procedure. While HSV titres did not show any association, a positive EBV titre was positively associated with the abundance of the Neisseria genus in the buccal microbiome (FDR < 0.05, Figure 6A). However, anti-CMV antibody titres showed a very strong association with the OM. In total, 10 genera had a significant association with CMV titres. Nine genera from the Buccal community, eight derived from the Proteobacteria phylum (Acinetobacter, Bradyrhizobium, Comamonadaceae_unclassified, Methylobacterium-Methylorubrum, Oxalobacteraceae_unclassified, and Sphingomonas genera) unveiled a negative association, whereas the genera of Alysiella and Neisseria demonstrated a positive association. One genus from Bacteroidetes phylum, Flavobacterium, also appeared to be negatively associated with positive CMV titres (FDR < 0.05, Figure 6B–J). Two genera of the Salivary community from the Proteobacteria phylum, Pseudomonas and Oxalobaceraceae_unclassified, exhibited a negative association (FDR < 0.05, Figure 7A,B). The full results for the buccal and the saliva microbiome antibody titre associations can be found in Supplementary Table S3. In a further step of sensitivity analysis, we included a positive antibody titre as a covariate into the regression models to investigate the differences in the genus abundances between ELA and controls. However, the results virtually remained the same, indicating either insufficient statistical power to detect potential mediation or that CMV exposure does not mediate ELA-related changes in the OM. This suggests that, unlike increased immunosenescence, the changes we saw in the oral bacterial community are independent of prior exposure to Herpesviridae.

2.4. Fractional Regression Models of the Immune–Microbiome Interactions

In the next step, we fitted a series of fractional regression models integrating the relative abundance of the taxonomic levels in the salivary and buccal compositions with our previously published immune-system profiling. Among the full dataset of 48 immune cell populations, we identified 11 significant associations with genera, most importantly for T cells and NK cells (Table 1).

2.4.1. Association with CD4 T-Cell Immunosenescence

Immunosenescence is a common result of adversity. Thus, we decided to look for possible associations between adversity, microbiome, and accelerated ageing of immune cells (Table 1). For screening the OM associations with the share of CD57 -positive CD4 and CD8 cells, we used multivariable fractional regressions including the genus abundance as the response variable, and the share of CD57-positive CD4 and CD8 cells as a predictor of interest and the basic set of covariates. Additionally, we included the study-group variable as a covariate to control for potential confounding factors related to ELA status. CD8 T cells were previously reported to be significantly associated with CMV [22], but we found no associated taxonomic markers from the OM. From CD4 T-cells tests, we identified six strong taxonomic associations. Two genera from the salivary microbiome, Selenomonas from the Firmicutes phylum showed a positive association and Oxalobacteraceae_unclassified from the Proteobacteria phylum showed a negative association. Four genera from the buccal community: Selenomonas from Firmicutes, Capnocytophaga from the Bacteroidetes phylum, and Campylobacter and Lautropia from the Proteobacteria phylum, displayed a positive association (FDR < 0.05, Figure 8A–F). For further exploration, we fit additional fractional regressions using the number of T-helper cells and T-killer cells as predictors of interest using the same set of covariates as before, finding no additionally significant associations after correction for multiple testing. Summary statistics for the buccal and saliva microbiomes are given in Supplementary Table S4.

2.4.2. Association with NK Cell Activity

Innate immune cells such as natural killer (NK) cells are the first line of defence and often interact with commensal bacteria. Adoptees of this cohort showed increased cytotoxicity on their NK cells [52]; hence, we thought to assess for a potential link with the microbiome. Through screening the OM for associations with various types of NK cells, we found three genera associated with cell counts with an FDR < 0.05, while seven additional associations reached an FDR < 0.1 (Supplementary Table S4), hinting that a better-powered study may find a broader association pattern. In the buccal community, the Oribacterium genus showed a negative association with the total number of NK cells and the total number of mid-maturation NK cells (FDR < 0.05, Figure 9A,B). In parallel, within the salivary community, several genera were significantly associated with different stages of NK maturation. Pseudomonas was found to be positively associated with the total number of CD25 expressing NK cells, which reflects an association with the global activation of NK cells (FDR < 0.05, Figure 9C). The abundance of Alloprevotella was positively associated with the abundance of activated immature CD25CD56hi expressing NK cells (FDR < 0.05, Figure 9D).
Summary statistics for the buccal and saliva microbiome are given in Supplementary Table S4 and bacterial taxa are highlighted in Table 1.

3. Discussion

In this study, we identified taxonomic differences in the OM 24 years after adversity that were common throughout a cohort of diverse cultural and ethnic origins. We identified genera that had a significantly reduced abundance in the adoptees, which were significantly associated with smoking; immunosenescence of CD4 T cells; circulating number and activation status of NK cells; and anti-CMV and, to a lesser degree, anti-EBV titres. Importantly, we were able to see these differences in both the salivary and buccal microbiomes, both of which are readily accessible and both are regularly and easily sampled, even if the buccal microbiome is somewhat underexplored to date. Our data highlight the distinctness of the salivary and buccal microbiomes in distinct oral niches with unique microbial signatures.
Our findings from the EpiPath cohort closely mirror those of Reid et al. [53], although in significantly different microbial communities. We report differences in the abundance of taxa associated with early institutionalisation and CMV seropositivity. Considering that the gut microbiome (GM) is far more labile to lifestyle and environmental impact than the OM [53], our findings build upon those of Ried et al., opening the possibility of much longer-term studies, as the enhanced stability of the OM suggests that differences may be stable over many decades [53]. Expanding our analyses to associations with the immunosenescent CD4 T cells and the activation status of circulating NK cells strengthens the possible role of microbe–immune cross-talk in ELA and the potentially detrimental outcomes. Furthermore, at the family taxonomic level, we observed highly similar differences to those reported by Reid et al. (e.g., Prevotella vs. Alloprevotella, both from the Prevotellaceae family). This highlights the link between the oral and GMs, as numerous studies provide evidence of bacteria migrating from the oral cavity and colonising the gut, whereas there is no evidence of the opposite happening [54,55,56,57].
Our current findings show that institutionalised, genetically unrelated individuals share particular taxa, identifiable 24 years later, independent of the event of adoption. The buccal community, in contrast to the salivary community, appears to be more prone to lifestyle habits such as smoking, agreeing with previous reports that the salivary community remains stable despite lifestyle-hygiene-related mediations such as flossing [42,43,58]. This agrees with several prior reports of the stability and resilience of oral communities over time [39,41,43,46,47,59,60,61,62]. Although host genetics help to shape microbial communities, previous reports of low variance between twins suggest that the shared early environment is the key determinant of the long-term composition [43,60,62]. Longitudinal observations of twins revealed that the salivary microbiome has a stable core community at the genus level, and as twin lives diverge over time, environmental differences increase the diversity between the microbiomes of twins [60,61]. Furthermore, genetically unrelated people with a shared environment show similar environment-related effects on microbiome composition in the mouth as well as other communities [43,47,61]. Cohabitation appears to have a greater impact on the skin microbiome rather than gut and oral communities, persisting after the cohabitation is terminated [43,58,60,63], an effect that is thought to persist for the long term despite leaving or changing household [41,48,49].
The importance of the OM should not be underestimated. As for the GM, there is a direct interaction between the microbiome and both oral and systemic health. Multiple oral inflammatory microbiome-associated conditions such as periodontitis and carries have strong epidemiological and mechanistic associations to other systemic and gastrointestinal diseases [61,64]. Further associations over the years have identified oral marker links to systemic complications, including cardiovascular, immune, metabolic, respiratory, osteopathic, obstetric, and perinatal complications [64,65,66]. In both healthy and inflammatory statuses, viable oral bacteria are often found to travel from the mouth to the gut and are capable of achieving successful colonisation [54,55,56]. Schmidt et al. found that more than half of identified species often found residing in both mouth and gut exhibited signs of oral–gut transmission for all their study participants. Nearly one-third of these are taxa known to be highly dominant in oral communities [54,57,67]. Interestingly, this is a one-way observation: although oral strains can travel to and colonise the gut, the opposite is unlikely to occur [54,57]. Hence, as dental health research has been suggesting for years, oral microbial composition hinges on oral and dental health. In contrast to the prevailing GM, OM shows rising importance as an indicator of systemic health.
Although observational, we report numerous clear associations and correlations in our statistical model that demonstrate the crosstalk between the microbiome and the immune system. Microbial transmission across the gastrointestinal tract, direct microbial contact with tissue-resident innate cells, probable oral bacterial infection, circulating bacterial toxins, and molecular mimicry are all valid candidate pathways that may explain the observed relationship [65,66]. The ELM plays a crucial role in educating immune cells (immune tolerance) that are completely naïve at birth. As immune cells learn to recognise host cells, they are also programmed to recognise antigens from the developing beneficial endogenous microbiome [68]. Tissue-resident dendritic cells harvest microbial antigens from local microbial communities and present them to other immune cells [68]. In germ-free mice, the absence of a microbiome during the early-life period alters immune functions and induces structural defects in lymphoid tissues. In the presence of microbial communities, these tissue structures form normally. Despite many such observations, it is unclear how this acts mechanistically to alter the formation of epithelial barriers. Evidence from the gut suggests that bacteria can direct the glycosylation of luminally exposed surface proteins, a process whose outcome differs in germ-free mice [69,70,71]. Initially, Th17 cells are absent in germ-free mice and only appear upon microbial colonisation [68,69].
It is now well-established that the relationship between stress and chronic disease starts in utero, as susceptibility and occurrence of disease can be predefined by maternal stress [72]. During this period, the naïve, uneducated, immune system develops [73]. NK cells are part of the body’s first line of immune defence, interacting with other immune cells as well as pathogens. In the majority of chronic diseases associated with the early-life environment, NK cells appear to either have an impaired function or an exaggerated cytotoxic activity [74,75]. The most studied NK cell populations are the CD56brightCD16− and CD56dimCD16bright cells and the associated cytotoxic CD56dim and cytokine-producing CD56bright cells [76]. NK cell cytotoxicity is initiated by target cell contact and recognition, which leads to immune synapse formation, resulting in NK-cell-induced target-cell death. The proliferation and expansion of NK cells depend on CD4+ Th1 cells. Nevertheless, due to the bidirectional relationship between innate and adaptive immunity, NK cells impact CD4+ and CD8+ T cells through cytokine production [77]. In the absence of short-chain fatty acids (SCFAs), metabolites produced from fibre fermentation by the local microbiome communities, certain CD4 T-cell subsets do not differentiate. Furthermore, naïve CD8 T cells do not differentiate into memory cells in germ-free mice [68,78,79]. The activation of NK cells by pathogen-associated molecular patterns (PAMPs) may initiate an unwanted response in the microbiome and lead to a strong inflammatory response [80]. Similarly, pathogen-driven activation of NK cells can result in increased on-site cytotoxicity, which can also be harmful to local microbial communities. Correspondingly, microbiome members regulate homeostasis by inducing NK cell expansion and cytokine production or driving the proliferation of anti-inflammatory cytokine-producing NK cells, a common event observed with tissue-resident cells and microbiome crosstalk [81]. Furthermore, immunomodulatory properties of the bacterial community may drive antiviral defences regulating the outcome of viral infection [82].
The OM is intimately linked to oral health. Poor oral health is often approached in an eco-social framework, as it is known to be associated with psychosocial adversity [83]. Both epigenetic and behavioural pathways were linked to poor oral health [83]. One of the most studied causal routes is diet. Affordability and access to a nourishing diet” are strongly influenced by socio-economic status [83], which in turn is linked to the composition of the OM. Detrimental shifts in the microbial composition associated with poor immune responses and mental health were documented for both hospitalised and long-term care home residents [84]. The multidirectional interconnected relationship between the microbial composition, the host’s immunological status, and the resulting life-long health trajectory is most probably highly dependent on constant exposure to particular irritants [84].
Our observation that psychosocial adversity is associated with changes in the OM opens many possibilities for future research. The collection of oral samples, primarily saliva, has been the sampling media of choice for psychobiology, lifestyle, and other social to clinical research areas for many decades. Saliva has long been recognised as an accurate, noninvasive, and cost-effective diagnostic approach that can be tailored to personalised medicine strategies [65,85]. Here, we opened up the possibility of using standard salivary swabs previously collected for microbiome studies. Such studies have the potential to provide a more holistic view of host–microbe interactions and the role of the microbiome in health, which is a potential that can now be applied in nearly all areas of psychobiology (and further afield). Our data also provide preliminary mechanistic insights and the perspectives for future detailed mechanistic studies. We know that early oral microbial colonisation is associated with IL-17-producing cells [86], and subsequent chronic oral disease is often initiated by Th17 cells and IL-17 [87,88,89]. In our EpiPath cohort, there was a strong ELA-associated increase in immunosenescence-associated chronic inflammation, together with increased Th17 cell numbers, although this narrowly missed significance (p = 0.06, [51]). The ELA-associated immunophenotype is centred on immunosenescence [22,51]. Here, we saw clear associations between Selenomonas, Campylobacter, and Capnocytophaga with T-cell immunosenescence, and together with the activated immature NK cell-associated Alloprevotella, these genera were all associated with periodontitis, gingivitis, and T2D. Diseases such as periodontitis an gingivitis have long been associated with changes in both the local and peripheral immune systems. This may be mediated by IL-17 from Th17 cells, and it has been implicated in periodontitis-associated distal diseases in many disease contexts [90,91]. This is directly induced by microbial dysbiosis [92]. Furthermore, direct microbial interaction with immune cells may underlie this, as loss of Toll-like receptor-2 (TLR2) in antigen-presenting cells reduces IL-17 secretion from Th17 cells that dysregulate the host immune system in periodontitis [93]. A similar direct link from the microbiome to the induction of a Th17 cell response was previously reported for Streptococcus [94]. As such, it is interesting to hypothesise that innate immune signalling from TLRs on immune cells within the oral cavity may directly mediate microbiome–immune interactions, acting locally and distally.
As with all investigations, our study is not without limitations. Due to the limited quantity of the biobanked samples, 16S sequencing was favoured over shotgun sequencing to ensure good-quality data, leading to a limited taxonomic resolution in comparison to metagenomics studies. Future metagenomics studies are needed to refine the herein-presented association pattern, exploring potential differences within one genus. The EpiPath cohort consists of only 115 participants. This is a considerable number for a study on ELA, in which a full psychosocial stress test was performed, together with full immune and psychological profiling. However, this sample size is considered small for a microbiome studies, where statistical screening leads to multiple testing, reducing the statistical power for detecting individual associations. Similarly, the reported mediation analyses lack statistical power, and negative results should not be interpreted as the absence of effects. Similarly, as EpiPath is an adoption cohort, metadata such as the mode of birth, if they were ever collected, were never transferred to the adoptive parents. It is also possible that our data could be interpreted as the early inoculation with different microbiomes that simply persisted until 24 years later. The invasive nature of the ELA questions meant that compromise on microbiome-specific metadata, such as dietary habits and oral health status, was unavoidable if maximum participation in the study was to be ensured. Such information would have enhanced the mechanistic potential of our dataset. Knowledge of potential oral complications such as carries or periodontitis will be necessary in future studies to ensure that mechanistic pathways can be explored [95]. As the cohort consists of observational human data, causal interpretations of the reported associations should be treated with care. However, we demonstrated that 16S sequencing, despite its known limitations, provided clear insight into the long-term effect of ELA on the microbiome. Follow-up studies using shotgun metagenomics may refine the reported associations on the species and strain level.

4. Materials and Methods

Participants For this study, we used our previously reported EpiPath cohort of 115 adults aged 20 to 25 years [19,22,51,96]. A total of 75 control participants were brought up by their biological parents and 40 participants were adopted in Luxembourg from institutions worldwide. The median age at adoption was 4.3 months (IQR 0–15 months) [51]. Basic immunoprofiling was available for all cohort members [22,51]. Furthermore, detailed NK cell profiling was available for 76 participants (19 cases and 57 controls), and immunosenescence profiles were available for 79 participants (19 ELA and 60 controls) [22,52]. Biobanked oral swabs were available for 98 participants (33 ELA and 65 controls) and buccal swabs for all 115 participants (40 ELA and 75 controls). For one participant without immunosenescence profiling, the body mass index and sex were missing. This individual was excluded from statistical analyses, where the BMI and/or sex were used as covariates.
Oral samples Saliva samples were collected using Salimetrics Oral Swabs (Salimetrics, Cambridge, UK). Salivary cortisol levels have previously been reported from these samples [19,96]. Buccal swabs were collected with Isoxelix Buccal Swabs (Isohelix, Harrietsham, U.K.). Microbial DNA was extracted using Qiagen DNA from a body fluids kit (Qiagen, Venlo, The Netherlands) according to the manufacturer’s protocol. Samples were quantified with Qubit 1.2 (Invitrogen, Merelbeke, Belgium) and quality was assessed with a Nanodrop (Thermofisher, Merelbeke, Belgium). The V4 region of the 16S gene was amplified from bacterial DNA using 515F [97] and 806R [98] forward and reverse primers (Eurogentec, Seraing, Belgium). The amplification reagents and library preparation were performed using a Quick-16S kit and its equivalent dual indexes (BaseClear, Leiden, The Netherlands) using the manufacturer’s low microbial DNA concentration protocol. Libraries were quantified with Qubit, 1.2, 1.4 (Thermofisher, Merelbeke, Belgium); quality and size were assessed using a BioAnalyser (Agilent, Diegem, Belgium). Sequencing was performed on an Illumina MiSeq system with v2 sequencing chemistry and 500 bp paired-end reads, as well as 10% PHIX control according to the manufacturer’s protocol.
Bioinformatic analyses Fastq files were processed, aligned, and classified using mothur 1.41v [99]. Alpha (inverse Simpson diversity index and Shannon evenness index) and beta diversity (Jaccard Index) were further calculated in the same pipeline. Sequences classification was aligned based on the Silva v138 database [100]. Further integration of microbiome data into the immunophenotype and metadata as well as visualisations were performed with R.
Statistical analyses For descriptive statistics, nominal variables are described by proportions, while metric variables are described by means and standard deviations. Evenness and Shannon entropy metrics were calculated for the OM as measures of alpha diversity and compared between ELA cases and controls with Wilcoxon rank sum tests. Additionally, diversity measures were compared across the OM using rank correlations. For investigating statistical associations between taxonomical units and immune-cell numbers, relative abundances for all genera were checked for outliers. Observations that were outliers both in immune-cell numbers and relative abundances (more than four standard deviations away the mean) were excluded from the analyses, when analysing genus–immune-cell associations. Only genera, or phyla, detected in more than 50% of all cases, were analysed. The microbial abundance data were analysed using fractional regressions [101,102]. Fractional regressions are semiparametric methods not relying on distributional assumptions, and are specifically designed for the analyses of relative abundance data, making them suitable for the analysis of microbiome data, as different species abundances may not be sampled from the same class of distributions. Fractional regressions can be parametrised by odds ratios, allowing for easy interpretation of the regression coefficients in terms of the chance that a certain sequence read is assigned to a taxonomic unit [102]. All fractional regression models, if not specified otherwise, included age, BMI, and sex as covariates, and were performed separately for the OM communities. The basic covariates were included mainly to reduce residual variance and thereby increase statistical power to detect associations with the predictor of interests. Using fractional regressions, we screened the microbiome for associations with the study group variable, basic covariates (age, sex, body mass index (BMI), and smoking), antibody titres for Epstein–Barr virus (EBV), cytomegalovirus (CMV), and the herpes simplex virus (HSV), immunosenescence markers, as well as immune cell counts. All p-values are reported two-tailed. Statistical analyses were performed in STATA 16/MP (College Station, TX, USA), and correction for multiple testing was performed by applying the false discovery rate (FDR) [103]. An FDR < 0.05 was considered to be significant. Summary statistics of the performed analyses are given in Supplementary Tables S1–S4.

5. Conclusions

Our data show a clear link between ELA and the OM that was visible 24 years later. The two oral communities investigated were clearly associated but distinct. We previously reported that ELA induced higher activation and senescence of the immune system. The taxonomic differences in the oral composition were not only associated with ELA but also with the immunosenescence of CD4 T cells, circulating numbers and activation status of NK cells, and anti-CMV titres. Although we do not yet have a detailed mechanistic explanation, our data suggest the presence of multiple links between ELA, immunosenescence, and cytotoxicity that persist through long-term changes in the microbiome.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms222312682/s1.

Author Contributions

Conceptualisation, E.G.C. and C.P.M.; data curation, E.G.C. and J.H.; formal analysis, E.G.C., J.H. and J.D.T.; funding acquisition, C.P.M., I.T. and J.D.T.; investigation, E.G.C., S.B.M., P.G., F.A.D.L. and M.M.C.E.; methodology, E.G.C., S.B.M., P.G. and J.H.; project administration, J.D.T.; supervision, J.D.T.; visualisation, E.G.C.; writing—original draft, E.G.C., J.H. and J.D.T.; writing—review and editing, E.G.C., S.B.M., P.G., C.P.M., F.A.D.L., M.M.C.E., I.T. and J.D.T. All authors have read and agreed to the published version of the manuscript.

Funding

E.G.C. is currently funded by the FNR (PRIDE/11823097/MICROH). The work of J.D.T. on the long-term consequences of ELA was further funded by FNR-CORE (C16/BM/11342695 “MetCOEPs” and C12/BM/3985792 “EpiPath”) and FNR-INTER (INTER/ANR/16/11568350 “MADAM”). J.D.T. is a management committee member of the EU-funded COST action CA18211 focused on early-life and birth-associated trauma. J.H. and I.T. were funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 757922) to I.T. The APC was funded by the Fonds National de Recherche, Luxembourg.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki revised in 2013, and approved by the National Research Ethics Committee of Luxembourg (Comité National d’Ethique de Recherche, CNER, reference 201303/10 v1.4). All participants provided written informed consent, and to compensate for time, effort, and inconvenience, all participants were reimbursed up to €150.

Informed Consent Statement

Written informed consent was obtained from all study participants.

Data Availability Statement

All data from the EpiPath study are available upon reasonable request to J.D.T.

Acknowledgments

The authors would like to thank Paul Wilmes and the members of the FNR-funded doctoral training unit “MicrOH” without whom this study would not have been possible. Furthermore, the authors would like to thank Lorieza Neuberger-Castillo for help with 16s-rRNA sequencing.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Grova, N.; Schroeder, H.; Olivier, J.L.; Turner, J.D. Epigenetic and Neurological Impairments Associated with Early Life Exposure to Persistent Organic Pollutants. Int. J. Genom. 2019, 2019, 2085496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Turner, J.D. Holistic, personalized, immunology? The effects of socioeconomic status on the transcriptional milieu of immune cells. Pediatr. Pulmonol. 2018, 53, 696–697. [Google Scholar] [CrossRef] [PubMed]
  3. Barker, D.J.; Osmond, C. Infant mortality, childhood nutrition, and ischaemic heart disease in England and Wales. Lancet 1986, 1, 1077–1081. [Google Scholar] [CrossRef]
  4. Dube, S.R.; Fairweather, D.; Pearson, W.S.; Felitti, V.J.; Anda, R.F.; Croft, J.B. Cumulative childhood stress and autoimmune diseases in adults. Psychosom. Med. 2009, 71, 243–250. [Google Scholar] [CrossRef] [Green Version]
  5. Wampach, L.; Heintz-Buschart, A.; Fritz, J.V.; Ramiro-Garcia, J.; Habier, J.; Herold, M.; Narayanasamy, S.; Kaysen, A.; Hogan, A.H.; Bindl, L.; et al. Birth mode is associated with earliest strain-conferred gut microbiome functions and immunostimulatory potential. Nat. Commun. 2018, 9, 5091. [Google Scholar] [CrossRef] [Green Version]
  6. Shao, Y.; Forster, S.C.; Tsaliki, E.; Vervier, K.; Strang, A.; Simpson, N.; Kumar, N.; Stares, M.D.; Rodger, A.; Brocklehurst, P.; et al. Stunted microbiota and opportunistic pathogen colonization in caesarean-section birth. Nature 2019, 574, 117–121. [Google Scholar] [CrossRef]
  7. Reyman, M.; van Houten, M.A.; van Baarle, D.; Bosch, A.; Man, W.H.; Chu, M.; Arp, K.; Watson, R.L.; Sanders, E.A.M.; Fuentes, S.; et al. Impact of delivery mode-associated gut microbiota dynamics on health in the first year of life. Nat. Commun. 2019, 10, 4997. [Google Scholar] [CrossRef] [Green Version]
  8. Sarkar, A.; Yoo, J.Y.; Valeria Ozorio Dutra, S.; Morgan, K.H.; Groer, M. The Association between Early-Life Gut Microbiota and Long-Term Health and Diseases. J. Clin. Med. 2021, 10, 459. [Google Scholar] [CrossRef]
  9. Yang, I.; Corwin, E.J.; Brennan, P.A.; Jordan, S.; Murphy, J.R.; Dunlop, A. The Infant Microbiome: Implications for Infant Health and Neurocognitive Development. Nurs. Res. 2016, 65, 76–88. [Google Scholar] [CrossRef] [Green Version]
  10. Moore, R.E.; Townsend, S.D. Temporal development of the infant gut microbiome. Open Biol. 2019, 9, 190128. [Google Scholar] [CrossRef] [Green Version]
  11. Backhed, F.; Roswall, J.; Peng, Y.; Feng, Q.; Jia, H.; Kovatcheva-Datchary, P.; Li, Y.; Xia, Y.; Xie, H.; Zhong, H.; et al. Dynamics and Stabilization of the Human Gut Microbiome during the First Year of Life. Cell Host Microbe 2015, 17, 690–703. [Google Scholar] [CrossRef] [Green Version]
  12. Eriksson, M.; Raikkonen, K.; Eriksson, J.G. Early life stress and later health outcomes—Findings from the Helsinki Birth Cohort Study. Am. J. Hum. Biol. 2014, 26, 111–116. [Google Scholar] [CrossRef]
  13. Spitzer, C.; Wegert, S.; Wollenhaupt, J.; Wingenfeld, K.; Barnow, S.; Grabe, H.J. Gender-specific association between childhood trauma and rheumatoid arthritis: A case-control study. J. Psychosom. Res. 2013, 74, 296–300. [Google Scholar] [CrossRef]
  14. Tomasdottir, M.O.; Sigurdsson, J.A.; Petursson, H.; Kirkengen, A.L.; Krokstad, S.; McEwen, B.; Hetlevik, I.; Getz, L. Self Reported Childhood Difficulties, Adult Multimorbidity and Allostatic Load. A Cross-Sectional Analysis of the Norwegian HUNT Study. PLoS ONE 2015, 10, e0130591. [Google Scholar] [CrossRef] [Green Version]
  15. Gern, J.E.; Visness, C.M.; Gergen, P.J.; Wood, R.A.; Bloomberg, G.R.; O’Connor, G.T.; Kattan, M.; Sampson, H.A.; Witter, F.R.; Sandel, M.T.; et al. The Urban Environment and Childhood Asthma (URECA) birth cohort study: Design, methods, and study population. BMC Pulm. Med. 2009, 9, 17. [Google Scholar] [CrossRef] [Green Version]
  16. Herzog, J.I.; Schmahl, C. Adverse Childhood Experiences and the Consequences on Neurobiological, Psychosocial, and Somatic Conditions across the Lifespan. Front. Psychiatry 2018, 9, 420. [Google Scholar] [CrossRef]
  17. Mansuri, F.; Nash, M.C.; Bakour, C.; Kip, K. Adverse Childhood Experiences (ACEs) and Headaches among Children: A Cross-Sectional Analysis. Headache 2020, 60, 735–744. [Google Scholar] [CrossRef]
  18. Elwenspoek, M.M.C.; Kuehn, A.; Muller, C.P.; Turner, J.D. The effects of early life adversity on the immune system. Psychoneuroendocrinology 2017, 82, 140–154. [Google Scholar] [CrossRef]
  19. Elwenspoek, M.M.C.; Hengesch, X.; Leenen, F.A.D.; Sias, K.; Fernandes, S.B.; Schaan, V.K.; Meriaux, S.B.; Schmitz, S.; Bonnemberger, F.; Schachinger, H.; et al. Glucocorticoid receptor signaling in leukocytes after early life adversity. Dev. Psychopathol. 2020, 32, 853–863. [Google Scholar] [CrossRef]
  20. Holland, J.F.; Khandaker, G.M.; Dauvermann, M.R.; Morris, D.; Zammit, S.; Donohoe, G. Effects of early life adversity on immune function and cognitive performance: Results from the ALSPAC cohort. Soc. Psychiatry Psychiatr. Epidemiol. 2020, 55, 723–733. [Google Scholar] [CrossRef]
  21. Reid, B.M.; Coe, C.L.; Doyle, C.M.; Sheerar, D.; Slukvina, A.; Donzella, B.; Gunnar, M.R. Persistent skewing of the T-cell profile in adolescents adopted internationally from institutional care. Brain Behav. Immun. 2019, 77, 168–177. [Google Scholar] [CrossRef]
  22. Elwenspoek, M.M.C.; Sias, K.; Hengesch, X.; Schaan, V.K.; Leenen, F.A.D.; Adams, P.; Meriaux, S.B.; Schmitz, S.; Bonnemberger, F.; Ewen, A.; et al. T Cell Immunosenescence after Early Life Adversity: Association with Cytomegalovirus Infection. Front. Immunol. 2017, 8, 1263. [Google Scholar] [CrossRef] [Green Version]
  23. Nielsen, C.M.; White, M.J.; Goodier, M.R.; Riley, E.M. Functional Significance of CD57 Expression on Human NK Cells and Relevance to Disease. Front. Immunol. 2013, 4, 422. [Google Scholar] [CrossRef] [Green Version]
  24. Judge, S.J.; Murphy, W.J.; Canter, R.J. Characterizing the Dysfunctional NK Cell: Assessing the Clinical Relevance of Exhaustion, Anergy, and Senescence. Front. Cell Infect. Microbiol. 2020, 10, 49. [Google Scholar] [CrossRef] [Green Version]
  25. Della Chiesa, M.; Falco, M.; Podesta, M.; Locatelli, F.; Moretta, L.; Frassoni, F.; Moretta, A. Phenotypic and functional heterogeneity of human NK cells developing after umbilical cord blood transplantation: A role for human cytomegalovirus? Blood 2012, 119, 399–410. [Google Scholar] [CrossRef]
  26. Lopez-Verges, S.; Milush, J.M.; Schwartz, B.S.; Pando, M.J.; Jarjoura, J.; York, V.A.; Houchins, J.P.; Miller, S.; Kang, S.M.; Norris, P.J.; et al. Expansion of a unique CD57(+)NKG2Chi natural killer cell subset during acute human cytomegalovirus infection. Proc. Natl. Acad. Sci. USA 2011, 108, 14725–14732. [Google Scholar] [CrossRef] [Green Version]
  27. Yang, X.; Xie, L.; Li, Y.; Wei, C. More than 9,000,000 unique genes in human gut bacterial community: Estimating gene numbers inside a human body. PLoS ONE 2009, 4, e6074. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, L.; Alammar, N.; Singh, R.; Nanavati, J.; Song, Y.; Chaudhary, R.; Mullin, G.E. Gut Microbial Dysbiosis in the Irritable Bowel Syndrome: A Systematic Review and Meta-Analysis of Case-Control Studies. J. Acad. Nutr. Diet. 2020, 120, 565–586. [Google Scholar] [CrossRef] [Green Version]
  29. Rogers, G.B.; Keating, D.J.; Young, R.L.; Wong, M.L.; Licinio, J.; Wesselingh, S. From gut dysbiosis to altered brain function and mental illness: Mechanisms and pathways. Mol. Psychiatry 2016, 21, 738–748. [Google Scholar] [CrossRef] [Green Version]
  30. Miller, G.E.; Chen, E.; Shalowitz, M.U.; Story, R.E.; Leigh, A.K.K.; Ham, P.; Arevalo, J.M.G.; Cole, S.W. Divergent transcriptional profiles in pediatric asthma patients of low and high socioeconomic status. Pediatr. Pulmonol. 2018, 53, 710–719. [Google Scholar] [CrossRef]
  31. Zijlmans, M.A.; Korpela, K.; Riksen-Walraven, J.M.; de Vos, W.M.; de Weerth, C. Maternal prenatal stress is associated with the infant intestinal microbiota. Psychoneuroendocrinology 2015, 53, 233–245. [Google Scholar] [CrossRef] [PubMed]
  32. Alcon-Giner, C.; Dalby, M.J.; Caim, S.; Ketskemety, J.; Shaw, A.; Sim, K.; Lawson, M.A.E.; Kiu, R.; Leclaire, C.; Chalklen, L.; et al. Microbiota Supplementation with Bifidobacterium and Lactobacillus Modifies the Preterm Infant Gut Microbiota and Metabolome: An Observational Study. Cell Rep. Med. 2020, 1, 100077. [Google Scholar] [CrossRef] [PubMed]
  33. Stewart, C.J.; Ajami, N.J.; O’Brien, J.L.; Hutchinson, D.S.; Smith, D.P.; Wong, M.C.; Ross, M.C.; Lloyd, R.E.; Doddapaneni, H.; Metcalf, G.A.; et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature 2018, 562, 583–588. [Google Scholar] [CrossRef] [PubMed]
  34. Xu, M.; Wang, C.; Krolick, K.N.; Shi, H.; Zhu, J. Difference in post-stress recovery of the gut microbiome and its altered metabolism after chronic adolescent stress in rats. Sci. Rep. 2020, 10, 3950. [Google Scholar] [CrossRef] [Green Version]
  35. Toussaint, L.; Shields, G.S.; Dorn, G.; Slavich, G.M. Effects of lifetime stress exposure on mental and physical health in young adulthood: How stress degrades and forgiveness protects health. J. Health Psychol. 2016, 21, 1004–1014. [Google Scholar] [CrossRef] [Green Version]
  36. Duran-Pinedo, A.E.; Solbiati, J.; Frias-Lopez, J. The effect of the stress hormone cortisol on the metatranscriptome of the oral microbiome. NPJ Biofilms Microbiomes 2018, 4, 25. [Google Scholar] [CrossRef] [Green Version]
  37. Boustedt, K.; Roswall, J.; Dahlen, G.; Dahlgren, J.; Twetman, S. Salivary microflora and mode of delivery: A prospective case control study. BMC Oral Health 2015, 15, 155. [Google Scholar] [CrossRef] [Green Version]
  38. Chu, D.M.; Ma, J.; Prince, A.L.; Antony, K.M.; Seferovic, M.D.; Aagaard, K.M. Maturation of the infant microbiome community structure and function across multiple body sites and in relation to mode of delivery. Nat. Med. 2017, 23, 314–326. [Google Scholar] [CrossRef] [Green Version]
  39. Mason, M.R.; Chambers, S.; Dabdoub, S.M.; Thikkurissy, S.; Kumar, P.S. Characterizing oral microbial communities across dentition states and colonization niches. Microbiome 2018, 6, 67. [Google Scholar] [CrossRef]
  40. Mark Welch, J.L.; Dewhirst, F.E.; Borisy, G.G. Biogeography of the Oral Microbiome: The Site-Specialist Hypothesis. Annu. Rev. Microbiol. 2019, 73, 335–358. [Google Scholar] [CrossRef]
  41. Carpenter, G.H. Salivary Factors that Maintain the Normal Oral Commensal Microflora. J. Dent. Res. 2020, 99, 644–649. [Google Scholar] [CrossRef]
  42. David, L.A.; Materna, A.C.; Friedman, J.; Campos-Baptista, M.I.; Blackburn, M.C.; Perrotta, A.; Erdman, S.E.; Alm, E.J. Host lifestyle affects human microbiota on daily timescales. Genome Biol. 2014, 15, R89. [Google Scholar] [CrossRef] [Green Version]
  43. Shaw, L.; Ribeiro, A.L.R.; Levine, A.P.; Pontikos, N.; Balloux, F.; Segal, A.W.; Roberts, A.P.; Smith, A.M. The Human Salivary Microbiome Is Shaped by Shared Environment Rather than Genetics: Evidence from a Large Family of Closely Related Individuals. mBio 2017, 8, e01237-17. [Google Scholar] [CrossRef] [Green Version]
  44. Marsh, P.D. Dental plaque as a biofilm and a microbial community—Implications for health and disease. BMC Oral Health 2006, 6 (Suppl. S1), S14. [Google Scholar] [CrossRef] [Green Version]
  45. Zaura, E.; Brandt, B.W.; Teixeira de Mattos, M.J.; Buijs, M.J.; Caspers, M.P.; Rashid, M.U.; Weintraub, A.; Nord, C.E.; Savell, A.; Hu, Y.; et al. Same Exposure but Two Radically Different Responses to Antibiotics: Resilience of the Salivary Microbiome versus Long-Term Microbial Shifts in Feces. mBio 2015, 6, e01693-15. [Google Scholar] [CrossRef] [Green Version]
  46. Kennedy, B.; Peura, S.; Hammar, U.; Vicenzi, S.; Hedman, A.; Almqvist, C.; Andolf, E.; Pershagen, G.; Dicksved, J.; Bertilsson, S.; et al. Oral Microbiota Development in Early Childhood. Sci. Rep. 2019, 9, 19025. [Google Scholar] [CrossRef]
  47. Abeles, S.R.; Jones, M.B.; Santiago-Rodriguez, T.M.; Ly, M.; Klitgord, N.; Yooseph, S.; Nelson, K.E.; Pride, D.T. Microbial diversity in individuals and their household contacts following typical antibiotic courses. Microbiome 2016, 4, 39. [Google Scholar] [CrossRef] [Green Version]
  48. Almeida-Santos, A.; Martins-Mendes, D.; Gaya-Vidal, M.; Perez-Pardal, L.; Beja-Pereira, A. Characterization of the Oral Microbiome of Medicated Type-2 Diabetes Patients. Front. Microbiol. 2021, 12, 610370. [Google Scholar] [CrossRef]
  49. Mukherjee, C.; Moyer, C.O.; Steinkamp, H.M.; Hashmi, S.B.; Beall, C.J.; Guo, X.; Ni, A.; Leys, E.J.; Griffen, A.L. Acquisition of oral microbiota is driven by environment, not host genetics. Microbiome 2021, 9, 54. [Google Scholar] [CrossRef]
  50. Jakubovics, N.S. Saliva as the Sole Nutritional Source in the Development of Multispecies Communities in Dental Plaque. Microbiol. Spectr. 2015, 3, 263–277. [Google Scholar] [CrossRef] [Green Version]
  51. Elwenspoek, M.M.C.; Hengesch, X.; Leenen, F.A.D.; Schritz, A.; Sias, K.; Schaan, V.K.; Meriaux, S.B.; Schmitz, S.; Bonnemberger, F.; Schachinger, H.; et al. Proinflammatory T Cell Status Associated with Early Life Adversity. J. Immunol. 2017, 199, 4046–4055. [Google Scholar] [CrossRef] [Green Version]
  52. Fernandes, S.B.; Patil, N.D.; Meriaux, S.B.; Theresine, M.; Leenen, F.A.D.; Elwenspoek, M.M.C.; Zimmer, J.; Turner, J.D. Unbiased Screening Identifies Functional Differences in NK Cells after Early Life Psycho-Social Stress. Front. Immunol. 2021, 12, 674532. [Google Scholar] [CrossRef]
  53. Reid, B.M.; Horne, R.; Donzella, B.; Szamosi, J.C.; Coe, C.L.; Foster, J.A.; Gunnar, M.R. Microbiota-immune alterations in adolescents following early life adversity: A proof of concept study. Dev. Psychobiol. 2021, 63, 851–863. [Google Scholar] [CrossRef]
  54. Prodan, A.; Levin, E.; Nieuwdorp, M. Does disease start in the mouth, the gut or both? eLife 2019, 8, e45931. [Google Scholar] [CrossRef]
  55. Valdes, A.M.; Walter, J.; Segal, E.; Spector, T.D. Role of the gut microbiota in nutrition and health. BMJ 2018, 361, k2179. [Google Scholar] [CrossRef] [Green Version]
  56. Ridlon, J.M.; Kang, D.J.; Hylemon, P.B.; Bajaj, J.S. Bile acids and the gut microbiome. Curr. Opin. Gastroenterol. 2014, 30, 332–338. [Google Scholar] [CrossRef] [Green Version]
  57. Schmidt, T.S.; Hayward, M.R.; Coelho, L.P.; Li, S.S.; Costea, P.I.; Voigt, A.Y.; Wirbel, J.; Maistrenko, O.M.; Alves, R.J.; Bergsten, E.; et al. Extensive transmission of microbes along the gastrointestinal tract. Elife 2019, 8, e42693. [Google Scholar] [CrossRef]
  58. Utter, D.R.; Mark Welch, J.L.; Borisy, G.G. Individuality, Stability, and Variability of the Plaque Microbiome. Front. Microbiol. 2016, 7, 564. [Google Scholar] [CrossRef]
  59. Belstrom, D.; Holmstrup, P.; Bardow, A.; Kokaras, A.; Fiehn, N.E.; Paster, B.J. Temporal Stability of the Salivary Microbiota in Oral Health. PLoS ONE 2016, 11, e0147472. [Google Scholar] [CrossRef]
  60. Stahringer, S.S.; Clemente, J.C.; Corley, R.P.; Hewitt, J.; Knights, D.; Walters, W.A.; Knight, R.; Krauter, K.S. Nurture trumps nature in a longitudinal survey of salivary bacterial communities in twins from early adolescence to early adulthood. Genome Res. 2012, 22, 2146–2152. [Google Scholar] [CrossRef] [Green Version]
  61. Gomez, A.; Espinoza, J.L.; Harkins, D.M.; Leong, P.; Saffery, R.; Bockmann, M.; Torralba, M.; Kuelbs, C.; Kodukula, R.; Inman, J.; et al. Host Genetic Control of the Oral Microbiome in Health and Disease. Cell Host Microbe 2017, 22, 269–278.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Premaraj, T.S.; Vella, R.; Chung, J.; Lin, Q.; Panier, H.; Underwood, K.; Premaraj, S.; Zhou, Y. Ethnic variation of oral microbiota in children. Sci. Rep. 2020, 10, 14788. [Google Scholar] [CrossRef] [PubMed]
  63. Blekhman, R.; Goodrich, J.K.; Huang, K.; Sun, Q.; Bukowski, R.; Bell, J.T.; Spector, T.D.; Keinan, A.; Ley, R.E.; Gevers, D.; et al. Host genetic variation impacts microbiome composition across human body sites. Genome Biol. 2015, 16, 191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Byrd, K.M.; Gulati, A.S. The “Gum-Gut” Axis in Inflammatory Bowel Diseases: A Hypothesis-Driven Review of Associations and Advances. Front. Immunol. 2021, 12, 620124. [Google Scholar] [CrossRef]
  65. Gomez, A.; Nelson, K.E. The Oral Microbiome of Children: Development, Disease, and Implications beyond Oral Health. Microb. Ecol. 2017, 73, 492–503. [Google Scholar] [CrossRef] [Green Version]
  66. Li, X.; Kolltveit, K.M.; Tronstad, L.; Olsen, I. Systemic diseases caused by oral infection. Clin. Microbiol. Rev. 2000, 13, 547–558. [Google Scholar] [CrossRef]
  67. Li, S.S.; Zhu, A.; Benes, V.; Costea, P.I.; Hercog, R.; Hildebrand, F.; Huerta-Cepas, J.; Nieuwdorp, M.; Salojarvi, J.; Voigt, A.Y.; et al. Durable coexistence of donor and recipient strains after fecal microbiota transplantation. Science 2016, 352, 586–589. [Google Scholar] [CrossRef]
  68. Zheng, D.; Liwinski, T.; Elinav, E. Interaction between microbiota and immunity in health and disease. Cell Res. 2020, 30, 492–506. [Google Scholar] [CrossRef]
  69. Round, J.L.; Mazmanian, S.K. The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 2009, 9, 313–323. [Google Scholar] [CrossRef]
  70. Bouskra, D.; Brezillon, C.; Berard, M.; Werts, C.; Varona, R.; Boneca, I.G.; Eberl, G. Lymphoid tissue genesis induced by commensals through NOD1 regulates intestinal homeostasis. Nature 2008, 456, 507–510. [Google Scholar] [CrossRef]
  71. Bry, L.; Falk, P.G.; Midtvedt, T.; Gordon, J.I. A model of host-microbial interactions in an open mammalian ecosystem. Science 1996, 273, 1380–1383. [Google Scholar] [CrossRef]
  72. Henriksen, R.E.; Thuen, F. Marital Quality and Stress in Pregnancy Predict the Risk of Infectious Disease in the Offspring: The Norwegian Mother and Child Cohort Study. PLoS ONE 2015, 10, e0137304. [Google Scholar] [CrossRef] [Green Version]
  73. Fragkou, P.C.; Karaviti, D.; Zemlin, M.; Skevaki, C. Impact of Early Life Nutrition on Children’s Immune System and Noncommunicable Diseases Through Its Effects on the Bacterial Microbiome, Virome and Mycobiome. Front. Immunol. 2021, 12, 644269. [Google Scholar] [CrossRef]
  74. Ong, S.; Rose, N.R.; Cihakova, D. Natural killer cells in inflammatory heart disease. Clin. Immunol. 2017, 175, 26–33. [Google Scholar] [CrossRef] [Green Version]
  75. Yang, Y.; Day, J.; Souza-Fonseca Guimaraes, F.; Wicks, I.P.; Louis, C. Natural killer cells in inflammatory autoimmune diseases. Clin. Transl. Immunol. 2021, 10, e1250. [Google Scholar] [CrossRef]
  76. Poli, A.; Michel, T.; Patil, N.; Zimmer, J. Revisiting the Functional Impact of NK Cells. Trends Immunol. 2018, 39, 460–472. [Google Scholar] [CrossRef]
  77. Abel, A.M.; Yang, C.; Thakar, M.S.; Malarkannan, S. Natural Killer Cells: Development, Maturation, and Clinical Utilization. Front. Immunol. 2018, 9, 1869. [Google Scholar] [CrossRef] [Green Version]
  78. Smith, P.M.; Howitt, M.R.; Panikov, N.; Michaud, M.; Gallini, C.A.; Bohlooly, Y.M.; Glickman, J.N.; Garrett, W.S. The microbial metabolites, short-chain fatty acids, regulate colonic Treg cell homeostasis. Science 2013, 341, 569–573. [Google Scholar] [CrossRef] [Green Version]
  79. Bachem, A.; Makhlouf, C.; Binger, K.J.; de Souza, D.P.; Tull, D.; Hochheiser, K.; Whitney, P.G.; Fernandez-Ruiz, D.; Dahling, S.; Kastenmuller, W.; et al. Microbiota-Derived Short-Chain Fatty Acids Promote the Memory Potential of Antigen-Activated CD8(+) T Cells. Immunity 2019, 51, 285–297. [Google Scholar] [CrossRef]
  80. Souza-Fonseca-Guimaraes, F.; Adib-Conquy, M.; Cavaillon, J.M. Natural killer (NK) cells in antibacterial innate immunity: Angels or devils? Mol. Med. 2012, 18, 270–285. [Google Scholar] [CrossRef]
  81. Theresine, M.; Patil, N.D.; Zimmer, J. Airway Natural Killer Cells and Bacteria in Health and Disease. Front. Immunol. 2020, 11, 585048. [Google Scholar] [CrossRef]
  82. Przemska-Kosicka, A.; Childs, C.E.; Maidens, C.; Dong, H.; Todd, S.; Gosney, M.A.; Tuohy, K.M.; Yaqoob, P. Age-Related Changes in the Natural Killer Cell Response to Seasonal Influenza Vaccination Are Not Influenced by a Synbiotic: A Randomised Controlled Trial. Front. Immunol. 2018, 9, 591. [Google Scholar] [CrossRef]
  83. Lee, J.Y.; Divaris, K. The ethical imperative of addressing oral health disparities: A unifying framework. J. Dent. Res. 2014, 93, 224–230. [Google Scholar] [CrossRef] [Green Version]
  84. Coman, V.; Vodnar, D.C. Gut microbiota and old age: Modulating factors and interventions for healthy longevity. Exp. Gerontol. 2020, 141, 111095. [Google Scholar] [CrossRef]
  85. Verma, D.; Garg, P.K.; Dubey, A.K. Insights into the human oral microbiome. Arch. Microbiol. 2018, 200, 525–540. [Google Scholar] [CrossRef] [PubMed]
  86. Koren, N.; Zubeidat, K.; Saba, Y.; Horev, Y.; Barel, O.; Wilharm, A.; Heyman, O.; Wald, S.; Eli-Berchoer, L.; Shapiro, H.; et al. Maturation of the neonatal oral mucosa involves unique epithelium-microbiota interactions. Cell Host Microbe 2021, 29, 197–209.e5. [Google Scholar] [CrossRef] [PubMed]
  87. Bellando-Randone, S.; Russo, E.; Venerito, V.; Matucci-Cerinic, M.; Iannone, F.; Tangaro, S.; Amedei, A. Exploring the Oral Microbiome in Rheumatic Diseases, State of Art and Future Prospective in Personalized Medicine with an AI Approach. J. Pers. Med. 2021, 11, 625. [Google Scholar] [CrossRef] [PubMed]
  88. Abusleme, L.; Moutsopoulos, N.M. IL-17: Overview and role in oral immunity and microbiome. Oral Dis. 2017, 23, 854–865. [Google Scholar] [CrossRef] [PubMed]
  89. Gaffen, S.L.; Moutsopoulos, N.M. Regulation of host-microbe interactions at oral mucosal barriers by type 17 immunity. Sci. Immunol. 2020, 5, eaau4594. [Google Scholar] [CrossRef]
  90. Konkel, J.E.; O’Boyle, C.; Krishnan, S. Distal Consequences of Oral Inflammation. Front. Immunol. 2019, 10, 1403. [Google Scholar] [CrossRef] [Green Version]
  91. De Aquino, S.G.; Talbot, J.; Sonego, F.; Turato, W.M.; Grespan, R.; Avila-Campos, M.J.; Cunha, F.Q.; Cirelli, J.A. The aggravation of arthritis by periodontitis is dependent of IL-17 receptor A activation. J. Clin. Periodontol. 2017, 44, 881–891. [Google Scholar] [CrossRef]
  92. Dutzan, N.; Kajikawa, T.; Abusleme, L.; Greenwell-Wild, T.; Zuazo, C.E.; Ikeuchi, T.; Brenchley, L.; Abe, T.; Hurabielle, C.; Martin, D.; et al. A dysbiotic microbiome triggers TH17 cells to mediate oral mucosal immunopathology in mice and humans. Sci. Transl. Med. 2018, 10, eaat0797. [Google Scholar] [CrossRef] [Green Version]
  93. De Aquino, S.G.; Abdollahi-Roodsaz, S.; Koenders, M.I.; van de Loo, F.A.; Pruijn, G.J.; Marijnissen, R.J.; Walgreen, B.; Helsen, M.M.; van den Bersselaar, L.A.; de Molon, R.S.; et al. Periodontal pathogens directly promote autoimmune experimental arthritis by inducing a TLR2- and IL-1-driven Th17 response. J. Immunol. 2014, 192, 4103–4111. [Google Scholar] [CrossRef] [Green Version]
  94. Goncalves, M.T.; Mitchell, T.J.; Lord, J.M. Immune ageing and susceptibility to Streptococcus pneumoniae. Biogerontology 2016, 17, 449–465. [Google Scholar] [CrossRef]
  95. Dimitrov, D.V.; Hoeng, J. Systems approaches to computational modeling of the oral microbiome. Front. Physiol. 2013, 4, 172. [Google Scholar] [CrossRef] [Green Version]
  96. Hengesch, X.; Elwenspoek, M.M.C.; Schaan, V.K.; Larra, M.F.; Finke, J.B.; Zhang, X.; Bachmann, P.; Turner, J.D.; Vogele, C.; Muller, C.P.; et al. Blunted endocrine response to a combined physical-cognitive stressor in adults with early life adversity. Child. Abuse Negl. 2018, 85, 137–144. [Google Scholar] [CrossRef]
  97. Parada, A.E.; Needham, D.M.; Fuhrman, J.A. Every base matters: Assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 2016, 18, 1403–1414. [Google Scholar] [CrossRef]
  98. Apprill, A.; McNally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microbial. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef] [Green Version]
  99. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  100. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glockner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  101. Papke, L.E.; Wooldridge, J.M. Econometric methods for fractional response variables with an application to 401 (k) plan participation rates. J. Appl. Econom. 1996, 11, 619–632. [Google Scholar] [CrossRef] [Green Version]
  102. Baldini, F.; Hertel, J.; Sandt, E.; Thinnes, C.C.; Neuberger-Castillo, L.; Pavelka, L.; Betsou, F.; Krüger, R.; Thiele, I. Parkinson’s disease-associated alterations of the gut microbiome predict disease-relevant changes in metabolic functions. BMC Biol. 2020, 18, 62. [Google Scholar] [CrossRef]
  103. Benjamini, Y. Discovering the false discovery rate. J. R. Stat. Soc. Ser. B 2010, 72, 405–416. [Google Scholar] [CrossRef]
Ijms 22 12682 g001 550
Figure 1. Overall composition of salivary bacterial community. (A) Overall microbial composition displayed in stacked area bar plot with the percentage relative abundance of all phyla found in each participant in both study arms. (B) Top 20 most abundant genera by mean abundance arranged graphically by phyla. Vertical line = mean; rectangle = 1st to 3rd quartile; horizontal lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points. Blue, control group; red, ELA group.
Figure 1. Overall composition of salivary bacterial community. (A) Overall microbial composition displayed in stacked area bar plot with the percentage relative abundance of all phyla found in each participant in both study arms. (B) Top 20 most abundant genera by mean abundance arranged graphically by phyla. Vertical line = mean; rectangle = 1st to 3rd quartile; horizontal lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points. Blue, control group; red, ELA group.
Ijms 22 12682 g001
Ijms 22 12682 g002 550
Figure 2. Overall composition of the buccal bacterial community. (A) Overall microbial composition displayed in stacked area bar plot with the percentage relative abundance of all phyla found for each participant in both study arms. (B) Top 20 most abundant genera by mean abundance arranged graphically by phyla. Vertical line = mean; rectangle = 1st to 3rd quartile; horizontal lines = 2.5th to 97.5the percentile. Outliers are indicated as individual data points. Blue, control group; red, ELA group.
Figure 2. Overall composition of the buccal bacterial community. (A) Overall microbial composition displayed in stacked area bar plot with the percentage relative abundance of all phyla found for each participant in both study arms. (B) Top 20 most abundant genera by mean abundance arranged graphically by phyla. Vertical line = mean; rectangle = 1st to 3rd quartile; horizontal lines = 2.5th to 97.5the percentile. Outliers are indicated as individual data points. Blue, control group; red, ELA group.
Ijms 22 12682 g002
Ijms 22 12682 g003 550
Figure 3. Diversity and evenness of the salivary and buccal bacterial communities in both study groups. (A) Inverse Simpson diversity index of saliva against buccal communities for control (left) and ELA (right). (B) Shannon evenness index of saliva against buccal communities for control (left) and ELA (right). No correlation was found between either community for either measure or group (Spearman’s rank correlation test, p > 0.47). Grey shaded area: 95% confidence interval.
Figure 3. Diversity and evenness of the salivary and buccal bacterial communities in both study groups. (A) Inverse Simpson diversity index of saliva against buccal communities for control (left) and ELA (right). (B) Shannon evenness index of saliva against buccal communities for control (left) and ELA (right). No correlation was found between either community for either measure or group (Spearman’s rank correlation test, p > 0.47). Grey shaded area: 95% confidence interval.
Ijms 22 12682 g003
Ijms 22 12682 g004 550
Figure 4. Taxonomic differences between study groups in the salivary bacterial community. Box plots of two phyla, (A) Proteobacteria and (B) Verruomicrobiota, as well as four genera: (C) Aquabacterium, (D) Comamonas, (E) Comamonadaceae, unclassified and (F) Acinetobacter. All are significantly associated with study group; fractional regressions against study group were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Figure 4. Taxonomic differences between study groups in the salivary bacterial community. Box plots of two phyla, (A) Proteobacteria and (B) Verruomicrobiota, as well as four genera: (C) Aquabacterium, (D) Comamonas, (E) Comamonadaceae, unclassified and (F) Acinetobacter. All are significantly associated with study group; fractional regressions against study group were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Ijms 22 12682 g004
Ijms 22 12682 g005 550
Figure 5. Taxonomic differences in the buccal bacterial community associated with smoking. Box plots of five genera: (A) Pasteurellaceae unclassified, (B) Neisseria, (C) Capnocytophaga, (D) Neisseriaceae unclassified, and (E) Planococcaceae unclassified. All are significantly associated with study group; fractional regressions against smoking were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Figure 5. Taxonomic differences in the buccal bacterial community associated with smoking. Box plots of five genera: (A) Pasteurellaceae unclassified, (B) Neisseria, (C) Capnocytophaga, (D) Neisseriaceae unclassified, and (E) Planococcaceae unclassified. All are significantly associated with study group; fractional regressions against smoking were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Ijms 22 12682 g005
Ijms 22 12682 g006 550
Figure 6. Taxonomic differences in the buccal bacterial community associated with anti-herpesviridiae serological status. Box plots of one genus (A) Neisseria significantly associated with anti-EBV antibody titres. Nine genera, (B) Neisseria, (C) Sphingomonas, (D) Acinetobacter, (E) Oxalobacteraceae unclassified, (F) Bradyrhizobium, (G) Flavobacterium, (H) Methylorubrum, (I) Comamonadeceae unclassified, and (J) Alysiella, were significantly associated with anti-CMV antibody titres. Fractional regressions against the presence of anti-EBV and anti-CMV antibodies were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Figure 6. Taxonomic differences in the buccal bacterial community associated with anti-herpesviridiae serological status. Box plots of one genus (A) Neisseria significantly associated with anti-EBV antibody titres. Nine genera, (B) Neisseria, (C) Sphingomonas, (D) Acinetobacter, (E) Oxalobacteraceae unclassified, (F) Bradyrhizobium, (G) Flavobacterium, (H) Methylorubrum, (I) Comamonadeceae unclassified, and (J) Alysiella, were significantly associated with anti-CMV antibody titres. Fractional regressions against the presence of anti-EBV and anti-CMV antibodies were calculated to determine significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Ijms 22 12682 g006
Ijms 22 12682 g007 550
Figure 7. Taxonomic differences in the salivary bacterial community associated with anti-CMV serological status. Box plots of two genera, (A) Pseudomonas and (B) Oxalobacteraceae unclassified, which were significantly associated with CMV antibody titres. Fractional regressions against the presence of anti-CMV antibodies were calculated to determined significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Figure 7. Taxonomic differences in the salivary bacterial community associated with anti-CMV serological status. Box plots of two genera, (A) Pseudomonas and (B) Oxalobacteraceae unclassified, which were significantly associated with CMV antibody titres. Fractional regressions against the presence of anti-CMV antibodies were calculated to determined significance (FDR < 0.05). Horizontal line = mean; rectangle = 1st to 3rd quartile; vertical lines = 2.5th to 97.5th percentile. Outliers are indicated as individual data points.
Ijms 22 12682 g007
Ijms 22 12682 g008 550
Figure 8. Taxonomic associations in both communities with immunosenescence. Scatter plots with regression lines for six genera. From the salivary community, two genera, (A) Oxalobacteraceae unclassified and (B) Selenomonas, as well as four genera from the buccal community, (C) Selenomonas, (D) Capnocytophaga, (E) Campylobacter, and (F) Lautropia, were all significantly associated with CD4 CD57 cell counts. Fractional regressions against CD4 CD57 cell counts were calculated to determine significance (FDR < 0.05). Regression lines were derived from fractional regressions with logistic parametrisation of the conditional mean.
Figure 8. Taxonomic associations in both communities with immunosenescence. Scatter plots with regression lines for six genera. From the salivary community, two genera, (A) Oxalobacteraceae unclassified and (B) Selenomonas, as well as four genera from the buccal community, (C) Selenomonas, (D) Capnocytophaga, (E) Campylobacter, and (F) Lautropia, were all significantly associated with CD4 CD57 cell counts. Fractional regressions against CD4 CD57 cell counts were calculated to determine significance (FDR < 0.05). Regression lines were derived from fractional regressions with logistic parametrisation of the conditional mean.
Ijms 22 12682 g008
Ijms 22 12682 g009 550
Figure 9. Taxonomic associations in both communities with natural killer cell activity. Scatter plots with regression lines for one genus from the buccal community, Oribacterium, was significantly associated with (A) the total number of NK cells and (B) the total number of mid-maturation NK cells as well as two genera from the salivary community: (C) Pseudomonas, associated with CD25+ NK cell counts; and (D) Alloprevotella, associated with CD25CD56hi NK cell counts. Fractional regressions against NK cell counts were calculated to determine significance (FDR < 0.05). Regression lines were derived from fractional regressions with logistic parametrisation of the conditional mean.
Figure 9. Taxonomic associations in both communities with natural killer cell activity. Scatter plots with regression lines for one genus from the buccal community, Oribacterium, was significantly associated with (A) the total number of NK cells and (B) the total number of mid-maturation NK cells as well as two genera from the salivary community: (C) Pseudomonas, associated with CD25+ NK cell counts; and (D) Alloprevotella, associated with CD25CD56hi NK cell counts. Fractional regressions against NK cell counts were calculated to determine significance (FDR < 0.05). Regression lines were derived from fractional regressions with logistic parametrisation of the conditional mean.
Ijms 22 12682 g009
Table
Table 1. An overview of the bacterial taxa associated with all the tested covariates in both saliva and buccal microbiomes.
Table 1. An overview of the bacterial taxa associated with all the tested covariates in both saliva and buccal microbiomes.
SalivaBuccal
Number of
Associations (FDR < 0.05)		Positively
Associated Taxa		Negatively
Associated Taxa		Number of
Associations (FDR < 0.05)		Positively
Associated Taxa		Negatively
Associated Taxa
Smoking0--5-Pasteurellaceae (unclassified), Neisseria, Capnocytophaga, Neisseriaceae (unc 1), Planococcaceae (unc)
ELA vs. control4-Aquabacterium, Comamonas, Comamonadaceae (unc), Acinetobacter0--
Anti-CMV seropositive2-Pseudomonas, Oxalobaceraceae (unc)9Alysiella, NeisseriaSphingomonas, Acinetobacter, Oxalobacteraceae (unc), Bradyrhizobium, Flavobacterium, Methylorubrum, Comamonadaceae (unc)
Anti-EBV seropositive0--1Neisseria-
HSV0--0--
CD4+
CD57+		2SelenomonasOxalobaceraceae (unc)4Selenomonas, Capnocytophaga, Campylobacter, Lautropia-
CD8+
CD57+		0--0--
Total CTLs0--0--
Total Th cells0--0--
1 unc = unclassified.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Charalambous, E.G.; Mériaux, S.B.; Guebels, P.; Muller, C.P.; Leenen, F.A.D.; Elwenspoek, M.M.C.; Thiele, I.; Hertel, J.; Turner, J.D. Early-Life Adversity Leaves Its Imprint on the Oral Microbiome for More Than 20 Years and Is Associated with Long-Term Immune Changes. Int. J. Mol. Sci. 2021, 22, 12682. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222312682

AMA Style

Charalambous EG, Mériaux SB, Guebels P, Muller CP, Leenen FAD, Elwenspoek MMC, Thiele I, Hertel J, Turner JD. Early-Life Adversity Leaves Its Imprint on the Oral Microbiome for More Than 20 Years and Is Associated with Long-Term Immune Changes. International Journal of Molecular Sciences. 2021; 22(23):12682. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222312682

Chicago/Turabian Style

Charalambous, Eleftheria G., Sophie B. Mériaux, Pauline Guebels, Claude P. Muller, Fleur A. D. Leenen, Martha M. C. Elwenspoek, Ines Thiele, Johannes Hertel, and Jonathan D. Turner. 2021. "Early-Life Adversity Leaves Its Imprint on the Oral Microbiome for More Than 20 Years and Is Associated with Long-Term Immune Changes" International Journal of Molecular Sciences 22, no. 23: 12682. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222312682

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop