Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Occurrence of Anti-Drug Antibodies against Interferon-Beta and Natalizumab in Multiple Sclerosis: A Collaborative Cohort Analysis

Abstract

Immunogenicity of biopharmaceutical products in multiple sclerosis is a frequent side effect which has a multifactorial etiology. Here we study associations between anti-drug antibody (ADA) occurrence and demographic and clinical factors. Retrospective data from routine ADA test laboratories in Sweden, Denmark, Austria and Germany (Dusseldorf group) and from one research study in Germany (Munich group) were gathered to build a collaborative multi-cohort dataset within the framework of the ABIRISK project. A subset of 5638 interferon-beta (IFNβ)-treated and 3440 natalizumab-treated patients having data on at least the first two years of treatment were eligible for interval-censored time-to-event analysis. In multivariate Cox regression, IFNβ-1a subcutaneous and IFNβ-1b subcutaneous treated patients were at higher risk of ADA occurrence compared to IFNβ-1a intramuscular-treated patients (pooled HR = 6.4, 95% CI 4.9–8.4 and pooled HR = 8.7, 95% CI 6.6–11.4 respectively). Patients older than 50 years at start of IFNβ therapy developed ADA more frequently than adult patients younger than 30 (pooled HR = 1.8, 95% CI 1.4–2.3). Men developed ADA more frequently than women (pooled HR = 1.3, 95% CI 1.1–1.6). Interestingly we observed that in Sweden and Germany, patients who started IFNβ in April were at higher risk of developing ADA (HR = 1.6, 95% CI 1.1–2.4 and HR = 2.4, 95% CI 1.5–3.9 respectively). This result is not confirmed in the other cohorts and warrants further investigations. Concerning natalizumab, patients older than 45 years had a higher ADA rate (pooled HR = 1.4, 95% CI 1.0–1.8) and women developed ADA more frequently than men (pooled HR = 1.4, 95% CI 1.0–2.0). We confirmed previously reported differences in immunogenicity of the different types of IFNβ. Differences in ADA occurrence by sex and age are reported here for the first time. These findings should be further investigated taking into account other exposures and biomarkers.

Citation: Bachelet D, Hässler S, Mbogning C, Link J, Ryner M, Ramanujam R, et al. (2016) Occurrence of Anti-Drug Antibodies against Interferon-Beta and Natalizumab in Multiple Sclerosis: A Collaborative Cohort Analysis. PLoS ONE 11(11): e0162752. https://doi.org/10.1371/journal.pone.0162752

Editor: Sreeram V. Ramagopalan, University of Oxford, UNITED KINGDOM

Received: June 16, 2016; Accepted: August 26, 2016; Published: November 2, 2016

Copyright: © 2016 Bachelet et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: ABIRISK does not own the data analyzed in this study, which were collected in the context of other past and current projects by both ABIRISK partners and non-ABIRISK research groups in the different European countries and provided to ABIRISK exclusively for the purpose of the performed analyses upon official request and ethical approval. For this reason access to the minimal dataset underlying the findings cannot be made publicly available, but can be obtained by interested researchers upon request to the ABIRISK steering committee by submission of an analysis plan. The analysis plan should explain the purpose of the use of the data and confirm the intention to use the data only for replication studies concerning anti-drug antibodies, since this is the limitation of the ethical permission on how this data can be used. The contact persons of the ABIRISK steering committee to whom the requests should be sent are Pierre Dönnes (pierre@scicross.com) and Sophie Tourdot (sophie.tourdot@inserm.fr).

Funding: The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement no [115303], resources of which are composed of financial contribution from the European Union's Seventh Framework Programme (FP7/2007-2013) and EFPIA companies. A. Lawton, N. Donnellan, A. Hincelin-Mery, D. Sikkema, J. Davidson: belong to EFPIA (European Federation of Pharmaceutical Industries and Association) member companies in the IMI JU and costs related to their part in the research were carried by the respective companies as in kind contributions under the IMI JU scheme. Bernard Hemmer and Dorothea Buck were supported by the German Competence Network for multiple sclerosis.

Competing interests: D. Bachelet, S. Hässler, C. Mbogning, J. Link, R. Ramanujam, M. Auer, P.E. Hyldgaard Jensen, K. Ingenhoven, D. Buck, V. Grummel, M. Pallardy, P. Dönnes, P. Broët: have nothing to declare. N. Koch-Henriksen has received honoraria for lecturing and participation in advisory councils, travel expenses for attending congresses and meetings, and financial support for monitoring the Danish Multiple Sclerosis Treatment Register from Bayer-Schering, MerckSerono, Biogen Idec, Teva, Sanofi-Aventis and Novartis. A. Fogdell-Hahn has received funding speaking honoraria from Biogen Idec and Pfizer. F. Deisenhammer participated in meetings sponsored by or received honoraria for acting as an advisor/speaker for Bayer Healthcare, Biogen Idec, Genzyme-Sanofi, Merck, Novartis Pharma, and Roche. C. Warnke participated in meetings sponsored by, or received honoraria for acting as an advisor/speaker for Bayer, Biogen, Novartis, and Teva. D. Buck has received compensation for activities with Bayer HealthCare, BiogenIdec, MerckSerono, and Novartis and she is supported by the Abirisk Consortium. H. Hartung has received honoraria for consulting, serving on steering committees and speaking from Biogen, GeNeuro, Genzyme, Merck, Novartis, Opexa, Receptos, Roche, Sanofi, Teva with approval by the president of Heinrich-Heine University. B. Kieseier has received honoraria for lecturing, travel expenses for attending meetings, and financial support for research from Bayer Health Care, Biogen Idec, Genzyme/Sanofi Aventis, Grifols, Merck Serono, Mitsubishi Europe, Novartis, Roche, Talecris, and TEVA. B. Hemmer has served on scientific advisory boards for Roche, Novartis, Bayer Schering, Merck Serono, Biogen Idec, GSK, Chugai Pharmaceuticals, Micromet, Genentech and Genzyme Corporation; serves on the international advisory board of Archives of Neurology and Experimental Neurology; has received speaker honoraria from Bayer Schering, Novartis, Biogen Idec, Merck Serono, Roche, and Teva Pharmaceutical Industries Ltd; and has received research support from Biogen Idec, Bayer Schering, Merck Serono, Five prime, Metanomics, Chugai Pharmaceuticals and Novartis. He has been filed a patent for the detection of antibodies and T cells against KIR4.1 in a subpopulation of MS patients and genetic determinants of neutralizing antibodies to interferon-beta. P.S. Sørensen has served on scientific advisory boards for Biogen Idec, Merck Serono, Novartis, Genzyme, Teva Pharmaceutical Industries Ltd., GlaxoSmithKline, medDay Pharmaceuticals and Forward Pharma; has been on steering committees or independent data monitoring boards in clinical trials sponsored by Merck Serono, Teva Pharmaceutical Industries Ltd., and GlaxoSmithKline; and has received speaker honoraria from Biogen Idec, Merck Serono, Teva Pharmaceutical Industries Ltd., Genzyme, and Novartis. His department has received research support from Biogen Idec, Bayer Schering, Merck Serono, TEVA, Baxter, Sanofi-Aventis, BioMS, Novartis, Bayer, RoFAR, Roche, Genzyme, from the Danish Multiple Sclerosis Society, the Danish Medical Research Council, and the European Union Sixth Framework Programme: Life sciences, Genomics and Biotechnology for health. M. Ryner has received research support from Biogen Idec and Sanofi-Aventis, and received speaker honoraria from Biogen Idec. Andy Lawton and Dan Sikkema are employed by GlaxoSmithKline. At the time of writing Julie Davidson was employed by and held stocks/shares in GlaxoSmithKline. Naoimh Donnellan is employed by IPSEN. Agnès Hincelin-Mery is employed by Sanofi. Pierre Dönnes is an employee of SciCross AB and has received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement No. [115303], resources of which are composed of financial contribution from the European Union's Seventh Framework Program (FP7/2007–2013) and EFPIA companies' in kind contribution. All the company-employed authors declare that this does not alter their adherence to PLOS ONE policies on sharing data and materials.

