{'fraction_non_alphanumeric': 0.03260043087061209, 'fraction_numerical': 0.014573564820681789, 'mean_word_length': 5.236909701639992, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Midday and nadir salivary cortisol appear superior to cortisol awakening response in burnout assessment and monitoring

Burnout and work-related stress symptoms of anxiety disorder and depression cause prolonged work absenteeism and early retirement. Hence, reliable identification of patients under risk and monitoring of treatment success is highly warranted. We aimed to evaluate stress-specific biomarkers in a population-based, “real-world” cohort (burnouts: n = 40, healthy controls: n = 26), recruited at a preventive care ward, at baseline and after a four-month follow up, during which patients received medical and psychological treatment. At baseline, significantly higher levels of salivary cortisol were observed in the burnout group compared to the control group. This was even more pronounced in midday- (p < 0.001) and nadir samples (p < 0.001) than for total morning cortisol secretion (p < 0.01). The treatment program resulted in a significant reduction of stress, anxiety, and depression scores (all p < 0.001), with 60% of patients showing a clinically relevant improvement. This was accompanied by a ~30% drop in midday cortisol levels (p < 0.001), as well as a ~25% decrease in cortisol nadir (p < 0.05), although not directly correlating with score declines. Our data emphasize the potential usefulness of midday and nadir salivary cortisol as markers in the assessment and biomonitoring of burnout.

However, there is no doubt that psychiatric disorders are increasing in working patients, thereby leading to prolonged sick leaves (30-50%) 10 and premature retirement (32.8%) 11 in Austria and other welfare countries. Our cohort of municipal employees typically report "burnout symptoms" reflecting merely work-related stress symptoms of mild to moderate MDD and anxiety disorder within classical DSM-IV boundaries 12 . In contrast to psychiatric inpatients, these primary care patients are typically underdiagnosed and subsequently undertreated. Hence, a population-based model is required to provide information for initial screening and treatment to this often neglected subpopulation at early stages of disease 13 .
As mentioned above, reliable biomarkers for early diagnosis of burnout and monitoring of treatment success are still missing 14 . We recently identified a set of peripherally assessable biomarkers amongst professional orchestra musicians in a performance situation 15 .Those markers have been reported to exhibit immunological messenger or effector function and to act as proxy for psychological stress. Amongst them was cortisol, which has repeatedly been implicated in mechanistic models of major depression or anxiety disorders 16 . Briefly, psychological stress is thought to be processed by the hippocampus, which serves as a brain-environment interface and should optimally prepare the organism for environmental changes. On a hormone level, this was reflected by alterations in the release into the blood circulation of cortisol releasing hormone (CRH). Increased CRH levels impact on brain function resulting in adaptive behavior in healthy subjects or maladaptive behavior and psychopathological symptoms in patients 17 . Moreover, hormonal signaling across the HPA axis also leads to increases in cortisol levels during situations of acute psychological stress or a blunted response 18 ; especially the cortisol awakening response (CAR), which is the rapid increase in cortisol levels across the first 30-45 min following morning awakening, represents a promising and convenient parameter that is associated with depressive symptoms and psychological stress 19 . Additionally, in our musician cohort we found associations between increased psychological stress levels and interleukin 6 (IL6), myeloperoxidase (MPO), and homocysteine (HCYS) 15 . All of these compounds can promote oxidative stress, which may trigger neuropsychiatric disorders by damaging brain structures 20 . In brief, IL6 is a mainly pro-inflammatory cytokine, which also responds to mental stress 21 and is linked to mood disorders through the cytokine model of depression 22 . MPO is the main antimicrobial effector molecule of neutrophilic granulocytes 23 . When released into the surrounding tissue during host defense or in response to psychological stress 24 , its peroxidase activity promotes inflammation, and damages also endogenous structures via generation of reactive oxygen and nitrogen species 25 . HCYS is potentially toxic, especially when it occurs within the nervous system, and is generated by demethylation of methionine 26 . The substance is naturally eliminated by both folate and vitamin B 12 , however, accumulation of HCYS precedes vascular events 27 and correlates with the presence of neuropsychiatric disorders 28,29 . Further support for the utility of these latter three parameters stems from independent studies showing that acute psychological stress leads to an increase of circulating IL6 30 , MPO 31 , and HCYS levels 32 .
Despite the abundance of reports highlighting the association between stress-related psychopathological symptoms and changes in blood parameters, there is still a debate on how such indicators 14 could be utilized for diagnostic purposes. To improve disease management, we assessed cortisol day profiles, but also IL6, MPO, and HCYS plasma concentrations so as to complement clinical psychological measurements in a representative cohort of burnout patients during a four month period at our stress outpatient clinic. We hypothesized higher saliva cortisol levels in so called "burnout" patients prior to clinical response and relative to healthy controls. Moreover, in secondary hypotheses we expected peripheral markers of oxidative stress to be increased ahead of treatment in this patient group possibly useful as future disease marker.

Methods
Patient recruitment and Study design. Outpatients of the Health and Prevention Center, which is the primary Healthcare Institution for municipal employees of the City of Vienna, were screened for depressive, anxiety, and unspecific work-related stress symptoms. Data collection was conducted between 09/2014 and 11/2016 utilizing a comprehensive psychosocial screening questionnaire (PSS) as being detailed below. Patients of both sexes were included beyond classical DSM-IV boundaries in line with the recently proposed Research Domain Criteria (RDoC) approach 12 . For this, the following inclusion criteria were applied: (1) ability to provide written informed consent, (2) commitment to fulfill the study protocol including baseline and follow-up visits, (3) age older than 18 yrs., and (4) exceeding cut-off scores in HADS depression, HADS anxiety (≥8), or a T-value of PSS stress sum-score ≥60. In contrast, (1) pregnancy, (2) non-compliance with the study protocol, or (3) a current treatment with glucocorticoids were defined as exclusion criteria. Clinical axis I main diagnoses were assessed according to DMS-IV by experienced and board-certified psychiatrists after study inclusion (see sample characteristics). All participants underwent a thorough medical exam at baseline, as well as pre-and posttreatment blood withdrawal, PSS assessments, and evaluation of saliva cortisol day profiles. The group of patients was provided an open four-month lasting treatment comprising intensive medical and psychological assistance (see below). In addition, a control cohort, which did not differ in terms of age and sex, without any concurrent psychiatric diagnosis was consecutively recruited at the Health and Prevention Center. Recruitment procedures are depicted in Fig. 1.
This study as well as all applicable amendments were approved by the Ethics Committee of the Medical University of Vienna (071/2009) and was performed in accordance with the Declaration of Helsinki and other relevant guidelines/regulations. All participants gave written informed consent prior to study participation. care service at the Health and Prevention Center, Sanatorium Hera. Following data were collected via paper/ pencil or online (http://service.hera.co.at/psychoscreen/): Depression and anxiety items were derived from the Hospital-Anxiety and Depression Scale (HADS), whose 2-factorial structure is well established and both subscales can be considered highly reliable (retest reliability in the interval of two weeks > 0.8) 33 . Work-related stress was assessed based on the burn-out risk screening according to the 12-phase model of Freudenberger and North 9 . This screening instrument results in a validated factorial structure of twelve sub-scales corresponding to the phases and a sum score transformed in a T-value. It has proved to be highly consistent (Cronbach's α = 0.937) 34 . Finally, lifestyle parameters were assessed including the Personal-Life-Style-Questionnaire (PLQ), the Fagerström nicotine dependency 35 and the Audit-GMAT for alcohol dependency 36 .
Symptom severity (CGI-S) was assessed at each visit utilizing the Clinical Global Impressions (CGI) Scale by a board certified psychiatrist on a scale ranging from 1 ("normal") to 7 ("among the most extremely ill patients"), with 2 as a cut-off for marginally ill patients 37 .
Collection of samples. Participants were fully instructed regarding the preanalytical requirements of saliva collection and asked to adhere to the following procedure the day before baseline and follow-up examinations: At-home collection of salivary samples 1) immediately after awakening, as well as 2) 15 minutes, 3) 30 minutes and 4) 45 minutes after awakening, 5) between 11:30 a.m. and 12:30 p.m., and 6) directly before going to bed. They were not allowed breakfasting (incl. caffeinated or sugared drinks) or smoking during the post-awakening phase. To improve compliance with salivary cortisol sampling we specifically emphasized to the participants the importance of accurate sample collection and reporting their sampling times 38,39 . The protocol was explained in detail to participants, thereby focusing on what is meant by the 'moment of awakening' . In addition, subjects were informed that analysis of post-awakening cortisol levels requires strict adherence to the sampling schedule, since delays of 10 to 15 min might already affect results 40 . To increase compliance, participants also kept a diary to record collection time points and activities during the day. Samples were stored at approximately 4 °C at the patients' residences until they were transferred to the Department of Laboratory Medicine. Blood withdrawal Laboratory analyses. Plasma concentrations of myeloperoxidase (MPO), interleukin 6 (IL6), homocysteine (HCYS), as well as salivary cortisol levels were measured at the Department for Laboratory Medicine, Medical University of Vienna, which, as the central facility for laboratory diagnostics at the General Hospital of Vienna, harbors a certified (ISO 9001:2008) and accredited (ISO 15189:2012, not applicable for myeloperoxidase) quality management system (www.kimcl.at).
In detail, MPO was quantified from EDTA-anticoagulated plasma by means of enzyme-linked immunosorbent assays (ELISA), using commercially available, IVD-certified MPO ELISA kits (Immundiagnostik AG, Bensheim, Germany), featuring 4-5% intra-assay and 12-15% inter-assay reproducibility. IL6 was measured from EDTA-anticoagulated whole-blood using commercially available IL6 electrochemiluminescent immunoassays (ECLIA) kits (Roche Diagnostics GmbH, Mannheim, Germany) on a cobas e602 module (Roche Diagnostics GmbH) with an intra-assay precision between 1.4 and 6.0% and an inter-assay reproducibility between 2.7 and 8.5%. In order to prevent in-vitro production of HCYS by erythrocytes, EDTA-anticoagulated whole blood was cooled after withdrawal and plasma was at the earliest possible time point separated from blood cells. HCYS was quantified by enzyme-cycling assays, using the oxidative capacity of homocysteine on cobas c701 modules (Roche Diagnostics GmbH). Cortisol was assessed from saliva collected in Salivette ® cortisol saliva tubes (Sarstedt AG & Co., Nuembrecht, Germany) by means of ECLIA on cobas e602 modules (Roche Diagnostics GmbH). Treatment Period. The Health and Prevention Center, Sanatorium Hera, provides a collaborative 41 and population-based care 13 by a multidisciplinary team of mental health professionals including psychiatrists, psychologists and psychotherapists, and specialists from other medical fields, such as physiotherapy, sports science and occupational medicine. Collaborative management of outpatients was coordinated weekly. Individual treatment strategies were implemented for specific clinical conditions according to patient's need and illness severity to maximize treatment adherence 42 and to provide a cost-effective case management 13 . Treatment in face-to-face and group settings followed current S3 treatment guidelines of depression 43 in terms of following treatment elements: patient education, pharmacological therapy, psychological and psychotherapeutic treatment including debriefing and stress coping, mindfulness based stress reduction/cognitive therapy (MBSR/MBCT), biofeedback, physical activation and relaxation techniques (e.g. progressive muscle relaxation, PMR).
Statistics. Continuous data were presented as median (quartile 1 -quartile 3) or, if appropriate, as mean and 95% confidence interval of the mean. Continuous data were given as counts and percentages. Differences in continuous data between groups were assessed by non-parametric Mann-Withney U tests, differences in nominal data by cross-tabulation and Pearson's χ² tests. Correlations between continuous variables were calculated according to Spearman and given as Spearman's ρ. The time courses of continuous variables were evaluated by two-way ANOVA with repeated measurements design (SPSS general linear models), giving readings of a continuous quantity (dependent variable) at two levels of a within-subject factor, and a dichotomous characteristic (e.g. group assignment) as independent, between-subjects factor. Interactions between the within-subject factor and the between-subject factor indicate that the within-subject factor develops different between both time points in the categories of the between-subject factor (e.g. results of a biomarker decline between baseline and follow-up examinations in group A, whereas they stagnate or even rise in group B).
Areas under the cortisol curve (AUC) have been calculated from 4 subsequent measurements (time of awakening, 15, 30, and 45 minutes after awakening) according to Pruessner et al. 44 , and given as AUC g and AUC i . In detail, AUC g (area under the curve with respect to ground) represents an approximated integral of the total saliva cortisol concentration over the first 45 minutes after awakening and is given in µg × h/dl: When subtracting the basal level measured at the time of awakening, integrated over the 0.75 hours of measurement, from the aforementioned AUC g , only the cortisol changes occurring due to the awakening response are taken into account. The resulting value is defined as AUC i (area under the curve with respect to increment): Notably, the term "cortisol awakening response (CAR)" is exclusively reserved for measurements of the post-awakening dynamic, e.g. AUC i 45 . Results being below the analytical lower limit of quantification (LLOD) were set to their arithmetical expected value (E (X) ) at LLOQ/2: Effect sizes were given with respect to the statistical tests applied: Mann-Whitney U test, r; Pearson's χ² test, ϕ; two-way ANOVA with repeated measurements design, partial η². Multiple testing correction of p-values was conducted according to Benjamini and Hochberg, and adjusted p-values were given as p BH . p-values < 0.05 were considered statistically significant. All calculations were made in SPSS 23.0 (IBM, Armonk, USA) and MedCalc Statistical Software (MedCalc Software bvba, Ostend, Belgium), figures were mainly drawn using GraphPad Prism (GraphPad Software Inc., La Jolla, USA). Data availability. Since publication of raw data was not consented by the study participants, data will not be openly available. However, raw data sets can be requested from the corresponding author by interested researchers.

