An Introduction to the Calibration of Computer Models

In the context of computer models, calibration is the process of estimating unknown simulator parameters from observational data. Calibration is variously referred to as model fitting, parameter estimation/inference, an inverse problem, and model tuning. The need for calibration occurs in most areas of science and engineering, and has been used to estimate hard to measure parameters in climate, cardiology, drug therapy response, hydrology, and many other disciplines. Although the statistical method used for calibration can vary substantially, the underlying approach is essentially the same and can be considered abstractly. In this survey, we review the decisions that need to be taken when calibrating a model, and discuss a range of computational methods that can be used to compute Bayesian posterior distributions

A Model to Investigate the Impact of Farm Practice on Antimicrobial Resistance in UK Dairy Farms

© 2021, The Author(s). The ecological and human health impact of antibiotic use and the related antimicrobial resistance (AMR) in animal husbandry is poorly understood. In many countries, there has been considerable pressure to reduce overall antibiotic use in agriculture or to cease or minimise use of human critical antibiotics. However, a more nuanced approach would consider the differential impact of use of different antibiotic classes; for example, it is not known whether reduced use of bacteriostatic or bacteriolytic classes of antibiotics would be of greater value. We have developed an ordinary differential equation model to investigate the effects of farm practice on the spread and persistence of AMR in the dairy slurry tank environment. We model the chemical fate of bacteriolytic and bacteriostatic antibiotics within the slurry and their effect on a population of bacteria, which are capable of resistance to both types of antibiotic. Through our analysis, we find that changing the rate at which a slurry tank is emptied may delay the proliferation of multidrug-resistant bacteria by up to five years depending on conditions. This finding has implications for farming practice and the policies that influence waste management practices. We also find that, within our model, the development of multidrug resistance is particularly sensitive to the use of bacteriolytic antibiotics, rather than bacteriostatic antibiotics, and this may be cause for controlling the usage of bacteriolytic antibiotics in agriculture

Protein kinase Cα (PKCα) Regulates Bone Architecture and Osteoblast Activity*

Bones' strength is achieved and maintained through adaptation to load bearing. The role of the protein kinase PKCα in this process has not been previously reported. However, we observed a phenotype in the long bones of Prkca−/− female but not male mice, in which bone tissue progressively invades the medullary cavity in the mid-diaphysis. This bone deposition progresses with age and is prevented by disuse but unaffected by ovariectomy. Castration of male Prkca−/− but not WT mice results in the formation of small amounts of intramedullary bone. Osteoblast differentiation markers and Wnt target gene expression were up-regulated in osteoblast-like cells derived from cortical bone of female Prkca−/− mice compared with WT. Additionally, although osteoblastic cells derived from WT proliferate following exposure to estradiol or mechanical strain, those from Prkca−/− mice do not. Female Prkca−/− mice develop splenomegaly and reduced marrow GBA1 expression reminiscent of Gaucher disease, in which PKC involvement has been suggested previously. From these data, we infer that in female mice, PKCα normally serves to prevent endosteal bone formation stimulated by load bearing. This phenotype appears to be suppressed by testicular hormones in male Prkca−/− mice. Within osteoblastic cells, PKCα enhances proliferation and suppresses differentiation, and this regulation involves the Wnt pathway. These findings implicate PKCα as a target gene for therapeutic approaches in low bone mass conditions

Lung function and microbiota diversity in cystic fibrosis

