key: cord-0776195-3499ol9u
authors: Karampitsakos, Theodoros; Torrisi, Sebastiano; Antoniou, Katerina; Manali, Effrosyni; Korbila, Ioanna; Papaioannou, Ourania; Sampsonas, Fotios; Katsaras, Matthaios; Vasarmidi, Eirini; Papakosta, Despoina; Domvri, Kalliopi; Fouka, Eva; Organtzis, Ioannis; Daniil, Zoe; Dimeas, Ilias; Kirgou, Paraskevi; Gourgoulianis, Konstantinos I.; Papanikolaou, Ilias C.; Markopoulou, Katerina; Kounti, Georgia; Tsapakidou, Eirini; Papadopoulou, Efthymia; Tatsis, Konstantinos; Gogali, Athena; Kostikas, Konstantinos; Tzilas, Vasilios; Chrysikos, Serafeim; Papiris, Spyridon; Bouros, Demosthenes; Kreuter, Michael; Tzouvelekis, Argyrios
title: Increased monocyte count and red cell distribution width as prognostic biomarkers in patients with Idiopathic Pulmonary Fibrosis
date: 2021-05-05
journal: Respir Res
DOI: 10.1186/s12931-021-01725-9
sha: bbce8aba8d9be8b113c6194c3ced40ff236de455
doc_id: 776195
cord_uid: 3499ol9u

BACKGROUND: Idiopathic Pulmonary Fibrosis (IPF) represents a chronic lung disease with unpredictable course. METHODS: We aimed to investigate prognostic performance of complete blood count parameters in IPF. Treatment-naïve patients with IPF were retrospectively enrolled from two independent cohorts (derivation and validation) and split into subgroups (high and low) based on median baseline monocyte count and red cell distribution width (RDW). RESULTS: Overall, 489 patients (derivation cohort: 300, validation cohort: 189) were analyzed. In the derivation cohort, patients with monocyte count ≥ 0.60 K/μL had significantly lower median FVC%pred [75.0, (95% CI 71.3–76.7) vs. 80.9, (95% CI 77.5–83.1), (P = 0.01)] and DLCO%pred [47.5, (95% CI 44.3–52.3) vs. 53.0, (95% CI 48.0–56.7), (P = 0.02)] than patients with monocyte count < 0.60 K/μL. Patients with RDW ≥ 14.1% had significantly lower median FVC%pred [75.5, (95% CI 71.2–79.2) vs. 78.3, (95% CI 76.0–81.0), (P = 0.04)] and DLCO%pred [45.4, (95% CI 43.3–50.5) vs. 53.0, (95% CI 50.8–56.8), (P = 0.008)] than patients with RDW < 14.1%. Cut-off thresholds from the derivation cohort were applied to the validation cohort with similar discriminatory value, as indicated by significant differences in median DLCO%pred between patients with high vs. low monocyte count [37.8, (95% CI 35.5–41.1) vs. 45.5, (95% CI 41.9–49.4), (P < 0.001)] and RDW [37.9, (95% CI 33.4–40.7) vs. 44.4, (95% CI 41.5–48.9), (P < 0.001)]. Patients with high monocyte count and RDW of the validation cohort exhibited a trend towards lower median FVC%pred (P = 0.09) and significantly lower median FVC%pred (P = 0.001), respectively. Kaplan–Meier analysis in the derivation cohort demonstrated higher all-cause mortality in patients with high (≥ 0.60 K/μL) vs. low monocyte count (< 0.60 K/μL) [HR 2.05, (95% CI 1.19–3.53), (P = 0.01)]. CONCLUSIONS: Increased monocyte count and RDW may represent negative prognostic biomarkers in patients with IPF. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12931-021-01725-9.

Idiopathic pulmonary fibrosis (IPF) represents a chronic, progressive lung disease with dismal prognosis despite the advent of novel antifibrotic compounds [1] [2] [3] [4] . The disease course is highly variable and prognosis remains challenging [5] . A clinicians' friendly, easily applicable and cost-effective prognostic biomarker with uniform cut-off values that will guide disease stratification and tailoring of therapeutic approaches is missing [6] . Several biomarkers including Mucin 5b and Toll-interacting protein single nucleotide polymorphisms, telomere length and gene expression signatures have been suggested as reliable prognosticators in patients with IPF [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] ; yet, measurement of such biomarkers remains laborious and requires expensive and sophisticated infrastructure. In addition, lack of standardization of samples' collection protocols leading to non-reproducible cut-off thresholds, further limits their widespread clinical applicability [6, 19] .

