Steroid responsive nucleic acid expression and prediction of disease activity

Information

  • Patent Application
  • 20070248978
  • Publication Number
    20070248978
  • Date Filed
    April 09, 2007
    18 years ago
  • Date Published
    October 25, 2007
    17 years ago
Abstract
The invention relates to methods useful for diagnosing and monitoring the steroid responsiveness of a subject by detecting expression of steroid modulated genes and for predicting transplant rejection and non-rejection.
Description
TECHNICAL FIELD

The invention relates to methods for detecting nucleic acid and protein expression modulated by steroids and using steroid responsiveness of a subject in predicting and monitoring disease activity.


BACKGROUND OF THE INVENTION

Steroids are used to ameliorate disease activity associated with immune disorders such as graft rejection, systemic lupus erythematosis (SLE), multiple sclerosis (MS) and cytomegalovirus (CMV) infection. Although steroids are used clinically to treat hyperactivity of the immune system, prolonged treatment has deleterious effects including diabetes, osteoporosis and weight gain. Given these and other side effects, clinicians avoid prescribing high dosages of steroid any longer than necessary. Since flare of immune disorders and transplantation require the use of steroids as an ongoing treatment, it is desirable to determine the steroid responsiveness of a subject in order to optimize outcome. An essential component of providing effective immunosuppression is monitoring subject or transplant status. In transplant patients, this monitoring is organ, tissue or cell-specific. For example, monitoring a subject with a cardiac transplant involves taking a biopsy of heart muscle and having a pathologist examine it for cytological evidence of rejection. Such biopsies are expensive, invasive, and painful and interpretation can only be focused on the biopsied cells, not the whole organ.


Although glucocorticoid induction of genes correlated with immune response has been studied in vitro (Galon et al. (2002) FASEB Journal 16:61-71); there is a need for methods to detect in vivo expression of steroid modulated nucleic acids. The present invention addresses this need by diagnosing and monitoring steroid responsiveness or immunological status, predicting flares or graft rejection, and designing, evaluating or monitoring treatment efficacy.


SUMMARY OF THE INVENTION

The present invention provides methods for detecting in vivo expression of nucleic acids and proteins modulated by steroid administration and metabolism. The invention presents a method of diagnosing or monitoring steroid responsiveness of a subject comprising detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression is correlated with steroid administration or dosage and applying at least one statistical method to the expression of the diagnostic set to diagnose or monitor steroid responsiveness of the subject.


In one embodiment, the diagnostic set further comprises at least one steroid modulated nucleic acid selected from each of at least two of the clusters of Table 1. In a second embodiment, the diagnostic set further comprises two or more steroid modulated nucleic acids selected from Table 2. In a third embodiment, the diagnostic set further comprises two or more steroid modulated nucleic acids selected from Table 3. In one aspect, detecting the expression of the diagnostic set of steroid modulated nucleic acids further comprises using hybridization or quantitative real-time polymerase chain reaction (RT-PCR) and a sample obtained from the subject by any sampling means. In a second aspect, the sample is a blood sample, and RNA is isolated from the peripheral blood mononuclear cells (PMBC) of the blood sample. In a third aspect, the statistical method is K-means clustering that produces clusters of genes that are correlated by p-value and their expression in a cell type or pathway or a prediction algorithm selected from a linear algorithm, a logistic regression algorithm, and a voting algorithm that produces a single value or score.


In a fourth embodiment, the diagnostic set further comprises selecting at least two oligonucleotides or a probe set to detect the expression of each steroid modulated nucleic acid of the diagnostic set. The invention also presents a kit comprising the oligonucleotides or probe sets that detect the expression of each steroid modulated nucleic acid of the diagnostic set. The invention further presents a method for diagnosing or monitoring steroid responsiveness of a subject comprising detecting the expression of nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.


The invention additionally presents a method for predicting rejection or non-rejection in a subject with a transplant comprising detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression of the steroid modulated nucleic acids correlates with transplant rejection or non-rejection, and applying at least one statistical method to the expression of the diagnostic set of steroid modulated nucleic acids to predict rejection or non-rejection.


In one embodiment, the diagnostic set of steroid modulated nucleic acids further comprises two or more nucleic acids selected from Tables 1-3. In one aspect, detecting the expression of the diagnostic set of steroid modulated nucleic acids further comprises using RT-PCR and RNA isolated from PMBCs. In a third aspect, the statistical method is a prediction algorithm selected from a linear algorithm, a logistic regression algorithm, and a voting algorithm that produces a single value or score that correlates with rejection or non-rejection. In a fourth aspect, the score that correlates with non-rejection is≦20 and the score that correlates with rejection is≧30. The invention yet further presents a method of predicting rejection or non-rejection comprising detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.


The invention further presents a method of diagnosing or monitoring the status of a subject with a transplant comprising detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression is correlated with dysfunction or rejection of the transplant, and applying at least one statistical method to the expression of the nucleic acids to monitor the status of the transplant. In one embodiment, the diagnostic set further comprises two or more nucleic acids selected from Tables 1-3. In a second embodiment, RT-PCR is used with RNA isolated from PMBC to detect expression of the steroid modulated nucleic acids and the expression is analyzed using a prediction algorithm that produces single value or score that correlates with the status of the subject with the transplant. In a third embodiment, diagnosing and monitoring the status of a subject with a transplant further comprises detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.


The invention also presents method for designing and monitoring a treatment plan for a subject with a transplant or an immune disorder comprising detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression correlates with the steroid responsiveness of the subject, and using the expression of the diagnostic set of steroid modulated nucleic acids to design and monitor the treatment plan of the subject. In one embodiment, the diagnostic set of steroid modulated nucleic acids comprises two or more nucleic acids selected from Tables 1-3. In a second embodiment, RT-PCR is used with RNA isolated from PMBC to detect expression of the steroid modulated nucleic acids and the expression is analyzed using a prediction algorithm that produces single value or score that correlates with the steroid responsiveness of the subject. In a third embodiment, diagnosing and monitoring the status of a subject with a transplant or immune disorder further comprises detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3 whose expression correlates with steroid responsiveness of a subject. In one aspect, the transplant is selected from bone marrow, heart, kidney, liver, lung, pancreas, pancreatic islets, stem cells, xenotransplants, and artificial implants. In another aspect, the immune disorder is selected from cytomegalovirus infection, multiple sclerosis, and systemic lupus erythematosus.


The invention yet still further presents a method for using primers and probe sets to detect steroid responsiveness of a subject with a transplant or an immune disorder comprising designing and generating primers or probe sets for nucleic acids whose expression is modulated by steroid administration or dosage, and using RT-PCR and the primers or probe sets on a sample from the subject to detect steroid responsiveness. In one embodiment, the nucleic acids whose expression is modulated by steroid administration or dosage are selected from Tables 1-3. In a second embodiment, the primers and probe sets are used in a diagnostic kit.


BRIEF DESCRIPTION OF THE TABLES

Table 1 presents ten clusters of genes whose nucleic acid and protein expression is modulated by steroids. Column 1 shows cluster number; column 2, microarray probe ID from Human Genome CGH 44A Microarray (Agilent Technologies); column 3, gene symbol; column 4, average p-value for expression of the nucleic acid in CARGO and LARGO; column 5, average Pearson correlation for expression of the nucleic acid in CARGO and LARGO; column 6, p-value for the expression of the nucleic acid in CARGO, column 7, p-value for the expression of the nucleic acid in LARGO; and column 8, the name of the gene as it appears in the GenBank database (NCBI, Bethesda Md.).


Table 2 summarizes steroid modulated nucleic acid expression for 104 subject post-transplant samples and a subset of 74 samples≦180 days post-transplant. Column 1 shows the nucleic acids whose probe sets were used in RT-PCR to detect expression in post-transplant subject samples. The overall score refers to the single value produced from all scores using a linear discriminant algorithm. Columns 2-5 show the data for rejection (R) subjects, non-rejection (NR) subjects, the ratio, and p-values for all days post-transplant (index), respectively. Columns 6-9 show the data for rejection (R) subjects, non-rejection (NR) subjects, the ratio, and p-values for<180 days post transplant samples (subset), respectively. Significant p-values are shown in red typeface.


Table 3 presents RT-PCR data for 33 nucleic acids expressed in pathways having genes modulated by steroids or regulating T-cell homeostasis. Column 1 of Table 3 shows the gene symbol; columns 2 and 3, the fold change and p-value for R (n=38)/NR (n=55) at all times post-transplant; columns 4 and 5, the fold change and p-value for R (n=27)/NR (n=40) at≦180 days post-transplant; and column 6, the gene name.







DETAILED DESCRIPTION OF THE INVENTION

The present invention addresses needs in the art by providing methods for detecting the in vivo expression of nucleic acids modulated by steroid administration or metabolism. The invention also provides methods for diagnosing and monitoring steroid responsiveness of a subject by detecting the expression of nucleic acids modulated by steroids. The invention uses detection of nucleic acids modulated by steroids to predict disease activity or transplant non-rejection or rejection and to determine status of an immune disorder or transplant. Such methods can be used to fine-tune immunosuppressant therapy and, more importantly, to reduce the number of invasive and costly tests and procedures that a subject must undergo. In particular, the invention can be used to predict transplant non-rejection or rejection. For example the invention can be used to predict transplant non-rejection or rejection allowing a clinician to reduce the number of biopsies performed in the first 180 days post-transplant or to begin anti-rejection therapy before cytological evidence of rejection is detectable. The invention also provides methods for evaluating the need for post-transplant monitoring and treatment or determining a subject's near-term prognosis based on steroid modulated nucleic acid expression.


Definitions


Unless defined otherwise, all scientific and technical terms are understood to have the same meaning as commonly used in the art to which they pertain. For the purpose of the present invention, the following terms are defined.


“Amplification” refers to any device, method or technique that can make copies of a nucleic acid. It can be achieved using a thermal cycler or a thermal gradient device and a polymerase chain reaction (PCR) technique such as linear amplification (cf. U.S. Pat. No. 6,132,997), rolling circle amplification, and the like. Further, amplification and detection can be combined as in Real-Time PCR (RT-PCR) using TAQMAN protocols and the Prism 7900HT Sequence Detection system and software (Applied Biosystems (ABI), Foster City Calif.).


“Array” refers to an ordered arrangement of at least two samples—nucleic acids, proteins or antibodies—in solution or on a substrate where at least one of the samples represents a control and/or normal sample and the other, a sample of diagnostic or prognostic interest. The ordered arrangement ensures that the size and signal intensity of each labeled complex, formed between at least one reagent and at least one sample to which the reagent specifically binds is readily detectable.


“Clusters” refers to groups of nucleic acids with expression that is directly or indirectly regulated by and correlated with the administration or metabolism of a steroid.


“Diagnostic set” refers to at least two nucleic acids whose expression is modulated by steroids and whose nucleic acids, oligonucleotides, primers and probe sets can be used in nucleic acid technologies or encoded proteins and antibodies or affinity reagents thereto can be used in protein technologies.


“Expression” refers to differential expression—increased or decreased expression as detected by presence, absence, or change in the amount of nucleic acid or protein expressed in a sample—as presented in a gene expression profile. A “gene expression profile” refers to the identification, characterization, quantification, and representation of a plurality of nucleic acids expressed in a sample from a subject as measured using nucleic acid or protein technologies. Nucleic acid expression is detected using nucleic acid technologies and mature mRNA transcript and/or regulatory sequences such as promoters, enhancers, introns, mRNA-processing intermediates, and 3′ untranslated regions. A gene expression profile from a subject can be compared with reference gene expression profiles based on detection of nucleic acid expression in control or normal, diseased, or treated samples.


“Immune disorders” refers to conditions, disorders and diseases associated with immunological response including but not limited to acute respiratory distress syndrome, Addison's disease, allograft rejection, ankylosing spondylitis, Takayasu's arteritis, arteriosclerosis, asthma, atherosclerosis, congestive heart failure, primary sclerosing cholangitis, Churg-Strauss syndrome, CREST syndrome, Crohn's disease, ulcerative colitis, diabetes mellitus, emphysema, glomerulonephritis, Wegener's granulomatosis, Grave's disease, autoimmune hepatitis, Kawasaki's syndrome, systemic lupus erythematosus, multiple sclerosis, myasthenia gravis, myelofibrosis, pancreatitis, polyarteritis nodosa, polymyositis, psoriasis, Raynaud's disease, Reiter's syndrome, rheumatoid arthritis, scleroderma, primary biliary sclerosis, systemic sclerosis, sepsis, septic shock syndrome, Sjogren's disease, ankylosing spondylitis, primary thrombocythemia, Hashimoto's thyroiditis, systemic vasculitis, Whipple's disease, complications of cancer, viral infection including CMV infection, bacterial infection, fungal infection, parasitic infection, protozoal infection, helminthic infection, and trauma.


“Immunosuppressant” refers to any therapeutic agent that suppresses immune response in a subject such as anticoagulents, antimalarials, heart drugs, non-steroidal anti-inflammatory drugs, and steroids including but not limited to aspirin, azathioprine, chloroquine, corticosteroids, cyclophosphamide, cyclosporin A, dehydroepiandrosterone, deoxyspergualin, dexamethasone, everolimus, fenoprofen, hydralazine, hydroxychloroquine, immunoglobulin, ibuprofen, indomethacin, leflunomide, ketoprofen, meclophenamate, mepacrine, 6-mercaptopurine, methotrexate, mizoribine, mycophenolate mofetil, naproxen, prednisone, methyprenisone, rapamycin (sirolimus), solumedrol, tacrolimus (FK506), thymoglobulin, tolmetin, tresperimus, triamcinoline, and the like.


“Monitoring” refers to repetitive testing for and detection of nucleic acid expression that provides useful information about a subject's health or disease status. Monitoring can include determining prognosis, risk-stratification, and efficacy of a particular drug; detecting subject response to a drug or ongoing therapy; predicting susceptibility, rejection or non-rejection, or disease activity; diagnosing onset, flare or complication of a disease; following disease progression or providing information related to a subject's status over time; selecting subjects most likely to benefit from a particular drug or experimental therapy especially where administration of that drug works for a small subset of subjects or where the drug does not have a label for a particular immune disorder; and screening a subject population to decide to use a more or less invasive or costly test; for example, moving from a non-invasive blood test to a more invasive option such as biopsy.


“Nucleic acid technology” refers to any and all devices, methods and systems used to detect expression of nucleic acids and produce a gene expression profile including but not limited to methods using arrays for hybridization, amplification in PCR, quantitative RT-PCR, TAQMAN protocol RT-PCR, multiplex PCR, thermal gradient devices, and the like, or hybridization in solution or on a substrate containing cDNAs, genomic DNAs, locked nucleic acids (LNAs), oligonucleotides, primers, peptide nucleic acids, polynucleotides, probe sets, RNAs and the like.


“Prediction” or “predicting” refers to the use of gene expression profile to provide information about a subject's health or the status of a disease, patient or transplant and can include determination of prognosis, risk-stratification, prediction of outcomes, and the like.


