Introduction {#s1}
============

Recombinant interferon β (IFNβ) is one of the most widely used first line therapies in multiple sclerosis (MS), but up to 50% of patients under treatment continue to experience relapses and progression of the disease and are considered suboptimal responders.[@R1] It is, therefore, strongly desirable to count with biomarkers that allow prediction of the therapeutic response to IFNβ.

The Tumour necrosis factor-Related Apoptosis Inducing Ligand (TRAIL)/TRAIL Receptor System has been shown to be implicated in MS pathogenesis[@R2] as well as in the mechanisms of action of IFNβ.[@R5] [@R6] Expression of TRAIL and its four surface receptors in different cell types during IFNβ treatment has been studied,[@R7] [@R8] and TRAIL mRNA expression in leucocytes has been proposed as a biomarker for the therapeutic response to IFNβ in patients with MS.[@R5] Nevertheless, none of these studies has taken into account the different splicing variants of each gene, which could lead to the expression of non-functional proteins. Three different splice variants have been discovered for TRAIL: TRAIL-α, the functional isoform; TRAIL-β, which lacks exon 3; and TRAIL-γ, which lacks exon 2 and 3. The loss of exon 3 leads to the lack of the extracellular domain and, therefore, to the absence of apoptotic potential.[@R9] TRAILR-2 has two different apoptotic variants: TRAILR-2 isoform 1 (TRICK2b), which has an 87 nucleotides insertion and TRAILR-2 isoform 2 (TRICK2a), which is homologous to TRAILR-1.[@R10] Finally, there are two different splice variants for TRAILR-4: TRAILR-4 α, the long isoform, and TRAILR-4 β, which lacks exon 3.[@R11] The loss of this exon leads to the truncation of the first cysteine rich domain 1, involved in the ligand--receptor complex, which may alter its ability to bind TRAIL.

In this study, we monitored the expression kinetics of TRAIL and its receptors after induction with IFNβ. Although kinetics of TRAIL expression have previously been reported,[@R12] it remains unknown whether TRAIL receptors follow the same pattern.

The aim of the present work was to assess the effects of IFNβ treatment on the expression of the splice variants of TRAIL and its receptors in different cell subpopulations from patients with MS, and to determine whether these expressions discriminated responders from non-responders to IFNβ therapy.

Methods {#s2}
=======

Subjects {#s2a}
--------

Forty-two Spanish patients with MS prone to start IFNβ treatment were recruited from the Multiple Sclerosis Unit at Malaga Regional University Hospital (Malaga, Spain) and followed for a minimum of 2 years of treatment.

Study inclusion criteria included definite MS according to the McDonald criteria,[@R13] and treatment-naive at baseline for at least 6 months. Once they started therapy with IFNβ, patients had to be on a stable dose for at least 24 months.

Blood samples were obtained during remissions before the first IFNβ administration (baseline) and again after 12--14 months (serum and peripheral blood mononuclear cells (PBMC)) and 24 months (serum) on that therapy.

As controls (HC), twenty age-matched healthy unrelated subjects were selected.

Samples were processed following standard procedures and frozen immediately after they were received by the Málaga Regional Hospital Biobank, as part of Andalusian Public Health System Biobank. All individuals participating in the study gave their informed consent and protocols were approved by an institutional ethical committee (Comisión de Ética y de Investigación del Hospital Regional Universitario Carlos Haya).

Two patients initially taking part in the study were excluded due to the presence of high and permanent titres of neutralising antibodies (NAbs) against IFNβ. The demographic and clinical variables of the remaining 40 patients with MS and HC are listed in [table 1](#JNNP2014309932TB1){ref-type="table"}.

###### 

Demographic and clinical characteristics of the participants at baseline

  ---------------------------------------------------------------------------------------------------------------------
  Characteristics                              Healthy controls\   Responders\         Non-responders\     p Values\*
                                               n=20                n=24                n=16                
  -------------------------------------------- ------------------- ------------------- ------------------- ------------
  Female/male ratio (% women)                  14/6 (70)           17/7 (70.8)         12/4 (75)           0.867

  Age at baseline (years)                      34.0 (31.7--43.5)   34.5 (28.5--45.7)   34.0 (29.2--46.0)   0.950

  Duration of disease                                              4.0 (2.0--7.7)      5.0 (3.2--6.7)      0.774

  RR/SP ratio (% RR)                           22/2 (91.7)         15/1 (93.8)         1.000               

  EDSS at baseline                             1.0 (1.0--1.37)     1.25 (1.0--2.87)    0.101               

  Number of relapses in the 2 previous years   1.0 (1.0--2.0)      2.0 (1.0--2.0)      0.149               

  Type of IFNβ                                                                                             

   Intramuscular IFNβ-1a (%)                   4 (16.7)            1 (6.2)             0.422               

   Subcutaneous IFNβ-1a (%)                    19 (79.2)           13 (81.3)                               

   Subcutaneous IFNβ-1b (%)                    1 (4.1)             2 (12.5)                                
  ---------------------------------------------------------------------------------------------------------------------

Quantitative data are presented as median (IQR).

\*p Values: Refers to p values obtained following comparisons between healthy controls, responders and non-responders by means of a Kruskal-Wallis test (age) and χ^2^ test (gender), as well as p values obtained following comparisons between responders and non-responders by means of a Mann-Whitney U test (duration, EDSS and number of relapses) and χ^2^ or Fisher test (clinical course and type of IFNβ).

EDSS, Expanded Disability Status Scale; IFN, interferon; RR, relapsing--remitting; SP, secondary progressive; TRAIL, Tumour necrosis factor-Related Apoptosis Inducing Ligand.

Definition of response to IFNβ therapy {#s2b}
--------------------------------------

Patients were classified as optimal responders if they had neither relapses within the previous year nor progression of the disease over the same follow-up period, and as non-responders if they experienced one or more relapses or an increase of at least one point in the Expanded Disability Status Scale (EDSS) that persisted for a minimum of two consecutive visits separated by a 6-month interval.[@R14]

Detection of Neutralising Antibodies NABs {#s2c}
-----------------------------------------

The presence of NAbs in serum was checked after 12 and 24 months of treatment by the IFN-induced inhibition of virus cytopathic effect on human cells in culture, following the WHO recommendations as previously described.[@R15] The neutralisation titre of serum samples was calculated according to Kawade *et al*.[@R16] Titres ≥20 TRU/mL were considered as positive.

Determination of expression kinetics in the Jurkat T-cell leukaemia cell line {#s2d}
-----------------------------------------------------------------------------

Owing to the limited yield of PBMC obtained from each patient, gene expression kinetics were carried out in an immortalised human T lymphoblastoid cell line. Jurkat cells (Clone E.6, American Type Culture Collection) were cultured in 12-well plates, at 1.5×10^6^ cells/mL, in RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS), 1.3% Hepes, 50 μg/mL gentamicin and 2 mM L-glutamine. After incubation at 37°C and 5% CO~2~ for 24 h, cells were washed and incubated for 2 h with the former media without FBS. Cells for pre-stimulation condition (3 replicates) were collected in 1 mL of Tripure Isolation Reagent (Roche Diagnostic GmbH) and the remaining cells were stimulated with 200 UI/mL of IFNβ-1a (Avonex, Biogen Idec) and collected in the same way at 4, 8, 12 and 24 h after IFNβ induction.

Gene expression in patients with MS {#s2e}
-----------------------------------

Blood samples were taken before the first IFNβ administration and after 1 year of treatment. PBMC were isolated by a Ficoll-Hypaque density gradient. Monocytes (CD14+), CD4+ T cells and CD8+ T cells were isolated using immunomagnetic microbeads and LS MACS columns (Miltenyi Biotec GmbH), following the manufacturer\'s protocol. Purity of each cellular subset (95--97%) was assessed by flow cytometry. After isolation, the cellular subsets were cultured in 96 well plates with medium without FBS. Cells were induced with 200 UI/mL of IFNβ-1a and collected in Tripure at 4 and 24 h after induction.

### RNA isolation, primer design and quantitative reverse transcription PCR {#s2e1}

Total RNA was isolated from Jurkat cells and PBMC using a modification of the phenol--chloroform method.[@R17] Total RNA yield and quality of product was assessed with a Nanodrop 2000 Spectrophotometer.

For complementary DNA (cDNA) synthesis, 1 μg of total RNA was reverse transcribed with the M-MLV reverse transcriptase as described elsewhere.[@R18] The cDNA was stored at −80°C until use.

Primers for all amplified products are shown in [table 2](#JNNP2014309932TB2){ref-type="table"}. Primers for *TRAILR-1* and *TRAILR-3* were designed using Primer 3 software.[@R19] Primers for *TRAIL, TRAILR-2* and *TRAILR-4* were custom designed to make them complementary to the specific exon boundaries of each splice variant.

###### 

Primers used for amplification of the splice variants of TRAIL and TRAIL receptors

  Product name         Primer                 Sequence 5'→3'              Location           Product length (bp)
  -------------------- ---------------------- --------------------------- ------------------ ---------------------
  TRAIL-α              TRAIL-fwd1             TTCTACAGTTCAAGAAAAGC        Exons 3--4 joint   145
  TRAIL-rev1           GCCCAGAGCCTTTTCATTCT   Exons 4--5 joint                               
  TRAIL-β              TRAIL-fwd2             AGCTCGTTAGAAAGAAAAGC        Exons 2--4 joint   195
  TRAIL-rev2           CTCAGGAATGAATGCCCACT   Exon                                           
  TRAIL-γ              TRAIL-fwd3             ACGAGCTGAAGCAGAAAAGC        Exons 1--4 joint   145
  TRAIL-rev1           GCCCAGAGCCTTTTCATTCT   Exons 4--5 joint                               
  TRAILR-1             TRAILR1-fwd            ACCAGGAACACAGCATGTCA        Exon               132
  TRAILR1-rev          GTCACTCCAGGGCGTACAAT   Exon                                           
  TRAILR-2 isoform 1   TRAILR2-fwd            CAGTGCGAAGAAGGCACCTT        Exon 4             145
  TRAILR2-rev1         CCCCACTGTGCTTTGTACCT   87 bp insertion in exon 5                      
  TRAILR-2 isoform 2   TRAILR2-fwd            CAGTGCGAAGAAGGCACCTT        Exon 4             157
  TRAILR2-rev2         TGATGCCTGATTCTTTGTGG   Exon 5                                         
  TRAILR-3             TRAILR3-fwd            TGCACAGAGGGTGTGGATTA        Exon               108
  TRAILR3-rev          GGTCATGGTGCAGGAACTTT   Exon                                           
  TRAILR-4 α           TRAILR4-fwd1           CAATTTGCCTTCTTGCCTGC        Exon 3             164
  TRAILR4-rev          CTCTGGGACACCCTGTTCTA   Exons 4--5 joint                               
  TRAILR-4 β           TRAILR4-fwd2           GGAGTGTCCAGCAGGTCAAA        Exons 2--4 joint   137
  TRAILR4-rev          CTCTGGGACACCCTGTTCTA   Exons 4--5 joint                               
  GADPH                GADPH-fwd              GAAGATGGTGATGGGATTTC        Exon               225
  GADPH-rev            GAAGGTGAAGGTCGGAGT     Exon                                           

TRAIL, Tumour necrosis factor-Related Apoptosis Inducing Ligand.

Conventional PCR with temperatures ranging from 57°C to 63°C were performed to ensure specificity of primer pairs and to determine the optimal annealing temperature. Annealing at 57°C was optimal for TRAILR-1 and TRAIL variants, at 60°C for GAPDH and at 58°C for the remaining genes.

Quantitative PCR was performed in duplicate in a Rotor Gene Q Thermocycler (Qiagen GmbH) in a 20 μL reaction mix containing DEPC-treated water, 20 mM primers (forward and reverse), Quantitec SYBR Green PCR Master Mix (Qiagen) and cDNA. The programme consisted of a step of 15 min at 95°C, followed by 40 cycles of 95°C for 30 s, 57°C, 58°C or 60°C for 30 s and 72°C for 30 s. A melting step from 65°C to 95°C was run, increasing 0.5°C every 5 s.

The relative TRAIL and TRAIL receptor variants messenger RNA (mRNA) expression levels were calculated according to the 2^−ΔΔCT^ method, first by normalising to GAPDH and then to a calibrator sample.

Statistical analysis {#s2f}
--------------------

Comparisons of demographic characteristics at baseline between HC, responder and non-responder patients with MS were performed by Kruskal-Wallis test (age) and χ^2^ test (gender). Comparisons of clinical characteristics at baseline between responders and non-responders were performed by means of a Mann-Whitney U test (duration, EDSS and number of relapses) and χ^2^ or Fisher test (clinical course and type of IFNβ).

A Wilcoxon test was used to compare relative expression between pretreatment samples and samples from patients after 1 year of treatment with IFNβ.

A Mann-Whitney U test was used to compare relative expression of TRAIL and its receptors between responders and suboptimal responders to IFNβ therapy.

Receiver operating characteristic (ROC) analyses were performed to evaluate the predictive value of gene expression before treatment onset on therapeutic response to IFNβ.

Results {#s3}
=======

Gene expression kinetics of TRAIL and TRAIL receptors in Jurkat cells {#s3a}
---------------------------------------------------------------------

The expression kinetics of each gene after in vitro induction with IFNβ was determined in Jurkat cells ([figure 1](#JNNP2014309932F1){ref-type="fig"}). The three TRAIL variants showed a similar expression pattern: gene expression was upregulated during 1--4 h after IFNβ stimulation (peak at 4 h), and declined thereafter, followed by a second increment after 24 h, in agreement with previous data.[@R12]

![Gene expression kinetics of the Tumour necrosis factor-Related Apoptosis Inducing Ligand (TRAIL) and TRAIL receptors splicing variants in Jurkat cells on in vitro stimulation with interferon β (IFNβ). (A) TRAIL splicing variant genes are upregulated during 1-4 h of IFNβ treatment and decline thereafter. (B) Death receptor isoforms begin to upregulate at 8 h and reach a maximum at 24 h. (C) Decoy receptor isoforms begin to be induced at the 12 h time point and beyond. Expression levels are represented as relative expression compared with the reference gene GAPDH, using the ΔΔCt method. Mean and SD of each time point for each gene isoform is represented.](jnnp-2014-309932f01){#JNNP2014309932F1}

TRAILR-1 and TRAILR-2 shared the same expression pattern, characterised by a decreased expression in the first hours, and a peak 24 h after IFNβ induction. The decoy receptors also began to be induced at 8--12 h post-stimulation with IFNβ, showing a progressive increment with a marked peak at 24 h. Thus, in PBMC, the expression of TRAIL and its receptors was assessed at 4 and 24 h after in vitro stimulation with IFNβ, respectively.

Gene expression of TRAIL and TRAIL receptors in untreated patients with MS and HC {#s3b}
---------------------------------------------------------------------------------

On in vitro stimulation with IFNβ, mRNA expression of TRAIL-α was significantly increased in HC compared to untreated patients with MS in the three cellular subsets (p=0.023 in monocytes, p=0.00004 in CD4+ T cells and p= 0.021 in CD8+ T cells), as observed in [figure 2](#JNNP2014309932F2){ref-type="fig"}. We detected no other significant differences in the expression of TRAIL and receptor variants at the RNA level.

![The Tumour necrosis factor-Related Apoptosis Inducing Ligand (TRAIL) α differential expression in healthy controls (HC) and untreated patients with multiple sclerosis (MS). Expression levels are represented as relative expression compared with the reference gene GAPDH, using the ΔΔCt method. Figures show box plots: the horizontal bars are the median and lower and upper edges of the boxes represent the 25--75th centiles. Lines extending from the box are 10th and 90th centiles. The Mann-Whitney U test was used to determine statistical differences.](jnnp-2014-309932f02){#JNNP2014309932F2}

Basal gene expression of TRAIL and TRAIL receptors in patients with MS who will subsequently become responders or non-responders to IFNβ therapy {#s3c}
------------------------------------------------------------------------------------------------------------------------------------------------

Within the patient group, prior to IFNβ treatment onset, we found higher mRNA expression of TRAIL-α in monocytes, CD4+ and CD8+ T cells from patients who will subsequently be classified as non-responders than from responders (p=0.008, p=0.00006 and 0.003, respectively), as shown in [figure 3](#JNNP2014309932F3){ref-type="fig"}. Conversely, the mRNA expression of the TRAIL-β isoform was significantly higher in CD4+ and CD8+ T cells from patients who will become responders (p=0.03 and p=0.025), and showed the same trend of expression in monocytes (p=0.052). Additionally, expression of TRAILR-1 was also significantly higher in monocytes and CD4+ T cells from patients who will become responders.

![Relative expression of Tumour necrosis factor-Related Apoptosis Inducing Ligand (TRAIL) α, TRAIL-β and TRAILR-1, on in vitro stimulation with interferon β (IFNβ), before treatment onset. Those patients who will become non-responders (NR) to IFNβ therapy showed higher expression of TRAIL-α mRNA in the three cellular subsets, while those who will become responders (R) showed higher baseline expression of TRAIL-β and TRAILR-1. Expression level is represented as relative expression compared with the reference gene GAPDH, using the ΔΔCt method. Figures show box plots: the horizontal bars are the median and lower and upper edges of the boxes represent the 25--75th centiles. Lines extending from the box are 10th and 90th centiles.](jnnp-2014-309932f03){#JNNP2014309932F3}

In the search for biomarkers to predict the therapeutic response to IFNβ, ROC analyses were performed to evaluate the predictive discriminating value of TRAIL-α, TRAIL-β and TRAILR-1 expression on in vitro induction with IFNβ before treatment onset, on different cell subsets from patients with MS.

Values of areas under the curve (AUC), sensitivity, specificity, positive predictive value and negative predictive value are shown in [table 3](#JNNP2014309932TB3){ref-type="table"}. Expression of TRAIL-α in the three cellular subsets better predicted those patients who will subsequently become non-responders to treatment, with AUC \>0.7 and sensitivities greater than 80% in the three cases. Conversely, expression of TRAILR-1 on monocytes and CD4+ T cells better predicted those patients who will subsequently become responders to treatment, with AUC \>0.8, sensitivities greater than 90% and specificities greater than 75%. TRAIL-β expression did not seem a good biomarker to predict therapeutic response to IFNβ, as AUCs were under 0.7.

###### 

Sensitivity and specificity of biomarkers to predict the therapeutic response to IFNβ

  Biomarker                                                                                                                                         Cell type      AUC     Cut-off   Sensitivity %   Specificity %   PPV %   NPV %
  ------------------------------------------------------------------------------------------------------------------------------------------------- -------------- ------- --------- --------------- --------------- ------- -------
  Biomarkers that better predict suboptimal response: TRAIL-α expression on in vitro induction with IFNβ, before the onset of IFNβ therapy                                                                                   
   TRAIL-α                                                                                                                                          Monocytes      0.736   0.81      88.2            70              62.5    91.3
  CD4+ T cells                                                                                                                                      0.794          0.12    83.3      65.2            55.5            88.2    
  CD8+ T cells                                                                                                                                      0.733          0.055   83.3      61.8            60.6            87.5    
  Biomarkers that better predict good response: TRAIL-β and TRAILR-1 expression on in vitro induction with IFNβ, before the onset of IFNβ therapy                                                                            
                                                                                                                                                                                                                             
   TRAIL-β                                                                                                                                          CD4+ T cells   0.662   0.095     77.3            47.8            73.9    52.4
  CD8+ T cell                                                                                                                                       0.683          0.115   71.1      52.6            75              47.6    
   TRAILR-1                                                                                                                                         Monocytes      0.803   0.67      90.5            78.6            79.2    81.8
  CD4+ T cells                                                                                                                                      0.883          0.29    90        75              85.7            81.8    

AUC, area under the curve; IFN, interferon; NPV, negative predictive value; PPV, positive predictive value; TRAIL, Tumour necrosis factor-Related Apoptosis Inducing Ligand.

Effects of long-term IFNβ therapy on TRAIL and TRAIL receptors expression in PBMC from patients with MS {#s3d}
-------------------------------------------------------------------------------------------------------

Relative expression of TRAIL and TRAIL receptors was significantly modified along the first year of therapy only for two genes. Expression of the isoform 2 of the death receptor TRAILR-2 was significantly reduced in monocytes from first-year responder as well as non-responder patients (p=0.004 and 0.016, respectively). Expression of TRAIL-α was significantly increased in both T cell subsets but exclusively in patients who showed a clinical response to systemic IFNβ treatment (p=0.004 in CD4+ cells and p=0.002 in CD8+ cells from responders). Induction of TRAIL-α mRNA expression in monocytes by IFNβ therapy was achieved in 13 of 24 responders and, therefore, overall expression in this subset was not significantly increased ([figure 4](#JNNP2014309932F4){ref-type="fig"}).

![Effects of long-term treatment with interferon β (IFNβ) on Tumour necrosis factor-Related Apoptosis Inducing Ligand TRAILR-2 isoform 2 and TRAIL-α expression. Expression levels are represented as relative expression compared with the reference gene GAPDH, using the ΔΔCt method. Expression before and after 1 year of IFNβ therapy in the same individual patients is plotted. (A) Long-term therapy with IFNβ significantly decreased mRNA expression of TRAILR-2 isoform 2 in monocytes from patients who will become responders (R) as well as those who will become non-responders (NR) to treatment. (B) TRAIL-α mRNA expression in CD4+ and CD8+ T cells significantly increased along with IFNβ therapy, but only in responders.](jnnp-2014-309932f04){#JNNP2014309932F4}

Discussion {#s4}
==========

TRAIL expression in PBMC has been proposed as a biomarker for therapeutic response to IFNβ in MS,[@R5] but different splice variants were not taken into account. Several diseases have been related to alternative splice variant or aberrant splicing.[@R20] We hypothesise that alternative splice variants could be implicated in response to treatment and, thus, they should be measured independently.

In our study, we first determined the expression kinetics of TRAIL and its receptors' isoforms after in vitro stimulation with IFNβ. In agreement with the existing literature, TRAIL isoforms reached a peak 4 h after IFNβ stimulation.[@R12] Conversely, TRAIL receptors reached maximum expression at 24 h after induction. To the best of our knowledge, our group has been the first to assess this issue. It remains unknown whether in vivo gene expression of TRAIL receptors after induction with IFNβ shares the same pattern as in vitro, or whether they go through a more complex regulation. We are currently performing a longitudinal study, using blood samples drawn at different time points after injection with IFNβ, to address this question.

Contrary to previous studies,[@R12] [@R21] [@R22] we found that the expression of TRAIL-α was significantly lower in monocytes and T cells from untreated patients with MS than from controls. In agreement, Kurne *et al* found lower soluble TRAIL concentrations in patients with MS than in HC.[@R23] One possible explanation for this finding may be the lower levels of endogenous IFNβ reported in patients with MS,[@R24] which are unable to induce TRAIL synthesis to the same extent as they can in healthy subjects.

Cell type-specific induction of TRAIL by IFNβ in patients with MS has been reported. We have shown upregulation of TRAIL-α expression after IFNβ therapy in monocytes, and in CD4+ and CD8+ T lymphocytes. Earlier works found increased levels of TRAIL after IFNβ-1a injection only on the surfaces of monocytes and/or granulocytes.[@R6] [@R8] [@R25] However, upregulation of TRAIL expression in lymphocytes[@R7] [@R26] [@R27] and T cell lines[@R28] from patients with MS and HC has also been reported. These differences could be attributed to the different culture conditions, kinetics and activation stimulus of the cells.

Within the patient group, those patients with MS who will be subsequently considered as responders to IFNβ therapy could be clearly distinguished from non-responders by baseline expression of TRAIL-α and TRAILR-1, indicating underlying differences in transcriptional regulation of these genes. The in vitro addition of IFNβ to previously IFNβ naïve cell subsets triggered a higher expression of TRAIL-α and a lower expression of TRAIL-β in monocytes and T cells from patients who will continue to experience relapses or disease progression while treated with IFNβ, and a higher expression of TRAILR-1 in monocytes and CD4+ T cells from those who will have neither relapses nor progression of the disease. Wandinger *et al*[@R5] reported no differences in baseline expression of TRAIL regarding IFNβ responsiveness, but found higher baseline concentrations of soluble TRAIL in responders, findings that could not be replicated by others.[@R29] Conversely, another study found higher levels of TRAIL in monocytes from relapsing--remitting quiescent disease patients (responders).[@R8] Comparison among studies regarding biomarkers for therapeutic response to IFNβ in MS is difficult, due to the unique design and experimental approach of each single study.

It is known that TRAIL is upregulated on activation of monocytes[@R24] [@R30] and T cells,[@R28] [@R29] while TRAILR-1 and TRAILR-2, in turn, are downregulated on the RNA level when T cells are activated.[@R28] In our ROC analysis, baseline TRAIL-α expression in the three cellular subsets at the cut-off levels described in [table 3](#JNNP2014309932TB3){ref-type="table"} predicts suboptimal response with a sensitivity greater than 80% and a specificity around 65%, while baseline TRAILR-1 expression over 0.803 in monocytes and over 0.883 in CD4+ T cells predicts a good response to IFNβ therapy with a sensitivity higher than 90% and a specificity over 75%, promoting the idea of using the expression of TRAIL-α and TRAILR-1 as biomarkers to predict IFNβ treatment response. Probably, the changes in baseline expression of TRAIL-α and TRAILR-1 of such a low magnitude in patients that will become non-responders may not be able to induce the expected changes after a significant increase in TRAIL, such as a strong inhibition of activation of autoreactive T cells, promotion of Tregs, or induction of neurons and oligodendrocytes apoptosis after lymphocyte infiltration.[@R31] They may be more probably reflecting a compensatory mechanism to prevent, in part, the activation of autoreactive T cells in those untreated patients with an increased state of activation who will continue to experience relapses or progression. Another explanation for this increase in TRAIL-α expression at baseline in non-responders may be attributed to a comparatively higher pretreatment serum concentration of endogenous IFNβ in non-responders than responders, as shown by Axtell *et al*.[@R32]

Long-term treatment with IFNβ has been reported to increase TRAIL mRNA expression.[@R5] [@R7] [@R8] [@R21] In our study, IFNβ therapy significantly increases TRAIL-α expression in both T cell subsets, but exclusively in responder patients, similar to results reported by others,[@R5] [@R7] [@R8] reaching levels of expression that overcome those found in HC. In this sense, upregulation of TRAIL in responders contributes to the deletion of autoreactive lymphocytes, the inhibition of activation of autoreactive lymphocytes and limitation of the production of proinflammatory cytokines, mechanisms that may contribute to the clinical benefit of IFNβ therapy. Moreover, TRAIL-α promotes the induction of endogenous IFNβ,[@R33] amplifying its effects. Although TRAIL-α is also known to induce apoptosis of oligodendrocytes and neurons, as IFNβ enhances the integrity of the blood--brain barrier (BBB) in responders, the upregulation of TRAIL induced by IFNβ therapy probably comprises immunoregulatory mechanisms outside the central nervous system (CNS) more than apoptosis of brain cells.

However, although non-responders had higher basal TRAIL-α levels than responders, they showed an inability to increase TRAIL-α expression in T cells in response to IFNβ therapy, that is, these patients are already responding to endogenous IFNβ, probably at maximum level, but still show a disease activity that obtains no benefit from further administration of IFNβ. This inability to increase TRAIL-α expression could, on one hand, lead to a failure in the apoptosis of auto-aggressive T cells invading the CNS and, on the other hand, prevent inhibition of the activation of autoreactive T cells contributing to the perpetuation of the inflammatory response and the appearance of relapses and CNS damage.

Moreover, we found that long-term IFNβ therapy decreased the expression of TRAILR-2 isoform 2 in monocytes from patients with MS, regardless of the therapeutic responsiveness to this drug. Expression of TRAILR-2 is constitutive in peripheral blood leucocytes and may be involved in the apoptosis of activated lymphocytes[@R22] [@R34] and monocytes. In later stages of tissue pathology, macrophages outnumber T cells in CNS infiltrates and are needed to execute myelin phagocytosis in the demyelinating process. In this sense, downregulation of TRAILR-2 isoform 2 in monocytes could be one of the mechanisms responsible for the clinical benefits of IFNβ, as apoptosis of peripheral monocytes will be prevented, allowing these cells to differentiate into macrophages, cross the BBB and exert their functions.

Conclusions {#s5}
===========

The present study shows that long-term therapy with IFNβ increased the expression of TRAIL-α in both T cell subsets exclusively from responder patients and decreased the expression of the isoform 2 of the death receptor TRAILR-2 in monocytes from responders as well as non-responders.

ROC analysis revealed that expression of TRAIL-α in the three cellular subsets, before the onset of IFNβ therapy, better predicted those patients who will subsequently become non-responders to treatment while expression of TRAILR-1 on monocytes and CD4+ T cells better predicted those patients who will become responders.

Further studies are necessary in order to elucidate the pathways implicated in this process, which could reveal the precise mechanism by which IFNβ is exerting its beneficial effects on patients with MS, but TRAIL-α and TRAILR-1 should be assessed as predictors in future IFNβ responsiveness studies. However, inconsistent results among different studies regarding therapeutic biomarkers to IFNβ in MS make us evaluate our results with caution. To prove the effectiveness of these predictive markers, studies with larger sample sizes are warranted.

The authors thank all the participants of this study. They also thank the "Red Española de Esclerosis Múltiple REEM (RD12/0032/0006)" and Fundación Española de Esclerosis Múltiple (FEDEM).

**Contributors:** CL-G, BO-M, M-JP-M and MS performed the cellular subsets isolation, the gene expression analysis and the Nabs test. VR-G, PU and OF recruited the patients, evaluated the response to therapy and participated in the analyses. CL-G and LL designed the study, evaluated the results and wrote the manuscript. OF revised it critically for important intellectual content.

**Funding:** The authors acknowledge the support from Fondo de Investigación Sanitaria & Fondo Europeo de Desarrollo Regional (PS09/01764) and Consejería de Salud de la Junta de Andalucía (SAS07/0231) to LL. The authors would like to thank Biogen Idec Iberia SL for financial support.

**Competing interests:** OF has received honoraria as consultant in advisory boards and as chairmen or lecturer in meetings, and has also participated in clinical trials and other research projects promoted by Biogen-Idec, Bayer-Schering, Merck-Serono, Teva, Novartis, Almirall and Allergan.

**Patient consent:** Obtained.

**Ethics approval:** Comisión de Ética y de Investigación del Hospital Regional Universitario Carlos Haya.

**Provenance and peer review:** Not commissioned; externally peer reviewed.

[^1]: OF and LL contributed equally to this study.
