Effect of disease-modifying anti-rheumatic drugs on bone structure and strength in psoriatic arthritis patients
Arthritis Research & Therapy volume 21, Article number: 162 (2019)
To address whether the use of methotrexate (MTX) and biological disease-modifying anti-rheumatic drugs (bDMARDs) impacts bone structure and biomechanical properties in patients with psoriatic arthritis (PsA).
This is a cross-sectional study in PsA patients receiving no DMARDs, MTX, or bDMARDs. Volumetric bone mineral densities (vBMDs), microstructural parameters, and biomechanical properties (stiffness/failure load) were determined by high-resolution peripheral quantitative CT and micro-finite element analysis in the respective groups. Bone parameters were compared between PsA patients with no DMARDs and those receiving any DMARDs, MTX, or bDMARDs, respectively.
One hundred sixty-five PsA patients were analyzed, 79 received no DMARDs, 86 received DMARDs, of them 52 bDMARDs (TNF, IL-17- or IL-12/23 inhibitors) and 34 MTX. Groups were balanced for age, sex, comorbidities, functional index, and bone-active therapy, while disease duration was longest in the bDMARD group (7.8 ± 7.4 years), followed by the MTX group (4.6 ± 7.4) and the no-DMARD group (2.9 ± 5.2). No difference in bone parameters was found between the no-DMARD group and the MTX group. In contrast, the bDMARD group revealed significantly higher total (p = 0.001) and trabecular vBMD (p = 0.005) as well as failure load (p = 0.012) and stiffness (p = 0.012). In regression models, age and bDMARDs influenced total vBMD, while age, sex, and bDMARDs influenced failure load and stiffness.
Despite longer disease duration, bDMARD-treated PsA patients benefit from higher bone mass and better bone strength than PsA patients receiving MTX or no DMARDs. These data support the concept of better control of PsA-related bone disease by bDMARDs.
Psoriatic arthritis (PsA) is a chronic inflammatory joint disease associated with psoriasis characterized by peripheral arthritis and enthesitis leading to structural damage [1,2,3]. Bone erosions and enthesiophytes are hallmarks of local structural damage in PsA. More recently, the impact of PsA on systemic bone is increasingly appreciated. Hence, systemic bone loss has been documented to occur in PsA [4, 5] and increased prevalence of fractures in PsA patients is reported .
In contrast to RA, little is known about the effect of disease-modifying anti-rheumatic drugs (DMARDs) on bone structure in PsA. In rheumatoid arthritis (RA), biological DMARDs (bDMARDs) have shown to inhibit bone loss and thus may prevent pathological fractures [7,8,9]. To date, no such studies have been done in PsA; however, it can be assumed that effective control of inflammation may also impact secondary bone loss and bone biomechanics in PsA patients. In support of this notion, bDMARD treatment in PsA patients retards the progression of periarticular bone erosions [10,11,12] and periarticular bone loss [13,14,15]. Whether systemic bone mass and bone biomechanical properties are influenced by DMARD treatment is unclear to date. Furthermore, methotrexate treatment may not necessarily share potential beneficial effects of bDMARDs on bone in PsA patients.
While dual-energy X-ray absorptiometry (DXA) or digital X-ray radiogrammetry (DXR) can quantify bone loss, they do not allow separate assessment of changes in the cortical and trabecular bone compartment or the biomechanical properties of bone. High-resolution peripheral quantitative CT (HR-pQCT) enables the analysis of bone mass and microstructure and via integration of micro-finite element analysis (μFEA) the assessment of bone strength . An HR-pQCT study in RA patients has already shown that the biomechanical properties of the bone are reduced in RA patients .
To investigate whether DMARD treatment influences systemic bone structure and function in PsA patients, we investigated bone density, bone microstructure, and biomechanical properties in patients receiving either no DMARDs, methotrexate (MTX) treatment, or biologic DMARDs (bDMARDs) by HR-pQCT. In this cross-sectional study, we were specifically interested, whether PsA patients treated with either MTX or bDMARDs show a better bone structure and function than PsA patients receiving no DMARDs.
Psoriatic arthritis patients
PsA patients were part of the Erlangen Imaging Cohort (ERIC), which prospectively assesses bone composition in patients with inflammatory arthritis . All participants were consecutively recruited during routine diagnostic assessments at the Department of Internal Medicine 3 of the University of Erlangen-Nuremberg. Patients were examined by experienced rheumatologists (AK, JR, AJH) and had to fulfill the Classification criteria for Psoriatic Arthritis (CASPAR) . Age, sex, and smoking habits were recorded. With respect to psoriatic disease, duration of PsA, duration of psoriasis, minimal disease activity (MDA) state, disease activity in psoriatic arthritis (DAPSA), Psoriasis Area and Severity Index (PASI), dermatology life quality index (DLQI), health assessment questionnaire (HAQ), nail and/or scalp involvement, anti-rheumatic and bone-active medication, and laboratory parameters (rheumatoid factor, anti-citrullinated protein antibodies) were collected. Patients were categorized into three different treatment groups: (A) methotrexate monotherapy (over at least 6 months), (B) bDMARDs (without MTX) including TNF inhibitors, secukinumab and ustekinumab (over at least 6 months), and (C) a control group of PsA patients receiving no DMARDs over the last 6 months. The no-DMARD group included patients that for different reasons were not on DMARDs including (i) treatment-naïve patients at their first visit; (ii) treatment-naïve patients who had mild disease that was not considered to require immediate DMARD treatment; (iii) DMARD pre-exposed patients that had stopped treatment for compliance, intolerance, or lack of efficacy; and (iv) patients that were in drug-free remission. The study was conducted on approval of the local ethics committee of the University Clinic of Erlangen. Each individual provided informed consent.
HR-pQCT was performed at the distal radius (dominant hand) by an XtremeCT scanner (Scanco Medical, Brüttisellen, Switzerland). The following bone parameters were assessed: volumetric bone mineral density (vBMD) of the total, trabecular, meta-trabecular, inner trabecular, and cortical bone (all: mg HA/cm3) and ratio of meta-to-inner density (%) and cross-sectional bone area (mm2). Bone microstructure was assessed by determining trabecular bone volume fraction (%), trabecular number (1/mm), thickness (mm), separation (mm), network inhomogeneity (SD of 1/trabecular number (mm)), and cortical thickness (mm) [20,21,22].
Micro-finite element analysis
For micro-finite element analysis (μFEA), finite element analysis software (FAIM, version 8.0, Numerics88 solution, Calgary, Canada) was used. In order to generate micro-finite element models, the segmented trabecular network and cortex of the HR-pQCT images were used . Mesh size of the resulting models ranged from 1.5 to 3.5 million equally sized brick elements. Single linear isotropic tissue modeling was applied by assigning a tissue modulus of 6829 MPa and a Poisson’s ratio of 0.3 homogeneously to each element . A linear uniaxial compression test was simulated. Nodes on the proximal bone surface were fixed in z direction but unconstrained in x and y directions. Nodes on the distal bone surface were also free in the x and y directions but exposed to a displacement equivalent to 1% strain along the z axis . Axial bone stiffness (kN/mm) as reaction force (RFz) divided by average displacement of the distal surface (Uz) and bone strength as estimated failure load (N) based on the Pistoia criterion was calculated .
Statistical analyses compared PsA patients treated by bDMARDs or methotrexate to patients receiving no DMARD with respect to vBMD, bone microstructure, and biomechanical properties. Categorical variables are presented as numbers and percentages, and continuous variables as mean ± SD. Frequency distributions of categorical variables were compared using χ2 tests. Clinical, bone structural, and μFEA parameters were compared by using Kruskal-Wallis test (KW) with subsequent pairwise Mann-Whitney U tests, if KW test was significant. In order to account for multiple testing, we applied Bonferroni-Holm adjustment for pairwise comparisons. Critical p values for adjusted levels of significance are shown in the corresponding tables. Finally, we fitted linear regression models to contrast the differences in HR-pQCT and μFEA measurements with methotrexate and bDMARD use in comparison to no DMARDs. Models were adjusted for age, gender, and gender-treatment interaction. As a sensitivity analysis, we re-ran the models including MDA. All data manipulation and analyses were conducted using R (V3.5.1, R Foundation for Statistical Computing, Vienna, Austria.).
Characteristics of psoriatic arthritis patients
One hundred sixty-five PsA patients were included, 86 of them received DMARD treatment. Thirty-four PsA patients received MTX, and 52 bDMARDs for at least 6 months. Within the bDMARD group, 31 patients had TNF inhibitors (13 adalimumab, 6 infliximab, 9 etanercept, 2 certolizumab, and 1 golimumab), 16 patients received the IL-17A inhibitor secukinumab, and 5 patients the IL12/23 inhibitor ustekinumab. The mean duration of bDMARD treatment was 3.9 ± 3.3 years. Eleven patients in the bDMARD group had received previous MTX, while no bDMARD patient per definition received concurrent MTX. Seventy-nine PsA patients serving as the control group received no DMARDs. Detailed information on the demographic and disease-specific characteristics of the patients is shown in Table 1. Briefly, age and sex distribution, functional index, and comorbidities were not different among the no-DMARD, MTX, and bDMARD groups. The use of bone-active treatments such as vitamin D and anti-resorptive drugs was also balanced. A greater proportion of bDMARD-treated PsA patients received glucocorticoids compared to the no-DMARD group (p = 0.002). In addition, a larger proportion of patients under bDMARD therapy had received glucocorticoids in the past (p = 0.001) and hence had a longer duration of glucocorticoid therapy.
bDMARD-treated PsA patients had also the longest disease duration (7.6 ± 8.4 years; p < 0.001 compared to no DMARDs) followed by the MTX group (4.6 ± 7.4) and the no-DMARD group (2.9 ± 5.2). With respect to disease control, more of the patients receiving bDMARDs (48%) and MTX (65%) were in the MDA state than those in the no-DMARD group (35%).
Better bone microstructure and functional properties in PsA patients taking DMARDs
We first compared the no-DMARD control group with PsA patients taking any DMARDs (MTX or bDMARDs). DMARD patients had higher total vBMD (312 ± 53 vs. 290 ± 54, p = 0.004) and trabecular vBMD (171 ± 38 vs. 156 ± 39, p = 0.010) compared to no-DMARD controls. In addition, they had better bone microstructure indicated by higher number of trabeculae (2.09 ± 0.33 vs. 1.99 ± 0.35, p = 0.047), lower trabecular separation (0.43 ± 0.11 vs. 0.47 ± 0.18, p = 0.025), and higher cortical thickness (0.77 ± 0.17 vs. 0.71 ± 0.16, p = 0.0012). Regarding biomechanical properties, patients receiving DMARDs had higher stiffness and failure load (stiffness, 50.0 ± 15.0 vs. 45.2 ± 13.7, p = 0.034; failure load, 2385 ± 687 vs. 2154 ± 621, p = 0.026).
Better bone microstructure and functional properties is confined to PsA patients taking bDMARDs
To test whether the observed better bone status of DMARD-treated PsA patients is based on MTX or bDMARD treatment, we compared bone parameters between the no-DMARD control group and the MTX or the bDMARD group, respectively. The results demonstrated that MTX had no influence on bone microstructure and functional properties (Table 2). In contrast, and despite their longer disease duration, the bDMARD group exhibited significantly higher total and trabecular vBMD (320 ± 44 vs. 290 ± 54, p = 0.001; 174 ± 36 vs. 156 ± 39, p = 0.005) as compared to the no-DMARD group. In addition, higher cortical thickness (0.80 ± 0.15 vs. 0.71 ± 0.16, p = 0.001), with numerically higher numbers and thicker trabeculae (2.12 ± 0.32 vs. 1.99 ± 0.35, p = 0.022; 0.069 ± 0.010 vs. 0.065 ± 0.010, p = 0.030 (the adjusted significance was not met) were observed in the bDMARD group. Furthermore, also the biomechanical properties of bone (stiffness, 52.1 ± 15.0 vs. 45.2 ± 13.7, p = 0.012; failure load, 2473 ± 704 vs. 2154 ± 621, p = 0.012) were better in the bDMARD group than in the no-DMARD group (Table 2, Fig. 1).
Bone structure and function in PsA patients depend on age, sex, bDMARD therapy, and disease activity state
To test the impact of bDMARD treatment on bone structure and function, we set up three regression models with total vBMD, stiffness, and failure load as the respective outcome variable. We found that age and bDMARD treatment had a significant impact on vBMD, while age, sex, and bDMARD treatment influenced the model stiffness and failure load (Table 3). Use of bDMARD treatment, despite the longer disease duration and higher proportion of corticosteroid use, was associated with better bone density and higher stiffness and failure load estimates while no such association was observed with the use of methotrexate in comparison to no-DMARD use. When adding minimal disease activity (MDA) to the regression models, results remained robust (Additional file 1: Table S1).
The results of this study reveal that PsA patients receiving bDMARDs show better bone microstructure and biomechanical properties as compared to PsA patients receiving no DMARD treatment. In contrast, such differences are not found in PsA patients receiving MTX treatment suggesting that the beneficial bone effect of DMARD treatment in PsA is confined to the use of bDMARDs. These data are remarkable since bDMARD-treated PsA patients are a selectively more active and more resistant patient population, including a higher glucocorticoid use, which would reflect a higher burden of disease on the bone. However, we found the bone structure in bDMARD-treated PsA patients is not worse but even better that in the control population comprising patients with mild disease indicating a specific beneficial effect of bDMARDs on the bone. Notably, this effect was observed despite significantly longer disease duration in the bDMARD group.
These observations may be explained by previous functional data showing that the two central pro-inflammatory mediators in PsA, IL-17 and TNF, trigger an imbalance in bone homeostasis, increasing osteoclast-medicated bone resorption and inhibiting osteoblast-medicated bone formation [25,26,27,28]. This concept of cytokine-mediated bone loss in PsA is in fact supported by this study showing that PsA patients treated with either IL-17 inhibitor or TNFa inhibitors show better structural and functional bone data. We did not find significant differences in bone structure and function between IL-17 inhibitor- and TNFa inhibitor-treated PsA patients, suggesting a similar impact of the two main bDMARD treatment strategies on systemic bone in PsA.
Methotrexate monotherapy did not impact bone structure or function in PsA patients. Values for bone mass, microstructure, and function were consistently in the range of no-DMARD controls. This finding is interesting since the control of signs and symptoms of PsA was similar in the MTX- and bDMARD-treated groups. Hence, indirect effects such as better control of inflammation by bDMARDs are less likely to attribute for these differences. This notion is also supported by the fact that the inclusion of minimal disease activity in the regression models did not affect the results. On the other hand, methotrexate does not seem to share the positive effects of cytokine blockade on the bone. Hence, methotrexate has shown to inhibit osteoblast differentiation and bone formation , preventing that the anti-inflammatory effects of MTX are accompanied by sufficient rebalancing of the disturbed bone homeostasis.
A limitation of this study is the fact that it is cross-sectional and not longitudinal. Thus, we cannot conclude that bDMARDs increase bone mass and/or improve bone biomechanics in individual patients. Nonetheless, considering the longer disease duration of DMARD-treated than naïve PsA patients and the preferential use of bDMARDs in more severe PsA cases would suggest more severe rather than milder bone disease in bDMARD-treated patients. Since the exact opposite outcome was observed, the data support a direct effect of bDMARDs on bone structure and function in PsA. Furthermore, patient groups were balanced for several factors that could influence the bone including age and sex and comorbidities like diabetes mellitus, which have recently been shown to influence bone structure in PsA  as well as bone-active therapies. As a next step, a longitudinal study based on repeated HR-pQCT measurements will be necessary to better understand the effects of individual bDMARDs on bone homeostasis in PsA.
In summary, this study shows that the use of bDMARDs is associated with better bone structure and function in PsA patients. Given that PsA is associated with increased fracture risk, fast and adequate neutralization of the key pro-inflammatory and bone-destructive mediators seems to be important to restore bone health and to limit fracture risk in PsA patients.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Biological disease-modifying anti-rheumatic drugs
Classification criteria for Psoriatic Arthritis
Disease activity in psoriatic arthritis
Dermatology life quality index
Disease-modifying anti-rheumatic drugs
Dual-energy X-ray absorptiometry
Digital X-ray radiogrammetry
Health assessment questionnaire
High-resolution peripheral quantitative computed tomography
Minimal disease activity
- mg HA/cm3 :
Milligram hydroxyapatite/cubic centimeter
Psoriasis Area Severity Index
Tumor necrosis factor
Volumetric bone mineral density
Micro-finite element analysis
Schett G, Lories RJ, D'Agostino MA, Elewaut D, Kirkham B, Soriano ER, et al. Enthesitis: from pathophysiology to treatment. Nat Rev Rheumatol. 2017;13(12):731–41.
Ritchlin CT, Colbert RA, Gladman DD. Psoriatic arthritis. N Engl J Med. 2017;376(10):957–70.
Mease PJ. Biologic therapy for psoriatic arthritis. Rheum Dis Clin N Am. 2015;41(4):723–38.
Attia EA, Khafagy A, Abdel-Raheem S, Fathi S, Saad AA. Assessment of osteoporosis in psoriasis with and without arthritis: correlation with disease severity. Int J Dermatol. 2011;50(1):30–5.
Kocijan R, Englbrecht M, Haschka J, Simon D, Kleyer A, Finzel S, et al. Quantitative and qualitative changes of bone in psoriasis and psoriatic arthritis patients. J Bone Miner Res. 2015;30(10):1775–83.
Ogdie A, Harter L, Shin D, Baker J, Takeshita J, Choi HK, et al. The risk of fracture among patients with psoriatic arthritis and psoriasis: a population-based study. Ann Rheum Dis. 2017;76(5):882–5.
Hoff M, Kvien TK, Kalvesten J, Elden A, Kavanaugh A, Haugeberg G. Adalimumab reduces hand bone loss in rheumatoid arthritis independent of clinical response: subanalysis of the PREMIER study. BMC Musculoskelet Disord. 2011;12:54.
Vis M, Voskuyl AE, Wolbink GJ, Dijkmans BA, Lems WF. Bone mineral density in patients with rheumatoid arthritis treated with infliximab. Ann Rheum Dis. 2005;64(2):336–7.
Krieckaert CL, Nurmohamed MT, Wolbink G, Lems WF. Changes in bone mineral density during long-term treatment with adalimumab in patients with rheumatoid arthritis: a cohort study. Rheumatology (Oxford). 2013;52(3):547–53.
Kavanaugh A, Ritchlin C, Rahman P, Puig L, Gottlieb AB, Li S, et al. Ustekinumab, an anti-IL-12/23 p40 monoclonal antibody, inhibits radiographic progression in patients with active psoriatic arthritis: results of an integrated analysis of radiographic data from the phase 3, multicentre, randomised, double-blind, placebo-controlled PSUMMIT-1 and PSUMMIT-2 trials. Ann Rheum Dis. 2014;73(6):1000–6.
van der Heijde D, Kavanaugh A, Gladman DD, Antoni C, Krueger GG, Guzzo C, et al. Infliximab inhibits progression of radiographic damage in patients with active psoriatic arthritis through one year of treatment: results from the induction and maintenance psoriatic arthritis clinical trial 2. Arthritis Rheum. 2007;56(8):2698–707.
van der Heijde D, Landewe RB, Mease PJ, McInnes IB, Conaghan PG, Pricop L, et al. Brief report: secukinumab provides significant and sustained inhibition of joint structural damage in a phase III study of active psoriatic arthritis. Arthritis Rheumatol. 2016;68(8):1914–21.
Hoff M, Kavanaugh A, Haugeberg G. Hand bone loss in patients with psoriatic arthritis: posthoc analysis of IMPACT II data comparing infliximab and placebo. J Rheumatol. 2013;40(8):1344–8.
Kampylafka E, d'Oliveira I, Linz C, Lerchen V, Stemmler F, Simon D, et al. Resolution of synovitis and arrest of catabolic and anabolic bone changes in patients with psoriatic arthritis by IL-17A blockade with secukinumab: results from the prospective PSARTROS study. Arthritis Res Ther. 2018;20(1):153.
Szentpetery A, McKenna MJ, Murray BF, Ng CT, Brady JJ, Morrin M, et al. Periarticular bone gain at proximal interphalangeal joints and changes in bone turnover markers in response to tumor necrosis factor inhibitors in rheumatoid and psoriatic arthritis. J Rheumatol. 2013;40(5):653–62.
Macneil JA, Boyd SK. Bone strength at the distal radius can be estimated from high-resolution peripheral quantitative computed tomography and the finite element method. Bone. 2008;42(6):1203–13.
Stemmler F, Simon D, Liphardt AM, Englbrecht M, Rech J, Hueber AJ, et al. Biomechanical properties of bone are impaired in patients with ACPA-positive rheumatoid arthritis and associated with the occurrence of fractures. Ann Rheum Dis. 2018;77(7):973–80.
Simon D, Kleyer A, Englbrecht M, Stemmler F, Simon C, Berlin A, et al. A comparative analysis of articular bone in large cohort of patients with chronic inflammatory diseases of the joints, the gut and the skin. Bone. 2018;116:87–93.
Taylor W, Gladman D, Helliwell P, Marchesoni A, Mease P, Mielants H. Classification criteria for psoriatic arthritis: development of new criteria from a large international study. Arthritis Rheum. 2006;54(8):2665–73.
Simon D, Kleyer A, Stemmler F, Simon C, Berlin A, Hueber AJ, et al. Age- and sex-dependent changes of intra-articular cortical and trabecular bone structure and the effects of rheumatoid arthritis. J Bone Miner Res. 2017;32(4):722–30.
Khosla S, Riggs BL, Atkinson EJ, Oberg AL, McDaniel LJ, Holets M, et al. Effects of sex and age on bone microstructure at the ultradistal radius: a population-based noninvasive in vivo assessment. J Bone Miner Res. 2006;21(1):124–31.
Boutroy S, Bouxsein ML, Munoz F, Delmas PD. In vivo assessment of trabecular bone microarchitecture by high-resolution peripheral quantitative computed tomography. J Clin Endocrinol Metab. 2005;90(12):6508–15.
van Rietbergen B, Weinans H, Huiskes R, Odgaard A. A new method to determine trabecular bone elastic properties and loading using micromechanical finite-element models. J Biomech. 1995;28(1):69–81.
Pistoia W, van Rietbergen B, Lochmuller EM, Lill CA, Eckstein F, Ruegsegger P. Estimation of distal radius failure load with micro-finite element analysis models based on three-dimensional peripheral quantitative computed tomography images. Bone. 2002;30(6):842–8.
Adamopoulos IE, Chao CC, Geissler R, Laface D, Blumenschein W, Iwakura Y, et al. Interleukin-17A upregulates receptor activator of NF-kappaB on osteoclast precursors. Arthritis Res Ther. 2010;12(1):R29.
Gravallese EM, Schett G. Effects of the IL-23-IL-17 pathway on bone in spondyloarthritis. Nat Rev Rheumatol. 2018;14(11):631–40.
Lam J, Takeshita S, Barker JE, Kanagawa O, Ross FP, Teitelbaum SL. TNF-alpha induces osteoclastogenesis by direct stimulation of macrophages exposed to permissive levels of RANK ligand. J Clin Invest. 2000;106(12):1481–8.
Sato K, Suematsu A, Okamoto K, Yamaguchi A, Morishita Y, Kadono Y, et al. Th17 functions as an osteoclastogenic helper T cell subset that links T cell activation and bone destruction. J Exp Med. 2006;203(12):2673–82.
Uehara R, Suzuki Y, Ichikawa Y. Methotrexate (MTX) inhibits osteoblastic differentiation in vitro: possible mechanism of MTX osteopathy. J Rheumatol. 2001;28(2):251–6.
Zhu TY, Griffith JF, Qin L, Hung VW, Fong TN, Au SK, et al. Density, structure, and strength of the distal radius in patients with psoriatic arthritis: the role of inflammation and cardiovascular risk factors. Osteoporos Int. 2015;26(1):261–72.
This study was supported by the Deutsche Forschungsgemeinschaft (DFG-SFB1181-A01; GS, DFG - FOR2886 PANDORA Z-project ), the ERC Synergy Grant 4D Nanoscope (ERC-Syg 810316 4DnanoSCOPE), the Bundesministerium für Bildung und Forschung (BMBF; project METARTHROS TP 2), the EU IMI2 funded project RTCure, the EU IMI2 funded project PREFER and FOREUM (AH).
Ethics approval and consent to participate
Ethical approval was obtained from the Institutional Review Board of the Universitätsklinikum Erlangen. Patients provided written informed consents.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Simon, D., Kleyer, A., Bayat, S. et al. Effect of disease-modifying anti-rheumatic drugs on bone structure and strength in psoriatic arthritis patients. Arthritis Res Ther 21, 162 (2019). https://doi.org/10.1186/s13075-019-1938-3