A “probe set” refers to groups of oligonucleotides or primers that can be used with a nucleic acid technology to detect groups of two or more nucleic acids. Primers in a probe set can contain rare or artificial nucleotides, be of any size useful in a nucleic acid technology, designed to detect a particular region or splice variant of a gene, labeled with one or more detectable moieties, and used in solution or attached to a substrate.


“Protein technology” refers to any and all devices, methods and systems that can be used to detect a peptide, polypeptide or protein expressed by a steroid modulated nucleic acid or gene and produce a gene expression profile including but not limited to activity assays, affinity assays, antibody or protein arrays, chromatographic separation, calorimetric assays, two-dimensional gel electrophoresis, ELISA, fluorescent-activated cell sorting, mass spectrophotometric detection, protein-fusion reporter constructs, western analysis, and the like. Protein expression, although time delayed, is correlated with and mirrors nucleic acid expression.


“Sample” is used in its broadest sense and refers to any biological material used for cytological or histological evaluation or to measure nucleic acid expression and obtained from a subject by any sampling means known to those of skill in the art. A sample can comprise a bodily fluid such as ascites, bile, blood, cerebrospinal fluid, synovial fluid, lymph, pus, semen, sputum, urine; the soluble fraction of a cell preparation, an aliquot of media in which cells were grown; a chromosome, an organelle, or membrane isolated or extracted from a cell; cDNA, genomic DNA, or RNA including but not limited to hnRNA, mRNA, mRNA processing intermediates, rRNA, and tRNA in solution or bound to a substrate; a cell; a cell, tissue or organ biopsy, and the like. Preferred samples for diagnosis, prognosis, or monitoring of immunological status are leukocytes, peripheral blood mononuclear cells (PBMC), or serum derived from whole blood.


“Sampling means” refers to any instrumentation and protocols for obtaining a biological sample and includes but is not limited to aspiration of a body fluid, aspiration of fluid following lavage, a biopsy (bronchoscopy or endoscopy) of cells, a tissue or organ, drawing of central or peripheral blood, and the like.


A “statistical method” refers to methods including but not limited to analysis of variance, canonical analysis, classification algorithms, classification and regression trees, cluster analysis including K-means clustering, factor analysis, Fisher's Exact test, k-nearest neighbor, linear algorithm, linear discriminatory analysis, linear regression, logistic algorithm, multidimensional scaling analysis, multiple regression, nearest shrunken centroids classifier, Pearson correlation, prediction algorithm, significance analysis of microarrays, one-tailed T-tests, two-tailed T-tests, voting algorithm, Wilcoxon's signed ranks test, and the like.


“Status” refers to any and all aspects of immune response in a subject who has an immune disorder or transplant including deterioration, improvement, progression, remission, or stability as determined from analyzing one or more samples from that subject for nucleic acid or protein expression that correlates with the degree and nature of response, steroid treatment or related complications including autoimmune cellular destruction, acute rejection, chronic rejection, humoral rejection, vasculopathy, and the like.


“Steroid modulated” refers to any gene product, nucleic acid or protein, whose expression is correlated with and results directly or indirectly from the administration or metabolism of steroids. For example, genes that have a steroid dependent regulatory element (sdre) in their promoter region (Dillner and Sanders (2002) J Biol Chem 277:33890-33894) are steroid modulated, primary response genes regulated by the presence and/or dosage of steroids. These primary response genes are often transcription factors that activate one or more indirectly affected, secondary response genes or pathways.


“Steroid responsive” or “steroid responsiveness” refers to any aspect of the immunological response of a subject to the administration or metabolism of steroids including improvement or worsening of symptoms, adjustment in dosage, change to another immunosuppressant, and the like.


“Subject” refers to an individual or patient who develops an infection, has an immune disorder or has received any allograft that elicits an immune response.


“Substrate” refers to any rigid or semi-rigid support to which antibodies, nucleic acids or proteins are bound and includes magnetic or nonmagnetic beads, capillaries or other tubing, chips, fibers, filters, gels, membranes, microparticles, plates, polymers, slides, and wafers with a variety of surface forms including channels, columns, pins, pores, trenches, wells and the like made from any natural or synthetic material or combination thereof.


“Transplant” refers to a subject's own genetically modified cells, or tissues grown from those cells; cells, tissues or organs from another subject or from an animal of a different species; and artificial implants such as mechanical or partially mechanical replacement organs.


“Transplant rejection” as detected or predicted using the methods and materials of the present invention refers to the rejection of bone marrow, heart, kidney, liver, lung, pancreas, pancreatic islets, stem cells, xenotransplants, and artificial implants.


“Quiescence” refers to the absence of signs or symptoms of histological or immunological response. For example, a diagnosis of remission in a subject with an immune disorder or non-rejection in a transplant patient indicates successful repression of immunological response and/or treatment with an immunosuppressant.


Description of the Invention


The correlation between the administration of steroids and the differential expression of steroid modulated nucleic acids and proteins provides an opportunity for developing pharmacogenomic markers for diagnosing and monitoring subjects with transplants, immune disorders such as SLE or MS, and CMV infection. As described in the Examples, the present invention provides methods, diagnostic sets of steroid modulated nucleic acids selected from Tables 1-3, and reagents such as antibodies, affinity reagents, primers and probe sets that can be used for determining, diagnosing, evaluating, monitoring, or predicting disease activity, non-rejection, rejection, status of a transplant or of an immune disorder, steroid responsiveness, and treatment plan of a subject with a transplant or immune disorder. In one embodiment, the ability to predict acute rejection can be used to begin immediate anti-rejection therapy while the ability to predict non-rejection can be used to determine the need for and timing of costly and invasive procedures such as biopsies. The invention additionally provides methods for designing and monitoring a treatment plan for a subject with an immune disorder or transplant and for evaluating the need for post-diagnosis or post-transplant monitoring and treatment.


The methods of the invention used RNA isolated from PBMC samples obtained from subjects enrolled in the Cardiac Allograft Rejection Gene Expression Observational (CARGO) and the Lung Allograft Rejection Gene Expression Observational (LARGO) studies. The samples were processed as described in Example 8 and used to study gene expression using nucleic acid technologies.


Microarray studies of gene expression were performed using the protocols described in Examples 9 and 10. These studies identified steroid modulated nucleic acids in the CARGO and LARGO samples from subjects treated with 1-100 mg doses of steroid as described in Example 1. Iterative cluster analysis and similarity testing as described in Example 4 were used to identify the nucleic acids modulated by steroids presented in Table 1. An exemplary RT-PCR study, carried out using the protocols described in Examples 13, used probe sets for 20 informative genes to investigate steroid responsiveness CARGO samples. The results of this study, as described in Example 5 and presented in Table 2, revealed that the expression of the nucleic acids known to be modulated by steroids were important both in diagnosing and monitoring steroid responsiveness and in predicting transplant rejection and non-rejection.


When the data from the exemplary RT-PCR study showed that differential expression of nucleic acids encoding FLT3, IL1R2, ITGAM, and PDCD1 proteins or fragments thereof was highly predictive of non-rejection within 180 days of transplant, a second study was performed to test additional nucleic acids in related pathways. Table 3 presents the results of the pathways RT-PCR study on 33 genes in the IL-1 or PDCD1 pathway, the ligand for FLT3, and genes induced and expressed in T cells. Using a p-value <0.05, the genes encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3 protein and fragments thereof showed differential expression correlated with rejection.


Primers or probes sets that detect expression of at least one nucleic acid from the steroid modulated genes of Tables 1-3 can be used in a diagnostic set to carry out the methods of the invention. In one embodiment, the steroid modulated nucleic acids of the invention were used to design, select, and test primers and probe sets that can be used to detect steroid responsiveness in a sample from a subject as described in Examples 11 and 12. In another embodiment, antibodies or other affinity reagents specifically binding to a protein or a fragment thereof, expressed from steroid modulated genes of Tables 1-3, can be used in a diagnostic set to carry out the methods of the invention. Protein expression and antibody production and testing are described in Examples 15 and 16.


In a preferred embodiment, the methods and diagnostic sets of this invention can be used on clinically stable subjects, those showing no histological signs of rejection in endomyocardial biopsies (EMB) within 180 days of transplant to predict the probability that transplant rejection will occur within the subsequent 12 weeks. For example, a prediction algorithm was applied to the nucleic acid expression from exemplary RT-PCR studies to produce a single score for each subject. Then quartile analysis was applied to the single scores as described in Example 6. When used in longitudinal studies of≦180 days post transplant, the score produced by the algorithm distinguished clinically stable cardiac transplant subjects who did not reject their transplant in the subsequent 12 weeks, 98.9% with a score≦20, from those who progressed to acute cellular rejection, 58% with a score≧30.


Using a nucleic acid technology or a protein technology to generate a gene expression profile, one of skill in the art would select the appropriate devices and methods based upon such factors as the particular immune disorder or transplant, ease and needed accuracy of measurement of each particular nucleic acid or protein, the number of primers, probe sets or antibodies in the diagnostic set, and the like. It is contemplated that a gene expression profile based on a small diagnostic set of steroid modulated nucleic acids can be produced on a low density array or a thermal gradient chip in a clinic or a doctor's office.


Knowing steroid modulated in vivo expression of nucleic acids or proteins and establishing a correlation between their expression and steroid responsiveness, one of skill in the art can use diagnostic sets of primers, probe sets, antibodies and the like to determine the steroid responsiveness of a particular individual. To establish such correlations, nucleic acid or protein expression will be measured multiple times, and statistical methods or algorithms will be applied to determine the reliability of the measurement and to establish a threshold for the correlation. Correlations can be determined using samples from steroid responsive subjects. For example, knowing the steroid modulated in vivo expression levels of nucleic acids and an established correlation between the expression levels of such nucleic acids and the steroid responsiveness of a group of subjects being treated for transplant rejection, one of skill in the art can extrapolate the steroid responsiveness of a previously untested subject.


The responsiveness of a subject, based on nucleic acid expression, can be used design or to modify a treatment plan including types and amounts of immunosuppressants or steroids being administered; the dose, frequency and duration of treatment; weaning protocol, and the like. If a subject develops or shows resistance to a particular immunosuppressant, nucleic acid or protein expression and established correlations or profiles can be used to re-evaluate, the subject's responsiveness and to revise the subject's treatment plan.


Reagents used to establish a gene expression profile include but are not limited to genes and their splice variants, amplicons, LNAs, oligonucleotides, peptide nucleic acids, primers, and probe sets that can be used in nucleic acid technologies; and proteins and their fragments, antibodies, and affinity reagents that can be used in protein technologies. These reagents can be used in assays or diagnostic kits to determine or monitor steroid responsiveness of a subject, to screen or monitor subjects for the development or flare of immune disorder or for transplant rejection, to design or evaluate a treatment protocol, and the like.


Assays or diagnostic kits based on the reagents and nucleic acid and protein technologies described herein can be used with a sample from a subject to diagnose, classify or rule out an immune disorder such as SLE or MS; to select a clinical trial, to predict flare, to detect immunosuppressant or steroid responsiveness, to determine efficacy of a potential therapeutic agent, to design treatment regimes, to monitor the status of the subject or the treatment regime. In one alternative, the diagnostic kit comprises an array of reagents; in another, probe sets for use in RT-PCR.


Pharmacogenomics is the study of an individual's response to a particular therapeutic agent, immunosuppressant or combinations thereof. In this context, response refers to whether a particular drug will work better for a subject with a particular immune disorder or transplant. The methods disclosed provide for assigning a subject to a clinical trial or treatment regime based on disease or transplant status (quiescent or flare for immune disorder, rejection or non-rejection for transplant). Pharmacogenomics is also important in determining the dosage of a therapeutic agent based on age, classification and status of the subject. Individual steroid responsiveness, dosage and even timing of administration must be taken into account relative to side effects or potential interactions of various therapeutic agents. Some potentially useful therapeutic agents, immunosuppressants and steroids are listed in the definitions and/or claims.


All of the references cited are hereby incorporated by reference herein. This invention will be better understood by reference to the following non-limiting Examples which serve to demonstrate the use of nucleic acid and protein expression to evaluate steroid responsiveness in subjects, to optimize steroid dosage, to predict periods of non-rejection in subjects with transplants in order to reduce the number of invasive procedures, EMBs, TBBs, and the like.


Tables 1, 2 and 3 described in detail in the Examples are provided below.

