Article Text

Extended report
Risk of skin and soft tissue infections (including shingles) in patients exposed to anti-tumour necrosis factor therapy: results from the British Society for Rheumatology Biologics Register
  1. James B Galloway1,
  2. Louise K Mercer1,
  3. Alison Moseley1,2,
  4. William G Dixon1,
  5. Andrew P Ustianowski3,
  6. Matthew Helbert4,
  7. Kath D Watson1,
  8. Mark Lunt1,
  9. Kimme L Hyrich1,
  10. Deborah PM Symmons1
  1. 1Arthritis Research UK Epidemiology Unit, Manchester Academic Health Sciences Centre, The University of Manchester, Manchester, UK
  2. 2University of Birmingham, Birmingham, UK
  3. 3Department of Infectious Diseases, North Manchester General Hospital, Manchester, UK
  4. 4Department of Immunology, Central Manchester NHS Foundation Trust, Manchester, UK
  1. Correspondence to Dr Kimme Hyrich, Arthritis Research UK Epidemiology Unit, Manchester Academic Health Sciences Centre, Stopford Building, Oxford Road, The University of Manchester, Manchester, M13 9PT, UK; kimme.hyrich{at}manchester.ac.uk.

Abstract

Introduction Anti-tumour necrosis factor (TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). In 2001, BSRBR was established to evaluate the safety of these agents. This paper addresses the safety of anti-TNF therapy in RA with specific reference to serious skin and soft tissue infections (SSSI) and shingles.

Methods A cohort of anti-TNF-treated patients was recruited alongside a comparator group with active RA treated with non-biological disease-modifying antirheumatic drugs (nbDMARD). 11 881 anti-TNF and 3673 nbDMARD patients were analysed. Follow-up was by 6-monthly questionnaires to patients and clinicians. Analyses considered SSSI and shingles separately. Incidence rates (IR) were calculated and then compared using survival analyses.

Results The crude IR for SSSI were: anti-TNF 1.6/100 patient-years (95% CI 1.4 to 1.8); nbDMARD 0.7/100 patient-years (95% CI 0.5 to 1.0) and shingles: anti-TNF 1.6/100 patient-years (95% CI 1.3 to 2.0); nbDMARD 0.8/100 patient-years (95% CI 0.6 to 1.1). Adjusted HR were SSSI 1.4 (95% CI 0.9 to 2.4), shingles 1.8 (95% CI 1.2 to 2.8). For SSSI, no significant differences were seen between anti-TNF agents. For shingles, the lowest risk was observed for adalimumab (adjusted HR vs nbDMARD) 1.5 (95% CI 1.1 to 2.0) and highest for infliximab (HR 2.2; 95% CI 1.4 to 3.4)).

Conclusion A significantly increased risk of shingles was observed in the anti-TNF-treated cohort. The risk of SSSI tended towards being greater with anti-TNF treatment but was not statistically significant. As with any observational dataset cause and effect cannot be established with certainty as residual confounding may remain. This finding would support the evaluation of zoster vaccination in this population.

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/3.0/ and http://creativecommons.org/licenses/by-nc/3.0/legalcode

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.

Compared with the general population, skin and soft tissue infections occur around three times more frequently in patients with rheumatoid arthritis (RA).1 This is due to a combination of factors, including both the immunosuppressive treatments as well as a complication of the disease itself.

It is now over a decade since the introduction of a new class of biological disease-modifying antirheumatic therapies to the treatment armamentarium for RA: the anti-tumour necrosis factor (TNF) agents. While these treatments have dramatically changed our ability to control the disease, they have also been associated with an increased risk of infection, especially early on after commencing the therapy.2,,5

In addition, there is evidence suggesting that this risk of infection is not constant across anatomical sites.6 Very few data exist specifically examining the risk of skin and soft tissue infections. However, a subgroup of skin infections caused by herpes zoster (shingles), has been studied by other European registries.7 ,8 Data from a German biologics register (RABBIT) identified an increased risk of shingles in patients treated with monoclonal antibodies against TNF (eg, infliximab, adalimumab but not etanercept, ETNa recombinant TNF receptor fusion protein).9

