A METHOD OF PREDICTING RESPONSE TO TREATMENT WITH A DISEASE-MODIFYING ANTI-RHEUMATIC DRUG, AND/OR CLASSIFYING DISEASE ACTIVITY IN A SUBJECT WITH RHEUMATOID ARTHRITIS

Abstract

The present invention relates to a method of predicting response to treatment with a disease-modifying anti-rheumatic drug, and/or classifying disease activity in a subject with rheumatoid arthritis by determining the quantitative level of one or more biomarkers and either predicting response to treatment with a disease-modifying anti-rheumatic drug, classifying disease activity, and/el or predicting response to treatment with a disease-modifying anti-rheumatic drug and classifying disease activity.

Description

FIELD OF THE INVENTION

The present invention relates to a method of predicting response to treatment with a disease-modifying anti-rheumatic drug, and/or classifying disease activity in a subject with rheumatoid arthritis. The method comprises determining the quantitative level of one or more biomarkers in a biological sample from the subject and either predicting response to treatment with a disease-modifying anti-rheumatic drug, classifying disease activity, and/or predicting response to treatment with a disease-modifying anti-rheumatic drug and classifying disease activity; based on the quantitative level of the or each biomarker in the biological sample.

BACKGROUND TO THE INVENTION

Many studies have highlighted the need for robust biomarkers in rheumatoid arthritis (RA) in order to assist with diagnosis, prognosis and determining treatment response. The human plasma proteome is of particular interest as it contains proteins secreted or shed by circulating cells and tissue involved in the pathogenesis of RA. When the body is under stress, for example during chronic inflammation, the plasma proteome changes as part of the immune response. This study focusses on identifying and quantifying changes in circulating plasma proteins between patients with variable treatment responses, as this approach can also improve our understanding of disease pathogenesis.

Disease-modifying anti-rheumatic drugs (DMARDs) are immunosuppressive and immunomodulatory agents, and are classified as either conventional DMARDs (cDMARDs) or biologic DMARDs (bDMARDs). Commonly used conventional DMARDs include methotrexate, leflunomide, hydroxychloroquine, and sulfasalazine. Biologic DMARDs were introduced in the early 1990s and are usually prescribed after the failure of conventional DMARD therapy (ongoing disease activity, or clinical or radiographic disease progression). Some biologic agents include infliximab, adalimumab, etanercept, rituximab, abatacept, rituximab, tocilizumab, tofacitinib, among others. Biologic DMARDs are highly specific and target a specific pathway of the immune system.

Currently, clinicians consider C-reactive protein (CRP) a biomarker of inflammation. Although NICE guidelines recommend regular monitoring of CRP levels throughout treatment of RA (NICE 2009), some studies suggest it can be used as a substitute of erythrocyte sedimentation rate (ESR) in composite scores used to determine treatment response such as the Disease Activity Score in 28 joints (DAS28). However, the increase in plasma CRP levels observed in RA compared to health is caused in part by the secretion of IL-6 from T-cells and macrophages during inflammation and conflicting evidence brings into question CRP's reliability as a marker of joint inflammation or disease activity. The lack of specificity of CRP to RA pathology may be a contributing factor to this contradictory evidence.

The wide variety of proteomic techniques now available has enabled research groups to robustly analyse a number of biological sample types. The most obvious, pathologically proximal sample to study disease activity in RA patients is synovial fluid. For example, a study using isobaric tag labelling of synovial peptides identified over 500 proteins that were differentially expressed between osteoarthritis (OA) and RA patients.

Cytosolic proteins isolated from synovial tissue and analysed by 2-dimensional gel electrophoresis discovered protein expression patterns which differentiate RA from spondyloarthropathy, enabling earlier diagnosis. Similar studies have used the same method to analyse proteins from fibroblast-like synoviocyte (FLS) cells and have also identified potential novel diagnostic proteins for RA. The analysis of proteins from synovial fluid or tissue may be useful in understanding disease pathogenesis, as this is the target of the autoimmune response in RA. However, this sampling method is rather invasive, and it is an impractical sample in comparator healthy control groups. Although the study of proteins from tissue-specific cells is advantageous in understanding the cell biology of RA, the lengthy laboratory procedures are somewhat unrealistic for regular determination of treatment response.

An increasing number of studies have investigated potential circulating biomarkers of treatment response in RA. Some studies have looked at circulating cellular markers of response. Peripheral blood is a potentially rich source of biomarker candidates to determine treatment response in RA. However, previous studies have failed to identify robust proteins that are indicative of both disease activity and treatment response. This study focuses on the circulating plasma proteome in RA as an accessible source of potential response biomarkers.

SUMMARY OF THE INVENTION

According to a first aspect of the present invention, there is provided a method of predicting response to treatment with a disease-modifying anti-rheumatic drug (DMARD) in a subject, the method comprising the steps of:

• • (a) determining the quantitative level of one or more biomarkers in a biological sample from the subject; and
  • (b) predicting response to treatment with a DMARD in the subject based on the quantitative level of the or each biomarker in the biological sample;
  • wherein the biomarker is DECR1.

Optionally or additionally, the method is a method of classifying disease activity in a subject with rheumatoid arthritis.

Optionally, there is provided a method of predicting response to treatment with a disease-modifying anti-rheumatic drug (DMARD), and/or classifying disease activity in a subject with rheumatoid arthritis, the method comprising the steps of:

• • (a) determining the quantitative level of one or more biomarkers in a biological sample from the subject; and
  • (b) predicting response to treatment with a DMARD, and/or classifying disease activity based on the quantitative level of the or each biomarker in the biological sample;wherein the biomarker is DECR1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises CD40-L.