TABLE 1AverageP valueClusterProbe IdGene SymbolP valuePearsonCARGOLARGOGene Name1A_24_P146211HIST1H2BD0.0000020.851.71E−080.000257histone 1, H2bd1A_23_P59069HIST1H2BO0.0000060.861.77E−070.000175histone 1, H2bo1A_23_P366216HIST1H2BH0.0000080.838.15E−080.00073histone 1, H2bh1A_24_P55148HIST1H2BJ0.0000100.861.24E−068.55E−05histone 1, H2bj1A_23_P30776HIST1H2BE0.0000100.852.64E−070.000406histone 1, H2be1A_23_P42178HIST1H2BF0.0000120.857.09E−070.000195histone 1, H2bf1A_23_P402081HIST1H2BN0.0000130.845.66E−070.000315histone 1, H2bn1A_24_P156911HIST2H2BE0.0000150.868.42E−070.000269histone 2, H2be1A_23_P8013HIST1H2BL0.0000190.843.65E−069.88E−05histone 1, H2bl1A_23_P111054HIST1H2BB0.0000210.841.39E−060.000304histone 1, H2bb1A_23_P93180HIST1H2BC0.0000230.83.31E−070.00159histone 1, H2bc1A_24_P152345LOC3915660.0000260.762.73E−080.0247Histone H2B.n1A_32_P578540.0000260.841.35E−060.0005DKFZp586A07221A_24_P3783HIST1H2BM0.0000300.792.39E−070.00379histone 1, H2bm1A_23_P111041HIST1H2BI0.0000310.85.68E−070.00166histone 1, H2bi1A_23_P30020PLA2G12A0.0000500.861.45E−050.000173phospholipase A2, group XIIA1A_23_P218131C14orf1510.0000570.875.99E−055.34E−05chromosome 14 ORF 1511A_23_P332992HIST3H2BB0.0000590.849.35E−060.000368histone 3, H2bb1A_23_P256618C6orf790.0000740.83.98E−060.00136chromosome 6 open readingframe 791A_24_P10884GRAP20.0000740.762.09E−070.026GRB2-related adaptor protein 21A_23_P167997HIST1H2BG0.0000790.782.72E−060.00229histone 1, H2bg1A_32_P100439Ells10.0000860.883.41E−050.000215hypothetical protein Ells11A_24_P219785CALM30.0000860.798.90E−070.00837calmodulin 31A_24_P164718MARCH20.0000880.811.13E−050.000679membrane-associated ring finger(C3HC4) 21A_23_P154065TUBA10.0001170.832.21E−050.000615tubulin, alpha 11A_23_P258093AGPAT10.0001470.835.21E−050.0004151-acylglycerol-3-phosphate O-acyltransferase 11A_23_P410312FLJ401420.0001510.772.23E−060.0102FLJ40142 protein1A_23_P29124GP1BB0.0001560.845.15E−050.000473glycoprotein Ib, beta polypeptide1A_23_P62351ARMCX60.0001590.759.75E−070.0258armadillo repeat containing, X-linked 61A_23_P164047MMD0.0001750.81.52E−050.00201monocyte to macrophagedifferentiation-associated1A_23_P33683MARCH20.0002040.84.40E−050.00095membrane-associated ring finger(C3HC4) 21A_23_P154070TUBA10.0002100.783.37E−050.00131tubulin, alpha 11A_23_P502710GAS2L10.0002400.846.95E−050.000831growth arrest-specific 2 like 11A_32_P122754MGC173370.0002590.799.01E−060.00746chromosome 9 ORF 301A_23_P128598TUBA20.0002810.813.96E−050.00199tubulin, alpha 21A_23_P40470H2BFS0.0002940.774.14E−060.0209H2B histone family, member S1A_23_P2114FLJ206250.0003510.773.21E−050.00383hypothetical protein FLJ206251A_23_P366254SLC10A30.0004310.783.30E−050.00562solute carrier family 10 member 31A_23_P103981HIST2H2AA0.0004470.729.79E−060.0204histone 2, H2aa1A_23_P39684TLK10.0004850.747.46E−060.0315tousled-like kinase 11A_24_P9296500.0004990.683.37E−060.07391A_23_P151120ACRBP0.0005120.794.40E−050.00595acrosin binding protein1A_32_P619360.0005150.782.89E−050.00919clone IMAGE: 51733891A_24_P124957RAB11A0.0005750.786.15E−050.00537RAB11A1A_23_P156550TREML10.0006150.778.03E−050.00471triggering receptor expressed onmyeloid cells-like 11A_32_P4814FAM11A0.0006250.786.66E−050.00586family with sequence similarity11 member A1A_23_P116264NRGN0.0006640.763.83E−050.0115neurogranin1A_24_P10657CTL20.0007250.750.0001040.00505solute carrier family 44 member 21A_23_P330611WASPIP0.0008280.750.000170.00403Wiskott-Aldrich syndromeprotein interacting protein1A_23_P165840ODC10.0008400.726.98E−050.0101ornithine decarboxylase 11A_23_P40718PARVB0.0009120.681.68E−050.0495parvin, beta1A_23_P166677MFSD10.0009960.760.0001260.00788major facilitator superfamilydomain containing 11A_23_P54488BG10.0010280.755.34E−050.0198acyl-CoA synthetase bubblegumfamily member 11A_23_P78209MAFG0.0011080.740.0003070.004v-maf musculoaponeuroticfibrosarcoma oncogene homolog G1A_24_P74371PPGB0.0013550.740.0001570.0117protective protein for beta-galactosidase1A_24_P259490ARF10.0014650.69.37E−060.229ADP-ribosylation factor 11A_23_P118038NUTF20.0014690.730.0007850.00275nuclear transport factor 21A_23_P98900FLJ224710.0014710.730.0001250.0173limkain beta 21A_23_P31177FLJ110000.0014910.676.11E−050.0364hypothetical protein FLJ110001A_24_P74374PPGB0.0015000.732.25E−050.1protective protein for beta-galactosidase1A_24_P8259420.0015570.641.05E−050.231FLJ10934 fis1A_24_P107695ACTN10.0016900.760.001210.00236actinin, alpha 11A_23_P147098MTPN0.0017880.70.0002830.0113myotrophin1A_32_P194848TAGLN20.0026500.655.53E−050.127transgelin 21A_32_P751410.0030410.740.001080.008561A_23_P76364CD90.0030600.730.0004660.0201CD9 antigen1A_23_P255444DAPP10.0031140.613.20E−050.303dual adaptor of phosphotyrosineand 3-phosphoinositides1A_23_P102109TUBA40.0031440.730.0003080.0321tubulin, alpha 41A_24_P55465MTPN0.0035190.647.94E−050.156myotrophin1A_23_P502224DIA10.0035900.680.0006510.0198cytochrome b5 reductase 31A_23_P100469TXNL4B0.0037170.670.0002250.0614thioredoxin-like 4B1A_23_P138717RGS100.0039660.690.00110.0143regulator of G-protein signalling101A_23_P162559SPPL30.0046190.70.009830.00217signal peptide peptidase 31A_24_P137897IFRD10.0046580.660.0007160.0303interferon-related developmentalregulator 11A_24_P147263USP310.0050960.690.0007570.0343ubiquitin specific peptidase 311A_23_P361773CCND30.0059430.650.0006780.0521cyclin D31A_23_P305711RYBP0.0060420.680.0007230.0505RING1 and YY1 binding protein1A_32_P192545LOC1589310.0064420.660.0004970.0835transcription elongation factor A(SII)-like 61A_23_P141394WIPI490.0075880.610.02360.00244WD repeat domain,phosphoinositide interacting 11A_23_P341392MGC321240.0077970.670.00320.019hypothetical protein MGC321241A_23_P138881ACTN30.0078080.640.001270.048actinin, alpha 31A_23_P434442TCEAL30.0082600.670.000710.0961transcription elongation factor A(SII)-like 31A_23_P302550RGS180.0090140.740.005490.0148regulator of G-protein signalling181A_23_P30799HIST1H3F0.0091380.620.01130.00739histone 1, H3f1A_24_P80135PTPN180.0096580.620.0006910.135protein tyrosine phosphatase,non-receptor type 181A_24_P319736MEIS10.0100570.670.01650.00613myeloid ecotropic viralintegration site 1 homolog1A_23_P6321CLDN50.0117390.620.04320.00319claudin 51A_24_P186414TEX270.0136160.50.0005740.323zinc finger, AN1-type domain 31A_23_P215479CYLN20.0151690.580.001330.173cytoplasmic linker 21A_23_P360379EGLN30.0151790.570.001440.16egl nine homolog 31A_23_P95470CD1510.0164680.540.0005650.48CD151 antigen1A_24_P29733PFTK10.0211280.590.00360.124PFTAIRE protein kinase 11A_24_P333525RABGAP1L0.0216670.560.001440.326RAB GTPase activating protein1-like1A_23_P200325RABGAP1L0.0244290.550.00160.373RAB GTPase activating protein1-like1A_23_P502915WDR10.0248800.520.002590.239WD repeat domain 11A_23_P132226TPST20.0257880.560.00250.266tyrosylprotein sulfotransferase 21A_24_P922357LOC1289770.0270260.530.006580.111hypothetical protein LOC1289771A_23_P141688RAB310.0290410.580.04510.0187RAB311A_24_P134834DKFZp547E0520.0302310.490.00370.247hypothetical protein LOC842361A_24_P1793390.0306260.550.01110.0845humanin1A_23_P126135MFN20.0308280.540.03370.0282mitofusin 21A_24_P244916SERF20.0325810.470.003860.275small EDRK-rich factor 21A_23_P141974TPM40.0327900.540.00480.224tropomyosin 41A_23_P251825IFRD10.0328500.60.01880.0574interferon-related developmentalregulator 11A_23_P48175MGC55760.0358750.570.0550.0234transmembrane protein 106C1A_32_P1191650.0375320.570.01560.09031A_32_P427800.0417110.510.05780.03011A_23_P502913WDR10.0418570.540.020.0876WD repeat domain 11A_23_P27207SCGB1C10.0454250.490.1140.0181secretoglobin, family 1C member 11A_24_P403303PHF20L10.0463850.50.01320.163PHD finger protein 20-like 11A_24_P3160590.0512200.550.0450.05831A_23_P147199ZNF2710.0669150.40.01160.386zinc finger protein 2711A_32_P559790.0806770.450.05760.1136-pyruvoyltetrahydropterinsynthase1A_32_P163089LOC3878820.0954360.450.0120.759hypothetical protein1A_23_P2661RAP1B0.3829990.260.3840.382RAP1B1A_23_P154294MGC130050.440704−0.120.1950.996FLJ44010 fis2A_23_P2103300.0000310.882.72E−053.45E−05CS0DL009YB17 of B cells2A_24_P365901MGC508440.0000380.872.91E−060.000507tetraspanin 332A_24_P226322SH3BGRL20.0000480.856.74E−060.000337SH3 domain binding glutamicacid-rich protein like 22A_23_P152906ALOX120.0000710.868.47E−060.000587arachidonate 12-lipoxygenase2A_24_P148321HIST2H2BE0.0000710.776.88E−070.00738histone 2, H2be2A_23_P256205ABLIM30.0000740.838.81E−060.000629actin binding LIM protein familymember 32A_24_P209171SH3BGRL20.0000810.841.39E−050.000471SH3 domain binding glutamicacid-rich protein like 22A_23_P390006PCSK60.0001340.775.81E−060.00309proprotein convertasesubtilisin/kexin type 62A_23_P129221FAH0.0001370.71.83E−070.102fumarylacetoacetate hydrolase2A_23_P430818HSPC1590.0001520.820.0001140.000202HSPC159 protein2A_24_P218905NET-50.0001620.84.72E−050.000556tetraspanin 92A_32_P1454770.0002180.845.04E−050.000947BX3502562A_24_P2901880.0002330.810.0001520.0003572A_24_P7067520.0002440.834.95E−050.00122A_23_P143720GRAP20.0002760.774.55E−060.0168GRB2-related adaptor protein 22A_23_P77971ITGA2B0.0003150.782.23E−050.00445integrin, alpha 2b2A_23_P212436CTDSPL0.0003940.823.79E−050.0041carboxy-terminal domain, RNApolymerase II, polypeptide A2A_23_P152926GP1BA0.0003980.850.001440.00011glycoprotein Ib, alphapolypeptide2A_24_P189997PCSK60.0004030.694.78E−060.034proprotein convertasesubtilisin/kexin type 62A_23_P38519ITGB30.0004080.810.0001110.0015integrin, beta 32A_24_P64167PTGS10.0004110.783.30E−050.00513prostaglandin-endoperoxidesynthase 12A_24_P318656ITGB30.0004220.830.0003230.000552integrin, beta 32A_23_P24140.0004470.731.24E−050.0161PSEC0021 fis2A_23_P216966PTGS10.0004550.782.95E−050.00703prostaglandin-endoperoxidesynthase 12A_32_P177430.0004700.790.0001740.001272A_23_P79978SLC24A30.0005810.784.70E−050.00718solute carrier family 24 member 32A_23_P43810LTBP10.0005820.80.0001850.00183latent transforming growth factorbeta binding protein 12A_23_P6034TUBB10.0005830.84.83E−050.00704tubulin, beta 12A_24_P176079WASF30.0006200.760.005427.10E−05WAS protein family member 32A_23_P202823CTTN0.0006550.760.002450.000175cortactin2A_23_P210358LIMS10.0006900.790.0002290.00208LIM and senescent cell antigen-like domains 12A_24_P929003ITGB30.0007150.80.0001330.00384integrin, beta 32A_23_P389118TMEM16F0.0007660.721.73E−050.0339DKFZp313M07202A_23_P106042CKLFSF50.0007690.780.0002110.0028CKLF-like MARVELtransmembrane domaincontaining 52A_24_P160104TUBA80.0007970.760.0002320.00274tubulin, alpha 82A_23_P207507ABCC30.0008090.770.0001220.00536ATP-binding cassette, sub-familyC member 32A_23_P102731SMOX0.0008190.758.77E−050.00764spermine oxidase2A_32_P1376040.0008380.810.0003480.00202clone IMAGE: 38692762A_23_P104624KIAA08300.0009100.750.0002190.00378KIAA0830 protein, partial cds2A_23_P3592770.0009650.766.65E−050.014ELOVL family member 72A_23_P151133NET-50.0010030.750.0001150.00875tetraspanin 92A_23_P105957ACTN10.0010080.780.0008990.00113actinin, alpha 12A_23_P17095TFPI0.0010310.720.001360.000782tissue factor pathway inhibitor2A_23_P25974TTC7B0.0010710.810.0006340.00181tetratricopeptide repeat domain7B2A_32_P168342C6orf250.0011130.780.000160.00774FLJ35073 fis2A_23_P215913CLU0.0011470.80.0001290.0102clusterin2A_23_P416581GNAZ0.0011550.80.0003350.00398guanine nucleotide bindingprotein2A_24_P122337SYTL40.0011750.740.0001070.0129synaptotagmin-like 42A_23_P166633ITGB50.0012070.790.0007290.002integrin, beta 52A_24_P185186LOC2011910.0012570.712.60E−050.0608sterile alpha motif domaincontaining 142A_24_P3333720.0012950.720.0004920.00341FLJ35984 fis2A_23_P217998JAM30.0013100.735.66E−050.0303junctional adhesion molecule 32A_23_P81930C6orf250.0013570.617.97E−060.231chromosome 6 ORF 252A_23_P152160SNN0.0014280.750.0002940.00694stannin2A_23_P109974RAB6B0.0014750.710.00350.000622RAB6B2A_23_P45524NGFRAP10.0017610.750.0008860.0035nerve growth factor receptorassociated protein 12A_23_P7642SPARC0.0017690.740.000310.0101secreted protein, acidic, cysteine-rich2A_23_P73457RUFY10.0018500.770.00210.00163RUN and FYVE domaincontaining 12A_32_P1364500.0019050.582.20E−050.165AF220206 Nedd4 WW domain-binding protein 22A_23_P17724SEP50.0020390.62.52E−050.165septin 52A_23_P42975PRKAR2B0.0020980.760.0002840.0155protein kinase, cAMP-dependent,regulatory, type II, beta2A_23_P19987IMP-30.0023570.750.0001570.0354IGF-II mRNA-binding protein 32A_32_P162250ARHGAP180.0025100.750.003480.00181Rho GTPase activating protein 182A_24_P251534CTDSPL0.0027660.780.001260.00607carboxy-terminal domain, RNApolymerase II, polypeptide A2A_23_P3915860.0028330.730.0008180.00981tropomyosin 1 transcript variant 32A_24_P319923MYLK0.0030910.720.002810.0034myosin, light polypeptide kinase2A_24_P13190ESAM0.0032220.720.001450.00716endothelial cell adhesionmolecule2A_23_P105562VWF0.0032470.680.0001720.0613von Willebrand factor2A_23_P111701GNG110.0032490.670.005360.00197guanine nucleotide bindingprotein, gamma 112A_24_P254850KIAA04200.0037070.740.0003870.0355KIAA0420 mRNA2A_24_P79403PF40.0038370.690.001620.00909platelet factor 42A_23_P121596PPBP0.0041650.680.000640.0271pro-platelet basic protein2A_23_P143817MYLK0.0046430.70.00550.00392myosin, light polypeptide kinase2A_23_P217428ARHGAP60.0049340.720.004120.00591Rho GTPase activating protein 62A_23_P146584MGC173370.0050010.740.004210.00594chromosome 9 ORF 302A_23_P149992PDLIM10.0054160.577.66E−050.383PDZ and LIM domain 12A_23_P500844PDE5A0.0054270.670.0004250.0693phosphodiesterase 5A, cGMP-specific2A_23_P99906HOMER20.0060030.690.0006150.0586homer homolog 22A_24_P921366CALD10.0063140.680.0002910.137caldesmon 12A_23_P125233CNN10.0063680.630.0005240.0774calponin 1, basic2A_23_P8906LRP120.0068270.70.002950.0158low density lipoprotein-relatedprotein 122A_32_P140139F13A10.0076550.670.0020.0293coagulation factor XIII, A1polypeptide2A_23_P360804CPNE50.0081400.60.001620.0409copine V2A_24_P188071TUBA60.0121960.590.0370.00402tubulin, alpha 62A_23_P137697SELP0.0136490.630.007860.0237selectin P2A_24_P8926120.0162100.640.004570.0575DKFZp313A1372A_23_P48212CLEC1B0.0201060.660.006070.0666C-type lectin domain family 1,member B2A_23_P58396PDGFC0.0395970.590.05640.0278platelet derived growth factor C2A_23_P2095270.0404570.560.01860.088A31642 villin2A_23_P168556STX1A0.0505710.530.04440.0576syntaxin 1A2A_32_P18723DKFZp762C11120.0513550.520.02980.0885FLJ38153 fis2A_23_P52207BAMBI0.0570180.460.007910.411BMP and activin membrane-bound inhibitor homolog2A_23_P431388SPOCD10.0686620.520.09410.0501SPOC domain containing 12A_23_P371266DNM30.0725000.490.06880.0764dynamin 32A_32_P1791380.0879640.430.1860.0416clone IMAGE: 53021583A_23_P111267SH3BGRL20.0001450.845.79E−050.000364SH3 domain binding glutamicacid-rich protein like 23A_23_P219045HIST1H3D0.0001580.764.14E−060.00601histone 1, H3d3A_24_P3152560.0001950.642.00E−070.193A_23_P91423C20orf1120.0002060.772.21E−050.00192chromosome 20 ORF 1123A_23_P149545HIST2H2BE0.0002340.757.28E−060.00752histone 2, H2be3A_23_P84448TUBA40.0003290.658.46E−070.128tubulin, alpha 43A_23_P405295LCE3C0.0003330.88.47E−050.00131late cornified envelope 3C3A_23_P152909ALOX120.0003750.790.0001080.0013arachidonate 12-lipoxygenase3A_23_P210939ITGB4BP0.0004740.760.000220.00102integrin beta 4 binding protein3A_23_P4944CALM30.0004970.642.13E−060.116calmodulin 33A_32_P221799HIST1H2AM0.0005110.790.0001330.00196histone 1, H2am3A_23_P436138MAX0.0006090.665.84E−060.0636MYC associated factor X3A_24_P180680LAPTM4B0.0007320.722.30E−050.0233lysosomal associated proteintransmembrane 4 beta3A_24_P753476LOC3405080.0007580.80.0001110.00518LOC3405083A_24_P65373ITGA2B0.0009880.645.64E−060.173integrin, alpha 2b3A_24_P918032LOC3390050.0010070.749.75E−050.0104LOC3390053A_23_P160546FLJ112800.0012390.760.0009190.00167family with sequence similarity63, member A3A_23_P41280PAICS0.0014360.720.005620.000367phosphoribosylaminoimidazolecarboxylase3A_24_P258633TUBB30.0014630.080.0007230.00296tubulin, beta 33A_24_P308506CML20.0015030.730.007530.0003putative N-acetyltransferaseCamello 23A_23_P206212THBS10.0015650.720.0001230.0199thrombospondin 13A_24_P382637GTPBP50.0017410.70.005490.000552GTP binding protein 53A_32_P33850.0017980.720.003580.000903CS0DI060YD223A_23_P156708TNXB0.0019690.680.0439.02E−05tenascin XB3A_23_P74138TAGLN20.0020920.665.86E−050.0747transgelin 23A_23_P215735ST70.0020940.61.72E−050.255suppression of tumorigenicity 73A_23_P113701PDGFA0.0022360.740.0003650.0137platelet-derived growth factor3A_23_P121564GUCY1B30.0026400.622.63E−050.265guanylate cyclase 1, soluble, beta 33A_24_P189533KIAA08300.0027960.720.0007740.0101KIAA08303A_32_P897090.0028780.720.0009910.00836tropomyosin 13A_23_P15647NLK0.0030350.645.98E−050.154nemo-like kinase3A_23_P24616CSE-C0.0030700.670.000160.0589sialic acid acetylesterase3A_23_P73239NCKAP10.0035800.670.0005430.0236NCK-associated protein 13A_23_P3946NT5M0.0037330.482.66E−050.5245′,3′-nucleotidase, mitochondrial3A_23_P19624BMP60.0040050.690.002250.00713bone morphogenetic protein 63A_24_P9267090.0040130.610.001890.008523A_23_P90407CASP140.0041290.70.00840.00203caspase 143A_23_P167096VEGFC0.0042340.670.0001990.0901vascular endothelial growthfactor C3A_23_P421843LOC2011910.0044490.574.90E−050.404sterile alpha motif domaincontaining 143A_23_P501831C5orf40.0048480.650.0005690.0413chromosome 5 ORF 43A_23_P417942FNBP1L0.0051060.630.0003490.0747formin binding protein 1-like3A_23_P307525ANKRD90.0051360.526.28E−050.42ankyrin repeat domain 93A_32_P922120.0054010.620.0002040.143IMAGE: 32717273A_23_P156284DBN10.0055440.670.0004560.0674drebrin 13A_23_P18539MMRN10.0057810.570.0002050.163multimerin 13A_24_P38387NDRG10.005842−0.60.005550.00615N-myc downstream regulatedgene 13A_23_P155979EGF0.0063970.670.001320.031epidermal growth factor (beta-urogastrone)3A_23_P401361PITPNM20.0073480.620.000360.15phosphatidylinositol transferprotein, membrane-associated 23A_24_P385313PTPRF0.0083400.640.007780.00894protein tyrosine phosphatase,receptor type, F3A_23_P141055TGFB1I10.0086070.489.89E−050.749transforming growth factor beta 1induced transcript 13A_24_P2042570.0119180.650.004190.03393A_23_P369899RIS10.0136870.450.0004090.458Ras-induced senescence 13A_24_P167654SLC8A30.0162780.570.002070.128solute carrier family 8 member 33A_24_P4059810.0183860.580.003130.108CS0DD001YH153A_23_P4318530.0184310.60.01140.0298A-COL042173A_23_P367043MGC264840.0203590.610.01320.0314CDC14 cell division cycle 14homolog C3A_23_P135499CLIC40.0264100.540.00940.0742chloride intracellular channel 43A_24_P324730.0265870.510.001980.357ELOVL family member 73A_23_P81934C6orf250.0314980.530.004390.226chromosome 6 ORF 253A_24_P414999LAPTM4B0.0322560.550.004390.237lysosomal associated proteintransmembrane 4 beta3A_23_P207414MGC27440.0347820.360.001430.846alanyl-tRNA synthetase domaincontaining 13A_32_P1414370.0391000.510.01040.147FKSG733A_32_P592620.0570920.520.02050.159IMAGE: 31040773A_23_P61945MITF0.0619350.510.0280.137microphthalmia-associatedtranscription factor3A_23_P16866VIL10.0660770.480.1140.0383villin 13A_23_P127642ARHGEF120.0729000.480.01460.364Rho guanine nucleotide exchangefactor 123A_24_P7131850.0751920.460.02870.197IMAGE: 42715223A_23_P69573GUCY1A30.1061530.380.5720.0197guanylate cyclase 1, soluble,alpha 33A_23_P257871DAB20.1377680.430.1460.13disabled homolog 2, mitogen-responsive phosphoprotein3A_24_P331882KIAA12110.1886150.350.06040.589DKFZp434F1173A_23_P154526GRB140.2052320.320.4320.0975growth factor receptor-boundprotein 143A_23_P45304XK0.2759350.20.1620.47Kell blood group precursor3A_23_P104493PAPSS20.5883940.120.3490.9923′-phosphoadenosine 5′-phosphosulfate synthase 24A_24_P143440DNCL2A0.0000040.852.56E−080.000546dynein, light chain, roadblock-type 14A_23_P208788C19orf330.0000400.850.0001828.84E−06chromosome 19 ORF 334A_24_P68631HIST2H2AB0.0000480.831.11E−050.000204histone 2, H2ab4A_23_P120364C20orf1490.0000600.89.45E−070.00385chromosome 20 ORF 1494A_23_P202029SPFH10.0001050.810.0001070.000104SPFH domain family, member 14A_24_P287075MAP4K20.0001110.84.13E−060.00296mitogen-activated protein kinasekinase kinase kinase 24A_23_P149301HIST3H2A0.0001110.831.06E−050.00117histone 3, H2a4A_24_P69210.0002880.80.0001930.00043LOC541471 protein4A_24_P1358010.0003140.87.96E−050.00124CS0DF024YI144A_24_P45767FLJ218390.0004170.80.000740.000235FLJ218394A_23_P42375RAB320.0005260.769.01E−050.00307RAB324A_23_P354705ST8SIA10.000559−0.660.0001850.00169ST8 alpha-N-acetyl-neuraminidealpha-2,8-sialyltransferase 14A_23_P407565CX3CR10.000593−0.640.008833.98E−05chemokine receptor 14A_23_P59045HIST1H2AE0.0006070.770.0001230.003histone 1, H2ae4A_24_P9119600.0006420.777.30E−050.00565IMAGE: 16997324A_32_P1849370.0006550.770.0005850.000733BU6789414A_23_P138117CAMTA10.0006780.760.0001830.00251calmodulin binding transcriptionactivator 14A_32_P31132-Mar0.0008020.790.0001190.0054membrane-associated ring finger(C3HC4) 24A_32_P54137UQCRH0.0010630.740.005760.000196ubiquinol-cytochrome c reductasehinge protein4A_23_P396626AP1GBP10.001214−0.660.002010.000733AP1 gamma subunit bindingprotein 14A_24_P227927IL21R0.001223−0.620.005710.000262interleukin 21 receptor4A_23_P110167MGST20.0013710.780.001080.00174microsomal glutathione S-transferase 24A_23_P120933ATF40.0014160.720.001050.00191activating transcription factor 44A_23_P55706RELB0.001465−0.640.005110.00042v-rel reticuloendotheliosis viraloncogene homolog B4A_24_P323835H3F3A0.0015730.710.000290.00853H3 histone, family 3A4A_24_P273143MGC46770.0016560.740.002940.000933hypothetical protein MGC46774A_23_P111037HIST1H3A0.0017030.740.001070.00271histone 1, H3a4A_24_P223384HIST1H2AB0.0018430.720.0007240.00469histone 1, H2ab4A_23_P1322850.0019830.740.002570.00153mercaptopyruvatesulfurtransferase4A_23_P52101NQO3A20.0023280.720.004170.0013cytochrome b5 reductase 14A_24_P6087900.0026190.70.003810.00184A_24_P122732SLC41A10.002673−0.630.0009010.00793solute carrier family 41, member 14A_32_P1323170.003079−0.630.001390.006824A_23_P218817CPT1B0.003379−0.630.00660.00173carnitine palmitoyltransferase 1B4A_32_P1321690.0034190.70.001910.006124A_24_P102769UQCRH0.0037120.690.005010.00275ubiquinol-cytochrome c reductasehinge protein4A_23_P148410FTHL170.0052960.516.36E−050.441ferritin, heavy polypeptide-like174A_23_P214330SERPINB10.0056420.640.005840.00545serpin peptidase inhibitor, cladeB, member 14A_32_P945210.0056880.680.0007950.04074A_32_P593020.0060400.610.005190.00703IMAGE: 62540314A_23_P121082GBE10.0067060.670.00480.00937glucan branching enzyme 14A_24_P79340.0074720.620.008710.00641LOC3917694A_23_P250671GPX10.0445710.460.04130.0481glutathione peroxidase 14A_23_P69218LOC558310.0476630.480.02920.0778transmembrane protein 1114A_24_P9136290.0836240.430.02220.3155A_23_P122007LOC903550.0000130.785.58E−070.000291hypothetical gene supported byAF0381825A_23_P210060MGC130570.0000320.882.36E−060.000436DKFZp686I152105A_23_P138417RSU10.0000350.81.88E−070.0066Ras suppressor protein 15A_23_P350591CXorf200.0000380.837.91E−070.00181chromosome X ORF 205A_23_P114275PGRMC10.0000520.752.08E−070.0131progesterone receptor membranecomponent 15A_24_P362540DDEF20.0001370.841.92E−050.000971development and differentiationenhancing factor 25A_23_P167983HIST1H2AC0.0001930.785.63E−060.0066histone 1, H2ac5A_23_P103070YWHAH0.0002230.89.09E−060.00546tyrosine 3-monooxygenase5A_24_P273666GNAS0.0002380.744.23E−060.0134GNAS complex locus5A_23_P333484HIST1H3H0.0002750.790.0001390.000545histone 1, H3h5A_23_P4072030.0002860.664.61E−070.177FLJ42816 fis5A_23_P414273NID670.0003140.820.0008410.000117MSTP1505A_23_P102391SLC40A10.0003250.761.69E−050.00624solute carrier family 40 member 15A_23_P2060180.0004030.761.69E−050.00963tropomyosin 15A_23_P72668SDPR0.0005670.788.47E−050.00379serum deprivation response5A_24_P228550TUBB10.0006060.820.0001810.00203tubulin, beta 15A_23_P107612RAB27B0.0007440.80.0001510.00367RAB27B5A_23_P77145RAB11A0.0008930.760.0001140.007RAB11A5A_23_P502797WDFY10.000895−0.010.003930.000204WD repeat and FYVE domaincontaining 15A_23_P211910PLOD20.0009430.770.0005320.00167procollagen-lysine, 2-oxoglutarate 5-dioxygenase 25A_24_P44462TPM10.0009910.734.86E−050.0202tropomyosin 15A_32_P1259170.0010680.790.0001440.00792BF2388435A_23_P157128SCAP20.0012590.720.000130.0122src family associatedphosphoprotein 25A_32_P168349C6orf250.0013220.749.06E−050.0193FLJ35073 fis5A_23_P216679CDC14B0.0013720.661.61E−050.117CDC14 cell division cycle 14homolog B5A_23_P63371TAL10.0014780.80.002080.00105T-cell acute lymphocyticleukemia 15A_23_P12884GRK50.0015450.736.47E−050.0369G protein-coupled receptor kinase 55A_23_P126836TNFSF40.0016110.780.0002380.0109tumor necrosis factorsuperfamily, member 45A_23_P23221GADD45A0.0020360.70.0009820.00422growth arrest and DNA-damage-inducible, alpha5A_23_P115608ARHGAP210.0021710.720.0001720.0274Rho GTPase activating protein 215A_24_P135444AMFR0.0021790.670.0001610.0295autocrine motility factor receptor5A_24_P118376UNQ93660.0027900.613.09E−050.252carcinoembryonic antigen-relatedcell adhesion molecule 205A_23_P124476CLCN30.0034500.630.0001540.0773chloride channel 35A_32_P357510.0035930.690.0002640.04895A_32_P1035580.0036530.740.004980.00268FLJ37480 fis5A_23_P334123CDA080.0043560.587.53E−050.252T-cell immunomodulatory protein5A_23_P143902P2RY120.0050460.690.000380.067purinergic receptor P2Y5A_23_P1366930.0052530.680.02190.00126DKFZp686D05215A_23_P33947EFHC20.0071070.660.01730.00292EF-hand domain containing 25A_23_P139486CDK2AP10.0072180.620.0005490.0949CDK2-associated protein 15A_23_P217611ARMCX30.0072950.670.0006490.082armadillo repeat containing, X-linked 35A_23_P86424NCOA40.0078910.520.0002650.235nuclear receptor coactivator 45A_23_P115375H3/o0.0079560.620.030.00211histone H3/o5A_23_P91900SMC4L10.0084770.040.007710.00932SMC4 structural maintenance ofchromosomes 4-like 15A_23_P422083DKFZp762O0760.0095000.60.0005860.154transmembrane protein 55A5A_23_P69226LOC558310.0098510.640.002390.0406transmembrane protein 1115A_23_P59547NT5C30.0108600.620.008190.01445′-nucleotidase, cytosolic III5A_24_P5006210.0123770.650.003820.0401FLJ23711 fis5A_24_P26897INPP5A0.0127550.620.00330.0493inositol polyphosphate-5-phosphatase5A_23_P11025ZNF1850.0131070.640.0008720.197zinc finger protein 1855A_24_P349560EIF4E0.0140110.540.0005180.379eukaryotic translation initiationfactor 4E5A_24_P941699PCGF50.0157160.560.002330.106polycomb group ring finger 55A_24_P147927EFHC20.0164700.540.005080.0534EF-hand domain containing 25A_23_P8763PTPN120.0228230.60.003830.136protein tyrosine phosphatase,non-receptor type 125A_24_P81947CORO1C0.0242310.580.006560.0895coronin, actin binding protein, 1C5A_23_P371239CMIP0.0246700.460.00170.358c-Maf-inducing protein5A_23_P135494CLIC40.0275430.580.008750.0867chloride intracellular channel 45A_23_P72643ADAM90.0298990.540.006340.141metallopeptidase domain 95A_24_P5038660.0492920.420.009680.251CS0DL005YE025A_24_P23411ARMCX30.0528660.520.01370.204armadillo repeat containing, X-linked 35A_24_P633902ZNF3640.0548950.520.01230.245zinc finger protein 3645A_32_P96134KIAA08770.0654650.440.006980.614KIAA08775A_23_P201376SSX2IP0.0717890.440.009010.572synovial sarcoma, X breakpoint 2interacting protein5A_32_P61720.0759080.460.2150.0268IMAGE: 52868435A_24_P27373PLDN0.1017390.420.04770.217pallidin homolog5A_23_P96041FLJ226790.1129880.40.02020.632FLJ226795A_32_P393840.1616550.090.3060.0854IMAGE: 48234166A_23_P145965TPST10.0000100.93.05E−052.97E−06tyrosylprotein sulfotransferase 16A_23_P33723CD1630.0000770.85.09E−050.000116CD163 antigen6A_24_P38081FKBP50.0001380.863.82E−050.000498FK506 binding protein 56A_23_P111206FKBP50.0002280.832.71E−050.00191FK506 binding protein 56A_23_P121602SAP300.0002440.80.0004840.000123sin3-associated polypeptide6A_23_P328729KLHL80.0002730.820.0001110.000672kelch-like 86A_23_P104804ZBTB160.0006070.790.003899.47E−05zinc finger and BTB domaincontaining 166A_23_P99442FLT30.0007800.790.0001110.00548fms-related tyrosine kinase 36A_32_P806841ARL4A0.0014430.68.40E−050.0248ADP-ribosylation factor-like 4A6A_32_P223985LOC3887520.0019170.740.001570.00234LOC3887526A_24_P322150.0022880.577.92E−050.0661ADP-ribosylation factor-like 4B6A_23_P145761ARL4A0.0022890.550.0001390.0377ADP-ribosylation factor-like 4A6A_23_P53838IRS20.0025150.740.001110.0057insulin receptor substrate 26A_24_P213296dJ341D10.10.0034410.750.007640.00155dJ341D10.16A_23_P415401KLF90.0508500.460.01170.221Kruppel-like factor 97A_23_P113212TMEM45A0.0000270.824.53E−060.000165transmembrane protein 45A7A_32_P1140200.0000600.858.32E−054.26E−05T328247A_32_P291400.0000610.850.001822.04E−06AA3446327A_32_P1309680.0001370.84.08E−050.000461IMAGE: 48262407A_23_P57658HRASLS0.0001750.741.43E−060.0214HRAS-like suppressor7A_23_P381714CA130.0002590.778.82E−050.000762carbonic anhydrase XIII7A_32_P1314490.0002870.743.48E−060.02377A_23_P151662MAX0.0003100.813.69E−050.00261MYC associated factor X7A_23_P17130MGC130570.0003190.87.85E−050.0013hypothetical protein MGC130577A_24_P76675MFAP3L0.0003360.847.72E−050.00146microfibrillar-associated protein3-like7A_23_P331253XPNPEP10.0004790.762.30E−050.00998X-prolyl aminopeptidase 17A_24_P394510HIST1H2AJ0.0004970.721.75E−050.0141histone 1, H2aj7A_23_P200001NEXN0.0005870.760.0001040.00331nexilin7A_32_P387450.0006450.80.0006860.0006067A_24_P409971NEXN0.0007840.785.04E−050.0122nexilin7A_24_P363615MTPN0.0008580.720.0001020.00721myotrophin7A_32_P1961420.0009440.790.002240.0003987A_32_P808KIAA14580.0009990.744.66E−050.0214KIAA14587A_32_P790410.0013490.680.03465.26E−05IMAGE: 61792617A_23_P217938SPHAR0.0016570.730.000280.0098S-phase response7A_23_P132619OXTR0.0018590.730.0008980.00385oxytocin receptor7A_24_P4538190.0020840.710.001180.00368IMAGE: 303309557A_23_P363344TPM10.0023460.660.0001730.0318tropomyosin 17A_23_P365685LIMS30.0023800.770.0007590.00746LIM and senescent cell antigen-like domains 37A_24_P148094LEPROT0.0024160.70.0001110.0526leptin receptor overlappingtranscript7A_23_P131825TNNC20.0025680.740.008220.000802troponin C type 27A_23_P39202C19orf330.0028790.70.000860.00964chromosome 19 ORF337A_23_P16733RALB0.0031960.650.06190.000165v-ral simian leukemia viraloncogene homolog B7A_23_P160336LEFTY10.0037380.740.001020.0137left-right determination factor 17A_32_P1179080.0041570.640.0001630.1067A_24_P5146780.0047370.690.002910.007717A_23_P1126LEPROT0.0047840.660.00210.0109leptin receptor overlappingtranscript7A_23_P160582HT0360.0057850.640.006510.00514hydroxypyruvate isomerasehomolog7A_32_P278780.0065560.670.003070.014AA3996567A_23_P93282HIST1H3J0.0072060.650.007750.0067histone 1, H3j7A_24_P5708060.0084300.630.003450.0206IMAGE: 48144377A_32_P805320.0089280.489.84E−050.81BF7339087A_24_P35478PARD30.0122400.60.001650.0908par-3 partitioning defective 3homolog7A_23_P38876LIPE0.0124550.320.7680.000202lipase, hormone-sensitive7A_23_P89902RTN20.0131120.50.000330.521reticulon 27A_24_P8798950.0137370.620.002930.0644IMAGE: 38836597A_24_P231104LEPR0.0154910.560.001420.169leptin receptor7A_24_P5242620.0197900.580.06110.00641Q80YT07A_23_P38106SPHK10.0231660.450.0009050.593sphingosine kinase 17A_23_P137173TMSNB0.0233970.530.002380.23thymosin-like 87A_32_P25639BET3L0.0391220.510.1980.00773FLJ11180 fis7A_23_P426663MITF0.0449960.540.007030.288microphthalmia-associatedtranscription factor7A_23_P169756HIPK20.0456150.390.5490.00379homeodomain interacting proteinkinase 27A_23_P92025CIDEC0.0599780.340.003630.991cell death-inducing DFFA-likeeffector c7A_32_P1812970.0611800.440.4880.00767CS0DK012YG127A_23_P377214FLJ323840.0653830.480.02170.197hexamthylene bis-acetamideinducible 27A_32_P44330.0700700.460.0980.0501BU6024857A_23_P213050HPGD0.0890010.470.04830.164hydroxyprostaglandindehydrogenase 15-(NAD)7A_23_P328740LOC930820.1243780.420.08740.177BC0123177A_24_P347447DAAM10.1343650.390.1770.102dishevelled associated activatorof morphogenesis 17A_23_P54116DAAM10.1549320.390.09340.257dishevelled associated activatorof morphogenesis 17A_23_P65674TMOD30.2575640.320.2730.243tropomodulin 37A_32_P2251350.4830720.230.2590.901IMAGE: 52778598A_23_P46369RAB130.0000310.838.86E−060.000106RAB138A_23_P130961ELA20.0001320.841.56E−050.00111elastase 28A_23_P140384CTSG0.0002710.827.50E−050.000982cathepsin G8A_23_P86653PRG10.0003260.752.86E−050.00372proteoglycan 1, secretory granule8A_23_P141173MPO0.0006340.786.51E−050.00618myeloperoxidase8A_23_P167005GPR1600.0010610.720.009380.00012G protein-coupled receptor 1608A_23_P121716ANXA30.0013150.720.001440.0012annexin A38A_23_P326080DEFA40.0014670.70.0003330.00646defensin, alpha 4, corticostatin8A_24_P347378ALOX5AP0.0015410.710.001440.00165arachidonate 5-lipoxygenase-activating protein8A_23_P201193TSPAN20.0019210.710.0004890.00755tetraspanin 28A_23_P150903MLSTD10.0019620.710.001040.0037male sterility domain containing 18A_23_P131789BPI0.0028740.660.001130.00731bactericidal/permeability-increasing protein8A_23_P169437LCN20.0029060.670.003640.00232lipocalin 28A_23_P159952BEX10.0037540.660.00520.00271brain expressed, X-linked 18A_23_P69171SUCNR10.0045680.650.003460.00603succinate receptor 18A_23_P71981ERAL10.0090850.620.009070.0091Era G-protein-like 19A_24_P63019IL1R20.0000020.832.00E−063.01E−06interleukin 1 receptor, type II9A_23_P60627ALOX15B0.0000100.850.0001258.69E−07arachidonate 15-lipoxygenase,second type9A_23_P4036HT0080.0000150.891.98E−051.17E−05testis expressed sequence 29A_23_P117582JDP20.0000340.840.0002315.06E−06jun dimerization protein 29A_32_P224094ZNF1430.0000560.790.0005295.91E−06zinc finger protein 1439A_24_P202567ITPKC0.0000620.841.65E−050.000232inositol 1,4,5-trisphosphate 3-kinase C9A_23_P162668CPM0.0000820.80.0002962.25E−05carboxypeptidase M9A_23_P255104LHFPL20.0001010.792.49E−060.0041lipoma HMGIC fusion partner-like 29A_23_P155765HMGB20.0001130.822.57E−050.000497high-mobility group box 29A_23_P169529HRB0.0001390.85.28E−060.00365HIV-1 Rev binding protein9A_23_P1161950.0001620.810.004815.45E−06Q7PKG09A_23_P11201GPR340.0001670.844.19E−050.000666G protein-coupled receptor 349A_23_P388900SLC22A150.0002100.810.0001070.000414solute carrier family 22, member159A_24_P938352CPM0.0002690.810.0001250.000577carboxypeptidase M9A_23_P423864PHC20.0004020.766.02E−050.00268polyhomeotic-like 29A_23_P138725MARVELD10.0005640.790.0001750.00182MARVEL domain containing 19A_24_P269687TOR1A0.0005860.75.25E−050.00653torsin family 1, member A9A_24_P9131150.0008140.770.0001370.00484CS0DK002YE209A_23_P93562SESN10.0011840.740.0003560.00394sestrin 19A_23_P104798IL180.0011900.780.0003430.00413interleukin 189A_23_P8640GPR300.0015360.760.000680.00347G protein-coupled receptor 309A_24_P78531CLEC4E0.0021610.770.002510.00186C-type lectin domain family 4,member E9A_23_P215566AHR0.0024740.750.003580.00171aryl hydrocarbon receptor9A_23_P415021DKFZP586A05220.0032170.710.006390.00162DKFZP586A05229A_24_P154037IRS20.0036070.80.002150.00605insulin receptor substrate 29A_24_P750164LOC1514380.0043840.710.002230.00862\FLJ31315 fis9A_23_P98085PTEN0.0049270.680.002740.00886phosphatase and tensin homolog9A_24_P233995FLJ223900.0086450.690.00870.00859MOCO sulphurase C-terminaldomain containing 110A_24_P235266GRB100.0000440.82.84E−060.000697growth factor receptor-boundprotein 1010A_23_P122863GRB100.0002070.765.88E−060.0073growth factor receptor-boundprotein 1010A_24_P360674CDKN2B0.0020520.698.42E−050.05cyclin-dependent kinase inhibitor2B10A_24_P323084FLJ394210.0072510.680.001060.0496chromosome 17 ORF 5510A_23_P502470IL6ST0.0073300.670.007270.00739interleukin 6 signal transducer