Abundant evidence has highlighted the role of Complete Blood Count (CBC) in the prognostication of patients with various chronic lung diseases [6, 20, 21] . Asthma researchers are currently using peripheral eosinophils to implement anti-IL5/13 therapeutic regimens [22] . Three major studies, encompassing an overall of almost 10,000 patients with IPF and scleroderma-associated interstitial lung disease have recently identified elevated peripheral blood monocyte count as a biomarker of disease progression and mortality [6, 23] . Given that monocyte count is a clinically applicable and inexpensive biomarker, these findings warrant further investigation in real-life studies. In line with this concept another parameter of CBC, red cell distribution width (RDW), has been associated with worse clinical outcomes in several chronic lung diseases including IPF and chronic obstructive pulmonary disease [20, [24] [25] [26] . Increased RDW seems to represent a biomarker of early hypoxemia [20, 27, 28] .

To this end, our aim was to evaluate the prognostic role of parameters of CBC, including monocyte count and RDW in two independent cohorts (derivation and validation) of patients with IPF in a real-life clinical setting.

This was an observational, retrospective study. Between 01/11/2018 and 31/08/2020, we retrospectively enrolled patients with IPF and available CBC at baseline (prior to anti-fibrotic treatment), as well as 6 and 12 months posttreatment. Only treatment-naïve patients at the time of baseline CBC were included in the analysis. There were no patients receiving antifibrotics or corticosteroids at the time of baseline CBC. Epidemiological data were derived from two independent cohorts.

Derivation cohort The derivation cohort included patients from referral centers for Interstitial Lung Diseases in Greece including Department of Internal and Respiratory Medicine, University Hospital of Patras, 1st and 2nd Academic Department of Respiratory Medicine, "SOTIRIA" and "ATTIKON" General Hospital, National and Kapodistrian University of Athens, Laboratory of Molecular and Cellular Pneumonology, Department of Respiratory Medicine, Faculty of Medicine, University of Crete, Heraklion, Crete, Medical School, University of Thessaly, Larissa, Respiratory Medicine Department, "Corfu General Hospital", Department of Respiratory Medicine, "G. PAPANIKOLAOU" General Hospital, Thessaloniki, Aristotle University of Thessaloniki and Department of Respiratory Medicine, Medical School, University of Ioannina.

Validation cohort The validation cohort included patients from the Center for interstitial and rare lung diseases, Pneumology, Thoraxklinik, University of Heidelberg, Germany and German Center for Lung Research, Heidelberg, Germany.

The study was approved by the Institutional Review Board and the Local Ethics Committee (Protocol Number: 458/06- [12] [13] [14] [15] [16] [17] [18] [19] . Diagnosis of IPF was based on ATS/ ERS/JRS/ALAT guidelines [1] . We collected parameters of CBC including monocyte count and RDW prior antifibrotic treatment, as well as 6 and 12 months post-treatment. Baseline demographics and comorbid conditions were recorded. Pulmonary hypertension was defined as elevated right ventricular systolic pressure on echocardiographic assessment, as all patients performed a baseline echocardiography but not right heart catheterization.

Median values of CBC parameters were recorded. Median values were used, as Kolmogorov-Smirnov test for normal distribution rejected normality. Patients were divided in subgroups based on the median value of each CBC parameter in the derivation cohort (high and low). We used median values based on the fact that monocyte count and RDW did not have significant differences over Conclusions: Increased monocyte count and RDW may represent negative prognostic biomarkers in patients with IPF.

Keywords: Monocyte count, RDW, Idiopathic pulmonary fibrosis, Biomarkers, Mortality the 1-year period both in our cohorts and in other studies [25, 29] . Mann-Whitney test was applied to assess differences in Forced Vital Capacity %predicted (FVC%pred) and Diffusion capacity of lung for carbon monoxide %predicted (DLCO%pred) between subgroups of patients split by the median value of CBC parameters. Prognostic performance of these cut-off thresholds was also assessed in the validation cohort. Kaplan-Meier survival analysis was applied to investigate differences in survival probability between high and low subgroups. Kaplan-Meier was also used to dichotomize patients based on the previously published cut-off threshold of monocyte count (0.95 K/μL) [6] . Differences in parameters of CBC between patients in need of Long Term Oxygen Therapy (LTOT) and patients without LTOT were investigated with Mann-Whitney. P -values < 0.05 were considered statistically significant.

Overall, 489 patients (derivation cohort: N = 300, validation cohort: N = 189) were included in the analysis. Patient demographics and disease characteristics are summarized in Table 1 . Median age (95% CI) was 74 (73-75) years for the derivation cohort and 74 (72-75) years for the validation cohort. Patients with IPF were predominantly male both in the derivation (83.3%, N = 250) and validation cohort (78.8%, N = 149). Median monocyte count (95% CI) was 0.60 (0.57-0.62) and 0.52 (0.50-0.58) K/μL for the derivation and the validation cohort, respectively. Median RDW (95% CI) for the derivation cohort was 14.1% (13.9-14. 3) and 13.7% (13.6-13.8) for the validation cohort. Finally, median FVC%pred (95% CI) was 77.0 (75.0-79.8) and 76.2 (71.7-80.8), while median DLCO%pred (95% CI) was 51.0 (47.1-53.8) and 41.9 (40.3-44.9) for the derivation and validation cohort, respectively (Table 1) . Median follow-up (95% CI) was 24.3 (23.4-28.7) and 15.0 (12.0-19.0) months for the derivation and the validation cohort, respectively. (Fig. 1a, c) . In the validation cohort, a trend towards lower FVC%pred [median: 72.4, (95% CI 68.8-79.6)] in the high monocyte group (≥ 0.60 K/μL) compared to the low monocyte group [median: 79.5, (95% CI 72.9-82.9)] was observed; yet, no statistical significance was reached (P = 0.09), (Fig. 1b) . Importantly, patients in the high monocyte group of the Multiple regression analysis of the overall study population showed that baseline monocyte count was independently associated with baseline DLCO%pred (P = 0.005). With regards to other factors, gender was independently associated with baseline FVC%pred (P = 0.007), while presence of pulmonary hypertension was independently associated with baseline DLCO%pred (P = 0.002) ( Table 2) . Moreover, multiple regression analysis investigating the impact of comorbidities on baseline CBC values, showed that presence of pulmonary hypertension was independently associated with baseline RDW (P < 0.001) (Additional file 1: Table S1 ). (Table 3) .

There was no statistically significant difference in 1-year FVC%pred and DLCO%pred decline between patients with high and low monocyte count [median ΔFVC%pred, derivation cohort: 0.0 (95% CI − 2.8- (Table 4 ). There was no statistically Table S3 ).

In descriptive analysis, median monocyte count was similar over 1 Table S4 ).

In descriptive analysis, median RDW was similar over Table S5 ).

In the derivation cohort, patients in the high monocyte group (≥ 0.60 K/μL) experienced increased risk of all-cause mortality compared to the low group (monocyte count < 0.60 K/μL) [HR 2.05, (95% CI 1.19-3.53), (P = 0.01)] (Fig. 3a) . This finding was not confirmed in the validation cohort (P = 0.79) (Fig. 3b) . Moreover, pooled analysis of the study population demonstrated an increased risk in all-cause mortality for patients with baseline monocyte count ≥ 0.95 K/μL compared to patients with baseline monocyte count < 0.95 Κ/μL [HR 2.47, (95% CI 0.94-6.47), (P = 0.005)] (Additional file 1: Figure S1 ). There was no increased risk of allcause mortality for patients in the high RDW group compared to the low group, in both the derivation (P = 0.82) and the validation (P = 0.90) cohort, as well as in the pooled analysis (P = 0.41). Data for mortality were available for 218, 155 and 373 patients in the derivation cohort, validation cohort and pooled analysis, respectively.

This real-life retrospective study demonstrated that peripheral blood monocyte count was predictive of allcause mortality in the derivation cohort and in a pooled collective of highly characterized patients with IPF. We also showed that patients with elevated levels of monocyte count and RDW exhibited more advanced disease at initial assessment compared to patients with low levels. There was no association of high monocyte count or RDW with 1-year disease progression, as assessed by functional decline. No effects of anti-fibrotic treatment on monocyte count or RDW were observed over 1-year of follow-up. Differences in baseline monocyte count, RDW, DLCO% pred and LTOT use between the two cohorts might be partially attributed to divergent endotypes across the world and/or different baseline functional status.

Our findings are consistent with those of previous reports evaluating a possible link between monocyte count and prognosis in patients with IPF [6, 25, 30, 31] . A previous retrospective, multicenter cohort study showed that monocyte count ≥ 0.95 K/μL was significantly associated with all-cause mortality compared to monocyte count < 0.95 K/μL in 7459 patients with IPF [6] . Nonetheless, IPF diagnosis in this study was based on ICD-10 medical records posing limitations to the findings. Analysis of 231 patients with IPF from the Australian registry corroborated evidence that elevated monocyte count were associated with worse clinical outcomes [30] . Most recently, pooled retrospective analysis of 2067 highly characterized patients with IPF derived from the pirfenidone trials (ASCEND, CAPACITY and INSPIRE) showed that patients with IPF and monocyte count in the range of 0.60-0.95 K/μL or ≥ 0.95 K/μL had a higher 1-year risk of IPF progression, all-cause hospitalization and all-cause mortality compared to patients with monocyte count of < 0.60 K/μL [29] . Given the results from pirfenidone clinical trials and our real-life study, monocyte count ≥ 0.60 K/μL, appears to be a highly robust and reproducible cut-off threshold which could potentially enrich the population of clinical trials, as a marker associated with greater risk of mortality and/or disease progression. In addition, it may alert clinicians in the context of risk stratification for timely interventions. In our study, monocyte count was predictive of all-cause mortality in the derivation but not the validation cohort. This might be partially attributed to the worse baseline functional status of the validation cohort, as indicated by the lower DLCO%pred and increased use of LTOT at baseline. With regards to RDW, a previous study enrolling 319 patients with IPF reported lower median DLCO%pred and increased mortality risk for patients with RDW > 15% compared to patients with RDW ≤ 15% [25] . Our study yielded similar results for DLCO%pred. Subgroup analysis of our cohorts did not show a survival benefit for patients with RDW < 14.1%; yet, our study was designed to assess differences in subgroups based on the median RDW (14.1%) and not based on the previously published cut-off threshold of 15%.

In our study, patients with increased baseline monocyte count and RDW exhibited more advanced disease at initial assessment as indicated by baseline FVC%pred and DLCO%pred. However, monocyte count and RDW were not associated with 1-year FVC%pred and DLCO%pred decline. Similarly to our findings, recent evidence using pooled data from the TOMORROW and INPULSIS trials, showed that the adjusted rate of FVC decline was similar between patients with high and low monocyte count receiving nintedanib [32] . Nonetheless, there is still a major knowledge gap for the longitudinal prognostic and theragnostic role of these biomarkers. The prognostic role of monocyte count in FVC decline requires further investigation, as its prognostic accuracy might be associated with the baseline status, the selected treatment or the population investigated [29, 32] . Further large studies are needed to address this issue, focusing on subgroup of patients that have been subjected to different treatment modalities.

Importantly, monocyte count and RDW were similar over 1-year follow-up of antifibrotic treatment either with pirfenidone or nintedanib. RDW has been widely considered a reproducible marker, given the relatively prolonged lifespan of red blood cells [27] . Previous reports have shown that patients with a high monocyte count at diagnosis maintained their high count through the disease course [6, 17] . There was no correlation between change in monocyte count over time and survival [6] . Instead, monocyte count seemed to be relatively stable over time indicating that patients with IPF retained the same risk profile [6, 17] . To this end, monocyte count seems to have greater potential as a prognostic biomarker rather than as a predictive biomarker of treatment response; nonetheless, this requires further investigation in future cohorts applying subgroup analyses.

The prognostic role of monocyte count and RDW in patients with IPF could be partially explained from recently emerged experimental evidence suggesting migration of monocytes from the bone marrow to the injured lung and differentiation to pro-fibrotic macrophages or even fibroblasts [33] [34] [35] [36] [37] . Single-cell RNA sequencing characterized the heterogeneity of macrophages in bleomycin-induced pulmonary fibrosis and identified a pathological subgroup of transitional macrophages (CX3CR1 + SiglecF +) required for the fibrotic response to the injurious stimuli [38] . Recent evidence has shown that the compartmental imbalance of fractalkine mediated the migration of CX3CR1 + nonclassical monocytes into fibrotic lung tissues, while non-classical monocytes-derived cells presented with a M2-like and phagocytic phenotype in fibrotic lungs [39] . Translational studies have demonstrated that accumulation of distinct populations of alveolar macrophages and higher levels of circulating fibrocytes, derived from the monocyte cell lineage, may be predictive of pulmonary fibrosis progression [36, 37, 40, 41] . Finally, C-C motif chemokine ligand 18, produced in a considerable extent by alveolar macrophages, has been suggested as a promising, serum biomarker of disease progression and mortality in patients with IPF [42] .

On the other hand, a causal-effect relationship between IPF and elevated RDW is highly unlikely. Instead, it is more likely that increased RDW is indicative of patients' hypoxemia and/or comorbidome in a similar way with other chronic lung diseases [20, [24] [25] [26] [27] . It has been proposed that arterial hypoxemia leads to increased erythropoietin secretion and thus to increased RDW through mechanisms involving regulation of erythrocyte maturation and survival [20, 27] . Elevated RDW might have a role in the early identification of patients with IPF and intermittent hypoxemia [25] . Patients with profound hypoxemia are easily diagnosed; intermittent hypoxemia might escape routine examination and there is still a need for biomarkers contributing to their identification.

Our study has some limitations. First of all, our study has the inherent weaknesses of a retrospective study. Nonetheless, the nature of this study enabled us to report longitudinal outcomes of patients with IPF. Secondly, our sample size is moderate compared to previous reports for the prognostic role of monocyte count; yet, the size is acceptable for a real-life study. Thirdly, we had data for LTOT, but not for po2 levels; thus, we could not further investigate the association of hypoxemia with monocyte count and RDW. Moreover, data for all-cause mortality was available; yet, the specific cause of death was not available for all patients. Finally, our results should be interpreted in the context of a real-life study that may be in part influenced by other factors including steroid use prior to admission at a referral center. To this end, multiple regression analysis was performed to adjust for these covariates.

This was the first real-life study of highly characterized patients with IPF showing that patients with IPF and high monocyte count (≥ 0.60 K/μL) exhibited more advanced disease at initial assessment and had a higher risk of allcause mortality compared to patients with low monocyte count (< 0.60 K/μL). RDW failed to predict disease progression and all-cause mortality. Our study coupled with previous reports demonstrating that peripheral blood monocytes can be easily incorporated into the routine clinical assessment of patients with IPF as a reliable prognostic biomarker considering its reproducibility, costeffectiveness and simplicity. Future prospective studies investigating the association of baseline and serial measurements of monocyte count with disease outcomes and treatment response are greatly anticipated.

Abbreviations CBC: Complete blood count; DLCO%pred: Diffusion capacity of lung for carbon monoxide %predicted; FVC%pred: Forced vital capacity %predicted; IPF: Idiopathic pulmonary fibrosis; LTOT: Long term oxygen therapy; RDW: Red cell distribution width.

Diagnosis of idiopathic pulmonary fibrosis. An official ATS/ ERS/JRS/ALAT clinical practice guideline

A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis

Pirfenidone in patients with idiopathic pulmonary fibrosis (CAPAC-ITY): two randomised trials

Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis

Clinical course and prediction of survival in idiopathic pulmonary fibrosis

Increased monocyte count as a cellular biomarker for poor outcomes in fibrotic diseases: a retrospective, multicentre cohort study

Prognostic and predictive biomarkers for patients with idiopathic pulmonary fibrosis treated with pirfenidone: post-hoc assessment of the CAPACITY and ASCEND trials

Validation of the prognostic value of MMP-7 in idiopathic pulmonary fibrosis

An epithelial biomarker signature for idiopathic pulmonary fibrosis: an analysis from the multicentre PROFILE cohort study

Longitudinal change in collagen degradation biomarkers in idiopathic pulmonary fibrosis: an analysis from the prospective, multicentre PROFILE study

Idiopathic pulmonary fibrosis: prospective, case-controlled study of natural history and circulating biomarkers

Effect of telomere length on survival in patients with idiopathic pulmonary • fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year submit your research ? Choose BMC and benefit from: ? Choose BMC and benefit from: fibrosis: an observational cohort study with independent validation

Peripheral blood proteins predict mortality in idiopathic pulmonary fibrosis

Biomarkers of collagen synthesis predict progression in the PROFILE idiopathic pulmonary fibrosis cohort

Association between the MUC5B promoter polymorphism and survival in patients with idiopathic pulmonary fibrosis

MUC5B, and the response to N-acetylcysteine among individuals with idiopathic pulmonary fibrosis

Validation of a 52-gene risk profile for outcome prediction in patients with idiopathic pulmonary fibrosis: an international, multicentre, cohort study

Biomarkers in the evaluation and management of idiopathic pulmonary fibrosis

Can monocytes predict prognosis of idiopathic pulmonary fibrosis?

The role of increased red cell distribution width as a negative prognostic marker in patients with COPD

Comparative effectiveness of LABA-ICS versus LAMA as initial treatment in COPD targeted by blood eosinophils: a population-based cohort study

Asthma biologics: comparing trial designs, patient cohorts and study results

Blood monocyte counts as a potential prognostic marker for idiopathic pulmonary fibrosis: analysis from the Australian IPF registry

Usefulness of red cell distribution width as a prognostic marker in pulmonary hypertension

The red cell distribution width as a prognostic indicator in idiopathic pulmonary fibrosis

Elevated red cell distribution width predicts poor outcome in young patients with community acquired pneumonia

Increased red cell distribution width: a novel predictor of adverse outcome in patients hospitalized due to acute exacerbation of chronic obstructive pulmonary disease

Increased red cell distribution width is associated with disease severity in hospitalized adults with SARS-CoV-2 infection: an observational multicentric study

Monocyte count as a prognostic biomarker in patients with idiopathic pulmonary fibrosis

Blood monocyte counts as a potential prognostic marker for IPF: analysis from the Australian IPF registry

Monocyte count and the risk for acute exacerbation of fibrosing interstitial lung disease: a retrospective cohort study

Monocyte count and decline in forced vital capacity (FVC) in patients with IPF

Idiopathic pulmonary fibrosis

A systematic review of the role of dysfunctional wound healing in the pathogenesis and treatment of idiopathic pulmonary fibrosis

Thyroid hormone inhibits lung fibrosis in mice by improving epithelial mitochondrial function

The transferrin receptor CD71 delineates functionally distinct airway macrophage subsets during idiopathic pulmonary fibrosis

Fibrocytes are increased in lung and peripheral blood of patients with idiopathic pulmonary fibrosis

Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage

CX3CR1-fractalkine axis drives kinetic changes of monocytes in fibrotic interstitial lung diseases

Elevation of pulmonary CD163(+) and CD204(+) macrophages is associated with the clinical course of idiopathic pulmonary fibrosis patients

Circulating fibrocytes are an indicator of poor prognosis in idiopathic pulmonary fibrosis

Serum CC-chemokine ligand 18 concentration predicts outcome in idiopathic pulmonary fibrosis

Publisher's Note

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

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12931-021-01725-9.Additional file1: Additional Tables and Figure. Authors' contributions Initial study conception and design: TK, DB, MK (Michael Kreuter), AT. Substantial contributions to the conception or design of the work or the acquisition, analysis, or interpretation of data for the work: All authors. Drafting initial version of the manuscript: TK, MK (Michael Kreuter), AT. Drafting the work or revising it critically for important intellectual content: All authors. Final approval of the version to be published: All authors. Agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved: All authors. Main guarantors of the paper: TK and AT. All authors read and approved the final manuscript.

No funding.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

The study was approved by the Institutional Review Board and the Local Ethics Committee (Protocol Number: 458/06-12-19).

All authors consent.