Preferably, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1 and CD40-L.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises ITGB1 BP2.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L and ITGB1 BP2.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises CASP-3.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2 and CASP-3.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises PDFG subunit A.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3 and PDFG subunit A.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises Dkk-1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A and Dkk-1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises SELP.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1 and SELP.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises SOD2.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP and SOD2.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises SORT1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2 and SORT1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity further comprises CD84.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Optionally, there is provided a method of predicting response to treatment with a disease-modifying anti-rheumatic drug (DMARD), and/or classifying disease activity in a subject with rheumatoid arthritis, the method comprising the steps of:

• • (a) determining the quantitative level of one or more biomarkers in a biological sample from the subject; and
  • (b) predicting response to treatment with a DMARD, and/or classifying disease activity based on the quantitative level of the or each biomarker in the biological sample;wherein the or each biomarker is selected from DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Optionally, the or each biomarker is a gene.

Optionally, the or each biomarker is a nucleic acid.

Optionally, the or each biomarker is a deoxyribonucleic acid.

Optionally, the or each biomarker is a ribonucleic acid.

Preferably, the or each biomarker is a protein.

Optionally, the or each biomarker is a peptide.

Optionally, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is a gene having a GenBank Accession Version Number selected from one or more of AK300069.1, L07414.1, BC108901.2, AJ413269.1, AH002928.2, AH009834.2, AL022146.1, X07834.1, AL390252.9 and AF054816.1.

Preferably, the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is a protein having a UniProt ID selected from one or more of Q16698, P29965, Q9UKP3, P42574, P04085, 094907, P16109, P04179, Q99523 and Q9U1B8.

Optionally, the or each biomarker for predicting response to treatment with a DMARD further comprises CASP-3, SELP, JAM-A, PDFG subunit A, ITGB1 BP2, Dkk-1, PECAM-1, ANG-1, SOD2, DAPP1, DCTN1, PAI and CD40-L.

Optionally, the or each biomarker for predicting response to treatment with a DMARD is selected from one or more of DECR1, CASP-3, SELP, JAM-A, PDGF subunit A, ITGB1 BP2, Dkk-1, PECAM-1, ANG-1, SOD2, DAPP1, DCTN1, PAI and CD40-L.

Optionally, the or each biomarker for predicting response to treatment with a DMARD further comprises HB-EGF, HCLS1, IRAK4, STK4, SH2B3, NEMO, CD84, EIF4G1, PRDX1, TANK, CXCL1, SIRT2, ZBTB16, NF2, PDGF subunit B, SPRY2, HEXIM1, PPP1R9B, PRDX3, PAR-1, ICA1, CCL17, STAMPB, MGMT, SRPK2, BML hydrolase, PLXNA4 and SORT1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD is selected from one or more of DECR1, CASP-3, SELP, JAM-A, PDGF subunit A, ITGB1 BP2, Dkk-1, PECAM-1, ANG-1, SOD2, DAPP1, DCTN1, PAI, CD40-L, HB-EGF, HCLS1, IRAK4, STK4, SH2B3, NEMO, CD84, EIF4G1, PRDX1, TANK, CXCL1, SIRT2, ZBTB16, NF2, PDGF subunit B, SPRY2, HEXIM1, PPP1 R9B, PRDX3, PAR-1, ICA1, CCL17, STAMPB, MGMT, SRPK2, BML hydrolase, PLXNA4 and SORT1.

Optionally, the or each biomarker for predicting response to treatment with a DMARD is a gene having a GenBank Accession Version Number selected from one or more of AK300069.1, AJ413269.1, AL022146.1, AL136649.1, AH002928.2, BC108901.2, BC107046.2, JQ287500.1, AY124380.1, X07834.1, AF178987.1 AK314352.1, X04731.1, AK222896.1, AC004634.1, AK312750.1, AY340965.1, Z93016.2, BC136451.1, AJ271718.1, AF054816.1, AK131407.1, CR407652.1, AC009299.5, BT006880.1, KF032391.1, BC029812.1, Y18000.1, M12783.1, AL354668.13, BC006460.1, KF495721.1, BT020007.1, BC002464.2, U37183.1, CH471092.1, CH471053.2, M60761.1, AC004884.1, AK312896.1, AC009365.9 and AL390252.9.

Preferably, the or each biomarker for predicting response to treatment with a DMARD is a protein having a UniProt ID selected from one or more of Q16698, P42574, P16109, Q9Y624, P04085, Q9UKP3, Q9UBT3, P16284, Q15389, P04179, Q9UN19, Q14203, P05121, P25942, Q99075, P14317, Q9NWZ3, Q13043, Q9UQQ2, Q9Y6K9, Q9U1B8, Q04637, Q06830, Q92844, P09341, Q8IXJ6, Q05516, P35240, P01127, 043597, 094992, Q96SB3, P30048, P25116, Q05084, Q92583, 095630, P16455, P78362, Q13867, Q9HCM2 and Q99523.

Preferably, the or each biomarker for predicting treatment response in a subject treated with a DMARD is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Optionally, the method is a method of classifying disease activity in a subject with rheumatoid arthritis.

Optionally, there is provided a method of classifying disease activity in a subject with rheumatoid arthritis, the method comprising the steps of:

• • (a) determining the quantitative level of one or more biomarkers in a biological sample from the subject; and
  • (b) classifying disease activity based on the quantitative level of the or each biomarker in the biological sample;
  • wherein the biomarker is DECR1.

