Article Text

Original research
Prognosis of light chain amyloidosis: a multivariable analysis for survival prediction in patients with cardiac involvement proven by endomyocardial biopsy
  1. Matthias Aurich1,
  2. Julian Bucur1,
  3. Johannes A Vey2,
  4. Sebastian Greiner1,
  5. Fabian aus dem Siepen1,
  6. Ute Hegenbart3,
  7. Stefan Schönland3,
  8. Hugo A Katus1,
  9. Norbert Frey1 and
  10. Derliz Mereles1
  1. 1Department of Internal Medicine III, Cardiology, Angiology and Pneumology, University Hospital Heidelberg, Heidelberg, Germany
  2. 2Institute of Medical Biometry, University of Heidelberg, Heidelberg, Germany
  3. 3Department of Internal Medicine V, Hematology, Oncology and Rheumatology, University Hospital Heidelberg, Heidelberg, Germany
  1. Correspondence to Dr Derliz Mereles; derliz.mereles{at}med.uni-heidelberg.de

Abstract

Background Cardiac involvement is a main determinant of mortality in light chain (AL) amyloidosis but data on survival of patients with cardiac AL amyloidosis proven by endomyocardial biopsy (EMB) are sparse.

Methods This study analysed clinical, laboratory, electrocardiography and echocardiographic parameters for their prognostic value in the assessment of patients with AL amyloidosis and cardiac involvement. Patients with AL amyloidosis who had their first visit to the amyloidosis centre at the University Hospital Heidelberg between 2006 and 2017 (n=1628) were filtered for cardiac involvement proven by EMB. In the final cohort, mortality-associated markers were analysed by univariate and multivariable Cox regression. Cut-off values for each parameter were calculated using the survival time.

Results One-hundred and seventy-four patients could be identified. Median overall survival time was 1.5 years and median follow-up time was 5.2 years. At the end of the investigation period, 115 patients had died. In multivariable analysis, New York Heart Association-functional class >II (HR 1.65; 95% CI 1.09 to 2.50; p=0.019), left ventricular global longitudinal strain (HR 1.12; 95% CI 1.03 to 1.22; p=0.007), left ventricular end-systolic volume (HR 1.02; 95% CI 1.01 to 1.03; p=0.001), systolic pulmonary artery pressure (HR 0.98; 95% CI 0.96 to 0.99; p=0.027), N-terminal pro-B-type natriuretic peptide (HR 1.57; 95% CI 1.17 to 2.11; p=0.003) and difference in free light chains (HR 1.30; 95% CI 1.05 to 1.62; p=0.017) were independently predictive.

Conclusion Among all patients with AL amyloidosis those with cardiac involvement represent a high-risk population with limited therapy options. Therefore, accurate risk stratification is necessary to identify cardiac amyloidosis patients with favourable prognosis. Incorporation of modern imaging techniques into existing or newly developed scoring systems is a promising option that might enable the implementation of risk-adapted therapeutic strategies.

  • Cardiomyopathy, Restrictive
  • Heart Failure, Systolic
  • Echocardiography

Data availability statement

Data are available on reasonable request.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/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.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Among all patients with light chain (AL) amyloidosis those with cardiac involvement represent a high-risk population with limited therapy options. Although current staging systems are powerful prognostic tools for AL amyloidosis in general, the prediction of survival in patients with cardiac involvement is not equally well investigated and patient data for endomyocardial biopsy-proven cardiac AL amyloidosis is sparse.

WHAT THIS STUDY ADDS

  • This study could demonstrate that a relevant number of patients with cardiac amyloidosis shows long-term survival over many years. To separate patients according to prognosis, we identified different clinical, laboratory and echocardiography parameters that are independent predictors of survival in cardiac AL amyloidosis.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • Strain imaging holds great potential to identify cardiac amyloidosis patients with favourable prognosis as it provides direct and objective read-outs on cardiac structure and function. Modern imaging techniques might, therefore, be incorporated into existing or newly developed scoring systems to enable the implementation of risk-adapted therapeutic strategies.

Introduction