TABLE 2













All Days post-transplant
≦180 days post-transplant

















Mean



Mean




Gene/
Mean R
NR
Ratio

Mean R
NR
Ratio


Protein
(n = 39)
(n = 65)
R/NR
p-value*
(n = 28)
(n = 46)
R/NR
p-value*



















27.4
23.9
NA
0.01
28.4
22.4
NA
0.0004


IL1R2
34.3
33.6
0.62
0.009
34.4
33.2
0.44
0.0003


PDCD1
32
32.4
1.32
0.03
32
32.4
1.32
0.06


FLT3
32
31.6
0.76
0.11
32.2
31.5
0.62
0.02


PF4
25
24.8
0.87
0.18
25
24.8
0.87
0.27


ITGAM
26.9
26.8
0.93
0.22
27
26.7
0.81
0.07


SEMA7A
34.3
34.4
1.07
0.31
34.3
34.5
1.15
0.16


RHOU
29.8
29.9
1.07
0.41
29.8
29.9
1.07
0.24


G6b
26.7
26.5
0.87
0.46
26.6
26.5
0.93
0.72


ITGA4
27.6
27.6
1
0.47
27.6
27.7
1.07
0.31


WDR40A
28.9
28.8
0.93
0.68
28.7
28.8
1.07
0.88