The primary aim of this research was to explore the whether anti-TNF therapy increases the risk of skin and soft tissue infections (including shingles) above that experienced by RA patients treated with more traditional immunosuppressive regimens. The secondary aims were to compare the risks in patients exposed to different anti-TNF agents and to examine whether the risk changed with the duration of exposure.

Methods

Ethics approval for this study was obtained from the Multicentre Research Ethics Committee for the northwest of England.

The British Society for Rheumatology Biologics Register (BSRBR) is a prospective observational cohort study. Details of this study including methods have been published previously.10 It was established in 2001, and is currently ongoing, with ethical approval to follow patients until at least 2013. It was initiated alongside national recommendations in the UK that all RA patients prescribed anti-TNF therapy should be enrolled on the register.11 Only etanercept, infliximab and adalimumab are considered in this study as the other anti-TNF therapies have only recently been introduced into the UK market. Recruitment to the infliximab and etanercept cohorts began from the start of the study, while recruitment to the adalimumab cohort began in 2003.

A comparison cohort of patients with active RA receiving non-biological disease-modifying antirheumatic drugs (nbDMARD) was recruited in parallel. Active RA was defined as having a 28 joint count disease activity score (DAS28)12 greater than 4.2. All patients in both cohorts were biological naive at entry. At the time of this analysis, the BSRBR data had over 90% power to detect a doubling in the rate of both skin infections and shingles in the anti-TNF cohort compared with the comparison cohort.

Patients prescribed biological agents were recruited from across the UK (over 250 hospitals), whereas the comparison cohort was recruited from 29 centres, reflecting a combination of secondary and tertiary care rheumatology centres distributed across the UK.

Baseline assessment

All patients in this study had a physician diagnosis of RA. Baseline information included demographics, disease duration, a measure of RA disease activity using the DAS28, a measure of self-reported physical function (the health assessment questionnaire; HAQ),13 baseline steroid use, smoking history and baseline co-morbidity.

Follow-up

Follow-up information was collected by postal questionnaire every 6 months from the treating rheumatology team as well as the patients. Patients were also provided with a diary to record details of all hospital attendances and new prescriptions. All patients were flagged with the UK National Health Service Information Centre that reports to the BSRBR on deaths of patients within the study. Adverse events were coded using the Medical Dictionary for Regulatory Activities terminology. All follow-up was censored at 31 December 2009. Follow-up was limited to the first 3 years per patient as information was available from all three sources of reporting (physician, patient and UK National Health Service Information Centre) during this time.

Case definition and verification

Serious skin and soft tissue infections (SSSI) were analysed separately from shingles. Skin and soft tissue infections were only included if they were serious in nature (defined as resulting in hospitalisation, requiring intravenous antibiotics or causing death). Postoperative infections (ie, any infection occurring within 30 days of surgery) were excluded. An additional proforma requesting further details was sent to the treating rheumatologist in all cases of reported serious infection. As many cases of shingles would not satisfy the definition of serious, all shingles cases were analysed irrespective of the source of reporting or the seriousness of the event. JBG reviewed all soft tissue infections (SSSI) reported to the BSRBR to ensure correct coding. JBG and AM reviewed and additionally categorised shingles cases independently. ‘Severe’ shingles was defined as being multidermatomal, requiring intravenous antiviral agents, or being a primary reason for hospital admission.

Statistical methods

Baseline characteristics are presented with p values to assess significant differences between cohorts. Baseline values for the individual anti-TNF cohorts refer to the first biological agent a patient received. For dichotomous variables, a Pearson's χ2 was used. For comparison of continuous variables Wilcoxon rank sum tests were used allowing for non-parametric data. When more than two groups were present the Kruskal–Wallis test was used.

Crude incidence rates were calculated as the number of events per 100 patient-years of follow-up. Patients contributed follow-up in the nbDMARD cohort from the date of registration until their first event, 3 years after registration, 31 December 2009, last returned follow-up or death, whichever came first. Follow-up was censored at the point of switching for those patients who subsequently started a biological agent. In the anti-TNF cohort, patients contributed follow-up time from the first dose of the drug until the first skin infection, 90 days after their first missed dose of anti-TNF agent, 3 years after their first dose, 31 December 2009, last returned follow-up or death, whichever came first. Patients on anti-TNF who switched anti-TNF agents remained under follow-up and contributed follow-up time to their second or subsequent anti-TNF drug. Patients who switched to a non-anti-TNF biological agent had their follow-up censored at the point of switching.