Amyloidosis comprises a group of diseases characterised by extracellular deposition of amyloid, a fibrillar aggregation of at least 30 different misfolded precursor proteins.1 Light chain (AL) amyloidosis accounts for more than three-quarters of cases in patients with systemic amyloidosis and among those, another 75% suffer from cardiac involvement.2 With an incidence of 1.2 cases per 100 000 person-years, AL amyloidosis is a rare disease, but a trend towards higher incidence rates with time might reflect an increased awareness of the disease and diagnostic improvements.3 4 Nevertheless, up to 37% of patients are diagnosed more than 1 year after the onset of initial symptoms.5

In the case of AL amyloidosis, cardiac deposition of misfolded immunoglobulin light chain-derived amyloid fibrils disrupts tissue integrity and causes stiffening of the myocardium eventually leading to restrictive cardiomyopathy and congestive heart failure (HF).6 Prefibril oligomers, in addition, have direct toxic effects on cardiomyocytes and thus amplify organ dysfunction even further.7–12

Amyloid deposits and myocardial injury can cause elevation of cardiac biomarkers such as N-terminal pro-brain natriuretic peptide (NT-proBNP) or cardiac troponins.13–15 These biomarkers have accordingly been incorporated into different risk classification scores.16–19 While these staging systems are powerful prognostic tools for AL amyloidosis in general, the prediction of survival in patients with endomyocardial biopsy (EMB) proven cardiac involvement is not equally well investigated.

Cardiac involvement is commonly associated with a poor prognosis of AL amyloidosis. Although new therapy options have improved overall survival (OS) over the last two decades,20 mortality is still high among patients with confirmed or suspected cardiac involvement and after onset of HF, outcome continues to deteriorate.21 However, some high-risk patients survive even for years,22 23 raising the question whether variables not included in the established scoring systems may indicate a potentially superior therapy response and survival. Therefore, we investigated the predictive value of different clinical, laboratory, electrocardiography and echocardiograpic parameters in an initially therapy-naïve cohort of patients with AL amyloidosis and cardiac involvement proven by EMB.

Methods

Study design and population

The present study is a single-centre investigation based on register datasets collected at the amyloidosis centre of the University Hospital Heidelberg, Germany. The records of all patients with AL amyloidosis who had their initial visit from 27 July 2006 to 27 December 2017 were filtered for cardiac involvement that was confirmed by an EMB. In our cardiology department EMB was performed in patients referred for further invasive HF diagnostics when coronary artery disease could be excluded during heart catheterisation or in case of uncertainty of cardiac amyloidosis based on basic diagnostics (eg, imaging, biomarkers) collected in the haematology or cardiology department.

To exclude a potential bias caused by former amyloidosis treatment or by diseases other than amyloidosis that would lead to altered functional imaging parameters and elevated cardiac biomarkers, patients who were already on chemotherapy for AL amyloidosis at the time of presentation and patients with significant cardiovascular or renal comorbidities were excluded from further analysis. To ensure optimal conditions for echocardiographic measurements, in particular strain imaging, only patients in sinus rhythm and with sufficient image quality were included.

Diagnosis of AL amyloidosis and cardiac involvement

Monoclonal gammopathy was confirmed by immunoelectrophoresis and immunofixation of serum and urine. Systemic AL amyloidosis with cardiac involvement was verified by light-microscopy immunohistochemistry and Congo red staining of myocardial biopsies showing green birefringence under polarised light.

Clinical, laboratory, electrocardiography and echocardiography assessment

Data on demographics, medical history, physical examination, blood samples and urine samples were collected on the day of first presentation at the Heidelberg University Hospital. Additionally, a 12-lead ECG and a comprehensive echocardiography examination were performed in every patient.

Electrocardiography was analysed for heart rhythm and heart rate, cardiac axis, time intervals and low-QRS-voltage pattern (QRS complex amplitude <5 mm in any limb lead).

Echocardiography examinations were conducted using commercially available ultrasound machines (Vivid E9 BT 11, GE Vingmed Ultrasound, Horten, Norway and iE33, Philips Medical Systems, Andover, Massachusetts, USA). Three consecutive heart cycles were acquired for each two-dimensional image, with a sampling rate of 55–60 frames/s. Offline analysis was conducted according to the recommendation of the American Society of Echocardiography and the European Association of Cardiovascular Imaging24 using Image Arena including the 2D cardiac performance analysis software (TOMTEC Imaging Systems, Unterschleissheim, Germany).