Optionally, the or each biomarker for classifying disease activity further comprises CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Optionally, the or each biomarker for classifying disease activity is selected from one or more of DECR1, CD40-L, ITGB1 BP2, CASP-3, PDFG subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Optionally, the or each biomarker for classifying disease activity is a gene having a GenBank Accession Version Number selected from one or more of AK300069.1, L07414.1, BC108901.2, AJ413269.1, AH002928.2, AH009834.2, AL022146.1, X07834.1, AL390252.9 and AF054816.1.

Preferably, the or each biomarker for classifying disease activity is a protein having a UniProt ID selected from one or more of Q16698, P29965, Q9UKP3, P42574, P04085, 094907, P16109, P04179, Q99523 and Q9UIB8.

Optionally, the determining step (a) comprises determining the quantitative level of one or more biomarkers in the biological sample from the subject.

Further optionally, the determining step (a) comprises determining the quantitative level of two, three, four, five, six, seven, eight, nine, ten, eleven, twelve, thirteen, fourteen, fifteen, twenty, twenty five, thirty, thirty five, forty, forty two or more biomarkers in the biological sample from the subject.

Optionally, the determining step (a) comprises determining the quantitative level of all of the biomarkers in the biological sample from the subject.

Optionally or additionally, the determining step (a) comprises determining the quantitative level of each of the biomarkers in the biological sample from the subject.

Optionally, the biological sample is selected from whole blood, serum, plasma, urine, interstitial fluid, peritoneal fluid, cervical swab, tears, saliva, buccal swab, skin, brain tissue, and cerebrospinal fluid.

Further optionally, the biological sample is selected from whole blood, serum, and plasma.

Preferably, the biological sample is whole blood.

Further preferably, the biological sample is plasma.

Optionally, the predicting and/or classifying step (b) comprises comparing the quantitative level of the or each biomarker in the biological sample from the subject with the quantitative level of the or each respective biomarkers in a normal sample.

Optionally, the normal sample is a biological sample from a subject not treated with a DMARD.

Optionally, the normal sample is a biological sample from a subject not treated with a DMARD and not suffering from rheumatoid arthritis.

Optionally, the subject is treated with one or more DMARDs Optionally, the subject is treated with one or more DMARDs, including any one or more of infliximab, adalimumab, etanercept, rituximab, abatacept, rituximab, tocilizumab, certolizumab, golimumab, anakinra, baricitinib, upadacitinab and tofacitinib.

Optionally, the subject is treated with one or more DMARDs, including any one or more of infliximab, adalimumab, etanercept, rituximab, abatacept, rituximab, tocilizumab, certolizumab, golimumab, anakinra, baricitinib, upadacitinab and tofacitinib, and will respond to treatment (responder).

Optionally, in the method of predicting response to treatment with a DMARD, the subject is treated with one or more DMARDs including any one or more of infliximab, adalimumab, etanercept, rituximab, abatacept, rituximab, tocilizumab, certolizumab, golimumab, anakinra, baricitinib, upadacitinab and tofacitinib, and will not respond to treatment (non-responder).

Optionally, the method is a method of predicting response to treatment with a conventional disease-modifying anti-rheumatic drug (cDMARD).

Optionally, the subject is treated with one or more cDMARDs.

Optionally, the subject is treated with one or more cDMARDs, including any one or more of apremilast, azathioprine, ciclosporin, cyclophosphamide, sodium aurothiomalate, mycophenolate mofetil, penicillamine, tacrolimus, methotrexate, sulfasalazine, leflunomide and hydroxychloroquine.

Optionally, the subject is treated with one or more cDMARDs, including any one or more of apremilast, azathioprine, ciclosporin, cyclophosphamide, sodium aurothiomalate, mycophenolate mofetil, penicillamine, tacrolimus, methotrexate, sulfasalazine, leflunomide and hydroxychloroquine, and will respond to treatment (responder).

Optionally, in the method of predicting response to treatment with a cDMARD, the subject is treated with one or more cDMARDs including apremilast, azathioprine, ciclosporin, cyclophosphamide, sodium aurothiomalate, mycophenolate mofetil, penicillamine, tacrolimus, methotrexate, sulfasalazine, leflunomide and hydroxychloroquine, and will not respond to treatment (non-responder).

Optionally, the determining step (a) further comprises determining the quantitative level in a first set of respective biomarkers(s).

Optionally, the determining step (a) further comprises determining the quantitative level in a first set of respective biomarkers(s) wherein the quantitative level of the first set of biomarkers is greater than the quantitative level of the respective biomarkers(s) in a normal sample, which is indicative of a non-responder.

Optionally, the first set of biomarkers is selected from any one or more of DECR1, CASP-3, SELP, JAM-A, PDFG subunit A, ITGB1BP2, Dkk-1, PECAM-1, ANG-1, SOD2, DAPP1, DCTN1, PAI, CD40-L, HB-EGF, HCLS1, IRAK4, STK4, SH2B3, NEMO, CD84, EIF4G1, PRDX1, TANK, CXCL1, SIRT2, ZBTB16, NF2, PDGF subunit B, SPRY2, HEXIM1, PPP1R9B, PRDX3, PAR-1, ICA1, CCL17, STAMPB, MGMT, SRPK2, BML hydrolase, PLXNA4 and SORT1.

Optionally, the determining step (a) further comprises determining the quantitative level in a second set of respective biomarkers(s).

Optionally, the determining step (a) further comprises determining the quantitative level in a second set of respective biomarkers(s) wherein the quantitative level of the second set of biomarkers is less than the quantitative level of the respective biomarkers(s) in a normal sample, which is indicative of a non-responder.

