Article Text

Original research
Prognostic value of C-reactive protein in adults with congenital heart disease
  1. Laurie W Geenen,
  2. Vivan J M Baggen,
  3. Annemien E van den Bosch,
  4. Jannet A Eindhoven,
  5. Robert M Kauling,
  6. Judith A A E Cuypers,
  7. Jolien W Roos-Hesselink,
  8. Eric Boersma
  1. Cardiology, Erasmus University Medical Center Rotterdam, Rotterdam, The Netherlands
  1. Correspondence to Dr Jolien W Roos-Hesselink, Cardiology, Erasmus MC, Rotterdam 3000 CA, The Netherlands; j.roos{at}erasmusmc.nl

Abstract

Background High-sensitivity C reactive protein (hs-CRP) has been associated with outcomes in adult congenital heart disease (ACHD). However, its prognostic value beyond N-terminal pro B type natriuretic peptide (NT-proBNP) or troponin T remains unknown. We studied the temporal evolution of hs-CRP, as well as the relation between hs-CRP and adverse clinical outcomes independent of NT-proBNP and troponin T in patients with ACHD.

Methods In this prospective cohort study, we enrolled 602 patients with ACHD (2011–2013) who underwent baseline and thereafter annual blood sampling during 4 years. Hs-CRP, hs-troponin T and NT-proBNP were measured. The primary endpoint was composed of death or heart failure (HF). Cox regression and Joint Modelling was used to relate 2log hs-CRP levels with the endpoint, with adjustment for baseline characteristics and (repeated) hs-troponin T and NT-proBNP measurements.

Results Hs-CRP was measured at baseline in 591 patients, median age 33 years, 58% men, 90% New York Heart Association I with an average of 4.3 measurements per patient. Median follow-up was 5.9 (IQR 5.3–6.3) years (99.2% complete) and 69 patients met the endpoint. Higher baseline hs-CRP was independently associated with higher risk of death or HF (HR 1.36, 95% CI 1.19 to 1.55). Hs-CRP increased over time prior to death or HF, and repeated hs-CRP measurements were associated with the endpoint, independent of repeated NT-proBNP and hs-troponin T (HR 1.54, 95% CI 1.24 to 1.98).

Conclusions Hs-CRP carries incremental prognostic value for the risk of death or HF, beyond NT-proBNP and hs-troponin T. Hs-CRP increased prior to the occurrence of HF or death, supporting the role of inflammation in the clinical deterioration of patients with ACHD.

  • complex congenital heart disease
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.

Introduction

Adults with congenital heart disease (ACHD) are a growing patient population with a high healthcare utilisation.1 2 Despite the increased survival prospects, patients with ACHD carry an increased risk of heart failure (HF), arrhythmias, reinterventions and early demise.3–5 Although mechanisms such as chronic pressure or volume overload that can lead to cardiac deterioration are well described in patients with ACHD,6 the potential role of chronic inflammation has been investigated to a lesser extent. A recent study by Opotowsky et al demonstrated that higher high-sensitivity C-reactive protein (hs-CRP) was associated with higher risk of cardiovascular events and mortality in patients with ACHD.7 These findings suggest that chronic inflammation may contribute to the pathophysiology of deterioration in ACHD, as hs-CRP is a sensitive marker of systemic inflammation.

Currently, it is unknown how hs-CRP relates to other biomarkers and echocardiographic measures in ACHD. Besides the great interest in the pathophysiological mechanism of hs-CRP, the additive prognostic value of hs-CRP beyond conventional biomarkers, such as N-terminal pro B type natriuretic peptide (NT-proBNP) and high-sensitivity troponin-T (hs-TnT), has not yet been assessed. It is therefore unclear how hs-CRP could enhance risk stratification. Moreover, it is unknown how hs-CRP evolves over time. Studying these temporal patterns may contribute to better understanding of hs-CRP release in patients with ACHD.

The aim of this study was to investigate the association of hs-CRP with other biomarkers and echocardiographic measurements, as well as its association with long-term cardiovascular outcomes in patients with ACHD. Moreover, we studied the temporal evolution of hs-CRP and the relation of these evolutions with the risk of cardiovascular events.