Follow-up data were retrieved from digital clinical records and/or telephone interviews with the patients, their relatives or their family doctors.

Statistical analysis

The study cohort was described using summary measures of the empirical distribution. Nominal variables are reported as absolute and relative frequencies and continuous variables are given as medians with 25% and 75% quantiles. In the comparison of survivors and nonsurvivors, for nominal variables, the χ2-test was used, and for continuous variables, the t-test for independent groups or the Mann-Whitney-Wilcoxon-test was used, as appropriate.

OS was investigated using the Kaplan-Meier estimator and the median follow-up was calculated with the reverse Kaplan-Meier method. The variables’ association with OS was examined by calculating HRs applying univariate Cox regression models, whereby NT-proBNP, high sensitivity Troponin T (hsTnT) and difference of free light chains (dFLC) were log-transformed due to their highly skewed distributions. For multivariable Cox regression modelling, the missing values (occurring in total in 1.82% of the investigated variables) were multiple (n=50) imputed. On each imputed dataset, forward selection using Cox regression was performed. In each step, the coefficients were pooled using Rubins’s rules and variable selection was done by using a criterion of p<0.1 for the pooled regression coefficients.25

Optimal cut-off values for risk stratification were calculated by optimising the log-rank test statistic and a 95% bootstrap CI was computed as certainty measure.26

Since this was a retrospective data analysis, all p values should be interpreted in a descriptive manner; p values<0.05 were accepted as significant. All statistical analyses were performed using the software R V.≥4.10.2.

Results

From 27 July 2006 to 27 December 2017, 1628 patients with AL amyloidosis had initial contact with the amyloidosis centre at the University Hospital Heidelberg. Figure 1 depicts a flow chart of the inclusion and exclusion criteria for this study. Only patients who were chemotherapy naïve for AL amyloidosis were included. Cardiac involvement was confirmed by EMB in 216 cases. To minimise the impact of confounders, patients with known cardiovascular (atrial fibrillation n=2, cardiac decompensation n=2, complete heart block n=2, pacemaker n=10, former heart surgery n=2, heart transplantation during follow-up n=8) or severe kidney disease (renal artery stenosis n=1, acute kidney failure n=2, dialysis n=2) as well as patients with suboptimal echocardiography image quality (n=11), were excluded. The final study population consisted of 174 patients. At the end of the follow-up period in August 2021, 115 patients had died.

Figure 1

Flow chart of the inclusion and exclusion criteria of the study. AL, amyloid light chain.

Characteristics of the study population

The clinical, laboratory, electrocardiography and echocardiograpic parameters of the entire study cohort, of survivors and non-survivors are summarised in table 1. Male subjects predominated. Survivors had higher systolic blood pressure and estimated glomerular filtration rate (eGFR) whereas New York Heart Association functional class (NYHA-FC), serum levels of dFLC, NT-proBNP, hsTnT, left ventricular (LV) global longitudinal strain and the ratio of E/e’ were lower. Furthermore, pericardial effusion was less often detected in survivors. Online supplemental table 1 compares patients who survived at least 1 year with patients who died within 1 year after initial presentation, excluding censored cases. Cardiac staging according to the MAYO 3b18 score is presented in table 2. Among survivors most patients were in stages II and III and among non-survivors most patients were in stages III and IIIb. None of the patients met the criteria for stage I.

Table 1

Characteristics of the study population

Table 2

Cardiac staging according to the MAYO 3b Score18

In addition to cardiac involvement, at least one other organ was affected in 41% of the patients, and at least two other organs were affected in 10%. Soft tissue was the second most often involved organ system (24%), followed by gastrointestinal system and kidneys (17% and 6%, respectively). In our study, no difference was observed in the type or number of organs involved between survivors and non-survivors (online supplemental table 2). There was also no difference between both groups with regard to the underlying haematological disease. However, survivors had lower bone marrow plasma cell infiltration (online supplemental table 3).

