The invention relates to a method for the prognosis of multiple sclerosis and, more specifically, to a method for predicting the clinical progress of patients diagnosed with multiple sclerosis by means of analyzing the expression levels of a series of genes.
Multiple sclerosis (MS) has a prevalence of around 70 cases every 100,000 inhabitants in Spain and in Western civilization it is the most common cause of chronic neurological disability in young adults after traffic accidents. Approximately 70% of cases starts between 20 and 40 years of age, with a peak in the age of onset around 25 or 30 years, so the huge impact it has on the professional, family and social life of those affected, as well as the enormous economic and social expense it generates, which is similar to that of Alzheimer's disease, is easily understandable.
Multiple sclerosis is a heterogeneous disease in its presentation and progress in which 80-85% of the patients present a clinical course which progresses to auto-limited flare-ups which, as they are repeated, cause a functional residual deficit (relapsing-remitting form; RR-MS). After 10-15 years of progress, 50% of them will pass to a secondary progressive course of increment of the disability not related to the flare-ups (secondary progressive form; SP-MS), and after 25 years the percentage reaches 90% of the patients. In 10% of the cases the course is progressive from the onset (primary progressive form; PP-MS). 10 to 20% of the patients will remain without significant sequelae 15 years after the onset of the disease (benign forms) and in 1-3% of the cases, however, the patients will progress accumulating a great disability in a few months after the start of the disease (aggressive or fulminant forms).
Interferon beta (Betaferon®, Rebif®, Avonex®) and glatiramer acetate or copolymer (Copaxone®) are the first medicines that have shown beneficial effects in the relapsing-remitting form of multiple sclerosis. These medicines reduce the formation of plaques and the number of flare-ups by one third compared with patients without treatment. Treatment with Natalizumab (Tysabri®) has recently been included in the therapy of multiple sclerosis, having greater effectiveness than prior treatments (they prevent flare-ups by approximately 60%), although with potential serious side effects. However, the individual response to treatment is unpredictable and ranges from excellent to complete ineffectiveness. The lack of biological tests which predict the activity and aggressiveness of the disease prevents prescribing the best treatment to each patient and forces administering a preventive treatment for life with the subsequent economic cost and the effect on the quality of life. Being able to have a predictive test of the aggressiveness and activity of the disease would allow a customized treatment.
Standard methods of diagnosis of multiple sclerosis include determining levels of IgG in CSF, brain imaging by means of magnetic resonance and spinal cord imaging and the exclusion of other autoimmune diseases by means of serum determinations. Although the usefulness of said tests in predicting the course of multiple sclerosis has recently been studied, their predictive capacity is very limited and in clinical practice they are used for diagnostic purposes but not for prognostic purposes or for deciding on or monitoring therapy. Therefore, reaching a diagnosis and a suitable prognosis continues to be a problem.
Different studies of the natural history of multiple sclerosis have allowed identifying some clinical variables associated with the progress of multiple sclerosis. The factors which best predict a relatively benign course are belonging to the female sex, the onset of the disease at an early age (less than 30 years), uncommon attacks, a relapsing-remitting pattern and a mild nature of the disease in the studies by means of magnetic resonance of the central nervous system. In contrast, the factors which predict a more aggressive course are the male sex and late onset (because it is associated with the progressive forms), the recurrence of the second flare-up in the first year after the first flare-up, accumulating disability early on, the clinical onset with motor or coordination symptoms and the persistence of sequela after the first flare-up, and especially reaching certain disability levels (Kurtzke Expanded Disability Status Scale EDSS: 3.0, 4.0, 6.0) at early ages.
Until now, the attempt has been made to develop methods of diagnosis and prognosis based on the detection of autoantibodies in serum (Bielekova, B. and, Martin R. Brain. 2004 July; 127(Pt 7):1463-78.; Berger, T. and Reindl, M., 2006, Disease Markers, 22:207-212). For example, Berger, T. et al. (New England J. Medicine, 2003, 146:181-197) have described that the presence of anti-myelin antibodies (anti-MOG) are capable of predicting the risk of first relapse in patients suffering a clinically isolated syndrome, suggesting multiple sclerosis. However, subsequent well-designed studies have not been able to confirm their predictive usefulness (Kuhle J. et al., N Engl J Med. 2007, 356:371-8.; Pelayo R. et al. N. Engl. J. Med. 2007, 356:426-8.)
However, until now no antibody has been identified which meets the requirements of a diagnostic or prognostic biomarker of multiple sclerosis. Nor is the simultaneous determination of several autoantibodies of use, and it furthermore involves great difficulties and a high cost.
In addition, gene expression studies (transcriptome) by means of DNA chips allow having a global vision of the genes that are participating in a process, so this type of analysis could become a valuable clinical tool for the diagnosis or prognosis of MS.
Until now differences have been described in the expression levels of various genes in multiple sclerosis, which could be candidates as biomarkers of multiple sclerosis (reviewed in Goertsches, R. et al., 2006, Current Pharmaceutical Design, 12:3761-3779). These studies compared the expression patterns between patients with multiple sclerosis and controls, but the association between said patterns and the progressive course and prognosis of the disease were never specifically studied.
Ramanathan et al (J. Neuroimmunol., 2001, 116:213-219) have described 34 genes differentially expressed in RR-MS patients in comparison with healthy subjects, most of them related to inflammatory and immune processes.
Bomprezzi et al. (Human Molecular Genetics, 2003, 12:2191-2199) identified a series of genes the expression levels of which in PBMCs allows distinguishing RR-MS and SP-MS patients and healthy volunteers. By means of this approach, over a thousand genes which allowed distinguishing samples of multiple sclerosis from controls could be identified. The strongly dominant genes included HSP70 and the CDC28 protein kinase (CKS) 2 which, combined with histone HI of the (HIF) 2 family and the PAFAH1B1, respectively, allowed a good discrimination between multiple sclerosis and controls. These pairs had a prediction value of 80% for classifying an independent sample in the right class. A correlation between the most discriminatory pairs of genes and relevant biological pathways of multiple sclerosis was also observed. Such molecules, which were highly expressed in multiple sclerosis included CD27, TNF receptor, the alpha locus of the T cell receptor and its associated chain ZAP70, and the zinc finger protein (ZNF) 148. Furthermore, the IL-7 receptor (IL-7R), which is required for the development of T and B cells, was also strongly overexpressed. The repressed genes in multiple sclerosis were HSP70 and CKS2, both involved in apoptosis regulation. It has previously been suggested that HSP70 can be an autoantigen in multiple sclerosis, but it can also be involved in the degradation of mRNA in the ubiquitin-proteasome pathway. The activation of the remodeling process of the extracellular matrix was evident due to the overexpression of the matrix metalloproteinase (MMP)-19 and the repression of a TIMP1 inhibitor.
The expression patterns for multiple sclerosis and the pathophysiology of the disease have been analyzed in several studies. Iglesias et al. (J. Neuroimmunol., 2004 150:163-77) identified a system of 553 genes differentially expressed in RR-MS compared with the healthy controls, 87 of which were highly significant. Among the genes differentially expressed, some involved in the activation and co-stimulation of T cells could be identified, which included several interferon response genes, such as IL-12, CD40, cytotoxic antigen 4 (CTLA4), T cell receptors, immunoglobulins, the IL-6 receptor, the IL-8 receptor, and integrins, for example VLA4 and VLA6, as well as different genes of the E2F pathway (E2F2, E2F3, CDC25A, CDK2), the thymopoietin (TMPO), and PRIM1. The importance of the E2F pathway in multiple sclerosis was validated in experimental autoimmune encephalitis (EAE). E2F1-deficient mice showed only a mild course of the EAE disease.
The gene patterns for the activity of multiple sclerosis have also been studied. International patent application WO03081201 identified a pattern of 1109 genes in PBMCs from 26 multiple sclerosis patients compared with healthy volunteers regardless of the state of the activation of the disease. The pattern was validated with the LOOCV method, which gave only two errors in the classification, proving that the patterns observed represent a true biological phenomenon. These genes included those related to activation and expansion, inflammatory stimuli of the T cell (cytokines and integrins), spreading epitope, and apoptosis. Comparison of the profiles of the expression in PBMCs of multiple sclerosis patients in a flare-up showed a pattern of 721 genes. Protease L, CTSLI and the MCP1 and MCP2 proteins were overexpressed during the flare-up. In contrast, several genes related to apoptosis such as cyclin G1 and the caspases (CASP) 2, 8 and 10 were repressed.
WO03023056 describes methods for the diagnosis of and/or the susceptibility to multiple sclerosis by means of determining variations in the expression levels of 25 genes.
Individually, a gene (CX3CR1) which has been identified in expression analysis in sub-populations of T cells has been proposed as a marker for the activity of the disease. CX3CR1 is repressed in RR-MS and PP-MS patients compared with healthy volunteers. This finding has been validated by real-time PCR and by flow cytometry in independent cohorts. The NK cells are responsible for the phenotype, whereas the expression of CX3CR1 is not altered in CD8 cytotoxic cells in multiple sclerosis patients with respect to healthy controls.
US2004/0091915 describes a method for predicting the survival rate of patients diagnosed with multiple sclerosis by means of detecting a deletion in the CCR5 gene.
WO2005054810 describes a method for predicting the survival rate of patients diagnosed with multiple sclerosis by means of detecting a deletion in the gene CD24.
In WO03001212 describes a method for the diagnosis of multiple sclerosis based on detecting in a sample isolated from the subject the absence of the wt-SARG-1 protein or of the mRNA which encodes it.
US20050064483 describes a method for monitoring the response of a multiple sclerosis patient to treatment with interferon-beta or with glatiramer acetate by means of detecting variations in the expression of at least 4 genes selected from a group of 34 genes.
US20050089919 describes a method for detecting multiple sclerosis which comprises detecting variations in the expression of at least one gene selected from a series of 31 genes.
US20050164253 describes a method for detecting multiple sclerosis which comprises detecting variations in the expression of at least one gene selected from the group of RIPK2, NFKBIE, TNFAIP3, DAXX, TNFSF10, BAG1, TOP1, ADPRT, CREB1, MYC, BAG4, RBBP4, GZMA, BCL2 and E2F5.
US20060115826 describes a method for the diagnosis of multiple sclerosis which comprises detecting variations in the expression of at least two genes selected from a set of 107 genes associated with inflammatory processes.
WO02079218 describes a method for the diagnosis of multiple sclerosis which comprises detecting variations in the expression of a selected gene panel in that they show variations in their expression level in an experimental animal model of autoimmune encephalitis. In this study, the analysis of the expression of the different genes was conducted by means of a human DNA chip in which about 14000 different genes were represented.
WO03081201 describes a method for the diagnosis of multiple sclerosis based on detecting variations in the expression of a gene panel represented in a human DNA chip which contained 12625 human genes.
WO03095618 describes methods for the diagnosis of multiple sclerosis, for the differential diagnosis of multiple sclerosis with respect to lateral amyotrophic sclerosis, for predicting the response of a subject diagnosed with multiple sclerosis to a treatment with Avonex by means of detecting variations of the expression of a series of genes involved in different signaling pathways.
However, all the methods described until now have been aimed at detecting differences between patients suspected of presenting multiple sclerosis and control subjects, whereby they have an essentially diagnostic use, but they do not allow predicting the progress of the disease in patients who have already been diagnosed with multiple sclerosis. Therefore, there is a need for methods which allow predicting the progress of the disease in patients already diagnosed.
In a first aspect, the invention relates to an in vitro method for determining the clinical prognosis of a patient who has multiple sclerosis which comprises
In another aspect, the invention relates to an in vitro method for determining the clinical prognosis of a patient who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for diagnosing multiple sclerosis in a subject which comprises
In another aspect, the invention relates to a kit comprising a set of probes wherein said set comprises a probe specific for each of the genes indicated in at least one table selected from the group of Tables 3, 5-8 and 11.
In another aspect, the invention relates to the use of a kit of the invention for determining the prognosis of a patient diagnosed with multiple sclerosis, for determining the effectiveness of a treatment for multiple sclerosis or for diagnosing multiple sclerosis in a patient.
The authors of the present invention, using DNA microarrays, have identified a series of genes which are differentially expressed in patients diagnosed with multiple sclerosis in which the disease has a good prognosis with respect to patients in which the disease shows a bad prognosis or to control subjects. Likewise, the authors of the invention have identified a series of genes the expression of which is modified in patients diagnosed with multiple sclerosis in which the disease has a bad prognosis. From an initial set of identified genes, a subset of 13 genes was validated by means of real-time PCR the expression variations of which allowed predicting the type of prognosis of the patients in a significant manner (p<0.05).
Thus, in a first aspect, the invention relates to an in vitro method (hereinafter the first method of the invention) for determining the clinical prognosis of a patient who has multiple sclerosis which comprises
Additionally, from the expression levels of these 13 genes and by using two clinical variables (ESSS and MSFC), a classifier was obtained which allowed predicting the progress of the disease with a precision of 95%. Thus, in another aspect, the invention relates to an in vitro method (hereinafter the first method of the invention) for determining the clinical prognosis of a patient who has sclerosis which comprises
According to the present invention, “determining the clinical prognosis” is understood as giving an opinion as to the future condition of the patient (clinical (physical and cognitive) disability) after a determined number of years (e.g. 2, 5, 10 years from the moment of the opinion). The clinical prognosis can be performed in recently diagnosed patients or after the first flare-up, as well as at any time of the course of his disease. The condition of the patient can be defined based on the symptoms of multiple sclerosis, including a reduced capacity for controlling small movements, reduced attention span, reduced coordination, reduced discerning capacity, reduced memory, depression, difficulty in speaking or understanding language, dizziness, double vision, eye discomfort, facial pain, fatigue, loss of balance, problems with movement that are slowly progressive and begin in the legs, muscle atrophy, muscle spasms (especially in the legs), muscle spasticity (uncontrollable spasm of muscle groups), numbing or abnormal sensation in any area, pain in the arms and legs, paralysis of one or both arms or legs, bad pronunciation, tingling sensation, shaking in one or both arms or legs, uncontrollable and rapid eye movements, increased urinary frequency, difficulty in urinating, urinary urgency, urinary incontinence, vertigo, loss of vision, walk/gait anomalies and weakness in one or both arms or legs.
“Modal value” is understood in the context of the present invention as a value of the variable (in this case, of the expression levels) which partitions the range of values of said variable into two or more sub-ranges. Suitable methods for determining said value have been described in Dougherty, J. et al., (Proc. of the 12th International Conference on Machine Learning; 1995. p. 194-202) and by Liu H. et al. (Data Mining and Knowledge Discovery, 2002, 6:393-423), the content of which is incorporated in the present application in its entirety. In a preferred embodiment, in order to obtain said modal value the variable is discretized by means of a supervised learning algorithm (computational) of the more informative discretization thresholds with respect to a reference variable (in the present invention, the diagnosis). To calculate the discretization ranges, if starting from the ordered sequence of values Ai={v1, v2, . . . , vn}, the information gain is evaluated with respect to a reference variable for all the possible partitions (n−1). The partition with the greatest information gain is the one that is used for comparing with the remaining attributes. The decisions of the nodes will be [A[x]<vi] and [A[x]≧vi].
Possible discretization algorithms suitable for their use in the present invention include the decision tree, the equal frequency algorithm and the equal distance algorithm. In a preferred embodiment, the discretization algorithm is the decision tree.
“Tables of conditional probabilities” is understood in the context of the present invention as a table in which the possible modal values of the expression of a determined gene or clinical variable are represented, and in which each of the modal values is correlated with a determined probability that the disease of the patient will follow a positive or negative prognosis. In a preferred embodiment, the tables of conditional probabilities between the modal values of the expression of each of the genes and the probability values that the multiple sclerosis has a good or bad prognosis and between the modal values of each of the clinical variables and the probability values that the multiple sclerosis has a good or bad prognosis are those indicated in Table 14.
“KLHDC5 gene” is understood as the gene encoding the kelch domain containing 5 protein the human variant of which is described in the GenEMBL database under accession number BC 108669.
“CASP2 gene” is understood as the gene encoding the precursor of caspase 2. The human form of said gene is described in the GenEMBL database under accession numbers U 13021 and U13022.
“EMID1 gene” is understood as the gene encoding the precursor of the EMI domain-containing protein 1. The human form of said gene is described in the GenEMBL database under accession number AJ416090.
“PRO1073 gene” is understood as the gene described in the GenEMBL database under accession number AF113016.
“BTBD7 gene” is understood as the gene encoding the BTB/POZ domain-containing protein 7. The human form of said gene is described in the GenEMBL database under accession number BX538231.
“MGC25181 gene” is understood as the gene encoding the hypothetical MGC25181 protein. The human form of said gene is described in the GenEMBL database under accession number AC114730.
“WDR20bis gene” is understood as the gene encoding the WD repeat-containing protein 20. The human form of said gene is described in the GenEMBL database under accession number BCO28387.
“NEK4 gene” is understood as the gene encoding the serine/threonine kinase Nek4. The human form of said gene is described in the GenEMBL database under accession number L20321.
“SYLT2 gene” is understood as the gene encoding the Synaptotagmin-like protein 2. The human form of said gene is described in the GenEMBL database under accession number AK000170.
“DOCK10 gene” is understood as the gene encoding the dedicator of cytokinesis protein 10. The human form of said gene is described in the GenEMBL database under accession number AB014594.
“TTC 10 gene” is understood as the gene encoding the tetratricopeptide repeat protein 10. The human form of said gene is described in the GenEMBL database under accession number U20362.
“PTPRC gene” is understood as the gene encoding the protein tyrosine phosphatase receptor type C. The human form of said gene is described in the GenEMBL database under accession number BC017863.
“CTLA4 gene” is understood as the gene encoding the precursor of cytotoxic T-lymphocyte protein 4. The human form of said gene is described in the GenEMBL database under accession number AF414120.
“EDSS” is understood as the Kurtzke Expanded Disability Status Scale, as it is defined in Kurtzke, J. F. (Neurology, 1983, 33:1444-1452).
“MSFC” is understood as the Multiple Sclerosis Functional Composite, as is defined in Fischer, J. S. et al. (National MS Society Clinical Prognoses Assessment Task Force. Mult. Scler. 1999, 5:244-250).
The determination of the expression values of a nucleic acid is performed by means of the relative measurement of the expression levels of a gene of interest compared to the expression levels of a reference nucleic acid. Said measurements can be carried out by any method known by the person skilled in the art, such as those included in Sambrook, J. et al. (Molecular Cloning: A Laboratory Manual. 2nd ed., Cold Spring Harbor Laboratory, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y. (1989)) and Ausubel et al. (Current Protocols in Molecular Biology, eds. Ausubel et al, John Wiley & Sons (1992)). Typical processes for detecting the polynucleotide resulting from the transcription of a gene of interest include the extraction of RNA from a cell or tissue sample, hybridization of said sample with a labeled probe, i.e., with a nucleic acid fragment having a sequence complementary to the molecule of nucleic acid to be detected, and detection of the probe (for example, by means of Northern blotting). The invention also contemplates the detection of the expression levels of a determined gene by means of using primers in a polymerase chain reaction (PCR), such as anchor PCR, RACE PCR, ligase chain reaction (LCR). In a preferred embodiment, the determination of the modal values of the expression is carried out by means of real-time PCR.
These methods include the steps of collecting a cell sample from a subject, isolating the mRNA from said samples, converting the mRNA present in the sample into complementary DNA (cDNA), contacting the cDNA preparation with one or several primers which specifically hybridize with the target gene in suitable conditions for the hybridization and amplification of said nucleic acid followed by the detection of the presence of an amplification product. Alternative amplification methods include self-sustained sequence replication (Guatelli, J C. et al., (1990) Proc. Natl. Acad. Sci. USA 87:1874-1878), transcriptional amplification system (Kwoh, D. E. et al., (1989) Proc. Natl. Acad. Sci. USA 86:1173-1177), Q-Beta replicase (Lizardi, P. M. et al. (1988) BioTechnology 6:1197) or any other known nucleic acid amplification method, followed by the detection of the amplified molecules using techniques that are well known by the person skilled in the art. These methods of detection are particularly suitable for detecting nucleic acids when said molecules are present in a very reduced number of copies.
In other embodiments, the genes per se can be used as markers of multiple sclerosis. For example, the increase of the expression of a determined gene can be due to the duplication of the corresponding gene, such that the duplication can be used as a diagnosis of the disease. The detection of the number of copies of a target gene can be carried out using methods that are well known by the person skilled in the art. The determination of the number of copies of a determined gene is typically carried out by means of Southern blot in which the complete DNA of a cell or of a tissue sample is extracted, hybridized with a labeled probe and said probe is detected. The labeling of the probe can be by means of a fluorescent compound, by means of an enzyme or an enzymatic cofactor. Other typical methods for the detection/quantification of DNA include direct sequencing, column chromatography and quantitative PCR using standard protocols.
The determination of the expression levels of a gene can be carried out in any biological sample from a subject, including different types of biological fluids, such as blood, serum, plasma, cerebrospinal fluid, peritoneal fluid, feces, urine and saliva, as well as tissue samples. The biological fluid samples can be obtained by any conventional method as can the tissue samples; by way of illustration said tissue samples can be biopsy samples obtained by surgical resection.
The second method of the invention contemplates the simultaneous determination of the expression values of a larger number of genes. Thus, the second method of the invention can include the determination of the expression values of at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10 and at least 11 genes.
In a preferred embodiment, the second method of the invention requires the determination of the expression values of the KLHDC5 gene and of the EDSS clinical variable. In another preferred embodiment, the second method of the invention requires additionally determining the expression value of the CASP2 gene. In another preferred embodiment, the expression value of the EMID1 gene is additionally determined. In an even more preferred embodiment, the value of the MSFC clinical variable is additionally determined. In an even more preferred embodiment, the method additionally comprises determining the expression value of the PRO1073 gene. In another preferred embodiment, the second method of the invention includes additionally determining the expression value of the BTBD7 gene. In an even more preferred embodiment, the method involves additionally determining the expression value of the MGC2518 gene. In another embodiment, the method of the invention involves additionally determining the expression value of the WDR20bis gene. In an even more preferred embodiment, the method of the invention involves additionally determining the expression value of the NEK4 gene. In another embodiment, the method of the invention involves additionally determining the expression value of the DOCK10 gene. In another preferred embodiment, the method of the invention involves additionally determining the expression value of the TTC10 gene. In an even more preferred embodiment, the method of the invention involves additionally determining the expression value of the PTPRC gene. In another preferred embodiment, the method of the invention involves additionally determining the expression value of the CTLA4 gene.
The inventors have additionally shown the existence of different genes which are differentially expressed in patients diagnosed with multiple sclerosis with a bad prognosis with respect to patients diagnosed with multiple sclerosis with a good prognosis and with respect to control subjects, which allows the development of prognostic methods for predicting the development of the disease. The authors of the present invention have additionally shown the existence of genes which are differentially expressed in subjects diagnosed with multiple sclerosis with respect to healthy patients, which allows the use of said genes for diagnostic purposes.
Thus, in another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for identifying the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for determining the clinical prognosis of a subject who has multiple sclerosis, for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a customized therapy to a subject who has sclerosis which comprises
In another aspect, the invention relates to a method for diagnosing multiple sclerosis in a subject which comprises
The definition of “determination of the clinical prognosis” has been described above.
“Monitoring the effect of the therapy administered to a subject who has multiple sclerosis” is understood according to the present invention as determining if a therapy has any incidence on the prognosis.
“Assigning a customized therapy to a subject who has multiple sclerosis” is understood as deciding, based on the prognosis of an individual, on the most suitable type of therapy for preventing the occurrence of the previously indicated symptoms. In cases of worse prognosis, a more aggressive therapy is applied from the time that said worse prognosis is detected. Thus, more aggressive therapies include immune modulators to aid in controlling the immune system, including interferons (Avonex, Betaseron or Rebif), monoclonal antibodies (Tysabri) and glatiramer acetate (Copaxone) and chemotherapy.
“Reference value” is understood as a measurement of the expression of a determined gene or polypeptide that can be calculated or established from one or several control samples. These can come from a healthy subject, from a subject with multiple sclerosis, or from subjects with a good or a bad prognosis, according to the objective of the method.
The person skilled in the art will understand that the determination of the expression levels of the genes included in Tables 3, 5, 6, 7 and 8 can be carried out using techniques known by the person skilled in the art.
The determination of the expression levels of a nucleic acid relating to the levels of a reference nucleic acid can be carried out by any method known by the person skilled in the art, as has been described above.
In other embodiments, the genes per se can be used as markers of multiple sclerosis in those cases in which the increase of the expression of a determined gene can be due to the duplication of the corresponding gene, such that the duplication can be used as a diagnosis of the disease. The detection of the number of copies of a target gene can be carried out using the methods described above.
Alternatively, the invention contemplates methods for determining the clinical prognosis of a subject who has multiple sclerosis or for monitoring the effect of the therapy administered to a subject who has multiple sclerosis or for assigning a therapy to a subject who has multiple sclerosis in which the expression level of one or several proteins encoded by the genes which are indicated in Tables 1 to 4 is determined. In this aspect, the invention requires the extraction of a protein sample from a cell or tissue sample followed by the incubation of said sample with a labeled reagent capable of binding specifically to said sample (for example, an antibody) and detecting said reagent, wherein the marker which includes the reagent is selected from the group of a radioisotope, a fluorescent compound, an enzyme or an enzymatic cofactor.
Typical immunodetection methods include ELISA, RIA, immunoradiometric assay, fluoroimmunoassay, chemoluminescent assays, bioluminescent assays and Western blot assays.
Generally all the immunoassays include a step in which a sample suspected of containing a determined antigen or in which the concentration of said antigen is to be known is contacted with a first antibody in suitable conditions for the formation of the immune complexes. Suitable samples for the determination include a tissue section or biopsy, a tissue extract or a biological fluid. Once the antigen-antibody complexes have been formed, the preparation is subjected to one or several washings to remove antibodies that have not specifically bound.
Then, the detection of the immune complexes is performed by means of methods that are well known by the person skilled in the art, such as radioactive, fluorescent, or biological methods or methods based on the determination of an enzymatic activity.
For the purpose of increasing sensitivity, it is possible to use an additional ligand, such as a second antibody or a ligand coupled to biotin, for example. In this situation, an additional incubation step for incubating the ligand-antibody complexes obtained in the first step with the second antibody in suitable conditions for the formation of the secondary immune complexes is necessary. The secondary complexes are subjected to a washing cycle to remove secondary antibodies which have non-specifically bound, and then the amount of secondary immune complex is determined by means of determining the signal emitted by the secondary antibody.
Additional methods include the detection of the primary immune complexes by means of a two-step process. In this process, a secondary ligand (an antibody), which has binding affinity for the antibody forming part of the immune complexes, is contacted with said complexes to form secondary immune complexes, as was described above. After a washing step, the secondary immune complexes are contacted with a tertiary ligand or antibody which binds with high affinity to the secondary antibody to give rise to the formation of the tertiary complexes. The third ligand or antibody is bound to a detectable marker which allows the detection of the tertiary immune complexes.
Other detection methods include Western blotting, dot blotting, FACS analysis and the like. In one embodiment, the antibodies directed against the antigens of the invention are immobilized on a surface showing affinity for the proteins (for example polystyrene). Then a composition in which the antigen to be detected is present is added. After a washing step to remove the non-specifically bound complexes, the bound antigen can be detected by means of a second antibody which is coupled to a detectable marker. This type of ELISA is referred to as sandwich ELISA. The detection can also be carried out by means of adding a second antibody and a third antibody having affinity for the second antibody and which is bound to a detectable marker.
In another type of ELISA, the samples containing the antigen are immobilized and are detected by means of a competitive method in which the sample in which the antigen to be detected is present is mixed with antibodies labeled for said antigen and is added on the surface in which the antibody is immobilized. The presence of antigen in the sample prevents the binding of the antibody to the immobilized antigen such that the amount of antibody that binds to the immobilized antigen is present in an inverse proportion with respect to the amount of antigen in the sample to be analyzed.
It is also possible to detect the antigen by means of immunohistochemistry and confocal microscopy in tissue sections obtained from frozen samples, fixed in formaldehyde or embedded in paraffin using techniques that are widely known by the person skilled in the art.
The reference sample which is used for determining the variation of the expression levels of the genes and proteins used in the present invention. In one embodiment, the reference value is obtained from the signal provided using a tissue sample obtained from a healthy individual. Samples are preferably taken from the same tissue of several healthy individuals and are pooled, such that the amount of mRNA or of polypeptides in the sample reflects the mean value of said molecules in the population.
The method of the present invention can be combined with other diagnostic methods (e.g. oligoclonal bands in the CSF, neuroimaging (MR, OCT), clinical variables (disability scales, rate of flare-ups, age, sex) or biological markers: a) genetic markers (polymorphisms, haplotypes); b) pathological patterns in biopsy; c) antibodies, etc.
The methods of the present invention are particularly useful for establishing the prognosis in patients who have suffered a single flare-up of multiple sclerosis, in a patient suffering RR-MS or in a patient suffering PP-MS. This method would therefore be applied once during task of diagnosing the disease. It could also be applied to patients with the disease already diagnosed but in which, given the great variability of the disease, it is unknown if the disease is stable or if it will progress, again with a prognostic nature and to decide on the treatment. Therefore, the vast majority of patients would take the test at least once, except those with the disease in a very advanced stage in which the bad progress is already obvious and in which there are not possibilities of choosing between treatments. A sub-group of patients could take the test on several occasions if, over the years, the clinical course of the disease seems to change and the prognosis is to be re-assured.
In principle, any sample isolated from a patient can be used in the methods of the invention. Thus, the determination of the mRNA or polypeptide levels can be performed in a tissue biopsy or in a biological fluid (serum, saliva, semen, sputum, CSF, tears, mucous, sweat, milk and the like). The determination can be carried out in tissue homogenates or in more or less clarified fractions thereof. In a preferred embodiment, the determination of the mRNA and polypeptide levels of the invention is carried out from mononuclear cell extracts obtained from peripheral blood.
In the event that the expression levels of several of the genes identified in the present invention are to be determined simultaneously, compositions containing at least one copy of a probe specific for each of the genes indicated in at least one table selected from the group of Tables 3, 5-8 and 11 are useful.
Thus, in another aspect, the invention relates to a kit comprising a probe specific for each of the genes indicated in at least one table selected from the group of Tables 3, 5-8 and 11.
“Kit” is understood in the context of the present invention as a product containing the different reagents necessary for carrying out the methods of the invention packaged so as to allow their transport and storage. Materials suitable for packaging the components of the kit include glass, plastic (polyethylene, polypropylene, polycarbonate and the like), bottles, vials, paper, sachets and the like. Additionally, the kits of the invention can contain instructions for the simultaneous, sequential or separate use of the different components in the kit. Said instructions can be in the form of printed material or in the form of an electronic medium capable of storing instructions such that they can be read by a subject, such as electronic storage media (magnetic discs, tapes and the like), optical media (CD-ROM, DVD) and the like. The media can additional or alternatively contain Internet addresses which provide said instructions.
In a preferred embodiment, the kit of the invention consists of a probe specific for each of the genes indicated in at least one table selected from the group of Tables 3, 5-8 and 11.
In a preferred embodiment, the genes forming part of the array are the genes indicated in Table 11 and at least one reference gene.
In another preferred embodiment, the probes or the antibodies forming the kit of the invention are coupled to an array.
In the event that the expression levels of several of the genes identified in the present invention are to be determined simultaneously, the inclusion of probes for all the genes the expression of which is to be determined in a hybridization microarray is useful.
The microarrays comprise a plurality of nucleic acids spatially distributed and stably associated with a support (for example, a biochip). The nucleic acids have a sequence complementary to particular subsequences of the genes the expression of which is to be detected, so they are capable of hybridizing with said nucleic acids. In the methods of the invention, a microarray comprising an array of nucleic acids is contacted with a nucleic acid preparation isolated from the patient object of study. The incubation of the microarray with the nucleic acid preparation is carried out in suitable conditions for hybridization. Subsequently, after the elimination of the nucleic acids that have not been retained in the support, the hybridization pattern is detected, which provides information about the genetic profile of the analyzed sample. Although the microarrays are capable of providing both qualitative and quantitative information of the nucleic acids present in a sample, the invention requires the use of arrays and methodologies capable of providing quantitative information.
The invention contemplates a variety of arrays in terms of type of probes and in terms of type of support used. The probes included in the arrays which are capable of hybridizing with the nucleic acids can be nucleic acids or analogs thereof which maintain the hybridization capacity, such as, for example, nucleic acids in which the phosphodiester bond has been replaced with a phosphorothioate, methylimino, methylphosphonate, forforamidate, guanidine bond and the like, nucleic acids in which the ribose of the nucleotides has been replaced with another hexose, peptide nucleic acids (PNA). The length of the probes can be 7, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 100 nucleotides and vary in the range of 10 to 1000 nucleotides, preferably in the range of 15 to 150 nucleotides, more preferably in the range of 15 to 100 nucleotides and they can be single-stranded or double-stranded nucleic acids.
The selection of the probes specific for the different target genes is carried out such that they bind specifically to the target nucleic acid with minimum hybridization to unrelated genes. However, there are probes 20 nucleotides in length which are not unique for a determined mRNA. Therefore, probes directed against said sequences will show cross-hybridization with identical sequences appearing in mRNA of unrelated genes. In addition, there are probes which do not specifically hybridize with the target genes in the conditions used (due to secondary structures or interactions with the substrate of the array). Probes of this type should not be included in the array. Therefore, the person skilled in the art will note that the probes that are going to be incorporated in a determined array must be optimized before their incorporation in the array. The optimization of the probes is generally carried out by generating an array containing a plurality of probes directed against the different regions of a determined target polynucleotide. This array is contacted first of all with a sample containing the target nucleic acid in an isolated manner and, second of all, with a complex mixture of nucleic acids. Probes showing a highly specific hybridization with the target nucleic acid but not a low or any hybridization with the complex sample are thus selected for their incorporation in the arrays of the invention. It is additionally possible to include in the array hybridization controls for each of the probes that is going to be studied. In a preferred embodiment, the hybridization controls contain an altered position in the central region of the probe. In the event that high levels of hybridization are observed between the studied probe and its hybridization control, the probe is not included in the array.
In a preferred embodiment, the array contains a plurality of probes complementary to subsequences of the target nucleic acid of a constant length or of a variable length in a range of 5 to 50 nucleotides. The array can contain all the specific probes of a determined mRNA of a determined length or it can contain probes selected from different regions of an mRNA. Each probe is assayed in parallel with a probe with a changed base, preferably in the central position of the probe. The array is contacted with a sample containing nucleic acids with sequences complementary to the probes of the array and the hybridization signal is determined with each of the probes and with the corresponding hybridization controls. Those probes in which a greater difference between the hybridization signal with the probe and its hybridization control is observed are selected. The optimization process can include a second optimization round in which the hybridization array is hybridized with a sample which does not contain sequences complementary to the probes of the array. After the second round of selection, those probes presenting hybridization signals less than a threshold level will be selected. Probes which exceed both controls, i.e., that show a minimum level of non-specific hybridization and a maximum level of specific hybridization with the target nucleic acid, are thus selected.
The microarrays of the invention contain not only probes specific for the polynucleotides indicative of a determined pathophysiological situation, but they also contain a series of control probes, which can be of three types: normalization controls, expression level controls and hybridization controls.
Normalization controls are oligonucleotides which are perfectly complementary to labeled reference sequences which are added to the nucleic acid preparation to be analyzed. The signals derived from the normalization controls after hybridization provide an indication of the variations in the hybridization conditions, intensity of the label, efficiency of the detection and another series of factors that can result in a variation of the hybridization signal between different microarrays. The signals detected from the remaining probes of the array are preferably divided by the signal emitted by the control probes thus normalizing the measurements. Virtually any probe can be used as a normalization control. However, it is known that the effectiveness of the hybridization ranges according to the nucleotide composition and the length of the probe. Therefore, preferred normalization probes are those which represent the mean length of the probes present in the array, although they can be selected such that they comprise a range of lengths which reflect the remaining probes present in the array. The normalization probes can be designed such that they reflect the mean nucleotide composition of the remaining probes present in the array. A limited number of normalization probes is preferably used, selected such that they suitably hybridize, i.e., they do not present a secondary structure and they do not show sequence similarity with any of the probes of the array. The normalization probes can be located in any position in the array or in multiple positions in the array to efficiently control variations in the hybridization efficiency related to the structure of the array. The normalization controls are preferably located in the corners of the array and/or in the center thereof.
The expression level controls are probes which specifically hybridize with genes which are constitutively expressed in the sample which is analyzed. The expression level controls are designed to control the physiological condition and the metabolic activity of the cell. The analysis of the covariance of the expression level control with the expression level of the target nucleic acid indicates if the variations in the expression levels are due to changes in the expression levels or if they are due to changes in the global transcription rate in the cell or in its general metabolic activity. Thus, in the case of cells presenting deficiencies in a determined metabolite essential for cell viability, it is expected that a reduction in both the expression levels of the target gene and in the expression levels of the control will be observed. In addition, if an increase in the expression of the target gene and of the control gene is observed, it is probable that it is due to an increase of the metabolic activity of the cell and not to a differential increase in the expression of the target gene. Any probe which corresponds to a constitutively expressed gene can be used, such as genes encoding proteins which perform essential functions of the cells, such as β-2-microglobulin, ubiquitin, ribosomal 18S protein, cyclophilin A, transferrin receptor, actin, GAPDH and the like. In a preferred embodiment, the expression level controls are GAPDH, tyrosine 3-monooxygenase/triptophan 5-monooxygenase activation protein (YWHAZ), ubiquitin, beta-actin and β-2-microglobulin.
The hybridization controls can be included for the probes directed against the target genes and for the probes directed against the expression level or against the normalization controls. Error controls are oligonucleotide probes identical to the probes directed against the target genes but they contain mutations in one or several nucleotides, i.e., they contain nucleotides in certain positions which do not hybridize with the corresponding nucleotide in the target gene. The hybridization controls are selected such that, by applying the suitable hybridization conditions, the target gene should hybridize with the specific probe but not with the hybridization control, or with a reduced efficiency. The hybridization controls preferably contain one or several modified positions in the center of the probe. The hybridization controls therefore provide an indication of the degree of non-specific hybridization or of cross-hybridization to a nucleic acid in the sample to a probe different from the one containing the exactly complementary sequence.
The arrays of the invention can also contain amplification and sample preparation controls which are probes complementary to subsequences of control genes selected because they typically do not appear in the biological sample object of study, such as probes for bacterial genes. The RNA sample is supplemented with a known amount of a nucleic acid which hybridizes with the selected control probe. The determination of the hybridization to said probe indicates the degree of recovery of the nucleic acids during its preparation as well as an estimate of the alteration caused in the nucleic acids during the processing of the sample.
Once a set of probes which show the suitable specificity and a set of control probes are available, the latter are arranged in the array in a known position such that, after the hybridization and detection steps, it is possible to establish a correlation between a positive hybridization signal and the particular gene from the coordinates of the array in which the positive hybridization signal is detected.
The microarrays can be high density arrays with thousands of oligonucleotides by means of in situ photolithographic synthesis methods (Fodor et al., 1991, Science, 767-773). Probes of this type are typically redundant, i.e., they include several probes for each mRNA to be detected.
In a preferred embodiment, the arrays are low density arrays, or LDA, containing less than 10000 probes in each one per square centimeter. In said low density arrays, the different probes are manually applied with the aid of a pipette in different locations of a solid support (for example, a glass surface, a membrane). The supports used for fixing the probes can be obtained from a wide variety of materials, including plastic, ceramic, metals, gels, membranes, glass and the like. The microarrays can be obtained using any methodology known by the person skilled in the art.
After hybridization, in the cases in which the non-hybridized nucleic acid is capable of emitting a signal in the detection step, a washing step is necessary to remove said non-hybridized nucleic acid. The washing step is carried out using methods and solutions known by the person skilled in the art.
In the event that the labeling in the nucleic acid is not directly detectable, it is possible to connect the microarray comprising the target nucleic acids bound to the array with the other components of the system necessary for producing the reaction which gives rise to a detectable signal. For example, if the target nucleic acids are labeled with biotin, the array is contacted with streptavidin conjugated with a fluorescent reagent in suitable conditions so that the binding occurs between the biotin and the streptavidin. After the incubation of the microarray with the system which generates the detectable signal, it is necessary to perform a washing step to remove all the molecules which have non-specifically bound to the array. The washing conditions will be determined by the person skilled in the art using suitable conditions according to the system which generates the detectable signal which has been used and which are well known by the person skilled in the art.
The resulting hybridization pattern can be viewed or detected in different ways, said detection being determined by the type of system used in the microarray. Thus, the detection of the hybridization pattern can be carried out by means of scintillation counting, autoradiography, determination of a fluorescent signal, calorimetric determinations, detection of a light signal and the like.
Before the detection step, it is possible to treat the microarrays with an endonuclease specific for single-stranded DNA, such that the DNA which has non-specifically bound to the array is removed, whereas the double-stranded DNA resulting from the hybridization of the probes of the array with the nucleic acids of the sample object of study remains unchanged. Endonucleases suitable for this treatment include S1 nuclease, Mung bean nuclease and the like. In the event that the treatment with endonuclease is carried out in an assay in which the target nucleic acid is not labeled with a directly detectable molecule (for example, in an assay in which the target nucleic acid is biotinylated), the treatment with endonuclease will be performed before contacting the microarray with the other members of the system which produces the detectable signal.
After hybridization and the possible subsequent washing and treatment processes, the hybridization pattern is detected and quantified, for which the signal corresponding to each hybridization point in the array is compared to a reference value corresponding to the signal emitted by a known number of terminal-labeled nucleic acids to thus obtain an absolute value of the number of copies of each nucleic acid that is hybridized at a determined point of the microarray.
In the event that the expression levels of several of the proteins identified in the present invention are to be determined simultaneously, compositions containing at least one antibody specific for each of the genes indicated in at least one table selected from the group of Tables 1 to 4 are useful. The antibody arrays such as those described by De Wildt et al. (2000) Nat. Biotechnol. 18:989-994; Lueking et al. (1999) Anal. Biochem. 270:103-111; Ge et al. (2000) Nucleic Acids Res. 28, e3, I-VII; MacBeath and Schreiber (2000) Science 289:1760-1763; WO 01/40803 and WO 99/51773A1, are useful for this purpose. The antibodies of the array include any immunological agent capable of binding to a ligand with high affinity, including IgG, IgM, IgA, IgD and IgE, as well as molecules similar to antibodies which have an antigen binding site, such as Fab′, Fab, F(ab′)2, single domain antibodies, or DABS, Fv, scFv and the like. The techniques for preparing said antibodies are well known by the person skilled in the art and include the methods described by Ausubel et al. (Current Protocols in Molecular Biology, eds. Ausubel et al., John Wiley & Sons (1992)).
The antibodies of the array can be applied at high speed, for example, using commercially available robotic systems (for example, those produced by Genetic Microsystems or Biorobotics). The substrate of the array can be nitrocellulose, plastic, glass or it can be made from a porous material such as, for example, archylamide, agarose or another polymer. In another embodiment, it is possible to use cells which produce the antibodies specific for detecting the proteins of the invention by means of their culture in array filters. After the inducing the expression of the antibodies, the latter are immobilized in the filter in the position of the array where the producer cell is arranged.
An antibody array can be contacted with a labeled target and the level of binding of the target to the immobilized antibodies can be determined. If the target is not labeled, a sandwich-type assay can be used which uses a second labeled antibody specific for the polypeptide which binds to the polypeptide which is immobilized in the support.
The quantification of the amount of polypeptide present in the sample in each point of the array can be stored in a database as an expression profile. The antibody array can be produced in duplicate and used for comparing the binding profiles of two different samples.
In another aspect, the invention relates to the use of a kit of the invention for determining the prognosis of a patient diagnosed with multiple sclerosis, for determining the effectiveness of a treatment for multiple sclerosis or for diagnosing multiple sclerosis in a patient.
The invention is described below by means of the following examples which must be considered as illustrative and non-limiting of the scope of the invention.
1. Screening with DNA Chips
6 multiple sclerosis patients were recruited, 3 of them were diagnosed as having a bad prognosis and 3 as having a good prognosis, and 3 healthy controls without any history of any autoimmune disease. The prognosis of the patients was determined by means of clinical data associated with the progress of multiple sclerosis in studies of the natural history of multiple sclerosis as a first flare-up type, time until the second flare-up, number of flare-ups in the first 2 to 5 years and initial sequelae (Table 1).
The purification of total RNA was performed from peripheral blood using the PAXgene™ Blood RNA Kit of PreAnalytiX. The use of this kit allows preserving the RNA expression profile after performing the blood extraction. During the purification of the total RNA, a treatment with DNase was performed to remove the possible DNA contamination. The samples were concentrated by means of Speed-vac and the quality and amount of the RNA purified was estimated by means of testing an aliquot in agarose gel and spectrophotometric measurement.
cDNA was synthesized from 6 μg of total RNA with the SuperScript Choice System Kit of Life Technologies, following the protocol of the Expression Analysis Technical Manual of Affymetrix. cRNA was synthesized from this cDNA following the protocol of the BioArray HighYield RNA Transcript Labeling Kit (T7) of Enzo. The cRNA thus synthesized was purified with the Clean-up module Kit of Affymetrix, being recovered in a final volume of 22 μl of water. Once synthesized and purified, the cRNA was fragmented (15 μg of each sample) to prepare the hybridization mixtures.
The hybridization and the development and scanning of the chips were performed following the protocols and equipment officially recommended by Affymetrix Inc. The chip used was HG-U133 Plus 2.0. The results of the chip were analyzed using the Microarray Suite 5.0 software (MAS 5.0; Affymetrix®) and Biometric Research Branch (BRB) Array Tools 3.2.3 (Dr. Richard Simon and Amy Peng Lam).
40 multiple sclerosis patients were recruited, 20 of them were diagnosed as having a bad prognosis and 20 as having a good prognosis, and 20 healthy controls without any history of any autoimmune disease using the clinical criteria described above.
The purification of total RNA was performed from peripheral blood mononuclear cells. By means of centrifuging with a density gradient using Ficoll-Paque of Pharmacia Biotech, the lymphocytes and monocytes were purified and immediately immersed in an RNAlater RNA Stabilization Reagent of Qiagen to preserve the gene expression patterns. The total RNA was purified using the RNeasy Mini Kit of Qiagen and during purification DNA residue was removed by means of treatment with DNase using the RNase-Free DNase Set of Qiagen. The synthesis of cDNA from total RNA was performed using the High-Capacity cDNA Archive Kit of Applied Biosystems.
The gene validation analysis was performed using the Low Density Arrays (LDAs; Applied Biosystems) technology. The LDAs contain 384 wells. The wells contain TaqMan assays validated by Applied Biosystems and the distribution of the assays is configurable by the user. In this project, the chosen plate design is 95 genes+1 control analyzed in duplicate and with two samples studied in each plate.
Taking into account that for LDAs only those assays which are inventoried by Applied Biosystems can be selected, the process for selecting the most suitable assays from the genes to be validated was according to the following criteria:
As a first step of the analysis, those samples presenting a standard deviation between replicas of the same PCR assay greater than 0.38 were ruled out. 0.38 is used as the limit value because it indicates that there is a difference of 0.75 between the minimum and the maximum Ct. Since each Ct is equivalent to a PCR cycle and in each cycle the amount of DNA is duplicated, the standard deviation of 0.38 indicates that there is almost twice the amount of DNA in one well than in the other. Then, and to calculate the expression values of each gene, the formula 2exp(Ctmin−Ctsample) is applied, where Ctmin is the minimum Ct value of each gene in all the samples and Ctsample is the Ct value of that gene in that sample.
By using the expression values of the constitutive genes after the processing (5 in this case; GAPDH, YWHAZ, UBC, ACTB, B2M) the normalization factor of each sample was calculated by means of the geNorm program (http://medgen.ugent.be/˜jvdesomp/genorm/), which will calculate the geometric mean of the expression value of a number of constitutive or internal control genes. These internal control genes were chosen according to those genes in which it was observed that there was less variation in their expression between the studied conditions in the gene expression analysis experiment in DNA chips previously performed. Once the normalization factor of each sample was obtained, the data of each gene in each sample was normalized with respect to this normalization factor obtained for said sample of the geNorm program.
For the statistical analysis, the data normalized with respect to the control genes were transformed to logarithmic scale (base 2). For the calculation of significant genes, a non-parametric test was applied, which will depend on if 2 (Mann-Whitney U for 2 independent samples) or 3 conditions (Kruskal Wallis H for 3 independent samples) are compared. In all the cases, the p values<0.05 were considered significant differences whereas the p values<0.01 were considered very significant differences. The statistical analysis was performed using the SPSS 11.0 program (SPSS Inc., Chicago, USA).
The Bayesian classifier was constructed using the Bayesian analysis software BayesiaLab 3.2 (Bayesia SA. Laval Cedex, France). To that end, the variables were previously discretized into a maximum of four ranges using the Decision Tree discretization algorithm taking the variable diagnosis as a reference. The learning was performed using the Augmented Naive Bayes algorithm.
1. Screening with DNA Chips
Table 2 shows the demographic and clinical characteristics of the multiple sclerosis patients and of the healthy controls used to perform the screening with DNA chips.
After normalizing the expression levels using MAS 5.0, the probes (genes) were filtered using the BRB Array Tools 3.2.3 according to the following criteria:
After filtering, 4,705 genes of the initial 54,675 complied with the criteria. A class comparison test was performed with these genes and 45 of them which presented significant differences between the three classes (control, bad and good prognosis) with a p value<0.001 (Table 3) were identified.
The distribution of these genes in Gene Ontology (GO) was not significantly different from that expected randomly (
A cluster analysis of the samples was performed with these 45 genes and it was observed that 3 highly diagnostic reproducible clusters with a mean value of correlation between each cluster of approximately 0.6 (
A cluster analysis was also performed with these 45 genes and it was observed that 4 clusters with a mean value of correlation of approximately 0.65 (
For the purpose of complementing the number of genes which allow differentiating between each diagnosis, a class comparison test with a p value<0.001 between the bad and good prognosis, control and bad prognosis, control and good prognosis and control and multiple sclerosis classes (Tables 5, 6, 7 and 8), as well as a class comparison test between the three diagnoses with a p value<0.005 (Table 9), were performed.
Homo sapiens,
The results of the cluster analysis both for samples and for genes obtained after these class comparisons are shown in
Table 10 includes the list of the 95+1 genes and assays selected to configure the LDAs from the results obtained in the screening with DNA chip.
The statistical analysis identified 25 genes which presented significant differences (p<0.01) in the expression levels between the three classes (
Using these 25 genes and the EDSS and MSFC clinical variables at the onset of the disease, a Bayesian classifier was constructed which showed a precision of 91.66% between the three diagnoses (
For the purpose of increasing precision when distinguishing between good and bad prognosis, a new classifier was constructed using the clinical variables and only those genes which presented significant differences (p<0.05) in the expression levels of both classes (
Table 12 shows the analysis of the information provided by each variable for establishing the prognosis. The weight represents the relative amount of information provided by each variable to the classification of the prognosis. The list of variables is arranged in descending order according to the information provided by each one. The a priori modal value describes the most probable value of each variable when the prognosis is unknown, whereas the modal value for each prognosis describes the most probable value for that prognosis. All the modal values are accompanied by their probability. The variation for each prognosis is a measure indicating the difference of probability between the a priori modal value and the modal value for the prognosis when the latter is known. The formula used for calculating it is: −log 2(P(modal value for the prognosis))+log2(P (modal value for the prognosis|value observed)). The simple underlined values simple indicate positive variations (the probability of the modal value for the prognosis is greater than that of the a priori modal value) whereas the values in italics indicate negative variations. Obviously no variation is indicated if the modal value for the prognosis is different from the a priori modal value. The modal value for the prognosis is then represented in bold print.
<=−3.145
0.6130
>2.000
0.7655
<=−1.920
<=−3.295
0.5536
−0.4948
−0.4948
0.3677
<=−2.825
0.4406
−0.3440
0.2775
−0.3440
0.2775
>−2.335
>−2.650
−0.2801
0.2345
−0.3692
0.2937
−0.2224
0.1926
<=−3.720
0.4285
Table 13 presents the precision of the classifier as the number of variables making it up increases. The variables are introduced in order according to the information provided by each one for establishing the prognosis. Within the prognoses the values are presented as correctly classified individuals (correct) with respect to the incorrectly classified individuals (incorrect).
Table 14 presents the conditional probabilities for the prognosis of the disease for each variable used by the classifier.
Number | Date | Country | Kind |
---|---|---|---|
P200703415 | Dec 2007 | ES | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/ES2008/000785 | 12/18/2008 | WO | 00 | 11/15/2010 |