Methods

Patient and public involvement

Patients or the public were not involved in the design, or conduct, or reporting, or dissemination plans of our research.

Study design and population

This is a prospective observational single-centre cohort study, including 602 clinically stable patients with moderately or severely complex ACHD. Detailed aspects of the study protocol have been described previously.8 9 Patients were enrolled during a routine visit at the outpatient clinic (2011–2013). Exclusion criteria were age <18 years, mild ACHD, renal impairment (creatinine >200 μmol/L) or pregnancy. All patients underwent physical examination by a cardiologist, ECG, echocardiography and venous blood sampling at the day of enrolment. Thereafter, yearly outpatient visits were scheduled for 4 years and venous blood sampling was repeated. Patient management was according to the discretion of the treating physician, based on prevailing guidelines.10 11

The study was performed according to the principles outlined in the Declaration of Helsinki. All patients gave written informed consent for their participation in the study.

Laboratory testing

During baseline and follow-up study visits, peripheral venous blood sampling was performed. Blood samples were directly used, or stored at −80°C in the biobank until batch analysis was performed. NT-proBNP and estimated glomerular filtration rate (eGFR) were directly measured in fresh blood samples in our clinical chemistry laboratory. Serum hs-CRP and hs-TnT were measured batch wise from the same samples (first batches of 591 baseline measurements in 2012; second batches of 1943 follow-up samples in 2018) in thawed serum samples in the same laboratory using Roche immunoturbidimetric assays and electrochemiluminescence immunoassays, respectively (Roche Diagnostics, Basel, Switzerland). Samples had not been exposed to a prior freeze-thaw cycle. On average 4.3 measurements per patient were available. The lower limits of detection of hs-CRP, hs-TnT and NT-proBNP were 0.30 mg/L, 5 ng/L and 0.6 pmol/L, respectively. Lower limit of blank of hs-TnT was 3 ng/L. For analytical purposes, hs-CRP levels<0.30 mg/L were substituted with 0.15 mg/L, and hs-TnT levels<3 ng/L were substituted with 1.5 ng/L. Elevated levels of hs-CRP and NT-proBNP were defined as; >3.0 mg/L and >14 pmol/L, respectively.

Definition and assessment of endpoints

The primary endpoint was defined as all-cause mortality or HF (defined as HF requiring initiation/intensification of HF medication or requiring hospitalisation). The secondary endpoint was defined as any major adverse cardiovascular event (MACE) and was composed of; all-cause mortality, HF, hospitalisation for cardiac reasons (eg, endocarditis), arrhythmia (symptomatic and recorded, or requiring treatment), thromboembolic events or cardiac (re)interventions (surgical or percutaneous). Survival status was checked in the Municipal Population Register. Endpoint events were adjudicated by two investigators who were blinded for any study related information.

Statistical analysis

Sample size calculation was performed and has been previously been described.9 In the current study, the composite endpoint of death or HF was considered the primary endpoint because of its greater clinical relevance than any MACE.

Continuous data are presented as mean±SD or median (IQR). Trends between clinical characteristics and hs-CRP tertiles were tested with the χ2 Mantel-Haenszel test for trend or linear regression. Skewed biomarker distributions were transformed by taking logarithms to the base 2. The variability of hs-CRP and reference change value (RCV) were assessed in primary endpoint-free patients. The methods are specified in online supplemental file 1.

Supplemental material

Survival curves were derived using the Kaplan-Meier method and were stratified based on elevated NT-proBNP and hs-CRP levels. The log-rank test was used to compare survival curves. Cox regression analyses were performed to investigate the association between baseline hs-CRP and the endpoints. We present crude HRs, and HRs that are adjusted for NT-proBNP, hs-TnT and baseline characteristics including age, sex, congenital diagnosis, New York Heart Association (NYHA) class, cardiac medication and systemic ventricular function.

