Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Risk factors for cutaneous reactions to allopurinol in Kinh Vietnamese: results from a case-control study

Arthritis Research & Therapy volume 22, Article number: 182 (2020) Cite this article

Abstract

Objective

The aim of this study was to investigate risk factors for cutaneous adverse reactions (CARs) in Kinh Vietnamese.

Methods

All patients were prospectively recruited in Ho Chi Minh City. Presence of the HLA-B*58:01 allele was determined by real-time PCR-sequence-specific amplification by using the PG5801 Detection Kit (Pharmigene, Taipei). Patients with severe (SCARs) and mild (MCARs) CARs and controls were compared for differences in features prospectively collected, and odds ratios (ORs) with 95% confidence intervals (CIs) were estimated.

Results

On comparing 32 patients with SCARs and 395 tolerant controls, we identified eight strong risk factors: increased age (OR 15.1 [95% CI 5.8–40.1], P < 0.0001), female sex (OR 333 [40–43,453], P < 0.0001), allopurinol for asymptomatic hyperuricemia (OR 955 [120–125,847], P < 0.0001), allopurinol starting dose > 150 mg (OR 316 [101–122], P < 0.0001), diuretics intake (OR 304 [35–40,018], P < 0.0001), eGFR < 60 ml/min/1.73 m2 (OR 100 [32–353], P < 0.0001), history of allopurinol-induced skin reaction (OR 78 [6–10,808], P = 0.004), and HLA-B*58:01 carriage (OR 147 [45–746], P < 0.0001). HLA-B*58:01 allele frequency in controls was 7.3%. For MCARs (n = 74), risk factors were eGFR < 60 ml/min/1.73 m2 (OR 4.9 [1.61–14.6], P = 0.006), history of allopurinol-induced skin reaction (OR 27 [2–3777], P = 0.01), and asymptomatic hyperuricemia (OR 27 [2–3777], P = 0.01).

Conclusion

This study confirmed 8 risk factors, including HLA-B*58:01, for SCARs and identified 3 risk factors for MCARs in Kinh Vietnamese. HLA-B*58:01 genotyping could guide the indication for allopurinol in Kinh Vietnamese patients with gout.

Introduction

After more than 60 years of discovery and medical use, allopurinol remains the most commonly used urate-lowering agent worldwide, and the first choice of urate-lowering therapy in an overwhelming majority of gout patients [1, 2]. Allopurinol is well tolerated in most cases, but cutaneous adverse reactions (CARs) may occur a few months after allopurinol initiation [3,4,5,6]. Mild CARs (MCARs) are the most common. Their incidence has been estimated at 2 to 4% in Caucasian allopurinol initiators [7, 8] and more than 15% in Taiwanese [9]. MCARs feature sparse maco-papular rash and no other organ involvement and quick resolution after drug discontinuation. In contrast, severe CARs (SCARs), characterized by widespread skin lesions, fever, organ involvement, and a mortality rate of 10 to 32% [3, 10, 11], are rare, although their frequency seems increased in Asians and some other non-Caucasian ethnicities [3, 9, 12, 13]. Allopurinol-induced SCARs include Stevens-Johnson syndrome/toxic epidermal necrosis (SJS/TEN), drug reaction with eosinophilia and systemic symptoms (DRESS), and acute generalized exanthematous pustulosis. Allopurinol was identified as the leading cause of SCARs [5, 14] and the second cause of DRESS [15], whereas acute generalized exanthematous pustulosis seems extremely rare.

The risk of hospitalized allopurinol-induced SCARs has been recently estimated at 0.051% (95% CI 0.045–0.059) in US whites and Hispanics, 0.070% (0.039 to 0.115) in Asians, and 0.082% (0.066 to 0.10) in Blacks [13]. Previous studies found that the incidence of SCARs due to allopurinol was higher in Asian than Caucasian populations. The incidence was estimated at 0.16% per year in Thailand [16], and for patients with chronic renal insufficiency (CRI), 2% and 2.3% per year in Korea [17] and Taiwan [9], respectively. More recently, the incidence of SCARs between 2008 and 2015 in Taiwan has been estimated at 0.10% per year; older estimations based on 2001–2004 data were close to 0.30% per year [18]. With this heterogeneity between and within populations, a search for genetic and clinical risk factors is needed.

Several strong risk factors for allopurinol-induced SCARs include recent (< 3 months) allopurinol introduction, high allopurinol starting dose, impaired renal function, and HLA-B*58:01 positivity [19,20,21,22]. The link between SCARs and the HLA-B*58:01 allele was first identified in Han Chinese in Taiwan [23]. The strength of the association was reported to vary across ethnicities. Strong associations were found in Han Chinese, Korean, and Thai people [24,25,26], whereas studies of Japanese and Caucasian populations reported lower allele prevalence and weaker associations [27,28,29]. Variations in the frequency of the at-risk allele across ethnicities appeared to explain at least in part the racial variations in the incidence of SCARs [12, 13]. To reduce the incidence of SCARs, routine HLA-B*58:01 screening in candidates for allopurinol was found cost-effective and recommended for certain ethnicities with high allele prevalence and a strong association [1, 16, 18, 30, 31]. In contrast to SCARs, for allopurinol-induced MCARs, risk factors have been little investigated.

Regarding its high therapeutic effectiveness and affordable cost, allopurinol is commonly prescribed to treat gouty arthritis and hyperuricemia in both public healthcare centers and private clinics in Vietnam. Kinh Vietnamese, the major ethnic group in Vietnam, present genomic profile and HLA allele distribution similar to surrounding populations including Thai and Han Chinese [32, 33]. However, little is known about allopurinol-induced skin reactions in Kinh Vietnamese population, particularly the clinical risk factors and the association with HLA-B*58:01. The aim of this study was to investigate the risk factors, including HLA-B*58:01 carriage, for MCARs and SCARs in Kinh Vietnamese, who account for up to 90% of the nearly 100 million inhabitants in Vietnam.

Methods

Patient recruitment