Optionally, the second set of biomarkers is selected from any one or more of DECR1, CD40-L, ITGB1BP2, CASP-3, PDGF subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

Further optionally, the quantitative level of the or each biomarker is higher than a threshold value and is indicative of a non-responder, wherein the respective threshold value of the or each biomarker is:

Biomarker      Threshold value
DECR1          7.1320
CASP-3         10.4740
SELP           10.5279
PDGF subunit A 5.5889
ITGB1BP2       6.2200
Dkk-1          10.4467
SOD2           9.9025
CD40-L         7.4976
CD84           6.0718
SORT1          8.5585

Optionally, the quantitative level of the or each biomarker is lower than a threshold value and is indicative of a responder, wherein the respective threshold value of the or each biomarker is:

Biomarker      Threshold value
DECR1          7.1320
CASP-3         10.4740
SELP           10.5279
PDGF subunit A 5.5889
ITGB1BP2       6.2200
Dkk-1          10.4467
SOD2           9.9025
CD40-L         7.4976
CD84           6.0718
SORT1          8.5585

Optionally, the quantitative level of the or each biomarker is between two threshold values and is indicative of high disease activity, wherein the respective threshold values of the or each biomarker are:

Biomarker      Threshold value
DECR1          7.8984-11.6893
CD40-L         8.4728-12.7503
ITGB1BP2       6.7802-11.4759
CASP-3         11.4114-13.6768
PDFG subunit A 5.7334-8.4577
Dkk-1          11.2034-12.7568
SELP           11.3715-12.6136
SOD2           10.0471-10.5467
SORT1          8.7990-9.7006
CD84           6.1830-7.8683

Optionally, the quantitative level of the or each biomarker is between two threshold values and is indicative of moderate disease activity, wherein the respective threshold values of the or each biomarker are:

Biomarker      Threshold values
DECR1          5.6510-10.5039
CD40-L         6.1622-10.1159
ITGB1BP2       3.1210-8.5301
CASP-3         8.0945-12.2662
PDFG subunit A 3.9517-7.0510
Dkk-1          9.1861-12.0037
SELP           9.2287-11.9099
SOD2           9.5039-10.2607
SORT1          7.5225-9.2595
CD84           4.9694-6.8399

Optionally, the quantitative level of the or each biomarker is between two threshold values and is indicative of low disease activity, wherein the respective threshold values of the or each biomarker are:

Biomarker      Threshold values
DECR1          6.6490-6.7309
CD40-L         6.5871-6.9401
ITGB1BP2       4.1040-4.5746
CASP-3         8.3368-8.7173
PDFG subunit A 3.8283-4.0907
Dkk-1          9.6015-9.8410
SELP           8.2174-9.3414
SOD2           9.6542-9.9702
SORT1          7.7980-8.5272
CD84           5.2689-5.5881

Optionally, the quantitative level of the or each biomarker is between two threshold values and is indicative of in remission, wherein the respective threshold values of the or each biomarker are:

Biomarker      Threshold value
DECR1          5.6064-6.9173
CD40-L         6.2510-7.3891
ITGB1BP2       3.4508-5.5040
CASP-3         7.9543-9.2353
PDFG subunit A 2.7648-4.9645
Dkk-1          9.2598-9.9684
SELP           8.8272-9.9551
SOD2           9.5195-9.7582
SORT1          7.6829-8.4218
CD84           4.6188-5.8750

Optionally, the method is a method of classifying disease activity, wherein the subject has high disease activity, moderate disease activity, low disease activity or is in remission.

Optionally, in a subject with high disease activity, the quantitative level of the or each biomarker is higher compared to a subject with moderate disease activity, low disease activity and in remission.

Optionally, in a subject with moderate disease activity, the quantitative level of the or each biomarker is higher compared to a subject with low disease activity and in remission, but the quantitative level of the or each biomarker is lower compared to a subject with high disease activity.

Optionally, in a subject with low disease activity, the quantitative level of the or each biomarker is higher compared to a subject with in remission, but the quantitative level of the or each biomarker is lower compared to a subject with high disease activity and moderate disease activity.

Optionally, in a subject in remission, the quantitative level of the or each biomarker is lower compared to a subject with high disease activity, moderate disease activity and low disease activity.

Optionally, the method of classifying disease activity is a method of prognosing disease activity.

Optionally, the method of predicting response to treatment with a DMARD, and/or classifying disease activity is an in vitro method.

Optionally, the method of predicting response to treatment with a cDMARD, and/or classifying disease activity is an in vitro method.

Optionally, the method of predicting response to treatment with a DMARD is an in vitro method.

Optionally, the method of predicting response to treatment with a cDMARD is an in vitro method.

Optionally, the method of classifying disease activity is an in vitro method.

Optionally, the method of prognosing disease activity is an in vitro method.

According to a further aspect of the present invention, there is provided a method of treating rheumatoid arthritis by predicting response to treatment with a DMARD; and treating the subject.

Optionally, there is provided a method of treating rheumatoid arthritis by predicting response to treatment with a cDMARD; and treating the subject.

Optionally, there is provided a method of treating rheumatoid arthritis by classifying disease activity; and treating the subject.

Optionally, there is provided a method of treating rheumatoid arthritis by predicting response to treatment with a DMARD, and/or classifying disease activity; and treating the subject.

Optionally, there is provided a method of treating rheumatoid arthritis by predicting response to treatment with a cDMARD, and/or classifying disease activity; and treating the subject.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the present invention will now be described with reference to the accompanying drawings and following non-limiting examples, in which:

FIG. 1 illustrates unsupervised clustering of proteins calculated to be significantly different between DMARD naïve (normal) (N), DMARD responder (DR) and DMARD non-responder (DNR) RA patients;