Abstract: Background: Chronic infection and concomitant airway inflammation is the leading cause of morbidity and mortality for people living with cystic fibrosis (CF). Although chronic infection in CF is undeniably polymicrobial, involving a lung microbiota, infection surveillance and control approaches remain underpinned by classical aerobic culture-based microbiology. How to use microbiomics to direct clinical management of CF airway infections remains a crucial challenge. A pivotal step towards leveraging microbiome approaches in CF clinical care is to understand the ecology of the CF lung microbiome and identify ecological patterns of CF microbiota across a wide spectrum of lung disease. Assessing sputum samples from 299 patients attending 13 CF centres in Europe and the USA, we determined whether the emerging relationship of decreasing microbiota diversity with worsening lung function could be considered a generalised pattern of CF lung microbiota and explored its potential as an informative indicator of lung disease state in CF. Results: We tested and found decreasing microbiota diversity with a reduction in lung function to be a significant ecological pattern. Moreover, the loss of diversity was accompanied by an increase in microbiota dominance. Subsequently, we stratified patients into lung disease categories of increasing disease severity to further investigate relationships between microbiota characteristics and lung function, and the factors contributing to microbiota variance. Core taxa group composition became highly conserved within the severe disease category, while the rarer satellite taxa underpinned the high variability observed in the microbiota diversity. Further, the lung microbiota of individual patient were increasingly dominated by recognised CF pathogens as lung function decreased. Conversely, other bacteria, especially obligate anaerobes, increasingly dominated in those with better lung function. Ordination analyses revealed lung function and antibiotics to be main explanators of compositional variance in the microbiota and the core and satellite taxa. Biogeography was found to influence acquisition of the rarer satellite taxa. Conclusions: Our findings demonstrate that microbiota diversity and dominance, as well as the identity of the dominant bacterial species, in combination with measures of lung function, can be used as informative indicators of disease state in CF. BBFJdPr3cu-jH3LTAhe361Video Abstrac

Disability, fatigue, pain and their associates in early diffuse cutaneous systemic sclerosis: the European Scleroderma Observational Study

Objectives; Our aim was to describe the burden of early dcSSc in terms of disability, fatigue and pain in the European Scleroderma Observational Study cohort, and to explore associated clinical features. Methods; Patients completed questionnaires at study entry, 12 and 24 months, including the HAQ disability index (HAQ-DI), the Cochin Hand Function Scale (CHFS), the Functional Assessment of Chronic Illness Therapy-fatigue and the Short Form 36 (SF36). Associates examined included the modified Rodnan skin score (mRSS), current digital ulcers and internal organ involvement. Correlations between 12-month changes were also examined. Results; The 326 patients recruited (median disease duration 11.9 months) displayed high levels of disability [mean (S.D.) HAQ-DI 1.1 (0.83)], with ‘grip’ and ‘activity’ being most affected. Of the 18 activities assessed in the CHFS, those involving fine finger movements were most affected. High HAQ-DI and CHFS scores were both associated with high mRSS (ρ = 0.34, P < 0.0001 and ρ = 0.35, P < 0.0001, respectively). HAQ-DI was higher in patients with digital ulcers (P = 0.004), pulmonary fibrosis (P = 0.005), cardiac (P = 0.005) and muscle involvement (P = 0.002). As anticipated, HAQ-DI, CHFS, the Functional Assessment of Chronic Illness Therapy and SF36 scores were all highly correlated, in particular the HAQ-DI with the CHFS (ρ = 0.84, P < 0.0001). Worsening HAQ-DI over 12 months was strongly associated with increasing mRSS (ρ = 0.40, P < 0.0001), decreasing hand function (ρ = 0.57, P < 0.0001) and increasing fatigue (ρ = −0.53, P < 0.0001). Conclusion; The European Scleroderma Observational Study highlights the burden of disability in early dcSSc, with high levels of disability and fatigue, associating with the degree of skin thickening (mRSS). Impaired hand function is a major contributor to overall disability

Treatment outcome in early diffuse cutaneous systemic sclerosis: the European Scleroderma Observational Study (ESOS).

