Article Text

Original research
Allostatic load as predictor of mortality: a cohort study from Lolland-Falster, Denmark
  1. Neda Esmailzadeh Bruun-Rasmussen1,
  2. George Napolitano2,
  3. Christian Christiansen3,
  4. Stig Egil Bojesen4,
  5. Christina Ellervik5,6,
  6. Randi Jepsen1,
  7. Knud Rasmussen5,
  8. Elsebeth Lynge1
  1. 1 Centre for Epidemiological Research, Nykobing Falster Hospital, Nykobing, Denmark
  2. 2 Department of Public Health, University of Copenhagen, Copenhagen, Denmark
  3. 3 Department of Internal Medicine, Nykobing Falster Hospital, Nykobing, Denmark
  4. 4 Department of Clinical Biochemistry, Herlev Hospital, Herlev, Denmark
  5. 5 Department of Data and Development Support, Region Sjaelland, Soro, Denmark
  6. 6 Department of Clinical Medicine, University of Copenhagen, Copenhagen, Denmark
  1. Correspondence to Dr Neda Esmailzadeh Bruun-Rasmussen; neebruun{at}gmail.com

Abstract

Objectives The purposes of the present study were to determine the association between (1) 10 individual biomarkers and all-cause mortality; and between (2) allostatic load (AL), across three physiological systems (cardiovascular, inflammatory, metabolic) and all-cause mortality.

Design Prospective cohort study.

Setting We used data from the Lolland-Falster Health Study undertaken in Denmark in 2016–2020 and used data on systolic blood pressure (SBP) and diastolic blood pressure (DBP), pulse rate (PR), waist–hip ratio (WHR) and levels of low-density lipoprotein cholesterol (LDL-c), high-density lipoprotein cholesterol (HDL-c), triglycerides, glycated haemoglobin A1c (HbA1c), C-reactive protein (CRP) and serum albumin. All biomarkers were divided into quartiles with high-risk values defined as those in the highest (PR, WHR, triglycerides, HbA1c, CRP) or lowest (HDL-c, albumin) quartile, or a combination hereof (LDL-c, SBP, DBP). The 10 biomarkers were combined into a summary measure of AL index. Participants were followed-up for death for an average of 2.6 years.

Participants We examined a total of 13 725 individuals aged 18+ years.

Primary outcome measure Cox proportional hazard regression (HR) analysis were performed to examine the association between AL index and mortality in men and women.

Results All-cause mortality increased with increasing AL index. With low AL index as reference, the HR was 1.33 (95% CI: 0.89 to 1.98) for mid AL, and HR 2.37 (95% CI: 1.58 to 3.54) for high AL.

Conclusions Elevated physiological burden measured by mid and high AL index was associated with a steeper increase of mortality than individual biomarkers.

  • Biomarker
  • Allostatic Load
  • Blood
  • Mortality
  • population-based
  • LOFUS

Data availability statement

Data are available upon reasonable request. Data from the study can be made available via Region Sjaelland following the Danish Data Protection Regulation.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Analysis based on a large population-based health study.

  • Complete follow-up for death via linkage with Danish Civil Registration System.

  • Biomarkers from only one point in time.

  • No biomarker from neuroendocrine system available.

Introduction

Biological markers (biomarkers) were originally defined as ‘cellular, biochemical or molecular alterations that are measurable in biological media such as cells, human tissues or fluids’.1 Later the definition was extended to include ‘indicators of normal biological processes, pathogenic processes and pharmacological responses to therapeutic interventions’.2 In clinical settings, measurement of biomarkers in blood samples is used to detect and diagnose medical conditions. Biomarkers as independent predictors of all-cause mortality are therefore of considerable clinical and research interest;3 dyslipidaemia including high levels of triglycerides and low-density lipoprotein cholesterol (LDL-c), and low levels of high-density lipoprotein cholesterol (HDL-c), have been reported to be independent risk factors for all-cause mortality.4–6 Lower levels of albumin7 and higher levels of C-reactive protein (CRP),8 and glycated haemoglobin A1c (HbA1c)9 have likewise been linked to mortality. Also, there is some evidence that the relationship between some of these biomarkers and all-cause mortality varies across sex and age groups.10 11

