Skip to main content

Oral ecological environment modifications by hard-cheese: from pH to microbiome: a prospective cohort study based on 16S rRNA metabarcoding approach



The oral ecosystem conditions dental health, and is known to be positively modified by oral hygiene which cannot always be performed between meals, especially outside home. It is therefore important to identify the practices to be adopted to influence the oral environment in an anticariogenic direction. Milk and cheese are considered functional foods and have a role on oral health. There are several mechanisms by which cheese exerts its beneficial effects on teeth. The aim of the present study was to examine whether short term consumption of hard cheese would affect the oral pH and microbial flora of healthy adults modifying ecological oral environment. The Next Generation Sequencing (NGS) approach was applied to study the effect of Italian Grana Padano (GP), as a prototype of typical hard cheese, on the oral microbiota composition. Finally, we explored Streptococcus mutans/sanguinis ratio as a marker of protective biofilm composition.


Nine oral-healthy adults were instructed to eat 25 gr of GP cheese for 5 consecutive days. Three time points were chosen for supragingival samples collection and pH measurement. 16S rRNA-gene sequences were obtained both from oral samples and GP cheese using the MiSeq platform and analyzed against the expanded Human Oral Microbiome Database (eHOMD). ProgPerm was used to perform statistical analyses to investigate strain differential representation after cheese consumption.


Taxonomic analyses of the oral microbiota revealed that Firmicutes was the most abundant phylum, followed by Proteobacteria and Actinobacteria. GP cheese significantly modifies oral pH, causing a shift toward basic conditions which are kept for a few hours. The Streptococcus mutans/Streptococcus sanguinis ratio lowers in the last observed timepoint.


Our results reveal that a portion of GP cheese eaten after dinner provides important micronutrients (i.e. calcium, vitamins and some aminoacids such as arginine) and changes oral pH toward basic conditions, resulting in a light modification of the oral microbiome towards the reduction of the overall amount of acidophilic bacteria. Furthermore, the S. mutans/S. sanguinis ratio is reduced, contributing to obtain a more protecting environment towards caries establishment and evolution.


Dental caries is the most prevalent noncommunicable disease (NCD), that affects all age groups across the life course all over the world (49) with increasing prevalence in many low-income and middle-income countries [1]. It affects about 95% of the world child population, while periodontal disease, common in the elderly, affects about 15% of the population [2]. Untreated caries in permanent teeth was the most prevalent health condition in 2010, affecting 35% of the global population, or 2.4 billion people worldwide [3]. Changed eating habits, lifestyle and oral hygiene during the COVID-19 pandemic, due to lockdown, wearing mask and the reduced access to dentistry care, increased the risk of developing tooth decay [4].

The human oral cavity is colonized by a microbial community that plays an important role in regulating the oral health status [1, 3, 5,6,7,8]. Oral diseases, in particular dental caries and periodontitis, develop as a result of complex perturbations of oral homeostasis. Environmental changes in the oral cavity, including the predominantly acidic or basic pH to which the oral cavity is exposed, affect the microbial community, consequently, it is fundamental to understand the composition of the oral microbial community and the conditions that maintain it in ecological balance to prevent and manage the progression of disease. Food influence is of paramount importance, and dairy foods have long been well known for their role in prevention of dental caries [8]. Several studies have shown anticariogenic properties of milk and cheese [9]. Moreover, recent epidemiological studies in children and adolescents have associated the frequency of milk and cheese consumption to a better oral health. Some of the anticariogenic properties of these dairy foods may be attributable to calcium, phosphate, casein, and postulated mechanisms involve buffering, salivary stimulation, reduction of bacterial adhesion, reduction of enamel demineralization, and/or promotion of remineralization by casein and ionizable Ca, P and casein phosphopeptides [10].