FIG. 2 illustrates unsupervised clustering of proteins calculated to be significantly different between RA patients with high, moderate or low disease activity, or who are in remission. Disease activity classifications are based on EULAR guidelines;

FIG. 3 illustrates ROC curve analysis of each protein when a cut-off (y-intercept) value is applied to distinguish RA cDMARD responders and non-responders; and

FIG. 4 illustrates correlation of NPX values of top ten proteins of interest with DAS28-ESR score.

EXAMPLES

Materials & Methods

The research team at Ulster University collaborated with rheumatologists from the Western Health and Social Care Trust and the Belfast Health and Social Care Trust to design, conduct and recruit patients to the study. Informed consent to participate was obtained from all RA patients enrolled in the study (Table 1).

TABLE 1
Demographics of patients treated with conventional disease-modifying anti-rheumatic drugs
DN = DMARD naïve (normal), DR = DMARD responder, DNR = DMARD non-responder.
	DN	DR	DNR
	(n = 6)	(n = 24)	(n = 32)
Female, n (%)	5 (83.3)	22	(91.7)	27 (81.8)
Age, mean (SD), years	42.5 (16.7)	59.7	(12.2)	57.7 (10.1)
Treatment duration since diagnosis, mean (SD), years	N/A	5.0	(3.8)	10.8 (12.1)
Erythrocyte Sedimentation Rate (ESR), mean (SD), mm in 1 hr	23.5 (15.2), n = 6	12.2 (7.9), n = 24	24.1 (20.0), n = 29
C-reactive Protein (CRP), mean (SD), mg/L	6.1 (2.9), n = 6	5.3 (7.3), n = 22	10.6 (9.4), n = 28
Disease Activity Score in 28 joints (DAS28-ESR), mean (SD)	4.3 (0.4), n = 4	2.6 (1.0), n = 24	4.9 (1.1), n = 19

Diagnosis for all patients fulfilled the 2010 American College of Rheumatology (ACR) criteria (American College of Rheumatology and European League Against Rheumatism 2010). Patients were treated with one or more conventional disease-modifying anti-rheumatic drugs (cDMARD), including methotrexate, sulfasalazine, leflunomide and hydroxychloroquine.

Patients were classified as ‘responders’ or ‘non-responders’ according to the National Institute for Health and Clinical Excellence (NICE) guidelines, where responders have exhibited a change in DAS28-ESR of >1.2 following treatment (NICE 2009). Responders and non-responders had undergone treatment for an average of 5 and 11 years respectively.

Disease activity classifications were according to the European League Against Rheumatism (EULAR) criteria, where a DAS28-ESR score of ≥5.1 was classified as ‘high’ disease activity (n=9), >3.2 and ≤5.1 was ‘moderate’ (n=21), >2.6 and ≤3.2 was ‘low’ (n=2), and <2.6 was classified as ‘remission’ (n=11).

Office for Research Ethics Committees Northern Ireland (ORECNI) (11/NI/0188) and Ulster University Research Ethics Committee (UREC) (REC/14/0053) approvals were obtained for the study.

Blood Processing

Peripheral whole blood used for plasma extraction was collected in an EDTA coated tube (Aquilant Scientific) using a venepuncture technique. Blood tubes were stored at 4° C. for a maximum of 2 hours until processing. The blood tube was centrifuged in an Eppendorf 5804 centrifuge at 300 rcf for 10 minutes at 18° C. The majority of the upper plasma layer was removed and stored at −80° C. for later analysis.

Proximity Extension Assay (PEA)

Plasma from each patient sample was thawed and a small volume of ˜40 μl was sent to OLINK (Analysis service Uppsala, Sweden). Each sample was analysed across OLINK's CVD II, CVD III, immune response and inflammatory panels. Each high throughput multiplex immunoassay consisted of 92 proteins.

The technology used is known as ‘Proximity Extension Assay’ (PEA), where matched antibody pairs containing unique DNA sequences were incubated with the sample. Sequences that were exactly matched hybridized and were extended using DNA polymerase. The elongated DNA was then amplified and a qPCR readout using Fluidigm® BioMARK™ was used to determine the results of each protein target. The assay included internal controls at each stage of the reaction, negative controls to determine background values, and inter-plate controls to account for variability between plates.

ΔCt values were log normalised to what is called ‘normalised protein expression’ values or NPX. Perseus software (version 1.6.0.7) was used to analyse the vast range of NPX values across patient groups. A Z-score was performed across the entire data set in order to standardise the NPX values.

In order to filter data based on significance of differentially expressed proteins across all patient groups, an analysis of variance (ANOVA) test was carried out with permutation-based false discovery rate (FDR). These proteins were expressed in a heat map with unsupervised clustering.

Example 1

Proteins that are Significantly Differentially Expressed Between Each Patient Group when Classified by Response, as Calculated by ANOVA

NPX values from all four OLINK panels across each patient group were analysed using Perseus software.