The study was approved by the Ethics Committee of University of Medicine and Pharmacy at Ho Chi Minh City (HCMC) (HEC/IRB no. 198/ĐHYD-HĐ). Oral and written informed consent was obtained for all patients. All patients were self-identified as Kinh Vietnamese and were prospectively recruited in HCMC. Patients diagnosed with allopurinol-induced SCARs were prospectively recruited in hospital departments (mostly Tropical Diseases, Cho Ray Hospital), where they had been admitted because of the severity of their skin reaction. Patients with mild skin reactions were recruited at the Vien Gut Clinic and at the outpatient clinic of the same hospitals. Patients who had been taking allopurinol and did not show any skin reactions 3 months after the latest escalated dose were prospectively and consecutively recruited at the Vien Gut Clinic as allopurinol-tolerant controls.

Clinical diagnosis and data collection

Physicians in charge diagnosed the skin reactions to allopurinol. Allopurinol had to be used for < 90 consecutive days, and patients with concomitant use of another drug with high allergy potential were excluded. Photographs of skin manifestations were reviewed and confirmed by certified dermatologists. SCARs were diagnosed according to RegiSCAR criteria [34]. Patients with a combination of high fever, maculopapular rash, lymph node enlargement, blood count abnormalities, and internal organ involvement were classified as having DRESS syndrome. Liver and kidney involvements were defined as more than a two-fold increase of serum liver transaminase levels from the upper limit of the reference range and a 1.5-fold increase of serum creatinine level from the pre-eruption value, respectively. Unspecific skin eruptions not classified as any type of SCAR were considered MCARs.

Clinical history, laboratory test results, and treatment regimen were documented. Concomitant diuretics intake was defined as thiazide or loop diuretics taken continuously during the allopurinol use. A baseline (before skin reaction) estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73 m2 was considered CRI.

DNA extraction and HLA-B*58:01 genotyping

For each patient, 2 ml peripheral venous blood was collected in a tube coated with ethylenediaminetetraacetic acid, and genomic DNA of mononuclear blood cells was extracted by using the GeneJET Whole Blood Genomic DNA purification mini kit (Thermo Fisher Scientific, Waltham, MA, USA). The DNA concentration was measured by using the NanoDrop 2000 Spectrophotometer.

The presence of HLA-B*58:01 allele was first determined by real-time polymerase chain reaction (PCR) sequence-specific amplification with the PG5801 Detection Kit (Pharmigene, Taipei). All PCR reactions were performed with Mastercycler Realplex4 Epgradient S (Eppendorf, Hamburg, Germany). The zygosity of HLA-B*58:01-positive samples was identified by DNA sequencing. In brief, PCR for HLA-B exon 2, 3 and the surrounding region was performed with 5-UT-F and Bin3M13-R primers (Supplementary table 1). Amplification involved using Mastercycler@proS (Eppendorf, Hamburg, Germany) with the following conditions: one cycle at 98 °C for 3 min and 40 cycles of denaturation at 98 °C for 20 s, annealing at 62 °C for 20 s, extension at 72 °C for 1 min, and a final extension at 72 °C for 2 min. The length of amplified products was 1073 bp and was confirmed by electrophoresis on 1% agarose gels with Diamond Nucleic Acid Dye (Promega, Madison, WI, USA). PCR products were cleaned with Exosap-IT glycerol solution (Thermo Fisher Scientific) and sequenced by using the BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA, USA) with Seq-BIn2-R and Bin3M13-R primers (Supplementary Table 1). Sequencing reactions were analyzed with the use of ABI 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Results were compared with the reference sequence of HLA-B*58:01 from the IMGT database [35]. Only samples that 100% matched the reference sequence and with unambiguity were considered homozygous. In the case of ambiguity, the results were identified as heterozygous.

All techniques were performed in the Center for Molecular Biomedicine of the University of Medicine and Pharmacy at Ho Chi Minh City according to the manufacturer’s protocol.

Statistical analysis

Descriptive data are presented as mean (standard deviation, SD) for quantitative variables or number (percentage, %) for qualitative variables. Differences were first tested between the 3 groups, then between controls and SCAR and MACR groups for items with a global significant difference. Global differences were tested by Kruskal-Wallis rank-sum test and Fisher exact test for quantitative and categorical variables, respectively. Two-by-two group comparisons were performed with the Wilcoxon rank-sum and Fisher exact tests for quantitative and categorical variables, respectively. The strength of the association was estimated by odds ratios (ORs) with 95% confidence intervals (95% CIs). Considering that some contingency tables include null counts, Firth’s bias-reduced penalized-likelihood logistic regression was systematically applied to obtain valid OR estimates [36, 37]. Due to the small number of cases, we cannot adjust on all potential confounders without being at risk of model over-parametrization. Thus, we proposed to stratify the impact of HLA-B*5801 on eGR, diuretic intake, and age above 65 years. Sensitivity and specificity of HLA-B*58:01 allele presence and positive and negative likelihood ratios with 95% CIs were estimated. Predictive and negative predictive values (PPV and NPV) depend on the incidence of SCARs in the studied population [38] and were estimated from previously estimated incidences of SCARs in Asian populations. All tests were two-sided, with P < 0.05 considered statistically significant. All statistical analyses were performed with R v3.6.1 [39].

Results

In total, 105 patients with allopurinol-induced skin reactions were recruited from October 2017 to May 31, 2018; 74 had MCARs, which developed after a mean allopurinol exposure of 18.5 days (range 1 to 89 days); one additional patient was excluded because he took celecoxib and omeprazole in addition to allopurinol when he developed his skin reaction. Thirty-one patients were hospitalized for SCARs (29 SJS/TEN and 2 DRESS), which developed after a mean exposure of 15.4 days (range 1 to 40 days). None of the patients died. We included 395 tolerant controls from October 2017 to May 31, 2018.

Study population characteristics

