{'fraction_non_alphanumeric': 0.03145674300254453, 'fraction_numerical': 0.01612595419847328, 'mean_word_length': 5.018568147013783, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Physical activity and Mediterranean diet as potential modulators of osteoprotegerin and soluble RANKL in gBRCA1/2 mutation carriers: results of the lifestyle intervention pilot study LIBRE-1

Purpose Emerging evidence suggests that the progesterone-mediated receptor activator of nuclear factor κB (RANK)/soluble RANK ligand (sRANKL)/osteoprotegerin (OPG) pathway plays an important role in mammary carcinogenesis and is hyperactivated in germline (g)BRCA1/2 mutation carriers. We analyzed the effects of a 3-month intensive lifestyle intervention within the LIBRE-1 study on the serum levels of OPG and sRANKL and hypothesized that the intervention program provides a beneficial impact on the biomarkers by increasing OPG and reducing sRANKL serum concentrations. Methods Serum levels of OPG and sRANKL of 49 gBRCA1/2 mutation carriers were quantified using enzyme-linked immunosorbent assays. We used previously collected blood samples from participants of the prospective LIBRE-1 study, who were randomized into an intervention group (IG), increasing physical activity and adherence to the Mediterranean diet (MedD) through supervised sessions from study entry to the first study visit after 3 months and a usual-care control group (CG). Differences in biomarker levels before and after the 3-month intervention were tested within and between study groups. Results The lifestyle intervention resulted in a significant increase in OPG for participants in both the IG (q = 0.022) and CG (q = 0.002). sRANKL decreased significantly in the IG (q = 0.0464) and seemed to decrease in the CG (q = 0.5584). An increase in the intake of Omega-3 polyunsaturated fatty acids was significantly associated with an increase in OPG (r = 0.579, q = 0.045). Baseline serum levels of sRANKL were a strong predictor for the change of sRANKL in the course of the intervention (ß-estimate = − 0.70; q = 0.0018). Baseline physical fitness (assessed as VO2peak) might predict the change of OPG in the course of the intervention program (ß-estimate = 0.133 pg/ml/ml/min/kg; p = 0.0319; q = 0.2871). Conclusion Findings from this pilot study seem to confirm our hypothesis by showing an increase in OPG and decrease in sRANKL over a 3-month lifestyle intervention and suggest that increased physical activity and adherence to the MedD are potent modulators of the biomarkers OPG and potentially sRANKL. Supplementary Information The online version contains supplementary material available at 10.1007/s10549-021-06400-7.

Study visit 3 months after study entry in LIBRE-1 VO 2 peak Maximal oxygen uptake in ml/min/kg Background Women with a germline mutation in the tumor suppressor genes BRCA1/2 have a high lifetime risk of developing breast cancer (BC) or ovarian cancer (OC) (69-72% and 16-59%, respectively) [1]. Since the penetrance of cancer disease is high but incomplete, the existence of risk-modulating factors has been postulated, some of which could be lifestyleassociated [2,3]. While scientific research has demonstrated that lifestyle factors, such as physical activity and nutrition, are modifiers of risk for sporadic BC [4][5][6], the impact on genetically predisposed women, such as gBRCA1/2 mutation carriers, is still unclear. At present, prophylactic mastectomy and/or salpingoovarectomy are the only available, but rather drastic primary prevention options for these women, presenting an urgent need for developing less invasive strategies [7].
These abnormal molecular patterns in gBRCA1/2 mutation carriers seem highly relevant in developing disease prevention strategies, as physical activity has been found to affect serum levels of OPG and possibly sRANKL. Studies in the general population suggest that endurance training increases circulating levels of OPG [27,28] and might diminish serum levels of progesterone [29,30] and RANKL [27].
Additionally, there is a growing body of evidence that the Mediterranean diet (MedD) might modulate these biomarkers [31,32]. The MedD describes a nutritional concept that is based on vegetables, grains, and nuts and is characterized by a proportionally high intake of olive oil and fish, which are rich in the n-3 and n-9 polyunsaturated fatty acids (PUFA), as well as a moderate intake of dairy products and low consumption of red meat, which contain primarily n-6 PUFA [33]. Several studies have shown that adherence to the MedD can reduce the risk of cardiovascular disease [34], inflammation [35], as well as mortality due to BC [36].
Using data from the LIBRE-1 study cohort, we aimed to explore the effects of a structured lifestyle intervention program on biomarkers that seem to be dysregulated in gBRCA1/2 mutation carriers [37,38].

Study cohort and participant recruitment
Women aged 18-69 years with a confirmed germline mutation in the BRCA1 or BRCA2 genes were eligible to participate in the LIBRE-1 study, a multicenter, prospective, two-armed randomized (1:1) controlled lifestyle intervention study. The LIBRE-1 study started in February 2014 [37] and was planned as the pilot study of the subsequently performed larger LIBRE-2 study, which started in 2015 [38]. The LIBRE-2 study is currently still recruiting (n = 484, as of May 2021). The primary endpoint of the LIBRE-1 study was to investigate the feasibility of a structured lifestyle intervention and its acceptance by the participants, which were both confirmed positive. [37,39]. In the present secondary analysis of the LIBRE-1 cohort, we analyzed OPG and sRANKL serum levels from all participants of the LIBRE-1 study (total n = 68) of whom both blood samples and all clinical data, including questionnaires, anthropometry, cardiopulmonary exercise testing (CPET), and measures on fatty acid profiles in red blood cell membranes (RBCM), were available for both time points of study entry (SE) as well as 3 months later (V1), resulting in 49 individuals (see Fig. 1). We focused on the 3-month period of the LIBRE-1 intervention as more pronounced lifestyle-associated effects on the biomarkers were expected during the intensive phase of the program [40]. Also, because of the small sample size and limited detection capacities of our quantitative tests, smaller changes in OPG and sRANKL serum levels might be more difficult to detect. Both gBRCA1/2 mutation carriers with a history of BC or OC prior to SE (diseased participants) and non-diseased participants, without previous BC or OC, were included in the LIBRE-1 study. Informed written consent was obtained from all participants. The study was approved by the ethics committees of the three participating centers (

Intervention program
The intervention included a structured endurance training and nutrition education based on the MedD. Over 3 months, the IG received two supervised and one homebased physical training session per week with the objective of increasing physical activity to at least 18 metabolic equivalents of task/-hour, mainly through endurance-based exercise to improve cardiopulmonary fitness. Nutrition education followed a group-based approach on fortnightly basis. The CG attended only one information session at SE on the benefits of physical activity and nutritional recommendations in accordance with the German Nutrition Society [37].

Assessment of clinical data and questionnaires
Data were collected by means of clinical interviews and questionnaires, CPET, and blood sampling during the study visits at SE and at 3 months (V1) and secured by each study center in a central database using OpenClinica (Waltham, MA, USA). To gain objective data on the effects of the intervention program on cardiopulmonary fitness, the maximal oxygen uptake (VO2 peak in ml/min/ kg) as a validated surrogate parameter for endurance performance [41] was measured using CPET. CPET describes a ramp protocol, which measures respiratory gases under gradually increasing resistance to achieve physical exhaustion [42]. Body Mass Index (BMI, kg/m 2 ) was determined during physical examination. The Mediterranean Diet Adherence Screener (MEDAS) [43] was used to score the nutritional adherence to the MedD with 14 items on food consumption and MedD habits. To reduce the bias of unanswered questions, the score was calculated as a percentage of positively answered questions (1 = stronger adherence to MedD, 0 = less adherence to MedD) to all answered questions [44]. Participants fasted a minimum of 12 h prior to each study visit. Blood samples were processed to aliquots and either analyzed within 48 h or stored at − 80 °C until usage. Furthermore, we used data of PUFA in the RBCM as well as in plasma that were conducted on the same LIBRE-1 blood samples, as described elsewhere [44]. ), we used enzyme-linked immunosorbent assay (ELISA) kits from Immundiagnostik AG (Bensheim, Germany). All tests were conducted according to the manufacturer's protocol. All samples were analyzed in double measurement on the same microtiter plate. As provided by the manufacturers, 6 standards and 2 control samples for sRANKL and 6 standards and 1 control sample for OPG were run within each kit. In total we used 6 kits for sRANKL and 4 kits for OPG. For quality control, we assessed the intra-and inter-assay coefficients of variation (CV) for each kit, being between 5.6 and 6.5% (inter-assay CV), as well as between 3.5 and 12.6% (intra-assay CV), which is generally accepted as good [45]. All standards and controls were within the default range, as stated by the manufacturer, with the exception that in one sRANKL ELISA run the low control was 6% above the reference range.

Statistical analysis
Statistical analysis was conducted using GraphPad Prism version 8.4.3. for MacOS (GraphPad Software, San Diego, California, USA). For sRANKL, we excluded data from 13 participants due to technical measurement errors, as serum levels were below the detection rate of the ELISA test. For OPG, we excluded data from 2 participants due to a technical measurement error (no concentrations were provided by the ELISA kits). We thus had data available from 20 participants in the IG and 27 participants in the CG for OPG as well as from 17 participants in the IG and 19 participants in the CG for sRANKL ( Fig. 1).
To test for normal distribution, we used the Shapiro-Wilk test [46] as well as histograms for graphic estimates. sRANKL as one of the key parameters in this analysis did not show normal distribution; therefore, to remain consistent, we used non-parametric tests and presented all metric data as median and IQR. To test whether two independent groups were statistically different, we used the Mann-Whitney U test. For categorical variables, we used the Fisher's exact test. To test for differences within each study group before and after the 3-month intervention phase (ΔV1-SE), we used the Wilcoxon matched-pairs signed-rank test. Associations between two variables were tested using Spearman's rank correlation coefficient (r); associations between more than two variables were tested using multiple linear regressions. Multiple testing was adjusted with a false discovery rate of 5% according to Benjamini and Hochberg [47]. An adjusted p value (q) < 0.05 was considered as statistically significant.
Some sRANKL concentrations yielded a strong positive deviation from the median, resulting in large variances. We conducted statistical outlier analyses using the Tukey method for Boxplots [48], identifying 8 data sets that could be regarded as outliers. However, we decided not to exclude these data for several reasons: firstly, no clinical reference values are yet available for the biomarkers in the context of hereditary BC, not allowing assumptions with regard to clinical implausibility. Secondly, we did not find an indication for measurement errors; hence, by including all measurable concentrations, we aimed at preventing arbitrariness. Thirdly, the comparison of data sets when excluding outliers achieved similar results.

Study population
Regarding all participant characteristics included in the present analysis, the study groups IG and CG did not differ significantly at SE (Table 1). To test whether our subgroup of 49 individuals differed from the total LIBRE-1 cohort (n = 68) in any of the characteristics shown in Table 1, we performed Mann-Whitney U tests between these cohorts, which showed no statistically significant results (data not shown).
Individuals from the IG and CG showed a similar significant increase in OPG between SE and V1 ( Table 2, Fig. 2A). Furthermore, the IG showed a significant decrease in sRANKL, while the decrease in this biomarker was nonsignificant for the CG ( Table 2, Fig. 2B). The extent to which OPG increased and sRANKL decreased was similar in the IG and CG.
For VO 2 peak, while non-significant within and between groups, the IG showed an improvement in VO 2 peak while the parameter in the CG stagnated from SE to V1 (Fig. 3A). Both groups (IG: p = 0.0001; q = 0.0002), CG: p = q = 0.0165) similarly increased their MEDAS scores significantly from SE to V1 (Fig. 3B).
The correlation analyses of PUFA with the biomarker included n = 26 participants from both study groups, as PUFA analyses were only available for these participants. The correlation analyses showed that an increase in total n-3 PUFA in RBCM was significantly correlated with an increase in OPG (Fig. 4), as well as an increase in the n-3 EPA with an increase in OPG (Supplementary Table 1). When looking at the study groups individually, an increase of n-3 PUFA correlated with an increase in OPG in the IG, yet the correlation was not statistically significant after FDR correction (r = 0.6176; p = 0.0212; q = 0.1201). There was no  There was no significant correlation between the changes in n-3, n-6, and n-9 PUFA measured in RBCM or plasma and the changes in sRANKL (Supplementary Table 1). Before FDR correction, both the total n-6 PUFA and oleic acid showed significant associations to sRANKL and OPG when measured in plasma. While an increase in total n-6 PUFA was positively correlated with an increase in sRANKL and negatively correlated with OPG, oleic acid showed the opposite trend (Supplementary Table 1).
For sRANKL, the explorative multivariate regression models showed a significant association between the baseline serum level and the changes in the biomarker from SE to V1 (Table 3). Furthermore, the analysis indicated that past smoking might be related to a stronger increase in sRANKL, yet this association was only significant before the p value was adjusted for multiple testing. For OPG, the multivariate analysis showed a significantly positive relationship between baseline VO 2 peak and changes in OPG, yet this association was no longer statistically significant after adjusting for multiple testing.