TABLE 2
Statistically significant differences in detected levels
of plasma proteins were calculated between each patient
group using an ANOVA test with permutation based FDR, where
FDR = 0.05 (with the number of randomizations set at 250).
	Difference
	(Non-	−log ANOVA
	responder −	p value (all	ANOVA	FDR
Protein	Responder)	patient groups)	p value	q value
CASP-3	1.654	12.924	<0.0001	<0.001
SELP	1.571	11.523	<0.0001	<0.001
JAM-A	1.589	11.496	<0.0001	<0.001
PDGF subunit A	1.591	10.616	<0.0001	<0.001
ITGB1BP2	1.391	8.522	<0.0001	<0.001
Dkk-1	1.358	8.010	<0.0001	<0.001
PECAM-1	1.363	7.441	<0.0001	<0.001
ANG-1	1.326	7.229	<0.0001	<0.001
DECR1	1.267	6.782	<0.0001	<0.001
SOD2	1.269	6.529	<0.0001	<0.001
DAPP1	1.353	6.492	<0.0001	<0.001
DCTN1	1.288	6.278	<0.0001	<0.001
PAI	1.300	6.022	<0.0001	<0.001
CD40-L	1.220	5.899	<0.0001	<0.001
HB-EGF	1.169	5.446	<0.0001	0.001
HCLS1	1.273	5.285	<0.0001	0.002
IRAK4	1.173	5.008	<0.0001	0.002
STK4	1.214	4.892	<0.0001	0.002
SH2B3	1.196	4.743	<0.0001	0.002
NEMO	1.034	4.686	<0.0001	0.002
CD84	1.065	4.605	<0.0001	0.002
EIF4G1	1.175	4.584	<0.0001	0.002
PRDX1	−0.945	4.435	<0.0001	0.002
TANK	1.073	4.254	0.0001	0.004
CXCL1	1.036	4.251	0.0001	0.004
SIRT2	1.122	4.168	0.0001	0.005
ZBTB16	1.028	3.645	0.0002	0.010
NF2	1.112	3.639	0.0002	0.009
PDGF subunit B	1.042	3.536	0.0003	0.010
SPRY2	1.098	3.498	0.0003	0.011
HEXIM1	1.063	3.444	0.0004	0.011
PPP1R9B	1.048	3.424	0.0004	0.012
PRDX3	1.081	3.382	0.0004	0.013
PAR-1	0.866	3.241	0.0006	0.016
ICA1	0.985	3.217	0.0006	0.017
CCL17	0.837	3.183	0.0007	0.017
STAMPB	0.951	3.023	0.0009	0.021
MGMT	0.943	2.835	0.0015	0.026
SRPK2	0.972	2.757	0.0017	0.028
BLM hydrolase	0.945	2.678	0.0021	0.033
PLXNA4	0.951	2.630	0.0023	0.035
SORT1	0.790	2.612	0.0024	0.035

Example 2

Unsupervised Clustering of Proteins Calculated to be Significantly Different Between DMARD Naïve (Normal) (N), DMARD Responder (DR) and DMARD Non-Responder (DNR) RA Patients

FIG. 1 illustrates a heat map showing only the proteins that were significantly different between groups, having an FDR q value of <0.05. Of the 42 proteins calculated as significantly different between patient groups, the majority were upregulated in cDMARD non-responders, compared to responders and cDMARD naïve (normal) patients. The unsupervised clustering of the heat map shows a split between a) cDMARD non-responders b) cDMARD naïve (normal) and c) cDMARD responders, where just one cDMARD non-responder is clustered among the responders. Furthermore, one cDMARD non-responder shows a reduced protein expression for the majority of the 42 proteins graphed, compared to all other patient samples. The reason for this outlier is unclear from available clinical data.

Example 3

Proteins that are Significantly Differentially Expressed Between Each Patient Group when Classified by Disease Activity Score (DAS), as Calculated by ANOVA.

DAS28-ESR scores from clinical data recorded on the same day as sampling were classified according to EULAR criteria (NICE 2009). A DAS28-ESR score of ≤5.1 was classified as ‘high’ disease activity, >3.2 and ≤5.1 was ‘moderate’, >2.6 and ≤3.2 was ‘low’, and <2.6 was classified as ‘remission’.

TABLE 3
Each group was analysed by ANOVA
with permutation based FDR.
	Difference
	(High disease	−log ANOVA
	activity −	p value (all	ANOVA	FDR
Protein	Remission)	patient groups)	p value	q value
ITGB1BP2	2.015	5.514	<0.0001	0.200
CD84	2.005	5.375	<0.0001	0.110
Dkk-1	1.990	5.355	<0.0001	0.073
SOD2	2.039	5.200	<0.0001	0.057
SORT1	1.961	5.158	<0.0001	0.046
DECR1	1.959	4.988	<0.0001	0.051
CASP-3	1.813	4.985	<0.0001	0.044
CD40-L	1.882	4.724	<0.0001	0.044
PDGF subunit A	1.737	4.616	<0.0001	0.043
SELP	1.611	4.570	<0.0001	0.040

Example 4

Unsupervised Clustering of Proteins Calculated to be Significantly Different Between RA Patients with High, Moderate or Low Disease Activity, or Who are in Remission

FIG. 2 illustrates a heat map with unsupervised clustering showing only the proteins that were significantly different between groups. Of the four protein panels analysed, only ten were statistically significantly different in normalised levels between patient groups.

Example 5

Analysis of the Top 10 Significantly Different Proteins

Each of the ten proteins were also significantly different between patient groups when classified by response. Statistical data of the ten proteins is summarised in Table 4 (see also FIG. 3 and FIG. 4). Proteins were ranked from left to right with the best performing protein on the left, based on AUC score. The benchmark false discovery AUC for this analysis was 0.6251.