Baseline cortisol levels differ between burnout patients and controls.
To validate the proposed stress biomarkers, baseline levels of blood-or saliva-borne molecules were measured in samples obtained from the 40 burnout patients, including cortisol, IL6, MPO, and HCYS. Marker levels were determined using commercially available immunoassays. Healthy volunteers (N = 26) were considered as controls. Results are summarized in Table 1.
Patients showed higher saliva cortisol concentrations for all awakening samples, with strongest effects for sample five and with the exception of sample three (30 minutes after awakening) and, accordingly, presented with a significantly higher AUC g value (p = 0.002, p BH < 0.01) when compared to the control group. However, differences in cortisol levels measured at later points in time (noon/evening) were even more pronounced, with patient levels failing to decrease during the day or at nighttime (both p/p BH < 0.001).
When comparing blood-borne biomarkers between stress patients and controls, we could not detect any significant differences. Though, IL6 showed trend-wise elevated levels within stress patients (p = 0.054, p BH = n.s.). Apart from that, we found an association between baseline IL6 and BMI among patients (ρ = 0.371, p = 0.019, p BH < 0.05), but not in controls (p = 0.237, p BH = n.s.). At follow up, this association was no longer detectable (ρ = 0.289, p = 0.070, p BH = n.s.; controls: p = 0.840, p BH = n.s.). Hence, the proposed parameters showed no discriminatory capabilities that could be used in identification of stress patients. Based on these data, it would appear that for the use of cortisol as a marker of burnout, midday and evening cortisol levels were at least as efficient as total salivary cortisol after awakening. In a next step, it was examined whether the investigated parameters might also allow for longitudinal monitoring of psychological stress in response to treatment.