The concept of allostatic load (AL) refers to the ‘wear and tear’ of the body resulting from repeated stimulation of stress responses via the hypothalamic–pituitary–adrenal axis and the sympathetic-adrenal-medullary system.12 As a latent variable, AL cannot be directly measured but it can be estimated using an AL index, which is composite of biomarkers from multiple organ systems integrated into a single score. The first AL developed by Seeman et al in 1997 included 10 biomarkers monitoring various physiological systems.13 However, the type and number of biomarkers used in published studies have ranged from 6 to 24.14 The most frequently used Al construct, originally proposed by Gruenewald et al in 2012,15 includes 24 biomarkers. It has been suggested that in the calculation of AL, the threshold of risk for each biomarker should be obtained by the quartiles or quintiles of the values of the biomarker.16 AL has been reported to be a better predictor of mortality than individual biomarkers, however, there are still gaps in the understanding of the associations.17 18 AL has been suggested also as a tool for allocation of nursing resources.19

This study provides data from the Lolland-Falster Health Study (LOFUS),20 a population-based survey undertaken in 2016–2020 in Lolland-Falster, a rural–provincial region in Denmark with a life expectancy much below the national average,21 and with health problems reported more frequently than in the rest of the country.22 Using the LOFUS data, the purposes of the present study were (1) to determine the association between 10 individual biomarkers and all-cause mortality; and (2) to examine the association between AL, across three physiological systems (cardiovascular, inflammatory, metabolic system) and all-cause mortality. The hypothesis is that AL can be used as an informative tool in predicting future risk of death in the general adult population.

Methods

Study population

We undertook a prospective cohort study of participants from LOFUS; a household-based population study with data collected between February 2016 and February 2020. Persons aged 18 years and above were randomly sampled from the Danish Civil Registration System and invited to participate together with the rest of their households. Participation required informed consent. A detailed description of the study protocol20 and information on the socioeconomic determinants of participation23 have been published previously. Persons below 18 years, and pregnant women were excluded from the present study.

Patient and public involvement

Patients were not actively involved in any stage of the present study. Once the paper has been published in the international literature, the key results will be reported also in the local press.

Self-reported data

From questionnaires, we used data on smoking (never, former, current), and presence of chronic conditions (cardiovascular disease, diabetes, cancer) at the time of participation in LOFUS.

Biomarkers

Non-fasting blood samples were collected in vacutainer blood collection tubes (Becton, Dickinson and Company; Franklin Lakes, New Jersey, USA) and kept at room temperature until same day analysis at the Department of Clinical Biochemistry at Nykøbing Falster Hospital, accredited by the International Organization for Standardization 15189. We used data on HDL-c, LDL-c, triglycerides, albumin, CRP and HbA1c. LDL-c was calculated by using Friedewald formula24 when the plasma triglyceride concentration was below 4.5 mmol/L. Systolic (SBP) and diastolic (DBP) blood pressure were based on three consecutive digital measurements on the upper left arm (apparatus type Welch Allyn Connex proBPO 3400). The mean values of the second and third measurement were used in this study (only one measurement was used if the other was missing). Waist-hip ratio (WHR) was calculated by waist-circumference divided by hip circumference. Body mass index (BMI) was defined as weight in kilograms divided by height in metres squared (kg/m2).

In the calculation of AL, biomarkers are most often dichotomised into low and high values based on either a percentile or a predetermined cut-off value.16 However, before doing so, we mapped for each biomarker the association between level of the marker and all-cause mortality, see Method below. For most biomarkers the association was monotonic (see online supplemental figure 1). These biomarkers were then dichotomised according to the sex-specific and age-specific quartiles. For age, we dichotomised at age 60. Some previous studies focused on AL in people aged 60 and above25 26 and we intuitively found it reasonable to distinguish in the same way between ‘young’ and ‘old’ people in our data; age 60 was furthermore the median age of our study population; and with this age-dichotomisation we avoided violations of the model assumption in the statistical analysis. We dichotomised biomarkers with high-risk values defined as those in the highest quartile of the sex-specific and age-specific distribution, except for HDL-c and albumin, where the lowest quartile was the high-risk value. For LDL-c, SBP and DBP the associations were U-shaped, and the high-risk values for these biomarkers were therefore defined as including both the lower and the upper quartiles (see online supplemental table 1). For biomarkers with U-shaped associations, we tested out also using octiles as cut-off points. However, this resulted in some violations of the model assumptions in the statistical analysis, and for SBP the upper octile cut-off was from a clinical point of view very high. On this basis we used the quartile cut-offs also for the biomarkers with the U-shaped association. For all biomarkers, the highest and lowest quartile of risk scores were either lower or similar to clinical cut-points.27–31