Survival modelling with Cox proportional hazards was used to compare event rates between the groups. An overall HR is presented for the entire follow-up period as well as by 6-month time windows to explore the variation of risk with time. A sensitivity analysis was performed for shingles that compared the rates of serious shingles between the two treatment groups.

Potential confounders were identified before the analysis: age, gender, disease duration, disease severity (DAS28), disability (HAQ score), baseline steroid exposure, smoking status, relevant co-morbidity (diabetes, chronic obstructive pulmonary disease) and year of entry into the study. Adjusting for potential confounders was performed using a propensity model. This involved generating a score for every patient, which reflected the likelihood of them being in the anti-TNF cohort or the nbDMARD cohort based on their baseline information (including all the identified potential confounders). This score was then used as a weight in the survival model. The use of this method in observational studies has been described in detail previously.14

Missing baseline data were replaced using multiple imputations with 20 cycles. Patterns of missing follow-up data were compared between cohorts. Different numbers of patients dropped out of follow-up in the nbDMARD and anti-TNF cohorts due to switching therapy (switching from nbDMARD to anti-TNF or other biological agent, or switching from anti-TNF to an alternative biological agent). A sensitivity analysis was performed to investigate the extent of bias that may be present due to this differential dropout using a marginal structural model, re-weighting the propensity model every 6 months for the likelihood of dropping out due to switching.

Further details of the statistical methods are available in the supplementary attachment (available online only). All analyses were performed using Stata 10.1.

Results

Participants

This analysis included 15 554 patients (11 881 anti-TNF and 3673 nbDMARD-treated patients). Figure 1 shows the numbers of patients continuing under follow-up at each time point. A breakdown of the reasons for dropping out from follow-up is shown in the supplementary attachment (available online only).

Figure 1

Flowchart showing follow-up of patients. Flowchart showing cumulative numbers of patients in follow-up at each time point. Details of the reasons for dropout are presented in the supplementary data (available online only). DMARD, disease-modifying antirheumatic drug; TNF, tumour necrosis factor.

Baseline characteristics

Anti-TNF-treated patients were younger and a higher percentage was women compared with nbDMARD (table 1). The anti-TNF-treated patients had more severe disease. Differences between patients treated with the three anti-TNF agents were less marked.

Table 1

Baseline characteristics of study population

Event rates and hazard estimates comparing nbDMARD with anti-TNF: SSSI

Three hundred and nine SSSI were reported to the registry (269 anti-TNF; 39 nbDMARD). The majority of these were cellulitis (227/309; 73%) with the remainder being coded as abscesses or infected ulcers. Four cases of SSSI were necrotising fasciitis (all in the anti-TNF cohort). Culture information was available for 130/269 (48%) in the anti-TNF cohort and 15/39 (38%) in the nbDMARD cohort. Staphylococci were the most frequently reported organisms in both cohorts (see table 2). There were proportionally more cases of pseudomonal infection in the anti-TNF cohort. The median hospital stay for SSSI was 5.5 days in the anti-TNF-treated cohort (IQR 3–10) and 5 days in the nbDMARD-treated cohort (IQR 4–15) (p value for χ2=0.459). The incidence rate for SSSI was 1.6/100 patient-years in the anti-TNF cohort, compared with 0.7/100 patient-years in the nbDMARD cohort (table 3). The unadjusted HR for SSSI for anti-TNF was 2.2 (95% CI 1.5 to 3.0).

Table 2

Organisms responsible for skin and soft tissue infections

Table 3

Incidence and HR of SSSI

Univariable adjusted results are available in supplementary table 1 (available online only). The following variables altered the unadjusted result by more than 10%: age, disability (HAQ score) and year of entry into the study. Year of entry into the study reflects the fact that patients recruited earlier in the study had more severe disease (higher DAS28 score, HAQ score, disease duration). The fully adjusted HR using the propensity model for SSSI was 1.3 (95% CI 0.8 to 2.2).

Comparing the individual anti-TNF agents: SSSI