Discussion
Experimental evidence suggests that the hyperactivation of the RANK/RANKL/OPG pathway plays an important role in the differentiation of the mammary gland [9], stem cell expansion [49], and ultimately in mammary carcinogenesis [10] in gBRCA1/2 mutation carriers. In the present analysis, we explored the potential effects of a structured lifestyle intervention on OPG and sRANKL serum concentrations and demonstrated that the MedD in combination with increased physical activity may be biomarker modulators. To the best of our knowledge, this is the first time that lifestyle intervention-based alterations in OPG and sRANKL serum  concentrations were shown in the context of gBRCA1/2 mutation carriers. Several studies propose that OPG can be increased directly through physical activity by the mechanism of mechanical bone stress [40,50,51]. While Ziegler et al. showed that OPG concentrations increased immediately after an endurance exercise [27], Bergström et al. found significantly elevated OPG levels after a structured 1-year endurance training [28]. Hence, the RANK/RANKL/OPG pathway seems to be responsive to acute as well as chronic exercise [52]. Our data support these findings by showing a significant increase in OPG serum concentrations over a time span of 3 months, indicating a direct responsiveness of the biomarker to changes in health behavior related to physical activity.
In contrast to OPG, the influence of physical activity on sRANKL concentrations is still unclear [28]. Some studies show that exercise can lower sRANKL serum concentrations, but changes are likely to depend on the intensity of exercise [53,54]. Ziegler et al. demonstrated that sRANKL decreased significantly only when running a marathon, but not a 15 km distance [27]. Scott et al. reported that while endurance training instantly increased OPG, only highintensity endurance exercise caused an additional decrease in sRANKL [55]. Our results support these findings, where we showed diminishing sRANKL concentrations within both study groups, but only significantly in the IG (Fig. 2B). This might indicate that only a structured and supervised intervention program like in the LIBRE-1 IG provides the intensity needed to effectively change sRANKL concentrations. However, since the intergroup difference in sRANKL changes was not significant, the LIBRE-1 lifestyle intervention might not be sufficient to significantly lower sRANKL levels. While VO 2 peak increased in the IG and decreased in the CG (Fig. 3A), indeed indicating a measurable improvement in cardiopulmonary fitness through the lifestyle intervention, there is a need for further studies with a more intense exercise design with either more frequent training sessions or higher intensities. After all, it remains unclear whether the diminishing effect of exercise on sRANKL is a direct mechanism or whether it is a result of an increase of inhibitory OPG levels [27].
Apart from exercise, OPG and sRANKL seem to be responsive to dietary habits. Various nutritional studies show a modulating impact of PUFA on these biomarkers [56]. In the present analysis we showed a significant association between the intake of n-3 PUFA, specifically EPA, and OPG (Fig. 4). Our results indicate that a structured nutritional intervention might change PUFA intake patterns sufficiently, as only the IG showed a significant association to n-3 PUFA in subgroup comparison to the CG (albeit the non-corrected p value), which supports other PUFA analyses on the LIBRE-1 cohort [44]. Our results are supported by several studies, such as Martin-Bautista et al. who reported an increase in OPG by 18% and decrease in sRANKL by 7% after a 1-year consumption of n-3 PUFA [57]. The latter negative correlation between n-3 PUFA and sRANKL was also described by several other studies [58,59]; however, our data were not able to replicate these results. This might be due to the rather large variances in our sRANKL data and might be aggravated by our small sample size. For n-6 PUFA, which are consumed less when following the MedD, a negative association to OPG as well as a positive association to sRANKL is described [60]. By diminishing the OPG/sRANKL expression, the intake of n-6 PUFA has been associated with an increased inflammatory and osteoclastic activity and hence risk of cardiovascular diseases and osteoporosis [61]. While being non-significant after adjusting for multiple testing, our correlations show congruent results for n-6 PUFA as well as for the n-9 PUFA oleic acid, for which olive oil is a primary dietary source and is known to inhibit osteoclastogenesis. In summary, our correlation analyses seem to confirm the associations between the different PUFA and OPG and sRANKL that are described by various studies using prospective data of gBRCA1/2 mutation carriers and suggest that the MedD might modulate the biomarker serum levels through modification of PUFA intake.
Furthermore, sRANKL is reported to depend on the metabolic status and to diminish as a response to a reduction of adipose tissue and thus circulating sex hormones, insulin, proinflammatory cytokines, and the adipocytokine leptin [62,63]. It has been shown that the MedD can improve the metabolic status [64,65]; hence the nutritional intervention, measured as a significant increase in the MEDAS score (Fig. 3B), might also explain the observed decrease in sRANKL concentrations in the IG. Pasanisi et al. showed in a similar study on gBRCA1/2 mutation carriers that a dietary MedD intervention can lead to a reduction of sRANKL through reduced levels of proinflammatory IGF-1 [32].
We performed exploratory multiple linear regression models to examine whether participant characteristics at SE might influence biomarker changes during the study. To the best of our knowledge, this is the first time that the dynamics of changes in OPG and sRANKL are assessed with regards to their initial serum levels. For sRANKL, we identified the baseline level as a significant influence on the change of the serum levels, but not for OPG, suggesting baseline sRANKL levels as a potential candidate to predict the outcome of a 3-month lifestyle intervention. For gBRCA1/2 mutation carriers with initially high sRANKL concentrations, the biomarker seems to decrease stronger over the intervention phase compared to women with lower initial sRANKL serum levels. The fact that sRANKL is associated with a proinflammatory metabolic status [66], indicates that especially participants with a less active lifestyle might benefit from such an intervention program. In turn, there might be a saturation effect for participants with an already healthy lifestyle and initially low sRANKL levels, leading to a smaller decrease in the biomarker through an additional change in health habits. Hence, initial biomarker levels might also affect the modulative potential of lifestyle changes. Multivariate regression models also showed that cigarette smoking might influence the change of sRANKL in the course of the study. It has been shown that smoking can induce the expression of RANKL mRNA [67] and increases the sRANKL/OPG ratio by promoting inflammatory cytokines [68]. In a large case-control study it was demonstrated that cigarette smoking increased the risk for BC in gBRCA1/2 mutation carriers [69] for which sRANKL/OPG system might present a molecular mechanism of action. According to our data, baseline physical fitness (assessed as VO 2 peak in CPET) was associated with changes in OPG. The higher the participants' initial VO 2 peak, as a proxy for cardiopulmonary fitness, the more pronounced OPG levels increased during the study. This finding supports the notion that OPG is primarily responsive to physical activity [27] and might indicate a dose-response relationship between cardiopulmonary fitness and OPG serum concentrations. In all, while our results must be regarded as purely explorative, we aim to confirm these findings using prospective data from the larger LIBRE-2 study.
While biomarker changes were more pronounced in the IG, we did not find significant intergroup differences in biomarker changes in any of our analyses. This suggests that all participants, independent of the study group allocation, seem to have made efforts within the study program to adapt to healthier lifestyle choices. We assume that information exchange between participants due to the open design of the study, as well as possible disappointment within the CG with regards to group allocation, which would be consistent with previous findings from the LIBRE-1 study [38], might explain why some effects that were expected only in the IG also occurred in the CG.
Furthermore, our study cohort was likely to inherit a selection bias toward participants prone to above-average health habits, as it has been shown that a prior cancer disease is likely to have confounding effects on subsequent health behavior [70]. Especially the fact that a gBRCA1/2 mutation or cancer diagnosis often occurs within the familial surrounding, participants-whether themselves diseased or not-are likely to be sensitized for a more proactive engagement toward lifestyle changes [71]. We also did not control for special dietary habits or sport programs of the participants prior to SE, hence we could not account for potentially confounding and/or prolonging effects of individual health behavior on the biomarkers.
The small sample size of this evaluation represents its main limitation. However, this explorative approach aimed at generating novel knowledge regarding a possible association between lifestyle, the biomarkers OPG and sRANKL, and BC in at-risk individuals. The present analysis was conducted as a secondary analysis and was no initial endpoint of the LIBRE-1 study; therefore, study design and sample size were not established based on statistical measures. Hence, we cannot ensure representativeness of our cohort (n = 49) compared to the full LIBRE-1 cohort (n = 68) or compared to the combined LIBRE-1 and LIBRE-2 cohort after both trials will be finished (planned n = 660). These initial findings can be the basis for future studies comprising larger sample sizes. Specifically, by showing that short-term lifestyle modifications can alter serum levels of the observed biomarkers, the present findings will be the foundation for a subsequent analysis of OPG and sRANKL in the larger LIBRE-2 main study.
In the present analysis we only explored the dynamics of OPG and sRANKL in individuals with a gBRCA1/2 mutation. While serum levels of OPG and sRANKL in healthy individuals have been evaluated in a number of clinical studies, the physiological ranges differ greatly. While compared to some studies, the median concentrations for both biomarkers seem to be higher in our cohort than in agematched reference subjects [72,73], both biomarkers were within normal range according to other references [74,75]. Whether these differences compared to healthy women are due to the gBRCA1/2 mutation or other factors is currently unclear and subject to further research.
With our analysis, we focused on the potential influences of exercise and nutritional habits on biomarker serum levels. While highlighting the preventive potential of health habits in reaching favorable biomarker constellations, this approach does not provide prognostic value. So far, there are ambiguous findings whether OPG and sRANKL actually are prognostic markers for BC risk. While Vik et al. reported a significantly inverse association between OPG and BRCA1/2-associated BC risk [76], Kotsopoulos and colleagues did not find evidence for an association between plasma OPG [77] or sRANKL [78] levels and BC risk. However, the validity of their study might be limited as plasma samples were collected in the late 1990s and hence stored over 20 years, possibly altering protein concentrations. Furthermore, Kotsopoulos et al. postulate that a single measurement of OPG and sRANKL is able to predict the BC risk in the future. This stands in contrast to our findings, suggesting the biomarkers are dynamic and responsive to lifestyle changes also in the short term.