Supplemental material

Supplemental material

BMI was divided into underweight (BMI less than 18.5) normal (BMI 18.5–24.9), overweight (BMI 25.0–29.9) or obese (BMI 30.0 or greater); reported diseases into either present or not; and smoking status into never, former or current.

Allostatic load scores

The AL scores were computed using biomarkers from: the cardiovascular system (CVS) (SBP, DBP and pulse rate (PR)); the metabolic system (MS) (LDL-c, triglycerides, HDL-c, WHR, HbA1c); and the inflammatory system (IS) (CRP, serum albumin).

Each system-specific AL score was then defined as the number of biomarkers with a high-risk value, hence as an integer value between 0 and 3 for the CVS, 0 and 5 for the MS and 0 and 2 for IS. The AL index was defined as the sum of all scores and divided in three groups based on tertiles contrasting individuals with (AL: 0–2), mid (AL: 3–4) and high (AL: 5–10). Note that, all biomarkers were given equal weight in accordance with previous studies.16 18

All-cause mortality

LOFUS participants were followed-up for death with data obtained from the Danish Civil Registration System on 26 February 2021.

Data management and statistical analyses

Observations with missing values in any of the variables were excluded from the analyses (1989 out of 15 714, ie, 12.6%, see online supplemental table 2). Values below the lower limit of detection were replaced with random numbers sampled with replacement from the set Embedded Image , where n is the variable-specific number of decimals reported in the data and Embedded Image the limit of detection (see online supplemental table 3).

Participants were followed-up from date of participation in the LOFUS study until date of death or end of follow-up on 26 February 2021, whichever came first. In order to define the biomarkers’ high-risk values, we first studied the association between levels of each individual biomarker and mortality, allowing for possible non-linear relations. This analysis was carried out via Cox proportional-hazard models with biomarker levels as continuous covariates, modelled with natural cubic splines with 2 df (except for LDL-c, where 3 df were used), and further adjusting for sex and age. By graphical inspection, a U-shaped association was found for LDL-c, SBP and DBP (see online supplemental figure 1). Therefore, for these biomarkers both the sex and age-specific (ie, below or above age 60) lower and upper quartiles were defined as high risk, while only one quartile for the others (upper or lower, in accordance with the existing literature) (see online supplemental table 1).

Associations between all-cause mortality and dichotomised biomarkers levels (low/high risk), system-specific AL scores and total AL index, were modelled with Cox proportional-hazard models. Here, we present two models: Model 1, where HRs are adjusted for sex and age; Model 2, where results are further adjusted for BMI, prevalent diseases and smoking status. HRs for the individual biomarkers (table 1) and for system-specific AL scores (table 2) are mutually adjusted. Proportional hazards assumptions in the above models have been tested using Schoenfeld residuals. Numbers below 5 are not reported. In addition, we report HRs for a one-point increase in the AL index.

Table 2

Multivariate Cox proportional hazard regression of all-cause mortality for Lolland-Falster Health Study participants by allostatic load index

Table 1

Multivariate Cox proportional hazard regression of all-cause mortality for Lolland-Falster Health Study participants by individual biomarkers

Data management, statistical analyses and plots were done in R V.4.0.3,32 with packages splines,32 survival,33 tidyverse,34 ggrepel35 and ggpubr.36

Results

The LOFUS database used for this study included 13 725 persons, of whom 53% were women and 47% men. The median follow-up time was 2.6 years (IQR 1.5) and the median age was 57.6 in women and 59.9 in men. One-fourth of the participants were obese, and one-fifth were current smokers. Presence of cardiovascular disease at the time of LOFUS participation was reported by 28%, diabetes by 5% and cancer by 4%. On the value of total AL index, participants were divided approximately into tertiles; 32% low, 40% mid and 38% high. During the follow-up period, 198 participants died; of these 39% were women and 61% men (table 3).

Table 3