Using the nbDMARD cohort as the comparison, the lowest risk was seen in the adalimumab cohort (HR 1.1; 95% CI 0.6 to 2.1). The risks observed for etanercept and infliximab were identical (HR 1.5; 95% CI 0.9 to 2.5). Compared with etanercept, the HR for SSSI for the monoclonal antibodies combined was not significantly different (p=0.306). The differences did not reach statistical significance when adalimumab was compared with either etanercept or infliximab.

Event rates and hazard estimates comparing nbDMARD with anti-TNF: shingles

There were 275 cases of shingles in the anti-TNF cohort (incidence 1.6/100 patient-years) and 45 in the nbDMARD cohort (incidence 0.8/100 patient-years; table 4). Twenty-one cases of shingles were classified as severe (anti-TNF 20 (7%), nbDMARD 1 (2%), p value for χ2=0.206). There were five cases of multidermatomal shingles, and six cases of ophthalmic shingles in the anti-TNF cohort. The one severe case in the nbDMARD cohort was ophthalmic. Univariable adjusted results are available in supplementary table 1 (available online only). The propensity adjusted HR for shingles remained significantly elevated at 1.7 (95% CI 1.1 to 2.7).

Table 4

Incidence and HR of shingles

Comparing the individual anti-TNF agents: shingles

The adalimumab cohort had the lowest risk of shingles compared with nbDMARD (HR 1.5; 95% CI 0.9 to 2.4) followed by etanercept (HR 1.7; 95% CI 1.0 to 2.7) and then infliximab (HR 2.2; 95% CI 1.4 to 3.4). No significant difference was apparent when comparing the rates of shingles with etanercept with the monoclonal antibodies combined. However, within the monoclonal antibody class, the risk of shingles was significantly higher with infliximab when compared with adalimumab (HR 1.5; 95% CI 1.1 to 2.0). Excluding patients who switched anti-TNF treatment during follow-up did not alter these findings.

Event rate over time: SSSI and shingles

The hazard for both SSSI and shingles varied over time, with a higher risk noted early on in therapy (figure 2). An early peak in hazard was more marked for SSSI, in which in the first 6 months of therapy 71 events were reported in the anti-TNF cohort (IR 5.0/100 patient-years; 95% CI 3.9 to 6.4) compared with three in the nbDMARD cohort (IR 0.9/100 patient-years; 95% CI 0.2 to 2.6), resulting in an adjusted HR of 3.5 (95% CI 0.6 to 19.2). Review of the cases occurring within this early time period did not reveal any unusual patterns. The type of infection was predominantly cellulitis.

Figure 2

Hazard for soft tissue infections (SSI) and shingles over time in anti-tumour necrosis factor (TNF)-treated cohort. Demonstration of the time-varying risk of events in the anti-TNF cohort. Event rates were higher earlier on during treatment in both cohorts.

Differential dropout during follow-up

A significantly greater proportion of the nbDMARD cohort dropped out from follow-up because of switching to biological therapy (12.9% vs 0.7%). The patients who switch inevitably comprise a cohort with higher disease activity. As disease severity and the associated requirement for corticosteroid therapy are independent risk factors for infections, this will result in an effective depletion of susceptible individuals from the nbDMARD cohort. The HR for SSSI for anti-TNF adjusting for dropout using the marginal structural model was 0.9 (95% CI 0.5 to 1.4). For shingles, the hazard was influenced less by the adjustment for dropout (HR 1.5; 95% CI 1.0 to 2.4).

Discussion

To date, only limited information on the site-specific risk of infection with anti-TNF therapy has been available. This paper examined skin infections in detail and demonstrated a non-significant increase in SSSI, and a significantly increased risk of shingles.

The absolute risk of SSSI in the anti-TNF cohort was more than double that seen in the nbDMARD cohort. This finding should encourage clinicians to be vigilant for skin infections in patients starting anti-TNF therapy. It is likely that some of the additional risk of SSSI seen in the anti-TNF cohort is attributable to other risk factors. This is demonstrated both by the reduction in the HR from 2.2 to 1.4 after adjustment for confounders and the disappearance of any association after adjusting for differential dropout.