The temporal evolution of hs-CRP was described using linear mixed effect models (LME) and Cox models were used to assess the association between hs-CRP and the endpoints. The LME model is described in online supplemental file 2. Parameters of LME and Cox models were estimated jointly in Joint Models to relate biomarker trajectories with endpoints and avoid bias.12 The rational and assumptions of Joint Models are outlined in online supplemental file 3. Hs-CRP measurements that were taken after occurrence of endpoints were discarded in the analyses. We reported crude HRs, and HRs adjusted for repeatedly measured NT-proBNP and hs-TnT,13 and baseline characteristics. Subanalyses were performed to assess the association between hs-CRP and endpoints in patients with elevated baseline NT-proBNP. The C-index was calculated to assess the discriminative value of hs-CRP.

Missing baseline covariates were handled by single imputation. SPSS (V.24) and R statistical software (V.3.6.1, packages mice, Survival, nlme, JMbayes) were used. A two-tailed p<0.05 was considered statistically significant.

Results

Baseline characteristics

In 591 patients (98.2%), hs-CRP was measured at baseline and was therefore included in the current analysis. Median age was 33 (IQR 25–41) years, 58% of the patients were men and 90% were NYHA class I (table 1). The most prevalent congenital diagnoses were: Tetralogy of Fallot (n=174, 29%), aortic stenosis (n=136, 23%) and aortic coarctation (n=110, 19%). Median hs-CRP level was 1.50 (IQR 0.60–3.50) mg/L. Hs-CRP was below the limit of detection in 44 patients (7.4%), and 164 patients (28%) had elevated hs-CRP. Higher hs-CRP was associated with older age, female sex, higher body mass index (BMI), cardiac medication use and NYHA class. Highest hs-CRP levels were found in patients with a functionally univentricular heart and those with pulmonary arterial hypertension (PAH) (figure 1). None of the echocardiographic measurements were related to hs-CRP, except for measures of left ventricular diastolic function (E/A ratio, E’ wave and E/E’ ratio). Hs-CRP showed significant, but weak, correlations with NT-proBNP (Spearman r=0.28, p<0.001) and eGFR (Spearman r=−0.14, p=0.001, but not with hs-TnT (Spearman r=0.04, p=0.323).

Table 1

Baseline characteristics of the study cohort for all patients and categorised according to the tertile distribution of hs-CRP

Figure 1

Hs-CRP levels according to congenital diagnosis. The median hs-CRP level in each diagnosis group is indicated by the horizontal red line together with the 25th and 75th percentile presented by the blue bars. P value is given for comparison of hs-CRP levels between diagnoses (for the analysis univentricular heart was combined with Fontan, and Rastelli/REV with PAH). ccTGA, congenitally corrected transposition of the great arteries; hs-CRP, high-sensitivity C reactive protein; PAH, pulmonary arterial hypertension; TGA-ASO, transposition of the great arteries corrected by the arterial switch operation; TGA-Mustard, transposition of the great arteries corrected by the Mustard procedure.

Follow-up and associations between baseline hs-CRP and endpoints

Follow-up data were complete in 586 patients (99.2%). Median follow-up was 5.9 (IQR 5.3–6.3) years. Death or HF occurred in 69 patients (11.8%) and MACE in 228 patients (38.9%). Occurrence of individual endpoint events was as follows: 25 deaths, 59 HF, 127 arrhythmias, 180 hospitalisations (of which n=13 endocarditis), 137 (re)interventions and 29 thromboembolic events. Patients in the first hs-CRP tertile (<0.80 mg/L) had a significant better HF-free and MACE-free survival than patients in higher tertiles (online supplemental file 4). Analysed continuously, hs-CRP was associated with death or HF (HR 1.36, 95% CI 1.19 to 1.55, C-index 0.66), also after adjustment for baseline NT-proBNP, hs-TnT and baseline clinical characteristics (table 2). Hs-CRP was also associated with death and HF when analysed as separate events (online supplemental material 5).

Table 2

Association between baseline hs-CRP levels and the endpoints, adjusted for clinical characteristics, NT-proBNP and hs-TnT

Patients with both elevated baseline hs-CRP and NT-proBNP, had the worst HF-free survival (figure 2).

Figure 2

Heart failure-free and MACE-free survival according to normal or elevated baseline levels of hs-CRP and NT-proBNP. Elevated hs-CRP was defined as >3 mg/L and elevated NT-proBNP as >14 pmol/L. hs-CRP, high-sensitivity C reactive protein; MACE, major adverse cardiovascular event; NT-proBNP; N-terminal pro B type natriuretic peptide.

Associations between repeated hs-CRP and endpoints

In total, 2421 hs-CRP measurements were available prior to the primary endpoint. Baseline hs-CRP was higher in patients with the primary endpoint (2.53 mg/L) than those who remained endpoint-free (1.29 mg/L)(p<0.001). Hs-CRP further increased in anticipation to death or incident HF, whereas endpoint-free patients showed stable hs-CRP levels throughout follow-up (figure 3). The average 4-year hs-CRP increase in patients with the primary endpoint was 3.09 mg/L, vs 0.07 mg/L in stable patients. In relation to MACE, no difference was found in hs-CRP evolution.

Figure 3

Average evolution of hs-CRP levels in patients with ACHD with and without the endpoint. All measurements taken after the endpoints were discarded. Time point zero is denoted as the time of study inclusion. The average evolution is accompanied by the 95% prediction band (thin lines). ACHD, adult congenital heart disease; hs-CRP, high-sensitivity C reactive protein; MACE, major adverse cardiovascular event.

Any time during follow-up, each twofold higher hs-CRP was associated with a 2.14 times higher hazard of mortality or HF (C-index 0.82) (table 3). This remained significant after adjustment for repeatedly measured NT-proBNP and hs-TnT, and baseline clinical characteristics. Repeated hs-CRP showed a significant association with MACE, but not when adjusted for baseline characteristics, or other repeatedly measured biomarkers. Associations between repeated hs-CRP and individual endpoints are shown in online supplemental file 5.

Table 3

Association between repeated hs-CRP levels and endpoints, adjusted for clinical characteristics and (repeated) NT-proBNP and hs-TnT

As posthoc sensitivity analysis, we repeated the analysis after excluding patient with baseline hs-CRP levels>10 mg/L (n=34). This did not result in qualitative meaningful differences in strengths of the associations. Similarly, analysis excluding women who became pregnant during follow-up (n=31) did not lead to different results. Posthoc analysis stratifying for sex showed that repeated hs-CRP was more strongly associated with the risk of death or HF in men than in women (HR 3.95 (95% CI 2.39 to 8.25) vs 1.55 (95% CI 1.20 to 2.03), p=0.005). This difference was not present when adjusted for NT-proBNP (HR 1.81 (95% CI 1.24 to 2.78) vs 1.51 (95% CI 1.17 to 1.99), p=0.157). Subanalysis restricted to patients with elevated baseline NT-proBNP demonstrated that the association between repeatedly measured hs-CRP and the risk of death or HF, is also present in these patients (online supplemental file 6).

Biological variation of hs-CRP

Based on repeated hs-CRP measurements in primary endpoint-free patients, the within-subject and between-subject variation of hs-CRP were 52% and 113%, respectively. The index of individuality was 0.48, indicating that it may be preferred to have a patient-based reference value based on previously measured hs-CRP levels of the individual patient. The RCV was 149% (lognormal RCV limits −75% to 304%).

Discussion

This study demonstrated that in patients with ACHD higher baseline hs-CRP is associated with a higher risk of death of HF, independent of NT-proBNP and hs-TnT. Furthermore, annually repeated hs-CRP measured up to 4 years, showed that on average hs-CRP levels steadily increased prior to the occurrence of death or HF and remained stable in HF-free patients. As repeated hs-CRP measurements provided prognostic information beyond repeated NT-proBNP and hs-TnT measurements, hs-CRP is a valuable potential biomarker for risk stratification in patients with ACHD.

Comparison with previous studies

In 2018, Opotowsky et al described hs-CRP in a prospective cohort consisting of 707 patients with ACHD.7 The proportion of patients with elevated hs-CRP (>3 mg/L) in our study was comparable to their study (28% vs 25%), as were study settings on most aspects, except for our cohort not including mild ACHD. This might explain why we found slightly higher hs-CRP levels (median 1.50 vs 1.22 mg/L), despite that our patients were younger, more often NYHA class I and had a lower BMI.

Whereas in our study baseline hs-CRP primarily showed an independent association with all-cause mortality and HF, Opotowsky et al found associations with other adverse cardiac events such as arrhythmia and thromboembolic events. An explanation could be the lower event rate and longer follow-up in our cohort, creating a longer time interval between baseline hs-CRP and event occurrence, probably diluting associations. The fact that repeatedly measured hs-CRP was associated with arrhythmia, reintervention and hospitalisation in our study may support this, as hs-CRP will be measured more closely to the event and therefore better reflects a patients’ disease status.

Scognamiglio et al investigated serial CRP measurements retrospectively in 225 adults with PAH due to congenital heart disease. Higher CRP correlated weakly with higher BNP, similar to our study, and was associated a higher mortality risk. While they excluded CRP measurements in patients with signs of active infection, elevated CRP levels (>10 mg/L) were still measured in 54% of the patients at least once during follow-up.14 Four out of nine patients with PAH in our study had baseline hs-CRP >10 mg/L; one patient had a cold without fever; in the others, no signs of infection were present. Therefore, elevated hs-CRP seem not restricted to acute infections but may be triggered by their disease state linked to PAH.

Role of inflammation in ACHD

Hs-CRP is an acute phase protein and is commonly used blood test to assess presence of acute inflammation. Improved hs-CRP assays now allow for assessment of chronic, low inflammatory levels in not eminently inflammatory conditions, like ACHD. Findings of our study support evidence regarding chronic inflammation accompanying cardiac deterioration and mortality in patients with ACHD.7 14 15 The question remains whether inflammation induces HF or is rather an ‘innocent bystander’ and thereon reflects disease severity. Hs-CRP was not associated with measures of right or left ventricular systolic function in our study, while hs-CRP was associated with measures of diastolic dysfunction. This could be explained by age as contributing factor to both increasing hs-CRP and diastolic dysfunction.16 Nevertheless, other biomarkers that are considered more linked to the HF pathophysiology, such as NT-proBNP and hs-TnT, seem to relate to echocardiographic measurements of cardiac function more closely.9 17 18 In addition, hs-CRP did not correlate with hs-TnT and only weakly with NT-proBNP. Therefore, causal involvement of hs-CRP in the pathophysiology of cardiac deterioration seems less plausible.

Data on whether inflammation are a bystander or a mediator in cardiovascular diseases, in general, remain controversial.19 20 In pulmonary arterial smooth muscle cells, known for their role in PAH, it was shown that impediment of inflammation by atorvastatin, reduces hs-CRP-induced expression of inflammatory pathways,21 suggesting causality. Unfortunately, no specific data on ACHD exist. Longitudinal hs-CRP patterns now showed that inflammation increased in patients with ACHD prior to deceasing or HF. Similar patterns have been observed in patients with chronic HF.22 This at least suggests that progression of HF is associated with increasing inflammatory levels and that stable inflammatory levels are indicative of a stable cardiac status.

Value of hs-CRP as prognostic biomarker in ACHD

There may be several explanations for the independent prognostic value of hs-CRP. First, hs-CRP may reflect other processes of deterioration in ACHD that are not captured by clinical characteristics and biomarkers measured in our study, this makes causal inferences on inflammation and cardiac deterioration difficult, as this requires measurement and adjustment for all other potential confounders. Residual confounding could distort the association between hs-CRP and endpoints, resulting in false conclusions about the relation between inflammation and cardiac deterioration. Nevertheless, hs-CRP can then still be valuable for risk stratification. As hs-CRP is a non-specific inflammatory marker, hs-CRP may identify vulnerable patients in general, and therefore have prognostic value.

The biological variation of hs-CRP found in this study was considerable, which is in line with previous studies.23 24 However, both within-subject and between-subject variation were higher and thus the corresponding RCV. This can most likely be explained by differences in study design: less frequent sampling, less hs-CRP measurements per patient and a different study population.

Clinical implications

Hs-CRP is able to provide prognostic value not captured by other prognosticators. In the absence of another likely cause of inflammation, an elevated hs-CRP (>3 mg/L) may help to detect patients at higher risk of HF or mortality. Hs-CRP may specifically be useful to detect vulnerable patients in copresence of an elevated NT-proBNP. Measuring hs-CRP should therefore be considered in addition to the current practice of NT-proBNP in patients with ACHD. Repeated NT-proBNP can be combined with repeated hs-CRP to further enhance risk stratification. Increasing hs-CRP levels over time can identify high-risk patients who require closer monitoring, while stable levels can identify low-risk patients who can be reassured and less frequently monitored. Nonetheless, the substantial biological variation and high RCV should be kept in mind when interpreting repeated hs-CRP measurements.

If hs-CRP is induced by low-grade inflammation of an extracardiac origin, that subsequently leads to cardiac deterioration, inflammation in general may have substantial negative influences on the prognosis of patients with ACHD. For clinical practice, this would indicate that attention should be paid to inflammatory diseases in ACHD, as this might trigger HF.

Last, inflammation may be a potential therapeutic target; however, further research is needed to elucidate the inflammatory mechanism in ACHD.

Limitations

We were not able to definitely rule out concomitant inflammatory processes which may have led to increasing hs-CRP levels. The association found in our study might therefore be interfered by occurrence of acute infections in some patients. However, we included stable patients with ACHD in who hs-CRP was measured during regular study visits, making the interference of infections with their prognosis and biomarker measurement very unlikely. Also, documented baseline chronic infections and inflammatory conditions were rare; 10 patients carried a rheumatologic diagnosis, 2 patients inflammatory bowel disease, 1 latent tuberculosis and 1 cystic fibrosis. However, we had no data on depression, while depression has been associated with inflammation and poorer prognosis in ACHD.25

We could not adjust for oral contraceptive or hormonal replacement use. This may have biased the associations, as hormonal replacement therapy is known to be associated with higher hs-CRP levels26 and it can be questioned whether oral contraception use could have increased the risk of certain cardiovascular events in women. Peak VO2 and peak heart rate have been associated with increased hs-CRP levels,7 as we did not have exercise data, we could not assess the prognostic value of hs-CRP beyond these variables.

Conclusion

In clinically stable patients with ACHD, an elevated hs-CRP in combination with an elevated NT-proBNP is associated with an increased risk of death or HF. Hs-CRP levels increased over time prior to the occurrence of HF or death, and serially repeated hs-CRP measurements provided additional prognostic information over repeated NT-proBNP measurements. This study shows that hs-CRP is a potential biomarker for risk stratification in patients with ACHD.

Key messages

What is already known on this subject?

  • In clinically stable patients with adult congenital heart disease (ACHD), N-terminal pro B type natriuretic peptide (NT-proBNP) can be used to identify patients with a poor prognosis. Recently, C-reactive protein (CRP) has been associated with adverse cardiac outcomes in ACHD; however, its additive prognostic value beyond NT-proBNP remains unknown.

What might this study add?

  • CRP can help to further discriminate between patients at low or high risk of death or heart failure. Higher CRP is associated with an increased risk of death or heart failure, independent of clinical characteristics, NT-proBNP and troponin T. Moreover, annual CRP measurements show that CRP levels increase prior to the occurrence of death or heart failure and repeated CRP measurements provide prognostic information beyond repeated NT-proBNP and troponin T measurements.

How might this impact on clinical practice?

  • In the routine clinical assessment of patients with ACHD, an annual CRP measurement should be considered in addition to an annual NT-proBNP measurement to improve risk stratification. Further research should investigate how low-graded inflammation is related to worsening heart failure in patients with ACHD.

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

  • Contributors All of the authors contributed to one or more of the following: conception of design of the work (VJMB, AEvdB, JAE, JAAEC, JWR-H, EB), data collection (LWG, VJMB, AEvdB, JAE, RMK, JAAEC, JWR-H), data analysis and interpretation (LWG, JWR-H, EB), drafting the manuscript (LWG, VJMB, JWR-H, EB) and critical revision of the manuscript (all authors). All authors provided final approval of the manuscript. LWG and EB are responsible for the overall content as guarantors.

  • Funding Roche diagnostics provided C-reactive protein assays free of charge to support this investigator initiated study.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval The study protocol was approved by the local medical ethics committee.

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

  • Data availability statement Data are available on reasonable request.

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