TABLE 4
Statistical analysis of top ten proteins that show significant difference
across RA patients when grouped by response or DAS28-ESR.
Individual
Protein					PDGF
Performance	DECR1	CD40-L	ITGB1BP2	CASP-3	subunit A	Dkk-1	SELP	SOD2	SORT1	CD84
p-value	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	0.0062	0.0001
AUC	1.0000	1.0000	0.9758	0.9550	0.9481	0.9412	0.9412	0.9308	0.9308	0.9204
Sensitivity	1.0000	1.0000	0.9412	1.0000	1.0000	1.0000	1.0000	0.9412	0.9412	0.8824
Specificity	1.0000	1.0000	0.9412	0.9412	0.9412	0.9412	0.9412	0.8824	0.8824	0.8824
Threshold	7.1320	7.4976	6.2200	10.4740	5.5889	10.4467	10.5279	9.9025	8.5585	6.0718
Correlation	0.0004	0.0002	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001
with DAS28-
ESR (p-value)
Correlation	0.6010	0.6273	0.6577	0.7118	0.6971	0.6853	0.6700	0.6515	0.6566	0.6748
with DAS28-
ESR (r-value)

Proteins are ranked from left to right, with the best performing protein on the left, based on AUC score. ROC curve analysis of each protein when a cut-off (y-intercept) value is applied to distinguish RA cDMARD responders and non-responders (FIG. 3). Also shown is correlation data of each NPX value with DAS28-ESR score (FIG. 4). We considered the effect of using combinations of proteins for classification. Proteins acting in combination all achieved an AUC score of 1. In order to establish whether this was due to overfitting, we generated nonsense data by randomly swapping patients between the responder and non-responder groups sufficiently to obtain decorrelated nonsense data. We then optimised classifiers to the nonsense data using the same protein combinations and calculated the corresponding AUC score. This was repeated twenty times and the mean AUC was calculated. The results can be seen in Table 5 where the classification AUC scores obtained using true data significantly outperforms the overfitting rate which is equal to the mean classification AUCs obtained using the nonsense data.

TABLE 5
Statistical analysis of the collective protein performance of the top ten proteins that
show significant difference across RA patients when grouped by response or DAS28-ESR.
Collective
Protein					PDGF
Performance	DECR1	CD40-L	ITGB1BP2	CASP-3	subunit A	Dkk-1	SELP	SOD2	SORT1	CD84
AUC	1	1	1	1	1	1	1	1	1	1
Sensitivity	1	1	1	1	1	1	1	1	1	1
Specificity	1	1	1	1	1	1	1	1	1	1
Benchmark AUC	0.625
Benchmark SE	0.013
Benchmark	0.541	0.5411	0.5411	0.541	0.541	0.541	0.54	0.541	0.541	0.541
Sensitivity
Benchmark	0.617	0.6176	0.6176	0.617	0.6176	0.6176	0.617	0.6176	0.6176	0.617
Specificity
Threshold	−28.84	−3.120	−2.540	−1.31	−1.2572	−1.275	−1.28	−1.2760	−0.680	−0.605
	205.70	−4.427	−2.3013	−0.55	−0.439	−0.472	−0.46	−0.4470	−0.315	−0.325
	0	55.773	−2.0427	0.1429	0.3053	0.2642	0.283	0.3107	0.0233	−0.0574
	0	0	49.256	−0.913	−0.560	−0.494	−0.332	−0.3449	−0.3096	−0.3124
	0	0	0	22.535	−0.518	−0.495	−0.378	−0.3728	−0.3396	−0.3677
	0	0	0	0	19.477	−0.125	−0.024	−0.0450	−0.0575	−0.0381
	0	0	0	0	0	20.538	−0.292	−0.301	−0.309	−0.322
	0	0	0	0	0	0	20.122	0.0805	0.0724	0.1029
	0	0	0	0	0	0	0	19.349	0.6848	0.5951
	0	0	0	0	0	0	0	0	9.8592	0.2706
	0	0	0	0	0	0	0	0	0	8.8411

Each column presents the results for the classifier obtained for the protein labelled in the column in combination with any proteins to the left of the column. As we move to the right of the table, incrementally more proteins are included in the classifier. AUC, sensitivity and specificity are shown for the protein combinations applied to true data. Also shown are the overfitting rate (mean benchmark AUC) and the standard error in the benchmark AUC, obtained from twenty iterations of calculating nonsense data, the mean sensitivity and specificity of the classifiers obtained for nonsense data and the thresholds identified which grow in dimension with increasing numbers of proteins.

One of the key aims of this work was to discover circulating protein biomarkers that may aid in differentiating clinical response to cDMARD treatment in RA. Of the 4 protein panels tested from OLINK (total number=384), 42 proteins were significantly differentially expressed between cDMARD naïve (normal), cDMARD responder and cDMARD non-responder RA patients. When classifying patients based on disease activity as determined by EULAR criteria, 10 proteins were shown to be differentially expressed between patient groups. Interestingly, these 10 proteins were also included in the 42 proteins previously shown to be differentially expressed when patients were classified by response. This suggests that these 10 proteins have the potential of distinguishing both disease activity and response to cDMARD treatment.

When a threshold is applied, each of the individual 10 proteins demonstrates good sensitivity (>88%) and specificity (>88%) in distinguishing responders from non-responders.

Overall, considering the results of the individual protein analysis and the combined protein analysis, DECR1 and CD40-L are the most optimal individual proteins in distinguishing response. No incremental benefit was observed when combining these proteins with each other or the others on the list (Table 5).

Claims

1. An in vitro method of predicting response to treatment with a disease-modifying anti-rheumatic drug (DMARD), and/or classifying disease activity in a subject with rheumatoid arthritis, the method comprising the steps of: (a) determining the quantitative level of one or more biomarkers in a biological sample from the subject; and(b) predicting response to treatment with a DMARD, and/or classifying disease activity based on the quantitative level of the or each biomarker in the biological sample;