There was a higher proportion of pseudomonal infections in the anti-TNF cohort. While the clinical significance of skin culture results is not always clearcut (some of the reported pathogens could represent commensals) the difference is an interesting observation and suggests that TNF inhibition may modify the risk of SSSI. With this in mind, it is important to consider that while the reported HR for SSSI with anti-TNF was not statistically significant at 1.4, the upper limit of the 95% CI was 2.4. Therefore, it is important to acknowledge that a clinically relevant increased risk has not been excluded.

The increased risk of shingles with anti-TNF is important because shingles is associated with substantial morbidity.15 In BSRBR the incidence of shingles in the anti-TNF cohort was almost 2% per year. Adjusting for confounders or differential dropout made only a minor difference in the HR estimates, suggesting that the association between anti-TNF and shingles is less likely to reflect selection bias. The estimated number needed to treat for 1 year with anti-TNF to see one additional case of shingles was 128. In the USA, the Centers for Disease Control and Prevention (CDC) has recommended that all adults over the age of 60 years be offered vaccination against shingles (Zostavax) irrespective of previous exposure status.16 The vaccine is a live therapy and should be used with caution in people who are immunosuppressed. The CDC recommends that, while low dose corticosteroids or methotrexate are not contraindications, the vaccine should not be administered to individuals already established on anti-TNF therapy. These findings should encourage clinicians to consider vaccination before exposure to anti-TNF therapy. It is worth noting that current estimates of any vaccine uptake among RA patients are low.17 Further study into the utility of vaccination in RA cohorts is clearly warranted.

A hypothesis that etanercept carried a lower risk of shingles compared with the monoclonal antibodies (infliximab and adalimumab) arose from research published by the German biologics registry.9 In this study we observed an increased risk of shingles with all three drugs. The hazard estimates from the BSRBR lie within the CI limits reported by the German registry, and the finding of a significant association for all drugs may simply reflect the greater power of the BSRBR dataset.

The risk of SSSI and shingles was greatest early on after initiating anti-TNF therapy. This pattern has been noted previously and a number of explanations have been hypothesised.18 It is likely that there are a number of factors contributing to the pattern. Undoubtedly, a ‘healthy user’ effect is present, but in addition, there is also likely to be changes in other time-varying confounders such as steroid exposure.

The strengths of this study are the large sample size combined with a simultaneously recruited comparator cohort allowing for robust comparisons of rare outcomes that could not be studied in the context of clinical trials. The multiple reporting methods incorporating both information from patients and clinicians further strengthens the validity of the results.

There are also limitations to acknowledge. Although these data allow us to say that patients chosen to start on anti-TNF therapy within the BSRBR are at higher risk of shingles and SSSI than the comparator cohort, we cannot say with certainty that this increased risk is directly attributable to the anti-TNF treatment. There remains the possibility that residual confounding exists. For example, the lack of high quality information about steroid use during follow-up precluded the analysis of steroids as a time-varying covariate. Diagnostic suspicion bias may also have played a role as clinicians may have had a lower threshold for admitting patients on anti-TNF therapy. However, if this was the case we might have expected to see shorter hospital admission durations in the anti-TNF arm, and this was not the case. Conversely, we may have underestimated the risk of SSSI or shingles with anti-TNF as it is possible that clinicians were less likely to prescribe anti-TNF therapy to patients who were perceived to have the highest risk of infection. Also, other potential unmeasured confounders such as frailty are not captured in our baseline covariates.

The study population reflects a predominantly Caucasian cohort based in the UK, and so generalisations to cohorts of other ethnicities or in regions where the background rates of infection differ may not be appropriate. Although the nbDMARD-treated patients were recruited from a limited number of centres, both cohorts came from a mixture of teaching and non-teaching hospitals spread geographically across the UK. While the observational design has important limitations, the data presented here reflect real world experience of anti-TNF therapy in RA in the UK, without patients having to fit the strict selection criteria often applied by randomised controlled trials.

Conclusions

There is an increased incidence of skin infections (including shingles) in RA patients prescribed anti-TNF therapy. Some of this risk may be a reflection of other co-morbidities or the concurrent use of corticosteroids. A better understanding of the individual risk profiles of different biological agents will enable patients and clinicians to make more personalised decisions regarding care.

Acknowledgments

