Journal of Nutrition and Metabolism

Journal of Nutrition and Metabolism / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 6629951 |

Teresa Louro, Carla Simões, Wilmara Lima, Laura Carreira, Paula Midori Castelo, Henrique Luis, Pedro Moreira, Elsa Lamy, "Salivary Protein Profile and Food Intake: A Dietary Pattern Analysis", Journal of Nutrition and Metabolism, vol. 2021, Article ID 6629951, 10 pages, 2021.

Salivary Protein Profile and Food Intake: A Dietary Pattern Analysis

Academic Editor: Roxana Valdés-Ramos
Received07 Oct 2020
Revised25 Mar 2021
Accepted30 Mar 2021
Published16 Apr 2021


Saliva research has gained interest due to its potential as a source of biomarkers. One of the factors inducing changes in saliva, in the short term, is food intake, and evidence exist about changes in salivary proteome induced by some food components. Since this topic of research is in its early stages, it was hypothesized that saliva protein composition could be associated with different levels of adherence to dietary patterns that contain higher amounts of plant products. The aim of the present study was to test this hypothesis, in adults, by comparing salivary protein electrophoretic profiles of individuals with different diet characteristics, particularly dietary patterns (DP) that exhibit different proportions of animal and plant-based products. Dietary habits were assessed in 122 adults (61 from each sex, with ages ranging from 20 to 59 years) using Food Frequency Questionnaires. To identify the dietary patterns, a principal component analysis was used. Individual’s non-stimulated saliva was evaluated for flow rate, pH, protein concentration, α-amylase activity, and electrophoretic protein profiles. Seven dietary patterns (DP) were identified. Salivary amylase enzymatic activity was positively associated with animal-based and starchy foods DP, and with plant-based fatty foods without wine DP. At the same time, protein bands containing amylase and type S cystatins were positively associated with the cheese/yoghurt and wine DP. Our results support the association of salivary proteomics and different dietary patterns and highlight the need of considering food consumption habits in studies using saliva, since this is a factor associated with variations in the composition of this fluid.

1. Introduction

Nowadays, due to global climate changes and the rise of food-related diseases (e.g., obesity, cardiovascular diseases, etc.), dietary patterns are becoming a major issue. Different methodologies are used to access dietary habits, most of which rely on subjective reports by individuals, which implies engagement and memory capacities [1], which may result in increased error in the data obtained. Moreover, individuals often misreport dietary intake, expressing what they believe to be adequate eating amounts, this misreport being even more frequent for people with a history of dieting and being overweight [2, 3]. For this reason, objective biomarkers identification continues to be important to increase accuracy of dietary intake assessment. Some biological markers have been tentatively assessed, with reports of urinary biomarkers of total protein [4], coffee and tea [5] or garlic intake [6], as well as plasmatic biomarkers of fruit and vegetable intake [7], among others, having been used.

Saliva study gained interest due to the potential of this fluid as a non-invasive and stress-free source of biomarkers for different pathologies and physiological conditions [8, 9]. In the nutrition field of study, the relationship between saliva and diet starts now to be studied (e.g., [10]) but many aspects remain to be understood. There are evidences that saliva participates in both oral processing and food sensory perception (e.g., [1113]). Moreover, salivary gland secretion responds to diet, representing a potential source for intake biomarkers that needs to be explored. In animals, it was observed that species having different feeding niches differ in their salivary proteomes [14, 15]. Tannins (a type of polyphenols) are the major dietary compounds linked to the salivary differences that these animal species present, and it is possible that these differences are due to the tannin amounts in food. In fact, changes in salivary glands and in the amounts of particular salivary proteins (e.g., proline-rich proteins and alpha amylase) have been linked to the proportion of polyphenols/tannins in diet (e.g., [1618]). From our knowledge, for humans, limited information exists about potential variations in saliva composition according to polyphenol intake levels, although the perception of astringency induced by these compounds has been demonstrated to be related to saliva composition [19]. Another nutritional compound potentially linked to saliva composition is starch. Variations in starch levels in diet was proposed to be related to differences in amylase gene (Amy1) copy number variants, with higher number of copies in populations with diets rich in starch [20, 21]. Additionally, short-term changes in salivary proteome were observed for bread ingestion [22, 23].

According to what was stated above, it is possible to hypothesize that saliva protein composition can vary among adults with different dietary patterns. Moreover, a relationship between saliva and diet composition was recently proposed, in studies performed in children [10, 24], reinforcing this hypothesis.

The objective of this work was to evaluate the potential of saliva as a non-invasive and objective tool for assessing potential association with dietary intake in adults, which may be of valuable use, particularly in new epidemiological studies.

2. Material and Methods

2.1. Participants

A convenience sample of 122 participants (61 from each sex), 20 to 59 years old, were recruited from North Alentejo region of Portugal, which is an interior region from the country characterized by a close link to agriculture. Participants were randomly recruited from all the municipalities of the referred region, since some variability in food habits may exist between cities and villages. This allowed us to have higher representativity of the region, although comparisons between dietary habits between different regions were not the aim of the present study. Only participants without self-reported and visible signs of oral (e.g., caries) or nasal health problems (e.g., smell loss, obstruction, etc.) participated in this study. On the data collection day, participants were asked to arrive at test room between 10:00 am. and 11:00 am., at least 1 h 30 m after breakfast intake. Before the beginning of the study, all subjects read and signed an informed consent form. All procedures were performed according to the Declaration of Helsinki for Medical Research Involving Human Subjects and had ethical approval from the Ethical Committee of the University of Evora.

2.2. Anthropometric Data Collection