Throughout the follow-up period, 94% of the patients received specific AL amyloidosis therapy. The therapy regimens are shown in online supplemental table 4. Bortezomib plus dexamethasone was the most frequently used treatment (53%). Other drugs commonly used in various combinations were cyclophosphamide, lenalinomide and melphalan. Bendamustine, doxorubicin and rituximab were ordered less often. One patient received stem cell transplantation (SCT). All others were classified as unsuitable for SCT due to advanced cardiac involvement. Ten patients remained untreated, all of whom belonged to the non-survivor group (five died before the therapy decision, three refused therapy, two were too old and in a significantly reduced general condition).

Survival analysis

The median OS time was 1.47 years (95% CI 1.04 to 2.92) or 17.6 moths (95% CI 12.5 to 35) figure 2. The median follow-up time was 5.16 years (95% CI 4.35 to 7.22) or 61.9 months (95% CI 52.2 to 86.6). The results of univariate Cox regression analysis and the final multivariable Cox regression model resulting from variable selection are presented in table 3. Male sex, blood pressure, NYHA-FC, low-voltage pattern (LVP), LV end-systolic volume (ESV), LV ejection fraction (EF), longitudinal ventricular function, pericardial effusion, NT-proBNP, hsTnT, dFLC and eGFR were associated with poor outcome, but only LV global longitudinal strain (HR 1.12; 95% CI 1.03 to 1.22; p=0.007), LV-ESV (HR 1.02; 95% CI 1.01 to 1.03; p=0.001), systolic pulmonary artery pressure (HR 0.98; 95% CI 0.96 to 0.99; p=0.027), NYHA-FC>II (HR 1.65; 95% CI 1.09 to 2.50; p=0.019), NT-proBNP (HR 1.57; 95% CI 1.17 to 2.11; p=0.003) and dFLC (HR 1.30; 95% CI 1.05 to 1.62; p=0.017) were independently predictive in the final Cox regression model.

Figure 2

Kaplan-Meier curve displaying overall survival of the complete study cohort. Median survival time was 1.47 years.

Table 3

Univariate and multivariable Cox regression analysis

Cut-off calculation and risk stratification

For all continuous parameters that were independent predictors of mortality in the final Cox regression model, their individual optimal cut-off values were computed (table 4). Corresponding Kaplan-Meier curves using the identified thresholds of NT-proBNP (5762 ng/L), global longitudinal strain (−8.9%), LV-ESV (43 mL), systolic pulmonary artery pressure (28 mm Hg) and dFLC (439 mg/L) for stratification are presented in figure 3A–F. Based on the aforementioned cut-off values a staging system incorporating global longitudinal strain and NT-proBNP was developed. The resulting stratification of our study population into three patient cohorts with low (both parameters below the threshold), intermediate (one of the two parameters above the threshold) and high risk (both parameters above the threshold) is shown in figure 4.

Figure 3

Kaplan-Meier curves displaying the survival probability according to parameters that were independent predictors in the multivariable Cox regression analysis. Stratification was done using individually calculated thresholds for each parameter. dFLC, difference of free light chains; ESV, end-systolic volume; LV-GLS, left ventricular-global longitudinal strain; NT-proBNP, N-terminal pro-brain natriuretic peptide; NYHA-FC, New York Heart Association functional class; sPAP, systolic pulmonary artery pressure.

Figure 4

Risk stratification based on a scoring system incorporating N-terminal pro-B-type natriuretic peptide (NT-proBNP) and global longitudinal strain (GLS). The thresholds were 5762 ng/L for NT-proBNP and −8.9% for GLS. A combination of both parameters separated groups with lower (both parameters below the threshold)=stage I, intermediate (one parameter above the threshold)=stage II and higher mortality risk (both parameters above the threshold)=stage III.

Table 4

Cut-off values for risk stratification with corresponding 95% bootstrap percentile of continuous variables

Discussion

This study aimed to investigate clinical, laboratory, electrocardiography and echocardiographic parameters for usefulness to predict outcome in AL amyloidosis with cardiac involvement proven by EMB.