Grana Padano (GP) is one of the most popular Italian hard cheeses. Its typical features are the result of a specific processing technique, a long ripening, and the quality of milk produced by dairy cows fed a specific diet. GP contains a high concentration of proteins and essential aminoacids: leucine, lysine, phenylalanine and tyrosine, to name a few. A very high arginine content makes it a great resource for ammonia producing bacteria, which buffer the oral acid environment. Additionally, being a dairy product, GP supplies proteins and bioactive peptides in a concentrated form. According to CREA (Consiglio per la Ricerca in agricoltura e l'analisi dell'Economia Agraria) tables of food composition (, 100 gr of GP PDO contain 1165 mg of calcium and high phosphorus and magnesium concentrations which are important for the development of bone and teeth in children, as well as for their maintenance in all segments of the population. GP is a good source of casein phosphopeptides, known to be essential for teeth remineralization (for nutritional information see Additional file 1: Table S1). Cheese is also an appropriate delivery vehicle for probiotic bacteria [11]. Lastly, a very notable feature of GP as well as other aged hard cheeses resides in the fact that it is lactose-free. Hard cheese consumption is widely spread across the world [12].

Next Generation Sequencing (NGS) showed that oral microbiota is a complex habitat where hundreds or thousands of microbial species co-inhabit the oral cavity and functionally interact [13]. NGS approach has identified more than 700 different predominant bacterial species in the oral cavity, out of which 35% have not been cultivated, yet [14]. However, recently a novel workflow based on targeted phenotypic culturing linked to large-scale whole-genome sequencing, phylogenetic analysis and computational modeling has been described to study the functions and phenotypes of several bacteria believed uncultivable. The application of this model will improve the knowledge about the role of a larger number of bacteria in the whole microbiota community [15]. Several investigations have shown that dental caries is not caused by a single agent, but is associated to an imbalance of microbial species that synergistically cause enamel demineralization by their acidogenic activity [16, 17]. A specific symbiotic relationship between the resident oral microbiota and the host is believed to be essential for maintenance of oral health.

The effect of diet and/or functional food on oral microbiota remains largely unexplored. In the absence of a specific diet, the composition of the salivary microbiota and oral plaque shows no diurnal variation and tends to remain stable for long periods [18, 19]. The aim of this work is therefore to investigate the potential effect of a typical hard cheese consumption on oral health using NGS technologies, by observing relevant changes in microbiome composition across timepoints. We characterized both cheese and oral microbiota to check if cheese could directly transfer its bacteria to the mouth or induce modifications in oral microbiota through other indirect mechanisms.

Material and methods

Patients and samples

A prospective cohort study [20] was carried out on nine healthy adult volunteers selected among the residents attending the University of Milan, Department of Biomedical, Surgical and Dental Sciences, Fondazione IRCCS Cà Granda, Ospedale Maggiore Policlinico Milan (Italy), following the STORMS reporting guidelines for human microbiome research [21]. All the subjects were periodontal healthy, according to the International Workshop for a Classification of Periodontal Diseases and Conditions criteria [22], free of caries and without lesions of the mucosa. Exclusion criteria included: diabetes; pregnancy; human immunodeficiency virus (HIV) positivity; use of immunosuppressant medications; antibiotic therapy or oral prophylactic procedures within the past 3 months and during the study; use of probiotics in the period 3 months prior to the study and during it; need for antibiotic coverage before dental treatment; presence of periodontal or gingival diseases; fewer than 28 teeth present in the dentition.

All subjects were instructed to eat a single portion (25 gr) of Italian hard cheese GP after dinner for five consecutive days and to keep a food diary. After the first day (i.e. after the first dinner), participants were asked to abstain from oral hygiene. The next morning, before tooth brushing, they collected a sample of supragingival plaque from the first and second mandibular molar. During the remaining days of the study all subjects brushed regularly their teeth three times per day, using commercially available fluoride toothpaste. All subjects were instructed to refrain from smoking, drinking (except for water) for 8 h, tooth brushing for 3 h and drinking water for 1 h prior to sample collection.

Oral samples were collected using the ORAcollect•DNA kit (DNA-Genotek Inc.) according to the kit protocol. After sampling, samples were immediately stored at − 80 °C. Three samples (A, B and C) of supragingival plaque were obtained from each patient at three different times (Fig. 1): Sample A = T0: baseline with no specific pre-treatment; Sample B = T1: the morning after the first consumption of cheese (no teeth brush before sampling); Sample C = T2: the morning after the fifth day of regular consumption of cheese (no teeth brushing before sampling). pH levels were recorded at T0, and after six timepoints (10, 30, 60, 120, 240 and 480 min) after GP cheese assumption.

Fig. 1
figure 1

Scheme of sampling. On the top of the picture, the time of the cheese meals (lunch, dinner, etc.) and/or the oral hygiene actions (a red x on the toothbrush means that the teeth were not washed). Yellow arrows indicate the time of the day when cheese was eaten and when the teeth were washed. The Blue arrows indicate the buccal swab sampling collected at T = 0, T = 1 and T = 2. Samples (A), (B) and (C) are, respectively, the samples at different times

Anthropometric analysis and evaluation of food intake

Anthropometric measurements were collected on the nine healthy volunteers, according to the conventional criteria and measuring procedures proposed by Lohman [19]. Body Weight (BW, Kg) was measured by Column scale (Seca 700 balance, Seca Corporation, Hanover, MD, USA) to the nearest 100 g with subjects wearing only light underwear and after bladder emptying, Body Height (BH, cm) was measured to the nearest 0.1 cm using a vertical stadiometer. Body mass index (BMI) was calculated using the formula: BMI = Kg/m2. Volunteers were invited to record their food consumption for five consecutive days using a 24-h food diary to estimate the dietary energy, micro and macronutrients intakes. Full written instructions on how to complete the diary were provided. Diaries were analyzed by specific software that linked diet to a comprehensive database of macro- and micronutrients (Dietosystem, DS Medica, Milan). The estimated daily intake was compared to the dietary reference values for adult standard portions in Italian Diet according to the INRAN guideline (Guidelines for a healthy Italian diet. Italian National Institute for Research on Food and Nutrition) [23].

pH measurement before and after consumption of a standard portion of hard cheese

After determining the resting salivary pH using pH test strips, the subjects were asked to eat the test food (hard cheese) as previously described.

Salivary pH was measured at time points 0, 10, 30, 60, 120, 240 and 480 min, before eating or drinking or brushing teeth (Saliva-Check BUFFER: D0425-Caries Susceptibility Test-GC America 0.2 resolution).

Unstimulated saliva from each subject, at every single time point, was collected in a plastic disposable container. Immediately after collection, salivary pH was measured using pH strips by dipping the strip into test solution until color development was complete and comparing to the sequence chart on the package to read the pH.

DNA extraction, library preparation and sequencing

Bacterial DNA was extracted from plaque samples using prepIT®•L2P according to the manufacturer’s instructions. The V3-V4 region of the 16S rRNA was amplified using 16S Metagenomic Sequencing Library Preparation kit (Illumina Inc.). In brief, a PCR amplicon of about 550 bp was obtained with V3 and V4 primer pairs included in the kit. The PCR products were purified with AMPure XP beads and dual indices, and Illumina sequencing adapters were added using the Nextera XT Index Kit. AMPure XP beads were used to clean up the final library before quantification. Nextera libraries were quantified by qPCR using Brilliant III SYBR Green qPCR Master Mix (Agilent Tech.) and the size profile was analyzed on a 2200 TapeStation (Agilent Tech). Libraries were then adjusted to a concentration of 5 µg/ml. Equal volumes of normalized DNAs were then pooled and run on MiSeq (Illumina Inc.).

DNA was extracted from the Italian GP cheese using the kit NucleoSpin soil (Macherey–Nagel). Library preparation and sequencing were carried out as for oral samples.

Sequence processing and analysis

Raw reads quality was checked using FastQC v0.11.2. Illumina raw sequences were trimmed using Trimmomatic v0.32 with default parameters. Minimum base quality 20 (Phred-scale) over a 4 bases sliding window was required. Only sequences longer than 36 bps were included into downstream analysis. A Qiime formatted taxonomy version of the expanded Human Oral Microbiome database (eHOMD 16S rRNA RefSeq Version 15.22) [24], considering 97% identity to identify cluster, was used to identify the 7 main taxa ranks (superkingdom, phylum, class, order, family, genus, species). For each sample raw data were aligned against the reference database using BWA MEM v0.7. Quantification was performed using SAMTOOLS. The OTU (Operational Taxonomic Unit) table was created using custom scripts. OTU-table in BIOM file format was generated using Qiime Biom convert. Specific primers were built using reference sequences for Rothia dentocariosa, as a proxy for the most relevant species, in order to adjust for relative quantification in the meta barcoding analyses. Prior to quantitative analysis, OTUs from the same genus with unclassified species were grouped, and OTUs with counts in less than one third of the samples were removed.

For oral samples, multidimensional scaling (MDS) was performed with the EdgeR R package ( to explore similarity in microbiota composition in different groups. ProgPerm [25] was used to study both the overall and the individual associations, and to provide measures on the robustness of the discoveries and the reliability of the results. The method progressively permutes the grouping factor labels of microbiome samples and performs differential testing by applying a Wilcoxon rank-sum test to the permuted data in each scenario. The signal strength (− log10 p-values) of top hits from the observed data is compared with their testing performance in permuted data sets to discriminate between significant and nonsignificant hits.

Arginolytic bacteria identification in the oral microbiome

To identify potentially arginolytic bacteria present in our dataset, we checked for the presence of the arginine deiminase arcA (alias sagP) gene in organisms included in the AnnoTree database [26] and matched the retrieved bacteria with those present in our dataset.


Nutritional status

Five females and four males, with mean age 36 ± 4 years participated in the trial. Their macro- and micronutrients intake was adequate in relation to daily recommendation, apart from calcium below RDA, in agreement with the general population data (Table 1).

Table 1 nutritional status. Data are shown as daily mean ± standard deviation

Influence of hard cheese assumption on salivary pH levels

pH levels were recovered just before eating a standard portion of GP cheese (T0) and at six time points (10, 30, 60, 120, 240 and 480 min) after GP cheese assumption. In the first 10 min salivary pH showed a significant increase, reaching a peak from baseline. Thereafter, the pH showed a small progressive and slow decline to baseline after more than 4 h (Fig. 2).

Fig. 2
figure 2

pH measurements at the seven recovered timepoints with standard deviation among samples

GP cheese microbiota composition

A total of 342,166 reads were obtained from the GP cheese sample, representing 329 OTUs after grouping OTUs from the same genus with unclassified species (Additional file 2: Table S2). Out of these, 126 showed ≥ 10 reads. OTUs with abundance higher than 1% of the total represented OTUs were all Lactobacillus species, Lactobacillus paracasei, casei and rhamnosus accounting for more than 63%, in total (Fig. 3). At the Phylum level, Firmicutes represent the majority of the dataset (96,96%); followed a long way behind by Proteobacteria and Actinobacteria (1,66% and 1,20%, respectively).

Fig. 3
figure 3

Species composition of GP cheese. Species represented by less than 1% OTUs were grouped

Microbiota profiling in oral samples

A mean of 758,441 reads per sample (standard deviation 213,584, min 131,797, max 1,279,983) was obtained. eHOMD version 15.22 contains 1015 OTUs, 535 of which were recovered in this investigation, after grouping OTUs from the same genus with unclassified species (Additional file 3: Table S3). The percentage of different phyla, classes and orders in groups A, B and C is reported in Fig. 4.

Fig. 4
figure 4

Pie chart showing the percentage of the major Phyla, Class and Order identified in samples A (T0), B (T1), and C (T2) by 16S rRNA analysis

The most abundant phylum detected in all groups is Firmicutes (58%, 55%, 53% in A, B, and C, respectively), followed by Proteobacteria, Actinobacteria, Bacteroidetes, and Fusobacteria. Bacilli are the most abundant class (51% in A, 45% in B, and 47% in C), followed by Actinobacteria and Gammaproteobacteria. Lactobacillales is the most abundant order (44% in A, 37% in B, and 41% in C). The number of OTU in the full lineage at each timepoint is reported in Additional file 4: Table S4. As a control, specific primers for Rothia dentocariosa, one of the most represented species in all datasets, were used to validate the numeric representation of reads, confirming data obtained with the bioinformatic procedure leading to read counts.

The MDS Plot based on the distribution of all species (Fig. 5) shows that microbiota composition reflects a high heterogeneity across individuals. In fact, samples collected from a same person at the three time points cluster much better than samples from different individuals collected at the same time point.

Fig. 5
figure 5

Multidimensional Scaling (MDS) plot of microbiota of different samples taken at different time points A (black), B (red) and C (green)

Differential abundance analysis

Differential abundance analysis on microbiome data is challenging due to many factors, the main of which are clustering sequence reads into OTUs and data heterogeneity across subjects. The degree of heterogeneity across samples can lead to the false conclusion that the microbiome is not associated with an experimental condition, or that some species are differentially represented across biological datasets, as the effects of the applied condition are masked among differences due to the personal microbiome composition, which is highly variable. This is particularly notable when experimental conditions interfere with quantitatively small effects on microbiome composition, as is the case of our tests. It is thus fundamental to explore both the overall and the individual associations to highlight significant differences between groups of subjects who underwent the same treatment, and to apply statistical pipelines that test the robustness of the highlighted associations.

To analyze our data, we chose to apply the Qiime pipeline for sequence clustering according to OTUs collected into eHOMD v15.22, the most updated and specific repository for oral microbiome representation, and to use ProgPerm for data abundance analysis. This method performs progressive permutation and applies correlation tests to highlight significant differences in OTU representation between experimental conditions. We performed comparisons between groups using all possible combinations, i.e. T0 vs T1, T0 vs T2, and T1 vs T2. 530 out of the 535 OTUs recovered after sequence clustering passed the cutoff of exhibiting counts in at least one third of the samples and were analyzed to retrieve abundance differences.

The T0 vs T1 comparison is aimed at individuating short term microbiome modifications induced by GP cheese consumption. The U-Curve obtained with ProgPerm (Fig. 6A), which can be used as a global measure to depict the overall association between microbiome compositions and different experimental outcomes, shows a light association in the observed data, suggesting a moderate influence of GP cheese eating on oral microbiome composition. This was expected, as we’re observing the effects of a single assumption, which cannot be compared to the effects of more influential long-term treatments that recover differences due to microbiome composition recalibration after long periods of administration of a particular diet, therapy, or drug. In such a condition, the variable’s effect might be modified by other variables and distorted by potential systematic bias, thus, highlighting significant differences between timepoints becomes particularly challenging. A similar behavior can be observed for the T0 vs T2 comparison (Fig. 6B), which aims at studying a more extended effect of eating cheese on the oral microbiota composition (5 days of administration). The T1 vs T2 comparison (Fig. 6C) aims at highlighting differences in microbiome composition due to the effects of GP extended assumption, starting from a condition where GP assumption has already imposed its modifications to the native microbiome composition. In this comparison, the U-Curve denotes the absence of significant variations in microbiome composition between the two datasets.

Fig. 6
figure 6

ProgPerm U-Curve of proportion of significant features versus proportion of mixing of the T0 vs T1 (A), T0 vs T2 (B), and T1 vs T2 (C) comparisons. The x-axis describes the proportion of mixing the two groups of data. The y-axis describes the proportion of significant features. The red dot describes the observed data, black dots describe the permuted data. The vertical bars describe the 95% quantile confidence intervals

Despite the low association denoted by the U-Curves, OTUs with significant differential abundance were individuated in all three comparisons. ProgPerm achieves the identification of robust variables by observing whether the − log10 p-values of targeted features lie within the 95% confidence interval of median − log10 p-values of the full permutation scenario. According to this classification, the T0 vs T1 comparison shows 13 species which are significantly more represented in T1 (Fig. 7A). 12 hits are more represented in T2 within the T0 vs T2 comparison (Fig. 7B). Linear discriminant analysis effect size (LEfSe)—determining the features most likely to explain differences between classes by coupling standard tests for statistical significance with additional tests encoding biological consistency and effect relevance—is also computed by ProgPerm (Fig. 7D and E). In both comparisons, Parvimonas micra—an anaerobic coccus frequently isolated from dental plaque especially in patients with chronic periodontitis, known to harbour aminopeptidase activity—is the species showing the highest effect size, and Parvimonas sp. and Selenomonas flueggei are also recovered. These three bacteria exhibit a significant increment in representation moving from T0 to T1, and T2 levels are very similar to T1 ones. Among the other statistically significant variations, those related to the Prevotella genus, classified as a proteolytic/amino acid-degrading bacteria, are evident in the T0 vs T1 comparison, with four species (marshii, bivia, nigrescens and oris), as well as Alloprevotella sp., overrepresented in T1. Apart from the already mentioned organisms, species recovered in the T0 vs T2 comparison differ from those recovered in T0 vs T1, with Bulleidia extructa—another anaerobic bacterium associated with oral infections, mainly periodontal disease [27]—exhibiting differential representation between the two timepoints, being more represented in T2. Bulleidia extructa is the only species with a significant effect size in the T1 vs T2 comparison (Fig. 7C and F), thus showing that its increment in representation occurs between T1 and T2, after a more prolonged assumption of GP cheese. In none of the three comparisons species being more represented in the former of the timepoints are recovered, and no association is displayed for classical saccharolytic bacteria, such as Streptococcus, Actinomyces, and Lactobacillus species, known to promote dental caries.

Fig. 7
figure 7

AC: P-value coverage plot of the top 50 features in the T0 vs T1, T0 vs T2, and T1 vs T2 comparisons with decreasing order. The color dots denote the − log10 p-value of top 50 features in the original data (permutation proportion is 0). The horizontal bars describe the 95% quantile confidence intervals of the − log10 p-value in the full permutation scenario. DF Effect size of the identified features in the T0 vs T1, T0 vs T2, and T1 vs T2 comparisons

Trend of caries- or periodontitis-related bacteria across timepoints

Despite not having been individuated by ProgPerm as statistically significant with respect to their quantitative representation across timepoints, the relative abundance of known cariogenic bacteria such as Streptococcus oralis, Streptococcus mutans, Streptococcus sanguinis, Streptococcus parasanguinis, Lactobacillus crispatus and Propionibacterium acidifaciens was observed in the three sample collections (Fig. 8). None of these bacteria exhibits a defined trend across timepoints, not even assuming less stringent cutoffs for statistical significance: the only observable differences are mainly due to samples heterogeneity within timepoints (see Additional file 3: Table S3). In all timepoints, Streptococcus oralis subsp. dentisani and subsp. tigurinus are highly represented, followed by Streptococcus oralis subsp. oralis and Streptococcus sanguinis.

Fig. 8
figure 8

Quantification of known cariogenic bacteria within timepoints

Streptococcus mutans has a demonstrated positive association with the establishment of severe early childhood caries (S-ECC) [28]. Streptococcus sanguinis was proposed to play an antagonistic role against S. mutans colonization in the oral cavity, thus reducing the risk of caries establishment [29, 30]. The two bacteria representation within the oral flora is extremely different, being S. sanguinis far more represented than S. mutans. This hampers the interpretation of the effects of their representation with ordinary data analysis techniques, but the observation of their ratio in different clinical situations showed that lower ratios are associated with a lower incidence of dental caries, suggesting that this ratio could be one of the indicators prone to be used to predict caries establishment and evolution in children [30]. In our study, this ratio doesn’t undergo significant changes from T0 to T1, while T2 exhibits a lower ratio with respect to the former timepoints (Fig. 9).

Fig. 9
figure 9

Streptococcus mutans/Streptococcus sanguinis ratio at the three timepoints

The trend of relative abundance of periodontal bacteria such as Porphyromonas_gingivalis, Aggregatibacter actinomycetemcomitans, Tannerella forsythia and different strains of Treponema denticola was also verified across timepoints (see Additional file 3: Table S3). As for cariogenic bacteria, the relative abundance of these strains was not statistically significant, but the trend shows that all of them tend to be more represented from T0 to T1 and then to T2, except for Aggregatibacter actinomycetemcomitans, which has its lowest representation in T1.

Behaviour of arginolytic bacteria across timepoints

Being GP an arginine-rich substrate, we investigated the hypothesis that arginolytic bacteria present in the oral cavity may be influenced in their representation by its consumption. In our dataset, we found 24 species having the arcA gene, of which proportional representation in all samples, at the three timepoints, is displayed in Fig. 10. Unfortunately, we do not dispose of expression data, and this prevents making any assumption on the activity of the arcA gene of these bacteria at the different timepoints. No common trend can be evidenced for arginolytic bacteria as a response to increased arginine availability. None of them exhibits statistically significant differences in representation across timepoints.

Fig. 10
figure 10

Proportional representation of arginolytic bacteria in the T0, T1 and T2 datasets. Representation in the three datasets was normalized according to the total number of sequences obtained for each timepoint


Oral microbiota and pH have an important role on dental health. Hence, an appealing hypothesis is that some foods may change the oral microbiota composition favoring the presence of bacteria with a positive effect on oral health. The anticariogenic effect of some food, such as hard cheese, and their ability to modify oral microbiota was not previously studied especially by using a 16S rRNA metabarcoding approach.

The present study focuses on the exploration of GP cheese-dependent short term (1 and 5 days after eating cheese) variation of the oral environment both in terms of pH and microbiota. To extrapolate the effects of presence of cheese microbiota versus the ability of cheese to induce modification of the oral microbiota, microbiota composition of the GP cheese was also characterized by NGS. At a taxonomic level, the percentage of different supragingival plaque bacteria identified in this study before (T0) and after GP consumption (T1 and T2) agrees with previous studies [31,32,33]. There were small differences at phylum, class and order level in T0 vs T1 and T2, indicating that eating cheese does not cause a big shift in the oral microbiota in a relatively short time span. Several reports showed stable oral microbiome with little or no variation in profiles within subjects over time [19, 32,33,34]. The stability of the microbiota composition is indispensable to study the effect of functional food on oral microbiota variation.

In our datasets, Streptococcus oralis subsp. dentisani and subsp. tigurinus are the most abundant OTUs in T0, T1 and T2. Streptococcus sp., Haemophilus_parainfluenzae, Rothia_dentocariosa, Gemella_haemolysans, and Streptococcus_infantis are also recovered within the 10 most represented OTUs at all the three timepoints, while Rothia aeria and Veillonella dispar are present within the top 10 most represented OTUs only in two out of three datasets (T0 and T2 the former, and T1 and T2 the latter). Lactobacillus paracasei, casei, rhamnosus and ultunensis, which represent almost 73% of the total OTUs detected in GP cheese, do not increase their representation moving from T0 to T2, suggesting that no direct bacteria transfer from cheese to the oral microbiota can be postulated.

Other abundant OTUs within the top 30 belong to the Streptococcus genus. In healthy individuals, Streptococci constitute more than 50% of the oral microbiota [35]. These bacteria generally possess low pathogenic potential. Streptococci such as Streptococcus sanguinis seem to have a positive effect on oral health: S. sanguinis is often found in subgingival biofilm that includes periodonto-pathogens and is correlated with a delay in colonization by these deleterious bacteria [36]. Streptococcus mutans is known as a cariogenic microorganism with the ability to colonize teeth, to cause a marked reduction in pH in the presence of a sugar substrate and consequently induce caries [37]. The ratio between S. mutans/S. sanguinis has been suggested to be helpful in determining caries risk. A low ratio would imply a low risk of caries and conversely, a high ratio would indicate high caries risk [28,29,30]. In this study few OTUs of Streptococcus mutans were detected, this is expected as this organism is always present in a limited proportion in oral samples, furthermore, all participants to the study were cavities' free. No conserved trend was observed among samples for Streptococcus mutans representation during timepoints, suggesting that GP assumption doesn’t influence its participation to the oral microbiome, nonetheless, the S. mutans/S. sanguinis ratio decreases in T2 when compared to T0 and T1, suggesting the establishment of a more protective environment with respect to caries development. In our experimental trial Parvimonas micra (P. micra) shows the highest effect size both in T0 vs T1 and T0 vs T2, due to a significant increase in representation from T0 to T1, while T1 and T2 show comparable levels. P. micra is normally found in the oral cavity and is one of the bacterial species most frequently isolated from infected root canals of teeth with chronic apical periodontitis.

A previous study showed a disparity between microbial communities localized in acidic versus neutral pH strata in dentinal caries [38]. Acidic conditions were associated with low diversity microbial populations, dominated by Lactobacillus species. In our study, 25 Lactobacillus species are recovered, which do not show a common trend during timepoints, being some of them incremented and others reduced moving from T0 to T2. Nonetheless, the total amount of Lactobacillus decreases during timepoints, passing from 0.266% of the total sequences mapped to known OTUs in T0, to 0.064% in T1 and 0.042% in T2, suggesting a general effect of the increase of pH due to GP cheese assumption on acidophilic bacteria. In our study, almost all Prevotella and Megasphaera species increased their representation after GP cheese assumption mostly in T1. Treponema species maintained this trend up to T2. Prevotella gram negative species were associated to dental caries [39], and some Prevotella spp. such as Prevotella intermedia, Prevotella nigrescens and Prevotella melaninogenica predominate in periodontal disease and periodontal abscesses [40]. However, Prevotella intermedia was also associated to oral health [31]. Treponema sp was differentially detected in healthy controls vs patient cases [36]. Moreover, Treponema was associated with dental health, whereas Lactobacillus, Streptococcus and Megasphaera to dentin lesions [41]. On the contrary, in other studies Treponema denticola was associated with periodontal diseases such as early-onset periodontitis, necrotizing ulcerative gingivitis, and acute pericoronitis [38]. Megasphaera was less frequently found in the caries group [42].

The overall information available on caries establishment and progression, as well as susceptibility to periodontitis, clearly shows that both pathologies cannot be ascribed to a few well identified pathogens, but depend on the interaction of multiple factors and organisms. Oral microbiota is a mirror of the state of health of an individual, and as it can be strongly affected by diet, it is important to determine which foods can be of help in maintaining a balanced and healthy composition. Much attention has been focused on the production of organic acids by oral pathogenic bacteria through the metabolism of carbohydrates. Emerging evidence indicates that alkali generation by bacteria, particularly through ammonia production from arginine and urea, plays a major role in pH homeostasis of oral biofilms and may moderate the initiation and progression of dental caries. Dental caries occurs when the acidification phases outweigh the alkalization phases, hence allowing the establishment of a more acidogenic, less alkalinogenic flora, which in turn results in lower plaque pH values with enhanced and prolonged enamel demineralization. As demonstrated by Stephan’s Curve [43], described by Robert Stephan in 1943 and reflecting plaque pH change following intake of fermentable carboydrates, acid liquids, or sugars in presence of acidogenic bacteria, sugar assumption causes a fall of the pH level under the threshold of 5.5, defined as critical because corresponding to the condition in which the demineralization of the enamel takes place. To reduce decalcification of tooth surfaces, patients should be encouraged to consume foods as hard cheese that do not result in a drop in plaque pH. A very high arginine content makes GP cheese a great source for ammonia producing bacteria which can buffer the oral acid environment. The pH curve shown in this study, in fact, has an opposite trend with respect to Stephan’s curve, showing a significant pH increase in the short term. Additionally, being a dairy product, GP supplies proteins and bioactive peptides as milk does, but in a concentrated form.


This study examines, for the first time by NGS technology, the short-time effects of eating a portion of hard cheese on the oral microbiome composition and the associated pH variations. NGS analysis of the bacterial composition of GP cheese, conducted in parallel, reveals that there is no direct transfer of microorganisms from cheese to mouth, and that other molecular mechanisms are responsible for microbiota modifications. Even if no direct evidence of a significant modification in the representation of bacteria known to be cariogenic or directly related to periodontitis can be provided, small changes in the amounts of acidophilic bacteria on their whole was observed, together with a decrease in the S. Mutans/S. sanguinis ratio, which mainly occurs as an effect of eating GP routinely. Although these trends require to be further investigated, they suggest a positive effect on the overall composition of the oral flora towards a healthy status which can lead to a major resistance to caries establishment. Moreover, we have observed that the variations in the alkaline sense that occur in the mouth immediately after the consumption of cheese are maintained over time, up to over four hours. The alkaline environment, together with calcium, phosphorus, caseinphosphopeptides and arginine, contained in large quantities in this cheese, not only inhibits demineralization, but also promotes the process of remineralization, making this food ideal for determining an oral environment that hinders the formation of caries in the periods of time between one meal and another, especially in the conditions of impossibility of correct oral hygiene practices.

Availability of data and materials

The datasets generated and analysed during the current study are available in the ENA Metagenomics repository, Project PRJEB51725, Samples ERS11067207- ERS11067233 (oral samples), and ERS11068245 (GP cheese).



Grana Padano


Next Generation Sequencing


Operational Taxonomic Unit


  1. Peres MA, Macpherson LM, et al. Oral diseases: a global public health challenge. Lancet. 2019;394(10194):P249-260.

    Article  Google Scholar 

  2. Whelton H. Nutrition and health-nutritional and oral health-promoting properties of dairy products: caries prevention and oral health. Encycl Dairy Sci. 2011: 1034–40.

  3. Kassebaum NJ, Bernabé E, Dahiya M, Bhandari B, Murray CJ, Marcenes W. Global burden of untreated caries: a systematic review and metaregression. J Dent Res. 2015;94(5):650–8.

    CAS  Article  Google Scholar 

  4. Docimo R, Costacurta M, Gualtieri P, Pujia A, Leggeri C, Attinà A, Cinelli G, Giannattasio S, Rampello T, Di Renzo L. Cariogenic risk and COVID-19 lockdown in a paediatric population. Int J Environ Res Public Health. 2021;18(14):7558–76.

    CAS  Article  Google Scholar 

  5. Marsh PD. Dental plaque as a biofilm and a microbial community -implications for health and disease. BMC Oral Health. 2006;6(Suppl 1):S14.

    Article  Google Scholar 

  6. Haffajee AD, Socransky SS. Microbiology of periodontal diseases: introduction. Periodontology. 2000;2005(38):9–12.

    Google Scholar 

  7. Haffajee AD, Socransky SS. Introduction to microbial aspects of periodontal biofilm communities, development and treatment. Periodontology. 2000;2006(42):7–12.

    Google Scholar 

  8. Paster BJ, Olsen I, Aas JA, Dewhirst FE. The breadth of bacterial diversity in the human periodontal pocket and other oral sites. Periodontology. 2000;2006(42):80–7.

    Google Scholar 

  9. Johansson I. Milk and dairy products: possible effects on dental health. Scand J Nutr. 2002;46:119–22.

    Article  Google Scholar 

  10. Aimutis WR. Bioactive properties of milk proteins with particular focus on anticariogenesis. J Nutr. 2004;134(4):989S-S995.

    CAS  Article  Google Scholar 

  11. Biziulevicius GA, Kazlauskaite J. Cheese: everyday, gourmet or medicinal food? Med Hypotheses. 2008;70(2):454–5.

    Article  Google Scholar 

  12. Accessed 08 Mar 2022.

  13. Bik EM, Long CD, Armitage GC, Loomer P, Emerson J, Mongodin EF, Nelson KE, Gill SR, Fraser-Liggett CM, Relman DA. Bacterial diversity in the oral cavity of 10 healthy individuals. ISME J. 2010;4(8):962–74.

    Article  Google Scholar 

  14. Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, et al. The human oral microbiome. J Bacteriol. 2010;192:5002–17.

    CAS  Article  Google Scholar 

  15. Browne HP, Forster SC, Anonye BO, Kumar N, Neville BA, Stares MD, Goulding D, Lawley TD. Culturing of “unculturable” human microbiota reveals novel taxa and extensive sporulation. Nature. 2016.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Marsh PD. Dental plaque as a biofilm and a microbial community - implications for health and disease. BMC Oral Health. 2006;6(Suppl 1):S14.

    Article  Google Scholar 

  17. Wilson M. The oral cavity and its indigenous microbiota. In: Wilson M, editor. Microbial inhabitants of humans. New York: Cambridge University Press; 2005. p. 318–74.

    Google Scholar 

  18. Belstrøm D, Holmstrup P, Bardow A, Kokaras A, Fiehn NE, Paster BJ. Temporal stability of the salivary microbiota in oral health. PLoS ONE. 2016;11(1): e0147472.

    Article  Google Scholar 

  19. Sato Y, Yamagishi J, Yamashita R, Shinozaki N, Ye B, Yamada T, Yamamoto M, Nagasaki M, Tsuboi A. Inter-individual differences in the oral bacteriome are greater than intra-day fluctuations in individuals. PLoS ONE. 2015;10(6): e0131607.

    Article  Google Scholar 

  20. Thiese MS. Observational and interventional study design types; an overview. Biochem Med. 2014;24(2):199–210.

    Article  Google Scholar 

  21. Mirzayi C, Renson A, Genomic Standards Consortium, et al. Reporting guidelines for human microbiome research: the STORMS checklist. Nat Med. 2021;27:1885–92.

    CAS  Article  Google Scholar 

  22. Armitage GC. Development of a classification system for periodontal diseases and conditions. Ann Periodontol. 1999;4:1–6.

    CAS  Article  Google Scholar 

  23. Accessed 08 Mar 2022.

  24. Fernández Escapa I, Chen T, Huang Y, Gajare P, Dewhirst FE, Lemon KP. New insights into human nostril microbiome from the expanded Human Oral Microbiome Database (eHOMD): a resource for species-level identification of microbiome data from the aerodigestive tract. mSystems. 2018.

    Article  Google Scholar 

  25. Zhang L, Shi Y, Do KA, et al. ProgPerm: progressive permutation for a dynamic representation of the robustness of microbiome discoveries. BMC Bioinformatics. 2021;22:126.

    CAS  Article  Google Scholar 

  26. Mendler K, Chen H, Parks DH, Hug LA, Doxey AC. AnnoTree: visualization and exploration of a functionally annotated microbial tree of life. Nucleic Acids Res. 2019;47:4442–8.

    CAS  Article  Google Scholar 

  27. Kloesel B, Beliveau M, Patel R, Trousdale RT, Sia IG. Bulleidia extructa periprosthetic hip joint infection, United States. Emerg Infect Dis. 2013;19(7):1170–1.

    Article  Google Scholar 

  28. Laitala M, Alanen P, Isokangas P, Söderling E, Pienihäkkinen K. A cohort study on the association of early mutans streptococci colonisation and dental decay. Caries Res. 2012;46(3):228–33.

    CAS  Article  Google Scholar 

  29. Lozano CP, Díaz-Garrido N, Kreth J, Giacaman RA. Streptococcus mutans and Streptococcus sanguinis expression of competition-related genes, under sucrose. Caries Res. 2019;53:194–203.

    CAS  Article  Google Scholar 

  30. Ge Y, Caufield PW, Fisch GS, Li Y. Streptococcus mutans and Streptococcus sanguinis colonization correlated with caries experience in children. Caries Res. 2008;42:444–8.

    CAS  Article  Google Scholar 

  31. Chen L, Qin B, Du M, Zhong H, Xu Q, Li Y, Zhang P, Fan M. Extensive description and comparison of human supra-gingival microbiome in root caries and health. PLoS One. 2015;10(2): e0117064.

    Article  Google Scholar 

  32. Utter DR, Mark Welch JL, Borisy GG. Individuality, stability, and variability of the plaque microbiome. Front Microbiol. 2016;22(7):564.

    Google Scholar 

  33. Lazarevic V, Whiteson K, Hernandez D, François P, Schrenzel J. Study of inter-and intra individual variations in the salivary microbiota. BMC Genomics. 2010;28(11):523.

    Article  Google Scholar 

  34. Belstrøm D, Paster BJ, Fiehn NE, Bardow A, Holmstrup P. Salivary bacterial fingerprints of established oral disease revealed by the Human Oral Microbe Identification using Next Generation Sequencing (HOMINGS) technique. J Oral Microbiol. 2016;14(8):30170.

    Article  Google Scholar 

  35. Human Microbiome Project. C Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.

    Article  Google Scholar 

  36. Lee SH. Antagonistic effect of peptidoglycan of Streptococcus sanguinis on lipopolysaccharide of major periodontal pathogens. J Microbiol. 2015;53(8):553–60.

    CAS  Article  Google Scholar 

  37. Loesche WJ. Role of Streptococcus mutans in human dental decay. Microbiol Rev. 1986;50:353–80.

    CAS  Article  Google Scholar 

  38. Kianoush N, Adler CJ, Nguyen KT, Browne GV, Simonian M, Hunter N. Bacterial profile of dentine caries and the impact of pHon bacterial population diversity. Plos ONE. 2014;9: e92940.

    Article  Google Scholar 

  39. Kianoush N, Adler CJ, Nguyen KA, Browne GV, Simonian M, Hunter N. Bacterial profile of dentine caries and the impact of ph on bacterial population diversity. Plos ONE. 2014;9(3): e92940.

    Article  Google Scholar 

  40. Tanaka S, Yoshida M, Murakami Y, et al. The relationship of Prevotella intermedia, Prevotella nigrescens and Prevotella melaninogenica in the supragingival plaque of children, caries and oral malodor. J Clin Pediatr Dent. 2008;32(3):195–200.

    Article  Google Scholar 

  41. Jiang W, Ling Z, Lin X, Chen Y, Zhang J, Yu J, Xiang C, Chen H. Pyrosequencing analysis of oral microbiota shifting in various caries states in childhood. Micob Ecol. 2014;67(4):962–9.

    Article  Google Scholar 

  42. Belstrøm D, Fiehn NE, Nielsen CH, Holmstrup P, Kirkby N, Klepac-Ceraj V, Paster BJ, Twetman S. Altered bacterial profiles in saliva from adults with caries lesions: a case-cohort study. Caries Res. 2014;48(5):368–75.

    Article  Google Scholar 

  43. Bowen WH. The Stephan curve revisited. Odontology. 2013;101(1):2–8.

    Article  Google Scholar 

Download references


The authors thank Paola Cremonesi for her technical support.


This project was financed by The Flagship Project InterOmics (MIUR).

Author information

Authors and Affiliations



ECL, GF, and AS conceived and designed the study. BL analyzed the data, participated in the design of the study and wrote the paper with ECL and AS. FB, SB, GT collaborated to write the manuscript. VL helped in acquisition of samples. AS, PC performed the statistical analysis and collaborated to writing and revising the manuscript. AS, ECL and SB coordinated the work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Barbara Lazzari.

Ethics declarations

Ethics approval and consent to participate

This study was conducted according to the Good Clinical Practice guidelines and was approved by the Human Ethic Committees of Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico (Registration number: 2081). All participants provided written informed consent.

Consent for publication

Consent to publication was obtained by all the participants.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

GP cheese nutritional information.

Additional file 2: Table S2.

GP cheese OTU classification. The table reports the counts for each OTU found in the GP cheese sample. OTUs from the same genus with unclassified species were grouped.

Additional file 3: Table S3.

OTU counts in all samples. Total counts for all samples A (T0), B (T1) and C (T2) are given, together with counts of single samples. OTUs from the same genus with unclassified species were grouped.

Additional file 4: Table S4.

Full lineage of the OTUs recovered in this study. Counts for each species are given as a total for all samples at each timepoint. OTUs from the same genus with unclassified species were grouped.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lorenzini, E.C., Lazzari, B., Tartaglia, G.M. et al. Oral ecological environment modifications by hard-cheese: from pH to microbiome: a prospective cohort study based on 16S rRNA metabarcoding approach. J Transl Med 20, 312 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • pH
  • Oral microbiota
  • 16S rRNA metabarcoding
  • Next generation sequencing
  • Dental caries
  • Periodontitis
  • Cheese