2. The method according to claim 1, wherein the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is a gene having a GenBank Accession Version Number selected from one or more of AK300069.1, L07414.1, BC108901.2, AJ413269.1, AH002928.2, AH009834.2, AL022146.1, X07834.1, AL390252.9 and AF054816.1.

3. The method according to claim 1, wherein the or each biomarker for predicting response to treatment with a DMARD, and/or classifying disease activity is a protein having a UniProt ID selected from one or more of Q16698, P29965, Q9UKP3, P42574, P04085, 094907, P16109, P04179, Q99523 and 09UIB8.

4. The method according to claim 1, wherein the method is a method for predicting response to treatment with a DMARD, and the or each biomarker further comprises JAM-A, PECAM-1, ANG-1, DAPP1, DCTN1 and PAI.

5. The method according to claim 4, wherein the or each biomarker for predicting response to treatment with a DMARD further comprises HB-EGF, HCLS1, IRAK4, STK4, SH2B3, NEMO, EIF4G1, PRDX1, TANK, CXCL1, SIRT2, ZBTB16, NF2, PDGF subunit B, SPRY2, HEXIM1, PPP1 R9B, PRDX3, PAR-1, ICA1, CCL17, STAMPB, MGMT, SRPK2, BML hydrolase and PLXNA4.

6. The method according to claim 5, wherein the or each biomarker for predicting response to treatment with a DMARD is a gene having a GenBank Accession Version Number selected from one or more of AK300069.1, AJ413269.1, AL022146.1, AL136649.1, AH002928.2, BC108901.2, BC107046.2, JQ287500.1, AY124380.1, X07834.1, AF178987.1 AK314352.1, X04731.1, AK222896.1, AC004634.1, AK312750.1, AY340965.1, Z93016.2, BC136451.1, AJ271718.1, AF054816.1, AK131407.1, CR407652.1, AC009299.5, BT006880.1, KF032391.1, BC029812.1, Y18000.1, M12783.1, AL354668.13, BC006460.1, KF495721.1, BT020007.1, BC002464.2, U37183.1, CH471092.1, CH471053.2, M60761.1, AC004884.1, AK312896.1, AC009365.9 and AL390252.9.

7. The method according to claim 5, wherein the or each biomarker for predicting response to treatment with a DMARD is a protein having a UniProt ID selected from one or more of Q16698, P42574, P16109, Q9Y624, P04085, Q9UKP3, Q9UBT3, P16284, Q15389, P04179, Q9UN19, Q14203, P05121, P25942, Q99075, P14317, Q9NWZ3, Q13043, Q9UQQ2, Q9Y6K9, Q9UIB8, Q04637, Q06830, Q92844, P09341, Q8IXJ6, Q05516, P35240, P01127, 043597, 094992, Q96SB3, P30048, P25116, Q05084, Q92583, 095630, P16455, P78362, Q13867, Q9HCM2 and Q99523.

8. The method according to claim 5, wherein the determining step (a) further comprises determining the quantitative level in a first set of respective biomarkers(s), wherein the quantitative level of the first set of biomarkers is greater than the quantitative level of the same respective biomarkers(s) in a normal sample, which is indicative of a non-responder, the first set of biomarkers is selected from any one or more of DECR1, CASP-3, SELP, JAM-A, PDFG subunit A, ITGB1BP2, Dkk-1, PECAM-1, ANG-1, SOD2, DAPP1, DCTN1, PAI, CD40-L, HB-EGF, HCLS1, IRAK4, STK4, SH2B3, NEMO, CD84, EIF4G1, PRDX1, TANK, CXCL1, SIRT2, ZBTB16, NF2, PDGF subunit B, SPRY2, HEXIM1, PPP1 R9B, PRDX3, PAR-1, ICA1, CCL17, STAMPB, MGMT, SRPK2, BML hydrolase, PLXNA4 and SORT1.

9. The method according to claim 5, wherein the determining step (a) further comprises determining the quantitative level in a second set of respective biomarkers(s), wherein the quantitative level of the second set of biomarkers is less than the quantitative level of the same respective biomarkers(s) in a normal sample, which is indicative of a non-responder, the second set of biomarkers is selected from any one or more of DECR1, CD40-L, ITGB1BP2, CASP-3, PDGF subunit A, Dkk-1, SELP, SOD2, SORT1 and CD84.

10. The method according claim 8, wherein the quantitative level of the or each biomarker is higher than a threshold value and is indicative of a non-responder, wherein the respective threshold value of the or each biomarker is:

11. The method according to claim 9, wherein the quantitative level of the or each biomarker is lower than a threshold value and is indicative of a non-responder, wherein the respective threshold value of the or each biomarker is:

12. The method according to any one of claims 1-3, wherein the quantitative level of the or each biomarker is between two threshold values and is indicative of high disease activity, wherein the respective threshold values of the or each biomarker are:

13. The method according to any one of claims 1-3, wherein the quantitative level of the or each biomarker is between two threshold values and is indicative of moderate disease activity, wherein the respective threshold values of the or each biomarker are:

14. The method according to any one of claims 1-3, wherein the quantitative level of the or each biomarker is between two threshold values and is indicative of low disease activity, wherein the respective threshold values of the or each biomarker are:

15. The method according to any one of claims 1-3, wherein the quantitative level of the or each biomarker is between two threshold values and is indicative of in remission, wherein the respective threshold values of the or each biomarker are:

16. The method according to any one of claims 1-15, wherein the biological sample is selected from whole blood, serum, plasma, urine, interstitial fluid, peritoneal fluid, cervical swab, tears, saliva, buccal swab, skin, brain tissue, and cerebrospinal fluid.