The mean (SD) age of the SCAR group was 60 (15.9) years, significantly higher than both the MCAR and control groups (45.7 [11.7] and 45.4 [10.3] years, respectively). All MCAR and control patients were male, but 9 of 31 SCAR patients were female. The control group had the highest mean allopurinol dose at the time of DNA sampling. Other factors, including history of allopurinol-induced skin reaction, asymptomatic hyperuricemia, and concomitant diuretic intake, were more commonly found in the SCAR than MCAR group and were totally absent in controls (Table 1). The proportion of patients with a starting dose of allopurinol > 150 mg/day was higher in the SCAR than MCAR and control groups (87.1% vs 4.1% and 1.8%). Baseline CRI (eGFR < 60 ml/min/1.73 m2) was noted in 66.7%, 8.1%, and 1.8% of SCAR, MCAR, and control groups, respectively. Of note, serum creatinine had not been measured before the toxic reaction in 13 SCAR patients. Overall, 29 (93.6%) of all SCAR patients, including 27 of the 29 SJS/TEN patients and the 2 DRESS patients, were HLA-B*58:01 positive, whereas only a few MCAR patients (6/74, 8.1%) and controls (29/395, 7.3%) carried the allele (Table 1). The proportion of HLA-B*58:01-positive patients who were homozygote for the allele was similar in the SCAR, MCAR, and control groups (25.9%, 16.4%, and 24.1%, respectively).

Table 1 Characteristics of Vietnamese Kinh people with cutaneous adverse reactions (MCARs), rare severe CARs (SCARs), and controls (no reaction to allopurinol)
Full size table

Identified risk factors for allopurinol-induced MCARs

For allopurinol-induced MCARs, we found only three risk factors: CRI, history of allopurinol-induced skin reaction, and asymptomatic hyperuricemia (Table 2). HLA-B*58:01 positivity was not associated with MCARs.

Table 2 Identified risk factors for allopurinol-induced MCARs
Full size table

Identified risk factors for allopurinol-induced SCARs

We identified seven risk factors for allopurinol-induced SCARs, with significant difference and high ORs between SCAR and control groups: female sex, age > 65 years, history of allopurinol-induced skin reaction, high starting dose of allopurinol, concomitant diuretics intake, eGFR < 60 ml/min/1.73 m2, and HLA-B*58:01 positivity (Table 3). More than half of the SCAR patients received allopurinol to treat asymptomatic hyperuricemia in contrast to our controls who all had gout, as the Vien Gut Medical Centre, where controls were recruited, follows the international guideline not to treat asymptomatic hyperuricemia. SCAR patients also more frequently had hypertension (P = 0.03) and dyslipidemia (P = 0.02) than controls. Mean allopurinol starting dose and ratio of starting dose to eGFR were higher with SCARs. An allopurinol starting dose > 150 mg/day was strongly associated with SCARs. SCARs developed only within patients with a ratio of starting dose to eGFR > 2.74. A similar strong association was found using a threshold ratio of 4 (corresponding to the first quartile of the starting dose/eGFR ratio in SCAR patients) with the advantage of greater specificity. The probability of HLA-B*58:01 positivity with the occurrence of allopurinol-induced SCARs was high (OR 147 [95% CI 38–564]) and was stronger with the concomitant risk factors age > 65 years (3372 [136–83,458]), eGFR < 60 ml/min/1.73 m2 (524 [91–3002]), concomitant diuretics intake (2199 [81–59,866]), and starting dose > 150 mg/day (7770 [346–174,435]) as compared with HLA-B*58:01 negativity (Fig. 1). Of note, all patients with SCARs and with renal insufficiency or > 65 years old were HLA-B*58:01 positive.

Table 3 Identified risk factors for allopurinol-induced SCARs
Full size table
Fig. 1
figure 1

Probability of severe cutaneous adverse reactions (SCARs) in the total population of Vietnamese Kinh people and in those with starting allopurinol dose > 150 mg/day, renal failure, concomitant diuretic intake, and age ≥ 65 years. Data are odds ratios (95% confidence intervals)

Full size image

The interest of HLA-B*58:01 testing to predict SCARs in Kinh Vietnamese was determined by several approaches. Sensitivity and specificity were estimated at 93.5% (95% CI 78.6–99.2) and 92.5% (89.8–94.7), respectively. Positive and negative likelihood ratios were estimated at 12.5 (9.0–17.5) and 0.07 (0.02–0.27). Therefore, both tests strongly supported the value of HLA-B*58:01 genotyping, because likelihood ratios > 10 and < 0.1 are considered strong evidence to rule in or rule out diagnoses, respectively [40]. Considering incidence figures observed in South Asia populations (see the “Introduction” section), a SCAR incidence < 0.30% per year led to estimate a PPV < 3.6% and NPV > 99.9%, respectively. With retained prevalence < 1/1000, the estimated PPV decreased to 1.2% and the NPV remained > 99.9%. In our patients with age > 65 years or eGFR < 60 ml/min/1.73 m2 in whom the SCAR prevalence could increase to 2/100, estimated PPVs reached 22.2% and 12.2%, respectively, and NPVs were still > 99.9%.

Discussion

The most common adverse reactions to allopurinol are MCARs [11, 41], which, in our study, associated with CRI, previous skin reaction to allopurinol, and management of asymptomatic hyperuricemia. Only 8.1% of MCAR patients were positive for HLA-B*58:01, for an insignificant strength of association. This finding agrees with reports from Portugal, Australia, and Korea [17, 29, 41]. In contrast, studies of Han Chinese and Thai people with MCARs revealed high HLA-B*58:01 carriage (85–100%) and a strong association similar to that of SCARs [24, 42].

SCARs are extremely severe and life-threatening and are the most important hazard for patients with allopurinol treatment. Many studies have focused on the identification of risk factors to minimize the risk of SCARs [3, 21, 23, 25, 29, 43], but these have been little investigated in Vietnamese people. Our study found several non-genetic and genetic risk factors for allopurinol-induced SCARs in the Kinh Vietnamese population, which confirmed the findings for other ethnicities [3, 10, 22].