Besides response to therapy, cardiac involvement is the main prognostic determinant in AL amyloidosis.27 Patients with cardiac involvement have a shorter survival than those without cardiac involvement, and mortality is even higher when congestive HF occurs.28 29 The search for prognostic markers therefore plays a crucial role, as they might aid in selecting therapeutic strategies, not only targeting the clonal disease but also regarding a timely and appropriate HF management.

Serum levels of NT-proBNP and cardiac troponins have emerged as important prognostic parameters and are therefore commonly used for risk stratification.16 18 In an unselected population of patients with AL amyloidosis the prognostic value of these biomarkers results from separating patients with a low (MAYO stage I) from those with a high probability of cardiac involvement (MAYO stage>I). Furthermore, they subdivide patients with suspected or proven cardiac involvement into a group with milder (MAYO stage II) and more advanced HF (MAYO stage III and IIIb). Beyond serological testing, functional cardiac imaging, particularly strain analysis, plays a fundamental role in risk stratification and prognostication in AL amyloidosis.30 However, the aforementioned parameters must be interpreted with caution because increased levels of biomarkers and pathological imaging findings can also be found in the context of other diseases and do not necessarily confirm or reflect cardiac involvement in AL amyloidosis.

The gold standard for the diagnosis of cardiac involvement is EMB.31 32 Due to its invasiveness, EMB is not common in most studies investigating survival criteria in AL amyloidosis in general23 33 34 and few studies have investigated the prognosis of patients with cardiac involvement proven by EMB in particular. These studies identified sex, Karnofsky index, diastolic blood pressure, NYHA-FC, NT-proBNP, GFR, mineralocorticoid receptor antagonist use, lack of β-blocker use, E-wave deceleration time, LV-EF, LV mass, LV wall thickness, LVP, intraventricular conduction delay, amyloid load, response to chemotherapy and SCT as markers associated with survival.22 35–37 Only LV-EF, LVP,35 NYHA-FC,36 NT-proBNP, GFR, response to chemotherapy and amyloid load37 were independent predictors of mortality in the multivariable model. However, cardiac deformation assessed by strain imaging had not been investigated for its prognostic value in patients with AL amyloidosis and EMB-proven cardiac involvement so far.

To close this knowledge gap, we screened the database of the amyloidosis centre at Heidelberg University Hospital, one of the largest of its kind worldwide. Among 1628 patients who had their initial visit from July 2006 to December 2017, we identified 216 patients with EMB proven cardiac involvement. To eliminate the effect of chemotherapy and to reduce the influence of comorbidities on our analysis, only therapy-naïve patients without a history of significant cardiac or renal disease at the time of EMB were further investigated. At the end of the screening process 174 patients met the criteria for inclusion in this study.

The median OS of the complete cohort was 1.47 years (95% CI 1.04 to 2.92 years) which is comparable with the literature on patients with AL amyloidosis.6 28 29 One hundred and fifteen patients (66%) had died after the 12 years of the investigation. The survival probability at 10 years from the time of first contact with the amyloidosis centre was 0.195 (95% CI 0.120 to 0.316).

Comparing survivors with non-survivors male sex, blood pressure, NYHA-FC, LVP, LV- ESV, LV-EF, longitudinal ventricular function, pericardial effusion, NT-proBNP, hsTNT, dFLC and eGFR were associated with poor outcome but only NYHA-FC >II, LV global longitudinal strain, LV-ESV, systolic pulmonary artery pressure, NT-proBNP and dFLC were independently predictive in multivariable analysis. Interestingly, from the time of EMB proven cardiac involvement, neither the serum level of cardiac hsTNT, as an indicator of myocardial damage and integral component of current amyloidosis staging systems, nor LV-EF, as a standard measure of cardiac systolic function and crucial criterion for HF classification, had independent prognostic significance in our study cohort. However, this does not mean that both parameters have no predictive value at all. Rather our variable selection was subject to some restrictions and the circumstances of missing values might have introduced additional uncertainty. Due to the number of predictors in relation to the effective sample size for multivariable Cox regression, the better approach of backward elimination was not feasible and so forward selection was performed for variable selection.