MIR
29.4
29.3
0.93
0.82
29.3
29.3
1
0.85







*Significant values in larger red typeface



















TABLE 3













All times post

<180 da post




transplant

transplant



R (n = 38)/

R (n = 27)/



NR (n = 55)

NR (n = 40)












Fold

Fold



Gene/Protein
Change
p-value*
Change
p-value














IL1R1
0.67
0.01
0.55
0.0008


TSC22D3
0.8
0.01
0.72
0.0009


FKBP5
0.85
0.18
0.68
0.007


THBS1
0.73
0.04
0.68
0.03


CD163
0.85
0.2
0.72
0.03


ABCB1
1.1
0.41
1.28
0.07


ANXA1
0.89
0.1
0.86
0.1


IL1B
1.29
0.19
1.45
0.11


EPOR
0.9
0.06
0.91
0.17


DUSP1
0.88
0.39
0.79
0.21


SGK
1.08
0.5
1.16
0.27


TGFB1
0.94
0.19
0.94
0.3


IL7R
1.08
0.54
1.19
0.3


NFKBIA
0.92
0.41
0.9
0.43


NR3C1
1.01
0.76
1.02
0.52


IL4R
0.98
0.75
0.97
0.56


SELP
0.88
0.36
0.93
0.62


IL1RN
0.97
0.73
0.97
0.78


