animals Article Evaluation of the Effect of Sampling Time on Biomarkers of Stress, Immune System, Redox Status and Other Biochemistry Analytes in Saliva of Finishing Pigs Alba Ortín-Bustillo 1 , María D. Contreras-Aguilar 1 , Camila P. Rubio 2 , María Botia 1 , José J. Cerón 1 , Marina López-Arjona 1 , Silvia Martínez-Subiela 1 , Damián Escribano 1,3 and Fernando Tecles 1, * 1 Interdisciplinary Laboratory of Clinical Analysis of the University of Murcia (Interlab-UMU), Regional Campus of International Excellence ‘Campus Mare Nostrum’, University of Murcia, Campus de Espinardo s/n, 30100 Murcia, Spain 2 Department of Animal and Food Science, Universitat Autònoma de Barcelona, 08193 Barcelona, Spain 3 Department of Animal Production, Regional Campus of International Excellence ‘Campus Mare Nostrum’, University of Murcia, Campus de Espinardo s/n, 30100 Murcia, Spain * Correspondence: ftecles@um.es; Tel.: +34-868887082 Simple Summary: Saliva is a sample with a high potential in pigs since it is usually easy to obtain and its collection from animals causes less stress than blood sampling. However, the possible effects of daily variations in many salivary biomarkers are still unknown in this species. In our report, the possible variations depending on the sampling time in the day in a panel of 26 salivary biomarkers related to stress, immune system, redox status and other physiological functions in the saliva of pigs Citation: Ortín-Bustillo, A.; were evaluated. In our experimental conditions, daily variations were observed in cortisol, α-amylase, Contreras-Aguilar, M.D.; Rubio, C.P.; total esterase, butyrylcholinesterase, lipase, adenosine deaminase isoenzyme 1, uric acid, superoxide Botia, M.; Cerón, J.J.; López-Arjona, dismutase, alanine aminotransferase, aspartate aminotransferase, alkaline phosphatase, creatine M.; Martínez-Subiela, S.; Escribano, kinase, lactate dehydrogenase, triglycerides and lactate. In some analytes, these differences appeared D.; Tecles, F. Evaluation of the Effect in both sexes, whereas others only showed differences in one sex. These variations should be taken of Sampling Time on Biomarkers of into consideration for an appropriate interpretation of these analytes in the saliva of healthy pigs. Stress, Immune System, Redox Status and Other Biochemistry Analytes in Abstract: This study aims to evaluate the possible variations due to the sampling time in the day in Saliva of Finishing Pigs. Animals 2022, 12, 2127. https://doi.org/ 26 analytes of pigs’ saliva, related to stress, the immune system, redox status and other biomarkers 10.3390/ani12162127 related to metabolism and selected tissues and organs, in order to know the possible effects of the hour of the day in their interpretation. These analytes were measured in saliva obtained from a Academic Editors: Arianna Miglio population of 40 clinically healthy pigs from 8 a.m. to 8 p.m., every 4 h in the same day. In our and Angelica Stranieri experimental conditions, daily variations were observed in cortisol, salivary α-amylase, total esterase Received: 22 July 2022 activity, butyrylcholinesterase, lipase, adenosine deaminase isoenzyme 1, uric acid, superoxide Accepted: 16 August 2022 dismutase, alanine aminotransferase, aspartate aminotransferase, alkaline phosphatase, creatine Published: 19 August 2022 kinase, lactate dehydrogenase, lactate and triglycerides. These changes appeared in both sexes, except Publisher’s Note: MDPI stays neutral for adenosine deaminase isoenzyme 1 and superoxide dismutase which only showed differences with regard to jurisdictional claims in in females. In conclusion, this report indicates that, in the experimental conditions of this trial, published maps and institutional affil- the time of the day and sex can influence the values obtained in various salivary analytes in pigs. iations. These variations should thus be taken into consideration for an adequate interpretation of these analytes when used for the evaluation of health and welfare in this species. Keywords: analytes; biomarkers; fattening pigs; hour of the day; saliva Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and 1. Introduction conditions of the Creative Commons Attribution (CC BY) license (https:// Saliva is gaining attention and interest in veterinary science and, in particular, in pig creativecommons.org/licenses/by/ production as a useful sample to evaluate different biomarkers related to health and 4.0/). welfare. This is because its collection is much easier than blood and does not induce Animals 2022, 12, 2127. https://doi.org/10.3390/ani12162127 https://www.mdpi.com/journal/animals Animals 2022, 12, 2127 2 of 12 animal stress or discomfort [1]. In pigs, a comprehensive panel of analytes including biomarkers of stress, inflammation, immune system, redox status and other analytes related to different physiological aspects and various tissues and organs, such as hepatic and muscle integrity, can currently be measured in this fluid [2]. This is included under the concept of sialochemistry [3,4] which has previously been applied to evaluate changes in analytes during pregnancy, farrowing, lactation [5] and fattening [6] in pigs and also in other species such as horses and cows [7,8]. The knowledge of the possible influence of the sample collection time during the day in salivary analytes is of high importance for an adequate interpretation of the biomarkers [9]. For example, in healthy horses, the values found at certain times of the day for some salivary analytes such as total esterase activity (TEA) or butyrylcholinesterase (BChE) reached those that were found in psychological acute stress situations or diseases [10]. In pigs, cortisol or immunoglobulin A (IgA) showed significant changes in saliva depending on the hour of the day [11–13]. However, no differences have been found in certain salivary analytes, such as chromogranin A [14]. In addition, C-reactive protein and haptoglobin can be measured in the saliva of pigs at any time of the day without significant differences [15]. Despite these previous reports, to our knowledge there are no studies about the variations depending on the hour of the day in other salivary biomarkers of stress different than cortisol, such as α-amylase (sAA), TEA, BChE or lipase (Lip). In addition, no studies about this topic were found for adenosine deaminase (ADA), which is a biomarker of the immune system in pigs [16], for biomarkers of redox status that can be measured in pigs [17], and other analytes included in the sialochemistry profile [2]. Therefore, this study aims to evaluate: (1) if the values from a panel of salivary an- alytes can change when measured at different times of the day in pigs; and (2) if these changes could be related to sex. This panel will be integrated by a total of 26 analytes including biomarkers of stress such as cortisol, sAA, TEA, BChE, and Lip; biomarkers of immune system such as ADA; biomarkers related to redox status such as the cupric reduc- ing antioxidant capacity (CUPRAC), the Trolox equivalent antioxidant capacity (TEAC), the ferric reducing ability of saliva (FRAS), superoxide dismutase (SOD), uric acid (UA), hydrogen peroxide (peroxidase activity Pox-Act), reactive oxygen-derived compounds (d-ROMs) and the advanced oxidation protein products (AOPP); and other analytes re- lated to metabolism and physiological status of different tissue and organs such as alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), γ-glutamine transferase (GGT), lactate dehydrogenase (LDH), creatine kinase (CK), urea, creatinine, triglycerides, total protein, glucose, lactate, calcium and phosphorus. 2. Materials and Methods 2.1. Animals The experiment was conducted on a commercial farm in the southeast of Spain in May 2022. All pigs included in the study were LargeWhite. A total of 40 animals were randomly selected for inclusion in the experimental procedure. They were 20 males and 20 females with an age of 5 months and 1 week, approximately. Animals were kept under general commercial housing and conditions conforming to European Union guidelines (Directive 2010/63/EU1), and had ad libitum access to a nutritionally balanced diet and water. Saliva was sampled from the pigs at 0800 h, 1200 h, 1600 h and 2000 h (8 a.m., 12 a.m., 4 p.m. and 8 p.m., respectively) within a same day. The temperature of the pens at these sampling time-points varied from 19 ◦ C in the morning to 24 ◦ C at midday. Saliva was collected using Salivette tubes (Salivette, Sarstedt, Aktiengesellschaft & Co, Nümbrecht, Germany) as reported previously [6]. Pigs were allowed to chew a sponge attached to a flexible thin metal rod for 1 min. Each sponge was then placed individually in the Salivette tubes and stored on ice until arrival at the laboratory (less than 1 h). At the laboratory, the tubes were centrifuged at 3000× g for 10 min. Saliva samples were transferred into 1.5-mL Eppendorf tubes and frozen at −80 ◦ C for the analysis. Animals 2022, 12, 2127 3 of 12 2.2. Stress Biomarkers All the assays used for stress biomarkers have been previously validated in porcine saliva. Salivary cortisol was measured by an in-house indirect competitive AlphaLisa method [18]. sAA activity was measured by an spectrophotometric assay (a-Amylase, OSR6182, Beckman Coulter Inc., Fullerton, CA, USA) [19]. TEA was analysed by an automated assay using 4-nitrophenyl acetate (4-NA) as substrate [20]. BChE was measured with 5,50 -dithiobis-2- nitrobenzoic acid (DTNB, Sigma-Aldrich, Merck KGaA, Darmstadt, Germany; final concentra- tion, 0.29 mM) as chromophore and butyrylthiocholine iodide (BTCI, Sigma-Aldrich; final concentration, 16.18 mM) as substrate [21]. Lip was determined by a colorimetric assay (Lipase, Beckman Coulter, Brea, CA, USA). Salivary sAA, BChE and Lip were measured in an automated analyser (Olympus AU600, Olympus Diagnostica GmbH, Beckman Coulter). 2.3. Immune System Biomarkers ADA activity was measured by a kinetic spectrophotometric assay (Adenosine Deami- nase assay kit, Diazyme Laboratories, Poway, CA, USA), validated in pigs [22]. Total ADA and ADA2 isoenzyme were measured, respectively, in absence and presence of the specific ADA1 inhibitor erythro-9-(2-hydroxy-3-nonyl) adenine (EHNA, Sigma-Aldrich, City of Saint Louis, MO, USA) in the reaction media at a final concentration of 0.47 mM, as previ- ously described [16]. The isoenzyme ADA1 was calculated from the difference between total ADA and ADA2. ADA assays were carried out in the Olympus AU600 analyser. 2.4. Redox Biomarkers CUPRAC was determined based in the reduction of Cu2+ to Cu1+ by the non-enzymatic antioxidants of the specimen [23]. TEAC was measured by the reduction of 2,20-azino-bis(3 ethylbenzothiazoline-6-sulfonic acid) radical [24]. FRAS was determined by the reduction of the ferric-tripyridyltriazine (Fe3+ -TPTZ) complex [25]. UA was measured by a commercially available colorimetric assay (Beckman Coulter). SOD was measured by a commercially available assay (SOD, Randox Laboratories Ltd., Crumlin, UK). The Pox-Act determination was based on the assay described by Tatzber et al. in which the oxidation of 3,5,30 ,50 - Tetramethylbenzidine (TMB) by peroxides in the sample is monitored. The results were also expressed in µmol of peroxide hydrogen per L of the sample (µmol/L) [26]. Salivary d- ROMs levels were determined based on monitoring the N,N-Diethyl-p-phenylenediamine radical cation concentration as previously described with results expressed in Carratelli units (Carr units) [27]. AOPP were assessed based on a previous reported assay [28]. All the assays used for redox biomarkers measurements were performed in an automated analyser (Olympus AU600), and had been validated in saliva of pigs [17]. 2.5. Analytes Related with Metabolism and Different Tissues and Organs AST, ALP, GGT, LDH, CK, urea, creatinine, triglycerides, glucose, lactate, calcium and phosphorous were measured by spectrophotometric methods (Beckman Coulter). Total protein (TP) was analysed by an assay for urinary proteins quantification (Protein in Urine and CSF, Spinreact, Barcelona, Spain). All the assays showed an imprecision lower than 15%, were linear after serial sample dilution and were performed in an Olympus AU600 biochemistry autoanalyzer. 2.6. Statistical Analysis A Shapiro–Wilk test was used to evaluate data distribution. All analytes showed a non-normal distribution. Therefore, they were log-transformed by applying the formula ln x = ln(x + 1) for further analysis [10]. After transformation, significant daily changes in the salivary biomarkers measured were evaluated by a mixed linear model of repeated measures considering fixed factors of time (from 0800 h to 2000 h) and gender, being the individual considered as a random factor. When any of the fixed factors was significant, it was further evaluated with a univariate analysis and Bonferroni post-hoc test. Changes were considered significant when p < 0.05. The statistical analyses were calculated using Animals 2022, 12, 2127 4 of 12 Graph Pad Prism 8 (GraphPad Software, San Diego, CA, USA) and SPSS statistics package (IBM SPSS Statistics for Windows, Version 28.0.1. IBM Corp, Armonk, NY, USA). 3. Results The results for the stress and immune system biomarkers appear in Figure 1. Cortisol, sAA, TEA, BChE and Lip showed a significant time effect, with an increase at 1600 h for cortisol, sAA, TEA and BChE, whereas Lip showed a significant increase at 2000 h. ADA1 isoenzyme showed significant changes for time and gender fixed effects, with a significant increase at 1600 h only in females. ADA2 isoenzyme values showed a significant gender effect, with higher values at 2000 h in females than in males. Figure 1. Results for the stress (cortisol, salivary α-amylase (sAA), total esterase activity (TEA), Figure 1. Results for the stress (cortisol, salivary α-amylase (sAA), total esterase activity (TEA), butyrylcholinesterase butyrylcholinesterase (BChE) and lipase (BChE) and lipase(Lip)) (Lip))and and immune immune system system (adenosine (adenosine deaminase deaminase (ADA) (ADA) isoenzymes isoenzymes11andand 2) biomarkersobtained 2) biomarkers obtained in in saliva saliva fromfrom 40 healthy 40 healthy fattening fattening pigs (20pigs (20 females, females, 20 20 males) males)obtained obtainedatat different hours different of the hours day.day. of the pB: ppvalue of theof B : p value Bonferroni post-hoc the Bonferroni analysis.analysis. post-hoc Animals 2022, 12, 2127 5 of 12 Results of the Redox status biomarkers are shown in Figure 2. No significant effects were seen in most of the biomarkers between the different hours of the day; with the Animals 2022, 12, x exception of UA which showed a significant time effect with an elevation 7atof2000 16 h, and SOD that showed a significant time and gender interaction with values showing significant increases at 1600 h only in females. Figure 2. Results for the Redox statuss biomarkers (cupric reducing antioxidant capacity (CUPRAC), Figure 2. Results for the Redox statuss biomarkers (cupric reducing antioxidant capacity (CU- the Trolox PRAC),equivalent antioxidant the Trolox equivalent capacity antioxidant (TEAC), capacity the the (TEAC), ferric reducing ferric reducingability of saliva ability of saliva (FRAS), (FRAS), uric acid uric superoxide (UA), acid (UA), superoxide dismutase dismutase (SOD),(SOD), the advanced the advanced oxidationprotein oxidation protein products products (AOPP), reactive oxygen-derived compounds (d-ROMs) and hydrogen peroxide (peroxidase activity Pox-Act)) obtained in saliva from 40 healthy fattening pigs (20 females, 20 males) obtained at different hours of the day. pB : p value of the Bonferroni post-hoc analysis. Animals 2022, 12, 2127 6 of 12 Results of enzymatic markers related to the physiological status of different tissue and organs are shown in Figure 3. A significant time effect was detected in ALT, AST, ALP and Animals 2022, 12, x CK, whose values were higher at 1600 h if compared with those obtained at 1200 9 of 16 h, and also for LDH, with higher values at 0800 h if compared with values at 2000 h. No gender effect was observed in any of those enzymatic activities. Figure 3. Results for the enzymatic activities (alanine aminotransferase (ALT), aspartate aminotrans- Figure 3. Results for the enzymatic activities (alanine aminotransferase (ALT), aspartate ami- ferase (AST), alkaline notransferase (AST), phosphatase (ALP), alkaline phosphatase γ-glutamine (ALP), γ-glutamine transferase transferase(GGT), lactatedehydro- (GGT), lactate dehydrogenase (LDH)genase (LDH) and creatine kinase (CK)) obtained in saliva from 40 healthy fattening pigs (20 (20 and creatine kinase (CK)) obtained in saliva from 40 healthy fattening pigs fe- females, males, 20 males) obtained at different hours of the day. p B: p value of the Bonferroni post-hoc 20 males) obtained at different hours of the day. pB : p value of the Bonferroni post-hoc analysis. anal- ysis. Regarding metabolites (Figure 4), no significant effects were observed for urea, crea- tinine, glucose, Ca and P. Lactate and triglycerides showed a significant time effect, with lactate values being higher at 0800 h compared to 1200 h and at 1600 h, whereas triglycerides values at 0800 h were lower than at 1600 h. Animals 2022, 12, x 11 of 16 Animals 2022, 12, 2127 7 of 12 Figure 4. Results for the metabolites (urea, creatinine, glucose, lactate, triglycerides, total protein (TP), calcium (Ca) and phosphorous (P)) obtained in saliva from 40 healthy fattening pigs (20 females, 20 males) obtained at different hours of the day. pB : p value of the Bonferroni post-hoc analysis. 4. Discussion This study aimed to evaluate if the values of a comprehensive panel of analytes measured in saliva of pigs can vary at different hours of the day and whether or not these differences could depend on the sex. These data could be of practical use in order to select Animals 2022, 12, 2127 8 of 12 the most appropriate saliva sampling times when measuring analytes for the evaluation of pig health and welfare. In this report, we performed four different samplings during the day in order to cover the different times in which sampling could be done in routine in farm conditions and according to previous studies made in pigs [13] in which four time points of the day were also selected. From the 26 analytes measured, to the authors’ knowledge, only cortisol has been previously analysed during different times of the day in saliva of pigs. Similarly to our results, cortisol has been described to peak in porcine saliva at 1600 h in other studies [12,29]. In a report in which pigs were fed two times per day, it was postulated that the increase in salivary cortisol at 1600 h could be explained by increased gluconeogenesis during the fasting periods (11). However, in our case, the pigs were feed ad libitum, so we were not able to determine when the pigs ingested food. In any case, this increase in cortisol at this time has been also described with pigs feeding ad libitum [29] and further studies should be made to elucidate the reason. Overall, the increase of cortisol found in our study at 1600 h, as previously described, supports the validity of our experimental model. Interestingly, in addition to cortisol, other analytes related to stress such as sAA, TEA and BChE showed a significant increase at 1600 h in our experimental conditions. Although sAA and BChE are more related with the autonomous nervous system [21], whereas cortisol is related with the hypothalamus-hypophysis axis, maybe there could be some mechanism that could potentially activate both axes at this time of the day, and lead to the increase in these analytes. ADA1 isoenzyme showed a significant increase at 1600 h in females. In human erythro- cytes, higher values of ADA were found in samples obtained in the afternoon in comparison with the morning [30]. The reason for this change should be elucidated, and possibly it could have relation with the lymphocyte dynamics, since ADA is an enzyme related with lymphocyte function [31]. In this line, in humans, numbers of total lymphocytes, T and B cells showed higher levels in the afternoon–night than in the morning [32]. Regarding markers of redox homeostasis, UA was the only analyte showing significant differences in both sexes, which a significant elevation at 2000 h. This would be in line with a previous study in humans reporting that the maximal UA concentration, more than two times greater than the minimal values, were measured during the night [33]. A high antioxidant salivary potential could be necessary during the night, when the concentrations of compounds with oxidant potential seem to be at their maximum [34,35]. In metabolites and markers of different tissues and organs, various enzymes such as ALT, AST, ALP and CK showed higher values at 1600 h. To the authors’ knowledge there are not previous studies about daily variations in these analytes and further research should be performed to evaluate the reasons for these changes. Regarding the possible practical applications of these biomarkers, in other species such as dogs, increases in urea and creatinine in saliva have been described in cases of renal failure [36] and increases in CK in muscle damage [37], with theses analytes showing high correlations between saliva and serum. These are examples indicating that in some situations, saliva could be an alternative to serum as biological sample for the measurement of selected metabolites. However, specific studies in the case of the porcine species should be performed to evaluate this. The gender of the animals did not affect most of the results, except for ADA isoenzymes and SOD activities, which were significantly higher in females than in males at some of the sampling points. Overall, it could be indicated that in general sex does not have a major influence in the variations of salivary analytes during the day with some selected exceptions such as ADA or SOD. Some analytes such as sAA and ADA showed large individual differences specially at 1600 h. In the case of salivary analytes, it has been reported that those with lower in- terindividual variability and the greatest magnitude of change with absence of overlapping between the conditions to be evaluated would be preferred [5]. For those biomarkers showing a high interindividual variability, each animal should ideally act as its own control. Animals 2022, 12, 2127 9 of 12 In addition, based on the findings of this report it could be added that the absence of variations depending on the hour of the day would be a positive factor for a biomarker, since it could be used in samples collected at any time of the day. This would be the case for ADA2 within the immune system biomarkers, GGT within the enzymatic activity markers, and several of the biomarkers for redox homeostasis, which will have the advantage of a lack of changes along the day. ADA has been proposed as a biomarker of inflammation and pain in pigs, showing increases in lameness and prolapses [38] and also increased in sows after parturition, potentially as a response to the postpartum inflammation [39]. GGT has been proven to be a good pain biomarker in species such as the horse [7] and also increases in sows after parturition [5]. In addition, some of the biomarkers of Redox status studied in the manuscript, such as CUPRAC, FRAS, AOPP, Pox-Act and d-ROMs, showed higher values in pigs with sepsis compared to healthy pigs [40]. It is important to point out that in our experimental conditions, in some cases such as cortisol, BChE, Lip and uric acid, the increases found at their peak hours of the day reached similar values than described in pigs with diseases and clinical signs such as lameness and prolapses [38], and this could influence their practical interpretation. Some of analytes that were used in this panel were significantly correlated with coeffi- cients of correlation higher than 0.7 (data not shown). For example, urea and creatinine which are biomarkers of renal function were correlated. Furthermore, significant correla- tions were found between the metabolites glucose, triglycerides and lactate; and between sAA and BChE, which are markers of acute stress [21]. In addition, among the antioxidant markers, FRAS was correlated with TEAC, and the oxidants Pox-Act and d-ROMS were correlated. However, it would be of interest to evaluate how these analytes correlate also in situations of stress and disease. This report has various limitations. First, it should be indicated that the results reported here have been obtained in our experimental conditions and with our analytical methods, and care should be taken when applied to other conditions or different assays. In addition, it should be considered as a pilot study to evaluate the possibility of the presence of changes in salivary analytes depending on the hour of the day and sex. Therefore, further studies should be made in the future at different times of the year and also with more sampling times, in order to obtain a more complete picture of the changes occurring and refine our findings and their possible applications. Studying the possible differences due to season would be of particular importance, since in humans, for example, differences have been found in ADA in winter versus summer [30]. In addition, the possible variations depending on the breed, age, season and productive conditions should be evaluated. Moreover, studying the possible changes depending on different situations of disease or stress would be of interest. Finally, the establishment of proper reference values considering all these sources of variation would be desirable. For this purpose, further studies involving a higher number of animals, at least 120 in case of non-parametric reference intervals [41], should be performed. 5. Conclusions In our experimental conditions there were changes in the values of selected analytes in saliva of pigs depending on the hour of the day in which the sample is collected and sex. Namely, daily variations were observed in cortisol, sAA, TEA, BChE, Lip, ADA1, UA, SOD, ALT, AST, ALP, CK, LDH, lactate and triglycerides. In some analytes, the differences appeared in both sexes, whereas others only showed differences in one sex. Additional studies should be made to elucidate the mechanisms underlying these changes and evaluate if other possible sources of variation such as breed, age, season or productive condition could influence them. However, based on our results, for the measurement of the analytes that showed changes in our study, it would be recommended to obtain saliva samples in the morning; and if afternoon samples should be used, different reference ranges should be established. Animals 2022, 12, 2127 10 of 12 Author Contributions: Conceptualization and design, J.J.C., D.E., S.M.-S. and F.T.; formal analysis, A.O.-B., M.D.C.-A., M.L.-A., S.M.-S., D.E., M.B. and C.P.R.; writing—original draft preparation, A.O.-B., J.J.C. and F.T.; writing—review and editing, J.J.C., D.E. and F.T.; supervision, J.J.C., D.E. and F.T.; project administration, J.J.C.; funding acquisition, J.J.C. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by MCIN/AEI/10.13039/501100011033 (Grant reference PCI2020- 120712-2) and European Union “NextGenerationEU”/PRTR (1st ICRAD Joint Cofund Call). A.O.-B. was funded by Fundación Séneca (21603/FPI/21), Región de Murcia (Spain). D.E. was funded by the postdoctoral contract “Generational renewal to promote research” of the University of Murcia. C.P.R. has a post-doctoral fellowship “Juan de la Cierva Formación” supported by the “Ministerio de Economía y Competitividad” (FJC2019-042475-I), Spain. Institutional Review Board Statement: The animal study protocol was Bioethical Commission of Murcia University according to the European Council Directives regarding the protection of animals used for experimental purposes (Approval number, 235/2016; Approval date, 25 April 2016). Data Availability Statement: The data presented in this study are available on request from the corresponding author. Conflicts of Interest: The authors declare no conflict of interest. References 1. Cerón, J. Acute Phase Proteins, Saliva and Education in Laboratory Science: An Update and Some Reflections. BMC Vet. Res. 2019, 15, 197. [CrossRef] 2. Cerón, J.; Contreras-Aguilar, M.; Escribano, D.; Martínez-Miró, S.; López-Martínez, M.; Ortín-Bustillo, A.; Franco-Martínez, L.; Rubio, C.; Muñoz-Prieto, A.; Tvarijonaviciute, A.; et al. Basics for the Potential Use of Saliva to Evaluate Stress, Inflammation, Immune System, and Redox Homeostasis in Pigs. BMC Vet. Res. 2022, 18, 81. [CrossRef] 3. Giri, U.; Nagaraj, V.; Sankaran, A.; Ramesh, R.; Rajaram, S.; Muthanandam, S.; Arumugam, S.; Santhanam, V. Sialochemical Profile in Depressive Individuals under Antidepressant Therapy: An Observational Study. J. Clin. Diagn. Res. 2018, 12, VC06–VC09. [CrossRef] 4. Miller, B.; Deutsch, O.; Redlich, M.; Konttinen, Y.; Benoliel, R.; Zaks, B.; Davidovich, E.; Palmon, A.; Aframian, D. Sialochemistry and Cortisol Levels in Patients with Sjogren’s Syndrome. Oral Dis. 2012, 18, 255–259. [CrossRef] [PubMed] 5. Contreras-Aguilar, M.; López-Arjona, M.; Martínez-Miró, S.; Escribano, D.; Hernández-Ruipérez, F.; Cerón, J.; Tecles, F. Changes in Saliva Analytes during Pregnancy, Farrowing and Lactation in Sows: A Sialochemistry Approach. Vet. J. 2021, 273, 105679. [CrossRef] [PubMed] 6. Ortín-Bustillo, A.; Escribano, D.; López-Arjona, M.; Botia, M.; Fuentes, P.; Martínez-Miró, S.; Rubio, C.P.; García-Manzanilla, E.; Franco-Martínez, L.; Pardo-Marín, L.; et al. Changes in a Comprehensive Profile of Saliva Analytes in Fattening Pigs during a Complete Productive Cycle: A Longitudinal Study. Animals 2022, 12, 1865. [CrossRef] 7. Contreras-Aguilar, M.; Escribano, D.; Martínez-Subiela, S.; Martín-Cuervo, M.; Lamy, E.; Tecles, F.; Cerón, J. Changes in Saliva Analytes in Equine Acute Abdominal Disease: A Sialochemistry Approach. BMC Vet. Res. 2019, 15, 187. [CrossRef] 8. Contreras-Aguilar, M.; Vallejo-Mateo, P.; Želvytė, R.; Tecles, F.; Rubio, C. Changes in Saliva Analytes Associated with Lameness in Cows: A Pilot Study. Animals 2020, 10, 2078. [CrossRef] 9. Humphrey, S.P.; Williamson, R.T. A Review of Saliva: Normal Composition, Flow, and Function. J. Prosthet. Dent. 2001, 85, 162–169. [CrossRef] 10. Contreras-Aguilar, M.D.; Lamy, E.; Escribano, D.; Cerón, J.J.; Tecles, F.; Quiles, A.J.; Hevia, M.L. Changes in Salivary Analytes of Horses Due to Circadian Rhythm and Season: A Pilot Study. Animals 2020, 10, 1486. [CrossRef] 11. Muneta, Y.; Yoshikawa, T.; Minagawa, Y.; Shibahara, T.; Maeda, R.; Omata, Y. Salivary IgA as a Useful Non-Invasive Marker for Restraint Stress in Pigs. J. Vet. Med. Sci. 2010, 72, 1295–1300. [CrossRef] [PubMed] 12. Hillmann, E.; Schrader, L.; Mayer, C.; Gygax, L. Effects of Weight, Temperature and Behaviour on the Circadian Rhythm of Salivary Cortisol in Growing Pigs. Animal 2008, 2, 405–409. [CrossRef] [PubMed] 13. Gallagher, N.L.; Giles, L.R.; Wynn, P.C. The Development of a Circadian Pattern of Salivary Cortisol Secretion in the Neonatal Piglet. Biol. Neonate 2002, 81, 113–118. [CrossRef] 14. Escribano, D.; Gutiérrez, A.; Fuentes-Rubio, M.; Cerón, J. Saliva Chromogranin A in Growing Pigs: A Study of Circadian Patterns during Daytime and Stability under Different Storage Conditions. Vet. J. 2014, 199, 355–359. [CrossRef] [PubMed] 15. Gutiérrez, A.M.; Escribano, D.; Fuentes, M.; Cerón, J.J. Circadian Pattern of Acute Phase Proteins in the Saliva of Growing Pigs. Vet. J. 2013, 196, 167–170. [CrossRef] 16. Contreras-Aguilar, M.; Tvarijonaviciute, A.; Monkeviciene, I.; Martín-Cuervo, M.; González-Arostegui, L.; Franco-Martínez, L.; Cerón, J.; Tecles, F.; Escribano, D. Characterization of Total Adenosine Deaminase Activity (ADA) and Its Isoenzymes in Saliva and Serum in Health and Inflammatory Conditions in Four Different Species: An Analytical and Clinical Validation Pilot Study. BMC Vet. Res. 2020, 16, 384. [CrossRef] Animals 2022, 12, 2127 11 of 12 17. Rubio, C.; Mainau, E.; Cerón, J.; Contreras-Aguilar, M.; Martínez-Subiela, S.; Navarro, E.; Tecles, F.; Manteca, X.; Escribano, D. Biomarkers of Oxidative Stress in Saliva in Pigs: Analytical Validation and Changes in Lactation. BMC Vet. Res. 2019, 15, 144. [CrossRef] 18. López-Arjona, M.; Escribano, D.; Mateo, S.V.; Contreras-Aguilar, M.D.; Rubio, C.P.; Tecles, F.; Cerón, J.J.; Martínez-Subiela, S. Changes in Oxytocin Concentrations in Saliva of Pigs after a Transport and during Lairage at Slaughterhouse. Res. Vet. Sci. 2020, 133, 26–30. [CrossRef] 19. Fuentes, M.; Tecles, F.; Gutiérrez, A.; Otal, J.; Martínez-Subiela, S.; Cerón, J. Validation of an Automated Method for Salivary Alpha-Amylase Measurements in Pigs (Sus Scrofa Domesticus) and Its Application as a Stress Biomarker. J. Vet. Diagn. Investig. 2011, 23, 282–287. [CrossRef] 20. Tecles, F.; Contreras-Aguilar, M.D.; Martínez-Miró, S.; Tvarijonaviciute, A.; Martínez-Subiela, S.; Escribano, D.; Cerón, J.J. Total Esterase Measurement in Saliva of Pigs: Validation of an Automated Assay, Characterization and Changes in Stress and Disease Conditions. Res. Vet. Sci. 2017, 114, 170–176. [CrossRef] 21. Tecles, F.; Escribano, D.; Martínez-Miró, S.; Hernández, F.; Contreras, M.D.; Cerón, J.J. Cholinesterase in Porcine Saliva: Analytical Characterization and Behavior after Experimental Stress. Res. Vet. Sci. 2016, 106, 23–28. [CrossRef] [PubMed] 22. Tecles, F.; Rubio, C.; Contreras-Aguilar, M.; López-Arjona, M.; Martínez-Miró, S.; Martínez-Subiela, S.; Cerón, J. Adenosine Deaminase Activity in Pig Saliva: Analytical Validation of Two Spectrophotometric Assays. J. Vet. Diagn. Investig. 2018, 30, 175–179. [CrossRef] [PubMed] 23. Campos, C.; Guzmán, R.; López-Fernández, E.; Casado, A. Evaluation of the Copper(II) Reduction Assay Using Bathocuproinedisulfonic Acid Disodium Salt for the Total Antioxidant Capacity Assessment: The CUPRAC-BCS Assay. Anal. Biochem. 2009, 392, 37–44. [CrossRef] 24. Arnao, M.; Cano, A.; Hernández-Ruiz, J.; García-Cánovas, F.; Acosta, M. Inhibition by L-Ascorbic Acid and Other Antioxidants of the 2.2’-Azino-Bis(3-Ethylbenzthiazoline-6-Sulfonic Acid) Oxidation Catalyzed by Peroxidase: A New Approach for Determining Total Antioxidant Status of Foods. Anal. Biochem. 1996, 236, 255–261. [CrossRef] 25. Benzie, I.; Strain, J. The Ferric Reducing Ability of Plasma (FRAP) as a Measure of “Antioxidant Power”: The FRAP Assay. Anal. Biochem. 1996, 239, 70–76. [CrossRef] [PubMed] 26. Tatzber, F.; Griebenow, S.; Wonisch, W.; Winkler, R. Dual Method for the Determination of Peroxidase Activity and Total Peroxides-Iodide Leads to a Significant Increase of Peroxidase Activity in Human Sera. Anal. Biochem. 2003, 316, 147–153. [CrossRef] 27. Cesarone, M.; Belcaro, G.; Carratelli, M.; Cornelli, U.; De Sanctis, M.; Incandela, L.; Barsotti, A.; Terranova, R.; Nicolaides, A. A Simple Test to Monitor Oxidative Stress. Int. Angiol. 1999, 18, 127–130. 28. Witko-Sarsat, V.; Friedlander, M.; Capeillère-Blandin, C.; Nguyen-Khoa, T.; Nguyen, A.; Zingraff, J.; Jungers, P.; Descamps-Latscha, B. Advanced Oxidation Protein Products as a Novel Marker of Oxidative Stress in Uremia. Kidney Int. 1996, 49, 1304–1313. [CrossRef] 29. Ruis, M.; Te Brake, J.; Engel, B.; Ekkel, E.; Buist, W.; Blokhuis, H.; Koolhaas, J. The Circadian Rhythm of Salivary Cortisol in Growing Pigs: Effects of Age, Gender, and Stress. Physiol. Behav. 1997, 62, 623–630. [CrossRef] 30. Cornélissen, G.; Touitou, Y.; Tritsch, G.; Bogdan, A.; Auzéby, A.; Reinberg, A.; Halberg, F. Circadian Rhythms of Adenosine Deaminase Activity in Human Erythrocytes: A Transverse Study on Young, Elderly and Senile Demented Subjects. Ric. Clin. Lab. 1985, 15, 365–374. [CrossRef] 31. Macdermott, R.P.; Tritsch, G.L.; Formeister, J.F.; York, N. Adenosine Deaminase and Nucleoside Phosphorylase Activities in Normal Human Blood Mononuclear Cell Subpopulations. Clin. Exp. Immunol. 1980, 42, 303–307. [PubMed] 32. Bertouch, J.V.; Roberts Thomson, P.J.; Bradley, J. Diurnal Variation of Lymphocyte Subsets Identified by Monoclonal Antibodies. Br. Med. J. (Clin. Res. Ed.) 1983, 286, 1171–1172. [CrossRef] [PubMed] 33. Castagnola, M.; Cabras, T.; Denotti, G.; Fadda, M.B.; Gambarini, G.; Lupi, A.; Manca, I.; Onnis, G.; Piras, V.; Soro, V.; et al. Circadian Rhythms of Histatin 1, Histatin 3, Histatin 5, Statherin and Uric Acid in Whole Human Saliva Secretion. Biol. Rhythm Res. 2010, 33, 213–222. [CrossRef] 34. Ferguson, D.B.; Fort, A. Circadian Variations in Human Resting Submandibular Saliva Flow Rate and Composition. Arch. Oral Biol. 1974, 19, 47–55. [CrossRef] 35. Ferguson, D.B.; Botchway, C.A. Circadian Variations in the Flow Rate and Composition of Whole Saliva Stimulated by Mastication. Arch. Oral Biol. 1979, 24, 877–881. [CrossRef] 36. Tvarijonaviciute, A.; Pardo-Marin, L.; Tecles, F.; Carrillo, J.; Garcia-Martinez, J.; Bernal, L.; Pastor, J.; Cerón, J.; Martinez-Subiela, S. Measurement of Urea and Creatinine in Saliva of Dogs: A Pilot Study. BMC Vet. Res. 2018, 14, 223. [CrossRef] 37. Tvarijonaviciute, A.; Barranco, T.; Rubio, M.; Carrillo, J.; Martinez-Subiela, S.; Tecles, F.; Carrillo, J.; Cerón, J. Measurement of Creatine Kinase and Aspartate Aminotransferase in Saliva of Dogs: A Pilot Study. BMC Vet. Res. 2017, 13, 168. [CrossRef] 38. Contreras-Aguilar, M.; Escribano, D.; Martínez-Miró, S.; López-Arjona, L.; Rubio, C.; Martínez-Subiela, S.; Cerón, J.; Tecles, F. Application of a Score for Evaluation of Pain, Distress and Discomfort in Pigs with Lameness and Prolapses: Correlation with Saliva Biomarkers and Severity of the Disease. Res. Vet. Sci. 2019, 126, 155–163. [CrossRef] Animals 2022, 12, 2127 12 of 12 39. Kaiser, M.; Jacobson, M.; Andersen, P.; Bækbo, P.; Cerón, J.; Dahl, J.; Escribano, D.; Jacobsen, S. Inflammatory Markers before and after Farrowing in Healthy Sows and in Sows Affected with Postpartum Dysgalactia Syndrome. BMC Vet. Res. 2018, 14, 83. [CrossRef] 40. López-Martínez, M.J.; Escribano, D.; Ortín-Bustillo, A.; Franco-Martínez, L.; González-Arostegui, L.G.; Cerón, J.J.; Rubio, C.P. Changes in Biomarkers of Redox Status in Saliva of Pigs after an Experimental Sepsis Induction. Antioxidants 2022, 11, 1380. [CrossRef] 41. Horn, P.S.; Pesce, A.J. Reference Intervals: An Update. Clin. Chim. Acta 2003, 334, 5–23. [CrossRef]