Due to the potential association between Body Mass Index (BMI) and saliva composition [13, 25], height and weight were assessed and measured according to the European Health Examination Survey procedures [26]. Since the relationship between food intake and obesity was not a major aim of this study, only BMI was assessed (rather than other anthropometric measures) for controlling a possible effect in saliva composition. The participants were in a stand position wearing light cloths and barefoot. A portable stadiometer (Seca 214) and a digital scale accurate to the nearest 0.1 kg (Seca 803) were used to access height and weight, respectively. BMI was calculated by dividing the weight (kg) by the square of the height (m2). Normal weight, pre-obese, and obese were considered for the respective values: 18.5 < BMI < 25, 25 < BMI < 30, and BMI > 30 [27]. Since only one individual presented a BMI lower than 18.5, she was considered in normal weight group.

2.3. Food Intake Assessment

Participants were asked to report their frequency of food consumption by completing a self-administered, semi-quantitative food frequency questionnaire (FFQ), validated for Portuguese adults [28]. Before questionnaire completion, the researcher responsible for data collection explained how the questionnaires should be filled. The FFQ is an 86-item questionnaire, that includes individual food items and/or food groups and beverages and assesses usual dietary intake frequency over the previous 12 months. A blank space to introduce other foods, not correspondent to any of the 86 items, was provided. In that blank space, participants wrote the food name, portion, and intake frequency, using the same frequency possibilities. Food intake was calculated by multiply one of the nine possibilities of frequency of consumption (from “never or less than once per month” to “six or more times a day”), by the weight of the standard portion size of the food item, resulting in the daily intake of each food item. Participants were instructed to fill the questionnaire taking into consideration the food portion indicated. For example, if meat is consumed one time per day, in an amount representing the double of the portion indicated, then participant filled as consuming 2 times per day. Energy and nutritional intake were estimated using an adapted Portuguese version of the nutritional analysis software Food Processor Plus (ESHA Research Inc., Salem, OR, USA). When other foods or culinary food preparations, not contained in the 86 items of the FFQ, were reported, their nutritional value was estimated and added to capture the person’s total nutrient intake.

For further constitution of dietary patterns, foods with similar nutritional characteristics were grouped, according to a previous study [29], resulting in 29 food/food groups that were further used in a multivariate analysis, as further described. The daily intake of each food group was obtained by summing the individual daily intakes of the foods that constituted it (for example, the daily intake of fresh fruit represents the sum of individual fruits daily intakes). Since that previous work was in children [29] and did not consider alcoholic beverages in that age group, we added wine and alcoholic beverages different from wine as two variables to include in the constitution of the dietary patterns.

2.4. Saliva Collection and Cleaning

Saliva collection occurred with a minimum of 1 h 30 m from the last food or beverage intake (breakfast), to avoid any influence of previous ingestion in salivary protein composition. Before saliva collection, individuals were asked to rinse with water to remove any food debris and “old” saliva. Saliva was collected without stimulation. Participants were requested to accumulate all saliva produced in the mouth and spitting it to a tube, maintained on ice, each time they need it. This collection occurred for a period of 3 minutes. Saliva was maintained on ice until laboratory arrival, where it was stored for one day at −20°C. To remove mucins and cell and/or food residues, saliva samples were thawed on ice and centrifuged at 13,000  for 30 minutes at 4°C. The supernatant was recovered and stored at −80°C until subsequent analysis.

2.5. Saliva Composition Analysis
2.5.1. Saliva Flow Rate, pH, and Total Protein Concentration

Saliva flow rate was assessed by assuming that saliva density is 1.0. Tubes containing saliva were weighed, with the weight of the empty tube being subtracted. The final value was divided by the number of minutes during which saliva was collected. The pH of saliva samples was measured using a calibrated pH meter (Hanna Instruments) and recording to two decimal places. Total protein concentration was determined by the Bradford method, using bovine serum albumin (BSA) as standard, and plates were read at 600 nm in a microplate reader (Glomax, Promega).

2.5.2. SDS-PAGE Separation and Protein Profile Analysis