Lolland-Falster Health Study (LOFUS). Baseline characteristics of study population and deaths in follow-up period, n(%). For definition of cut-off values, see online supplemental table 1

The multivariate Cox proportional hazard regression for individual biomarker and all-cause mortality, adjusted for sex and age and additionally for BMI, reported diseases and smoking, are listed in table 1. For all biomarkers, apart from triglycerides, a high-risk value was associated with an increased mortality level. However, only the HRs for low albumin and high CRP were statistically significantly elevated; HR 1.54 (95% CI: 1.16 to 2.06) and 1.41 (95% CI: 1.04 to 1.91), respectively.

The HR for all-cause mortality increased with increasing level of the AL from low as the reference over mid to high, table 2 and figure 1. For the IS AL score, the HR was 1.02 (95% CI: 0.73 to 1.42) for mid AL, and 2.38 (95% CI: 1.67 to 3.39) for high AL. For the MS AL score, the HRs were 1.18 (95% CI: 0.75 to 1.85) and 1.54 (95% CI: 1.00 to 2.38), respectively. For the CVS AL score, the HRs were 1.65 (95% CI: 1.02 to 2.65) and 1.89 (95% CI: 1.20 to 2.99), respectively. The gradient for the total AL index was a HR of 1.33 (95% CI: 0.89 to 1.98) for mid AL, and 2.37 (95% CI: 1.58 to 3.54) for high AL. HRs for one unit increase in AL (continuous AL) was 1.23 (1.14–1.32) when adjusted for age and sex, and 1.22 (1.13–1.32), when additionally adjusted for BMI, reported diseases and smoking status.

Figure 1

All-cause mortality by level of allostatic load index, as HR (95% CI). AL, allostatic load.

Discussion

In this population-based study from a rural–provincial area of Denmark, we followed the adult population up for a median period of 2.6 years. High levels of individual biomarkers were overall associated with increased mortality, but most of them at a modest level of 20%–30%, and statistically significantly elevated for only CRP and albumin. High levels of physiological system-specific AL scores were associated with increased mortality at the level of 50%–140%; statistically significantly for the IS and CVS, and at borderline of significance for the MS. The composite measure of total AL index was a strong predictor of all-cause mortality. Persons with a high versus low total AL index had about 2.5 times the mortality. The total AL index was thus a better predictor of all-cause mortality than individual biomarkers and the MS and CVS AL scores, a pattern consistent with previous studies.16 18 37

The most comprehensive studies on AL and mortality all used data from the National Health and Nutrition Examination Survey (NHANES). Borrell et al 38 examined 12-year mortality by using data from 13 715 adults aged 25+ years of whom 2491 had died. They calculated AL based on nine biomarkers; albumin, CRP, total cholesterol, HDL-c, HbA1c, WHR, SBP, DBP and PR. Using a clinical cut-off AL score, they found that, compared with persons with an AL score of ≤1, those with AL scores of 2 and 3+ had adjusted HRs of 1.40 (95% CI: 1.11 to 1.76) and 1.88 (95% CI: 1.56 to 2.26), respectively.

Levine and Crimmins39 examined 10-year all-cause and disease-specific mortality. In total, 15 042 persons were eligible, but biomarker data were available for only 9942 adults aged 30+, of whom 1076 had died. They included data on albumin, CRP, WHR, total cholesterol, HDL-c, HbA1c, PR, SBP and DBP. For each of the nine biomarkers, a person was classified as high or low based on clinical cut-off points, and the AL score was the number of biomarkers classified as high. In addition, an expanded AL score included 5 additional biomarkers defined by quintiles; and a continuous AL score used a continuous z-score measure for all 14 biomarkers. For the first AL score, a HR of 2.75 (p<0.001) was found for all-cause mortality when persons with the highest quintile of AL were compared with those with the lowest. Somewhat stronger gradients were found for the expanded; 3.62 (p<0.0001) and continuous; 6.97 (p<0.0001), ALs.

Howard and Sparks40 studied 11 733 participants from NHANES. Imputation was used to estimate missing values. Their AL measure was based on DBP, SBP, PR, total cholesterol, HDL-c, triglycerides, HbA1c, BMI, albumin and CRP. They found that a one-unit increase in AL represented a 7% increase in risk of death when adjusted for age, sex, ethnicity, socioeconomic status and health behaviour.