Conclusion
The aim of our analysis was to objectively demonstrate the effects of a controlled lifestyle intervention, including physical activity and adherence to a MedD, on serum concentrations of OPG and sRANKL.
Showing that changes in health habits can lead to favorable changes in biomarker levels, especially OPG may provide a first step in research for establishing clinical risk factors as well as potential preventative options for gBRCA1/2 mutation carriers with regards to the incidence of hereditary BC disease.
Acknowledgements We would like to thank all staff members involved in the LIBRE-1 study and especially Alexandra Sturmheit for the great support during the implementation of the laboratory analyses.
Author contributions LN participated in the data collection/laboratory analyses (OPG, sRANKL), statistical analysis, and manuscript writing. MY-D contributed to project development, data preparation and analysis, and manuscript writing. SG performed supervision as project manager and contributed to project conceptualization and development, methodology, and manuscript writing. CE contributed to the conception and design of the LIBRE-1 study and to the acquisition of the data and he revised the manuscript critically for important intellectual content. MH and SCB contributed to conceptualization of the LIBRE-1 study and manuscript editing. MK contributed to PI of the LIBRE-1: conceptualization of the LIBRE-1 study design, project administration, project resources, and manuscript editing. BS contributed to data collection/laboratory analyses (RBCM PUFA), manuscript editing. MB contributed to project development, data collection, and manuscript editing. JL contributed to conceptualization and manuscript editing. AB-E performed CPET data collection and manuscript editing. JR participated in the methodology of laboratory analyses and resources. ASQ contributed to project resources and manuscript editing. TS, UN, KR, and RS contributed to LIBRE-1 project administration and manuscript editing. All authors read and approved the final manuscript to be published. All the authors agreed 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.
Funding Open Access funding enabled and organized by Projekt DEAL. The LIBRE-1 study is funded by the German Cancer Aid (Deutsche Krebshilfe, http:// www. krebs hilfe. de) within the Priority Program "Primary Prevention of Cancer" (Grant No. 110013). The funder has no authority and is not involved in the following activities: study design; collection, management, analysis, and interpretation of data; writing of the report; and the decision to submit the report for publication.

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. Individual participant data will be available after de-identification (text, tables, figures). The study protocol has been published elsewhere.
Code availability Statistical analysis was conducted using GraphPad Prism version 8.4.3. for MacOS (GraphPad Software, San Diego, California, USA).

Conflict of interest
The authors declare no potential conflicts of interest. Informed consent Informed written consent was obtained from all participants prior to study entry. They were informed that they can withdraw their consent and stop participation at any time without disclosing the reasons and without negative consequences for their future medical care. All participants have consented to publishing their data anonymously.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.