Treatment of burnout patients corresponds to lower cortisol levels.
To determine whether cortisol measurement might be useful in monitoring therapy success in burnout patients, these were first assessed on a behavioral level. The aforementioned 40 burnout patients underwent a four month period of combined treatment, with a median of five (2-10) personal sessions, whereby the actual number of visits was very specifically oriented towards the individual patients' needs and requirements. By filling out appropriate psychosocial screening questionnaires (PSS, HADS anxiety and depression scores), they were evaluated by trained psychological personnel. Healthy volunteers, which did not receive any neuro-psychological treatment in this period of time, were assessed in analogy to patients. General linear models with repeated measurements design were used to arrive at statistical significance. Results are given in Fig. 3 and Table 2.
The results indicated that biopsychosocial treatment was beneficial to the burnout patients. Subsequently, it was evaluated whether this improvement would also be seen in reductions in biomarker levels over the treatment period. Therefore, the same patient group was asked to provide a saliva day profile and blood samples, once at the onset of therapy and once at its termination, and these were subjected to the aforementioned analyses. The control group was similarly assessed. Test statistics are presented in Table 3, temporal developments of biomarkers are visually depicted in Fig. 4.
Of the assessed stress associated biomarkers, only the temporal developments of midday cortisol and saliva cortisol nadir showed a significant statistical interaction with the group variable (although no longer statistically significant after p-value adjustment), indicating that both values decreased, as expected, in the burnout group, but not in controls. In detail, investigated patients presented with a midday cortisol decline between baseline and follow-up of −0.10 ± 0.14 µg/dl (p/p BH < 0.001) and also with a slightly less pronounced evening cortisol drop of −0.05 ± 0.14 (p = 0.022, p BH < 0.05). Interestingly, biomarker dynamics during the therapeutic phase highly depended on baseline levels, as higher baseline concentrations correlated with a less pronounced decrease at follow-up (see Fig. 5) for most parameters (all p/p BH < 0.001).
However, changes of biomarkers after the treatment phase did in neither case correlate with any improvement in one of the stress, mood or anxiety-related questionnaires (all p/p BH > 0.05). The significance of the utility of saliva cortisol, sampled conveniently at noon or at the evenings, to diagnose and monitor burnout is discussed.