OBJECTIVES: The rarity of early diffuse cutaneous systemic sclerosis (dcSSc) makes randomised controlled trials very difficult. We aimed to use an observational approach to compare effectiveness of currently used treatment approaches. METHODS: This was a prospective, observational cohort study of early dcSSc (within three years of onset of skin thickening). Clinicians selected one of four protocols for each patient: methotrexate, mycophenolate mofetil (MMF), cyclophosphamide or 'no immunosuppressant'. Patients were assessed three-monthly for up to 24 months. The primary outcome was the change in modified Rodnan skin score (mRSS). Confounding by indication at baseline was accounted for using inverse probability of treatment (IPT) weights. As a secondary outcome, an IPT-weighted Cox model was used to test for differences in survival. RESULTS: Of 326 patients recruited from 50 centres, 65 were prescribed methotrexate, 118 MMF, 87 cyclophosphamide and 56 no immunosuppressant. 276 (84.7%) patients completed 12 and 234 (71.7%) 24 months follow-up (or reached last visit date). There were statistically significant reductions in mRSS at 12 months in all groups: -4.0 (-5.2 to -2.7) units for methotrexate, -4.1 (-5.3 to -2.9) for MMF, -3.3 (-4.9 to -1.7) for cyclophosphamide and -2.2 (-4.0 to -0.3) for no immunosuppressant (p value for between-group differences=0.346). There were no statistically significant differences in survival between protocols before (p=0.389) or after weighting (p=0.440), but survival was poorest in the no immunosuppressant group (84.0%) at 24 months. CONCLUSIONS: These findings may support using immunosuppressants for early dcSSc but suggest that overall benefit is modest over 12 months and that better treatments are needed. TRIAL REGISTRATION NUMBER: NCT02339441

Oral abstracts 3: RA Treatment and outcomesO13. Validation of jadas in all subtypes of juvenile idiopathic arthritis in a clinical setting

Background: Juvenile Arthritis Disease Activity Score (JADAS) is a 4 variable composite disease activity (DA) score for JIA (including active 10, 27 or 71 joint count (AJC), physician global (PGA), parent/child global (PGE) and ESR). The validity of JADAS for all ILAR subtypes in the routine clinical setting is unknown. We investigated the construct validity of JADAS in the clinical setting in all subtypes of JIA through application to a prospective inception cohort of UK children presenting with new onset inflammatory arthritis. Methods: JADAS 10, 27 and 71 were determined for all children in the Childhood Arthritis Prospective Study (CAPS) with complete data available at baseline. Correlation of JADAS 10, 27 and 71 with single DA markers was determined for all subtypes. All correlations were calculated using Spearman's rank statistic. Results: 262/1238 visits had sufficient data for calculation of JADAS (1028 (83%) AJC, 744 (60%) PGA, 843 (68%) PGE and 459 (37%) ESR). Median age at disease onset was 6.0 years (IQR 2.6-10.4) and 64% were female. Correlation between JADAS 10, 27 and 71 approached 1 for all subtypes. Median JADAS 71 was 5.3 (IQR 2.2-10.1) with a significant difference between median JADAS scores between subtypes (p < 0.01). Correlation of JADAS 71 with each single marker of DA was moderate to high in the total cohort (see Table 1). Overall, correlation with AJC, PGA and PGE was moderate to high and correlation with ESR, limited JC, parental pain and CHAQ was low to moderate in the individual subtypes. Correlation coefficients in the extended oligoarticular, rheumatoid factor negative and enthesitis related subtypes were interpreted with caution in view of low numbers. Conclusions: This study adds to the body of evidence supporting the construct validity of JADAS. JADAS correlates with other measures of DA in all ILAR subtypes in the routine clinical setting. Given the high frequency of missing ESR data, it would be useful to assess the validity of JADAS without inclusion of the ESR. Disclosure statement: All authors have declared no conflicts of interest. Table 1Spearman's correlation between JADAS 71 and single markers DA by ILAR subtype ILAR Subtype Systemic onset JIA Persistent oligo JIA Extended oligo JIA Rheumatoid factor neg JIA Rheumatoid factor pos JIA Enthesitis related JIA Psoriatic JIA Undifferentiated JIA Unknown subtype Total cohort Number of children 23 111 12 57 7 9 19 7 17 262 AJC 0.54 0.67 0.53 0.75 0.53 0.34 0.59 0.81 0.37 0.59 PGA 0.63 0.69 0.25 0.73 0.14 0.05 0.50 0.83 0.56 0.64 PGE 0.51 0.68 0.83 0.61 0.41 0.69 0.71 0.9 0.48 0.61 ESR 0.28 0.31 0.35 0.4 0.6 0.85 0.43 0.7 0.5 0.53 Limited 71 JC 0.29 0.51 0.23 0.37 0.14 -0.12 0.4 0.81 0.45 0.41 Parental pain 0.23 0.62 0.03 0.57 0.41 0.69 0.7 0.79 0.42 0.53 Childhood health assessment questionnaire 0.25 0.57 -0.07 0.36 -0.47 0.84 0.37 0.8 0.66 0.4