The National Child Development Study was followed-up for deaths from birth in 1958 to 1 December 2013, that is, to the age of 55 years.18 AL based on 10 biomarkers was calculated and divided into three levels. All-cause mortality for persons with mid or high AL was compared with that of persons with low AL, and adjusted for early life, childhood, young and adulthood confounders. The HR of death was 1.71 (95% CI: 1.07 to 2.72) for persons with mid AL, and 2.57 (95% CI: 1.59 to 4.15) for those with high AL. The association between AL and all-cause mortality was stronger than the associations between of the individual 10 biomarkers and all-cause mortality.

The NHANES studies vary in number of participants included in the studies, in length of follow-up for mortality, in biomarkers included, in the definition of AL and in methods used for AL calculation. Nevertheless, all the studies indicated that all-cause mortality increased with increasing AL. The study by Borell et al 38 is the one methodologically most similar to our study and the gradient of 1.88 (95% CI: 1.56 to 2.26) is compatible with the 1 of 2.37 (95% CI: 1.58 to 3.54) found in our study, and so is the gradient of 2.57 (95% CI: 1.59 to 4.15) found in the National Child Development Study.

For individual biomarkers in our study, HRs were highest for CRP and albumin. CRP is the prototypical acute-phase response protein that increases during systemic inflammation,41 while albumin is a major component of plasma protein, required for transportation and to maintain oncotic pressure, acid–base function, microvascular permeability and to prevent platelet aggregation.42 Inflammation increases capillary permeability and thereby escape of serum albumin, leading to expansion of interstitial space and increasing the distribution volume of albumin causing lower serum albumin concentrations. High level of CRP and low level of albumin have thus previously been linked with a variety of health outcomes including morbidity and mortality.7 8 43

We found a U-shaped association between LDL-c and mortality. Elevated LDL-c is a well-established risk factor of atherosclerosis and cardiovascular disease, and the general perception is that high level of LDL-c is associated with an increased risk of morbidity and mortality.44 45 Nevertheless, studies on the association between LDL-c levels and mortality have provided conflicting results. Some studies found increasing level of LDL-c to be associated with lower mortality,5 ,46 and some studies found no association.45 47 48 However, most studies were conducted in elderly people often with an intake of lipid-lowering agents. A more recent study in young Koreans found an association between low level of LDL-c and an increased risk of cancer, cardiovascular and all-cause mortality.49 These findings were supported by a Chinese study of participants aged 40+ years.50 A recent Danish study among 108 243 individuals aged 20–100 years found the lowest all-cause mortality at an LDL-c concentration of 3.6 mmol/L (140 mg/dL), and higher mortality at both lower and higher levels.51 Our findings for LDL-c were thus in accordance with these recent observations. Seplaki et al suggested that both high and low ends of the risk continuum for the construct of AL could be more informative than simply using high-risk quartiles. They assigned a value of ‘1’ for values above the 75th percentile and below the 25th percentile of the distribution, and a value of ‘0’ for intermediate values.52

We found both higher and lower levels of DBP to be associated with an increased mortality, and a similar tendency was indicated for SBP. The association between lower BP and mortality is still of discussion.53–55 Most studies have found this association among elderly people and linked it to chronic disease, for example, cardiovascular disease (cardiac failure or ischaemic heart disease), cancer, poor functional status or frailty. Low BP has also been associated with poor function and low quality of life,56 57 but in previous studies only the highest quartile or the clinical cut-off value have been used as predictor of all-cause mortality.

Several methods have been used to define an AL composite index, including the count-based, canonical correlation, z-score and grade of membership method.58 59 The most commonly used method is the count-based method, where a summary index is calculated by summing the number of biomarkers falling within the high-risk category, either defined by the percentile (ie, upper or lower 25th percentile of the sample’s distribution) or by the clinical cut-off value. In our analysis with the two-tail cut-off points, we found HRs for LDL-c of 1.13 (95% CI: 0.85 to 1.51); for SBP of 1.17 (95% CI: 0.88 to 1.57; and for DBP of 1.28 (95% CI: 0.95 to 1.72). If we have used instead the single high-risk quartile cut-off point, we would have found HRs for LDL-c of 0.71 (95% CI: 0.49 to 1.03); for SBP of 0.96 (95% CI: 0.68 to 1.35) vs), and for DBP of 1.24 (95% CI: 0.86 to 1.81). The two-tail cut-off points thus provided a better identification of persons with high mortality than the one-tail cut-off points.

