The present invention relates to the identification and use of diagnostic markers related to venous thromboembolic disease (“VTED”). In a various aspects, the invention relates to methods and compositions for use in the diagnosis of VETD, pulmonary embolism, and deep venous thrombosis, and in the stratification of risk in such patients.
The following discussion of the background of the invention is merely provided to aid the reader in understanding the invention and is not admitted to describe or constitute prior art to the present invention.
Venous thromboembolic disease (“VTED”) represents a spectrum of conditions that includes deep venous thrombosis (DVT) and pulmonary embolism (PE). The estimated annual incidence of VETD is 117 cases per 100,000 persons. The incidence rises markedly in persons 60 years and older and may be as high as 900 cases per 100,000 by the age of 85 years. Silverstein et al., Arch. Intern. Med. 158: 585-93, 1998. Risk factors for VETD include increasing age, prolonged immobility, surgery, trauma, malignancy, pregnancy, estrogenic medications (e.g., oral contraceptive pills, hormone therapy, tamoxifen), congestive heart failure, hyperhomocystinemia, diseases that alter blood viscosity (e.g., polycythemia, sickle cell disease, multiple myeloma), and inherited thrombophilias. About 75 percent of patients with VETD have at least one established risk factor. Heit et al., Arch. Intern. Med. 162: 1245-48, 2002.
Most clinically important PEs originate from proximal DVT of the leg, particularly the popliteal, femoral, or iliac veins. Moser and LeMoine, Ann. Intern. Med. 94(4 pt 1): 439-44, 1981. Upper extremity DVT is less common but also may lead to PE. A much less common cause of upper extremity DVT is Paget-Schroetter syndrome (idiopathic upper extremity DVT in young athletes). Classic symptoms of DVT include swelling, pain, and discoloration in the affected extremity. Physical examination may reveal the palpable cord of a thrombosed vein, unilateral edema, warmth, and superficial venous dilation. Such classic signs of DVT are of low predictive value and can occur in other conditions such as musculoskeletal injury, cellulitis, and venous insufficiency. Like DVT, PE is also characterized by a constellation of nonspecific signs and symptoms that are associated with other diseases. The most common symptoms in individuals without preexisting cardiopulmonary disease are dyspnea, pleuritic chest pain, cough, leg edema, leg pain, hemoptysis, and palpitations.
The fibrin degradation polypeptide known as “D-dimer” is a marker of endogenous fibrinolysis and should therefore be detectable in the blood of patients with venous thromboembolic disease. One recent study found that the combination of clinical assessment and a negative D-dimer assay result effectively rules out DVT. Wells et al., N. Engl. J. Med. 349: 1227-1235, 2003. Similarly, use of the D-dimer tests in combination with the clinical assessment is effective in ruling out PE in patients who present to the emergency department. Wells et al., Ann. Intern. Med. 135: 98-107, 2001.
The present invention relates to the identification and use of markers for the detection of venous thromboembolic disease and in the stratification of risk in VTED patients. The methods and compositions described herein can meet the need in the art for rapid, sensitive and specific diagnostic assay to be used in the diagnosis and differentiation of various forms of VTED. Moreover, the methods and compositions of the present invention can also be used to facilitate the treatment of patients and the development of additional diagnostic and/or prognostic indicators.
In various aspects, the invention relates to materials and procedures for identifying markers that are associated with the diagnosis, prognosis, or differentiation of VTED in a patient; to using such markers in diagnosing and treating a patient and/or to monitor the course of a treatment regimen; to using such markers to identify subjects at risk for one or more adverse outcomes related to VTED; and for screening compounds and pharmaceutical compositions that might provide a benefit in treating or preventing such conditions, e.g., for efficacy.
In a first aspect of the present invention, methods for diagnosing VTED, PE, and/or DVT are described. Such methods comprise performing one or more assays on one or more test samples obtained from the subject, such assay(s) being configured to detect one or more markers as described herein, and using results of the assays performed, typically in the form of a presence or amount of the markers, to assign the presence or absence of VTED, PE, and/or DVT to the subject. A plurality of different assays detecting different biochemical markers are preferably used togetherin a diagnostic “panel.” Such panels may comprise 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, or more or individual assays, and optionally comprise one or more values related to physical characteristics of the subject. These panels most preferably comprise performing a plurality of assays on one or more test samples obtained from the subject, such assays being configured to detect a plurality of markers selected from the group consisting of markers related to blood pressure regulation, markers related to inflammation, markers related to apoptosis, markers related to reactive oxygen species, markers related to myocardial injury, markers related to pulmonary injury, and markers related to coagulation and hemostasis.
In a related aspect, the invention features methods of assigning a risk of one or more clinical outcomes to a subject suffering from VTED, PE, and/or DVT. Such methods comprise performing one or more assays on one or more test samples obtained from the subject, such assay(s) being configured to detect one or more markers as described herein, and using results of the assays performed, typically in the form of a presence or amount of the markers, to assign a risk of one or more clinical outcomes to the subject. As in the case of diagnosis, a plurality of different assays detecting different biochemical markers are preferably used togetherin a prognostic “panel.” Such panels may comprise 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, or more or individual assays, and optionally comprise one or more values related to physical characteristics of the subject.
Suitable markers for use in such methods are described in detail hereinafter. Preferably, such methods comprise performing one or more assays configured to detect one or more markers selected from the group consisting of acidic calponin, adrenomedullin, angiopoietin-4, basic calponin, BMP-4, BNP, BNP1-108, BNP3-108, BNP79-108, caspase-3, CCL11, CGRP, CK-BB, CK-MB, CRP, D-dimer, sELAF, endothelin-1, GSTP, hFABP, IL-1ra, IL-25, leptin, sLTBR, MCP-1, MMP-9, MPO, NDKA, neuropilin-2, NGAL, PLGF-1, PLGF 1+2, activated protein C, total protein C, PSAP-A, PSAP-B, PSAP-C, PSAP-D, sRAGE, sPECAM-1, spectrin 145, TIE-2, tissue factor, sTNFR1a, sTNFRSF7, TNFsR14, TpP, UFDP1H, UPA, VCAM-1, VE-cadherin, VEGF, sVEGF-R1, VWF-integrin, and ANP28-151. These preferred markers may be used individually, or in panels that comprise 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, or more or individual assays, and optionally comprise one or more values related to physical characteristics of the subject.
In certain embodiments, result(s) of assay(s) can each be compared to a level (a “threshold”) that is associated with the diagnosis or prognosis of VTED, PE, and/or DVT. By correlating each of the subject's selected marker results to thresholds for each marker of interest, the subject may be assigned to a diagnostic group (e.g., suffering from one of these conditions, or not suffering from one of these conditions. Similarly, by correlating the subject's marker results to prognostic thresholds for each marker, the probability or risk that the subject will suffer one or more future adverse outcomes may be determined.
In other embodiments, particular thresholds for one or more markers in a panel are not relied upon to determine if a profile of marker levels obtained from a subject are indicative of a particular diagnosis or prognosis. Rather, the present invention may utilize an evaluation of the entire profile of markers. For example, by plotting ROC curves for the sensitivity of a particular panel of markers versus 1-(specificity) for the panel at various marker levels, a profile of marker measurements from a subject may be considered together to provide a global probability (a “panel response” expressed either as a numeric score or as a percentage risk) that the symptom(s) observed in an individual are caused by a particular underlying disease. In such embodiments, an increase (or decrease) in a certain subset of markers may be sufficient to indicate a particular diagnosis in one patient, while an increase (or decrease) in a different subset of markers may be sufficient to indicate the same or a different diagnosis in another patient. Methods for performing such analyses are described hereinafter.
In yet other embodiments, multiple determinations of one or more markers can be made, and a temporal change in the markers can be used to rule in or out one or more particular etiologies for observed symptom(s). For example, one or more markers may be determined at an initial time, and again at a second time, and the change (or lack thereof) in the marker level(s) over time determined. In such embodiments, an increase in the marker from the initial time to the second time may be diagnostic of a particular disease underlying one or more symptoms, a particular prognosis, etc. Likewise, a decrease in the marker from the initial time to the second time may be indicative of a particular disease underlying one or more symptoms, a particular prognosis, etc. Temporal changes in one or more markers may also be used together with single time point marker levels to increase the discriminating power of marker panels. In yet another alternative, a “panel response” may be treated as a marker, and temporal changes in the panel response may be indicative of a particular disease underlying one or more symptoms, a particular prognosis, etc.
Receiver Operating Characteristic curves, or “ROC” curves, are typically calculated by plotting the value of a variable versus its relative frequency in “normal” and “disease” populations, where “normal” and “disease” simply indicates the absence and presence of some characteristic of interest. For any particular marker, a distribution of marker levels for subjects with and without a “disease” will likely overlap. Under such conditions, a test does not absolutely distinguish normal from disease with 100% accuracy, and the area of overlap indicates where the test cannot distinguish normal from disease. A threshold is selected, above which (or below which, depending on how a marker changes with the disease) the test is considered to be abnormal and below which the test is considered to be normal. The area under the ROC curve is a measure of the probability that the perceived measurement will allow correct identification of a condition. ROC curves can be used even when test results don't necessarily give an accurate number. As long as one can rank results, one can create an ROC curve. For example, results of a test on “disease” samples might be ranked according to degree (say 1=low, 2=normal, and 3=high). This ranking can be correlated to results in the “normal” population, and a ROC curve created. These methods are well known in the art. See, e.g., Hanley et al., Radiology 143: 29-36 (1982).
Preferably, a plurality of marker assay results are combined to increase the predictive value of the analysis in comparison to that obtained from the markers individually. The skilled artisan will also understand that diagnostic markers, differential diagnostic markers, prognostic markers, time of onset markers, etc., may be combined in a single assay or device. For example, certain marker assay results measured by a device or instrument may be used to diagnose VTED, PE, and/or DVT. Each condition may be diagnosed with sets of markers that may comprise unique markers, or may include markers that overlap with one or both of the other sets. Markers may also be commonly used for multiple purposes by, for example, applying a different set of analysis parameters (e.g., a threshold or a different weighting factor) to the marker(s) for the different purpose(s). For example, a marker at one concentration or weighting may be used, alone or as part of a larger panel, to indicate a diagnosis of VTED, PE, and/or DVT, and the same marker at a different concentration or weighting may be used, alone or as part of a larger panel, to indicate prognosis associated with the diagnosis.
The sensitivity and specificity of a diagnostic and/or prognostic test depends on more than just the analytical “quality” of the test—they also depend on the definition of what constitutes an abnormal result. In practice, In preferred embodiments, markers and/or marker panels are selected to exhibit at least 75% sensitivity, more preferably at least 80% sensitivity, even more preferably at least 85% sensitivity, still more preferably at least 90% sensitivity, and most preferably at least 95% sensitivity, combined with at least 75% specificity, more preferably at least 80% specificity, even more preferably at least 85% specificity, still more preferably at least 90% specificity, and most preferably at least 95% specificity. In particularly preferred embodiments, both the sensitivity and specificity are at least 75%, more preferably at least 80%, even more preferably at least 85%, still more preferably at least 90%, and most preferably at least 95%.
In other embodiments, a positive likelihood ratio, negative likelihood ratio, odds ratio, or hazard ratio is used as a measure of a test's ability to predict risk or diagnose a disease. In the case of a positive likelihood ratio, a value of 1 indicates that a positive result is equally likely among subjects in both the “diseased” and “control” groups; a value greater than 1 indicates that a positive result is more likely in the diseased group; and a value less than 1 indicates that a positive result is more likely in the control group. In the case of a negative likelihood ratio, a value of 1 indicates that a negative result is equally likely among subjects in both the “diseased” and “control” groups; a value greater than 1 indicates that a negative result is more likely in the test group; and a value less than 1 indicates that a negative result is more likely in the control group. In certain preferred embodiments, markers and/or marker panels are preferably selected to exhibit a positive or negative likelihood ratio of at least about 1.5 or more or about 0.67 or less, more preferably at least about 2 or more or about 0.5 or less, still more preferably at least about 5 or more or about 0.2 or less, even more preferably at least about 10 or more or about 0.1 or less, and most preferably at least about 20 or more or about 0.05 or less. The term “about” in this context refers to +/−5% of a given measurement.
In the case of an odds ratio, a value of 1 indicates that a positive result is equally likely among subjects in both the “diseased” and “control” groups; a value greater than 1 indicates that a positive result is more likely in the diseased group; and a value less than 1 indicates that a positive result is more likely in the control group. In certain preferred embodiments, markers and/or marker panels are preferably selected to exhibit an odds ratio of at least about 2 or more or about 0.5 or less, more preferably at least about 3 or more or about 0.33 or less, still more preferably at least about 4 or more or about 0.25 or less, even more preferably at least about 5 or more or about 0.2 or less, and most preferably at least about 10 or more or about 0.1 or less. The term “about” in this context refers to +/−5% of a given measurement.
In the case of a hazard ratio, a value of 1 indicates that the relative risk of an endpoint (e.g., death) is equal in both the “diseased” and “control” groups; a value greater than 1 indicates that the risk is greater in the diseased group; and a value less than 1 indicates that the risk is greater in the control group. In certain preferred embodiments, markers and/or marker panels are preferably selected to exhibit a hazard ratio of at least about 1.1 or more or about 0.91 or less, more preferably at least about 1.25 or more or about 0.8 or less, still more preferably at least about 1.5 or more or about 0.67 or less, even more preferably at least about 2 or more or about 0.5 or less, and most preferably at least about 2.5 or more or about 0.4 or less. The term “about” in this context refers to +/−5% of a given measurement.
One or more markers may lack predictive value when considered alone, but when used as part of a panel, such markers may be of great value in determining a particular diagnosis/prognosis. Weighting factors may also be applied to one or more markers in a panel, for example, when a marker is of particularly high utility in identifying a particular diagnosis/prognosis, it may be weighted so that at a given level it alone is sufficient to signal a positive result. Likewise, a weighting factor may provide that no given level of a particular marker alone is sufficient to signal a positive result, but only signals a result in combination with one or more other markers in the panel.
While exemplary panels are described herein, assays to detect one or more particular markers may be replaced, added, or subtracted from these exemplary panels while still providing clinically useful results. Panels may comprise both specific markers of VTED, PE, and/or DVT; and/or non-specific markers (e.g., markers that are increased or decreased due to inflammation, regardless of the cause; markers that are increased or decreased due to changes in hemostasis, regardless of the cause, etc.). While non-specific (and/or specific) markers may not individually be diagnostic of VTED, PE, and/or DVT, a particular “fingerprint” pattern of changes may, in effect, act as a specific indicator of a disease or condition. As discussed above, that pattern of changes may be obtained from a single sample, or may optionally consider temporal changes in one or more members of the panel (or temporal changes in a panel response value).
In a particularly preferred embodiment, the plurality of assays used in the diagnostic and prognostic methods and compositions described herein comprises at least one, and preferably two or more, assays configured to detect marker(s) related to coagulation and hemostasis. Particularly preferred markers related to coagulation and hemostasis include those selected from the group consisting of plasmin, thrombin, antithrombin-III, fibrinogen, one or more forms of von Willebrand factor, D-dimer, PAI-1, soluble urokinase plasminogen activator surface receptor (uPAR), Protein C, soluble endothelial protein C receptor (EPCR), TAFI, fibrinopeptide A, plasmin alpha 2 antiplasmin complex, platelet factor 4, platelet-derived growth factor, P-selectin, prothrombin fragment 1+2, B-thromboglobulin, thrombin antithrombin III complex, thrombomodulin, thrombus precursor protein, tissue factor, tissue factor pathway inhibitor-α, and tissue factor pathway inhibitor-β, or markers related thereto. Most preferred markers are selected from the group consisting of D-dimer, uPAR, TpP, and one or more forms of von Willebrand factor, or markers related thereto.
In another particularly preferred embodiment, the plurality of assays used in the diagnostic and prognostic methods and compositions described herein comprises at least one, and preferably two or more, assays configured to detect marker(s) related to blood pressure regulation. Particularly preferred markers related to blood pressure regulation are selected from the group consisting of atrial natriuretic peptide (“ANP), pro-ANP, B-type natriuretic peptide (“BNP”), NT-pro BNP, pro-BNP C-type natriuretic peptide, urotensin II, urocortin I, urocortin II, urocortin III, arginine vasopressin, aldosterone, angiotensin I, angiotensin II, angiotensin III, bradykinin, calcitonin, procalcitonin, calcitonin gene related peptide, adrenomedullin, calcyphosine, endothelin-2, endothelin-3, renin, and urodilatin, or markers related thereto. Most preferred markers are selected from the group consisting of BNP, proBNP, NT-proBNP, BNP79-108, and BNP3-108, or markers related thereto.
In still another particularly preferred embodiment, the plurality of assays used in the diagnostic and prognostic methods and compositions described herein comprises at least one, and preferably two or more, assays configured to detect marker(s) related to inflammation. Particularly preferred markers are selected from the group consisting of acute phase reactants, cell adhesion molecules such as vascular cell adhesion molecule (“VCAM”), intercellular adhesion molecule-1 (“ICAM-1”), intercellular adhesion molecule-2 (“ICAM-2”), and intercellular adhesion molecule-3 (“ICAM-3”), C-reactive protein, HMG-1 (also known as HMGB1), interleukins such as IL-1β, IL-6, IL-8, interleukin-1 receptor agonist, monocyte chemotactic protein-1, caspase-3, lipocalin-type prostaglandin D synthase, mast cell tryptase, eosinophil cationic protein, KL-6, haptoglobin, tumor necrosis factor α, tumor necrosis factor β, Fas ligand, soluble Fas (Apo-1), TRAIL, TWEAK, fibronectin, macrophage migration inhibitory factor (MIF), and vascular endothelial growth factor (“VEGF”), or markers related thereto.
Acute phase reactants may be selected from the group consisting of hepcidin, HSP-60, HSP-65, HSP-70, asymmetric dimethylarginine (an endogenous inhibitor of nitric oxide synthase), matrix metalloproteins 11, 3, and 9, defensin HBD 1, defensin HBD 2, serum amyloid A, oxidized LDL, insulin like growth factor, transforming growth factor β, inter-α-inhibitors, e-selectin, glutathione-S-transferase, hypoxia-inducible factor-1α, inducible nitric oxide synthase (“I-NOS”), intracellular adhesion molecule, lactate dehydrogenase, matrix metalloproteinase-9 (“MMP-9”), monocyte chemoattractant peptide-1 (“MCP-1”), n-acetyl aspartate, prostaglandin E2, receptor activator of nuclear factor (“RANK”) ligand, TNF receptor superfamily member 1A, and cystatin C, or markers related thereto.
Likewise, the plurality of assays used in the diagnostic and prognostic methods and compositions described herein may comprise at least one, and preferably two or more, assays configured to detect marker(s) related to reactive oxygen species. The marker(s) may be selected from the group consisting of superoxide dismutase, glutathione, α-tocopherol, ascorbate, inducible nitric oxide synthase, lipid peroxidation products, nitric oxide, myeloperoxidase, and breath hydrocarbons (preferably ethane), or markers related thereto.
Additional markers and/or marker classes may be added to such panels to provide further ability to discriminate amongst diseases. For example, the inflammatory response and resulting effects on capillaries and reduced oxygenation of tissues implicate one or more markers related to the acute phase response, one or more markers related to vascular tissues, and one or more tissue-specific (e.g., neural-specific) markers, the levels of which are increased in ischemic conditions. Preferably, one or more markers selected from the group consisting of α-2 actin, basic calponin 1, β-1 integrin, acidic calponin, caldesmon, cysteine rich protein-2 (“CRP 2” or “CSRP 2”), elastin, fibrillin 1, latent transforming growth factor beta binding protein 4 (“LTBP 4”), smooth muscle myosin, smooth muscle myosin heavy chain, and transgelin, or markers related thereto (referred to collectively as “markers related to vascular tissue”) may be included in such a panel. Additional marker classes, such as markers related to myocardial injury, markers related to neural tissue injury, markers related to pulmonary injury, etc., are described hereinafter.
Preferred markers related to pulmonary injury may be selected from the group consisting of neutrophil elastase, KL-6, LAMP 3, LAMP3, pulmonary surfactant protein A, pulmonary surfactant protein B, pulmonary surfactant protein C, pulmonary surfactant protein D, phospholipase D, PLA2G5, SFTPC, HT156, and HTII280, or markers related thereto.
Preferred markers related to myocardial injury may be selected from the group consisting of cardiac troponin I (free and/or complexed), cardiac troponin T (free and/or complexed), annexin V, B-enolase, CK-MB, glycogen phosphorylase-BB, heart type fatty acid binding protein, phosphoglyceric acid mutase, and S-100ao, or markers related thereto.
The markers and marker assays described herein may be variously combined to provide suitable marker panels. In various embodiments, the plurality of markers comprises at least one marker related to coagulation and hemostasis and at least one marker related to inflammation; the plurality of markers comprises at least one marker related to blood pressure regulation and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to inflammation, and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to apoptosis and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to reactive oxygen species and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to myocardial injury and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to pulmonary injury and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to myocardial injury, and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to pulmonary injury, and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to apoptosis, and at least one marker related to coagulation and hemostasis; the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to apoptosis, at least one marker related to inflammation, and at least one marker related to coagulation and hemostasis; or the plurality of markers comprises at least one marker related to blood pressure regulation, at least one marker related to myocardial injury, and at least one marker related to coagulation and hemostasis. These combinations are not meant to be limiting.
These markers may be combined in various combinations. For example, preferred methods may include 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more assays configured to detect one or more markers selected from the group consisting of B-type natriuretic peptide, NT-proBNP, proBNP, BNP79-108, BNP3-108, pulmonary surfactant protein A, B, C, and/or D, caspase-3, CRP, D-dimer, TpP, MCP-1, MMP-9, myeloperoxidase, free cardiac troponin I, complexed cardiac troponin I, free and complexed cardiac troponin I, free cardiac troponin T, complexed cardiac troponin T, free and complexed cardiac troponin T, total cardiac troponin, or marker(s) related thereto.
Particularly preferred panels comprise performing assays configured to detect D-dimer, one or more markers related to blood pressure regulation, preferably selected from the group consisting of B-type natriuretic peptide, NT-proBNP, proBNP, BNP79-108, or BNP3-108, and/or one or more markers related to myocardial injury, preferably selected from the group consisting of free cardiac troponin I, complexed cardiac troponin I, free and complexed cardiac troponin I, free cardiac troponin T, complexed cardiac troponin T, free and complexed cardiac troponin T, total cardiac troponin, annexin V, B-enolase, CK-MB, glycogen phosphorylase-BB, heart type fatty acid binding protein, phosphoglyceric acid mutase, and S-100ao, or marker(s) related thereto.
Other particularly preferred panels comprise performing assays configured to detect D-dimer, one or more markers related to blood pressure regulation, preferably selected from the group consisting of B-type natriuretic peptide, NT-proBNP, proBNP, BNP79-108, or BNP3-108, and/or one or more markers related to pulmonary injury, preferably selected from the group consisting of neutrophil elastase, KL-6, LAMP 3, LAMP3, pulmonary surfactant protein A, pulmonary surfactant protein B, pulmonary surfactant protein C, pulmonary surfactant protein D, phospholipase D, PLA2G5, SFTPC, HTI56, and HTII280, or markers related thereto.
Still other particularly preferred panels comprise performing assays configured to detect D-dimer, one or more markers related to blood pressure regulation, preferably selected from the group consisting of B-type natriuretic peptide, NT-proBNP, proBNP, BNP79-108, or BNP3-108, one or more markers related to myocardial injury, preferably selected from the group consisting of free cardiac troponin I, complexed cardiac troponin I, free and complexed cardiac troponin I, free cardiac troponin T, complexed cardiac troponin T, free and complexed cardiac troponin T, total cardiac troponin, annexin V, B-enolase, CK-MB, glycogen phosphorylase-BB, heart type fatty acid binding protein, phosphoglyceric acid mutase, and S-100ao, and one or more markers related to pulmonary injury, preferably selected from the group consisting of neutrophil elastase, KL-6, LAMP 3, LAMP3, pulmonary surfactant protein A, pulmonary surfactant protein B, pulmonary surfactant protein C, pulmonary surfactant protein D, phospholipase D, PLA2G5, SFTPC, HTI56, and HTII280, or markers related thereto.
As discussed herein, these markers may be measured at a single time point, and/or may be measured at multiple time points for calculation of a change in the marker level(s) over time.
In related aspects, the present invention relates to methods for identifying marker panels for use in the foregoing methods. In developing a panel of markers useful in diagnosis and/or prognosis, data for a number of potential markers may be obtained from a group of subjects by testing for the presence or level of certain markers. The group of subjects may then be divided into sets. For example, a first set includes subjects who have been confirmed as having a disease or, more generally, being in a first condition state. The confirmation of this condition state may be made through a more rigorous and/or expensive testing. A second set of subjects is selected from those who do not fall within the first set.
The data obtained from subjects in these sets includes levels of a plurality of markers. Preferably, data for the same set of markers is available for each patient. Exemplary markers are described herein. Actual known relevance of the marker(s) to the disease of interest is not required. Methods for comparing these subject sets for relevance of one or more markers is described hereinafter. Embodiments of the methods and systems described herein may be used to determine which of the candidate markers are most relevant to the diagnosis of the disease or condition or of a given prognosis.
In another aspect, the invention relates to methods for determining a treatment regimen for use in a subject exhibiting VTED, PE, and/or DVT. The methods preferably comprise performing the methods described herein to rule in or out VTED, PE, and/or DVT; and/or to assign a prognosis to a subject diagnosed with VTED, PE, and/or DVT. One or more treatment regimens can then be selected, based on the condition and/or prognosis assigned to the subject.
In a further aspect, the invention relates to kits to rule in or out VTED, PE, and/or DVT; and or to assign a prognosis to a subject diagnosed with VTED, PE, and/or DVT. These kits preferably comprise devices and reagents for measuring a plurality of marker levels in a patient sample, and instructions for performing the assay. Optionally, the kits may contain one or more means for correlating marker level(s) in order to provide a diagnosis and/or a prognosis. Such kits preferably contain sufficient reagents to perform at least one, and preferably two or more assays described above, and/or Food and Drug Administration (FDA)-approved labeling.
In yet a further aspect, the invention relates to devices to rule in or out VTED, PE, and/or DVT; and/or to assign a prognosis to a subject diagnosed with VTED, PE, and/or DVT. Such devices preferably contain a plurality of diagnostic zones, each of which is configured to provide a signal indicative of one or more of the assay results described above. Such devices may be referred to as “arrays” or “microarrays.” Following reaction of a sample with the devices, a signal is generated from the diagnostic zone(s), which may then be correlated to the presence or amount of the markers of interest. Numerous suitable devices are known to the skilled artisan.
The present invention relates in part to methods and compositions for diagnosis, prognosis, and determination of treatment regimens in subjects. In particular, the invention relates to methods and compositions selected to rule in or out VTED, PE, and/or DVT; and or to assign a prognosis to a subject diagnosed with VTED, PE, and/or DVT.
Differential diagnosis refers to methods for diagnosing the particular disease(s) underlying the symptoms in a particular subject, based on a comparison of the characteristic features observable from the subject to the characteristic features of those potential diseases. Depending on the breadth of diseases that must be considered in the differential diagnosis, the types and number of tests that must be ordered by a clinician can be quite large. The clinician must then integrate information obtained from a battery of tests, leading to a clinical diagnosis that best represents the range of symptoms and/or diagnostic test results obtained for the subject.
Patients presenting for medical treatment often exhibit one or a few primary observable changes in bodily characteristics or functions that are indicative of disease. Often, as in the case of VTED, these “symptoms” are nonspecific, in that a number of potential diseases can present the same observable symptom or symptoms.
The present invention describes methods and compositions that can assist in the differential diagnosis of one or more nonspecific symptoms by providing diagnostic markers that are designed to rule in or out one, and preferably a plurality, of possible etiologies for the observed symptoms. Symptom-based differential diagnosis described herein can be achieved using panels of diagnostic markers designed to distinguish between possible diseases that underlie a nonspecific symptom observed in a patient.
Definitions
The term “marker” as used herein refers to proteins, polypeptides, glycoproteins, proteoglycans, lipids, lipoproteins, glycolipids, phospholipids, nucleic acids, carbohydrates, etc. or small molecules to be used as targets for screening test samples obtained from subjects. “Proteins or polypeptides” used as markers in the present invention are contemplated to include any fragments thereof, in particular, immunologically detectable fragments. Markers can also include other measurable physical characteristics including results from blood pressure measurements, temperature measurements, pulse oximetry measurements, patient history, radiography, electrocardiogram, exercise treadmill testing, blood chemistry analysis, echocardiography, bronchoprovocation testing, spirometry, pulse oximetry, esophageal pH monitoring, laryngoscopy, computed tomography, histology, cytology, magnetic resonance imaging, etc. Similarly, markers can also include clinical “scores” such as a pre-test probability assignment, a pulmonary hypertension “Daniel” score, an NIH stroke score, a Sepsis Score of Elebute and Stoner, a Duke Criteria for Infective Endocarditis, a Mannheim Peritonitis Index, an “Apache” score, etc.
Preferably, the methods described hereinafter utilize one or more markers that are derived from the subject. The term “subject-derived marker” as used herein refers to protein, polypeptide, phospholipid, nucleic acid, prion, glycoprotein, proteoglycan, glycolipid, lipid, lipoprotein, carbohydrate, or small molecule markers that are expressed or produced by one or more cells of the subject. The presence, absence, amount, or change in amount of one or more markers may indicate that a particular disease is present, or may indicate that a particular disease is absent. In the case of markers that are known to exist in membrane-bound form, such as type I membrane proteins, cell surface receptors, etc., soluble forms which may be measured in body fluids typically are generated by alternative splicing and/or cleavage and release from the plasma membrane.
Additionally, “non-subject-derived markers” may also be used. Such markers are not derived from the subject as defined herein, but rather that are characteristics of the subject observable by the artisan. Such markers are discussed above, and can include various measureable characteristics and/or clinical scores. This list is not meant to be limiting.
The term “related marker” as used herein refers to one or more fragments of a particular subject-derived marker or its biosynthetic parent that may be detected as a surrogate for the marker itself or as independent markers. For example, human BNP is derived by proteolysis of a 108 amino acid precursor molecule, referred to hereinafter as BNP1-108. Mature BNP, or “the BNP natriuretic peptide,” or “BNP-32” is a 32 amino acid molecule representing amino acids 77-108 of this precursor, which may be referred to as BNP77-108. The remaining residues 1-76 are referred to hereinafter as BNP1-76. Additionally, related markers may be the result of covalent modification of the parent marker, for example by oxidation of methionine residues, ubiquitination, cysteinylation, phosphorylation, nitrosylation, glycosylation, etc.
The sequence of the 108 amino acid BNP precursor pro-BNP (BNP1-108) is as follows, with mature BNP (BNP77-108) underlined:
BNP1-108 is synthesized as a larger precursor pre-pro-BNP having the following sequence (with the “pre” sequence shown in bold):
While mature BNP itself may be used as a marker in the present invention, the prepro-BNP, BNP1-108 and BNP1-76 molecules represent BNP-related markers that may be measured either as surrogates for mature BNP or as markers in and of themselves. In addition, one or more fragments of these molecules, including BNP-related polypeptides selected from the group consisting of BNP77-106, BNP79-106, BNP76-107, BNP69-108, BNP79-108, BNP80-108, BNP81-108, BNP83-108, BNP39-86, BNP53-85, BNP66-98, BNP30-103, BNP11-107, BNP9-106, and BNP3-108 may also be present in circulation. In addition, natriuretic peptide fragments, including BNP fragments, may comprise one or more oxidizable methionines, the oxidation of which to methionine sulfoxide or methionine sulfone produces additional BNP-related markers. See, e.g., U.S. patent Ser. No. 10/419,059, filed Apr. 17, 2003, which is hereby incorporated by reference in its entirety including all tables, figures and claims. Preferred BNP-related molecules are proBNP, NT-proBNP, BNP79-108, and BNP3-108.
Because production of marker fragments is an ongoing process that may be a function of, inter alia, the elapsed time between onset of an event triggering marker release into the tissues and the time the sample is obtained or analyzed; the elapsed time between sample acquisition and the time the sample is analyzed; the type of tissue sample at issue; the storage conditions; the quantity of proteolytic enzymes present; etc., it may be necessary to consider this degradation when both designing an assay for one or more markers, and when performing such an assay, in order to provide an accurate prognostic or diagnostic result. In addition, individual antibodies that distinguish amongst a plurality of marker fragments may be individually employed to separately detect the presence or amount of different fragments. The results of this individual detection may provide a more accurate prognostic or diagnostic result than detecting the plurality of fragments in a single assay. For example, different weighting factors may be applied to the various fragment measurements to provide a more accurate estimate of the amount of natriuretic peptide originally present in the sample.
In a similar fashion, many of the markers described herein are synthesized as larger precursor molecules, which are then processed to provide mature marker; and/or are present in circulation in the form of fragments of the marker. Thus, “related markers” to each of the markers described herein may be identified and used in an analogous fashion to that described for BNP.
Removal of polypeptide markers from the circulation often involves degradation pathways. Moreover, inhibitors of such degradation pathways may hold promise in treatment of certain diseases. See, e.g., Trindade and Rouleau, Heart Fail. Monit. 2: 2-7, 2001. However, the measurement of the polypeptide markers has focused generally upon measurement of the intact form without consideration of the degradation state of the molecules. Assays may be designed with an understanding of the degradation pathways of the polypeptide markers and the products formed during this degradation, in order to accurately measure the biologically active forms of a particular polypeptide marker in a sample. The unintended measurement of both the biologically active polypeptide marker(s) of interest and inactive fragments derived from the markers may result in an overestimation of the concentration of biologically active form(s) in a sample.
The failure to consider the degradation fragments that may be present in a clinical sample may have serious consequences for the accuracy of any diagnostic or prognostic method. Consider for example a simple case, where a sandwich immunoassay is provided for BNP, and a significant amount (e.g., 50%) of the biologically active BNP that had been present has now been degraded into an inactive form. An immunoassay formulated with antibodies that bind a region common to the biologically active BNP and the inactive fragment(s) will overestimate the amount of biologically active BNP present in the sample by 2-fold, potentially resulting in a “false positive” result. Overestimation of the biologically active form(s) present in a sample may also have serious consequences for patient management. Considering the BNP example again, the BNP concentration may be used to determine if therapy is effective (e.g., by monitoring BNP to see if an elevated level is returning to normal upon treatment). The same “false positive” BNP result discussed above may lead the physician to continue, increase, or modify treatment because of the false impression that current therapy is ineffective.
Likewise, it may be necessary to consider the complex state of one or more markers described herein. For example, troponin exists in muscle mainly as a “ternary complex” comprising three troponin polypeptides (T, I and C). But troponin I and troponin T circulate in the blood in forms other than the I/T/C ternery complex. Rather, each of (i) free cardiac-specific troponin I, (ii) binary complexes (e.g., troponin I/C complex), and (iii) ternary complexes all circulate in the blood. Furthermore, the “complex state” of troponin I and T may change over time in a patient, e.g., due to binding of free troponin polypeptides to other circulating troponin polypeptides. Immunoassays that fail to consider the “complex state” of troponin may not detect all of the cardiac-specific isoform of interest.
The term “test sample” as used herein refers to a sample of bodily fluid obtained for the purpose of diagnosis, prognosis, or evaluation of a subject of interest, such as a patient. In certain embodiments, such a sample may be obtained for the purpose of determining the outcome of an ongoing condition or the effect of a treatment regimen on a condition. Preferred test samples include blood, serum, plasma, cerebrospinal fluid, urine, saliva, sputum, and pleural effusions. In addition, one of skill in the art would realize that some test samples would be more readily analyzed following a fractionation or purification procedure, for example, separation of whole blood into serum or plasma components.
As used herein, a “plurality” as used herein refers to at least two. Preferably, a plurality refers to at least 3, more preferably at least 5, even more preferably at least 10, even more preferably at least 15, and most preferably at least 20. In particularly preferred embodiments, a plurality is a large number, i.e., at least 100.
The term “subject” as used herein refers to a human or non-human organism. Thus, the methods and compositions described herein are applicable to both human and veterinary disease. Further, while a subject is preferably a living organism, the invention described herein may be used in post-mortem analysis as well. Preferred subjects are “patients,” i.e., living humans that are receiving medical care or are being evaluated in a medical setting. This includes persons with no defined illness who are being investigated for signs of pathology.
The term “diagnosis” as used herein refers to methods by which the skilled artisan can estimate and/or determine whether or not a patient is suffering from a given disease or condition. The skilled artisan often makes a diagnosis on the basis of one or more diagnostic indicators, i.e., a marker, the presence, absence, amount, or change in amount of which is indicative of the presence, severity, or absence of the condition.
Similarly, a prognosis is often determined by examining one or more “prognostic indicators.” These are markers, the presence or amount of which in a patient (or a sample obtained from the patient) signal a probability that a given course or outcome will occur. For example, when one or more prognostic indicators reach a sufficiently high (or low) level in samples obtained from such patients, the level may signal that the patient is at an increased probability for experiencing a future outcome in comparison to a similar patient exhibiting a lower (or higher) marker level. A level or a change in level of a prognostic indicator, which in turn is associated with an increased probability of morbidity or death, is referred to as being “associated with an increased predisposition to an adverse outcome” in a patient.
The term “correlating,” as used herein in reference to the use of diagnostic and prognostic markers, refers to comparing the presence or amount of the marker(s) in a patient to its presence or amount in persons known to suffer from, or known to be at risk of, a given condition; or in persons known to be free of a given condition. As discussed above, a marker level in a patient sample can be compared to a level known to be associated with a specific diagnosis. The sample's marker level is said to have been correlated with a diagnosis; that is, the skilled artisan can use the marker level to determine whether the patient suffers from a specific type diagnosis, and respond accordingly. Alternatively, the sample's marker level can be compared to a marker level known to be associated with a good outcome (e.g., the absence of disease, etc.). In preferred embodiments, a profile of marker levels is correlated to a global probability or a particular outcome using ROC curves.
The phrase “determining the diagnosis” as used herein refers to methods by which the skilled artisan can determine the presence or absence of a particular disease in a patient. The term “diagnosis” does not refer to the ability to determine the presence or absence of a particular disease with 100% accuracy, or even that a given course or outcome is more likely to occur than not. Instead, the skilled artisan will understand that the term “diagnosis” refers to an increased probability that a certain disease is present in the subject. In preferred embodiments, a diagnosis indicates about a 5% increased chance that a disease is present, about a 10% chance, about a 15% chance, about a 20% chance, about a 25% chance, about a 30% chance, about a 40% chance, about a 50% chance, about a 60% chance, about a 75% chance, about a 90% chance, and about a 95% chance. The term “about” in this context refers to +/−2%.
The term “discrete” as used herein refers to areas of a surface that are non-contiguous. That is, two areas are discrete from one another if a border that is not part of either area completely surrounds each of the two areas.
The term “independently addressable” as used herein refers to discrete areas of a surface from which a specific signal may be obtained.
The term “antibody” as used herein refers to a peptide or polypeptide derived from, modeled after or substantially encoded by an immunoglobulin gene or immunoglobulin genes, or fragments thereof, capable of specifically binding an antigen or epitope. See, e.g. Fundamental Immunology, 3rd Edition, W. E. Paul, ed., Raven Press, N.Y. (1993); Wilson (1994) J. Immunol. Methods 175:267-273; Yarmush (1992) J. Biochem. Biophys. Methods 25:85-97. The term antibody includes antigen-binding portions, i.e., “antigen binding sites,” (e.g., fragments, subsequences, complementarity determining regions (CDRs)) that retain capacity to bind antigen, including (i) a Fab fragment, a monovalent fragment consisting of the VL, VH, CL and CHI domains; (ii) a F(ab′)2 fragment, a bivalent fragment comprising two Fab fragments linked by a disulfide bridge at the hinge region; (iii) a Fd fragment consisting of the VH and CHI domains; (iv) a Fv fragment consisting of the VL and VH domains of a single arm of an antibody, (v) a dAb fragment (Ward et al., (1989) Nature 341:544-546), which consists of a VH domain; and (vi) an isolated complementarity determining region (CDR). Single chain antibodies are also included by reference in the term “antibody.”
Preferred assays are immunoassays that provide an assay result indicative of polypeptides binding to an antibody, where such antibody specifically binds a target molecule of interest, together with those polypeptides that are related markers and that contain the epitope(s) necessary to bind to the antibody used in the assay. Assays that are “configured to detect” a target molecule of interest may, but need not, specifically detect a particular target (that is, detect detect one form of a molecule (e.g., BNP) but not a related molecule (e.g., proBNP)). Because an antibody epitope is on the order of 8 amino acids, an immunoassay will detect other polypeptides (e.g., related markers) so long as the other polypeptides contain the epitope(s) necessary to bind to the antibody used in the assay. Thus, an assay configured to detect a particular molecule (e.g., detects BNP) may also detect other related polypeptides (e.g., may also detect proBNP, together with one or more fragments of BNP or proBNP) that may exist in the sample, to the extent that such molecules contain the necessary epitopes to be detected in the assay.
The term “specifically binds” is not intended to indicate that an antibody binds exclusively to its intended target and appropriate related markers. Rather, an antibody “specifically binds” if its affinity for its intended target is about 5-fold greater when compared to its affinity for a non-target molecule. Preferably the affinity of the antibody will be at least about 5 fold, preferably 10 fold, more preferably 25-fold, even more preferably 50-fold, and most preferably 100-fold or more, greater for a target molecule than its affinity for a non-target molecule. In preferred embodiments, Specific binding between an antibody or other binding agent and an antigen means a binding affinity of at least 106 M−1. Preferred antibodies bind with affinities of at least about 107 M−1, and preferably between about 108 M−1 to about 109 M−1, about 109 M−1 to about 1010 M−1, or about 1010 M−1 to about 1011 M−1.
Affinity is calculated as Kd=koff/kon (koff is the dissociation rate constant, kon is the association rate constant and Kd is the equilibrium constant. Affinity can be determined at equilibrium by measuring the fraction bound (r) of labeled ligand at various concentrations (c). The data are graphed using the Scatchard equation: r/c=K(n−r):
where
r=moles of bound ligand/mole of receptor at equilibrium;
c=free ligand concentration at equilibrium;
K=equilibrium association constant; and
n=number of ligand binding sites per receptor molecule
By graphical analysis, r/c is plotted on the Y-axis versus r on the X-axis thus producing a Scatchard plot. The affinity is the negative slope of the line. koff can be determined by competing bound labeled ligand with unlabeled excess ligand (see, e.g., U.S. Pat. No. 6,316,409). The affinity of a targeting agent for its target molecule is preferably at least about 1×10−6 moles/liter, is more preferably at least about 1×10−7 moles/liter, is even more preferably at least about 1×10−8 moles/liter, is yet even more preferably at least about 1×10−9 moles/liter, and is most preferably at least about 1×10−10 moles/liter. Antibody affinity measurement by Scatchard analysis is well known in the art. See, e.g., van Erp et al., J. Immunoassay 12: 425-43, 1991; Nelson and Griswold, Comput. Methods Programs Biomed. 27: 65-8, 1988.
Identification of Marker Panels
In accordance with the present invention, there are provided methods and systems for the identification of one or more markers for differential diagnosis and/or risk stratification of a subject. Suitable methods for identifying markers useful for the diagnosis of disease states are described in detail in U.S. Provisional Patent Application No. 60/436,392 filed Dec. 24, 2002, PCT application US03/41426 filed Dec. 23, 2003, U.S. patent application Ser. No. 10/331,127 filed Dec. 27, 2002, and PCT application No. US03/41453, each of which is hereby incorporated by reference in its entirety, including all tables, figures, and claims.
One skilled in the art will also recognize that univariate analysis of markers can be performed and the data from the univariate analyses of multiple markers can be combined to form panels of markers to differentiate different disease conditions. Such methods include multiple linear regression, determining interaction terms, stepwise regression, neural net methods, etc.
In developing a panel of markers useful in differential diagnosis, data for a number of potential markers may be obtained from a group of subjects by testing for the presence or level of certain markers. The group of subjects is divided into two sets. The first set includes subjects who have been confirmed as having a disease or, more generally, being in a first condition state. For example, this first set of patients may be those diagnosed with VTED, PE, and/or DVT. The confirmation of this condition state may be made through a more rigorous and/or expensive testing to confirm the condition state. Hereinafter, subjects in this first set will be referred to as “diseased.”
The second set of subjects is simply those who do not fall within the first set. Subjects in this second set will hereinafter be referred to as “non-diseased”. Preferably, the first set and the second set each have an approximately equal number of subjects. This set may be normal patients, and/or patients suffering from another cardiovascular disease.
The data obtained from subjects in these sets includes levels of a plurality of markers. Preferably, data for the same set of markers is available for each patient. This set of markers may include all candidate markers that may be suspected as being relevant to the detection of a particular disease or condition. Actual known relevance is not required. Embodiments of the methods and systems described herein may be used to determine which of the candidate markers are most relevant to the diagnosis of the disease or condition. The levels of each marker in the two sets of subjects may be distributed across a broad range, e.g., as a Gaussian distribution. However, no distribution fit is required.
As noted above, a marker often is incapable of definitively identifying a patient as either diseased or non-diseased. For example, if a patient is measured as having a marker level that falls within the overlapping region, the results of the test will be useless in diagnosing the patient. An artificial cutoff may be used to distinguish between a positive and a negative test result for the detection of the disease or condition. Regardless of where the cutoff is selected, the effectiveness of the single marker as a diagnosis tool is unaffected. Changing the cutoff merely trades off between the number of false positives and the number of false negatives resulting from the use of the single marker. The effectiveness of a test having such an overlap is often expressed using a ROC (Receiver Operating Characteristic) curve. ROC curves are well known to those skilled in the art.
The horizontal axis of the ROC curve represents (1-specificity), which increases with the rate of false positives. The vertical axis of the curve represents sensitivity, which increases with the rate of true positives. Thus, for a particular cutoff selected, the value of (1-specificity) may be determined, and a corresponding sensitivity may be obtained. The area under the ROC curve is a measure of the probability that the measured marker level will allow correct identification of a disease or condition. Thus, the area under the ROC curve can be used to determine the effectiveness of the test.
As discussed above, the measurement of the level of a single marker may have limited usefulness, e.g., it may be non-specifically increased due to inflammation. The measurement of additional markers provides additional information, but the difficulty lies in properly combining the levels of two potentially unrelated measurements. In the methods and systems according to embodiments of the present invention, data relating to levels of various markers for the sets of diseased and non-diseased patients may be used to develop a panel of markers to provide a useful panel response. The data may be provided in a database such as Microsoft Access, Oracle, other SQL databases or simply in a data file. The database or data file may contain, for example, a patient identifier such as a name or number, the levels of the various markers present, and whether the patient is diseased or non-diseased.
Next, an artificial active region may be initially selected for each marker. The location of the active region may initially be selected at any point, but the selection may affect the optimization process described below. In this regard, selection near a suspected optimal location may facilitate faster convergence of the optimizer. In a preferred method, the active region is initially centered about the center of the overlap region of the two sets of patients. In one embodiment, the active region may simply be a cutoff point. In other embodiments, the active region may have a length of greater than zero. In this regard, the active region may be defined by a center value and a magnitude of length. In practice, the initial selection of the limits of the active region may be determined according to a pre-selected percentile of each set of subjects. For example, a point above which a pre-selected percentile of diseased patients are measured may be used as the right (upper) end of the cutoff range.
Each marker value for each patient may then be mapped to an indicator. The indicator is assigned one value below the active region and another value above the active region. For example, if a marker generally has a lower value for non-diseased patients and a higher value for diseased patients, a zero indicator will be assigned to a low value for a particular marker, indicating a potentially low likelihood of a positive diagnosis. In other embodiments, the indicator may be calculated based on a polynomial. The coefficients of the polynomial may be determined based on the distributions of the marker values among the diseased and non-diseased subjects.
The relative importance of the various markers may be indicated by a weighting factor. The weighting factor may initially be assigned as a coefficient for each marker. As with the active region, the initial selection of the weighting factor may be selected at any acceptable value, but the selection may affect the optimization process. In this regard, selection near a suspected optimal location may facilitate faster convergence of the optimizer. In a preferred method, acceptable weighting coefficients may range between zero and one, and an initial weighting coefficient for each marker may be assigned as 0.5. In a preferred embodiment, the initial weighting coefficient for each marker may be associated with the effectiveness of that marker by itself. For example, a ROC curve may be generated for the single marker, and the area under the ROC curve may be used as the initial weighting coefficient for that marker.
Next, a panel response may be calculated for each subject in each of the two sets. The panel response is a function of the indicators to which each marker level is mapped and the weighting coefficients for each marker. In a preferred embodiment, the panel response (R) for each subject (j) is expressed as:
Rj=ΣwiIi,j,
where i is the marker index, j is the subject index, wi is the weighting coefficient for marker i, I is the indicator value to which the marker level for marker i is mapped for subject j, and Σ is the summation over all candidate markers i. The value “R” may be referred to as a “panel index.”
One advantage of using an indicator value rather than the marker value is that an extraordinarily high or low marker levels do not change the probability of a diagnosis of diseased or non-diseased for that particular marker. Typically, a marker value above a certain level generally indicates a certain condition state. Marker values above that level indicate the condition state with the same certainty. Thus, an extraordinarily high marker value may not indicate an extraordinarily high probability of that condition state. The use of an indicator which is constant on one side of the active region eliminates this concern.
The panel response may also be a general function of several parameters including the marker levels and other factors including, for example, race and gender of the patient, various clinical characteristics, clinical “scores” such as a pulmonary hypertension “Daniel” score, an NIH stroke score, a Sepsis Score of Elebute and Stoner, a Duke Criteria for Infective Endocarditis, a Mannheim Peritonitis Index, an “Apache” score, etc. Other factors contributing to the panel response may include the slope of the value of a particular marker over time. For example, a patient may be measured when first arriving at the hospital for a particular marker. The same marker may be measured again an hour later, and the level of change may be reflected in the panel response. Further, additional markers may be derived from other markers and may contribute to the value of the panel response. For example, the ratio of values of two markers may be a factor in calculating the panel response.
Having obtained panel responses for each subject in each set of subjects, the distribution of the panel responses for each set may now be analyzed. An objective function may be defined to facilitate the selection of an effective panel. The objective function should generally be indicative of the effectiveness of the panel, as may be expressed by, for example, overlap of the panel responses of the diseased set of subjects and the panel responses of the non-diseased set of subjects. In this manner, the objective function may be optimized to maximize the effectiveness of the panel by, for example, minimizing the overlap.
In a preferred embodiment, the ROC curve representing the panel responses of the two sets of subjects may be used to define the objective function. For example, the objective function may reflect the area under the ROC curve. By maximizing the area under the curve, one may maximize the effectiveness of the panel of markers. In other embodiments, other features of the ROC curve may be used to define the objective function. For example, the point at which the slope of the ROC curve is equal to one may be a useful feature. In other embodiments, the point at which the product of sensitivity and specificity is a maximum, sometimes referred to as the “knee,” may be used. In an embodiment, the sensitivity at the knee may be maximized. In further embodiments, the sensitivity at a predetermined specificity level may be used to define the objective function. Other embodiments may use the specificity at a predetermined sensitivity level may be used. In still other embodiments, combinations of two or more of these ROC-curve features may be used.
It is possible that one of the markers in the panel is specific to the disease or condition being diagnosed. When such markers are present at above or below a certain threshold, the panel response may be set to return a “positive” test result. When the threshold is not satisfied, however, the levels of the marker may nevertheless be used as possible contributors to the objective function.
An optimization algorithm may be used to maximize or minimize the objective function. Optimization algorithms are well-known to those skilled in the art and include several commonly available minimizing or maximizing functions including the Simplex method and other constrained optimization techniques. It is understood by those skilled in the art that some minimization functions are better than others at searching for global minimums, rather than local minimums. In the optimization process, the location and size of the active region for each marker may be allowed to vary to provide at least two degrees of freedom per marker. Such variable parameters are referred to herein as independent variables. In a preferred embodiment, the weighting coefficient for each marker is also allowed to vary across iterations of the optimization algorithm. In various embodiments, any permutation of these parameters may be used as independent variables.
In addition to the above-described parameters, the sense of each marker may also be used as an independent variable. For example, in many cases, it may not be known whether a higher level for a certain marker is generally indicative of a diseased state or a non-diseased state. In such a case, it may be useful to allow the optimization process to search on both sides. In practice, this may be implemented in several ways. For example, in one embodiment, the sense may be a truly separate independent variable that may be flipped between positive and negative by the optimization process. Alternatively, the sense may be implemented by allowing the weighting coefficient to be negative.
The optimization algorithm may be provided with certain constraints as well. For example, the resulting ROC curve may be constrained to provide an area-under-curve of greater than a particular value. ROC curves having an area under the curve of 0.5 indicate complete randomness, while an area under the curve of 1.0 reflects perfect separation of the two sets. Thus, a minimum acceptable value, such as 0.75, may be used as a constraint, particularly if the objective function does not incorporate the area under the curve. Other constraints may include limitations on the weighting coefficients of particular markers. Additional constraints may limit the sum of all the weighting coefficients to a particular value, such as 1.0.
The iterations of the optimization algorithm generally vary the independent parameters to satisfy the constraints while minimizing or maximizing the objective function. The number of iterations may be limited in the optimization process. Further, the optimization process may be terminated when the difference in the objective function between two consecutive iterations is below a predetermined threshold, thereby indicating that the optimization algorithm has reached a region of a local minimum or a maximum.
Thus, the optimization process may provide a panel of markers including weighting coefficients for each marker and active regions for the mapping of marker values to indicators. Certain markers may be then be changed or even eliminated from the panel, and the process repeated until a satisfactory result is obtained. The effective contribution of each marker in the panel may be determined to identify the relative importance of the markers. In one embodiment, the weighting coefficients resulting from the optimization process may be used to determine the relative importance of each marker. The markers with the lowest coefficients may be eliminated or replaced.
In certain cases, the lower weighting coefficients may not be indicative of a low importance. Similarly, a higher weighting coefficient may not be indicative of a high importance. For example, the optimization process may result in a high coefficient if the associated marker is irrelevant to the diagnosis. In this instance, there may not be any advantage that will drive the coefficient lower. Varying this coefficient may not affect the value of the objective function.
To allow a determination of test accuracy, a “gold standard” test criterion may be selected which allows selection of subjects into two or more groups for comparison by the foregoing methods. Measures of test accuracy may be obtained as described in Fischer et al., Intensive Care Med. 29: 1043-51, 2003, and used to determine the effectiveness of a given marker or panel of markers. These measures include sensitivity and specificity, predictive values, likelihood ratios, diagnostic odds ratios, and ROC curve areas. As discussed above, suitable tests may exhibit one or more of the following results on these various measures:
at least 75% sensitivity, combined with at least 75% specificity;
ROC curve area of at least 0.7, more preferably at least 0.8, even more preferably at least 0.9, and most preferably at least 0.95; and/or
a positive likelihood ratio (calculated as sensitivity/(1-specificity)) of at least 5, more preferably at least 10, and most preferably at least 20, and a negative likelihood ratio (calculated as (1-sensitivity)/specificity) of less than or equal to 0.3, more preferably less than or equal to 0.2, and most preferably less than or equal to 0.1.
Exemplary Markers
In a preferred embodiment, the following discussion provides Swiss-Prot accession numbers for the human precursor of certain markers described herein. Additional markers are described by their common names hereinafter.
In addition to the use of markers individually (“univariately), a panel consisting of the markers referenced herein may be constructed to provide relevant information related to the differential diagnosis of interest. Such a panel may be constructed using 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 or 20 individual markers. The analysis of a single marker or subsets of markers comprising a larger panel of markers could be carried out by one skilled in the art to optimize clinical sensitivity or specificity in various clinical settings. These include, but are not limited to ambulatory, urgent care, critical care, intensive care, monitoring unit, inpatient, outpatient, physician office, medical clinic, and health screening settings. Furthermore, one skilled in the art can use a single marker or a subset of markers comprising a larger panel of markers in combination with an adjustment of the diagnostic threshold in each of the aforementioned settings to optimize clinical sensitivity and specificity. The following provides a brief discussion of additional exemplary markers for use in identifying suitable marker panels by the methods described herein.
Additional Markers
A panel consisting of the markers referenced herein and/or their related markers may be constructed to provide relevant information related to the diagnosis of interest. Such a panel may be constructed using 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, or more individual markers. The analysis of a single marker or subsets of markers comprising a larger panel of markers could be carried out by one skilled in the art to optimize clinical sensitivity or specificity in various clinical settings. These include, but are not limited to ambulatory, urgent care, critical care, intensive care, monitoring unit, inpatient, outpatient, physician office, medical clinic, and health screening settings. Furthermore, one skilled in the art can use a single marker or a subset of markers comprising a larger panel of markers in combination with an adjustment of the diagnostic threshold in each of the aforementioned settings to optimize clinical sensitivity and specificity.
The following table provides a list of additional preferred markers for use in the present invention. Further detail is provided in US2005/0148029, which is hereby incorporated by reference in its entirety. As described herein, markers related to each of these markers are also encompassed by the present invention. Classification of markers refers to the following categories: Blood pressure regulation: markers related to blood pressure regulation; Inflammation: markers related to inflammation; Apoptosis: markers related to apoptosis; Reactive oxygen: markers related to reactive oxygen species; Myocardial injury: markers related to myocardial injury; Pulmonary injury: markers related to pulmonary injury; Coagulation and hemostasis: markers related to coagulation and hemostasis; Vascular tissue: markers related to vascular tissue injury; Neural tissue injury: markers related to neural tissue injury; Collagen synthesis & degradation: markers related to collagen synthesis and degradation; Tissue injury: markers related to general tissue injury.
Ubiquitination of Markers
Ubiquitin-mediated degradation of proteins plays an important role in the control of numerous processes, such as the way in which extracellular materials are incorporated into a cell, the movement of biochemical signals from the cell membrane, and the regulation of cellular functions such as transcriptional on-off switches. The ubiquitin system has been implicated in the immune response and development. Ubiquitin is a 76-amino acid polypeptide that is conjugated to proteins targeted for degradation. The ubiquitin-protein conjugate is recognized by a 26S proteolytic complex that splits ubiquitin from the protein, which is subsequently degraded.
It has been reported that sepsis stimulates protein breakdown in skeletal muscle by a nonlysosomal energy-dependent proteolytic pathway, and because muscle levels of ubiquitin mRNA were also increased, the results were interpreted as indicating that sepsis-induced muscle protein breakdown is caused by upregulated activity of the energy-ubiquitin-dependent proteolytic pathway. The same proteolytic pathway has been implicated in muscle breakdown caused by denervation, fasting, acidosis, cancer, and burn injury. Thus, levels of ubiquitinated proteins generally, or of specific ubiquitin-protein conjugates or fragments thereof, can be measured as additional markers of the invention. See, Tiao et al., J. Clin. Invest. 99: 163-168, 1997. Moreover, circulating levels of ubiquitin itself can be a useful marker in the methods described herein. See, e.g., Majetschak et al., Blood 101: 1882-90, 2003.
The skilled artisan will recognize that an assay for ubiquitin may be designed that recognizes ubiquitin itself, ubiquitin-protein conjugates, or both ubiquitin and ubiquitin-protein conjugates. For example, antibodies used in a sandwich immunoassay may be selected so that both the solid phase antibody and the labeled antibody recognize a portion of ubiquitin that is available for binding in both unconjugated ubiquitin and ubiquitin conjugates. Alternatively, an assay specific for ubiquitin conjugates of the muscle protein troponin could use one antibody (on a solid phase or label) that recognizes ubiquitin, and a second antibody (the other of the solid phase or label) that recognizes troponin.
The present invention contemplates measuring ubiquitin conjugates of any marker described herein. Preferred ubiquitin-muscle protein conjugates for detection as markers include, but are not limited to, troponin I-ubiquitin, troponin T-ubiquitin, troponin C-ubiquitin, binary and ternary troponin complex-ubiquitin, actin-ubiquitin, myosin-ubiquitin, tropomyosin-ubiquitin, and α-actinin-ubiquitin.
Assay Measurement Strategies
Numerous methods and devices are well known to the skilled artisan for the detection and analysis of the markers of the instant invention. With regard to polypeptides or proteins in patient test samples, immunoassay devices and methods are often used. See, e.g., U.S. Pat. Nos. 6,143,576; 6,113,855; 6,019,944; 5,985,579; 5,947,124; 5,939,272; 5,922,615; 5,885,527; 5,851,776; 5,824,799; 5,679,526; 5,525,524; and 5,480,792, each of which is hereby incorporated by reference in its entirety, including all tables, figures and claims. These devices and methods can utilize labeled molecules in various sandwitch, competitive, or non-competitive assay formats, to generate a signal that is related to the presence or amount of an analyte of interest. Additionally, certain methods and devices, such as biosensors and optical immunoassays, may be employed to determine the presence or amount of analytes without the need for a labeled molecule. See, e.g., U.S. Pat. Nos. 5,631,171; and 5,955,377, each of which is hereby incorporated by reference in its entirety, including all tables, figures and claims. One skilled in the art also recognizes that robotic instrumentation including but not limited to Beckman Access, Abbott AxSym, Roche ElecSys, Dade Behring Stratus systems are among the immunoassay analyzers that are capable of performing the immunoassays taught herein.
Preferably the markers are analyzed using an immunoassay, and most preferably sandwich immunoassay, although other methods are well known to those skilled in the art (for example, the measurement of marker RNA levels). The presence or amount of a marker is generally determined using antibodies specific for each marker and detecting specific binding. Any suitable immunoassay may be utilized, for example, enzyme-linked immunoassays (ELISA), radioimmunoassays (RIAs), competitive binding assays, and the like. Specific immunological binding of the antibody to the marker can be detected directly or indirectly. Direct labels include fluorescent or luminescent tags, metals, dyes, radionuclides, and the like, attached to the antibody. Indirect labels include various enzymes well known in the art, such as alkaline phosphatase, horseradish peroxidase and the like.
The use of immobilized antibodies specific for the markers is also contemplated by the present invention. The antibodies could be immobilized onto a variety of solid supports, such as magnetic or chromatographic matrix particles, the surface of an assay place (such as microtiter wells), pieces of a solid substrate material or membrane (such as plastic, nylon, paper), and the like. An assay strip could be prepared by coating the antibody or a plurality of antibodies in an array on solid support. This strip could then be dipped into the test sample and then processed quickly through washes and detection steps to generate a measurable signal, such as a colored spot.
For separate or sequential assay of markers, suitable apparatuses include clinical laboratory analyzers such as the ElecSys (Roche), the AxSym (Abbott), the Access (Beckman), the ADVIA® CENTAUR® (Bayer) immunoassay systems, the NICHOLS ADVANTAGE® (Nichols Institute) immunoassay system, etc. Preferred apparatuses perform simultaneous assays of a plurality of markers using a single test device. Particularly useful physical formats comprise surfaces having a plurality of discrete, adressable locations for the detection of a plurality of different analytes. Such formats include protein microarrays, or “protein chips” (see, e.g., Ng and Ilag, J. Cell Mol. Med. 6: 329-340 (2002)) and certain capillary devices (see, e.g., U.S. Pat. No. 6,019,944). In these embodiments, each discrete surface location may comprise antibodies to immobilize one or more analyte(s) (e.g., a marker) for detection at each location. Surfaces may alternatively comprise one or more discrete particles (e.g., microparticles or nanoparticles) immobilized at discrete locations of a surface, where the microparticles comprise antibodies to immobilize one analyte (e.g., a marker) for detection.
Preferred assay devices of the present invention will comprise, for one or more assays, a first antibody conjugated to a solid phase and a second antibody conjugated to a signal development element. Such assay devices are configured to perform a sandwich immunoassay for one or more analytes. These assay devices will preferably further comprise a sample application zone, and a flow path from the sample application zone to a second device region comprising the first antibody conjugated to a solid phase.
Flow of a sample along the flow path may be driven passively (e.g., by capillary, hydrostatic, or other forces that do not require further manipulation of the device once sample is applied), actively (e.g., by application of force generated via mechanical pumps, electroosmotic pumps, centrifugal force, increased air pressure, etc.), or by a combination of active and passive driving forces. Most preferably, sample applied to the sample application zone will contact both a first antibody conjugated to a solid phase and a second antibody conjugated to a signal development element along the flow path (sandwich assay format). Additional elements, such as filters to separate plasma or serum from blood, mixing chambers, etc., may be included as required by the artisan. Exemplary devices are described in Chapter 41, entitled “Near Patient Tests: Triage® Cardiac System,” in The Immunoassay Handbook, 2nd ed., David Wild, ed., Nature Publishing Group, 2001, which is hereby incorporated by reference in its entirety.
A panel consisting of the markers referenced above may be constructed to provide relevant information related to differential diagnosis. Such a panel may be constructed using 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, or more or individual markers. The analysis of a single marker or subsets of markers comprising a larger panel of markers could be carried out by one skilled in the art to optimize clinical sensitivity or specificity in various clinical settings. These include, but are not limited to ambulatory, urgent care, critical care, intensive care, monitoring unit, inpatient, outpatient, physician office, medical clinic, and health screening settings. Furthermore, one skilled in the art can use a single marker or a subset of markers comprising a larger panel of markers in combination with an adjustment of the diagnostic threshold in each of the aforementioned settings to optimize clinical sensitivity and specificity. The clinical sensitivity of an assay is defined as the percentage of those with the disease that the assay correctly predicts, and the specificity of an assay is defined as the percentage of those without the disease that the assay correctly predicts (Tietz Textbook of Clinical Chemistry, 2nd edition, Carl Burtis and Edward Ashwood eds., W.B. Saunders and Company, p. 496).
The analysis of markers could be carried out in a variety of physical formats as well. For example, the use of microtiter plates or automation could be used to facilitate the processing of large numbers of test samples. Alternatively, single sample formats could be developed to facilitate immediate treatment and diagnosis in a timely fashion, for example, in ambulatory transport or emergency room settings.
In another embodiment, the present invention provides a kit for the analysis of markers. Such a kit preferably comprises devises and reagents for the analysis of at least one test sample and instructions for performing the assay. Optionally the kits may contain one or more means for using information obtained from immunoassays performed for a marker panel to rule in or out certain diagnoses. Other measurement strategies applicable to the methods described herein include chromatography (e.g., HPLC), mass spectrometry, receptor-based assays, and combinations of the foregoing.
Selection of Antibodies
The generation and selection of antibodies may be accomplished several ways. For example, one way is to purify polypeptides of interest or to synthesize the polypeptides of interest using, e.g., solid phase peptide synthesis methods well known in the art. See, e.g., Guide to Protein Purification, Murray P. Deutcher, ed., Meth. Enzymol. Vol 182 (1990); Solid Phase Peptide Synthesis, Greg B. Fields ed., Meth. Enzymol. Vol 289 (1997); Kiso et al., Chem. Pharm. Bull. (Tokyo) 38: 1192-99, 1990; Mostafavi et al., Biomed. Pept. Proteins Nucleic Acids 1: 255-60, 1995; Fujiwara et al., Chem. Pharm. Bull. (Tokyo) 44: 1326-31, 1996. The selected polypeptides may then be injected, for example, into mice or rabbits, to generate polyclonal or monoclonal antibodies. One skilled in the art will recognize that many procedures are available for the production of antibodies, for example, as described in Antibodies, A Laboratory Manual, Ed Harlow and David Lane, Cold Spring Harbor Laboratory (1988), Cold Spring Harbor, N.Y. One skilled in the art will also appreciate that binding fragments or Fab fragments which mimic antibodies can also be prepared from genetic information by various procedures (Antibody Engineering: A Practical Approach (Borrebaeck, C., ed.), 1995, Oxford University Press, Oxford; J. Immunol. 149, 3914-3920 (1992)).
In addition, numerous publications have reported the use of phage display technology to produce and screen libraries of polypeptides for binding to a selected target. See, e.g., Cwirla et al., Proc. Natl. Acad. Sci. USA 87, 6378-82, 1990; Devlin et al., Science 249, 404-6, 1990, Scott and Smith, Science 249, 386-88, 1990; and Ladner et al., U.S. Pat. No. 5,571,698. A basic concept of phage display methods is the establishment of a physical association between DNA encoding a polypeptide to be screened and the polypeptide. This physical association is provided by the phage particle, which displays a polypeptide as part of a capsid enclosing the phage genome which encodes the polypeptide. The establishment of a physical association between polypeptides and their genetic material allows simultaneous mass screening of very large numbers of phage bearing different polypeptides. Phage displaying a polypeptide with affinity to a target bind to the target and these phage are enriched by affinity screening to the target. The identity of polypeptides displayed from these phage can be determined from their respective genomes. Using these methods a polypeptide identified as having a binding affinity for a desired target can then be synthesized in bulk by conventional means. See, e.g., U.S. Pat. No. 6,057,098, which is hereby incorporated in its entirety, including all tables, figures, and claims.
The antibodies that are generated by these methods may then be selected by first screening for affinity and specificity with the purified polypeptide of interest and, if required, comparing the results to the affinity and specificity of the antibodies with polypeptides that are desired to be excluded from binding. The screening procedure can involve immobilization of the purified polypeptides in separate wells of microtiter plates. The solution containing a potential antibody or groups of antibodies is then placed into the respective microtiter wells and incubated for about 30 min to 2 h. The microtiter wells are then washed and a labeled secondary antibody (for example, an anti-mouse antibody conjugated to alkaline phosphatase if the raised antibodies are mouse antibodies) is added to the wells and incubated for about 30 min and then washed. Substrate is added to the wells and a color reaction will appear where antibody to the immobilized polypeptide(s) are present.
The antibodies so identified may then be further analyzed for affinity and specificity in the assay design selected. In the development of immunoassays for a target protein, the purified target protein acts as a standard with which to judge the sensitivity and specificity of the immunoassay using the antibodies that have been selected. Because the binding affinity of various antibodies may differ; certain antibody pairs (e.g., in sandwich assays) may interfere with one another sterically, etc., assay performance of an antibody may be a more important measure than absolute affinity and specificity of an antibody.
Those skilled in the art will recognize that many approaches can be taken in producing antibodies or binding fragments and screening and selecting for affinity and specificity for the various polypeptides, but these approaches do not change the scope of the invention.
Selecting a Treatment Regimen
Just as the potential causes of any particular nonspecific symptom may be a large and diverse set of conditions, the appropriate treatments for these potential causes may be equally large and diverse. However, once a diagnosis is obtained, the clinician can readily select a treatment regimen that is compatible with the diagnosis. The skilled artisan is aware of appropriate treatments for numerous diseases discussed in relation to the methods of diagnosis described herein. See, e.g., Merck Manual of Diagnosis and Therapy, 17th Ed. Merck Research Laboratories, Whitehouse Station, N.J., 1999.
In addition, since the methods and compositions described herein provide prognostic information, the panels and markers of the present invention may be used to monitor a course of treatment. For example, inproved or worsened prognostic state may indicate that a particular treatment is or is not efficacious.
The following examples serve to illustrate the present invention. These examples are in no way intended to limit the scope of the invention.
Blood specimens are collected by trained study personnel using EDTA as the anticoagulant and centrifuged for greater than or equal to 10 minutes. The plasma component is transferred into a sterile cryovial and frozen at −20° C. or colder. Clinical histories are available for each of the patients to aid in the statistical analysis of the assay data.
For Examples 6, 7, 8, and 9 below, samples were obtained from 120 normal healthy controls. A study population consisted of prospectively enrolled inpatients and emergency department patients with confirmed pulmonary embolism. The inclusion criterion that triggered a review for eligibility was the presumed diagnosis of pulmonary embolism based upon the initial interpretation of either the computerized tomography (CT) chest angiography scan or ventilation-perfusion (V/Q) lung scan. Exclusion was based upon the following criteria: (1)>12 hours since start of heparin therapy, (2) systolic hypotension (<100 mm Hg for two consecutive measurements obtained greater than 15 min apart, (3) initial treatment with fibrinolytic therapy, catheter fragmentation, or surgical embolectomy, (4) illnesses with a predicted 6-month mortality >50% (e.g., metastatic cancer, end-stage AIDs, end-stage heart or renal failure with no plan for transplantation or hemodialysis therapy), (5) clinical plan to not treat the patient, (6) anatomy or clinical condition that precluded the ability to measure right ventricular function via echocardiography, (7). permanent inability to walk for any reason, (8) prisoners with expected incarceration time >6 months, (9) personal physicians' refusal, (10) inability or unwillingness to participate in the informed consent process. Exclusions that were applied after informed consent were: (1) reinterpretation of CT angiography as negative for pulmonary embolism and no deep venous thrombosis observed on both CT venography and venous ultrasonography, (2) inability to perform echocardiography, (3) patient withdrawal from the study, or loss to follow-up in survivors.
Samples were obtained upon enrollment for 100 subjects diagnosed with PE and/or DVT, and 37 subjects diagnosed with PE and/or DVT and later experiencing an adverse outcome (defined as death, shock, intubation, surgical thrombectomy during hospitalization, or severe disability defined as dilated right ventricle, moderate or severe hypokinesis, sPAP>35 and either rest dyspnea or upon walking <330 m, over 6 months of follow up).
Individual assays were performed as described below in Examples 4 and 5. Assays were configured to bind the following markers, and results are reported in the following examples using the following units (alternate names and units of measurement are in parenthesis): acidic calponin (ng/mL); adrenomedullin (pg/ml); angiopoietin-4 (ng/ml); basic calponin (ng/ml); bone morphogenetic protein 4 (BMP4) (ng/ml); B-type natriuretic peptide (BNP, BNP77-108) (pg/ml); BNP1-108 (proBNP) (pg/ml); BNP3-108 (pg/ml); BNP79-108 (pg/ml); caspase-3 (ng/mL); CCL-11 (pg/ml); calcitonin gene-related peptide (CGRP) (pg/ml); creatine kinase-BB (CK-BB) (ng/ml); creatine kinase-MB (CK-MB) (ng/ml); C-reactive protein (CRP) (μg/ml); D-dimer (ng/ml); soluble elastin fragments (sELAF) (ng/ml); endothelin-1 (pg/ml); glutathione-S-transferase 3 (GSTP) (ng/ml); heart-type fatty acid binding protein (hFABP) (ng/ml); interleukin 1 receptor antagonist (IL-1ra) (pg/ml); interleukin-25 (IL-25) (pg/ml); leptin (ng/ml); soluble lymphotoxin b receptor (sLTBR, sTNFRSF3) (ng/ml); monocyte chemotactic protein-1 (MCP-1) (pg/ml); matrix metalloproteinase-9 (MMP9) (ng/ml); myloperoxidase (MPO) (ng/ml); myoglobin (MYO) (ng/ml); nucleoside diphosphate kinase A (NDKA) (ng/ml); neuropilin-2 (ng/ml); neutrophil gelatinase-associated lipocalin (NGAL) (ng/ml); placental growth factor-1 (PLGF-1) (pg/ml); placental growth factor 1 and 2 (PLGF-1/2) (pg/ml); ANP28-151 (des-Asn-Pro proANP) (ng/ml); activated protein C (ng/ml); total protein C (latent+active PC) (ng/ml); pulmonary surfactant protein A (PSAP-A) (ng/ml); pulmonary surfactant protein B (PSAP-B) (ng/ml); pulmonary surfactant protein C (PSAP-C) (ng/ml); pulmonary surfactant protein D (PSAP-D) (ng/ml); soluble receptor for advanced glycosylation end products (sRAGE) (ng/ml); soluble platelet endothelial cell adhesion molecule-1 (sPECAM-1) (ng/ml); spectrin alpha chain, 120 kda (spectrin 120) (pg/ml); spectrin alpha chain, 145 kda (spectrin 145) (ng/ml); angiopoietin-1 receptor (TIE-2) (ng/ml); tissue factor (pg/ml); soluble tumor necrosis factor receptor 1a (sTNFR1a) (ng/ml); soluble tumor necrosis factor receptor superfamily member 7 (sTNFRSF7, cd27) (ng/ml); soluble tumor necrosis factor receptor superfamily member 14 (TNFsR14) (ng/ml); cardiac troponin I, free and complexed (cTNI) (ng/ml); thrombus precursor protein (TpP) (μg/ml); ubiquitin fusion degradation protein 1 homolog (UFDP1H) (ng/ml); urokinase-type plasminogen activator (UPA) (ng/ml); vascular cell adhesion protein 1 (VCAM-1) (ng/ml); VE-cadherin (cd144 antigen) (ng/ml); vascular endothelial growth factor (VEGF) (pg/ml); soluble flt-1 (sVEGF-R1) (pg/ml); soluble kdr (sVEGF-R2) ng/ml; von Willebrand factor containing a VWF integrin domain (VWF-integrin) (ng/ml).
For the sandwich immunoassay in microtiter plates, a monoclonal antibody directed against a selected analyte was biotinylated using N-hydroxysuccinimide biotin (NHS-biotin) at a ratio of about 5 NHS-biotin moieties per antibody. The antibody-biotin conjugate was then added to wells of a standard avidin 384 well microtiter plate, and antibody conjugate not bound to the plate was removed. This formed the “anti-marker” in the microtiter plate. Another monoclonal antibody directed against the same analyte was conjugated to alkaline phosphatase, for example using succinimidyl 4-[N-maleimidomethyl]-cyclohexane-1-carboxylate (SMCC) and N-succinimidyl 3-[2-pyridyldithio]propionate (SPDP) (Pierce, Rockford, Ill.).
Biotinylated antibodies were pipetted into microtiter plate wells previously coated with avidin and incubated for 60 min. The solution containing unbound antibody was removed, and the wells washed with a wash buffer, consisting of 20 mM borate (pH 7.42) containing 150 mM NaCl, 0.1% sodium azide, and 0.02% Tween-20. The plasma samples (10 μL, or 20 μL for CCL4) containing added HAMA inhibitors were pipeted into the microtiter plate wells, and incubated for 60 min. The sample was then removed and the wells washed with a wash buffer. The antibody-alkaline phosphatase conjugate was then added to the wells and incubated for an additional 60 min, after which time, the antibody conjugate was removed and the wells washed with a wash buffer. A substrate, (AttoPhos®, Promega, Madison, Wis.) was added to the wells, and the rate of formation of the fluorescent product is related to the concentration of the analyte in the sample tested.
For competitive immunoassays in microtiter plates, a murine monoclonal antibody directed against a selected analyte was added to the wells of a microtiter plate and immobilized by binding to goat anti-mouse antibody that is pre-absorbed to the surface of the microtiter plate wells (Pierce, Rockford, Ill.). Any unbound murine monoclonal antibody was removed after a 60 minute incubation. This forms the “anti-marker” in the microtiter plate. A purified polypeptide that is either the same as or related to the selected analyte, and that can be bound by the monoclonal antibody, was biotinylated as described above for the biotinylation of antibodies. This biotinylated polypeptide was mixed with the sample in the presence of HAMA inhibitors, forming a mixture containing both exogenously added biotinylated polypeptide and any unlabeled analyte molecules endogenous to the sample. The amount of the monoclonal antibody and biotinylated marker added depends on various factors and was titrated empirically to obtain a satisfactory dose-response curve for the selected analyte.
This mixture was added to the microtiter plate and allowed to react with the murine monoclonal antibody for 120 minutes. After the 120 minute incubation, the unbound material was removed, and Neutralite-Alkaline Phosphatase (Southern Biotechnology; Birmingham, Ala.) was added to bind to any immobilized biotinylated polypeptide. Substrate (as described above) was added to the wells, and the rate of formation of the fluorescent product was related to the amount of biotinylated polypeptide bound, and therefore was inversely related to the endogenous amount of the analyte in the specimen.
Immunoassays were performed using microfluidic devices essentially as described in Chapter 41, entitled “Near Patient Tests: Triage® Cardiac System,” in The Immunoassay Handbook, 2nd ed., David Wild, ed., Nature Publishing Group, 2001.
For sandwich immunoassays, a plasma sample is added to the microfluidic device that contains all the necessary assay reagents, including HAMA inhibitors, in dried form. The plasma passes through a filter to remove particulate matter. Plasma enters a “reaction chamber” by capillary action. This reaction chamber contains fluorescent latex particle-antibody conjugates (hereafter called FETL-antibody conjugates) appropriate to an analyte of interest, and may contain FETL-antibody conjugates to several selected analytes. The FETL-antibody conjugates dissolve into the plasma to form a reaction mixture, which is held in the reaction chamber for an incubation period (about a minute) to allow the analyte(s) of interest in the plasma to bind to the antibodies. After the incubation period, the reaction mixture moves down the detection lane by capillary action. Antibodies to the analyte(s) of interest are immobilized in discrete capture zones on the surface of a “detection lane.” Analyte/antibody-FETL complexes formed in the reaction chamber are captured on an appropriate detection zone to form a sandwich complex, while unbound FETL-antibody conjugates are washed from the detection lane into a waste chamber by excess plasma. The amount of analyte/antibody-FETL complex bound on a capture zone is quantified with a fluorometer (Triage® MeterPlus, Biosite Incorporated) and is related to the amount of the selected analyte in the plasma specimen.
For competitive immunoassays, the procedure and process is similar to that described for sandwich immunoassays, with the following exceptions. In one configuration, fluorescent latex particle-marker (FETL-marker) conjugates are provided in the reaction chamber, and are dissolved in the plasma to form a reaction mixture. This reaction mixture contains both the unlabeled analyte endogenous to the sample, and the FETL-marker conjugates. When the reaction mixture contacts the capture zone for a analyte of interest, the unlabeled endogenous analyte and the FETL-marker conjugates compete for the limited number of antibody binding sites. Thus, the amount of FETL-marker conjugate bound to the capture zone is inversely related to the amount of analyte endogenously present in the plasma specimen. In another configuration, antibody-FETL conjugates are provided in the reaction chamber as described above for sandwich assays. In this configuration, the capture zone contains immobilized marker on the surface of the detection lane. Free antibody-FETL conjugates bind to this immobilized marker on the capture zone, while antibody-FETL conjugates bound to an analyte of interest do not bind as readily or at all to this immobilized marker. Again, the amount of FETL captured in the zone is inversely related to the amount of the selected analyte in the plasma specimen. One skilled in the art will recognize that either configuration may be used depending on the characteristics and concentrations of the selected analyte(s).
The assays were calibrated using purified proteins (that is either the same as or related to the selected analyte, and that can be detected in the assay) diluted gravimetrically into EDTA plasma treated in the same manner as the sample population specimens. Endogenous levels of the analyte present in the plasma prior to addition of the purified marker protein was measured and taken into account in assigning the marker values in the calibrators. When necessary to reduce endogenous levels in the calibrators, the endogenous analyte was stripped from the plasma using standard immunoaffinity methods. Calibrators were assayed in the same manner as the sample population specimens, and the resulting data used to construct a “dose-response” curve (assay signal as a function of analyte concentration), which may be used to determine analyte concentrations from assay signals obtained from subject specimens.
Significantly higher levels of the following biomarkers were observed in the PE and/or DVT disease group, compared to normal health controls (Non-parametric Kruskal-Wallis test for median, all p<0.001, except for MCP-1, where p<0.007). The results are presented in the following table:
The markers were also subjected to ROC analysis, with the following results:
There were also significantly higher levels of BNP, BNP3-108, and MCP-1 found in PE and/or DVT subjects who experienced an adverse outcome, compared to PE and/or DVT subjects without adverse outcome. Wilcoxon test p-values for the comparisons are 0.002, 0.004, and 0.0312, respectively.
Table 3: PE and/or DVT; normal=no adverse outcome; disease=adverse outcome as defined above in Example 2.
The markers were also subjected to ROC analysis, with the following results:
The prevalence of adverse outcome in the study group was 28%. If a patient had a history of prior cardiopulmonary disease, a positive medical panel (defined as either persistent hypoxemia defined as a pulse oximetry reading <95% while the patient breathed room air, or evidence of pulmonary hypertension from 12-lead electrocardiography defined as a Daniel score >8, or evidence of cardiac cell necrosis from the serum troponin T measurement concentration defined as >0.1 ng/mL), a D-dimer >8 ug/mL, B-type natriuretic peptide >90 pg/mL, the probability of an adverse outcome (i.e., the predictive value positive) increased to 57%, 42%, 50%, 57%, respectively. A three marker panel consisting of a combination of troponin T>0.1 ng/mL, D-dimer >8 μg/mL, or B-type natriuretic peptide >90 pg/mL yielded sensitivity, specificity and accuracy of 72% (95% CI: 57 to 84%), 68% (58 to 76%) and 69% (61 to 76%) respectively.
Using the methods described in PCT application no. US03/41426, filed Dec. 23, 2003, exemplary panels for risk stratification is VETD were identified. Starting with a large number of potential markers, an iterative procedure was applied. In this procedure, individual threshold concentrations for the markers are not used as cutoffs per se, but are used as values to which the assay values for each patient are compared and normalized. A window factor was used to calculate the minimum and maximum values above and below the cutoff. Assay values above the maximum are set to the maximum and assay values below the minimum are set to the minimum. The absolute values of the weights for the individual markers adds up to 1. A negative weight for a marker implies that the assay values for the control group are higher than those for the diseased group.
A “panel response” is calculated using the cutoff, window, and weighting factors. The panel responses for the entire population of patients and controls are subjected to ROC analysis as is commonly performed for individual markers, and a “panel response” cutoff is selected to yield the desired sensitivity and specificity for the panel. After each set of iterations, the weakest contributors to the equation may be eliminated and the iterative process started again with the reduced number of markers. This process is continued until a minimum number of markers that will still result in acceptable sensitivity and specificity of the panel is obtained.
In the present examples, the “diseased” dataset represents a population of subjects diagnosed as having PE and/or DVT and an adverse outcome as defined above in Example 3, while the “normal” dataset represents a population of subjects diagnosed as having PE and/or DVT and no adverse outcome. Samples were obtained for these subjects at hospital admission.
For Examples 11 and 12 below, samples were collected from 439 human subjects enrolled in one of two studies on presentation at an emergency medical facility. Inclusion criteria required that the subjects were suspected of having pulmonary embolism or deep vein thrombosis. Exclusion criteria included those already treated for the current episode of PE or DVT, those with active heart failure, a history of renal disease requiring dialysis, or experiencing an MI, PTCA, CABG or Unstable Angina within the past month. A population of apparently healthy indivduals not enrolled in the study (n=49) were selected to serve as normal control samples.
The diagnoses of the 439 enrolled subjects, provided by standard clinical methods in use at the medical facility, were as follows:
57 PE+ 382 PE−
51 DVT+ 381 DVT− 7 DVT status uncertain
For assessment of prognosis, subjects were followed for 45 days. The following outcomes were observed:
Among the 57 PE+ patients: 4 deaths in the follow up period; 23 combined adverse outcomes defined as ICU admissions on presentation, recurring PE in the follow up period and/or deaths in the follow up period.
Among the 439 enrolled subjects: 11 deaths in the follow up period; 34 combined adverse outcomes by the same definition.
The subject population was subdivided into the following groups (1) through (4) for diagnostic comparisons: (1) normal subjects gathered outside of study enrollment criteria; (2) PE+ (positive for pulmonary embolism, and positive, negative, or undetermined for DVT); (3) DVT+ (positive for DVT, and positive or negative for pulmonary embolism), and (4) PE− (negative for pulmonary embolism, and positive, negative, or undetermined for DVT). The subject population was further subdivided into the following groups (5) through (12) for prognostic comparisons: (5) PE+ and alive through the followup period; (6) PE+, and death within the follow up period; (7) all enrolled subjects, and alive at 45 days; (8) all enrolled subjects, and death within the follow up period; (9) PE+, and no combined adverse outcome through the follow up period; (10) PE+, and one or more of the combined adverse outcomes within the follow up period; (11) all enrolled subjects, and no combined adverse outcome through the follow up period; and (12) all enrolled subjects, and one or more of the combined adverse outcomes within the follow up period. Marker assays were performed as described above in Examples 3, 4, and 5. The results of the assays in each of these groups are presented in the following Table 7:
The following Tables 8 and 9 present results obtained by analyzing the univariate ability of each of the individual marker assays to distinguish between various subject groups. Data are presented as area under the ROC curve for the particular group comparison, calculated according to standard receiver operator characteristic methods, and include p values calculated using a Kruskal-Wallis Test assessing whether two groups are statistically different.
Preferred assays for diagnosis of VETD, PE, and/or DVT include those configured to detect one or more markers providing a ROC area of ≧0.6 in any of the three comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: acidic calponin, adrenomedullin, basic calponin, BMP-4, BNP, BNP1-108, BNP3-108, BNP79-108, CCL11, CK-BB, CRP, D-dimer, sELAF, endothelin-1, GSTP, IL-1ra, IL-25, leptin, sLTBR, MMP-9, MPO, NDKA, neuropilin-2, NGAL, PLGF-1, PLGF 1+2, activated protein C, total protein C, PSAP-A, PSAP-C, sRAGE, sPECAM-1, TIE-2, tissue factor, sTNFR1a, sTNFRSF7, TNFsR14, TpP, UFDP1H, UPA, VCAM-1, VE-cadherin, VEGF, sVEGF-R1, VWF-integrin, and ANP28-151. Particularly preferred assays include those detecting markers providing a ROC area of ≧0.7 in any of the three comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: acidic calponin, adrenomedullin, basic calponin, BNP1-108, BNP3-108, CK-BB, CRP, D-dimer, sELAF, GSTP, IL-1ra, MMP-9, MPO, NDKA, neuropilin-2, NGAL, PLGF 1+2, total protein C, sTNFR1a, UFDP1H, VCAM-1, VEGF, sVEGF-R1, and VWF-integrin. Most preferred assays include those detecting markers providing a ROC area of ≧0.8 in any of the three comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: adrenomedullin, basic calponin, D-dimer, GSTP, MMP-9, NDKA, neuropilin-2, total protein C, and UFDP1H.
Certain assays, such as those detecting one or more of BNP, BNP3-108, BNP79-108, CRP, D-dimer, sELAF, GSTP, sLTBR, MPO, NDKA, sTNFR1a, TNFsR14, UFDP1H, sVEGF-R1, and VWF-integrin, which provide a ROC area of ≧0.6 in distinguishing group (4) from group (2), are preferred for distinguishing amongst causes of VETD, and particularly in distinguishing PE from DVT.
Each of these preferred markers and marker assays may be used individually, or in panels comprising 1, 2, 3, 4, 5, 6, or more of these preferred markers, and optionally comprising one or more additional markers other than these preferred markers. The diagnostic methods of the present invention may be used to rule in or out VETD, PE, and/or DVT, most preferably when used together with other clinical signs and symptoms related to these conditions.
Preferred assays for prognosis of subjects diagnosed VETD, PE, and/or DVT include those configured to detect one or more markers providing a ROC area of ≧0.6 in any of the four comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: acidic calponin, adrenomedullin, angiopoietin-4, basic calponin, BMP-4, BNP, BNP1-108, BNP3-108, BNP79-108, CCL11, CGRP, CK-BB, CK-MB, CRP, D-dimer, sELAF, hFABP, IL-1ra, IL-25, leptin, sLTBR, MCP-1, MMP-9, MPO, NDKA, neuropilin-2, NGAL, PLGF 1+2, activated protein C, PSAP-A, PSAP-B, PSAP-C, PSAP-D, sRAGE, sPECAM-1, spectrin 145, tissue factor, sTNFR1a, sTNFRSF7, TNFsR14, TpP, UFDP1H, UPA, VCAM-1, VE-cadherin, VEGF, sVEGF-R1, VWF-integrin, and ANP28-151. Particularly preferred assays include those detecting markers providing a ROC area of ≧0.7 in any of the four comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: acidic calponin, adrenomedullin, basic calponin, BMP-4, BNP, BNP3-108, CK-BB, CK-MB, CRP, D-dimer, sELAF, IL-25, leptin, MCP-1, MMP-9, MPO, NDKA, NGAL, PSAP-A, PSAP-B, sPECAM-1, sTNFRSF7, UFDP1H, UPA, VEGF, and sVEGF-R1. Most preferred assays include those detecting markers providing a ROC area of ≧0.8 in any of the four comparison groups in the preceeding table. These include one or more assays detecting one or more of the following markers: CRP, D-dimer, IL-25, leptin, MCP-1, PSAP-A, PSAP-B, sPECAM-1, UFDP1H, and VEGF.
Each of these preferred markers and marker assays may be used individually, or in panels comprising 1, 2, 3, 4, 5, 6, or more of these preferred markers, and optionally comprising one or more additional markers other than these preferred markers. The prognostic methods of the present invention may be used to assign a prognosis to a subject diagnosed as having VETD, PE, and/or DVT, most preferably when used together with other clinical signs and symptoms related to these conditions.
In the following example, the “diseased” dataset represents a population of subjects from the subject collection described in Example 10 diagnosed as having PE and an adverse outcome as defined above in Example 3, while the “normal” dataset represents a population of subjects from the subject collection described in Example 10 for which a PE diagnosis was excluded. Samples were obtained for these subjects at presentation. Panels of markers were identified using the methods described in PCT application no. US03/41426, filed Dec. 23, 2003, and Example 8.
The following markers were selected for use in panels: D-Dimer (univariate ROC 0.887); MPO (univariate ROC 0.789); NDKA (univariate ROC 0.682); Neuropilin-2 (univariate ROC 0.52); CRP (univariate ROC 0.809); VEGF-r1 (univariate ROC 0.651); TNF-sR14 (univariate ROC 0.614); sPECAM-1 (univariate ROC 0.51); Tissue Factor (univariate ROC 0.588).
As is demonstrated by the foregoing data, panels comprising assays that detect the preferred markers from Example 11 can be used to assign a diagnosis of pulmonary embolism to a subject. Moreover, as measured by ROC area, such panels can provide diagnostic methods superior to methods in which a marker is used individually. Such panels preferably comprise 2, 3, 4, 5, 6, 7, 8, 9, 10, or more of the preferred marker assays from Example 11. Most preferred panels comprise an assay that detects D-dimer, and one or more assays that detect one or more of the preferred markers from Example 11 other than D-dimer.
One skilled in the art readily appreciates that the present invention is well adapted to carry out the objects and obtain the ends and advantages mentioned, as well as those inherent therein. The examples provided herein are representative of preferred embodiments, are exemplary, and are not intended as limitations on the scope of the invention.
It will be readily apparent to a person skilled in the art that varying substitutions and modifications may be made to the invention disclosed herein without departing from the scope and spirit of the invention.
All patents and publications mentioned in the specification are indicative of the levels of those of ordinary skill in the art to which the invention pertains. All patents and publications are herein incorporated by reference to the same extent as if each individual publication was specifically and individually indicated to be incorporated by reference.
The invention illustratively described herein suitably may be practiced in the absence of any element or elements, limitation or limitations which is not specifically disclosed herein. Thus, for example, in each instance herein any of the terms “comprising”, “consisting essentially of” and “consisting of” may be replaced with either of the other two terms. The terms and expressions which have been employed are used as terms of description and not of limitation, and there is no intention that in the use of such terms and expressions of excluding any equivalents of the features shown and described or portions thereof, but it is recognized that various modifications are possible within the scope of the invention claimed. Thus, it should be understood that although the present invention has been specifically disclosed by preferred embodiments and optional features, modification and variation of the concepts herein disclosed may be resorted to by those skilled in the art, and that such modifications and variations are considered to be within the scope of this invention as defined by the appended claims.
Other embodiments are set forth within the following claims.
Number | Date | Country | |
---|---|---|---|
60688942 | Jun 2005 | US |