Spondylarthropathies (including psoriatic arthritis): 244. Validity of Colour Doppler and Spectral Doppler Ultrasound of Sacroilicac Joints Againts Physical Examination as Gold Standard

Background: Sacroiliac joints (SJ) involvement is a distinctive and charasteristic feature of Spondyloarthritis (SpA) and x-ray is the test routinely used to make a diagnosis. However, x-ray reveals late structural damage but cannot detect active inflammation. The objective of this study was to assess the validity of Doppler ultrasound in SJ. Methods: Prospective blinded and controlled study of SJ, in which three populations were compared. We studied 106 consecutive cases, who were divided into three groups: a) 53 patients diagnosed with SpA who had inflammatory lumbar and gluteal pain assessed by a rheumatologist; b) 26 patients diagnosed with SpA who didn't have SJ tenderness and had normal physical examination; c) control group of 27 subjects (healthy subjetcs or with mechanical lumbar pain). All patients included that were diagnosed with SpA met almost the European Spondyloarthropathy Study Group (ESSG) classification criteria. Physical examination of the SJ included: sacral sulcus tenderness, iliac gapping, iliac compression, midline sacral thrust test, Gaenslen's test, and Patrick s test were used as gold standard. Both SJ were examined with Doppler ultrasound (General Electric Logiq 9, Wauwatosa WI, USA) fitted with a 9-14 Mhz lineal probe. The ultrasonographer was blinded to clinical data. Doppler in SJ was assessed as positive when both Doppler colour and resistance index (RI) < 0.75 within the SJ area were present. Statistical analysis was performed estimating sensitivity and specificity against gold standard. The Kappa correlation coefficient was used for reliability study. Results: 106 cases (53 female, 55 male; mean age 36 10 years) were studied. There were no statistical differences between groups related to age or sex. Physical examination of SJ was positive in 38 patients (59 sacroiliac joints). US detected Doppler signal within SJ in 37 patients (58 SJ): 33 of them were symptomatic SpA (52 SJ), one of them were asymptomatic SpA (1 SJ) and one was a healthy control (1 SJ). The accuracy of US when compared to clinical data as gold standard at subject level in the overall group was: sensitivity of 68.6% and specificity of 85.7%, positive predictive value of 70.5% and negative predictive value of 84.5%. A positive likelihood ratio of 4.8, a negative likelihood ratio of 0.36 and a kappa coefficient of 0.55 were achieved. Conclusions: Doppler US of SJ seems to be a valid method to detect active SJ inflammation. Disclosure statement: The authors have declared no conflicts of interes

Clinical Utility of Random Anti–Tumor Necrosis Factor Drug–Level Testing and Measurement of Antidrug Antibodies on the Long-Term Treatment Response in Rheumatoid Arthritis