The issue of whether a clinical or sample-based cut-off criteria should be used is still of discussion,17 however, studies comparing distinct measurement approaches have found only modest differences in their predictive utility.15 60 61

Strengths and limitations

The strengths of our study included the size of the cohort in terms of the large number of individuals recruited from a general adult population, and the complete follow-up for death by linkage with the Danish Civil Registration System.

Our study also had some limitations. First, the choice of biomarkers used to construct the AL index. The AL theory emphasises the importance of measuring dysregulation across different physiological systems, including biomarkers from the neuroendocrine, CVS, MS and IS.13 The neuroendocrine system (stress response) is believed to play a key role in allostasis and subsequent AL, as a series of physiological changes takes place before initial stress responses occur (such as rapid increases in blood sugar and BP that supply the body with additional energy). However, biomarkers from the neuroendocrine system are difficult to measure, as repeated measurements over 1–2 days are recommended. These requirements cannot be fulfilled in population studies, where participants are examined only once, and biomarkers from the neuroendocrine system were therefore not available for our study.

Second, the initial stress responses are followed by secondary outcomes from the MS, IS and CVS, and these markers were all available in our data. Nevertheless, greater sensitivity could have been achieved by studying the dynamic changes over time in these markers to fully capture the flexibility of stress response mechanisms across the lifespan.

Finally, differences across studies in construction of AL indices could influence the comparison of results. We used the shape of the association between level of a given biomarker and all-cause mortality as the basis for the categorisation of the biomarker into low and high values. One can argue therefore that our analysis was circular in the way that we used outcome on the dependent variable to categorise levels of the independent variable. We believe that this was justifiable in the context here where the purpose was to optimise the predictive power of the AL index. However, validation in other data sets are needed before our approach can be recommended for research in general and for eventual clinical use.

Conclusion

Our findings demonstrated that an optimally constructed AL index was a strong predictor of all-cause mortality. This supported the conceptual validity of AL as an effective marker of the cumulative physiological burden on the body. These findings can contribute to the evidence for the use of an AL index as a basis for targeted efforts to bring down continued stress exposures, and in this way prevent the potential detrimental effect of these exposures on health. Our findings on the U-shaped association with LDL-c, DBP and SBP and all-cause mortality suggested that AL measures incorporating risks at both the low and the high end of biomarkers may yield the best prediction of all-cause mortality.

Data availability statement

Data are available upon reasonable request. Data from the study can be made available via Region Sjaelland following the Danish Data Protection Regulation.

Ethics statements

Patient consent for publication

Ethics approval

Approval was granted by Region Zealand’s Ethical Committee on Health Research (Reg: SJ-421). All data storage and management were approved by the Regional Data Protection Agency of Zealand (REG-024-2019 and REG-24-2015). LOFUS is registered in ClinicalTrials.gov (NCT02482896). Participants gave informed consent to participate in the study before taking part.

Acknowledgments

We thank inhabitants of Lolland-Falster for their participation in LOFUS and for their permission to use their data for research purposes. The Lolland-Falster Health Study (LOFUS), Nykøbing Falster Hospital, Denmark, is a collaboration between Region Zealand, Nykøbing Falster Hospital, and Lolland and Guldborgsund Municipalities. The authors are grateful to LOFUS for making the LOFUS research data available. However, LOFUS bears no responsibility for the analysis or the interpretation conducted within this study.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Twitter @RandiJepsen

  • Contributors All authors contributed significantly to the study. RJ provided the LOFUS data. NEB-R, EL and GN designed the study, interpreted the data and drafted the manuscript. GN performed the statistical analysis. CE, CC, RJ, KR and SEB contributed to the interpretation and writing of the manuscript. All authors critically revised and approved the final manuscript. EL is the guarantor.

  • Funding (1) Region Zealand/University of Copenhagen, (2) Nykøbing Falster Hospital, (3) Professor grants for Elsebeth Lynge and (4) Danish Health Insurance Foundation 19-B-0188. The funding bodies had no role in the design of the study, neither in the collection, analysis and interpretation of data, nor in the writing of the manuscript.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.