New treatment options have improved the survival of patients with AL amyloidosis over the past few decades.20 A basic requirement before starting therapy is optimal risk stratification. However, only patients in MAYO stage I had significant survival improvement over time, while the survival improvement in MAYO stage II and III/IIIb was marginal or nil, respectively.38 Risk stratification tools in patients with both, AL amyloidosis and cardiac involvement, should therefore provide a better predictive power. A first step to improve the original MAYO score was to implement a more appropriate cut-off level for NT-proBNP to identify a very poor-prognosis subgroup of patients in MAYO stage III.18 From the cardiologist’s point of view additional options exist, that could supplement a prognostic evaluation based purely on blood tests. In this context, cardiac imaging holds great potential as it provides direct and objective structural and functional read-outs. A focused analysis of myocardial performance may be superior to clinical classification systems like NYHA-FC, the assessment of which may vary substantially, or biomarkers such as dFLC that primarily reflects haematological disease burden.17

Global longitudinal strain could be a promising candidate to optimise the risk stratification of cardiac AL amyloidosis. Compared with other echocardiographic parameters that were independently predictive in our study, global longitudinal strain has lower intraobserver and interobserver variability than LV-ESV39 40 and, unlike systolic pulmonary artery pressure, should be measurable in most patients. Furthermore, recent studies have demonstrated that strain imaging can predict survival and therapy response in AL amyloidosis in general41–45 while our data further suggest that global longitudinal strain in combination with NT-proBNP may offer a useful survival risk stratification for AL amyloidosis patients with cardiac involvement in particular. However, a prognostic benefit for patients with cardiac AL amyloidosis needs to be prospectively validated and tested against other accepted staging systems in the field.

The inclusion criteria of our study resulted in a highly selected patient cohort especially since some comorbidities such as atrial fibrillation might have been directly related to AL amyloidosis. Furthermore, the rates of atrial fibrillation, additional kidney amyloidosis and the use of SCT were very low compared with a general population with AL amyloidosis. Therefore, our findings are not transferable to all AL amyloidosis patients, especially to those without cardiac involvement.

Limitations

Although the amyloidosis centre at Heidelberg University Hospital has a nationwide catchment area, this is a single-centre study. The retrospective nature of our study and the requirement for EMB proven cardiac involvement may contribute to a selection bias. We minimised any bias by including all patients with AL amyloidosis and cardiac involvement during the above-mentioned time interval. Nevertheless, we want to point out that the cohort examined in this study is not representative of the total population of patients diagnosed with AL amyloidosis. This is reflected, among other things, in a low rate of concomitant renal involvement or atrial fibrillation.

The onset and therefore the true duration of AL amyloidosis are difficult to assess because patients often show no clinical signs even at the beginning of cardiac involvement. The fact that all patients in our study underwent EMB indicates a high suspicion for an advanced stage of the disease in our cohort.

While the multivariable Cox regression analysis identified predictors for mortality, other important indicators may have been missed. Additionally, the variable selection might have suffered from the uncertainty due to the multiple imputation.

Conclusion

Among all patients with AL amyloidosis those with cardiac involvement represent a high-risk population with limited therapy options and poor prognosis. However, the present study suggests, that a relevant number of patients with cardiac AL amyloidosis can survive many years after the initial diagnosis. Therefore, accurate risk stratification is necessary to identify cardiac amyloidosis patients with favourable prognosis. Incorporation of modern imaging techniques into existing or newly developed scoring systems is a promising option that might enable the implementation of risk-adapted therapeutic strategies.

Data availability statement

Data are available on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study involves human participants and this study was carried out after approval by the ethics committee of the University of HeidelbergReference number: 123/2006. Participants gave informed consent to participate in the study before taking part.

Acknowledgments

For the publication fee, we acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding programme 'Open Access Publikationskosten' as well as by Heidelberg University.

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 @echobasics

  • Contributors Conception and design of data: MA and DM. Analysis and interpretation: MA, JB, JAV and SG. Drafting of the manuscript: MA. Criticall revision of the manuscript for important intellectual content: FADS, UH and SS. Final approval of the manuscript: HAK, NF and DM. DM is responsible for the overall content, accepts full responsibility for the finished work, had access to the data, and controlled the decision to publish.

  • Funding This study was funded by Deutsche Forschungsgemeinschaft, Heidelberg University.

  • Competing interests None declared.

  • 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.