Each saliva sample was run in duplicate. For each sample, a volume corresponding to 7.5 μg total protein was mixed with sample buffer and run on each lane of a 14% polyacrylamide mini-gel (Protean xi, Bio-Rad, CA, USA) using a Laemmli buffer system, as described elsewhere [30]. Each electrophoretic run was performed at a constant voltage of 140 V until front dye reached the end of the gel. Gels were fixed for 1 hour in 40% methanol/10% acetic acid, followed by staining for 2 hours with Coomassie Brilliant Blue (CBB) G-250. Gel images were acquired using a scanning Molecular Dynamics densitometer with internal calibration and LabScan software (GE Healthcare), and images were analysed using GelAnalyzer software (GelAnalyzer 2010a by Istvan Lazar, for the volume percentage of each protein band. Molecular masses were determined in accordance with molecular mass standards (Bio-Rad Precision Plus Protein Dual Colour 161–0394) run with protein samples. The identification of the proteins contained in the bands observed in SDS-PAGE salivary profiles was based on previous data [22].

2.5.3. Salivary Amylase Enzymatic Activity

A Salimetrics® kit was used to determine the enzymatic activity of salivary amylase according to the manufacturer’s recommendations. Briefly, saliva samples were diluted 200× and applied on the microplate in duplicate, followed by application of a substrate (2-chloro-p-nitrophenol) preheated to 37°C. The mixture was incubated at 37°C for 1 minute, and absorbance values were read at 405 nm in a plate reader spectrophotometer, followed by incubation for an additional 2 minutes at 37°C and a new reading at 405 nm. The enzymatic activity of amylase (U/ml) was calculated by the following formula: (ΔAbs./min × TV × DF)/(MMA × SV × LP), where ΔAbs./Min is absorbance variation per minute, TV is total test volume (0.287 mL), DF is dilution factor, MMA is millimolar absorbance of substrate 2-chloro-p-nitrophenol (12.9), SV is sample volume (0.007 mL), and LP is light path (0.97, specific for plate received with kit).

2.6. Statistical Analysis

The values ​​of total protein concentration, salivary secretion rate, protein band amount (volume percentage), and salivary amylase enzymatic activity were analysed statistically. Descriptive statistics was performed, and normality and homoscedasticity were tested through Shapiro–Wilk and Levene tests, respectively. Continuous variables were compared between sexes using Student’s t-test when assumptions were fulfilled and using Mann–Whitney when not. For BMI classes (normal weigh, pre-obese and obese) comparison, one-way ANOVA or the non-parametric Kruskal–Wallis test was used, considering a confidence interval of 95%.

To identify dietary patterns within the study population, multivariate statistical techniques were used. The 29 foods/food groups constituted from the 86 items of the FFQ, as described earlier and according to a previous study [29], were considered together with 2 additional groups, namely, wine, and alcoholic beverages excluding wine. These 31 variables (amount of food consumed per day for each food/food group) were reduced through principal component analysis (PCA). PCA with orthogonal rotation (Varimax) was used to estimate the latent factors emerging and to obtain optimal non-correlated components. The number of components (dietary patterns) to be retained was decided based on the observation of correlation matrix, on the total of variance explained by the components and on Kaiser–Meyer–Olkin (KMO) measure of sampling adequacy values above 0.600 [31, 32].

With the aim of explaining α-amylase enzymatic activity and salivary protein profile variation (namely, the expression levels of each SDS-PAGE protein band) as dependent variables, each dietary pattern extracted from the PCA analysis was used as independent variable in the regression model, adjusting for age, BMI, sex, and total energy consumed. A stepwise backward process was used to obtain the final model after looking for changes in adjusted R2 and F values for each retained variable; the assumptions of collinearity (VIF and tolerance), non-dependent errors (Durbin–Watson), and homoscedasticity (residual analysis) were also considered.

3. Results

3.1. Participants Characteristics and Dietary Patterns

Forty percent of the individuals participating in this study were normal weight and 60% overweight (36% pre-obese + 24% obese). When comparing the sexes, men presented statistically significantly higher BMI than women (Table 1), with 67.7% of men being overweight (38.7% pre-obese + 29.0% obese) whereas 29.8% of women being pre-obese and 17.5% obese (e.g., a total of 47.3% overweight women).

ParametersMenWomen-value (sex)NWPre-obeseObese-value (BMI)

Age (years)40.7 ± 10.938.3 ± 10.70.47737.1 ± 10.941.5 ± 1.742.6 ± 9.6.054
BMI (kg/m2)27.9 ± 4.225.9 ± 5.20.014

Nutritional parameters
Energy (kcal/day)2346 ± 7462280 ± 7340.6272304 ± 7232361 ± 7872211 ± 5380.650
Protein (g/day)102.8 ± 37.6107.8 ± 35.60.457103.1 ± 36.1105.5 ± 35.7106.5 ± 32.10.924
Protein (% TEV)17.6 ± 3.919.2 ± 3.70.02718.1 ± 4.018.1 ± 3.319.4 ± 4.30.198
Total carbohydrate (g/day)263.4 ± 101.7256.1 ± 92.20.911265.3 ± 95.0268.4 ± 106.5236.8 ± 75.00.472
Total carbohydrate (% TEV)44.9 ± 10.044.9 ± 8.80.97646.3 ± 11.045.4 ± 7.042.4 ± 8.10.143
Total sugars (g/day)119.0 ± 58.5113.5 ± 51.40.853124.9 ± 60.5111.5 ± 47.7103.6 ± 44.60.351
Total sugars (% TEV)20.5 ± 7.819.9 ± 6.40.65121.6 ± 8.619.4 ± 4.918.5 ± 5.90.226
SFA (g/day)26.6 ± 11.026.3 ± 9.50.95126.7 ± 10.126.8 ± 11.425.6 ± 8.90.866
SFA (% TEV)10.3 ± 2.610.5 ± 2.30.68410.5 ± 2.510.2 ± 2.410.5 ± 2.20.949
MUFA (g/day)41.2 ± 16.942.9 ± 17.80.53540.2 ± 16.343.0 ± 17.941.4 ± 13.90.685
MUFA (% TEV)15.8 ± 4.017.1 ± 4.00.07216.0 ± 4.016.3 ± 3.916.8 ± 3.30.458
PUFA (g/day)15.8 ± 6.015.8 ± 6.40.90114.6 ± 5.716.6 ± 6.215.9 ± 5.80.137
PUFA (% TEV)6.1 ± 1.56.3 ± 1.40.3515.9 ± 1.56.3 ± 1.36.4 ± 1.30.106
Cholesterol (mg/day)361.3 ± 167.7345.6 ± 134.00.578345.0 ± 148.3358.6 ± 165.9369.3 ± 141.70.774
Fibre (g/day)26.2 ± 13.028.4 ± 12.90.20627.8 ± 13.028.5 ± 13.425.7 ± 12.50.695

Foods/food groups (g/day)
Vegetables387.9 ± 268.0499.7 ± 255.40.038434.9 ± 283.1424.4 ± 230.6543.4 ± 265.50.106
Vegetable soup194.5 ± 167.4272.0 ± 154.60.003247.6 ± 237.2206.0 ± 130.1298.0 ± 136.90.016
Fresh fruits253.8 ± 177.1353.6 ± 239.30.027348.6 ± 270.3304.6 ± 196.8255.1 ± 188.30.281
Canned fruits10.7 ± 21.25.6 ± 12.80.0066.4 ± 13.59.3 ± 23.57.8 ± 11.80.300
Olives8.1 ± 15.413.4 ± 36.60.23014.3 ± 38.98.9 ± 18.27.7 ± 12.30.719
Nuts10.6 ± 15.017.2 ± 22.90.19313.9 ± 29.915.5 ± 19.317.2 ± 23.20.223
Fish72.7 ± 46.382.9 ± 6.40.16868.5 ± 43.682.7 ± 48.289.6 ± 55.60.156
White meat55.1 ± 47.072.9 ± 47.70.039
Red meat43.0 ± 34.330.8 ± 27.30.03232.7 ± 27.838.1 ± 35.438.8 ± 27.90.572
Processed meat131.9 ± 40.517.5 ± 14.30.01223.8 ± 34.424.5 ± 35.024.6 ± 17.10.347
Eggs14.1 ± 11.114.7 ± 9.80.39712.9 ± 9.515.4 ± 11.615.7 ± 10.50.325
Olive oil10.9 ± 9.512.9 ± 9.10.08411.4 ± 10.614.0 ± 12.911.6 ± 9.10.868
Vegetable oil2.1 ± 3.51.5 ± 3.10.2131.0 ± 2.42.3 ± 3.32.0 ± 3.80.006
Butter2.5 ± 2.23.5 ± 3.00.0592.9 ± 2.32.9 ± 2.72.9 ± 3.10.843
Margarine1.2 ± 1.61.3 ± 1.60.4381.7 ± 2.21.5 ± 2.31.2 ± 2.40.366
Milk167.7 ± 164.4231.5 ± 206.70.114199.4 ± 212.4216.4 ± 179.4186.9 ± 152.50.524
Yoghurt62.7 ± 85.393.6 ± 90.60.00984.6 ± 99.960.4 ± 72.686.0 ± 90.90.399
Cheese18.1 ± 21.315.7 ± 16.30.88916.7 ± 19.316.5 ± 19.219.0 ± 19.20.623
Milk-based puddings13.9 ± 16.69.2 ± 9.20.32311.1 ± 13.212.0 ± 15.310.2 ± 10.70.982
Ice-cream7.5 ± 8.75.6 ± 8.10.0836.3 ± 8.47.9 ± 10.44.9 ± 4.70.791
Starch-rich foods (rice, potatoes, pasta)2136.8 ± 72.0136.4 ± 71.11.000120.1 ± 66.7156.5 ± 53.8153.8 ± 109.40.021
Bread63.8 ± 51.777.0 ± 53.70.11076.6 ± 54.775.1 ± 57.159.4 ± 39.70.542
Ready-to-eat cereals11.6 ± 14.89.3 ± 14.10.3069.5 ± 14.212.5 ± 15.37.5 ± 12.90.357
Sugar sweetened beverages (SSB)3170.1 ± 216.384.6 ± 100.20.100134.6 ± 192.4107.2 ± 145.3128.1 ± 180.60.839
Coffee76.5 ± 57.663.4 ± 49.50.33064.5 ± 56.674.5 ± 49.578.8 ± 50.70.585
Tea32.4 ± 58.668.3 ± 88.30.24642.2 ± 79.745.2 ± 74.043.1 ± 72.10.994
Fast-food437.8 ± 27.829.3 ± 15.40.16332.7 ± 23.832.1 ± 20.337.5 ± 33.80.936
Pastry40.1 ± 36.430.6 ± 25.70.24536.0 ± 31.040.3 ± 38.125.7 ± 22.70.237
Crackers/cookies55.5 ± 7.610.9 ± 13.50.0299.8 ± 12.47.7 ± 10.46.8 ± 9.90.433
Pulses51.9 ± 57.635.5 ± 36.70.49844.2 ± 49.946.6 ± 50.142.9 ± 45.60.975
Wine36.6 ± 41.312.2 ± 25.3<0.00117.4 ± 30.825.4 ± 36.938.0 ± 42.60.072
Alcoholic beverages (apart from wine)661.9 ± 74.514.1 ± 30.4<0.00126.1 ± 43.643.1 ± 63.252.1 ± 83.60.173

The food/food groups according to previous work [29], which were further used to constitute the dietary patterns.1Processed meats include ham, sausages, and bacon; 2Starch-rich foods include pasta, rice, and potato. 3Sugar sweetened beverages include all types of juices with sugar, cola, and ice-tea. 4Fast-food includes hamburger, pizza, snacked fried food, and sausages. 5Cookies with sugar amounts lower that 20% of total ingredients. 6Beer + spirits.  < .05. SFA, saturated fatty acids; MUFA, monounsaturated fatty acids; PUFA, polyunsaturated fatty acids; TEV, total energy value.

When participants were analysed for their food intake, assessed by the FFQ, some statistically significant differences were observed between men and women in the energy contribution of protein (higher in women) and in the consumption of some food groups (Table 1). Women have significant higher consumption of vegetables, vegetable soup, yoghurts, and crackers/cookies with approximately 20% sugar or less, as well as white meat and lower consumption levels of canned fruit, red meat, processed meat, and alcoholic drinks. No significant differences between sexes were found in relation to mean total energy or total macronutrient intake.

It was observed a tendency for older individuals to have higher BMI. Despite no significant differences in total energy or total macronutrient intake, for foods like vegetable soup, vegetable oil, and starchy-rich foods, the daily intake was observed to be higher in obese individuals (Table 1).

To estimate the dietary patterns from the food/food groups considered [29], PCA was performed as described in material and methods section. PCA adequacy was evaluated before analysis. Matrix component inspection allowed choosing the variables (food/food groups) with coefficients higher than 0.30 in at least one component, resulting in 24 variables (food/food groups) including the final model. The general measure of KMO was 0.606, and Bartlett sphericity was statistically significant ( < 0.0001).

It was possible to obtain 7 components, which explained 59% of total variance (Table 2), corresponding to the following dietary patterns: DP1 with positive loadings for vegetables and soup, olives, fresh fruit, bread and starch-rich foods, and sugar sweetened beverages (SSB) (plant-based with SSB DP); DP2 presented positive loadings for fish, meat, eggs, and starch-rich foods (animal-based and starchy foods pattern); DP3 presented high positive scores for fast food, processed foods, and SSB (fast-food pattern); DP4 presented positive loadings for olive oil, margarine, and butter, as well as for pastry and negative loadings for pulses (dietary fats and pastry pattern); DP5 presented positive loadings for dairy desserts, pastry, canned fruits, and ready-to-eat cereals (sweet foods pattern); DP6 presented positive loadings for olives and nuts, but negative loadings for wine (plant-based fatty foods without wine pattern); and DP7 presented positive loadings for yoghurt, cheese, and wine (cheese/yoghurt and wine pattern).

% Var explainedComponent

Milk-based puddings0.599
Processed meat0.575
Olive oil0.665
Ready-to-eat cereals0.584
Starch-rich foods (rice, potatoes, pasta)0.4560.471
Fresh fruits0.521
Canned fruits0.727
Vegetable soup0.712

Rotation converted with 16 interactions. Values lower than 0.40 were omitted. 1Sugar sweetened beverages.
3.2. Saliva Composition

SDS-PAGE electrophoresis allowed the consistent separation and comparison of well resolved 12 protein bands, with apparent molecular masses between 14.1 and 88.0 kDa (Figure 1), containing proteins previously identified [22] (Table 3).

SDS-PAGE BandProtein identification#Apparent molecular mass (kDa)

APolymeric immunoglobulin receptor88.0
BSerum albumin71.0
Cα-Amylase 164.0
Dα-Amylase 158.2
Eα-Amylase 155.0
Carbonic anhydrase VI
GZymogen granule protein 16 homolog B31.5
Immunoglobulin kappa constant
HImmunoglobulin kappa constant27.0
I1Prolactin inducible protein (PIP)16.8
I2Prolactin inducible protein (PIP)16.2
KCystatin B14.1

Protein identification based on previous studies with human saliva collected in the same conditions.

By comparing the sexes for saliva composition, men presented higher expression levels of the protein bands G (4.39 ± 0.73 vs. 3.98 ± 0.85, for men and women, respectively;  = .005) and J (8.54 ± 2.74 vs. 7.41 ± 2.14, for men and women, respectively;  = .014).

3.3. Dietary Patterns and Saliva Composition

The significant associations found for linear multiple regression models between salivary parameters, namely, amylase enzymatic activity, protein bands E (identified as amylase) and J (identified as cystatins), and dietary patterns, are shown in Table 4. The protein bands A, B, I1, I2, and K did not show significant associations, and protein bands C, D, and G were only associated with age (positively), BMI (positively), and sex (high in men), respectively.

Dependent variableIndependent variableCoeffCI (95%)t valueFR2 adjDurbin− Watson

Band EConstant1.3900.686 to 2.0953.912<0.0012.7650.0721.953
Dietary pattern 70.1060.012 to 0.1992.2470.027
Energy6.208E−60.000 to 0.0000.1360.892
Sex0.023−0.159 to 0.2040.2480.804
Age−0.008−0.017 to 0.000−1.8680.064
BMI0.0270.008 to 0.0472.7660.007

Band JConstant9.0355.472 to 12.5985.027<0.0013.8680.1321.598
Dietary pattern 7.6910.224 to 1.1572.9340.004
Energy0.000−0.001 to 0.000−0.5130.609
Sex−0.625−1.574 to .323−1.3070.194
Age0.041−0.002 to 0.0851.8920.061
BMI−0.054−0.151 to 0.043−1.0990.274

α-Amylase (U/mL)Constant167.348−46.456 to 381.1521.5530.1244.1260.1501.569
Dietary pattern 242.66813.620 to 71.7152.9140.004
Dietary pattern 631.0122.395 to 59.6292.1500.034
Energy−0.007−0.037 to 0.023−.4590.647
Sex18.932−36.866 to 74.7300.6730.502
Age4.4741.905 to 7.0433.4560.001
BMI−5.641−11.433 to 0.150−1.9330.056

All the models were adjusted for sex, age, BMI, and energy intake. Significant for  < 0.05.

Looking at the linear multiple regression models between dietary patterns and salivary parameters, adjusting for confounders (sex, age, BMI, and energy intake), in Table 4, it was possible to observe that DP7 (fermented dairy and wine pattern) was positively associated with bands E (containing amylase) and J (containing S-type cystatins), and DP6 (plant-based fatty foods without wine pattern) and DP2 (animal-based and starchy foods pattern) were positively associated with amylase enzymatic activity.

4. Discussion

In the present study, PCA analysis was used to identify combinations of foods, named as dietary patterns and which constituted not-correlated variables, used as independent variables in the regression models. As such, these dietary patterns do not represent groups of individuals with particular food consumption patterns.

Dietary patterns constituted by high proportion of plant-based fatty foods (without wine, DP6) together with dietary patterns based on foods of animal origin and starch-rich foods (DP2) were associated with the salivary amylase enzymatic activity. This is in line with other studies associating salivary amylase with starchy foods intake [21] although, to the best of our knowledge, no previous studies related salivary amylase with olives and nuts intake. The different results obtained when amylase was considered in terms expression levels, comparatively to when it was considered in terms of its enzymatic activity, are not surprising, since a lack of association between the enzymatic activity of this protein and the expression levels of its different proteoforms has been previously reported [33].

In relation to a dietary pattern with high contribution of fermented dairy and wine (DP7), an association with the expression level of one band containing α-amylase and one band containing S-type cystatins was observed. Concerning wine, it may be a source of tannins, which are known to affect saliva proteome, in animal models [16, 34]. In fact, in mice, it was previously observed that increased levels of polyphenols (tannins) in the diet result in increased expression levels of salivary amylase bands, in SDS-PAGE profiles [18]. Concerning cystatins, these are inhibitors of cysteine proteases, the family of S-type cystatins being secreted by salivary glands. Different studies show a relationship between salivary cystatins and the bitterness or astringency (e.g., [13, 35, 36]) that is associated with the polyphenol content of foods and beverages, such as wines (particularly red wines). If higher levels of salivary S-type cystatins are present in individuals with low sensitivity to bitterness and astringency, as reported by some authors [11], it can be hypothesized that the positive association between S-type cystatins and DP7 results from the lower intensity with which these oral sensations are perceived and consequently the higher acceptance/preference for wine, in consequence of salivary cystatins. At the same time, another hypothesis is that a higher intake of wine can induce higher levels of salivary S-type cystatins. This hypothesis is supported by studies showing that exposing rats to tannins or bitter compounds results in the increase of S-type salivary cystatins secretion and in the increase of bitter taste compounds acceptance [35, 37]. Anyway, this hypothesis needs to be further tested. It is worthwhile to remember that DP7 has also a considerable contribution of fermented dairy. Although few information exists about effect of dairy in saliva composition, the levels of salivary cystatins were recently reported to increase in response to chocolate milk [38].

The association between the dietary patterns obtained by PCA analysis and saliva composition may be additionally confounded by sex, BMI, and age. As such, those parameters were included in the regression models used. Although we are unable to guarantee that the observed salivary differences between sexes are due exclusively to the sex factor, and not to the differences in the eating habits of women and men, the effect of sex in salivary proteome needs to be considered, since it has been already reported [13, 39, 40] and reinforced in recent publications (e.g., [41]). Minor differences were observed among BMI groups in food consumption habits, although starchy foods were highly consumed by overweight and obese individuals. As for the case of sex, we cannot guarantee that the effect of BMI in saliva was not influenced by this particular difference in food intake. As such, the positive association between BMI and the band containing one amylase form (band E) can result from this higher intake of starchy foods by individuals with high BMI. But, since band E was not related to the patterns with starch-rich foods (DP2), other factors apart from starch consumption may explain this positive association between band E and obesity. Differences in saliva composition among individuals with different BMI have been already observed, which included variations in the expression levels of salivary proteins like carbonic anhydrase VI and some forms of amylase [30]. Finally, the positive association between salivary amylase enzymatic activity and age goes in line with previous observations about increasing salivary amylolytic activity in people with higher age [42].

4.1. Our Study Has Some Limitations

(1)Firstly, although in the regression model for statistical analysis, α-amylase enzymatic activity and salivary protein profile were considered as dependent variables, while dietary patterns extracted from the PCA analysis were used as independent variables, the cross-sectional design of our study does not allow us to establish causal relationships between dietary intake and salivary parameters. Moreover, since the dietary patterns considered are, in fact, the components generated by PCA, which consist in the combinations of the different food/food groups, it is not possible to allocate each individual subject to one of these dietary patterns.(2)Secondly, whereas the salivary parameters represent the saliva collected at one time point, the intake reported represents the frequency of the food groups consumption in the last 12 months. Short-term effect of food intake was not controlled. Moreover, it is recognized that FFQ have the limitation of being difficult to carry out to complete, and suffer from difficulties in self-estimation of portion size and biases resulting from misreporting [43]. However, the FFQ used is validated for Portuguese adults and considered as a valid tool to estimate the usual food intake in large samples. Intervention-based studies must be further done to get more controlled results about the effect of specific dietary patterns in saliva composition.

Our study has also important strengths. According to our best knowledge, this is the first study that aimed to assess how dietary habits relate to saliva protein composition. One of the main findings is that some types of salivary proteins appear to be associated with the intake of particular foods, namely, animal-based and starch-rich foods, plant-based fatty foods, fermented dairy, and wine. It is a fact that the regression models obtained presented relatively low R2, which means that the model only explains a minor percentage of the variation in salivary proteins levels. But taking into consideration that we are working with biological models, where several other factors besides food frequency consumption can affect salivary protein composition [44, 45], this low R2 is not surprising and the statistical significance of the model reinforces an association between salivary proteins and food habits.

5. Conclusions

The results obtained in the present study emphasize the potential of saliva to reflect differences in food intake. Salivary amylase and S-type cystatins are associated with dietary patterns. It is interesting to note that some combinations of foods are more related to the total enzymatic activity of amylase whereas others with the expression levels of forms of this protein. Salivary amylase and cystatins are proteins previously observed to have a link to oral food perception and more studies are necessary to confirm if they can be the result of the type of food eaten or rather the cause of food choices.

To consider salivary proteins as non-invasive biomarkers of intake, more studies are required, to define criteria, such as sex, BMI, age, and hour of the day, among others. Moreover, it will be important, in further work, to get more detailed information using other proteomic approaches. Even so, this study evidences that some of the variation in salivary protein composition, among individuals, can be related to food habits and it is important to take into account foods as a source of variation, when using saliva for studies with different purposes.

Data Availability

All the data necessary to reach the objectives of the study are included. Raw data files, which generated the results, are with the researchers.


The FCT-Portuguese Science Foundation was not involved in carrying out this study or submitting it for publication.

Conflicts of Interest

The authors have no conflicts of interest.


The authors would like to thank all the subjects who took part in the study. This paper was funded by the FEDER as part of the Operational Programme for Competitiveness Factors (COMPETE) and received funding from the Science and Technology Foundation (FCT) under the Project UID/AGR/00115/2019 (ICAAM, University of Évora). Funding was additionally provided by the FCT-Portuguese Science Foundation, Research Contract CEECIND/04397/2017 to Elsa Lamy. Support was also received from the State of São Paulo Research Foundation (FAPESP), grant no. 2017/26400-6.


  1. J.-S. Shim, K. Oh, and H. C. Kim, “Dietary assessment methods in epidemiologic studies,” Epidemiology and Health, vol. 22, no. 36, Article ID e2014009, 2014. View at: Publisher Site | Google Scholar
  2. J. Maurer, D. L. Taren, P. J. Teixeira et al., “The psychosocial and behavioral characteristics related to energy misreporting,” Nutrition Reviews, vol. 64, no. 2, 2006. View at: Publisher Site | Google Scholar
  3. R. K. Johnson, R. P. Soultanakis, and D. E. Matthews, “Literacy and body fatness are associated with underreporting of energy intake in US low-income women using the multiple-pass 24-hour recall: a doubly labeled water study,” Journal of the American Dietetic Association, vol. 98, no. 10, pp. 1136–1140, 1998. View at: Publisher Site | Google Scholar
  4. S. A. Bingham, “Urine nitrogen as a biomarker for the validation of dietary protein intake,” Journal of Nutrition, vol. 24, pp. 921S–924S, 2003. View at: Publisher Site | Google Scholar
  5. J. A. Rothwell, F. Madrid-Gambin, M. Garcia-Aloy et al., “Biomarkers of intake for coffee, tea, and sweetened beverages,” Genes & Nutrition, vol. 15, 2018. View at: Publisher Site | Google Scholar
  6. L. Scheffler, Y. Sauermann, A. Heinlein, C. Sharapa, and A. Buettner, “Detection of volatile metabolites derived from garlic (Allium sativum) in human urine,” Metabolites, vol. 6, no. 4, p. 43, 2016. View at: Publisher Site | Google Scholar
  7. M. Pennant, M. Steur, C. Moore, A. Butterworth, and L. Johnson, “Comparative validity of Vitamin C and carotenoids as indicators of fruit and vegetable intake: a systematic review and meta-analysis of randomised controlled trials,” British Journal of Nutrition, vol. 231, pp. 14–18, 2013. View at: Publisher Site | Google Scholar
  8. E. Lamy, “Salivary proteomics in ingestive behaviour research: advances, potentialities and limitations,” Journal of Integrated OMICS, vol. 8, no. 1, 2018. View at: Publisher Site | Google Scholar
  9. M. Castagnola, E. Scarano, G. C. Passali et al., “Salivary biomarkers and proteomics: future diagnostic and clinical utilities,” ACTA Otorhinolaryngologica Italica, vol. 37, no. 2, pp. 94–101, 2017. View at: Publisher Site | Google Scholar
  10. M. Morzel, C. Truntzer, E. Neyraud et al., “Associations between food consumption patterns and saliva composition: specificities of eating difficulties children,” Physiology & Behavior, vol. 173, pp. 116–123, 2017. View at: Publisher Site | Google Scholar
  11. M. Dsamou, O. Palicki, C. Septier et al., “Salivary protein profiles and sensitivity to the bitter taste of caffeine,” Chemical Senses, vol. 37, pp. 87–95, 2012. View at: Publisher Site | Google Scholar
  12. T. Cabras, M. Melis, M. Castagnola et al., “Responsiveness to 6-n-propylthiouracil (PROP) is associated with salivary levels of two specific basic proline-rich proteins in humans,” PLoS One, vol. 7, Article ID e30962, 2012. View at: Publisher Site | Google Scholar
  13. L. Rodrigues, G. da Costa, C. Cordeiro, C. C. Pinheiro, F. Amado, and E. Lamy, “Relationship between saliva protein composition and 6-n-Propylthiouracil bitter taste responsiveness in young adults,” Journal of Sensory Studies, vol. 32, 2017. View at: Publisher Site | Google Scholar
  14. E. Lamy, G. da Costa, R. Santos et al., “Sheep and goat saliva proteome analysis: a useful tool for ingestive behavior research?” Physiology & Behavior, vol. 98, pp. 393–401, 2009. View at: Publisher Site | Google Scholar
  15. E. Sales Baptista, E. Lamy, M. Mau, F. Capela e Silva, and A. V. Coelho, “Variation in salivary protein composition related to feeding behavior and its ecological implications,” 2013. View at: Google Scholar
  16. E. Lamy, G. Graça, G. da Costa et al., “Changes in mouse whole saliva soluble proteome induced by tannin-enriched diet,” Proteome Science, vol. 8, p. 65, 2010. View at: Publisher Site | Google Scholar
  17. E. Lamy, E. S. Baptista, A. V. Coelho, and F. C. Silva, “Morphological alterations in salivary glands of mice (Mus musculus) submitted to tannin enriched diets: comparison with sialotrophic effects of sympathetic agonists stimulation,” The Arquivo Brasileiro de Medicina Veterinária e Zootecnia, vol. 62, 2010. View at: Publisher Site | Google Scholar
  18. G. da Costa, E. Lamy, F. Capela e Silva, J. Andersen, E. Sales Baptista, and A. V Coelho, “Salivary amylase induction by tannin-enriched diets as a possible countermeasure against tannins,” Journal of Chemical Ecology, vol. 34, pp. 376–387, 2008. View at: Publisher Site | Google Scholar
  19. C. Dinnella, A. Recchia, S. Vincenzi, H. Tuorila, and E. Monteleone, “Temporary modification of salivary protein profile and individual responses to repeated phenolic astringent stimuli,” Chemical Senses, vol. 35, pp. 75–85, 2010. View at: Publisher Site | Google Scholar
  20. A. L. Mandel, C. Peyrot des Gachons, K. L. Plank, S. Alarcon, and P. A. S. Breslin, “Individual differences in AMY1 gene copy number, salivary α-amylase levels, and the perception of oral starch,” PLoS One, vol. 5, Article ID e13352, 2010. View at: Publisher Site | Google Scholar
  21. J. L. Santos, E. Saus, S. V Smalley et al., “Copy number polymorphism of the salivary amylase gene: implications in human nutrition research,” Journal of Nutrigenetics and Nutrigenomics, vol. 5, pp. 117–131, 2012. View at: Publisher Site | Google Scholar
  22. L. Carreira, P. Midori Castelo, C. Simões, F. Capela e Silva, C. Viegas, and E. Lamy, “Changes in salivary proteome in response to bread odour,” Nutrients, vol. 12, p. 1002, 2020. View at: Publisher Site | Google Scholar
  23. E. Lamy, V. Santos, S. Barrambana et al., Saliva Protein Composition Relates with Interindividual Variations in Bread Sensory Ratings, Starch - Stärke, Wiley, Hoboken, NJ, USA, 2020.
  24. A. Tvarijonaviciute, N. Martinez-Lozano, R. Rios, M. C. Marcilla de Teruel, M. Garaulet, and J. J. Cerón, “Saliva as a non-invasive tool for assessment of metabolic and inflammatory biomarkers in children,” Clinical Nutrition, vol. 39, no. 8, pp. 2471–2478, 2019. View at: Publisher Site | Google Scholar
  25. L. Rodrigues, R. Espanca, A. R. Costa et al., “Comparison of salivary proteome of children with different sensitivities for bitter and sweet tastes: association with body mass index,” International Journal of Obesity, vol. 43, no. 4, pp. 701–712, 2019. View at: Publisher Site | Google Scholar
  26. B. Nunes, M. Barreto, A. P. Gil et al., “The first Portuguese National Health Examination Survey (2015): design, planning and implementation,” Journal of Public Health, vol. 41, 2019. View at: Publisher Site | Google Scholar
  27. WHO, “Preventing and Managing the Global Epidemic. World Health Organization: Technical Report Series,” WHO, Geneva, Switzerland, 2000. View at: Google Scholar
  28. C. Lopes, A. Aro, A. Azevedo, E. Ramos, and H. Barros, “Intake and adipose tissue composition of fatty acids and risk of myocardial infarction in a male Portuguese community sample,” Journal of the American Dietetic Association, vol. 107, pp. 276–286, 2007. View at: Publisher Site | Google Scholar
  29. P. Moreira, S. Santos, P. Padrão et al., “Food patterns according to Sociodemographics, physical activity, sleeping and obesity in Portuguese children,” International Journal of Environmental Research and Public Health, vol. 7, no. 3, pp. 1121–1138, 2010. View at: Publisher Site | Google Scholar
  30. E. Lamy, C. Simões, L. Rodrigues et al., “Changes in the salivary protein profile of morbidly obese women either previously subjected to bariatric surgery or not,” Journal of Physiology and Biochemistry, vol. 71, 2015. View at: Publisher Site | Google Scholar
  31. H. F. Kaiser, “An index of factorial simplicity,” Psychometrika, vol. 39, pp. 31–36, 1974. View at: Publisher Site | Google Scholar
  32. H. F. Kaiser, “The application of electronic computers to factor Analysis,” Educational and Psychological Measurement, vol. 20, pp. 141–151, 1960. View at: Publisher Site | Google Scholar
  33. M. D. Contreras-Aguilar, D. Escribano, S. Martínez-Subiela, S. Martínez-Miró, J. J. Cerón, and F. Tecles, “Changes in alpha-amylase activity, concentration and isoforms in pigs after an experimental acute stress model: an exploratory study,” BMC Veterinary Research, vol. 14, p. 256, 2018. View at: Publisher Site | Google Scholar
  34. E. Lamy, H. Rawel, F. J. Schweigert et al., “The effect of tannins on mediterranean ruminant ingestive behavior: the role of the oral cavity,” Molecules, vol. 16, 2011. View at: Publisher Site | Google Scholar
  35. L. E. Martin, L. V. Nikonova, K. Kay, A. B. Paedae, R. J. Contreras, and A. M. Torregrossa, “Salivary proteins alter taste-guided behaviors and taste nerve signaling in rat,” Physiology & Behavior, vol. 184, pp. 150–161, 2018. View at: Publisher Site | Google Scholar
  36. S. Soares, N. Mateus, and V. de Freitas, “Interaction of different classes of salivary proteins with food tannins,” Food Research International, vol. 49, pp. 807–813, 2012. View at: Publisher Site | Google Scholar
  37. L. E. Martin, L. V. Nikonova, K. E. Kay, and A. M. Torregrossa, “Altering salivary protein profile can increase acceptance of a novel bitter diet,” Appetite, vol. 136, pp. 8–17, 2019. View at: Publisher Site | Google Scholar
  38. C. R. Crawford and C. A. Running, “Addition of chocolate milk to diet corresponds to protein concentration changes in human saliva,” Physiology & Behavior, vol. 225, 2020. View at: Publisher Site | Google Scholar
  39. L. Rodrigues, R. Espanca, A. R. Costa et al., “Association between salivary leptin levels and taste perception in children,” Journal of Nutrition and Metabolism, vol. 2017, 2017. View at: Publisher Site | Google Scholar
  40. L. Rodrigues, G. Costa, C. Cordeiro, C. Pinheiro, F. Amado, and E. Lamy, “Salivary proteome and glucose levels are related with sweet taste sensitivity in young adults,” Food & Nutrition Research, vol. 61, Article ID 1389208, 2017. View at: Publisher Site | Google Scholar
  41. L. Franco-Martínez, J. M. González-Hernández, A. Horvatić et al., “Differences on salivary proteome at rest and in response to an acute exercise in men and women: a pilot study,” Journal of Proteomics, vol. 214, Article ID 103629, 2020. View at: Publisher Site | Google Scholar
  42. A. C. Mosca, M. Stieger, E. Neyraud, H. Brignot, A. van de Wiel, and J. Chen, “How are macronutrient intake, BMI, ethnicity, age, and gender related to the composition of unstimulated saliva? a case study,” Journal of Texture Studies, vol. 50, no. 1, 2019. View at: Publisher Site | Google Scholar
  43. A. Skinner, Z. Toumpakari, C. Stone, and L. Johnson, “Future directions for integrative objective assessment of eating using wearable sensing Technology,” Frontiers in Nutrition, vol. 7, 2020. View at: Publisher Site | Google Scholar
  44. C. V. Esteves, W. G. de Campos, M. M. de Souza, S. V. Lourenço, W. L. Siqueira, and C. A. Lemos-Júnior, “Diagnostic potential of saliva proteome analysis: a review and guide to clinical practice,” Brazilian Oral Research, vol. 33, p. e043, 2019. View at: Publisher Site | Google Scholar
  45. F. M. L. Amado, R. P. Ferreira, and R. Vitorino, “One decade of salivary proteomics: current approaches and outstanding challenges,” Clinical Biochemistry, vol. 46, pp. 506–517, 2013. View at: Publisher Site | Google Scholar

Copyright © 2021 Teresa Louro et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.