THBS2
0.97
0.74
1.03
0.79


ITGAX
1.02
0.8
0.96
0.86


TNFRSF1
0.94
0.61
1.02
0.89


ADA
1.26
0.002
1.35
0.0008


GZMA
1.19
0.15
1.4
0.01


TRBC1
1.27
0.8
1.5
0.02


FLT3LG
1.16
0.12
1.31
0.03


CD28
1.21
0.12
1.33
0.08


CD8A
1.15
0.37
1.32
0.1


PDCD1L
1.2
0.6
1.21
0.12


CTLA4
1.19
0.17
1.23
0.18


CD274
1.08
0.38
1.15
0.2


CD4
1.01
0.87
1.08
0.35


NFKB1
1.09
0.03
1.1
0.02


TNF
1.21
0.06
1.32
0.03









EXAMPLES
Example 1
Study Objectives and Subjects

Nucleic acid technologies were used to produce gene expression profiles for PBMC samples from subjects who had been treated with various dosages of steroid and were enrolled in the Cardiac Allograft Rejection Gene Expression Observational (CARGO) and the Lung Allograft Rejection Gene Expression Observational (LARGO) studies. All studies were approved by local Institutional Review Boards.


The CARGO study was initiated in 2001 to study gene expression in blood samples as a means for managing transplant rejection in cardiac patients. The eight transplant centers contributing to the studies handle more than 20% of cardiac transplants. The LARGO study was initiated in 2004 to collect blood samples and clinical data, including the results from TBB from lung transplant subjects, at fourteen centers in five different countries.


Microarrays as described in Example 10 were used to study gene expression in 95 samples from LARGO subjects being treated with 5-40 mg of steroid, 68 samples from CARGO subjects being treated with 1-100 mg of steroid, and 56 samples from CARGO or LARGO subjects being treated with 0-50 mg of steroid for CMV infection.


RT-PCR was used in exemplary and pathways studies with PBMC samples from CARGO subjects between 30 days and 12 months post-transplant whose transplants were graded as rejection or non-rejection. The principle inclusion criteria were: a) clinically stable defined as absence of signs or symptoms of acute cardiac allograft rejection, b) histologically stable defined as current EMB indicating non-rejection, c) absence of cardiac dysfunction by invasive hemodynamics and/or echocardiogram, and d) absence of ISHLT (International Society for Heart and Lung Transplant)≧3A rejection, graft dysfunction, or administration of rejection therapy within 30 days prior to enrollment. The demographic and treatment characteristics of the cardiac transplant subjects are shown in the following Table 4.

TABLE 4Subjects-all days postSubjects ≦180 days posttransplanttransplantGroups-No SubjectsR = 39NR = 65p-valueR = 28NR = 46p-valueMedian Age (Range) 60 (25-68) 59 (8-76)0.58  59 (25-68)  59 (8-76)0.73Sex-Male (%) 32 (82.1) 54 (83.1)1  22 (78.6)  41 (53.6)0.31Race-No (%)0.330.025White 23 (59.0) 47 (72.3)  15 (53.6)  38 (82.6)Black 10 (25.6) 10 (15.6)   8 (28.6)  5 (10.9)Other 6 (15.4) 8 (12.1)   5 (17.8)  3 (6.5)Immunosuppression Regimen-No (%)0.320.29Cyclosporine/Mycophenolate 20 (51.3) 37 (56.9)  15 (53.6)  28 (60.9)Cyclosporine/Sirolimus 1 (2.6) 2 (3.1)   1 (3.6)  2 (4.3)Tacrolimus/Mycophenolate 10 (25.6) 19 (29.2)   6 (21.4)  12 (26.1)Tacrolimus/Sirolimus 6 (15.4) 3 (4.6)   5 (17.9)  2 (4.3)Other 2 (5.1) 4 (6.2)   1 (3.6)  2 (4.3)Median Dose (Range)Index Sample 10 (2-30) 10 (1-60)0.6213.25 (2-30)12.5 (1-60)0.75R/NR Sample 7.5 (1-25) 7.5 (2-20)0.8  10 (2-25)  10 (2.5-20)0.6Post-recovery Sample 10 (1-80) 6 (1-20)0.003  10 (2-80) 7.5 (2-20)0.003Days Post-Transplant-Median (Range)Index Sample138 (32-491)133 (33-317)0.3  93 (32-180)  83 (33-177)0.54R/NR Sample180 (53-565)166 (56-342)0.33  130 (53-240) 124 (56-242)0.58Post-recovery Sample189 (62-579)228 (70-471)0.56  155 (62-249) 152 (70-304)0.35Days from Index to R/NR 35 (14-77) 34 (14-76)0.99  32 (14-63)  31 (14-70)0.89ISHLT Biopsy-No (%)0.00060.008Grade 0 12 (30.8) 43 (66.2)   9 (32.1)  30 (65.2)Grade 1A 27 (69.2) 22 (34.4)  19 (67.9)  16 (34.8)


Column 1 of the table characterizes the subjects, immunosuppression regimen, days post-transplant and ISHLT grades. Columns 2, 3 and 4 show the data for rejection (R) and non-rejection (NR) subjects and p-values for characteristics all days post-transplant. Columns 4, 5, and 6 show the data for rejection (R) and non-rejection (NR) subjects and p-values for characteristics <180 days post transplant.


Subjects in both the R and NR groups were on standard steroid weaning protocols with no significant difference (p=0.75) in steroid dose. A two-tailed independent t-test or a Fisher Exact test was used to compare quantitative characteristics, and a Wald (Mann Whitney) test was used to compare categorical characteristics. There was no significant difference in the distribution of characteristics between groups except that ISHLT 1A biopsies and African-Americans were more prevalent in the R group.


Example 2
Sample Collection, Transplant Protocol, and Immunosuppressive Therapy

A blood sample was collected from each subject at each clinical encounter, and clinical data including results of EMB or TBB, immunosuppressive regime, laboratory data, and clinical complications were obtained. Samples were processed as described in Example 8.


Standard cardiac transplant center protocols generally require invasive EMBs to be performed weekly in the first 30 days post transplant (4 biopsies), every two weeks between 31-90 days post transplant (4 biopsies), every 4 weeks between 91-180 days post transplant (3 biopsies), and every 8 weeks between 181-365 days post transplant (3 biopsies). Histology was graded by a local pathologist and two or three pathologists blinded to subject data and outcomes. Agreement of at least two of the pathologists was required to diagnose ISHLT≧3A rejection, and agreement of three pathologists was required for ISHLT 0/1A non-rejection.


Standard lung transplant center protocols generally require at least six invasive TBBs during the first six months post transplant. These tissue samples are examined by at least three pathologists for signs of rejection and rated on a five point ISHLT scale of increasing severity based on the extent of perivascular inflammation, A0=normal lung tissue, A1=minimal, A2=mild, A3=moderate, and A4=severe rejection. A TBB rated≧A2 generally requires therapeutic intervention.


All subjects received center-specific immunosuppressive therapy consisting of cyclosporine or tacrolimus in combination with either mycophenolate mofetil or sirolimus and corticosteroids. The cardiac rejection group (R) had 39 subjects who progressed to acute cellular rejection within 12 weeks. The control group (NR) had 65 subjects who remained rejection-free and were matched with subjects in the rejection group by demographic characteristics, time post-transplant, and immunosuppressive therapy.


Example 3
Steroid Modulated Nucleic Acids and Their Expression

Steroid modulated genes are described in the clusters of Table 1, in the diagnostic set of genes of Table 5, in the pathways genes of Table 3, and among the sequences listed in the published applications and patents incorporated by reference herein in their entirely and shown in the table below.

TABLE 5TitleApplication No; Filing DatePatent/Publication NoMethods And CompositionsUSSN 10/131,827; Apr. 24, 2002USPN 6,905,827For Diagnosing AndPCT/US03/13015; Apr. 24, 2003WO03/090694Monitoring AutoimmuneAnd Chronic InflammatoryDiseasesMethods And CompositionsUSSN 10/325,899;US2003/123086For Diagnosing AndDec. 20, 2002WO04/042346Monitoring TransplantPCT/US03/129456RejectionLeukocyte ExpressionPCT/US01/47856;WO02/057414ProfilingOct. 22, 2001


The steroid modulated genes were identified using at least one statistical method on nucleic acid expression from the microarray study as described in Example 4 and RT-PCR studies as described in Example 5. Primers and probe sets for use in a diagnostic set for detecting genes modulated by steroids can be generated as described in Examples 11 and 12.


Example 4
Microarray Study and Results

Protocols used with the microarrays are described in Examples 9 and 10. For the microarray studies, the manufacturer's software was used to download microarray data. To be included in the analysis, a probe had to be flagged as present (versus marginal or absent) and have a signal of at least 100 for at least 80% of the arrays.


Nucleic acids expressed on Human Genome CGH 44A microarrays (Agilent Technologies, Palo Alto Calif.) that correlated with steroid treatment were identified separately in the samples from the CARGO and LARGO projects. Feature Extraction and GeneSpring software (Agilent Technologies) were used to download microarray data. As shown in the first table in Example 1, the initial filtering flagged 28,997 out of 41,000 probes. Signals were normalized to the median expression of each chip to achieve chip-to-chip comparability.


K-means clustering was applied to the expression of 28,997 nucleic acids in 219 samples as shown in the table below. The parameters for clustering were the number of clusters (20), number of iterations (400), and similarity measure (p-value, Pearson correlation). In one alternative, similarity measure can be a t-test.


In the initial analysis, nucleic acid expression converged after 147 iterations. Using a p-value<0.01, CARGO samples showed expression in 3,604 genes; LARGO samples, in 699 genes. The CARGO and LARGO samples had 278 expressed nucleic acids in common, and cluster 14 (highlighted) was found to be highly enriched in steroid modulated (SM) genes (62.9%), with another 24.7% whose expression correlated with steroid dose (CSD).

TABLE 6ClusterNo. GenesNo. SM Genes% SM Genes% of CSD Genes 1190420.70.1 2156220.70.1 3221820.70.1 4323620.70.1 5221251.80.2 6130510.40.1 7202410.40 8803000 9117420.70.2102059248.61.21197510.40.112121920.70.2133360001470917562.924.715304207.26.616101531.10.317330362.20.21851531.10.619981000201143279.72.4Total2899727810037.2


Column one of Table 6 shows the cluster number; column two, the number of genes in that cluster; column 3, the number of steroid modulated genes; column four, the percent of steroid modulated genes; and column five, the percent of genes correlated with steroid dose.


Candidate steroid modulated nucleic acids (709 genes from cluster 14 and 278 steroid dose correlated genes) were subjected to additional rounds of K-means clustering. The parameters were number of clusters (40), number of iterations (100), and similarity measure (p-value, Pearson correlation). After each round, any cluster containing zero or one steroid modulated nucleic acid was eliminated. Clusters containing two or more steroid modulated nucleic acids were combined for next round of clustering. After four rounds of K-means clustering, 518 genes were in clusters that contained two or more steroid modulated nucleic acids and 262 (50.5%) were nucleic acids whose expression were correlated with steroid dose (data not shown). These 518 genes were subjected to further rounds of clustering with the parameters: number of clusters (10), number of iterations (100), similarity measure (p-value, Pearson correlation). As shown in the table below, all genes had converged into ten clusters after 14 iterations. The 518 steroid modulated genes are described in their respective clusters in Table 1.

TABLE 7ClusterNo. of SM genesNo. CSD Genes111646295553732144540567206151175822816169282810 53Total518262


Column one of Table 7 shows the cluster number; column two, the number of genes; and column three, the number of genes correlated with steroid dose (CSD).


Example 5
RT-PCR Studies and Results

An exemplary RT-PCR study demonstrated the utility of steroid modulated nucleic acids and proteins in diagnosing and monitoring steroid responsiveness. Genes were chosen for the diagnostic set, and nucleic acid expression was reported as threshold cycle (CT) as measured using RT-PCR. The ratios of expression are calculated from the Ct values as 2(Ct(Control)-Ct(Rejection).


Gene expression was processed into a single score using voting, logistic regression or linear algorithms as detailed in Examples 1-3 of U.S. Ser. No. 11/433,191 and in Example 5 of U.S. Pat. No. 6,905,827, both incorporated by reference herein in their entirety. The diagnostic set of the 20 genes (11 formative, six normalization, three control) contained probes that were designed and tested as described in Examples 11 and 12, and RT-PCR, as described in Example 13, was conducted in triplicate RT-PCR reactions on samples from subjects on standard weaning protocols.


Of 104 index subjects, longitudinal gene expression profiles including post rejection and matched post non-rejection samples were available for 34 R subjects and 56 matched NR subjects at similar time points. The findings of the index study were extended to include samples and expression from an additional 192 consecutive subject encounters satisfying the inclusion criteria stated above. This set included samples from 118 new subjects and from 74 previous subjects and was used to estimate the prevalence of non-rejection in any 12 week period following sampling.


Longitudinal changes in expression from the index group were compared to corresponding scores for the larger group of 192 using repeated measure ANOVA. The probability that the transplant would not be rejected (negative predictive value) was calculated using EMB, rejection and non-rejection data. The Wald test was used with multivariate analysis to determine if, after controlling for clinical variables, the gene expression score remained a significant predictor of rejection.


Gene expression score, as calculated using a prediction algorithm, was found to be an independent predictor of future rejection at p=0.0266 when the clinical variables of recipient age, gender and race, panel reactive antibody, CMV serology status, and immunosuppression regimen (Wald test) were included. In fact, independent predictive value at p=0.015, was further enhanced in subjects≦180 days post-transplant.


Table 2 showed the p-value, as calculated using a t-test, for gene expression score and subject nucleic acid expression for 104 index samples, and for the subset of 74 samples <180 days post-transplant. Several of the individual genes shown in Table 2 showed differential expression associated with acute transplant rejection. Expression of IL1R2 decreased significantly (p=0.009, 1.6 fold) and PDCD1, increased significantly (p=0.032, 1.3 fold). In addition, IL1R2 (p<0.001) and FLT3 (p=0.024) demonstrated greater significance during the ≦180 day period and significant decreases in expression (2.3 and 1.6 fold, respectively) in subjects who progressed to rejection. During acute rejection, erythropoiesis genes, MIR and WDR40A, were up-regulated (both p=0.02), and FLT3 was down-regulated (p=0.03). The overall score was also significant using a Wilcoxon test for all subjects who progressed to rejection, p=0.011, and for those who did not progress, p<0.001. Those subjects who showed evidence of incipient rejection were placed immediately on anti-rejection therapy and subsequently showed a significant decrease in gene expression score (p<0.01).


The first RT-PCR study using a diagnostic set corresponding to the genes shown in Table 2 concluded: a) treatment of rejection with high dose steroids led to a statistically significant change in expression, b) low expression scores or a low value derived from evaluating expression scores with a prediction algorithm identified a group of subjects at very low risk for current and future rejection, and c) expression can be used to stratify subjects as to their risk of future rejection and lead to reduced number of cardiac biopsies.