We found that starting allopurinol at > 150 mg/day was highly associated with SCARs (OR 316 [95% CI 101–1224]), which confirmed previous results about the importance of starting dose [4, 13, 21], even though the starting dose in our tolerant controls was higher than the dose of ≤ 100 mg/day usually recommended [1, 44]. Our tolerant controls were recruited at the Vien Gut Medical Centre, where the policy is to restrict allopurinol prescription to patients with eGFR > 60 ml/min/1.73 m2 and to start the drug at 150 mg/day (owing to the sole availability of 300 mg tablets in this center), with subsequent escalation of the dose by 150 mg/day every month until the desirable uricemia target is reached. Because no SCAR was observed in more than 10,000 gouty patients prescribed allopurinol according to this protocol at the Vien Gut Medical Centre over the last 4 years, 150 mg/day may be an acceptable starting dose in Vietnamese patients without renal insufficiency. This approach is supported by the fact that all the SCAR patients in the present study had a starting dose of at least 300 mg/day or eGFR < 60 ml/min/1.73 m2. Allopurinol was finally given at a significantly higher dose in our tolerant controls than in those with skin reactions, which confirms the involvement of starting dose and not maintenance dose in provoking skin reactions. Strikingly, more than 50% of patients with SCARs had received allopurinol for asymptomatic hyperuricemia. Other investigators also found that allopurinol prescription for asymptomatic hyperuricemia notably increased the risk of SCARs [4]. Thus, great caution is warranted in treating asymptomatic hyperuricemia with allopurinol.