Introduction

The introduction of biotechnology-derived proteins (BPs) has been a critical step forward in the treatment of multiple sclerosis (MS), with BPs such as interferon beta (IFNβ) registered since the 1990s and natalizumab available since 2006 in the EU. Natalizumab is mainly used as a second line treatment for patients with relapsing-remitting MS (RRMS) that do not respond to treatment with first line drugs such as IFNβ. Failure of response to therapy may occur due to an unwanted immune response against the BP. This immunological reaction to BPs is not fully understood and the fraction of patients developing anti-drug antibodies (ADA) has varied from 2 to 47%, depending upon either the product or the laboratory methods and cut-off used to define antibody positivity [1].

The consequences of this immunological reaction range from transient ADA to high titer ADA with clear clinical effects [2,3]. In particular, the occurrence of ADA may decrease the efficacy of BPs by neutralizing their activity or modifying their drug clearance and ADA may also be associated with BP-specific hypersensitivity reactions [4]. Clinically relevant neutralizing ADA (NAbs) against IFNβ usually develop between 9 and 18 months after start of therapy [5]. Patients who have been persistently NAb negative during the first 2 years of IFNβ therapy only rarely become NAb positive [6]. For natalizumab, the antibodies are detected in approximatively 9% of treated patients with 3% that are transient [7]. These ADA develop quickly, within 4 months following the first infusion of treatment [8]. In establishing the used ADA assay, the cut-point set for ADA positivity by enzyme-linked immunosorbent assay (ELISA) correspond to neutralizing antibody positive by flow cytometry blocking assay [7].

Many factors (patient-, disease- or product-related) may influence the risk of BP immunogenicity but the relative contribution of these factors to the development of ADA is not currently understood. The ability to identify patients at high risk for ADA and subsequent lack of treatment efficacy or hypersensitivity reactions would allow better personalized therapy. The immunogenic potential of BPs can only be definitively assessed in human studies. Previous investigations have identified oligoclonal band negativity as a protective factor [9], specific HLA haplotypes as risk factors for IFNβ ADA [10–13], and smoking as a risk factor for both IFNβ and natalizumab ADA [7–9].

In order to carry out an unbiased evaluation of the relationships between bio-clinical factors and the occurrence of ADAs, we selected a target population that was extracted from the whole ABIRISK tranSMART database based on five European MS cohorts and that is composed by BP naive adult patients with biological samples being taken within the first two years of treatment.

The present study gathers data from four European central laboratories and one research center routinely performing ADA tests in order to identify patient-related and treatment-related factors associated with ADA development in IFNβ- and natalizumab-treated MS patients using a time-to-event analysis.

Patients and Methods

Patients

The present study is a collaborative analysis of historical data collected from a number of European MS cohorts under the leadership of the ABIRISK EU-consortium, standardized into a common format and grouped into the ABIRISK tranSMART database including data from national routine ADA test laboratories in Sweden, Denmark, Austria and Germany (Dusseldorf group) and from one research study in Germany (Munich group).