Discussion
Burnout syndrome is gaining medical importance in the industrialized world, whereby people at their height of productivity might not be able to work any further, thereby meaning not only a major social-economic loss 46 , but additionally being exposed to serious suffering. Despite the ever increasing incidence of burnout 11 , reliable biomarkers to help with diagnosis or treatment monitoring are still missing 14 . Thus, the objective of this study was to determine whether cortisol day profiles, IL6, MPO, and HCYS plasma concentrations could be utilized as biomarkers in such patients with respect to clinical course and treatment outcome. To our knowledge, we showed here for the first time in a prospective setting that cortisol levels sampled at midday and in the evenings were not inferior to the more complex cortisol awakening response in distinguishing burnout patients and healthy controls.
In the present study, we measured saliva cortisol day profiles as well as plasma IL6, MPO and HCYS at baseline and in response to a four months treatment. At the clinical level, stress relief during the course of treatment was subjectively reported by anxiety, depression, and burnout-risk questionnaires, and objectively assessed by CGI ratings (Table 1). At the follow up visit the majority of patients (60%) showed a significant clinically improvement (indicated by a drop in the CGI score) that allowed patients to resume their work highlighting the socioeconomic value associated with such treatment programs. The results of this study demonstrated that the night time, midday (11:30-12:30) and, to a lesser extent, also morning cortisol concentrations were useful biomarkers for work-related stress symptoms when compared between patients and healthy controls, and within patients during recovery from major depression and anxiety disorders. Whereas blood-borne biomarkers did not exhibit any discriminatory or predictive usefulness, saliva cortisol levels presented with considerable diagnostic accuracy. Interestingly, cortisol samples taken either at midday or night time appeared to be superior to the more complex assessment of the cortisol awakening response. To our knowledge, this is the first study that prospectively monitored biomarkers of inflammation and oxidative stress during multiprofessional collaborative treatment 41 of work-related stress patients recruited beyond classical DSM-IV limitations according to RDoC principles 12 . The applied design was superior compared to cross-sectional studies regarding causal inference on longitudinal processes 47 . Hence, it is tempting to speculate that the observed symptom alleviation might be linked to the reduction of cortisol levels following therapy.
Various parameters of cortisol secretion have been applied in studies on psychiatric disorders. Regarding baseline levels, there are many studies reporting no difference in afternoon cortisol concentrations between controls and patients with current major depression or anxiety disorders 16 . By contrast, older adults suffering from  Table 2. Test statistics for general linear models (two-way ANOVA with repeated measurements design) assessing differences in relevant scores between baseline and follow-up examinations (main effects score), differences in neuropsychological scores between stress cohort and control group (main effect group), and, whether temporal developments of scores vary between groups (interaction group × score). Effect sizes are given as partial η². P-values given as p BH  depression were found to display significantly higher levels of basal cortisol than healthy controls during all phases of diurnal cycle but particularly during the evening and night hours 48 . Similar to these results, we have observed that, compared to controls, our burnout group had higher cortisol levels during the day, with the most significant value in the evening. Notably, degrees of depression, anxiety and stress were well reflected in AUC g values and even more pronounced in the nadir levels of cortisol. Treatment of the patients caused a decrease in midday and nadir cortisol only, whereas the concentrations of morning cortisol appeared to be unaffected. The cortisol awakening response is the most prominent component of the diurnal cortisol profile investigated in psychosocial and mental health studies 38 . Deviations from a regular CAR profile are assumed to indicate maladaptive neuroendocrine processes. Meta-analysis of the AUC i showed positive associations between the magnitude of the AUC i and occupational stress and stressful life-events 49 . By contrast, there was a tendency for negative correlations between AUC i and fatigue, exhaustion, or burnout, the latter representing a symptom cluster with quite inconsistent results [50][51][52][53][54] . However, the present study neither identified an altered level of AUC i in the burnout group nor any statistically significant association between AUC i and symptom scores. This is in line with previous studies showing no relationship between burnout questionnaire and CAR 51 , and no impact of depression, hopelessness and negative affect on CAR 49 . The fact that we found equal AUC i values in patients and controls whilst the burnout group exhibited an overall elevation of the cortisol profile may be explained by the well-known uncoupling of CAR from basal circadian cortisol secretion 38 . A similar effect has been reported when assessing differences in CAR in relation to severity of burnout 52 . Patients with low burnout displayed lower basal levels of cortisol than patients with moderate burnout whereas the relative increase in cortisol over the first 30 min after morning awakening appeared to be unchanged. Taken together, these results point to some limitation of morning AUC i for comparisons of cortisol patterns in psychiatric disorders.
In our study, midday and nadir cortisol showed the most pronounced correlations with HADS anxiety score. In addition, a significant correlation between AUC g and anxiety score was found, whereas no association between AUC i and degree of anxiety could be observed. A possible explanation could be that the more complex assessment of morning response curves produces a higher pre-analytical variability 38 . Indeed, measures of dispersion are slightly higher in AUC g (Standard error of the mean, SEM = 0.024) when compared to midday (SEM = 0.020) or nadir cortisol (SEM = 0.014). Above that, cortisol hypersecretion is clearly associated with depressive behavior 55 . Especially rumination, which is an aggravating symptom of burnout and depression 56 , was shown to delay the cortisol decline during the day 57 . Since especially nadir concentrations of healthy individuals were hardly measurable and presented with very low variability, the difference to burnout patients might become more obvious. Moreover, our results support the findings of Merswolken et al. 58 , who demonstrated on the one hand a  significant correlation between anxiety and AUC g , and on the other hand only a trend towards an association between anxiety and AUC i in CHD patients. A recent report by Hakamata et al. 59 provides a possible mechanistic model for the relationship between cortisol and anxiety by underlining the importance of the strength of amygdala-hippocampus functional connectivity. Estimated marginal means are calculated by general linear models (two-way ANOVA with repeated measurements design) and p-values were adjusted according to Benjamini and Hochberg. Data is given as estimated marginal mean and 95% confidence interval, and confidence intervals from stress patients are light grey shaded, whereas confidence intervals from controls are dark grey shaded. (b) Variation in cortisol day profiles at baseline (light grey) and before follow-up examinations (dark grey) within the stress (top) and the control cohort (bottom). AUC g , area under the curve with respect to ground; AUC i , cortisol awakening response -area under the curve with respect to increase; MPO, plasma myeloperoxidase; IL6, interleukin 6; n.s., not significant; ***p BH < 0.001; *p BH < 0.05. IL6, MPO, HCYS plasma concentrations did not change within subjects during recovery from stress, and did not differ between patients and controls. Even though IL6 was trend-wise higher in the patient cohort, these measures might, according to our findings, not be useful as biomarkers to monitor stress relief in patients.
This study comes with two limitations, which, however, did not markedly affect quality: First, the sample size appears prima facie to be relatively small compared to the number of covariates. However, correlations between changes in biomarkers and psychosocial variables presented with very low effect sizes (ρ < 0.25), which might probably not be diagnostically relevant. Therefore, the statistical detection of those relationships would have no clinical value. The differences in sample size between control and intervention group can be considered irrelevant as well, since relevant parameters (e.g. age, sex) are equally distributed in both groups. Taking this together, our sample size can be considered as sufficient for testing clinically relevant interdependencies. Indeed, inspired by previous clinical studies 60 , we preferred a naturalistic and population-based setting to ensure realistic conditions with higher external validity. The second limitation concerns saliva pre-analytics. Of course, subjects were informed about the necessity of accurate sample collection at home. However, temporal accuracy cannot be completely assured, since collections were not monitored by study personnel. To enhance compliance, subjects were required to document the exact time points in a day protocol.
In conclusion, the present results encompass that high levels of stress, depression, and anxiety are reflected by increased concentrations of salivary cortisol, especially just before going to bed at night (nadir). Although patients responded on a clinical level, neither of the assessed blood parameters showed to be a predictor of clinical treatment outcome. Hence, our study suggests that measurement of salivary cortisol might be a useful objective adjunct in addition to clinical diagnostics, and, to a lesser extent, as objective marker for disease course. Notably, single cortisol values from saliva collected at noon or bedtime exhibited considerably better diagnostic quality than more expensive and error-prone sequences, as e.g. the cortisol awakening response, a variable which has led to inconsistent results in the biomonitoring of burnout 14 .