Concomitant use of diuretics was also a risk factor for allopurinol-induced SCARs (OR 304 [95% CI 35–40,018) similar to other reports [13, 19, 20, 22]. This finding could be related to the decreased renal clearance of oxypurinol, which has been documented in patients taking thiazides [45] and furosemide [46]. However, the exact role of diuretics in allopurinol-induced skin reactions remains controversial, because studies of both Asians [4, 17, 23] and Caucasians [5] showed that taking diuretics with allopurinol, after adjustment for other factors, in particular comorbidities, did not significantly increase the risk of hypersensitivity. We did not use multiple adjustments in our study because of the limited number of included SCAR patients. Of note, our 7 SCAR patients who took diuretics took them for hypertension.

CRI has been highlighted as an independent risk factor for allopurinol-induced SCARs in many studies [13, 20, 22, 24, 43, 47], leading to recommend limiting the allopurinol dose based on baseline eGFR to minimize the risk [44]. Most regulatory agencies declared this recommendation, but it has recently been debated because it appeared to compromise the urate-lowering effect of allopurinol in renal patients [48]. Also, in some studies, SCARs seemed to occur regardless of adjustment of allopurinol dose to renal function [3, 5, 23, 48], whereas in others, plasma oxypurinol level was correlated with prognosis rather than SCARs [10, 20, 48]. In our study, CRI was associated with both SCARs (OR 100 [95% CI 32–353]) and MCARs (OR 4.9 [1.6–14.6]). However, this result must be taken with caution because 13 of our 31 SCAR patients had no creatinine measurement before their eruption, and patients with renal failure were not included in our tolerant control group because of the mode of recruitment.

Several in vitro studies shed some light on the idiosyncratic characteristics of allopurinol-induced SCARs. Delayed hypersensitivity via CD8+ T cells seems to play a central role in the pathogenesis of SCARs, involving the antigen presentation process. The response requires a sufficient amount of drug molecules binding to peptide–major histocompatibility complex and also a certain specificity between antigen-presenting cells and CD8+ T cells [49]. This finding led to two crucial categories of risk factors: the first includes drug dose, reduced clearance due to impaired renal function or diuretics, and increased oxypurinol plasma concentration, and the second, regulation of oxypurinol presentation to immune cells by immune receptors such as HLA-B*58:01 [49].

HLA-B*58:01 has emerged as a promising genetic marker to predict allopurinol-induced SCARs. However, substantial difference by ethnicity was an important determinant of clinical cost-effective analyses: HLA-B*58:01 screening was beneficial in certain countries with high-risk ethnic groups but was ineffective in others [50]. Our Ho Chi Minh study found a strong association between HLA-B*58:01 and allopurinol-induced SCARs (OR 147 [95% CI 45–745]) in accordance with a previous study in Hanoi [51]. Also, the frequency of the HLA-B*58:01 genotype was 7.3% in our tolerant gouty patients, similar to the results for the general Kinh population of both South and North Vietnam [33, 52]. Hence, our results suggest that Kinh Vietnamese are another high-risk group in whom HLA-B*58:01 would be worth genotyping before prescribing allopurinol.

Within the HLA-B locus, the number of 58:01 alleles can be one (heterozygous) or two (homozygous) and can result in a gene-dosage effect. Ng et al. reported a higher incidence of SCARs in homozygotes than heterozygotes (OR 72.45 vs 15.25) [47]. Shim et al. showed similar results, but the association was not statistically significant [53]. In our study, risk of SCARs was not higher in patients homozygous than heterozygous for HLA-B*58:01.

The PPVs and NPVs estimated in our work agree with previous reports finding a PPV for carriage of < 10% [47, 54], whereas with a non-carriage state, the probability of allopurinol-induced SCAR was extremely unlikely: the NPV was almost 100% [55]. Our results for patients with renal failure also agree with those of Jung et al., who found a PPV and NPV of 18% and 100%, respectively, in a Korean population with CRI [17].

Given both the cost-effectiveness of allopurinol therapy and the rare occurrence of SCARs, a single risk factor could be insufficient to justify absolute restraint from this drug. To date, the most widely recommended approach to risk evaluation of allopurinol-induced SCAR has been to combine multiple well-studied factors. Lee et al. proposed screening for HLA-B*58:01 in patients with one or more characteristics including baseline eGFR < 60 mL/min/1.73 m2, co-administration with multiple drugs especially diuretics, East Asian or Southeast Asian background, and history of allopurinol intolerance [41]. Jung et al. reported a 6-fold rate of SCARs in patients with HLA-B*58:01 carriage and CRI versus those with normal renal function (18% vs 2.7%) [17], which resulted in necessary avoidance of allopurinol in CRI patients carrying the at-risk allele. Ng et al. reported that a combination of HLA-B*58:01 and poor renal function significantly improved the ORs, PPV, and area under the receiver operating characteristic curve as compared with a single risk factor [47].

In our study of Kinh Vietnamese, with several risk factors (8 for SCARs and 3 for MCARs), a number should be avoided such as starting allopurinol at a high dose and treating asymptomatic hyperuricemia or a history of skin intolerance. Our results support that patients with no clinical risk factor could be started on allopurinol, with a starting dose not exceeding 150 mg/day. HLA typing could be advised in patients taking diuretics, older people, and those with CRI, to guide allopurinol prescription, which could be allowed in HLAB*58:01-negative patients, at a low starting dose (Fig. 2), as supported by the high ORs for SCARs associated with the at-risk allele in these patients (Fig. 1).

Fig. 2
figure 2

A proposed decision tree for allopurinol prescription in Kinh Vietnamese

Full size image

This study has limitations. Despite prospective enrolment, a number of data were missing, especially in the SCAR patient group as these patients went to the hospital on emergency and had not been followed-up as out-patients. The mode of recruitment of tolerant controls at the Vien Gut Medical Centre excluded asymptomatic hyperuricemia patients as this center follows the policy not to treat asymptomatic hyperuricemia. This made a valid comparison with cutaneous intolerance patients impossible. The case-control study design we used exposes to a higher risk of bias than a cohort study. The main risk of bias lies in the difficulty of collecting the exposure factor (recall bias), but for a genetic study, this risk is very low as exposure is assessed by PCR. Moreover, in the presence of rare disease or event, it is the only way to achieve results within a reasonable time and at a reasonable cost. Nevertheless, we cannot rule out that confounders such as advanced age or comorbidities could impact the results of our study. Due to the small number of cases, we could not adjust on all potential confounders without being at risk of model over-parametrization. Thus, we proposed to stratify the impact of HLA-B*5801 on eGFR, diuretic intake, and age above 65 years (Fig. 1). Results still showed an impact of the HLA-B*5801 and seem robust when age or renal insufficiency are taking into account. In logistic regression analysis of small or sparse data sets, results obtained by classical maximum likelihood methods cannot be generally trusted. This situation has been termed “separation,” and it typically occurs whenever no or few events are observed in one of the two groups defined by a dichotomous covariate and in the presence of highly predictive risk factors [37, 56]. In such situation, Heinze and Schemper have demonstrated that Firth’s modification of the logistic model unbiased estimates and valid confidence intervals. However, the width of confidence intervals should lead to relativize the value of some point estimates.

Conclusion

The multiple risk factors for allopurinol-induced SCARs in Kinh Vietnamese identified in this study should further raise awareness about the risk of hypersensitivity and the need for careful evaluation before allopurinol prescription. Considering that the cost of HLA*58:01 testing is gradually decreasing in Vietnam, identifying genetic risk factors could be extremely valuable to rule out candidates with inevitable risk factors for allopurinol treatment. As a result, a urate-lowering strategy could be wisely tailored to achieve both therapeutic safety and effectiveness.

Availability of data and materials

The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CAR:

Cutaneous adverse reaction

CI:

Confidence interval

CRI:

Chronic renal insufficiency

DRESS:

Drug reaction with eosinophilia and systemic symptoms

eGFR:

Estimated glomerular filtration rate

HCMC:

Ho Chi Minh City

HLA:

Human leucocyte antigen

MCAR:

Mild cutaneous adverse reaction

NA:

Not applicable

NPV:

Negative predictive value

OR:

Odds ratio

PCR:

Polymerase chain reaction

PPV:

Positive predictive value

SCAR:

Severe cutaneous adverse reaction

SD:

Standard deviation

SJS/TEN:

Stevens-Johnson syndrome/toxic epidermal necrosis

References

  1. Khanna D, FitzGerald JD, Khanna PP, Bae S, Singh M, Neogi T, et al. 2012 American College of Rheumatology guidelines for management of gout. Part I: systematic non-pharmacologic and pharmacologic therapeutic approaches to hyperuricemia. Arthritis Care Res. 2012;64(10):1431–46.

    CAS  Google Scholar 

  2. Rai SK, Aviña-Zubieta JA, McCormick N, De Vera MA, Shojania K, Sayre EC, et al. The rising prevalence and incidence of gout in British Columbia, Canada: population-based trends from 2000 to 2012. Semin Arthritis Rheum. 2017;46(4):451–6.

    PubMed  Google Scholar 

  3. Ramasamy SN, Korb-Wells CS, Kannangara DRW, Smith MWH, Wang N, Roberts DM, et al. Allopurinol hypersensitivity: a systematic review of all published cases, 1950-2012. Drug Saf. 2013;36(10):953–80.

    CAS  PubMed  Google Scholar 

  4. Yang C-Y, Chen C-H, Deng S-T, Huang C-S, Lin Y-J, Chen Y-J, et al. Allopurinol use and risk of fatal hypersensitivity reactions: a nationwide population-based study in Taiwan. JAMA Intern Med. 2015;175(9):1550–7.

    PubMed  Google Scholar 

  5. Halevy S, Ghislain P-D, Mockenhaupt M, Fagot J-P, Bouwes Bavinck JN, Sidoroff A, et al. Allopurinol is the most common cause of Stevens-Johnson syndrome and toxic epidermal necrolysis in Europe and Israel. J Am Acad Dermatol. 2008;58(1):25–32.

    PubMed  Google Scholar 

  6. Kim SC, Newcomb C, Margolis D, Roy J, Hennessy S. Severe cutaneous reactions requiring hospitalization in allopurinol initiators: a population-based cohort study. Arthritis Care Res. 2013;65(4):578–84.

    CAS  Google Scholar 

  7. McInnes GT, Lawson DH, Jick H. Acute adverse reactions attributed to allopurinol in hospitalised patients. Ann Rheum Dis. 1981;40(3):245–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Hoigné R, Sonntag MR, Zoppi M, Hess T, Maibach R, Fritschy D. Occurrence of exanthems in relation to aminopenicillin preparations and allopurinol. N Engl J Med. 1987;316(19):1217.

    PubMed  Google Scholar 

  9. Lin C-W, Huang W-I, Chao P-H, Chen W-W, Hsiao F-Y. Risk of cutaneous adverse reactions associated with allopurinol or febuxostat in real-world patients: a nationwide study. Int J Clin Pract. 2019;73(5):e13316.

    PubMed  Google Scholar 

  10. Stamp LK, Day RO, Yun J. Allopurinol hypersensitivity: investigating the cause and minimizing the risk. Nat Rev Rheumatol. 2016;12(4):235–42.

    CAS  PubMed  Google Scholar 

  11. Mockenhaupt M. Epidemiology of cutaneous adverse drug reactions. Chem Immunol Allergy. 2012;97:1–17.

    PubMed  Google Scholar 

  12. Lu N, Rai SK, Terkeltaub R, Kim SC, Menendez ME, Choi HK. Racial disparities in the risk of Stevens-Johnson Syndrome and toxic epidermal necrolysis as urate-lowering drug adverse events in the United States. Semin Arthritis Rheum. 2016;46(2):253–8.

    PubMed  PubMed Central  Google Scholar 

  13. Keller SF, Lu N, Blumenthal KG, Rai SK, Yokose C, Choi JWJ, et al. Racial/ethnic variation and risk factors for allopurinol-associated severe cutaneous adverse reactions: a cohort study. Ann Rheum Dis. 2018;77(8):1187–93.

    PubMed  PubMed Central  Google Scholar 

  14. Atzori L, Pinna AL, Mantovani L, Ferreli C, Pau M, Mulargia M, et al. Cutaneous adverse drug reactions to allopurinol: 10 year observational survey of the dermatology department--Cagliari University (Italy). J Eur Acad Dermatol Venereol. 2012;26(11):1424–30.

    CAS  PubMed  Google Scholar 

  15. Kardaun SH, Sekula P, Valeyrie-Allanore L, Liss Y, Chu CY, Creamer D, et al. Drug reaction with eosinophilia and systemic symptoms (DRESS): an original multisystem adverse drug reaction. Results from the prospective RegiSCAR study. Br J Dermatol. 2013;169(5):1071–80.

    CAS  PubMed  Google Scholar 

  16. Saokaew S, Tassaneeyakul W, Maenthaisong R, Chaiyakunapruk N. Cost-effectiveness analysis of HLA-B*5801 testing in preventing allopurinol-induced SJS/TEN in Thai population. PLoS One. 2014;9(4):e94294.

    PubMed  PubMed Central  Google Scholar 

  17. Jung J-W, Song W-J, Kim Y-S, Joo KW, Lee KW, Kim S-H, et al. HLA-B58 can help the clinical decision on starting allopurinol in patients with chronic renal insufficiency. Nephrol Dial Transplant. 2011;26(11):3567–72.

    CAS  PubMed  Google Scholar 

  18. Ko T-M, Tsai C-Y, Chen S-Y, Chen K-S, Yu K-H, Chu C-S, et al. Use of HLA-B*58:01 genotyping to prevent allopurinol induced severe cutaneous adverse reactions in Taiwan: national prospective cohort study. BMJ. 2015;351:h4848.

    PubMed  PubMed Central  Google Scholar 

  19. Day RO, Kannangara DRW, Stocker SL, Carland JE, Williams KM, Graham GG. Allopurinol: insights from studies of dose-response relationships. Expert Opin Drug Metab Toxicol. 2017;13(4):449–62.

    CAS  PubMed  Google Scholar 

  20. Chung W-H, Chang W-C, Stocker SL, Juo C-G, Graham GG, Lee M-HH, et al. Insights into the poor prognosis of allopurinol-induced severe cutaneous adverse reactions: the impact of renal insufficiency, high plasma levels of oxypurinol and granulysin. Ann Rheum Dis. 2015;74(12):2157–64.

    CAS  PubMed  Google Scholar 

  21. Stamp LK, Taylor WJ, Jones PB, Dockerty JL, Drake J, Frampton C, et al. Starting dose is a risk factor for allopurinol hypersensitivity syndrome: a proposed safe starting dose of allopurinol. Arthritis Rheum. 2012;64(8):2529–36.

    CAS  PubMed  Google Scholar 

  22. Saksit N, Tassaneeyakul W, Nakkam N, Konyoung P, Khunarkornsiri U, Chumworathayi P, et al. Risk factors of allopurinol-induced severe cutaneous adverse reactions in a Thai population. Pharmacogenet Genomics. 2017;27(7):255–63.

    CAS  PubMed  Google Scholar 

  23. Hung S-I, Chung W-H, Liou L-B, Chu C-C, Lin M, Huang H-P, et al. HLA-B*5801 allele as a genetic marker for severe cutaneous adverse reactions caused by allopurinol. Proc Natl Acad Sci U S A. 2005;102(11):4134–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Cao Z, Wei Z, Zhu Q, Zhang J, Yang L, Qin S, et al. HLA-B*58:01 allele is associated with augmented risk for both mild and severe cutaneous adverse reactions induced by allopurinol in Han Chinese. Pharmacogenomics. 2012;13(10):1193–201.

    CAS  PubMed  Google Scholar 

  25. Tassaneeyakul W, Jantararoungtong T, Chen P, Lin P-Y, Tiamkao S, Khunarkornsiri U, et al. Strong association between HLA-B*5801 and allopurinol-induced Stevens-Johnson syndrome and toxic epidermal necrolysis in a Thai population. Pharmacogenet Genomics. 2009;19(9):704–9.

    CAS  PubMed  Google Scholar 

  26. Kang H-R, Jee YK, Kim Y-S, Lee CH, Jung J-W, Kim SH, et al. Positive and negative associations of HLA class I alleles with allopurinol-induced SCARs in Koreans. Pharmacogenet Genomics. 2011;21(5):303–7.

    CAS  PubMed  Google Scholar 

  27. Tohkin M, Kaniwa N, Saito Y, Sugiyama E, Kurose K, Nishikawa J, et al. A whole-genome association study of major determinants for allopurinol-related Stevens-Johnson syndrome and toxic epidermal necrolysis in Japanese patients. Pharmacogenomics J. 2013;13(1):60–9.

    CAS  PubMed  Google Scholar 

  28. Lonjou C, Borot N, Sekula P, Ledger N, Thomas L, Halevy S, et al. A European study of HLA-B in Stevens-Johnson syndrome and toxic epidermal necrolysis related to five high-risk drugs. Pharmacogenet Genomics. 2008;18(2):99–107.

    CAS  PubMed  Google Scholar 

  29. Gonçalo M, Coutinho I, Teixeira V, Gameiro AR, Brites MM, Nunes R, et al. HLA-B*58:01 is a risk factor for allopurinol-induced DRESS and Stevens-Johnson syndrome/toxic epidermal necrolysis in a Portuguese population. Br J Dermatol. 2013;169(3):660–5.

    PubMed  Google Scholar 

  30. Park D-J, Kang J-H, Lee J-W, Lee K-E, Wen L, Kim T-J, et al. Cost-effectiveness analysis of HLA-B5801 genotyping in the treatment of gout patients with chronic renal insufficiency in Korea. Arthritis Care Res. 2015;67(2):280–7.

    Google Scholar 

  31. Jutkowitz E, Dubreuil M, Lu N, Kuntz KM, Choi HK. The cost-effectiveness of HLA-B*5801 screening to guide initial urate-lowering therapy for gout in the United States. Semin Arthritis Rheum. 2017;46(5):594–600.

    PubMed  Google Scholar 

  32. Le VS, Tran KT, Bui HTP, Le HTT, Nguyen CD, Do DH, et al. A Vietnamese human genetic variation database. Hum Mutat. 2019;40(10):1664–75.

    CAS  PubMed  Google Scholar 

  33. Do MD, Le LGH, Nguyen VT, Dang TN, Nguyen NH, Vu HA, et al. High-resolution HLA typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by using next-generation sequencing. Front Genet. 2020;11 Available from: https://www.frontiersin.org/articles/10.3389/fgene.2020.00383/full?&utm_source=Email_to_authors_&utm_medium=Email&utm_content=T1_11.5e1_author&utm_campaign=Email_publication&field=&journalName=Frontiers_in_Genetics&id=510718. [cited 2020 May 4].

  34. Roujeau JC. The spectrum of Stevens-Johnson syndrome and toxic epidermal necrolysis: a clinical classification. J Invest Dermatol. 1994;102(6):28S–30S.

    CAS  PubMed  Google Scholar 

  35. Robinson J, Barker DJ, Georgiou X, Cooper MA, Flicek P, Marsh SGE. IPD-IMGT/HLA database. Nucleic Acids Res. 2020;48(D1):D948–55.

    PubMed  Google Scholar 

  36. Firth D. Bias reduction of maximum likelihood estimates. Biometrika. 1993;80(1):27–38.

    Google Scholar 

  37. Heinze G, Schemper M. A solution to the problem of separation in logistic regression. Stat Med. 2002;21(16):2409–19.

    PubMed  Google Scholar 

  38. Altman DG, Bland JM. Diagnostic tests. 1: Sensitivity and specificity. BMJ. 1994;308(6943):1552.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation; 2019. Available from: https://www.R-project.org/.

    Google Scholar 

  40. Deeks JJ, Altman DG. Diagnostic tests 4: likelihood ratios. BMJ. 2004;329(7458):168–9.

    PubMed  PubMed Central  Google Scholar 

  41. Lee MH, Stocker SL, Anderson J, Phillips EJ, Nolan D, Williams KM, et al. Initiating allopurinol therapy: do we need to know the patient’s human leucocyte antigen status? Intern Med J. 2012;42(4):411–6.

    CAS  PubMed  Google Scholar 

  42. Sukasem C, Jantararoungtong T, Kuntawong P, Puangpetch A, Koomdee N, Satapornpong P, et al. HLA-B (*) 58:01 for allopurinol-induced cutaneous adverse drug reactions: implication for clinical interpretation in Thailand. Front Pharmacol. 2016;7:186.

    PubMed  PubMed Central  Google Scholar 

  43. Park HJ, Yun J, Kang DY, Park JW, Koh YI, Kim S, Kim SH, Nam YH, Jeong YY, Kim CW, Park HK, Kim SH, Kang HR, Jung JW. Unique Clinical Characteristics and Prognosis of Allopurinol-Induced Severe Cutaneous Adverse Reactions. Drug Allergy Work Group of KAAACI. J Allergy Clin Immunol Pract. 2019;7(8):2739-49.e3. httsp://doi.org/10.1016/j.jaip.2019.05.047. Epub 2019 Jun 13.

  44. Richette P, Doherty M, Pascual E, Barskova V, Becce F, Castañeda-Sanabria J, et al. 2016 updated EULAR evidence-based recommendations for the management of gout. Ann Rheum Dis. 2017;76(1):29–42.

    CAS  PubMed  Google Scholar 

  45. de Vries JX, Voss A, Ittensohn A, Walter-Sack I, Löffler W, Landthaler R, et al. Interaction of allopurinol and hydrochlorothiazide during prolonged oral administration of both drugs in normal subjects. II. Kinetics of allopurinol, oxipurinol, and hydrochlorothiazide. Clin Investig. 1994;72(12):1076–81.

    PubMed  Google Scholar 

  46. Stamp LK, Barclay ML, O’Donnell JL, Zhang M, Drake J, Frampton C, et al. Furosemide increases plasma oxypurinol without lowering serum urate--a complex drug interaction: implications for clinical practice. Rheumatol Oxf Engl. 2012;51(9):1670–6.

    CAS  Google Scholar 

  47. Ng CY, Yeh Y-T, Wang C-W, Hung S-I, Yang C-H, Chang Y-C, et al. Impact of the HLA-B(*)58:01 allele and renal impairment on allopurinol-induced cutaneous adverse reactions. J Invest Dermatol. 2016;136(7):1373–81.

    CAS  PubMed  Google Scholar 

  48. Dalbeth N, Kumar S, Stamp L, Gow P. Dose adjustment of allopurinol according to creatinine clearance does not provide adequate control of hyperuricemia in patients with gout. J Rheumatol. 2006;33(8):1646–50.

    CAS  PubMed  Google Scholar 

  49. Yun J, Marcaida MJ, Eriksson KK, Jamin H, Fontana S, Pichler WJ, et al. Oxypurinol directly and immediately activates the drug-specific T cells via the preferential use of HLA-B*58:01. J Immunol Baltim Md 1950. 2014;192(7):2984–93.

    CAS  Google Scholar 

  50. Saito Y, Stamp LK, Caudle KE, Hershfield MS, McDonagh EM, Callaghan JT, et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines for human leukocyte antigen B (HLA-B) genotype and allopurinol dosing: 2015 update. Clin Pharmacol Ther. 2016;99(1):36–7.

    CAS  PubMed  Google Scholar 

  51. Do T, Tran T, Vu T, Tran Q, Nguyen T, Chu C, et al. Association between HLA*B5801 and allopurinol-induced severe cutaneous adverse drug reaction in Bachmai hospital. Tap Chi Hoc Phong. 2015;25:396–400.

    Google Scholar 

  52. Hoa BK, Hang NTL, Kashiwase K, Ohashi J, Lien LT, Horie T, et al. HLA-A, -B, -C, -DRB1 and -DQB1 alleles and haplotypes in the Kinh population in Vietnam. Tissue Antigens. 2008;71(2):127–34.

    CAS  PubMed  Google Scholar 

  53. Shim J-S, Yun J, Kim M-Y, Chung SJ, Oh JH, Kang D-Y, et al. The presence of HLA-B75, DR13 homozygosity, or DR14 additionally increases the risk of allopurinol-induced severe cutaneous adverse reactions in HLA-B*58:01 carriers. J Allergy Clin Immunol Pract. 2019;7(4):1261–70.

    PubMed  Google Scholar 

  54. Somkrua R, Eickman EE, Saokaew S, Lohitnavy M, Chaiyakunapruk N. Association of HLA-B*5801 allele and allopurinol-induced Stevens Johnson syndrome and toxic epidermal necrolysis: a systematic review and meta-analysis. BMC Med Genet. 2011;12:118.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Yu K-H, Yu C-Y, Fang Y-F. Diagnostic utility of HLA-B*5801 screening in severe allopurinol hypersensitivity syndrome: an updated systematic review and meta-analysis. Int J Rheum Dis. 2017;20(9):1057–71.

    CAS  PubMed  Google Scholar 

  56. Heinze G. A comparative investigation of methods for logistic regression with separated or nearly separated data. Stat Med. 2006;25(24):4216–26.

    PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

ART-Viggo and Société Française de Rhumatologie

Author information

Author notes

  1. Minh Duc Do, Thao Phuong Mai and Anh Duy Do contributed equally to this work.

Authors and Affiliations

  1. Center for Molecular Biomedicine, University of Medicine and Pharmacy, at Ho Chi Minh City, Vietnam

    Minh Duc Do, Linh Gia Hoang Le & Vu Anh Hoang

  2. Department of Physiology, Pathophysiology and Immunology, Faculty of Medicine, University of Medicine and Pharmacy, Ho Chi Minh City, Vietnam

    Thao Phuong Mai

  3. Department of Physiology, Pathophysiology and Immunology, Pham Ngoc Thach University of Medicine, Ho Chi Minh City, Vietnam

    Anh Duy Do

  4. French Vietnamese Research Center on Gout and Chronic Diseases, Vien Gut Medical Centre, Ho Chi Minh City, Vietnam

    Quang Dinh Nguyen, Nghia Hieu Le & Thomas Bardin

  5. Department of Scientific Research, Cho Ray Hospital, Ho Chi Minh City, Vietnam

    Anh Ngoc Le

  6. Department of Tropical Disease, Cho Ray Hospital, Ho Chi Minh City, Vietnam

    Hung Quoc Le

  7. Université de Paris, U1132, INSERM, 75010, Paris, France

    Pascal Richette & Thomas Bardin

  8. Department of Rheumatology, AP-HP, Lariboisière hospital, 2 rue A. Paré, 75010, Paris, France

    Pascal Richette & Thomas Bardin

  9. Université de Paris, ECSTRRA Team U1153, INSERM, 75010, Paris, France

    Matthieu Resche-Rigon

  10. Department of Biostatistics, AP-HP, Saint-Louis hospital, 75010, Paris, France

    Matthieu Resche-Rigon

Authors
  1. Minh Duc Do
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Thao Phuong Mai
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Anh Duy Do
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Quang Dinh Nguyen
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nghia Hieu Le
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Linh Gia Hoang Le
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Vu Anh Hoang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Anh Ngoc Le
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Hung Quoc Le
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Pascal Richette
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Matthieu Resche-Rigon
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Thomas Bardin
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MDD, MRR, PR, and TB designed the study. MDD, ADD, ANL, QDN, TB, and HQL included the patients. NHL and TB collected the anonymized data. MDD, TPM, LGHL, and VAH performed the HLA genotyping. MRR performed the statistical analysis. ADD and TB wrote the first version of the manuscript, which was read and approved by all the authors.

Corresponding author

Correspondence to Thomas Bardin.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the University of Medicine and Pharmacy at Ho Chi Minh City (HEC/IRB no. 198/ĐHYD-HĐ). Oral and written informed consent was obtained for all patients.

Consent for publication

Not applicable

Competing interests

None

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Primers used for PCR and sequencing HLA-B.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Do, M.D., Mai, T.P., Do, A.D. et al. Risk factors for cutaneous reactions to allopurinol in Kinh Vietnamese: results from a case-control study. Arthritis Res Ther 22, 182 (2020). https://doi.org/10.1186/s13075-020-02273-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-020-02273-1

Keywords

Arthritis Research & Therapy

ISSN: 1478-6362

Contact us