The second RT-PCR study used PMBC samples from CARGO subjects and 33 nucleic acids/genes expressed in steroid modulated pathways. Analyses were based on all samples for which mRNA was available, 93 of 104 subjects in the all times post transplant group and 67 of 74 subjects in the ≦180 days post transplant group. Most of the nucleic acids came from the IL-1 and PDCD1 pathways and nucleic acids induced and expressed in T cells.


Table 3 shows the 33 genes grouped as to pathway, T cell associated, and other (TNF and NFKB1) and presented according to p-Value within the group. Differential expression of the genes is presented as fold change calculated as 2(mean controlCt-mean rejection Ct). Genes whose mRNA levels demonstrated a fold change >1 were up-regulated (increased) in subjects with rejection while those with a fold change <1 were down-regulated (decreased). P-value was based on t-test, and similar significance was obtained using the Mann-Whitney non-parametric test.


Using a p-value <0.05, five of the additional 33 genes tested supported the algorithm's steroid modulated constituents (IL1R2 and FLT3) while six, supported T-cell activation (PDCD1). Specifically, IL1R1, TSC22D3, FKBP5, THBS1 and CD163 showed significantly reduced expression; and ADA, GZMA, TRBC1, NFKB1, TNF and FLT3LG, significantly increased expression. Thus the methods of the invention and diagnostic sets of genes including but not limited to ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3 and selected from Tables 1-3 can be used for determining, diagnosing, evaluating, monitoring, or predicting disease activity, non-rejection, rejection, status of a transplant or of an immune disorder, steroid responsiveness, and treatment plan of a subject with a transplant or immune disorder.


Informative nucleic acids from the RT-PCR studies are listed in the table below as referenced to sequences in U.S. Pat. No. 6,905,827 or GenBank.

GENESEQ ID NOs in USPN 6,905,827CD1633857FKBP56299FLT3See GenBank sequence NM_004119IL1R24685ITGAM1981, 62THB14109, 264


Example 6
Prediction of Rejection or Non-Rejection

Quartile analysis was applied to the exemplary RT-PCR data for 74 subjects≦180 days post transplant. Subjects in the lowest quartile had expression scores≦20, and no subjects progressed to rejection in the subsequent 12 weeks (n=19). Subjects in the top quartile had expression scores≧30, and 58% of these subjects had rejection episodes (n=19) within 12 weeks of histological stability.


When this analysis was extended to the larger group of 192 representative consecutive samples, the incidence of subjects with expression scores≦20 were 33% of samples≦180 days post-transplant, and 98.9% of these remained rejection-free during the ensuing 12 weeks. Since the predictive value did not differ significantly by segmental time periods post transplant (30-60; 61-90; 90-180 days), a clinician can order 2-5 fewer EMBs for a subject with a low risk of rejection during the subsequent 12 weeks.


Example 7
Statistical Methods

The steroid modulated nucleic acids shown in the tables herein were identified in samples from subjects to whom steroids had been administered using at least one statistical method selected from various classification and prediction algorithms, software and programs. These methods include, but are not limited to, analysis of variance, classification and regression trees (Brieman et al. (1984) Classification and Regression Trees, Wadsworth, Belmont Calif.), cluster analysis including K-means clustering (MacQueen (1967) Proceedings of 5th Berkeley Symposium on Mathematical Statistics and Probability, University of California Press 1:281-297), Fisher Exact test, linear discriminatory analysis, logistic regression (Agresti (1996) An Introduction to Categorical Data Analysis. John Wiley and Sons Inc), multiple additive regression trees (Friedman (2002) Stanford University, Stanford Calif.), Mann-Whitney test, multivariate analysis, nearest shrunken centroids classifier (Tibshirani et al. (2002) PNAS 99:6567-6572), significance analysis of microarrays (Tusher et al. (2001) PNAS 98:5116-5121), one and two tailed T-tests, Wald test (Wald (1943) Trans Am Math Soc 54:426-482), Wilcoxon's signed ranks test, quartile analysis, and the like. Many of the above methods can be performed using SAS (SAS Institute, Cary N.C.) or Statistica (Statsoft, Tulsa Okla.). As noted in Example 1, the statistical methods applied to expression in order to chose a diagnostic set of nucleic acids or proteins are fully described in the Examples 1-3 of U.S. Ser. No. 11/433,191 and in Example 16 of U.S. Pat. No. 6,905,827, both incorporated by reference herein in their entirety.


Example 8
Preparation of Blood Samples, RNA Isolation from Lysate

Peripheral blood mononuclear cells (PBMC) were isolated from 8 mL venous blood using a VACUTAINER CPT tube (BD Biosciences (BD), San Jose Calif.) containing the anticoagulant sodium citrate, Ficoll Hypaque density fluid, and a thixotropic polyester gel. After the blood and tube components were mixed by inverting the tube 5-10 times, the tube was centrifuged, and mononuclear cells were collected from the fluid above the barrier layer. Approximately 2 mls of mononuclear cell suspension were transferred to a microfuge tube and centrifuged for 3 min at 16,000 rpm to pellet the cells. The pellet was resuspended and pipetted up and down in 1.8 ml of RLT lysis buffer (Qiagen, Chatsworth Calif.). Cell lysate was frozen and stored at −80 EC until total RNA was isolated.


After adding 5 ml of chloroform to the thawed lysate, the samples were vortexed and incubated at room temperature for 3 min. The aqueous layer was transferred to a new tube and purified using the RNeasy kit (Qiagen) according to the manufacturer's protocol. Isolated RNA was treated with DNAse on a QIASHREDDER column (Qiagen) and purified RNA was eluted in 50 μl of water. RNA purity was checked using the 2100 bioanalyzer and RNA 6000 microfluidics chips (Agilent Technologies, Palo Alto Calif.).


