Clinicogenomic factors of biotherapy immunogenicity in autoimmune disease: A prospective multicohort study of the ABIRISK consortium.

BACKGROUND: Biopharmaceutical products (BPs) are widely used to treat autoimmune diseases, but immunogenicity limits their efficacy for an important proportion of patients. Our knowledge of patient-related factors influencing the occurrence of antidrug antibodies (ADAs) is still limited. METHODS AND FINDINGS: The European consortium ABIRISK (Anti-Biopharmaceutical Immunization: prediction and analysis of clinical relevance to minimize the RISK) conducted a clinical and genomic multicohort prospective study of 560 patients with multiple sclerosis (MS, n = 147), rheumatoid arthritis (RA, n = 229), Crohn's disease (n = 148), or ulcerative colitis (n = 36) treated with 8 different biopharmaceuticals (etanercept, n = 84; infliximab, n = 101; adalimumab, n = 153; interferon [IFN]-beta-1a intramuscularly [IM], n = 38; IFN-beta-1a subcutaneously [SC], n = 68; IFN-beta-1b SC, n = 41; rituximab, n = 31; tocilizumab, n = 44) and followed during the first 12 months of therapy for time to ADA development. From the bioclinical data collected, we explored the relationships between patient-related factors and the occurrence of ADAs. Both baseline and time-dependent factors such as concomitant medications were analyzed using Cox proportional hazard regression models. Mean age and disease duration were 35.1 and 0.85 years, respectively, for MS; 54.2 and 3.17 years for RA; and 36.9 and 3.69 years for inflammatory bowel diseases (IBDs). In a multivariate Cox regression model including each of the clinical and genetic factors mentioned hereafter, among the clinical factors, immunosuppressants (adjusted hazard ratio [aHR] = 0.408 [95% confidence interval (CI) 0.253-0.657], p < 0.001) and antibiotics (aHR = 0.121 [0.0437-0.333], p < 0.0001) were independently negatively associated with time to ADA development, whereas infections during the study (aHR = 2.757 [1.616-4.704], p < 0.001) and tobacco smoking (aHR = 2.150 [1.319-3.503], p < 0.01) were positively associated. 351,824 Single-Nucleotide Polymorphisms (SNPs) and 38 imputed Human Leukocyte Antigen (HLA) alleles were analyzed through a genome-wide association study. We found that the HLA-DQA1*05 allele significantly increased the rate of immunogenicity (aHR = 3.9 [1.923-5.976], p < 0.0001 for the homozygotes). Among the 6 genetic variants selected at a 20% false discovery rate (FDR) threshold, the minor allele of rs10508884, which is situated in an intron of the CXCL12 gene, increased the rate of immunogenicity (aHR = 3.804 [2.139-6.764], p < 1 × 10-5 for patients homozygous for the minor allele) and was chosen for validation through a CXCL12 protein enzyme-linked immunosorbent assay (ELISA) on patient serum at baseline before therapy start. CXCL12 protein levels were higher for patients homozygous for the minor allele carrying higher ADA risk (mean: 2,693 pg/ml) than for the other genotypes (mean: 2,317 pg/ml; p = 0.014), and patients with CXCL12 levels above the median in serum were more prone to develop ADAs (aHR = 2.329 [1.106-4.90], p = 0.026). A limitation of the study is the lack of replication; therefore, other studies are required to confirm our findings. CONCLUSION: In our study, we found that immunosuppressants and antibiotics were associated with decreased risk of ADA development, whereas tobacco smoking and infections during the study were associated with increased risk. We found that the HLA-DQA1*05 allele was associated with an increased rate of immunogenicity. Moreover, our results suggest a relationship between CXCL12 production and ADA development independent of the disease, which is consistent with its known function in affinity maturation of antibodies and plasma cell survival. Our findings may help physicians in the management of patients receiving biotherapies.

Assessing the effect of genetic markers on drug immunogenicity from a mechanistic model-based approach.

BACKGROUND: With the growth in use of biotherapic drugs in various medical fields, the occurrence of anti-drug antibodies represents nowadays a serious issue. This immune response against a drug can be due either to pre-existing antibodies or to the novel production of antibodies from B-cell clones by a fraction of the exposed subjects. Identifying genetic markers associated with the immunogenicity of biotherapeutic drugs may provide new opportunities for risk stratification before the introduction of the drug. However, real-world investigations should take into account that the population under study is a mixture of pre-immune, immune-reactive and immune-tolerant subjects. METHOD: In this work, we propose a novel test for assessing the effect of genetic markers on drug immunogenicity taking into account that the population under study is a mixed one. This test statistic is derived from a novel two-part semiparametric improper survival model which relies on immunological mechanistic considerations. RESULTS: Simulation results show the good behavior of the proposed statistic as compared to a two-part logrank test. In a study on drug immunogenicity, our results highlighted findings that would have been discarded when considering classical tests. CONCLUSION: We propose a novel test that can be used for analyzing drug immunogenicity and is easy to implement with standard softwares. This test is also applicable for situations where one wants to test the equality of improper survival distributions of semi-continuous outcomes between two or more independent groups.

Quantifying the Predictive Accuracy of a Polygenic Risk Score for Predicting Incident Cancer Cases : Application to the CARTaGENE Cohort.

With the increasing use of polygenic risk scores (PRS) there is a need for adapted methods to evaluate the predictivity of these tools. In this work, we propose a new pseudo-R 2 criterion to evaluate PRS predictive accuracy for time-to-event data. This new criterion is related to the score statistic derived under a two-component mixture model. It evaluates the effect of the PRS on both the propensity to experience the event and on the dynamic of the event among the susceptible subjects. Simulation results show that our index has good properties. We compared our index to other implemented pseudo-R 2 for survival data. Along with our index, two other indices have comparable good behavior when the PRS has a non-null propensity effect, and our index is the only one to detect when the PRS has only a dynamic effect. We evaluated the 5-year predictivity of an 18-single-nucleotide-polymorphism PRS for incident breast cancer cases on the CARTaGENE cohort using several pseudo-R 2 indices. We report that our index, which summarizes both a propensity and a dynamic effect, had the highest predictive accuracy. In conclusion, our proposed pseudo-R 2 is easy to implement and well suited to evaluate PRS for predicting incident events in cohort studies.

A Machine Learning Approach for High-Dimensional Time-to-Event Prediction With Application to Immunogenicity of Biotherapies in the ABIRISK Cohort.

Predicting immunogenicity for biotherapies using patient and drug-related factors represents nowadays a challenging issue. With the growing ability to collect massive amount of data, machine learning algorithms can provide efficient predictive tools. From the bio-clinical data collected in the multi-cohort of autoimmune diseases treated with biotherapies from the ABIRISK consortium, we evaluated the predictive power of a custom-built random survival forest for predicting the occurrence of anti-drug antibodies. This procedure takes into account the existence of a population composed of immune-reactive and immune-tolerant subjects as well as the existence of a tiny expected proportion of relevant predictive variables. The practical application to the ABIRISK cohort shows that this approach provides a good predictive accuracy that outperforms the classical survival random forest procedure. Moreover, the individual predicted probabilities allow to separate high and low risk group of patients. To our best knowledge, this is the first study to evaluate the use of machine learning procedures to predict biotherapy immunogenicity based on bioclinical information. It seems that such approach may have potential to provide useful information for the clinical practice of stratifying patients before receiving a biotherapy.