In Sweden, Germany and Denmark an ethical permission for the analysis of anonymized patient retrospective data was required. For the Swedish cohort the study was approved by the Regional Ethics Committee in Stockholm, Sweden (http://www.epn.se/en/start/the-organisation/), and approval for export of data from the national Multiple Sclerosis Registry (MSreg) was given by the Research Board of the Swedish MS Society. For the German cohorts, the study was approved by the Ethics Committee of the School of Medicine of the Technical University of Munich and the Ethics Committee of the medical faculty of the University of Dusseldorf. The Danish MS-cohort collected in the Danish Multiple Sclerosis Biobank was approved by the Regional Scientific Ethics Committee in Copenhagen and Frederiksberg, Denmark, and approval for use of data was given by the Danish Data Protection Agency in the Capital Region, Denmark. In Austria the Ethics Committee judged that for retrospective analyses based on anonymized data that were collected during routine procedures, an ethical approval was not required.

In total the database contained 20,465 patients from the four countries and ADA data of 42,843 samples together with treatment history, type of ADA assay, age and sex. Data about anti-natalizumab ADA tests were available only in reference laboratories from Sweden, Austria and Denmark.

Populations eligible for inclusion in the IFNβ and natalizumab historical longitudinal analyses were selected based on the following criteria: (i) for IFNβ, naive to IFNβ prior to the IFNβ for which the first NAb test is available, (ii) for IFNβ and natalizumab, date of first BP dose administration known (at least with a month precision), (iii) for IFNβ, at least one NAb test result available after first dose and before 24 months of treatment, (iv) for natalizumab, at least one ADA test result available after baseline and before 12 months of treatment, and (v) for IFNβ and natalizumab, age at date of first BP dose more than 18 years.

Due to the low sample size obtained from the Swiss and Spanish cohorts, we focused on the cohorts from Sweden, Austria, Denmark, and Germany (Dusseldorf and Munich groups). IFNβ-treated MS patients tested for ADA from 2002 to 2013 in Sweden, from 1996 to 2014 in Austria, from 2009 to 2014 in Denmark, from 2008 to 2014 in the German reference laboratory (Dusseldorf group) and from 2005 to 2008 in the German research study (Munich group) were included in the time-to-event analysis. Eligible natalizumab-treated patients were tested from 2006 to 2013 in Sweden, from 2007 to 2014 in Austria and from 2007 to 2014 in Denmark.

The ABIRISK tranSMART database

The ABIRISK database integrates immunogenicity data from a number of European countries in a well-structured format. The data is harmonized over different cohorts and the format is based on the CDISC standards (http://www.cdisc.org/). Where no previous variable description can be found in CDISC, a local variable description is used. Data was prepared by data custodians using a data load plan describing the variable semantics and format. Anonymized data was uploaded into the ABIRISK database, which is based on the tranSMART platform [14]. The tranSMART platform is an open-source knowledge management platform for translational science, supported by a large number of organizations (http://www.transmartfoundation.org/). tranSMART provides tools to explore data, select specific cohorts and do basic statistics. It is also possible to export data for further analysis in tools like SAS and R.

Biological samples testing

The reference laboratories received serum samples for analysis of neutralizing antibodies.

According to the countries, ADA screening was either mandatory for all BP-treated MS patients as in Denmark since 1996, sponsored by the pharmaceutical companies and so widely performed as in Sweden from 2003 to 2006 (thereafter performed ad hoc at the request of the neurologists) or always performed ad hoc at the request of the neurologists as in Austria and Germany (Dusseldorf). In the Munich study, ADA testing was performed in one single instance for research purposes during an outpatient visit in the period 2003 to 2008.

In each laboratory, different methods to analyze ADA against IFNβ were used. Some laboratories search for both binding ADA and NAbs (Austria, Germany, Denmark) while others only tested for NAbs (Sweden). Since these latter are more frequently persistent and clinically relevant than binding ADA, in the present study only ADA determined to be NAbs against IFNβ were considered.

In Sweden, three different NAb assays were performed across time: patient samples sent to the laboratory in 2003–2006 were analyzed with the myxovirus resistance protein A (MxA) assay (MPA) [15] while later samples were analyzed with the MxA gene expression assay (MGA) [16,17]. Since 2012, samples were analyzed with the commercially available iLite anti-Human IFNβ-1a bioassay (Biomonitor/Eurodiagnostica). In a recent publication a strong correlation of NAb titers was found on a series of 118 serum samples titrated with the MPA and the MGA assay in the Swedish central laboratory [12]. There were also significant correlations between NAb titer levels measured with MGA and iLite [18]. For the three assays, titers were adjusted using the Kawade formula [19]. NAb-positive screening samples were further titrated and NAb titers were expressed as 10-fold reduction units per milliliter (TRU/ml) and the patients were classified according to the following categories: negative (<10 TRU/ml), low positive (≥ 10–50 TRU/ml), medium positive (>50–200 TRU/ml), and high positive (>200 TRU/ml).

In Austria, three different NAb assays were also performed across time: the MPA and the MGA assays, as well as the Luciferase assay (LUC) where NAb were measured using human fibrosarcoma cells that express IFN-β receptors on their surface and were stably transfected with a luciferase reporter gene cassette [20,21]. The cut-off for positivity for all NAb-tests (MPA, MGA and LUC) was 20 TRU/ml, whereby 20–99 TRU/ml were considered low positive and titers beginning from 100 TRU/ml high positive.

In Denmark, only results from the LUC method were applied. NAb-positivity was defined by a titer ≥ 20 TRU/ml.

In Dusseldorf, ADA were first screened in the ELISA. Then NAbs were usually only tested with LUC in case of a positive ADA testing. The patients were classified according to the following categories: negative (<20 TRU/ml), low positive (≥ 20–100 TRU/ml) and high positive (>100 TRU/ml).

In Munich, ADA were first tested using a capture ELISA as described by Pachner et al. [22]. Next antibodies were defined as NAb if MxA induction was decreased by more than 50% compared with newly treated antibody-negative control donors. A cross-validation study on 247 samples was performed between Austria and Munich and has shown an agreement of 96% for classifying samples in NAb positive (>30 TRU/ml) versus NAb negative (personal communication from F. Deisenhammer).

In our analyses, to define clinically relevant NAb-positivity against IFNβ, we kept the cut-off used in each reference laboratory (≥200 TRU/ml in Sweden, ≥100 TRU/ml in Austria, ≥20 TRU/ml in Denmark, ≥100 TRU/ml in Germany–Dusseldorf, and >50% inhibition in Munich).

Concerning natalizumab, the same assay method undergoing regular proficiency testing was used in Sweden, Austria and Denmark. Antibodies were measured using a bridging ELISA method as described in the Biogen Idec protocol [7].

For ease of reading, we use "ADA" as a generic term in the remainder of the text for either NAbs concerning IFNβ-treated patients or ELISA-detected ADA concerning natalizumab-treated patients.

ADA outcome

In this work, we focused on a two-year ADA occurrence window after the start of therapy. To take into account that ADA occurrence is a dynamic event which appears mainly within this two-year window but can be censored, we performed time-to-event analyses. The outcome was the time between the date of the first BP dose and the first time of ADA occurrence (positive ADA). Patients without ADA occurrence were censored at the date of their last clinic visit, the end of treatment or at the last visit before a drug switch.

In practice, the ADA occurrence process was not under continuous monitoring. As a consequence, a positive ADA result indicated that ADA have been produced during a particular time interval but the exact event time of the ADA occurrence was unknown. We only knew that it occurred between two time points: the time of the last negative laboratory test and the time of the first positive laboratory test. For the analyses, we considered interval-censored survival methods that explicitly took into account the information obtained for each patient.

Statistical analysis

For both therapies IFNβ and natalizumab, we present results of separate analyses in each cohort and pooled analysis.

The non-parametric estimator for the cumulative probability distribution function under interval censoring was computed as the complement of the survival estimator obtained using the self-consistency or Turnbull’s algorithm [23]. The cumulative probability distribution function curves were plotted for each BP (IFNβ and natalizumab separately) and according to categorical variables (type of therapy, age group, gender, month of first treatment and month of birth). In order to compare these curves, p-values for rank-based tests adapted for interval censoring [23–25] are presented.

For the univariate and multivariate analyses, we considered the Cox proportional hazards model via a multiple imputation strategy for the unobserved survival times as proposed by Pan [26].

Univariate Cox model analyses were carried out to assess the influence of clinical variables on the time of occurrence of ADA. Hazard ratios (HR) and 95% confidence intervals (CIs) were determined and reported. We investigated the impact of age, sex, month of therapy start and month of birth. Previous IFNβ therapy before natalizumab were also explored (last IFNβ taken, first IFNβ taken, last ADA status during IFNβ therapy, total duration of IFNβ therapy and number of different IFNβ therapies).

All variables associated with the time to ADA occurrence in univariate analysis (with p-value less than 0.20) were considered for the multivariate analysis using the Cox proportional hazards model. Hazard ratios (HR) and 95% confidence intervals (CIs) were determined and reported.

Analyses were performed using R software (version 3.0.2 software) with "interval", "Icens" and "MIICD" packages.

Results

The flowchart (Fig 1) gives the details of the patients eligible in each cohort according to the chosen inclusion criteria for IFNβ and natalizumab.

thumbnail
Download:
Fig 1. Inclusion flow chart.

Flow chart of patients by country for IFNβ (left side) and natalizumab (right side) treated patients.

https://doi.org/10.1371/journal.pone.0162752.g001

Among the 17,631 available IFNβ-treated patients in the database, 5,638 eligible patients (32%) were considered for the survival analysis. For natalizumab, the pooled eligible cohort consisted of 3,440 patients (85% of the patients available in the database).

We compared included patients treated with IFNβ to those not included with respect to age at therapy start, sex, type of first IFNβ and ADA status for the first test performed. Regarding the BPs, the fraction of included patients was higher for IFNβ-1a intramuscular (i.m.) as compared to the other type of treatment (S5 Table). There was also a significant difference in age at start of therapy in Germany-Dusseldorf (fewer patients older than 40 were included) and ADA status at first visit differed between included and excluded in Sweden, Austria and Germany-Dusseldorf (S5 Table). Concerning natalizumab-treated patients, a lower proportion was excluded. They did not differ from included patients, except for ADA status at first test available in Sweden: more excluded patients were ADA positive. In most of the cases, this first test is a pre-treatment sample. All the patients with baseline ADA positive tests were excluded from analysis.

In this selected historical cohort, the median number of tests for ADA against IFNβ done for one patient was 2 (max = 7) in Sweden, 2 (max = 41) in Austria, 2 (max = 7) in Denmark, 1 (max = 7) in Germany—Dusseldorf group, and 1 in Germany–Munich group (current status data). Concerning natalizumab, the median number of ADA tests done for one patient was 4 (max = 10) in Sweden, 2 (max = 50) in Austria, and 3 (max = 9) in Denmark.

In the reference laboratories, the first IFNβ ADA test was performed around one year after first dose of therapy in Sweden, Austria and Germany, and 6 months after first dose of therapy in Denmark. Thereafter tests were done approximatively every 6 months in each country.

Concerning natalizumab, time interval between two visits was approximatively from 3 to 6 months.

Univariate analyses

Results for IFNβ.

According to the univariate analyses, the type of IFNβ treatment was strongly associated with the rate of ADA occurrence (p < 0.0001 in each cohort) (Fig 2).

thumbnail
Download:
Fig 2. Cumulative probability of developing anti-IFNβ ADA.

Estimated cumulative probability of developing anti-IFNβ ADA by cohort and by categorical variables with p-values for rank-based tests adapted for interval censoring.

https://doi.org/10.1371/journal.pone.0162752.g002

As expected most of the ADA occurred within two years of IFNβ therapy. Of the two different routes of administration for IFNβ-1a, the one with a single weekly intramuscular (i.m.) dose was less immunogenic with 1.1 to 4.2% (95% confidence intervals (CI) ranging from 0 to 6.8%) of the patients who developed ADA before 18 months of therapy as compared to 11.1 to 23.0% (95% CI ranging from 1.0 to 30.7%) for the patients with three weekly subcutaneous (s.c.) doses (Table 1). Patients treated with IFNβ-1b s.c. every other day had an even higher rate of ADA development before 18 months of therapy: 17.4 to 33.3% (95% CI ranging from 0 to 51.7%).

thumbnail
Download:
Table 1. Estimated cumulative probabilities with 95% confidence interval of ADA development before 18 months of IFNβ therapy.

https://doi.org/10.1371/journal.pone.0162752.t001

Looking to all the cumulative probability curves across the countries, we observed very similar rates of ADA development by type of IFNβ-therapy.

Age at start of IFNβ-therapy was associated with immunogenicity in each country (p < 0.01) with the exception of Germany (Dusseldorf group: p = 0.45; Munich group: p = 0.08). Older patients had an increased risk of ADA. In Austria, Germany and Denmark, male patients had an immune response to IFNβ more frequently than females, this association between sex and ADA being significant in Austria and in Germany–Munich group (p < 0.01), and borderline significant for Denmark (p = 0.06). This relationship did not reach significance for Germany–Dusseldorf group (p = 0.09) and Sweden (p = 0.14). When comparing by period of treatment start, patients who started IFNβ therapy in April were at higher risk of developing ADA compared to patients who started at any other month in Germany (p < 0.0001 in Dusseldorf-group, and p = 0.07 in Munich-group). This association was not significant in Sweden (p = 0.15) and in Denmark (p = 0.48). Conversely, in Austria, patients who started IFNβ therapy in April were at lower risk of developing ADA (p = 0.04).

The month of birth was not associated with ADA development in any of the studied cohorts.

The ADA assay methods which have changed across time in Sweden and Austria were associated with ADA occurrence, suggesting a difference in the sensitivity of the assays as previously reported [18]. The MPA method detected a higher rate of ADA with 20% (95% CI 14.7–25.1) of the patients who developed ADA before 18 months of therapy in Austria and 30.3% (95% CI 18.6–36.1) in Sweden, while 7% (95% CI 1.3–10.7) developed ADA with the MGA method in Austria and 8.4% (95% CI 6.9–10.6) in Sweden. Concerning the more recent assay methods, 11% (95% CI 8.0–14.6) of the patients developed ADA before 18 months with the LUC method in Austria and 4% (95% CI 0.7–7.4) with the iLite method in Sweden.

Results for natalizumab.

In the univariate analysis of the pooled data from Sweden, Denmark and Austria, 6.0% [5.1–6.8] of the patients developed anti-natalizumab ADA. As expected most of the ADA occurred within six months of natalizumab therapy. The country was associated with occurrence of ADA (p = 0.004) (Fig 3).

thumbnail
Download:
Fig 3. Cumulative probability of developing anti-natalizumab ADA.

Estimated cumulative probability of developing anti-natalizumab ADA by categorical variables with p-values for rank-based tests adapted for interval censoring.

https://doi.org/10.1371/journal.pone.0162752.g003

In Sweden 4.6% (95% CI 3.7–5.9) of the patients developed ADA before six months of therapy, while 7.2% (95% CI 5.0–8.9) and 7.7% (95% CI 5.3–9.8) did in Denmark and Austria respectively (Table 2).

thumbnail
Download:
Table 2. Estimated cumulative probabilities of ADA development before six months of natalizumab therapy.

https://doi.org/10.1371/journal.pone.0162752.t002

Age at start of natalizumab therapy was associated with ADA occurrence (p = 0.03). Among patients older than 45 years, 7.5% (95% CI 4.9–9.5) had a detectable immunogenicity at six months of therapy compared to 5.5% (95% CI 4.6–6.5) of the younger patients.

Sex was significantly associated with ADA occurrence (p = 0.03). More female patients (6.6%, 95% CI 5.6–7.8) have measurable ADA than males (4.5%, 95% CI 3.2–5.9) at six months of therapy.

None of the months of natalizumab therapy start and none of the months of birth were associated with ADA production when the countries were analyzed separately. However in the pooled analysis, patients starting the therapy in November had a higher risk of developing ADA (S3 and S4 Tables).

None of the variables characterizing previous IFNβ therapy (last IFNβ taken, first IFNβ taken, last ADA status during IFNβ therapy, total duration of IFNβ therapy and number of different IFNβ therapies) was significantly associated with development of ADA against natalizumab (S1 Fig).

Multivariate analyses

Results for IFNβ.

Type of IFNβ, age at start of therapy, sex and month of therapy start were included in a multivariate Cox regression model for each cohort and results are reported in Table 3 and Fig 4.

thumbnail
Download:
Fig 4. Multivariate risk factors for anti-IFNβ ADA development.

Multivariate risk factors for anti-IFNβ ADA development in each cohort, adjusted on assay method at the last visit in Sweden and Austria.

https://doi.org/10.1371/journal.pone.0162752.g004

thumbnail
Download:
Table 3. Adjusted hazard ratios for risk factors of anti-IFNβ ADA development.

https://doi.org/10.1371/journal.pone.0162752.t003

As the ADA assay method was acting as a confounding factor for the Swedish and Austrian cohorts, we adjusted for this variable in the Cox model. The adjusted hazard ratios (HR) varied from 4.9 to 7.8 for IFNβ-1a s.c. compared to IFNβ-1a i.m. (95% CI ranging from 2.9 to 20.0) and from 6.4 to 10.4 for IFNβ-1b s.c. compared to IFNβ-1a i.m. (95% CI ranging from 3.2 to 25.9) (Table 3 and Fig 4).

HRs relative to age older than 50 compared to age between 18 and 30 varied from 1.3 to 1.9 across the cohorts (95% CI ranging from 0.8 to 3.8) (Table 3 and Fig 4).

The adjusted HRs for males compared to females were significant in Austria (1.4, 95% CI 1.1–1.9) and Denmark (1.6, 95% CI 1.0–2.7), at the limit of significance in Germany-Dusseldorf (1.4, 95% CI 0.9–2.0) and not significant in Sweden (1.0, 95% CI 0.7–1.3) (Table 3 and Fig 4).

After adjustment for other risk factors, start of therapy in April compared to start in any other month was still a risk factor for ADA development in Sweden (1.6, 95% CI 1.1–2.4) and Germany (Dusseldorf) (2.4, 95% CI 1.5–3.9), and was associated with a lower risk in Austria (0.6, 95% CI 0.3–1.0) (Table 3 and Fig 4).

According to the analysis of the pooled dataset of the five cohorts we found similar associations with type of treatment (HR = 6.4 [4.9–8.4] for IFNβ-1a s.c. and 8.7 [6.6–11.4] for IFNβ-1b s.c. both versus IFNβ-1a i.m.), age (HR = 1.8 [1.4–2.3] for patients older than 50 versus 18–30), sex (HR = 1.3 [1.1–1.6] for males versus females), month of therapy start (HR = 1.3 [1.0–1.6] for April compared to other months) (Table 3).

Results for natalizumab.

Age at start of natalizumab therapy and sex were analyzed together with a Cox regression model for each cohort and for the pooled dataset with adjustment on the country (Table 4 and Fig 5). The two variables were identified as risk factors for ADA occurrence in the pooled cohort, whereas the associations with ADA development were not statistically significant in each individual cohort.

thumbnail
Download:
Fig 5. Multivariate risk factors for anti-natalizumab ADA development.

(N = 3440).

https://doi.org/10.1371/journal.pone.0162752.g005

thumbnail
Download:
Table 4. Adjusted hazard ratios for risk factors of anti-natalizumab ADA development.

https://doi.org/10.1371/journal.pone.0162752.t004

HRs for age older than 45 at start of therapy compared to age 18–45 varied from 1.2 to 1.6 across the cohorts (95% CI ranging from 0.7 to 3.0) and was 1.4 (95% CI 1.0–1.8) over all the cohorts pooled together.

The HR for men compared to women was 0.7 (95% CI 0.4–1.3) in Sweden, 0.8 (95% CI 0.5–1.3) in Denmark and 0.6 (95% CI 0.3–1.2) in Austria. In the pooled dataset, it was 0.7 (95% CI 0.5–0.98) (Table 4 and Fig 5).

Discussion

To the best of our knowledge, this is the first collaborative cohort analysis investigating the impact of patient characteristics and type of BPs on the occurrence of ADA in MS. In order to carry out an unbiased evaluation of the relationships between bio-clinical factors and the occurrence of ADAs, we have selected a target population which was extracted from the whole ABIRISK database. This selected population is composed by naive adult patients with biological samples being taken within the first two years of treatment. From this selected population which combines data from four national European reference laboratories and one research center, we report that the occurrence of ADA against BPs used to treat MS (IFNβ and natalizumab) were associated with age and sex. Older patients were at higher risk of developing ADA against both IFNβ and natalizumab. In contrast, the sex effect differs between the two BPs, with male patients being more at risk of anti-IFNβ ADA and female patients more at risk of anti-natalizumab ADA.

We confirm the previously reported differences in the cumulative probability of ADA occurrence at 18 months of treatment for the different types of IFNβ (ranging between 1.1–4.2% for IFNβ-1a i.m., 11.1–23.0% for IFNβ-1a s.c. and 17.4–33.3% for IFNβ-1b s.c. in the present study) [1].

Interestingly, we also observe that patients who started IFNβ therapy in the month of April showed a higher risk of ADA development. This relationship is not seen in every cohort but stays significant in the German (Dusseldorf) and in the Swedish cohorts after adjustment for multiple testing. We do not find a relationship between ADA occurrence and the month of birth neither for IFNβ nor for natalizumab-treated patients. In the latter, no factor related to previous IFNβ exposure (last IFNβ taken, first IFNβ taken, last ADA status during IFNβ therapy, total duration of IFNβ therapy and number of different IFNβ therapies) is associated to anti-natalizumab ADA development.

Since this study relies on historical cohorts, we had to cope with various follow-up schemes. For example, in Sweden and Denmark the assays were done routinely, while in Germany the testing was left to the discretion of the treating physician. However, from a previous study, the characteristics of the patients who have been tested due to an indication of worsening disease and therefore a suspicion of ADA did not differ significantly in ADA frequency from those that were tested according to regular routine dates [27]. Moreover, in our analyses we observed similar follow-up schemes by age and sex.

In this study, we do not know the exact time of ADA occurrence but we know whether or not ADA has occurred within a monitoring interval. These so-called interval censored data have often been analyzed with various ad-hoc solutions such as the use of either the right-sided endpoint for subject with positive ADA or the midpoint in the interval. These solutions are known to be non-optimal and biased when the times between visits are long or the frequency of ADA assessment varies in the population which is the case in our study that included various cohorts with different follow-up scheme. Here, all the analyses were performed using methodologies designed for such interval-censored data.

In this study, it is worth noting that different types of assay were sequentially used over time and that these differences might have reflected changes in treatment choice over time or changes in treatment formulation as for the IFNβ-1a s.c. in 2007 [17]. We cannot assess the effect of change in IFNβ-1a s.c. formulation as its date is too close to the date of change in type of assay used. Nevertheless when we adjusted on the type of IFNβ treatment in the multivariate analyses, differences between the assay methods in terms of ADA development remained significant, suggesting that assay methods may differ in sensitivity of ADA detection as previously reported [18].

Due to the wide range of ADA detection process for IFNβ, the analyses were initially performed separately on each country. Access to reference laboratory data enabled us to have sufficient number of patients in each country that gave us sufficient power to detect the association of interest if it existed. Since we did not observe discrepancies between the results obtained for each country, we then also performed the multivariate analysis on the pooled cohort.

For natalizumab, since the same assay method was used and the laboratories had been tested yearly for concordant test results of 20 samples (personal communication from A. Fogdell-Hahn), the univariate analysis was applied on the pooled cohort. As we observed a country effect, multivariate analysis was performed both on separate and pooled cohort.

Our selection criteria implied that we included more patients who began their treatment more recently. This led for instance to the fact that our study population includes more IFNβ-1a i.m. treated patients than the global ABIRISK database. Because of the differences between included and excluded patients in terms of type of treatment, age at first treatment or first ADA status, our study population might not be representative of the general population of the IFNβ-treated or natalizumab-treated and ADA tested MS patients. This necessary selection in order to cope with the two-year ADA occurrence window could influence the global prevalence of ADA development, but not the tests of comparisons between groups and the measure of associations between potential risk factors and ADA development. Another consequence of the selection is that we might lose power in some less represented subcategories. For instance, older patients in Germany (Dusseldorf) were less included because it might have been more difficult to track their date of first dose, which may explain why the association with age could not be seen.

The differences among treatments that we observed were in the same range and in the same direction as those described in previous studies [1]. The higher immunogenicity of IFNβ-1b compared to IFNβ-1a has been attributed to its lack of glycosylation which makes it highly prone to form aggregates. This has been confirmed in an animal model, where aggregated IFNβ-1b and IFNβ-1a induced higher ADA seroprevalence than non-aggregated IFNβ-1a [28]. Animal studies have also suggested that both a T-independent and a T-dependent response might contribute to IFNβ-1b immunogenicity, where the T-independent response is characterized by low ADA titers and low ADA persistency due to lack of a memory response [29]. Conversely, IFNβ-1a immunogenicity, although less frequent, generates higher titer and persistent ADA responses, often with a neutralizing activity [15,21]. In our study as in previous ones, although risk of ADA occurrence was higher in IFNβ-1b s.c. treated MS-patients than in IFNβ-1a i.m. treated, the difference with IFNβ-1a s.c.-treated was not significant.

We found a significantly higher immunogenicity of both IFNβ and natalizumab in older patients. The impact of age was assessed in previous studies but was not found to be associated with ADA development [30,31]. However a higher median age was observed in ADA positive compared to ADA negative patients in one study (45 years versus 43 years, p<0.0001, n = 1115) [15]. The negative studies had lower number of patients (n = 301 and n = 541) and may have lacked in power to detect the association.

As thymic function declines after puberty, the main source of naïve T cells to counterbalance T-cell turnover in adult life comes from their homeostatic proliferation [32]. Mathematical models corroborated by experimental data postulate that in spite of T-cell turnover, homeostatic mechanisms keep a constant T-cell pool and the maintenance of T-cell repertoire diversity is IL7- and TCR-dependent [33–35]. Continuous loss of senescent T-cells and conversion of naïve T cells into memory T cells due to pathogen exposure progressively induce a reduction of the naïve T-cell repertoire diversity and lymphopenia with a consequent decrease in adaptive immune responses against novel antigens and increase in homeostatic proliferation in the elderly [36]. The lymphopenic condition is known to induce a stronger proliferation of high affinity autoreactive effector T cells to self-antigens, and if not compensated by an appropriate proliferation or novel thymic production of regulatory T cells it may easily lead to development of autoimmune responses [37]. Interestingly, more than 50% of alemtuzumab-treated MS-patients experience secondary autoimmune manifestations as a consequence of lymphopenia, and development of autoimmunity is strongly associated with insufficient recovery of thymopoiesis after a strong initial homeostatic proliferation [38]. Similarly, lymphopenia induced autoimmunity might explain the increased risk of anti-IFNβ ADA development in older MS-patients: chronic stimuli with IFNβ might induce homeostatic proliferation of effector T cells against a self-antigen not paralleled by an appropriate thymopoiesis of novel regulatory T cells.

In this study we identify a significant association of sex with ADA occurrence, with men being more susceptible to develop IFNβ ADA. In one study about immunogenicity of IFNβ in MS, it was observed that a higher proportion of female than male patients became ADA positive after 12 months of treatment with IFNβ-1a subcutaneously [39]. Later this result was not confirmed as sex was not a predictor of development of neutralizing antibodies in clinical data collected prospectively in Denmark [31]. Apart from these studies, the association between sex-related factors and ADA development has not been investigated so far. Their characterization is of utmost importance since some of these factors may be controlled (e.g. oral contraception, breastfeeding) for prevention of ADA development. A thorough knowledge of X-linked genetic factors, hormonal factors and behavioral factors is needed to distinguish which sex- (or gender-) specific factor mostly contributes to this difference. Among interesting gene candidates encoded on the X-chromosome are several transcription factors such as FOXP3 and immunity receptors such as TLR7 and TLR8. Sex hormones are important modulators of immune responses and several autoimmune diseases display a sex bias with a higher prevalence in women [40]. The hormonal status of the patient at therapy start is therefore likely to influence the development of ADA. Interestingly the sex difference in IFNβ ADA occurrence was observed in all the analyzed European countries except Sweden. A hypothesis to account for this peculiarity might be a local behavioral difference influencing a relevant exposure such that the exposure would be similar in men and women in Sweden, but higher in men in the other countries.

Concerning natalizumab, inversely to IFNβ, we observed that women developed ADA more often than men. IFNβ being a cytokine and natalizumab a monoclonal antibody, they are structurally different and might have distinct immunogenicity mechanisms. A slower pharmacokinetics in women than in men has been reported for several drugs, including the TNFα-blocking monoclonal antibody adalimumab [41], and might induce a prolonged exposure to the drug in women with possible effects on immunogenicity. The early pharmacokinetics study on natalizumab did not analyze the influence of sex [42], and a sex difference in natalizumab clearance warrants further investigation.

An association of month of birth with risk of developing MS has clearly been established in several case-control studies from the Northern and Southern hemisphere, albeit with a small effect size [43]. In the Northern hemisphere individuals born in April-May carry a higher risk of MS, whereas individuals born in October-November display a slight protection from MS. The month of birth effect is thought to derive from fetal exposure to a seasonal factor in the last trimester of pregnancy, a crucial period in the negative selection of autoreactive T cells in the fetal thymus [44]. As ADA (in particular against IFNβ) may be considered an autoimmune phenomenon, we investigated whether a month of birth effect could be linked also to ADA occurrence, but we were not able to detect an association. Considering the small effect size which has been reported for the association between month of birth and MS (OR = 1.08 for May in a UK case-control study [45]), we did not have the power to detect a similarly small effect in our study.

Instead we observed a significant association of the month of therapy start with IFNβ ADA in Sweden, Germany (Dusseldorf group), Austria, and borderline significant in Germany (Munich group). For all the cohorts the month associated with ADA was April, but surprisingly start in April conferred an increased risk in all the cohorts except Austria, where it conferred decreased risk. Further investigations are warranted to confirm this association and to identify the underlying seasonal factor. Candidate seasonal factors which might influence the immune status at start of therapy are allergy causing spring pollens or early spring vitamin D insufficiency.

Vitamin D is a modulator of T- and B-cell responses [46], and genetically and environmentally determined vitamin D insufficiency is a risk factor for MS [47,48]. B-lymphocytes express the VDR and vitamin D is able to suppress their proliferation, immunoglobulin production, class switch and differentiation to plasma cells [49], thus it might play a role in ADA development. It has been suggested that part of the effect of IFNβ is mediated by induction of vitamin D synthesis [50]. From these results and from the previous knowledge of the immunomodulatory activity of vitamin D it could be hypothesized that vitamin D levels at start of therapy, influenced by the seasonal UV-radiation and possibly also by early response to IFNβ therapy, is a factor associated with ADA development. Several randomized controlled clinical trials on vitamin D supplementation to IFNβ or glatiramer therapy in MS are ongoing [51,52]. Studies to correlate vitamin D levels before and during IFNβ therapy with later ADA development are warranted.

In conclusion, the findings obtained from this consortium cohort provide new insights on ADA occurrence. Age and sex were found to be risk factors of ADA development against IFNβ and natalizumab MS therapies in data from four national European reference laboratories and one research study. We also observed an association with the month of start of therapy which warrants further investigation. The next step of our research project, which is currently ongoing, is to include data from additional MS registries and research study questionnaires to assess other potential explanatory factors such as vitamin D serum levels, previous and concomitant medications, medical history, infections, vaccinations, tobacco consumption and genetic factors. The effect size of the associations with age and sex were not high enough to discourage clinicians in terms of drug prescription, but a more frequent ADA monitoring is recommended for the risk groups to allow earlier detection of ADA and an appropriate drug switch.

Supporting Information

S1 Fig. Association of variables characterizing previous IFNβ therapy with anti-natalizumab ADA occurrence risk.

https://doi.org/10.1371/journal.pone.0162752.s001

(TIF)

S1 File. List of ABIRISK Consortium partners and members.

https://doi.org/10.1371/journal.pone.0162752.s002

(DOCX)

S1 Table. Anti-IFNβ ADA development hazard ratios with 95% confidence intervals for each month of start of therapy (versus all the other months).

https://doi.org/10.1371/journal.pone.0162752.s003

(DOCX)

S2 Table. Anti-IFNβ ADA development hazard ratios with 95% confidence intervals for each month of birth (versus all the other months).

https://doi.org/10.1371/journal.pone.0162752.s004

(DOCX)

S3 Table. Anti-natalizumab ADA development hazard ratios with 95% confidence intervals for each month of start of therapy (versus all the other months).

https://doi.org/10.1371/journal.pone.0162752.s005

(DOCX)

S4 Table. Anti-natalizumab ADA development hazard ratios with 95% confidence intervals for each month of birth (versus all the other months).

https://doi.org/10.1371/journal.pone.0162752.s006

(DOCX)

S5 Table. Comparison of the of the included to the excluded IFNβ-treated patients.

https://doi.org/10.1371/journal.pone.0162752.s007

(DOCX)

Acknowledgments

The research leading to these results was conducted as part of the ABIRISK consortium (Anti-Biopharmaceutical).

Immunization: Prediction and analysis of clinical relevance to minimize the risk). For further information, please refer to www.abirisk.eu. All the authors of this article are members of the ABIRISK consortium.

We thank the eTRIKS consortium (https://www.etriks.org) for their support in the construction of the tranSMART database.

Author Contributions

  1. Conceptualization: D. Bachelet SH JL RR JD AF-H PB.
  2. Data curation: JL MA PEH-J NK-H CW KI D. Buck VG AL ND PD.
  3. Formal analysis: D. Bachelet SH.
  4. Methodology: CM PB.
  5. Project administration: AH-M DS MP JD AF-H PB.
  6. Resources: MA PEH-J NK-H CW D. Buck BK BH HPH PSS FD AF-H.
  7. Supervision: PB.
  8. Visualization: D. Bachelet SH.
  9. Writing – original draft: D. Bachelet SH PB.
  10. Writing – review & editing: CM JL MR RR PEH-J CW D. Buck BK BH HPH FD JD AF-H PD AH-M.

References

  1. 1. Deisenhammer F. Neutralizing antibodies to interferon-beta and other immunological treatments for multiple sclerosis: prevalence and impact on outcomes. CNS Drugs. 2009;23(5):379–96. pmid:19453200
  2. 2. Bertolotto A, Capobianco M, Amato MP, Capello E, Capra R, Centonze D, et al. Guidelines on the clinical use for the detection of neutralizing antibodies (NAbs) to IFN beta in multiple sclerosis therapy: report from the Italian Multiple Sclerosis Study group. Neurol Sci Off J Ital Neurol Soc Ital Soc Clin Neurophysiol. 2014 Feb;35(2):307–16.
  3. 3. Sørensen PS, Deisenhammer F, Duda P, Hohlfeld R, Myhr K-M, Palace J, et al. Guidelines on use of anti-IFN-beta antibody measurements in multiple sclerosis: report of an EFNS Task Force on IFN-beta antibodies in multiple sclerosis. Eur J Neurol. 2005 Nov;12(11):817–27. pmid:16241970
  4. 4. Deisenhammer F. Interferon-Beta: Neutralizing Antibodies, Binding Antibodies, Pharmacokinetics and Pharmacodynamics, and Clinical Outcomes. J Interferon Cytokine Res Off J Int Soc Interferon Cytokine Res. 2014 Dec;34(12):938–45.
  5. 5. Hesse D, Sørensen PS. Using measurements of neutralizing antibodies: the challenge of IFN-β therapy. Eur J Neurol. 2007;14(8):850–859. pmid:17662004
  6. 6. Sorensen PS, Koch-Henriksen N, Ross C, Clemmesen KM, Bendtzen K. Appearance and disappearance of neutralizing antibodies during interferon-beta therapy. Neurology. 2005 Jul 12;65(1):33–9. pmid:15888603
  7. 7. Calabresi PA, Giovannoni G, Confavreux C, Galetta SL, Havrdova E, Hutchinson M, et al. The incidence and significance of anti-natalizumab antibodies: results from AFFIRM and SENTINEL. Neurology. 2007 Oct 2;69(14):1391–403. pmid:17761550
  8. 8. Oliver-Martos B, Órpez-Zafra T, Urbaneja P, Maldonado-Sanchez R, Leyva L, Fernández O. Early development of anti-natalizumab antibodies in MS patients. J Neurol [Internet]. 2013 Jun 14 [cited 2013 Jul 31]; Available: http://link.springer.com.gate2.inist.fr/content/pdf/10.1007%2Fs00415-013-6991-2.pdf
  9. 9. Lundkvist M, Greiner E, Hillert J, Fogdell-Hahn A. Multiple sclerosis patients lacking oligoclonal bands in the cerebrospinal fluid are less likely to develop neutralizing antibodies against interferon beta. Mult Scler. 2010 Jan 7;16(7):796–800. pmid:20534645
  10. 10. Buck D, Cepok S, Hoffmann S, Grummel V, Jochim A, Berthele A, et al. Influence of the HLA-DRB1 genotype on antibody development to interferon beta in multiple sclerosis. Arch Neurol. 2011 Apr;68(4):480–7. pmid:21482927
  11. 11. Hoffmann S, Cepok S, Grummel V, Lehmann-Horn K, Hackermüller J, Hackermueller J, et al. HLA-DRB1*0401 and HLA-DRB1*0408 are strongly associated with the development of antibodies against interferon-beta therapy in multiple sclerosis. Am J Hum Genet. 2008 Aug;83(2):219–27. pmid:18656179
  12. 12. Link J, Lundkvist Ryner M, Fink K, Hermanrud C, Lima I, Brynedal B, et al. Human leukocyte antigen genes and interferon beta preparations influence risk of developing neutralizing anti-drug antibodies in multiple sclerosis. PloS One. 2014;9(3):e90479. pmid:24608124
  13. 13. Núñez C, Cénit MC, Alvarez-Lafuente R, Río J, Fernández-Arquero M, Arroyo R, et al. HLA alleles as biomarkers of high-titre neutralising antibodies to interferon-β therapy in multiple sclerosis. J Med Genet. 2014 Jun;51(6):395–400. pmid:24748646
  14. 14. Szalma S, Koka V, Khasanova T, Perakslis ED. Effective knowledge management in translational medicine. J Transl Med. 2010;8:68. pmid:20642836
  15. 15. Sominanda A, Rot U, Suoniemi M, Deisenhammer F, Hillert J, Fogdell-Hahn A. Interferon beta preparations for the treatment of multiple sclerosis patients differ in neutralizing antibody seroprevalence and immunogenicity. Mult Scler Houndmills Basingstoke Engl. 2007 Mar;13(2):208–14.
  16. 16. Bertolotto A, Sala A, Caldano M, Capobianco M, Malucchi S, Marnetto F, et al. Development and validation of a real time PCR-based bioassay for quantification of neutralizing antibodies against human interferon-beta. J Immunol Methods. 2007 Apr 10;321(1–2):19–31. pmid:17335844
  17. 17. Jungedal R, Lundkvist M, Engdahl E, Ramanujam R, Westerlind H, Sominanda A, et al. Prevalence of anti-drug antibodies against interferon beta has decreased since routine analysis of neutralizing antibodies became clinical practice. Mult Scler Houndmills Basingstoke Engl. 2012 Dec;18(12):1775–81.
  18. 18. Hermanrud C, Ryner ML, Engdahl E, Fogdell-Hahn A. Anti-interferon beta antibody titers strongly correlate between two bioassays and in vivo biomarker expression, and indicates that a titer of 150 TRU/mL is a biologically functional cut-point. J Interferon Cytokine Res Off J Int Soc Interferon Cytokine Res. 2014 Jul;34(7):498–504.
  19. 19. Kawade Y. Quantitation of neutralization of interferon by antibody. Methods Enzymol. 1986;119:558–73. pmid:2429157
  20. 20. Lam R, Farrell R, Aziz T, Gibbs E, Giovannoni G, Grossberg S, et al. Validating parameters of a luciferase reporter gene assay to measure neutralizing antibodies to IFNbeta in multiple sclerosis patients. J Immunol Methods. 2008 Jul 31;336(2):113–8. pmid:18511063
  21. 21. Hegen H, Millonig A, Bertolotto A, Comabella M, Giovanonni G, Guger M, et al. Early detection of neutralizing antibodies to interferon-beta in multiple sclerosis patients: binding antibodies predict neutralizing antibody development. Mult Scler Houndmills Basingstoke Engl. 2014 Apr;20(5):577–87.
  22. 22. Pachner AR. An improved ELISA for screening for neutralizing anti-IFN-beta antibodies in MS patients. Neurology. 2003 Nov 25;61(10):1444–6. pmid:14638976
  23. 23. Turnbull BW. The Empirical Distribution Function with Arbitrarily Grouped, Censored and Truncated Data. J R Stat Soc Ser B Methodol. 1976 Jan 1;38(3):290–5.
  24. 24. Fay MP, Shaw PA. Exact and Asymptotic Weighted Logrank Tests for Interval Censored Data: The interval R package. J Stat Softw. 2010 Aug;36(2). pmid:25285054
  25. 25. Gentleman R, Geyer CJ. Maximum likelihood for interval censored data: Consistency and computation. Biometrika. 1994 Jan 9;81(3):618–23.
  26. 26. Pan W. A multiple imputation approach to Cox regression with interval-censored data. Biometrics. 2000 Mar;56(1):199–203. pmid:10783796
  27. 27. Rot U, Sominanda A, Fogdell-Hahn A, Hillert J. Impression of clinical worsening fails to predict interferon-beta neutralizing antibody status. J Int Med Res. 2008 Dec;36(6):1418–25. pmid:19094453
  28. 28. van Beers MMC, Jiskoot W, Schellekens H. On the role of aggregates in the immunogenicity of recombinant human interferon beta in patients with multiple sclerosis. J Interferon Cytokine Res Off J Int Soc Interferon Cytokine Res. 2010 Oct;30(10):767–75.
  29. 29. Sauerborn M, van Beers MMC, Jiskoot W, Kijanka GM, Boon L, Schellekens H, et al. Antibody response against Betaferon® in immune tolerant mice: involvement of marginal zone B-cells and CD4+ T-cells and apparent lack of immunological memory. J Clin Immunol. 2013 Jan;33(1):255–63. pmid:22945588
  30. 30. Rudick RA, Simonian NA, Alam JA, Campion M, Scaramucci JO, Jones W, et al. Incidence and significance of neutralizing antibodies to interferon beta-1a in multiple sclerosis. Multiple Sclerosis Collaborative Research Group (MSCRG). Neurology. 1998 May;50(5):1266–72. pmid:9595973
  31. 31. Sorensen PS, Ross C, Clemmesen KM, Bendtzen K, Frederiksen JL, Jensen K, et al. Clinical importance of neutralising antibodies against interferon beta in patients with relapsing-remitting multiple sclerosis. Lancet. 2003 Oct 11;362(9391):1184–91. pmid:14568740
  32. 32. Kilpatrick RD, Rickabaugh T, Hultin LE, Hultin P, Hausner MA, Detels R, et al. Homeostasis of the naive CD4+ T cell compartment during aging. J Immunol Baltim Md 1950. 2008 Feb 1;180(3):1499–507.
  33. 33. Sportès C, Hakim FT, Memon SA, Zhang H, Chua KS, Brown MR, et al. Administration of rhIL-7 in humans increases in vivo TCR repertoire diversity by preferential expansion of naive T cell subsets. J Exp Med. 2008 Jul 7;205(7):1701–14. pmid:18573906
  34. 34. Mackall CL, Bare CV, Granger LA, Sharrow SO, Titus JA, Gress RE. Thymic-independent T cell regeneration occurs via antigen-driven expansion of peripheral T cells resulting in a repertoire that is limited in diversity and prone to skewing. J Immunol Baltim Md 1950. 1996 Jun 15;156(12):4609–16.
  35. 35. Johnson PLF, Goronzy JJ, Antia R. A population biological approach to understanding the maintenance and loss of the T-cell repertoire during aging. Immunology. 2014 Jun;142(2):167–75. pmid:24405293
  36. 36. Sauce D, Larsen M, Fastenackels S, Roux A, Gorochov G, Katlama C, et al. Lymphopenia-driven homeostatic regulation of naive T cells in elderly and thymectomized young adults. J Immunol Baltim Md 1950. 2012 Dec 15;189(12):5541–8.
  37. 37. Krupica T, Fry TJ, Mackall CL. Autoimmunity during lymphopenia: a two-hit model. Clin Immunol Orlando Fla. 2006 Aug;120(2):121–8.
  38. 38. Jones JL, Thompson SAJ, Loh P, Davies JL, Tuohy OC, Curry AJ, et al. Human autoimmunity after lymphocyte depletion is caused by homeostatic T-cell proliferation. Proc Natl Acad Sci U S A. 2013 Dec 10;110(50):20200–5. pmid:24282306
  39. 39. Ross C, Clemmesen KM, Svenson M, Sørensen PS, Koch-Henriksen N, Skovgaard GL, et al. Immunogenicity of interferon-beta in multiple sclerosis patients: influence of preparation, dosage, dose frequency, and route of administration. Danish Multiple Sclerosis Study Group. Ann Neurol. 2000 Nov;48(5):706–12. pmid:11079533
  40. 40. González DA, Díaz BB, Rodríguez Pérez M del C, Hernández AG, Chico BND, de León AC. Sex hormones and autoimmunity. Immunol Lett. 2010 Sep 6;133(1):6–13. pmid:20637236
  41. 41. Ternant D, Bejan-Angoulvant T, Passot C, Mulleman D, Paintaud G. Clinical Pharmacokinetics and Pharmacodynamics of Monoclonal Antibodies Approved to Treat Rheumatoid Arthritis. Clin Pharmacokinet. 2015 Nov;54(11):1107–23. pmid:26123705
  42. 42. Sheremata WA, Vollmer TL, Stone LA, Willmer-Hulme AJ, Koller M. A safety and pharmacokinetic study of intravenous natalizumab in patients with MS. Neurology. 1999 Mar 23;52(5):1072–4. pmid:10102433
  43. 43. Dobson R, Giovannoni G, Ramagopalan S. The month of birth effect in multiple sclerosis: systematic review, meta-analysis and effect of latitude. J Neurol Neurosurg Psychiatry. 2013 Apr;84(4):427–32. pmid:23152637
  44. 44. Handunnetthi L, Ramagopalan SV, Ebers GC. Multiple sclerosis, vitamin D, and HLA-DRB1*15. Neurology. 2010 Jun 8;74(23):1905–10. pmid:20530326
  45. 45. Disanto G, Chaplin G, Morahan JM, Giovannoni G, Hyppönen E, Ebers GC, et al. Month of birth, vitamin D and risk of immune-mediated disease: a case control study. BMC Med. 2012;10:69. pmid:22764877
  46. 46. Baeke F, Takiishi T, Korf H, Gysemans C, Mathieu C. Vitamin D: modulator of the immune system. Curr Opin Pharmacol. 2010 Aug;10(4):482–96. pmid:20427238
  47. 47. Mokry LE, Ross S, Ahmad OS, Forgetta V, Smith GD, Leong A, et al. Vitamin D and Risk of Multiple Sclerosis: A Mendelian Randomization Study. PLoS Med. 2015 Aug;12(8):e1001866. pmid:26305103
  48. 48. Ascherio A, Munger KL, Simon KC. Vitamin D and multiple sclerosis. Lancet Neurol. 2010 Jun;9(6):599–612. pmid:20494325
  49. 49. Chen S, Sims GP, Chen XX, Gu YY, Chen S, Lipsky PE. Modulatory effects of 1,25-dihydroxyvitamin D3 on human B cell differentiation. J Immunol Baltim Md 1950. 2007 Aug 1;179(3):1634–47.
  50. 50. Stewart N, Simpson S Jr, van der Mei I, Ponsonby A-L, Blizzard L, Dwyer T, et al. Interferon-β and serum 25-hydroxyvitamin D interact to modulate relapse risk in MS. Neurology. 2012 Jul 17;79(3):254–60. pmid:22700816
  51. 51. Dörr J, Ohlraun S, Skarabis H, Paul F. Efficacy of vitamin D supplementation in multiple sclerosis (EVIDIMS Trial): study protocol for a randomized controlled trial. Trials. 2012;13:15. pmid:22316314
  52. 52. Smolders J, Hupperts R, Barkhof F, Grimaldi LME, Holmoy T, Killestein J, et al. Efficacy of vitamin D3 as add-on therapy in patients with relapsing-remitting multiple sclerosis receiving subcutaneous interferon β-1a: a Phase II, multicenter, double-blind, randomized, placebo-controlled trial. J Neurol Sci. 2011 Dec 15;311(1–2):44–9. pmid:21620416