The authors acknowledge the enthusiastic collaboration of all consultant rheumatologists and their specialist nurses in the UK in providing the data. In addition, the authors acknowledge the support from Dr Ian Griffiths (past) and Professor David Isenberg (current), chairs of the BSRBR Management Committee; Professor Gabriel Panayi, Professor David GI Scott, Dr Andrew Bamji, Professor Deborah Bax and Professor DL Scott, presidents of the British Society for Rheumatology (BSR) during the period of data collection, for their active role in enabling the register to undertake its tasks, and to Samantha Peters (CEO of the BSR), Mervyn Hogg, Nia Taylor and members of the BSRBR Scientific Steering Committee. The authors also acknowledge the seminal role of the BSR Clinical Affairs Committee for establishing national biological guidelines and recommendations for such a register. Finally, the authors would like to acknowledge the substantial contribution of Andy Tracey, Katie McGrother and Dr Mark Lunt in database design and manipulation and Professor Alan Silman in his previous role as a principal investigator of the BSRBR and BSRBR Control Centre Consortium. The BSRBR Control Centre Consortium consists of the following institutions (all in the UK): Antrim Area Hospital, Antrim (Dr Nicola Maiden); Cannock Chase Hospital, Cannock Chase (Dr Tom Price); Christchurch Hospital, Christchurch (Dr Neil Hopkinson); Royal Derby Hospital, Derby (Dr Sheila O’Reilly); Dewsbury and District Hospital, Dewsbury (Dr Lesley Hordon); Freeman Hospital, Newcastle-upon-Tyne (Dr Ian Griffiths); Gartnavel General Hospital, Glasgow (Dr Duncan Porter); Glasgow Royal Infirmary, Glasgow (Professor Hilary Capell); Haywood Hospital, Stoke-on-Trent (Dr Andy Hassell); Hope Hospital, Salford (Dr Romela Benitha); King's College Hospital, London (Dr Ernest Choy); Kings Mill Centre, Sutton-In-Ashfield (Dr David Walsh); Leeds General Infirmary, Leeds (Professor Paul Emery); Macclesfield District General Hospital, Macclesfield (Dr Susan Knight); Manchester Royal Infirmary, Manchester (Professsor Ian Bruce); Musgrave Park Hospital, Belfast (Dr Allister Taggart); Norfolk and Norwich University Hospital, Norwich (Professor David Scott); Poole General Hospital, Poole (Dr Paul Thompson); Queen Alexandra Hospital, Portsmouth (Dr Fiona McCrae); Royal Glamorgan Hospital, Glamorgan (Dr Rhian Goodfellow); Russells Hall Hospital, Dudley (Professor George Kitas); Selly Oak Hospital, Selly Oak (Dr Ronald Jubb); St Helens Hospital, St Helens (Dr Rikki Abernethy); Weston General Hospital, Weston-super-Mare (Dr Shane Clarke/Dr Sandra Green); Withington Hospital, Manchester (Dr Paul Sanders); Withybush General Hospital, Haverfordwest (Dr Amanda Coulson); North Manchester General Hospital (Dr Bev Harrison); Royal Lancaster Infirmary (Dr Marwan Bukhari) and The Royal Oldham Hospital (Dr Peter Klimiuk).

References

Supplementary materials

Footnotes

  • Funding Funding for this project was provided by the British Society for Rheumatology.

  • Competing interests The BSR commissioned the Biologics Register (BSRBR) as a UK-wide national project to investigate the safety of biological agents in routine medical practice. DPMS and KLH are principal investigators on the BSRBR. BSR receives restricted income from UK pharmaceutical companies, presently Abbott Laboratories, Biovitrum, Merk Sharpe & Dohme, Pfizer and Roche. This income finances a wholly separate contract between the BSR and the University of Manchester. The principal investigators and their team have full academic freedom and are able to work independently of pharmaceutical industry influence. All decisions concerning analyses, interpretation and publication are made autonomously of any industrial contribution. Members of the Manchester team, BSR trustees, committee members and staff complete an annual declaration in relation to conflicts of interest. The authors declare no other competing interests.

  • Ethics approval Ethics approval for this study was obtained from the Multicentre Research Ethics Committee for the northwest of England.

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