In the alternative, blood samples were collected in PAXgene Blood RNA tubes (Qiagen, Valencia Calif. and total RNA was purified using the PAXgene Blood RNA kit (Qiagen).


Example 9
cDNA Synthesis

cDNA was synthesized from purified RNA using reverse transcription with OLIGO-dT primers/random hexamers (Invitrogen, Carlsbad Calif.) at a final concentration of 0.5 ng/μl and 3 ng/μ, respectively. For the first strand reaction, 0.5 μg of mononuclear RNA and 1 μl of the OLIGO-dT/random hexamers (Invitrogen) were added to water in a reaction tube to a final volume of 11.5 μl. The tube was incubated at 70° C. for 10 min, chilled on ice, centrifuged, and 88.5 μl of first strand buffer mix (Invitrogen) was added to the tube.


The first strand buffer mix contained 1× first strand buffer, 10 mM DTT (Invitrogen), 0.5 mM dATP (New England Biolabs (NEB), Beverly Mass.), 0.5 mM dGTP (NEB), 0.5 mM dTTP (NEB), 0.5 mM dCTP (NEB), 200 U of SUPERSCRIPT RNAse H reverse transcriptase (Invitrogen), and 18 U of RNAGuard inhibitor (GE Healthcare (GEH), Piscataway N.J.). After the reaction was incubated at 42° C. for 90 min, the enzyme was heat-inactivated at 70° C. for 15 min. After adding 2 U of RNAse H (NEB) to the reaction tube, it was incubated at 37° C. for 20 min.


For second strand synthesis, 40 U of E. coli DNA polymerase (Invitrogen) and 2 U RNaseH (Invitrogen) were added to the previous reaction to bring the final volume to 150 μl. Salts and nucleotides were added to a final concentration of 20 mM Tris-HCl (pH 7.0; Fisher Scientific, Pittsburgh Pa.), 90 mM KCl (Teknova, Half Moon Bay Calif.), 4.6 mM MgCl2 (Teknova), 10 mM(NH4)2SO4 (Fisher Scientific), 1× second strand buffer (Invitrogen), 0.266 mM dGTP, 0.266 mM dATP, 0.266 mM dTTP, and 0.266 mM dCTP.


After second strand synthesis for 150 min at 16° C., the cDNA was purified away from the enzymes, dNTPs, and buffers using phenol-chloroform extraction followed by ethanol precipitation in the presence of glycogen. Alternatively, the cDNA was purified on a QIAQUICK silica-gel column (Qiagen) followed by ethanol precipitation in the presence of glycogen. The cDNA was centrifuged at >10,000×g for 30 min. After the supernatant was aspirated, the pellet was washed with 150 μl of 70% ethanol and centrifuged. Following centrifugation, the supernatant was removed, and residual ethanol evaporated.


Example 10
Arrays

Arrays were used to identify steroid modulated genes in gene expression profiles from CARGO and LARGO subjects treated with steroids. In basic format, an array contains reagents specific for at least two nucleic acids or proteins, one that binds to a gene product of the invention and one that binds to a control gene product.


Nucleic Acid Arrays


Human Genome CGH 44A microarrays (Agilent Technologies) were used to determine differential expression. These Cy3/Cy5 chips contained 41,675 probes (60-mers) that represented most the genes found in REFSEQ database (NCBI); additional genes on the chip represented various controls. The chips were run as recommended by the manufacturer and scanned using an Agilent DNA microarray scanner. The data was extracted using Feature Extraction v 7.5 software (Agilent Technologies).


In the alternative, Affymetrix U133A Human GeneChips (Affymetrix, Santa Clara Calif.) with probe sets representing about 14,500 full length genes and 22,000 features were used according to the manuals and product inserts supplied by the manufacturer. Affymetrix Microarray Suite (MAS) v 5.0 software was used to generate expression values for each gene. To correct for slight differences in overall chip hybridization intensity and allow for comparison between samples, each chip was scaled to an overall intensity of 1500.


In another alternative, a low density array containing amplicons produced using probe sets for the nucleic acids selected from Tables 1-3 are harvested from PCR reactions, purified using Sephacryl-400 beads (GEH) and arrayed on a membrane. The membrane is UV irradiated, washed in 0.2% SDS at room temperature and rinsed three times in distilled water. Non-specific binding sites on the array are blocked by incubation in 0.2% casein in PBS for 30 min at 60° C., and the arrays are washed in 0.2% SDS and rinsed in distilled water prior to hybridization.


cDNAs are prepared from subject blood samples; diluted to a concentration of 40-50 ng in 45 μl TE buffer, denatured by heating to 100° C. for five min, and briefly centrifuged. The denatured cDNA is prepared using the Amersham CYSCRIBE first strand cDNA labeling kit (GEH) according to the manufacturer's instructions. The labeling reaction is stopped by adding 5 μl of 0.2M EDTA, and probe is purified from unincorporated nucleotides using a GFX Purification kit (GEH). The purified probe is heated to 100° C. for five min, cooled for two min on ice, and used in membrane-based hybridizations as described below.


Membranes are pre-hybridized in hybridization solution containing 1% Sarkosyl and 1× high phosphate buffer (0.5 M NaCl, 0.1 M Na2HPO4, 5 mM EDTA, pH 7) at 55° C. for two hr. The probe is diluted in 15 ml fresh hybridization solution and added to the membrane. The membrane is hybridized with the probe at 55° C. for 16 hr. Following hybridization, the membrane is washed once for 15 min at 25° C. in 1 mM Tris (pH 8.0) and 1% Sarkosyl and four times for 15 min each at 25° C. in 1 mM Tris (pH 8.0). To detect hybridization complexes, the membrane is exposed to x-ray film (Eastman Kodak) overnight at −70° C., developed, and examined visually or quantified using a scintillation counter (BeckmanCoulter, Fullerton Calif.).


Antibody arrays


Monoclonal antibodies are immobilized on a membrane, slide or dipstick or added to the wells of an ELISA plate using methods well known in the art. The array is incubated in the presence of serum or cell lysate until protein:antibody complexes are formed. The proteins encoded by genes or their splice variants are identified by the known position and labeling of the antibody that binds an epitope of that protein on the array. Quantification is normalized using the antibody:protein complex of various controls.


Example 11
Designing and Selecting Primers

Primers and probe sets were designed for the steroid modulated, normalization, and control genes using the PRIMER3 program (Whitehead Research Institute (WRI), Cambridge Mass.). Default values were used for all parameters but melting temperature (Tm). Tm was set between 71.7 and 73.7° C.; amplicon size, between 50 and 150 bases in length (optimum, about 100 bases); and primers or probes were allowed to be 36 nucleotides in length. Salt concentration, a critical parameter affecting the Tm of the probes and primers, was used at the default concentration, 50 mM.


The C source code for the PRIMER3 program was downloaded and compiled for use on machines running the Windows operating system (Microsoft, Redmond Wash.). To generate a number of potential primers, the program was run in batch mode from the command line using an input file that contained the sequences and the parameters for primer design. The first step was masked out repetitive sequences in the mRNA using the REPEATMASKER program (Institute for Systems Biology, University of Washington, Seattle Wash.). The second step masked out all known SNPs with allelic heterozygosity higher than 1% as annotated in the SNP database at NCBI (Bethesda Md.). The masked sequence was submitted to PRIMER3 using the parameters above, and the top pairs of primers were selected. Alternatively, the Primer3 program was used on the MIT website (Massachusetts Institute of Technology, Cambridge Mass.) to examine a specific region of the mRNA of a gene.


In the alternative, primer design software such as the web-based ProbeFinder software (Roche Diagnostics, Indianapolis Ind.), or provided by other suppliers of oligonucleotides, can be used to design primers and probes sets of the invention. The two step design process requires the name of the target organism and a sequence, gene name, or transcript ID number. The software will identify the Universal ProbeLibrary probes that will detect the nucleic acid. Primers were ordered from Roche Diagnostics, Integrated DNA Technologies (Coralville Iowa), or a similar commercial source.


Example 12
Testing of Primers and Probe Sets

Control genes: Experimental variation was monitored by adding one or more control genes to each array. β-actin, β-GUS, 18s ribosomal subunit, GAPDH, and β2-microglobulin were selected for low variability between samples and high expression across samples.


Primer Testing: Primers were tested at least once to see whether they produced an amplicon of the correct size and to determine their efficiency in a set of RT-PCR reactions using 5 serial dilutions of cDNA in water (1:10, 1:20, 1:40, 1:80, and 1:160). Each primer pair was tested on cDNA made from mononuclear cell RNA. The PCR reaction contained 1× RealTime-PCR buffer (Ambion, Austin Tex.), 2 mM MgCl2 (ABI), 0.2 mM dATP (NEB), 0.2 mM dTTP (NEB), 0.2 mM dCTP (NEB), 0.2 mM dGTP (NEB), 0.625 U AMPLITAQ Gold enzyme (ABI), 0.3 μM of each primer to be used (Sigma Genosys, The Woodlands Tex.), 5 μl of the reverse transcription reaction, and water added to a final volume of 19 μl.


Following 40 cycles of PCR, 10 μl of each PCR product was combined with Sybr Green dye at a final dilution of 1:72,000. Melt curves for each product were determined on a PRISM 7900HT Sequence detection system (ABI), and primer pairs yielding a product with one clean peak were chosen for further analysis. One μl of product from each probe set assay was examined by agarose gel electrophoresis or using a DNA 1000 chip kit and an Agilent 2100 bioanalyzer (Agilent Technologies). From primer design and the genomic sequence, the expected size of the amplicon was known. Only primer pairs showing amplification of the single desired product, and minimal amplification of contaminants, were used in assays.


Example 13
RT-PCR Assays and Analysis

CARGO: Ten μl RT-PCR reactions were performed to evaluate expression in the CARGO samples. TAQMAN Universal PCR Master mix (ABI) was aliquoted into light tight tubes, one for each gene. The primer pair for each gene was added to the tube of PCR master mix labeled for that gene. A FAM/TAMRA dual labeled TAQMAN probe (Biosearch Technologies, Novato Calif.) was added to each tube. Alternatively, different combinations of commercially available fluorescent reporter dyes and quenchers were used such that the absorption wavelength for the quencher matches the emission wavelength for the reporter. In the alternative, Universal ProbeLibrary probes (LNAs; Roche Diagnostics were substituted for TAQMAN probes.


Assays and Analysis: Each sample was dispensed into triplicate wells of a 384 well plate (ABI) for each primer pair. PCR reactions were run on the PRISM 7900HT Sequence Detection system (ABI) with the following conditions: 10 min at 95° C.; 40 cycles of 95° C. for 15 sec, 60° C. for 1 min. Sequence detection system v2.0 software (ABI) was used to analyze the fluorescent signal from each reaction. RT-PCR amplification product was measured as CT during the PCR reaction to observe amplification before any reagent became rate limiting. Threshold was set to a point where all of the reactions were in their linear phase of amplification. A lower CT indicated a higher amount of starting material (greater expression in the sample) since an earlier cycle number meant the threshold was crossed more quickly. A CT of less than 30 based on appropriate cDNA dilutions provided linear results for the blood samples from CARGO subjects. In the alternative, other technologies can be used to measure PCR product. Molecular beacons (Invitrogen) use FRET technology and disposable, microfluidic chip (Thermal Gradient, Pittsford N.Y.) employ silicon wafers to performed 30 cycle PCR in 4.4 min.


Example 14
Labeling Moieties

Labeling moieties can be used for detection of an antibody, nucleic acid or protein in any of the assays or diagnostic kits described herein. These labeling moieties include fluorescent, chemiluminescent, or chromogenic agents, cofactors, enzymes, inhibitors, magnetic particles, radionuclides, reporters/quenchers, substrates and the like that can be attached to or incorporated into the antibody, nucleic acid or protein. Visible labels and dyes include but are not limited to anthocyanins, avidin-biotin, β glucuronidase, biotin, BIODIPY, Coomassie blue, Cy3 and Cy5, 4,6-diamidino-2-phenylindole (DAPI), digoxigenin, ethidium bromide, FAM/TAMRA, FITC, fluorescein, gold, green fluorescent protein, horseradish peroxidase, lissamine, luciferase, phycoerythrin, reporter/quencher pairs (HEX/TAMRA, JOE/TAMRA, ROX/BHQ2, TAMRA/BHQ2, TET/BHQ1, VIC/BHQ1, and the like), rhodamine, spyro red, silver, streptavidin, and the like. Radioactive markers include radioactive forms of hydrogen, iodine, phosphorous, sulfur, and the like. They can be added to a primer or probe or to an antibody using standard protocols well know in the art and described in the specific nucleic acid and protein technologies described in Examples 9-14 and 16-17, respectively.


Example 15
Protein Expression

Adapter sequences for subcloning are added at either end of a coding region specific to a gene or a portion thereof and amplified using PCR. An epitope or affinity tag (6×his) or sequences for secretion from a cell can be added to the adapter sequence to facilitate purification and/or detection of the protein. The amplified cDNA is inserted into a shuttle or expression vector that can replicate in bacteria, insect, yeast, plant, or mammalian cells. Such vectors typically contain a promoter that operably links to the coding region, replication start sites, and antibiotic resistance or metabolite selection sequences.


The expression vector can be used in an in vitro translation system or to transfect cells. For example, Spodoptera frugiperda (Sf9) insect cells are infected with recombinant Autographica californica nuclear polyhedrosis virus (baculovirus). The polyhedrin gene is replaced with the cDNA by homologous recombination, and the polyhedrin promoter drives transcription. The protein is synthesized as a fusion protein with an affinity tag that enables purification.


Clones of transformed cells are analyzed to ensure that the inserted sequence is expressed. Once expression is verified, the cells are grown under selective conditions; and the protein is isolated from cells, or if secreted, from the growth media using chromatography, size exclusion chromatography, immunoaffinity chromatography, or other methods including cell fractionation, ion exchange, or selective precipitation.


The isolated and purified protein is then used as a reagent on an array or as an antigen to produce specific antibodies.


Example 16
Antibody Production and Testing

If antibodies are to be used as reagents, the sequence of the gene or splice variant is analyzed to determine regions of high immunogenicity (LASERGENE software; DNASTAR, Madison Wis.), and an appropriate oligopeptide is synthesized and conjugated to keyhole lympet hemocyanin (KLH; Sigma-Aldrich, St Louis Mo.).


Immunization


Rabbits are injected with the oligopeptide-KLH complexes in complete Freund=s adjuvant, and the resulting antisera is tested for specific recognition of the protein or fragments thereof. Antisera that react positively with the protein are affinity purified on a column containing beaded agarose resin to which the synthetic oligopeptide has been conjugated (SULFOLINK kit; Pierce Chemical, Rockford Ill.). The column is equilibrated using 12 ml IMMUNOPURE Gentle Binding buffer (Pierce Chemical). Three ml of rabbit antisera is combined with one ml of binding buffer and poured into the column. The column is capped (on the top and bottom), and antisera is allowed to bind with the oligopeptide by gentle shaking at room temperature for 30 min. The column is allowed to settle for 30 min, drained by gravity flow, and washed with 16 ml binding buffer (4×4 ml additions of buffer). The antibody is eluted in one ml fractions with IMMUNOPURE Gentle Elution buffer (Pierce Chemical), and absorbance at 280 nm is determined. Peak fractions are pooled and dialyzed against 50 mM Tris, pH 7.4, 100 mM NaCl, and 10% glycerol. After dialysis, the concentration of the purified antibody is determined using the BCA assay (Pierce Chemical), aliquoted, and frozen.


Electrophoresis and Blotting


Samples containing protein are mixed in 2× loading buffer, heated to 95° C. for 3-5 min and loaded on 4-12% NUPAGE Bis-Tris precast gel (Invitrogen). Unless indicated, equal amounts of total protein are loaded into each well. The gel is electrophoresed in 1× MES or MOPS running buffer (Invitrogen) at 200 V for approximately 45 min on an XCELL II apparatus (Invitrogen) until the RAINBOW marker (GEH) resolves and the dye front approaches the bottom of the gel. The gel is soaked in 1×transfer buffer (Invitrogen) with 10% methanol for a few minutes; and a PVDF membrane (Millipore, Billerica Mass.) is soaked in 100% methanol for a few seconds to activate it. The membrane, the gel, and supports are placed on the TRANSBLOT SD transfer apparatus (Biorad, Hercules Calif.) and a constant current of 350 mA is applied for 90 min.


Conjugation with Antibody and Visualization


After the proteins are transferred to the membrane, it is blocked in 5% (w/v) non-fat dry milk in 1× phosphate buffered saline (PBS) with 0.1% Tween 20 detergent (blocking buffer) on a rotary shaker for at least 1 hr at room temperature or at 4° C. overnight. After blocking, the buffer is removed, and 10 ml of primary antibody in blocking buffer is added and incubated on the rotary shaker for 1 hr at room temperature or overnight at 4° C. The membrane is washed 3 times for 10 min each with PBS-Tween (PBST), and secondary antibody, conjugated to horseradish peroxidase, is added at a 1:3000 dilution in 10 ml blocking buffer. The membrane and solution are shaken for 30 min at room temperature and washed three times for 10 min with PBST.


The wash solution is carefully removed, and the membrane is moistened with ECL+chemiluminescent detection system (GEH) and incubated for approximately 5 min. The membrane, protein side down, is placed on x-ray film (Eastman Kodak, Rochester N.Y.) and developed for approximately 30 seconds. Antibody:protein complexes are visualized and/or scanned and quantified.

Claims
  • 1. A method of diagnosing or monitoring steroid responsiveness of a subject comprising: a) detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression is correlated with steroid administration or dosage; and b) applying at least one statistical method to the expression of the diagnostic set to diagnose or monitor steroid responsiveness of the subject.
  • 2. The method of claim 1 wherein the diagnostic set further comprises at least one steroid modulated nucleic acid selected from each of at least two of the clusters of Table 1.
  • 3. The method of claim 1 wherein the diagnostic set further comprises two or more steroid modulated nucleic acids selected from Tables 2 and Table 3.
  • 4. The method of claim 1 wherein detecting the expression further comprises using hybridization or quantitative real-time polymerase chain reaction (RT-PCR).
  • 5. The method of claim 1 wherein the sample further comprises a fluid obtained from the subject by any sampling means.
  • 6. The method of claim 1 wherein the sample is blood containing peripheral blood mononuclear cells (PMBC).
  • 7. The method of claim 1 wherein detecting expression of the diagnostic set of steroid modulated nucleic acids further comprises isolating RNA from the sample.
  • 8. The method of claim 1 wherein the statistical method is K-means clustering or a prediction algorithm.
  • 9. The method of claim 8 wherein K-means clustering produces clusters of genes that are correlated by p-value and their expression in a cell type or pathway.
  • 10. The method of claim 8 wherein the prediction algorithm is selected from a linear algorithm, a logistic regression algorithm, and a voting algorithm and produces a single value or score.
  • 11. The method of claim 1 wherein detecting expression of a diagnostic set further comprises selecting at least two oligonucleotides or a probe set to detect the expression of each nucleic acid of the diagnostic set.
  • 12. A kit comprising the oligonucleotides or probe sets of claim 13.
  • 13. The method of claim 1 wherein diagnosing or monitoring steroid responsiveness further comprises detecting the expression of nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.
  • 14. A method for predicting rejection or non-rejection in a subject with a transplant comprising: a) detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression of the steroid modulated nucleic acids correlates with transplant rejection or non-rejection; and b) applying at least one statistical method to the expression of the diagnostic set of steroid modulated nucleic acids to predict rejection or non-rejection.
  • 15. The method of claim 14 wherein the diagnostic set of steroid modulated nucleic acids further comprises two or more nucleic acids selected from Tables 1-3.
  • 16. The method of claim 14 wherein the sample is PMBC.
  • 17. The method of claim 14 wherein detecting expression of the diagnostic set of steroid modulated nucleic acids further comprises isolating RNA from the sample.
  • 18. The method of claim 14 wherein detecting expression of the diagnostic set of steroid modulated nucleic acids further comprises using RT-PCR.
  • 19. The method of claim 14 wherein the statistical method is a prediction algorithm that produces single value or score that correlates with rejection or non-rejection.
  • 20. The method of claim 19 wherein the score that correlates with non-rejection is<20 and the score that correlates with rejection is>30.
  • 21. The method of claim 14 wherein predicting rejection or non-rejection further comprises detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.
  • 22. A method of diagnosing or monitoring the status of a subject with a transplant comprising: a) detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression is correlated with dysfunction or rejection of the transplant; and b) applying at least one statistical method to the expression of the nucleic acids to monitor the status of the transplant.
  • 23. The method of claim 22 wherein the diagnostic set further comprises two or more nucleic acids selected from Tables 1-3.
  • 24. The method of claim 22 wherein the sample is PMBC.
  • 25. The method of claim 22 wherein detecting expression further comprises isolating RNA from the sample.
  • 26. The method of claim 22 wherein detecting expression further comprises using RT-PCR.
  • 27. The method of claim 22 wherein the statistical method is a prediction algorithm that produces single value or score that correlates with the status of a subject with a transplant.
  • 28. The method of claim 22 wherein diagnosing and monitoring the status of a subject with a transplant further comprises detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3.
  • 29. A method for designing and monitoring a treatment plan for a subject with a transplant or an immune disorder comprising: a) detecting expression of a diagnostic set of at least two steroid modulated nucleic acids in a sample from the subject wherein the expression correlates with the steroid responsiveness of the subject; and b) using the expression of the diagnostic set of steroid modulated nucleic acids to design and monitor the treatment plan of the subject.
  • 30. The method of claim 29 wherein the diagnostic set of steroid modulated nucleic acids further comprises two or more nucleic acids selected from Tables 1-3.
  • 31. The method of claim 29 wherein the sample is PMBC.
  • 32. The method of claim 29 wherein detecting expression further comprises isolating RNA from the sample.
  • 33. The method of claim 29 wherein detecting expression further comprises using RT-PCR.
  • 34. The method of claim 29 wherein the statistical method is a prediction algorithm.
  • 35. The method of claim 29 wherein diagnosing and monitoring the treatment plan of a subject with a transplant or immune disorder further comprises detecting the expression of a diagnostic set of steroid modulated nucleic acids encoding ADA, CD163, FKBP5, FLT3, FLT3LG, GZMA, IL1R1, IL1R2, ITGAM, NFKB1, PDCD1, THBS1, TNF, TRBC1 and TSC22D3 whose expression correlates with steroid responsiveness of a subject.
  • 36. The method of claim 29 wherein the transplant is selected from bone marrow, heart, kidney, liver, lung, pancreas, pancreatic islets, stem cells, xenotransplants, and artificial implants.
  • 37. The method of claim 29 wherein the immune disorder is selected from cytomegalovirus infection, multiple sclerosis, and systemic lupus erythematosus.
  • 38. A method for using primers and probe sets to detect steroid responsiveness of a subject with a transplant or an immune disorder comprising: a) designing and generating primers or probe sets for nucleic acids whose expression is modulated by steroid administration or dosage; and b) using RT-PCR and the primers or probe sets on a sample from the subject to detect steroid responsiveness.
  • 39. The method of claim 38 wherein the nucleic acids whose expression is modulated by steroid administration or dosage are selected from Tables 1-3.
RELATED APPLICATION

This application claims priority to U.S. Patent Application No. 60/790,474, filed 7 Apr. 2006, which is incorporated by reference herein in its entirety.

Provisional Applications (1)
Number Date Country
60790474 Apr 2006 US