- Open Access
Drinking alcohol is associated with variation in the human oral microbiome in a large study of American adults
Microbiomevolume 6, Article number: 59 (2018)
Dysbiosis of the oral microbiome can lead to local oral disease and potentially to cancers of the head, neck, and digestive tract. However, little is known regarding exogenous factors contributing to such microbial imbalance.
We examined the impact of alcohol consumption on the oral microbiome in a cross-sectional study of 1044 US adults. Bacterial 16S rRNA genes from oral wash samples were amplified, sequenced, and assigned to bacterial taxa. We tested the association of alcohol drinking level (non-drinker, moderate drinker, or heavy drinker) and type (liquor, beer, or wine) with overall microbial composition and individual taxon abundance. The diversity of oral microbiota and overall bacterial profiles differed between heavy drinkers and non-drinkers (α-diversity richness p = 0.0059 and β-diversity unweighted UniFrac p = 0.0036), and abundance of commensal order Lactobacillales tends to be decreased with higher alcohol consumption (fold changes = 0.89 and 0.94 for heavy and moderate drinkers, p trend = 0.005 [q = 0.064]). Additionally, certain genera were enriched in subjects with higher alcohol consumption, including Actinomyces, Leptotrichia, Cardiobacterium, and Neisseria; some of these genera contain oral pathogens, while Neisseria can synthesize the human carcinogen acetaldehyde from ethanol. Wine drinkers may differ from non-drinkers in microbial diversity and profiles (α-diversity richness p = 0.048 and β-diversity unweighted UniFrac p = 0.059) after controlling for drinking amount, while liquor and beer drinkers did not. All significant differences between drinkers and non-drinkers remained after exclusion of current smokers.
Our results, from a large human study of alcohol consumption and the oral microbiome, indicate that alcohol consumption, and heavy drinking in particular, may influence the oral microbiome composition. These findings may have implications for better understanding the potential role that oral bacteria play in alcohol-related diseases.
More than 700 different bacterial species and a range of other microorganisms (archaea, fungi, and viruses) colonize the human oral cavity, known collectively as the oral microbiome [1, 2]. Oral microbiota play important roles in human health, including in immune response, carcinogen metabolism, and nutrient digestion [3,4,5]. Evidence indicates that oral microbiota dysbiosis is related to local oral diseases, such as periodontitis and dental caries  and potentially to systemic diseases, including gastrointestinal cancers [7, 8] and cardiovascular disease . However, little is known regarding exogenous exposures that cause dysbiosis of the oral microbiota.
We hypothesized that alcohol drinking habits are associated with changes in the oral microbial community. Alcohol intake may impact the oral microbiota in several ways: by direct cytotoxic effects on bacteria , by disturbing saliva-bacterium interactions [11,12,13,14], and by providing ethanol as a substrate for bacterial metabolism . Several previous studies, in both animals and humans, have also observed effects of alcohol consumption on oral bacteria. Animal studies showed that diet containing 20% ethanol increases colonization by Streptococcus mutans , a caries-related bacterium, and dramatically decreases the number of detectable bacterial species in the oral biofilms of rats . Similarly, the association of excessive drinking (> 3 drinks per day) with poor oral health was observed in a population study . In small-scale human studies, drinking at least one glass of red wine per day was associated with reduced species richness and reduction of certain anaerobic bacteria in sub- and supra-gingival plaque , while excessive co-use of tobacco and alcohol was associated with reduced species richness and decreased abundance of Neisseria, Aggregatibacter, and Fusobacteria, in oral mucosa biofilms . Aside from direct effects, alcohol may indirectly impact the oral microbiota through disturbing the host defense system [20,21,22,23], subsequently resulting in host-mediated periodontitis [24, 25]. Large population-based studies have consistently demonstrated that alcohol consumption is associated with increased risk of periodontal disease in a dose-dependent fashion [24, 26]. Evidence shows that the oral microbiome is closely tied to oral health status [27, 28]. Despite this evidence suggesting an impact of alcohol on the oral microbiome, no study has comprehensively investigated the relationship of the oral microbiome to alcohol drinking habits in terms of drinking amount and types of alcoholic beverages consumed.
We tested the relationship of level and types of alcoholic beverages, with the oral microbiome in 1044 individuals from two large US national cohorts. The oral microbiome was characterized by bacterial 16S rRNA gene sequencing. Comparisons of overall community composition and taxon abundance were conducted across groups with respect to drinking status.
Participants were drawn from the American Cancer Society (ACS) Cancer Prevention Study II (CPS II) Nutrition cohort  and the National Cancer Institute (NCI) Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (PLCO) cohort . As previously described , subjects included in the present analyses were originally selected from the CPS II and PLCO cohorts as cases or controls for collaborative nested case-control studies of the oral microbiome in relation to head and neck cancer and pancreatic cancer. Oral wash samples were collected by mail from 70,004 CPS II Nutrition cohort participants between 2000 and 2002 and in the PLCO control arm (n = 77,445) at recruitment from 1993 to 2001. After excluding participants without information on alcohol consumption status, 458 participants from CPS II (n = 169 from the head and neck study and n = 289 from the pancreas study) and 586 participants from PLCO (n = 231 from the head and neck study and n = 355 from the pancreas study) were included in the current study. All participants provided informed consent and all protocols were approved by the New York University School of Medicine Institutional Review Board.
Alcohol consumption and covariate assessment
Comprehensive demographic and lifestyle information was collected by baseline and follow-up questionnaires in the PLCO and CPS II Nutrition cohorts. Detailed information on alcohol consumption over the past year, including frequency of consumption, serving size, and type of alcoholic beverages consumed (wine, beer, and liquor), was ascertained via questionnaires that were most close in time to oral wash sample collection for both cohorts. The 137-item Food Frequency Questionnaire (FFQ) and Diet History Questionnaire (DHQ) were used in PLCO and CPS II Nutrition cohorts, respectively [32, 33]. According to the drinking level definition of the Dietary Guidelines for Americans 2010 , we defined moderate drinkers as > 0 but ≤ 1 drinks per day, on average, for women, and > 0 but ≤ 2 drinks per day, on average, for men. Women and men who had greater than one or two drinks per day, respectively, were considered heavy drinkers. This definition considers the gender differences in blood concentration of ethanol after drinking , which is equal to salivary ethanol for up to 5 h after consumption . Drinkers exclusively consuming wine, beer, or liquor were defined as wine drinkers, beer drinkers, and liquor drinkers, respectively.
Oral microbiota characterization using 16S rRNA gene amplification and sequencing
Participants in both cohorts were asked to swish vigorously with 10 mL Scope mouthwash (P&G) for 30 s and then to expectorate into a specimen tube. Samples were shipped to each cohort’s biorepository, pelleted, and stored at − 80 °C until use [30, 37]. Bacterial genomic DNA was extracted from the samples using the MoBio PowerSoil DNA Isolation Kit (Carlsbad, CA), with the bead-beating method in the MoBio Powerlyzer instrument. As reported previously , 16S rRNA gene sequencing on the extracted DNA was performed. 16S rRNA amplicon libraries were generated using primers incorporating FLX Titanium adapters and a sample barcode sequence, allowing unidirectional sequencing covering variable regions V3 to V4 (Primers: 347F- 5′GGAGGCAGCAGTAAGGAAT-3′ and 803R- 5′CTACCGGGGTATCTAATCC-3′). Five nanograms of genomic DNA was used as the template in 25 uL PCR reaction buffer for 16S rRNA amplicon preparation. Cycling conditions were one cycle of 94 °C for 3 min, followed by 25 cycles of 94 °C for 15 s, 52 °C for 45 s, and 72 °C for 1 min followed by a final extension of 72 °C for 8 min. The generated amplicons were then purified using Agencourt AMPure XP kit (Beckman Coulter, CA). Purified amplicons were quantified by fluorometry using the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, CA). Equimolar amounts (107 molecules/uL) of purified amplicons were pooled for sequencing. Pyrosequencing (Roche 454 GS FLX Titanium) was carried out according to the manufacturer’s instructions .
Upstream sequence analysis of microbiome data
16S rRNA gene amplicon sequences were processed and analyzed using the QIIME pipeline . Multiplexed libraries were deconvoluted based on the barcodes assigned to each sample. Poor-quality sequences were excluded using the default parameters of the QIIME script split_libraries.py (minimum average quality score = 25, minimum/maximum sequence length = 200/1000 base pairs, no ambiguous base calls, and no mismatches allowed in the primer sequence). From 1044 pre-diagnostic oral wash samples, we obtained 11,395,395 high-quality 16S rRNA gene sequence reads (mean 10,915 [SD = 3049] sequences per sample), with similar number of reads in all cohorts . Filtered sequence reads were clustered into operational taxonomic units (OTUs) and subsequently assigned to taxa by using the Human Oral Microbiome Database (HOMD) pre-defined taxonomy map of reference sequences with ≥ 97% identity . Summary of sequence reads per sample which were assigned to HOMD reference is shown in Additional file 1: Table S1.
Blinded positive quality control (QC) specimens were used across all sequencing batches. We previously reported good agreement of microbiome parameters in replicates from these QC subjects . Negative control samples (with Scope mouthwash only) were used to detect possible bacteria in Scope and reagent, as well as environment contamination in all sequencing batches. No DNA was detected from the negative control samples. We further conducted an independent experiment which compared the oral microbiota in paired mouthwash samples and immediately frozen whole saliva samples from 10 healthy subjects. The results showed that oral microbial profiles in mouthwash samples were similar to the profiles in saliva samples (Additional file 1: Figure S1). Thus, the mouthwash samples are suitable to test the hypothesis in this study.
The effect of drinking level was assessed by comparing subjects in each drinking category (moderate, heavy) to non-drinkers. For trend tests, drinking level was treated as a continuous variable by assigning the numbers 0, 1, and 2 to non-, moderate, and heavy drinkers, respectively. The effect of drinking type (wine, beer, or liquor) was examined by modeling the three dichotomous variables of any consumption of wine, beer or liquor, while additionally adjusting for multiple drinking types using cross-product terms (e.g., formula to test effect of wine drinking: ~any_wine + any_beer + any_liquor + any_wine*any_beer + any_wine*any_liquor). To control for the effects of potential confounders, all models were adjusted for age, sex, race, BMI category (normal weight, overweight, obese), smoking status (never, former, current), education level (no college, some college, college graduate), and study (CPSII-a, CPSII-b, PLCO-a, PLCO-b); drinking type models were additionally adjusted for drinking amount (grams of ethanol per day).
We assessed α-diversity (within-subject diversity) using numbers of observed species (richness) and the inverse Simpson’s Index (evenness). These α-diversity indices were calculated in 500 iterations of rarefied OTU tables with the minimum sequencing depth of 1325 among all study subjects. The average over the iterations was taken for each participant. Linear regression with covariate adjustment was used to examine the difference of α-diversity indices among drinking groups. β-diversity (between-subject diversity) was assessed using unweighted and weighted UniFrac distance matrices accounting for both presence or absence of observed organisms and their abundance, respectively . We performed Permutational Multivariate Analysis of Variance (PERMANOVA; adonis function, vegan package, R)  and partial constrained analysis of principal coordinates (partial CAP; “capscale” function, vegan package, R)  to examine statistically and visually whether bacterial community profiles differed by drinking level or type. Unlike the commonly used unconstrained principal coordinate analysis (PCoA) representation of UniFrac distances, partial CAP displays the variation attributable to alcohol consumption with covariate adjustment. We also conducted pairwise comparisons among drinking groups for each of the first three coordinates in unconstrained PCoA, using the Kruskal-Wallis post-hoc test (Dunn’s test). The community type analysis of the oral wash samples was performed with Dirichlet multinomial mixture (DMM) model using counts of sequencing reads at the genus level .
Operational taxonomic units (OTUs) were classified into 13 phyla, 23 classes, 37 orders, 76 families, 227 genera, and 341 species, according to their alignment with the HOMD reference database. We used the “DESeq” function within the DESeq2 package  in R to test for differentially abundant taxa by drinking level and types of alcoholic beverages, at each taxonomic level. This function models raw counts using a negative binomial distribution and adjusts internally for “size factors” which normalize for differences in sequencing depth between samples. We analyzed taxa within the major oral phyla (carried by more than 90% of study subjects). Additionally, we filtered to include only taxa with greater than 2 sequences in at least 100 participants using R command in DESeq “rowSums(counts(OTU.deseq) > 2) ≥ 100” to remove low-count taxa for the phylum through species level analysis. This resulted in inclusion of 5 phyla, 9 classes, 15 orders, 30 families, 96 genera, and 95 species in the analysis. DESeq2 default outlier replacement and filtering of count outliers were turned off. Taxa models with a maximum Cook’s distance > 10 were removed prior to p value adjustment for the false discovery rate (FDR). To account for the multiple comparisons at each taxonomic level, we considered an FDR-adjusted p value (q value) less than 0.10 as significant. Pearson’s linear correlation was used to explore the correlation among the log-transformed DESeq2-normalized taxa abundance. All statistical tests were two-sided, and all statistical analyses were carried out using R version 3.2.1.
To examine if smoking, future cancer case/control status, or oral health status (periodontal disease and caries) has confounding effects on the observed alcohol-microbiome associations, we conducted several sensitivity analyses: (a) excluding current smokers, (b) adjusting for future case/control status in models, (c) excluding subjects who had detectable periodontal pathogens Porphyromonas gingivalis and/or Aggregatibacter actinomycetemcomitans in their oral samples (surrogate markers for periodontal disease ), and (d) stratified analysis by the median abundance of Streptococcus mutans (surrogate marker for dental caries ). The latter two analyses used surrogate markers of oral disease as we lacked information on oral disease in our study.
Of the 1044 participants in this study, 25.9% (n = 270) were non-drinkers, 58.8% (n = 614) were moderate drinkers, and 15.3% (n = 160) were heavy drinkers. The study participants were predominantly above middle age (mean 67.7, range 55–87) and White (95%), with age and race distributions not differing significantly by drinking level. The alcohol drinking groups had higher percentages of men and smokers and tended to be leaner and more educated than the non-drinker group (Table 1). Among alcohol drinkers, 13.0% (n = 101) were wine-only drinkers, 5.0% (n = 39) were beer-only drinkers, and 3.4% (n = 26) were liquor-only drinkers (Table 2).
Drinking level analysis
We first examined the overall microbial composition according to drinking level. In the α-diversity analysis, richness (Fig. 1a, c) increased in heavy drinkers (HD) and moderate drinkers (MD) with statistical significance (meanHD = 100.3, p from linear regression = 0.0059 and meanMD = 100.1, p = 0.0073) relative to non-drinkers (meannon-drinker = 94.6). These differences in richness remained significant after further adjustment for case/control status (p = 0.0051 for HD and 0.0072 for MD). Evenness tended to be higher in drinkers compared to non-drinkers (Fig. 1b, d); however, the differences were not statistically significant (p = 0.092 and 0.062 for HD and MD, respectively). When assessing β-diversity according to the unweighted UniFrac distance, partial CAP revealed that HD and MD separated from non-drinkers on the first CAP axis, based on position of the group centroids (Fig. 2a). Similarly, HD and MD differed from non-drinkers in the first principal coordinate in PCoA (p = 0.029 and 0.038 respectively) (Fig. 2c). In PERMANOVA analysis of the unweighted UniFrac distance, HD significantly differed from non-drinkers after controlling for covariates (p = 0.0036), while MD did not (p = 0.22). The difference between HD and non-drinkers in PERMANOVA analysis persisted after further adjustment for case/control status (p = 0.0050). When assessing β-diversity according to the weighted UniFrac distance, HD, MD, and non-drinkers clustered together in partial CAP plot (Fig. 2b) and showed no differences in the first three principal coordinates in PCoA, with the exception of the third coordinate comparing MD and non-drinkers (p = 0.018) (Fig. 2d). In PERMANOVA analysis of the weighted UniFrac distance, neither HD nor MD were found to be different from non-drinkers (p = 0.22 and 0.24 for HD and MD, respectively). The differences in the microbial community composition according to alcohol drinking were confirmed by the results of community type classification (Additional file 1: Figure S2).
Associations between drinking level and the abundance of specific oral taxa were examined using negative binomial generalized linear models  (Table 3 and Fig. 3a). Abundance of class Bacilli in phylum Firmicutes was decreased with higher drinking levels (fold changes [FCs] = 0.92 and 0.94 for HD and MD, FDR-adjusted q trend = 0.050), as was its major order Lactobacillales (FCs = 0.89 and 0.94 for HD and MD, q = 0.064). Low-abundance genera Streptococcus and Lachnoanaerobaculum in Firmicutes were increased with higher drinking levels. Taxa in other phyla were also enriched with higher drinking level, including the genus Actinomyces (FCs = 1.41 and 1.15 for HD and MD, q = 0.041) and its constituent species A. graevenitzzi (FCs = 1.53 and 1.17 for HD and MD, q = 0.088). Leptotrichia (FCs = 1.61 and 1.24 for HD and MD, q = 0.0027) and low-abundance genera Cardiobacterium (FCs = 1.41 and 1.61 for HD and MD, q = 0.041) and Neisseria (FCs = 2.12 and 2.42 for HD and MD, q = 0.027) were increased in both heavy and moderate drinkers, although read counts tended to be low for these bacteria (Additional file 1: Figure S3). Additional taxa that were differentially abundant by drinking level are shown in Table 3 and Fig. 3a. The correlation matrix of the taxa associated with drinking level indicated that class Bacilli and order Lactobacillales, which were decreased with higher drinking level, generally had negative correlations with other taxa that were increased with higher drinking level (Fig. 3b). Results were not altered by further adjustment for case/control status.
We examined the homogeneity of our results across two independent cohorts. The partial CAP plots indicated a similar trend of HD separated from non-drinkers on the first CAP axis in both cohorts (Additional file 1: Figure S4a–d), which was further supported by the PERMANOVA analysis showing that HD differed from non-drinkers according to the unweighted UniFrac distance in both CPS II and PLCO cohorts (p = 0.097 and p = 0.064 respectively). Most of the taxa associated with drinking level showed similar trends in both cohorts (Additional file 1: Figure S3 and Table S2), with the exception of family Bacteroidales[F-2], its constituent genus Bacteroidales[G-2], and genus Corynebacterium, which all were shown to increase significantly with drinking level in the PLCO, but not CPS II, cohorts.
In a sensitivity analysis to determine potential influences of smoking and oral disease on study results, the observed alcohol-microbiota associations (i.e., in α-diversity, β-diversity, and taxon abundance) were unchanged after excluding current smokers or subjects carrying periodontal pathogens (Additional file 1: Table S3). Additionally, similar trends for alcohol-microbiota associations were observed in those with low or high abundance of S. mutans, a marker of dental caries (Additional file 1: Table S3).
Drinking type analysis
We next investigated the overall microbial composition and taxon abundances for exclusive wine (n = 101), beer (n = 39), and liquor (n = 26) drinkers, relative to non-drinkers, after controlling for drinking amount. In the α-diversity analysis (Additional file 1: Figure S5), beer drinkers and liquor drinkers did not differ from non-drinkers in richness (p from linear regression = 0.87 and 0.91 for beer and liquor drinkers, respectively), while wine drinkers (mean = 100.9, p = 0.048) had increased richness compared to non-drinkers (mean = 94.6). Significance remained after further adjustment for case/control status (p = 0.048) and exclusion of current smokers (p = 0.053). Similarly, in covariate-adjusted PERMANOVA analysis of the unweighted UniFrac distance, wine drinkers differed from non-drinkers with marginal significance (p = 0.059), but liquor drinkers or beer drinkers did not (p = 0.73 and p = 0.59). Evenness (Additional file 1: Figure S5) and β-diversity estimated by weighted UniFrac distance did not reveal any differences by drinking type. Compared to non-drinkers, wine drinking was associated with decreased abundance of genus Peptococcus (FC = 0.45, q = 0.073); beer drinking was associated with decreased abundance of Porphyromonas (FC = 0.32, q = 0.071) and increased abundance of genus Parascardovia (FC = 5.87, q = 0.016); and liquor drinking was associated with decreased abundance of Lachnospiraceae[G-2] (FC = 0.39, q = 0.090) (Fig. 4 & Additional file 1: Table S4). Other taxa associated with wine, beer, or liquor drinking were in genera Corynebacterium, Prevotella, and Aggregatibacter and Eikenella, which had been identified as also associated with drinking level (Fig. 4 and Additional file 1: Table S4).
In this first large comprehensive human study, we observed that overall microbiome community composition in the oral cavity differed by drinking level, with consistent results across two independent cohorts and after exclusion of current smokers. Drinkers had decreased abundance of order Lactobacillales, the major order in the Firmicutes phylum. Other taxa, some of which are potentially pathogenic, were enriched with higher alcohol consumption. Additionally, wine drinkers may have a shifted microbiome community composition relative to non-drinkers. These results suggest that alcohol drinking and beverage type may influence the oral microbiota.
We studied the oral microbiome in mouthwash samples, which contain bacteria shed from adhering microbial communities on various oral sites, including tooth surfaces, gingival crevices, tongue dorsum, and buccal mucosa [49,50,51]. We found an increase in within-subject diversity (richness) and qualitative variation  in microbial community profiles (unweighted UniFrac) by drinking level. The observed increased diversity and altered profiles in drinkers may be due to direct effects of alcohol or may reflect poor oral health conditions in drinkers. Alcohol intake has been associated with risk of periodontal diseases and caries [17, 24, 26]. Additionally, poor oral health conditions, including higher plaque index, presence of decayed teeth, gingival bleeding, and deeper periodontal pockets, have been associated with higher phylogenetic diversity of salivary microbiota . In a sensitivity analysis, we used periodontal pathogens (P. gingivalis and A. actinomycetemcomitans) and the abundance of S. mutans in saliva as surrogate markers for periodontitis and dental caries, respectively, and major findings remained unchanged. Thus, it is possible that alcohol consumption has a direct effect on oral bacteria composition independently of oral disease.
Our results indicate that alcohol consumption is associated with decreased abundance of Lactobacillales, a dominant order in class Bacilli and phylum Firmicutes. Bacteria in the order Lactobacillales (also named lactic acid bacteria or LAB), produce lactic acid as an end-product of carbohydrate fermentation  and are among the most common microbes employed as probiotics . Firmicutes and Lactobacillales abundance have also been shown to decrease in the intestines of mice fed with ethanol chronically [54, 55], while increases were observed in potentially pathogenic bacteria, Proteobacteria and Actinobacteria, in response to ethanol feeding . While the mechanism for ethanol-induced decreases in Lactobacillales are unclear, depletion of Lactobacillales may promote growth of alkaline-tolerant bacteria , which may explain the inverse correlations of Lactobacillales with the other taxa associated with alcohol drinking. In support of this, probiotic treatment of ethanol-fed mice with Lactobacillus rhamnosus GG decreased luminal pH and prevented ethanol-related pathogenic increases in the gut microbiome . Evidence shows that Lactobacillales have a beneficial effect on oral health . Some Lactobacillales can reduce the risk of caries development , possibly through competitive exclusion and displacement of pathogens, or by production of antibacterial compounds [58,59,60]. Additionally, oral Lactobacilli species can suppress the growth of periodontal pathogens , and oral intake of some Lactobacillales can promote reduction of gingival inflammation . Thus, the alcohol-related depletion of beneficial commensal Lactobacillales may lead to further oral health inflammation-related disturbances.
We also identified some taxa enriched in drinkers. Ethanol may indirectly increase certain bacterial taxa by decreasing Lactobacillales thereby increasing pH, as mentioned above, or by inhibiting the antimicrobial properties of saliva and by disturbing the host-microbial balance. Both acute and chronic ethanol exposure can lead to functional changes in saliva, including decreased flow rate and impaired output of total protein, amylase, and electrolytes [12, 14]. Additionally, alcohol could impair neutrophil function (contributing to bacterial overgrowth and increased bacterial penetration) , reduce monocyte production of inflammatory cytokines (allowing for microbial proliferation) , and have adverse effects on teeth (stimulate bone resorption and suppress bone formation)  and the periodontium . These detrimental effects of alcohol on host defense potentially lead to periodontitis . Among the taxa enriched in drinkers, the genera Aggregatibacter, Actinomyces, Kingella, Leptotrichia, Cardiobacterium, Bacteroidales[G-2], and Prevotella may contain human pathogenic and cardiogenic pathogens. Aggregatibacter, a genus in the Pasteurellaceae family of the Proteobacteria phylum, is related to periodontal disease [63,64,65,66] and infective endocarditis [67, 68], as are species in Actinomyces, a genus from the Actinobacteria phylum [69, 70], and Kingella, a genus from the Proteobacteria phylum [68, 71, 72]. Leptotrichia is one of the two major genera in phylum Fusobacteria, which tends to cause disease in the presence of local or general pre-disposing factors and has been reported to be involved in various human infections . The only two species in genus Cardiobacterium, C. hominis and C. valvarum, both can cause endocarditis [74,75,76]. The constituent species oral taxon 274 in Bacteroidales[G-2] was shown to be more prevalent in subgingival sites in periodontitis subjects than in healthy subjects . Species of genus Prevotella are part of the indigenous microbiota of mucous membranes from oral cavity, while some black-pigmented Prevotella, for instance P. intermedia and P. melaninogenica, associate with periodontal disease . Thus, the changes in the composition of oral microbiome related to alcohol consumption may potentially contribute to dental caries, periodontal diseases, and other health consequences.
In addition to these compositional changes related to alcohol use, alcohol consumption may result in the bacterial production from ethanol of acetaldehyde, a World Health Organization (WHO) group 1 human carcinogen . Acetaldehyde is toxic, mutagenic, and carcinogenic in human cells and animal models [80,81,82,83]. Certain bacteria are involved in producing acetaldehyde from ethanol, while others can beneficially metabolize acetaldehyde to less toxic compounds. Neisseria is mainly a non-pathogenic genus in the oral cavity, but it has extremely high alcohol dehydrogenase (ADH) activity and produces significant amounts of acetaldehyde from ethanol . The increased abundance of Neisseria with alcohol drinking that we have observed is consistent with the ADH activity of Neisseria and is supported by another human study . Conversely, strains of gastrointestinal Lactobacillus, decreased in alcohol consumers in our study, show a high capacity to metabolize acetaldehyde to less toxic forms . Thus, ethanol-related increase in abundance of Neisseria, which produces acetaldehyde, and decrease in abundance of LAB, involved in acetaldehyde degradation, suggest that oral bacteria may play a role in alcohol-related carcinogenesis through acetaldehyde.
Analyses by beverage type revealed that wine drinkers may have a higher richness and different microbial profiles from non-drinkers, with decreased abundance of phyla Bacteroidetes and Firmicutes, and family Peptostreptococcaceae. These findings are consistent with previous in vivo and in vitro studies investigating the antimicrobial properties of wine, showing that wine drinkers had decreased abundance of certain species in sub- and supra-gingival plaques, including species in Peptostreptococcaceae [18, 86]. The results on microbial richness in these studies are inconsistent with our findings, which may relate to different sample types used (subgingival plaque vs. mouthwash samples). Beer has a relatively low alcohol concentration but a complex chemical composition , containing a greater abundance of proteins, B type vitamins, fiber, minerals, antioxidants, and varied flavors than other alcoholic beverage types , which potentially may have differential effects on the growth of oral bacteria. Due to the high concentration of ethanol in liquor, the bacterial variations associated with liquor drinking may only more closely represent the effects of pure ethanol. Bacterial taxa similarly altered relative to non-drinkers in wine, beer, and liquor drinkers (e.g., Prevotella_gid_170) may simply reflect an effect of drinking alcohol vs. not drinking. Because our study had limited numbers of subjects who exclusively consumed beer, wine, or liquor, further study is required to disentangle the differential effect of each type of alcoholic beverage on oral microbial composition.
Our study has several strengths. First, the use of 16S rRNA gene sequencing for microbiome analysis allowed us to comprehensively study overall bacterial community composition and specific oral taxon abundances. Second, our very large sample size provided excellent statistical power to detect variation in the oral microbiome with alcohol drinking and allowed us to confirm our findings in multiple independent study groups. Finally, the detailed information on alcohol drinking amount and type allowed us to classify alcohol intake in detail, and the detailed demographic and lifestyle information allowed us to adjust for potential confounding factors. A limitation of our study is that it is observational, limiting our ability to establish a causal relationship. A trial where moderate or heavy drinkers were randomized to continue or stop drinking for a long enough time to influence the oral microbiome could provide more definitive information. Additionally, a few enriched bacteria in drinkers has overall low abundance. Further, the majority of this study population was White and above middle age, limiting the generalizability of our findings to other races and age groups. We lacked information on the oral health status of the study participants, though sensitivity analysis using bacterial markers as proxies of dental caries and periodontitis suggests that the observed alcohol-microbiota associations are independent of oral diseases. Lastly, the oral wash sample might not be representative of oral bacterial profiles in specific oral niches, such as the subgingival plaque, which is more directly involved in oral disease.
We found that alcohol consumption is related to overall oral microbiome community composition and to the abundance of specific oral taxa. Heavy drinking may influence bacterial composition, including potential depletion of beneficial commensal bacteria and increased colonization of potentially pathogenic bacteria. Such changes potentially contribute to alcohol-related diseases, including periodontal disease, head and neck cancer, and digestive tract cancers, but further research is needed to relate alcohol-related composition changes to disease phenotypes. The taxa we have identified can be further investigated to tease out their potential relationship with underlying oral health status and to elucidate their potential role in alcohol-related health consequences. Future studies should also investigate the impact of alcohol drinking on the metagenomic (functional) content of the oral microbiome. Improved understanding of the causes and health impacts of oral dysbiosis can lead to microbiome-targeted approaches for disease prevention.
American Cancer Society
Constrained analysis of principal coordinates
- CPS II:
Cancer Prevention Study II
False discovery rate
Human Oral Microbiome Database
Lactic acid bacteria
National Cancer Institute
Operational taxonomic units
Principal coordinate analysis
Permutational Multivariate Analysis of Variance
Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial
World Health Organization
Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, Lakshmanan A, Wade WG. The human oral microbiome. J Bacteriol. 2010;192(19):5002–17.
Baker JL, Bor B, Agnello M, Shi W, He X. Ecology of the oral microbiome: beyond bacteria. Trends Microbiol. 2017;25(5):362–74.
Slocum C, Kramer C, Genco CA. Immune dysregulation mediated by the oral microbiome: potential link to chronic inflammation and atherosclerosis. J Intern Med. 2016;280(1):114–28.
Homann N, Tillonen J, Meurman JH, Rintamaki H, Lindqvist C, Rautio M, Jousimies-Somer H, Salaspuro M. Increased salivary acetaldehyde levels in heavy drinkers and smokers: a microbiological approach to oral cavity cancer. Carcinogenesis. 2000;21(4):663–8.
Moye ZD, Zeng L, Burne RA. Fueling the caries process: carbohydrate metabolism and gene regulation by Streptococcus mutans. J Oral Microbiol. 2014;6(1):24878.
Wade WG. The oral microbiome in health and disease. Pharmacol Res. 2013;69(1):137–43.
Ahn J, Chen CY, Hayes RB. Oral microbiome and oral and gastrointestinal cancer risk. Cancer Causes Control. 2012;23(3):399–404.
Fan X, Alekseyenko AV, Wu J, Peters BA, Jacobs EJ, Gapstur SM, Purdue MP, Abnet CC, Stolzenberg-Solomon R, Miller G, et al. Human oral microbiome and prospective risk for pancreatic cancer: a population-based nested case-control study. Gut. 2018;67(1):120–7.
Koren O, Spor A, Felin J, Fak F, Stombaugh J, Tremaroli V, Behre CJ, Knight R, Fagerberg B, Ley RE, et al. Human oral, gut, and plaque microbiota in patients with atherosclerosis. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4592–8.
Ingram LO. Ethanol tolerance in bacteria. Crit Rev Biotechnol. 1990;9(4):305–19.
Scannapieco FA. Saliva-bacterium interactions in oral microbial ecology. Crit Rev Oral Biol Med. 1994;5(3–4):203–48.
Enberg N, Alho H, Loimaranta V, Lenander-Lumikari M. Saliva flow rate, amylase activity, and protein and electrolyte concentrations in saliva after acute alcohol consumption. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2001;92(3):292–8.
Tenovuo J. Antimicrobial function of human saliva—how important is it for oral health? Acta Odontol Scand. 1998;56(5):250–6.
Maier H, Born IA, Mall G. Effect of chronic ethanol and nicotine consumption on the function and morphology of the salivary glands. Klin Wochenschr. 1988;66(Suppl 11):140–50.
Kantorski KZ, de Souza DM, Yujra VQ, Junqueira JC, Jorge AO, da Rocha RF. Effect of an alcoholic diet on dental caries and on Streptococcus of the mutans group. Study in rats. Braz Oral Res. 2007;21(2):101–5.
Jabbour Z, Nascimento C, Kotake BGS, El-Hakim M, Henderson JE, de Albuquerque Junior RF. Assessing the oral microbiota of healthy and alcohol-treated rats using whole-genome DNA probes from human bacteria. Arch Oral Biol. 2013;58(3):317–23.
Jansson L. Association between alcohol consumption and dental health. J Clin Periodontol. 2008;35(5):379–84.
Signoretto C, Bianchi F, Burlacchini G, Sivieri F, Spratt D, Canepari P. Drinking habits are associated with changes in the dental plaque microbial community. J Clin Microbiol. 2010;48(2):347–56.
Thomas AM, Gleber-Netto FO, Fernandes GR, Amorim M, Barbosa LF, Francisco AL, de Andrade AG, Setubal JC, Kowalski LP, Nunes DN, et al. Alcohol and tobacco consumption affects bacterial richness in oral cavity mucosa biofilms. BMC Microbiol. 2014;14:250.
Szabo G. Consequences of alcohol consumption on host defence. Alcohol Alcohol. 1999;34(6):830–41.
Szabo G, Mandrekar P, Girouard L, Catalano D. Regulation of human monocyte functions by acute ethanol treatment: decreased tumor necrosis factor-alpha, interleukin-1 beta and elevated interleukin-10, and transforming growth factor-beta production. Alcohol Clin Exp Res. 1996;20(5):900–7.
Pepersack T, Fuss M, Otero J, Bergmann P, Valsamis J, Corvilain J. Longitudinal study of bone metabolism after ethanol withdrawal in alcoholic patients. J Bone Miner Res Off J Am Soc Bone Miner Res. 1992;7(4):383–7.
Ogden GR, Wight AJ, Rice P. Effect of alcohol on the oral mucosa assessed by quantitative cytomorphometry. J Oral Pathol Med. 1999;28(5):216–20.
Pitiphat W, Merchant AT, Rimm EB, Joshipura KJ. Alcohol consumption increases periodontitis risk. J Dent Res. 2003;82(7):509–13.
Bartold PM, Van Dyke TE. Periodontitis: a host-mediated disruption of microbial homeostasis. Unlearning learned concepts. Periodontol. 2013;62(1):203–17.
Tezal M, Grossi SG, Ho AW, Genco RJ. Alcohol consumption and periodontal disease. The third national health and nutrition examination survey. J Clin Periodontol. 2004;31(7):484–8.
Takeshita T, Kageyama S, Furuta M, Tsuboi H, Takeuchi K, Shibata Y, Shimazaki Y, Akifusa S, Ninomiya T, Kiyohara Y, et al. Bacterial diversity in saliva and oral health-related conditions: the Hisayama study. Sci Rep. 2016;6:22164.
Darveau RP. Periodontitis: a polymicrobial disruption of host homeostasis. Nat Rev Microbiol. 2010;8(7):481–90.
Calle EE, Rodriguez C, Jacobs EJ, Almon ML, Chao A, McCullough ML, Feigelson HS, Thun MJ. The American cancer society cancer prevention study II nutrition cohort: rationale, study design, and baseline characteristics. Cancer. 2002;94(9):2490–501.
Hayes RB, Reding D, Kopp W, Subar AF, Bhat N, Rothman N, Caporaso N, Ziegler RG, Johnson CC, Weissfeld JL, et al. Etiologic and early marker studies in the prostate, lung, colorectal and ovarian (PLCO) cancer screening trial. Control Clin Trials. 2000;21(6 Suppl):349S–55S.
Wu J, Peters BA, Dominianni C, Zhang Y, Pei Z, Yang L, Ma Y, Purdue MP, Jacobs EJ, Gapstur SM, et al. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016;10(10):2435–46.
Haftenberger M, Heuer T, Heidemann C, Kube F, Krems C, Mensink GB. Relative validation of a food frequency questionnaire for national health and nutrition monitoring. Nutr J. 2010;9:36.
Subar AF, Thompson FE, Kipnis V, Midthune D, Hurwitz P, McNutt S, McIntosh A, Rosenfeld S. Comparative validation of the Block, Willett, and National Cancer Institute food frequency questionnaires: the eating at America’s table study. Am J Epidemiol. 2001;154(12):1089–99.
United States. Department of Health and Human Services., United States. Department of Agriculture., United States. Dietary Guidelines Advisory Committee. Dietary guidelines for Americans, 2010. 7th ed. Washington, D.C: Goverment Printing Office; 2010.
Mumenthaler MS, Taylor JL, O'Hara R, Yesavage JA. Gender differences in moderate drinking effects. Alcohol Res Health. 1999;23(1):55–64.
Jones AW. Distribution of ethanol between saliva and blood in man. Clin Exp Pharmacol Physiol. 1979;6(1):53–9.
Feigelson HS, Rodriguez C, Welch R, Hutchinson A, Shao W, Jacobs K, Diver WR, Calle EE, Thun MJ, Hunter DJ, et al. Successful genome-wide scan in paired blood and buccal samples. Cancer Epidemiol Biomark Prev. 2007;16(5):1023–5.
Wu J, Lin I, Hayes RB, Ahn J. Comparison of DNA extraction methods for human oral microbiome research. Br J Med Med Res. 2014;4(10):1980–91.
Nossa CW, Oberdorf WE, Yang LY, Aas JA, Paster BJ, DeSantis TZ, Brodie EL, Malamud D, Poles MA, Pei ZH. Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J Gastroenterol. 2010;16(33):4135–44.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Chen T, Yu WH, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE. The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database. 2010;2010:baq013.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.
Anderson MJ, Willis TJ. Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology. 2003;84(2):511–25.
Holmes I, Harris K, Quince C. Dirichlet multinomial mixtures: generative models for microbial metagenomics. PLoS One. 2012;7(2):e30126.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Newman MG. Carranza’s clinical periodontology. 11th ed. Philadelphia: London: Saunders; 2012.
Loesche WJ. Role of Streptococcus mutans in human dental decay. Microbiol Rev. 1986;50(4):353–80.
Zhou Y, Gao H, Mihindukulasuriya KA, La Rosa PS, Wylie KM, Vishnivetskaya T, Podar M, Warner B, Tarr PI, Nelson DE, et al. Biogeography of the ecosystems of the healthy human body. Genome Biol. 2013;14(1):R1.
Segata N, Haake SK, Mannon P, Lemon KP, Waldron L, Gevers D, Huttenhower C, Izard J. Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples. Genome Biol. 2012;13(6):R42.
Mager DL, Ximenez-Fyvie LA, Haffajee AD, Socransky SS. Distribution of selected bacterial species on intraoral surfaces. J Clin Periodontol. 2003;30(7):644–54.
Wright AV, Lars A. Lactic acid bacteria. In: Lactic acid bacteria, microbiological and functional aspects. Boca Raton: CRC Press, Taylor & Francis Group; 2011. p. 1–16.
Ljungh A, Wadstrom T. Lactic acid bacteria as probiotics. Curr Issues Intest Microbiol. 2006;7(2):73–89.
Bull-Otterson L, Feng W, Kirpich I, Wang Y, Qin X, Liu Y, Gobejishvili L, Joshi-Barve S, Ayvaz T, Petrosino J, et al. Metagenomic analyses of alcohol induced pathogenic alterations in the intestinal microbiome and the effect of lactobacillus rhamnosus GG treatment. PLoS One. 2013;8(1):e53028.
Yan AW, Fouts DE, Brandl J, Starkel P, Torralba M, Schott E, Tsukamoto H, Nelson KE, Brenner DA, Schnabl B. Enteric dysbiosis associated with a mouse model of alcoholic liver disease. Hepatology. 2011;53(1):96–105.
Jukka HM, Iva S. Lactic acid bacteria in oral health. In: Lactic acid bacteria. Boca Raton: CRC Press; 2011. p. 403–21.
Nase L, Hatakka K, Savilahti E, Saxelin M, Ponka A, Poussa T, Korpela R, Meurman JH. Effect of long-term consumption of a probiotic bacterium, Lactobacillus rhamnosus GG, in milk on dental caries and caries risk in children. Caries Res. 2001;35(6):412–20.
Simark-Mattsson C, Emilson CG, Hakansson EG, Jacobsson C, Roos K, Holm S. Lactobacillus-mediated interference of mutans streptococci in caries-free vs. caries-active subjects. Eur J Oral Sci. 2007;115(4):308–14.
Strahinic I, Busarcevic M, Pavlica D, Milasin J, Golic N, Topisirovic L. Molecular and biochemical characterizations of human oral lactobacilli as putative probiotic candidates. Oral Microbiol Immunol. 2007;22(2):111–7.
Lang C, Bottner M, Holz C, Veen M, Ryser M, Reindl A, Pompejus M, Tanzer JM. Specific Lactobacillus/Mutans Streptococcus co-aggregation. J Dent Res. 2010;89(2):175–9.
Koll-Klais P, Mandar R, Leibur E, Marcotte H, Hammarstrom L, Mikelsaar M. Oral lactobacilli in chronic periodontitis and periodontal health: species composition and antimicrobial activity. Oral Microbiol Immunol. 2005;20(6):354–61.
Twetman S, Derawi B, Keller M, Ekstrand K, Yucel-Lindberg T, Stecksen-Blicks C. Short-term effect of chewing gums containing probiotic Lactobacillus reuteri on the levels of inflammatory mediators in gingival crevicular fluid. Acta Odontol Scand. 2009;67(1):19–24.
Zambon JJ. Actinobacillus-actinomycetemcomitans in human periodontal-disease. J Clin Periodontol. 1985;12(1):1–20.
Slots J, Reynolds HS, Genco RJ. Actinobacillus-actinomycetemcomitans in human periodontal-disease—a cross-sectional microbiological investigation. Infect Immun. 1980;29(3):1013–20.
Slots J, Zambon JJ, Rosling BG, Reynolds HS, Christersson LA, Genco RJ. Actinobacillus-actinomycetemcomitans in human periodontal-disease—association, serology, leukotoxicity, and treatment. J Periodontal Res. 1982;17(5):447–8.
Zambon JJ, Christersson LA, Slots J. Actinobacillus-actinomycetemcomitans in human periodontal-disease-prevalence in patient groups and distribution of biotypes and serotypes within families. J Periodontol. 1983;54(12):707–11.
Norskov-Lauritsen N. Classification, identification, and clinical significance of haemophilus and aggregatibacter species with host specificity for humans. Clin Microbiol Rev. 2014;27(2):214–40.
Revest M, Egmann G, Cattoir V, Tattevin P. HACEK endocarditis: state-of-the-art. Expert Rev Anti-Infect Ther. 2016;14(5):523–30.
Vielkind P, Jentsch H, Eschrich K, Rodloff AC, Stingu CS. Prevalence of Actinomyces spp. in patients with chronic periodontitis. Int J Med Microbiol. 2015;305(7):682–8.
Kottam A, Kaur R, Bhandare D, Zmily H, Bheemreddy S, Brar H, Herawi M, Afonso L. Actinomycotic endocarditis of the eustachian valve: a rare case and a review of the literature. Tex Heart Inst J. 2015;42(1):44–9.
Yagupsky P, Porsch E, St Geme JW 3rd. Kingella kingae: an emerging pathogen in young children. Pediatrics. 2011;127(3):557–65.
Chen C. Distribution of a newly described species, Kingella oralis, in the human oral cavity. Oral Microbiol Immunol. 1996;11(6):425–7.
Eribe ER, Olsen I. Leptotrichia species in human infections. Anaerobe. 2008;14(3):131–7.
Brouqui P, Raoult D. Endocarditis due to rare and fastidious bacteria. Clin Microbiol Rev. 2001;14(1):177–207.
Berbari EF, Cockerill FR, Steckelberg JM. Infective endocarditis due to unusual or fastidious microorganisms. Mayo Clin Proc. 1997;72(6):532–42.
Bothelo E, Gouriet F, Fournier PE, Roux V, Habib G, Thuny F, Metras D, Raoult D, Casalta JP. Endocarditis caused by Cardiobacterium valvarum. J Clin Microbiol. 2006;44(2):657–8.
Li CL, Liang JP, Jiang YT. Association of uncultivated oral phylotypes AU126 and X112 with periodontitis. Oral Dis. 2006;12(4):371–4.
Mario JAC, Maria RLS, Elerson G-J. Prevotella. In: Molecular detection of human bacterial pathogens. Boca Raton: CRC Press, Taylor & Francis Group; 2011. p. 585–97.
World Health Organization. Iarc monographs on the evaluation of carcinogenic risks to humans. Iarc Monog Eval Carc. 2010;95:9–38.
Higginson J, DeVita VT, Jr. IARC monographs on the evaluation of carcinogenic risk of chemicals to humans. Am Ind Hyg Assoc J 1980, 41(5):A26, A28, A30 passim.
Helander A, Lindahlkiessling K. Increased frequency of acetaldehyde-induced sister-chromatid exchanges in human-lymphocytes treated with an aldehyde dehydrogenase inhibitor. Mutat Res. 1991;264(3):103–7.
Dellarco VL. A mutagenicity assessment of acetaldehyde. Mutat Res. 1988;195(1):1–20.
Obe G, Jonas R, Schmidt S. Metabolism of ethanol in vitro produces a compound which induces sister-chromatid exchanges in human peripheral lymphocytes in vitro: acetaldehyde not ethanol is mutagenic. Mutat Res. 1986;174(1):47–51.
Muto M, Hitomi Y, Ohtsu A, Shimada H, Kashiwase Y, Sasaki H, Yoshida S, Esumi H. Acetaldehyde production by non-pathogenic Neisseria in human oral microflora: implications for carcinogenesis in upper aerodigestive tract. Int J Cancer. 2000;88(3):342–50.
Nosova T, Jousimies-Somer H, Jokelainen K, Heine R, Salaspuro M. Acetaldehyde production and metabolism by human indigenous and probiotic lactobacillus and Bifidobacterium strains. Alcohol Alcohol. 2000;35(6):561–8.
Munoz-Gonzalez I, Thurnheer T, Bartolome B, Moreno-Arribas MV. Red wine and oenological extracts display antimicrobial effects in an oral bacteria biofilm model. J Agric Food Chem. 2014;62(20):4731–7.
Gerhauser C. Beer constituents as potential cancer chemopreventive agents. Eur J Cancer. 2005;41(13):1941–54.
Preedy VR. Beer in health and disease prevention. Cambridge: Academic Press; 2011.
We thank our colleague, Dr. Yilong Zhang, from the Division of Biostatistics, Department of Population Health, NYU School of Medicine, who provided insight and expertise that greatly assisted the analysis of data.
Research reported in this publication was supported in part by the US National Cancer Institute under award numbers R01CA159036, U01CA182370, R01CA164964, R03CA159414, P30CA016087, and R21CA183887, and by AACR/Pancreas Cancer Action Network Career Development Award.
Availability of data and materials
The data are available from the Sequence Read Archive (SRA) with accession number SRP133146 and SRP133149.
The American Cancer Society (ACS) funds the creation, maintenance, and updating of the Cancer Prevention Study II cohort. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health, the US Department of Veterans Affairs, or the US Government.
Ethics approval and consent to participate
The studies described in this manuscript were approved by the institutional review boards (IRB) of the NYU School of Medicine, and all participants provided written informed consent (IRB study number: S12-00721 and 10-01834).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Tables S1-S4 and Figures S1-S5. Table S1. Number of sequence reads per sample assigned to reference taxonomy map. Table S2. Mean counts of taxa associated with drinking levels in all study participants and subjects in each cohort. Table S3. Associations between oral bacteria diversity and abundance of order Lactobacillales, genera Actinomyces and Neisseria with alcohol drinking levels in sensitivity analysis. Table S4. Mean counts and fold changes of taxa associated with type of alcoholic beverages in all study participants. Figure S1. Beta-diversity of oral bacterial communities in the paired mouthwash-saliva samples. Figure S2. Distribution of alcohol consumption in the enterotypes. Figure S3. Count boxplots of order Lactobacillales and selected genera that were associated with drinking level. Figure S4. Partial Constrained Analysis of Principal Coordinates plots. Figure S5. Richness and evenness of oral microbiome by alcohol drinking type. (DOCX 1148 kb)