Objective: To investigate whether antidrug antibodies and/or drug non-trough levels predict the long-term treatment response in a large cohort of patients with rheumatoid arthritis (RA) treated with adalimumab or etanercept and to identify factors influencing antidrug antibody and drug levels to optimize future treatment decisions.  Methods: A total of 331 patients from an observational prospective cohort were selected (160 patients treated with adalimumab and 171 treated with etanercept). Antidrug antibody levels were measured by radioimmunoassay, and drug levels were measured by enzyme-linked immunosorbent assay in 835 serial serum samples obtained 3, 6, and 12 months after initiation of therapy. The association between antidrug antibodies and drug non-trough levels and the treatment response (change in the Disease Activity Score in 28 joints) was evaluated.  Results: Among patients who completed 12 months of followup, antidrug antibodies were detected in 24.8% of those receiving adalimumab (31 of 125) and in none of those receiving etanercept. At 3 months, antidrug antibody formation and low adalimumab levels were significant predictors of no response according to the European League Against Rheumatism (EULAR) criteria at 12 months (area under the receiver operating characteristic curve 0.71 [95% confidence interval (95% CI) 0.57, 0.85]). Antidrug antibody–positive patients received lower median dosages of methotrexate compared with antidrug antibody–negative patients (15 mg/week versus 20 mg/week; P = 0.01) and had a longer disease duration (14.0 versus 7.7 years; P = 0.03). The adalimumab level was the best predictor of change in the DAS28 at 12 months, after adjustment for confounders (regression coefficient 0.060 [95% CI 0.015, 0.10], P = 0.009). Etanercept levels were associated with the EULAR response at 12 months (regression coefficient 0.088 [95% CI 0.019, 0.16], P = 0.012); however, this difference was not significant after adjustment. A body mass index of ≥30 kg/m2 and poor adherence were associated with lower drug levels.  Conclusion: Pharmacologic testing in anti–tumor necrosis factor–treated patients is clinically useful even in the absence of trough levels. At 3 months, antidrug antibodies and low adalimumab levels are significant predictors of no response according to the EULAR criteria at 12 months

Basic science232. Certolizumab pegol prevents pro-inflammatory alterations in endothelial cell function

Background: Cardiovascular disease is a major comorbidity of rheumatoid arthritis (RA) and a leading cause of death. Chronic systemic inflammation involving tumour necrosis factor alpha (TNF) could contribute to endothelial activation and atherogenesis. A number of anti-TNF therapies are in current use for the treatment of RA, including certolizumab pegol (CZP), (Cimzia ®; UCB, Belgium). Anti-TNF therapy has been associated with reduced clinical cardiovascular disease risk and ameliorated vascular function in RA patients. However, the specific effects of TNF inhibitors on endothelial cell function are largely unknown. Our aim was to investigate the mechanisms underpinning CZP effects on TNF-activated human endothelial cells. Methods: Human aortic endothelial cells (HAoECs) were cultured in vitro and exposed to a) TNF alone, b) TNF plus CZP, or c) neither agent. Microarray analysis was used to examine the transcriptional profile of cells treated for 6 hrs and quantitative polymerase chain reaction (qPCR) analysed gene expression at 1, 3, 6 and 24 hrs. NF-κB localization and IκB degradation were investigated using immunocytochemistry, high content analysis and western blotting. Flow cytometry was conducted to detect microparticle release from HAoECs. Results: Transcriptional profiling revealed that while TNF alone had strong effects on endothelial gene expression, TNF and CZP in combination produced a global gene expression pattern similar to untreated control. The two most highly up-regulated genes in response to TNF treatment were adhesion molecules E-selectin and VCAM-1 (q 0.2 compared to control; p > 0.05 compared to TNF alone). The NF-κB pathway was confirmed as a downstream target of TNF-induced HAoEC activation, via nuclear translocation of NF-κB and degradation of IκB, effects which were abolished by treatment with CZP. In addition, flow cytometry detected an increased production of endothelial microparticles in TNF-activated HAoECs, which was prevented by treatment with CZP. Conclusions: We have found at a cellular level that a clinically available TNF inhibitor, CZP reduces the expression of adhesion molecule expression, and prevents TNF-induced activation of the NF-κB pathway. Furthermore, CZP prevents the production of microparticles by activated endothelial cells. This could be central to the prevention of inflammatory environments underlying these conditions and measurement of microparticles has potential as a novel prognostic marker for future cardiovascular events in this patient group. Disclosure statement: Y.A. received a research grant from UCB. I.B. received a research grant from UCB. S.H. received a research grant from UCB. All other authors have declared no conflicts of interes