The present invention relates to the prognosis of patients, especially patients suffering from hepatic disorders. More particularly, the present invention relates to a non-invasive prognostic method for assessing the risk of death or of liver-related event in a subject, wherein the subject is preferably affected with a hepatic disorder.
Epidemiological data highlight the significant burden of liver diseases worldwide, including in the United States and in Europe. The most common chronic liver diseases are caused by excessive alcohol consumption, or infections by the hepatitis C and B viruses, and an increasing proportion is related to obesity and metabolic syndrome (non-alcoholic fatty liver disease). Chronic liver diseases are important causes of morbidity and mortality.
In the prior art, interest was mainly focused on the diagnosis of these patients. The Applicant, for example, filed several patent applications on methods for assessing the presence and/or the severity of liver diseases. These methods are diagnostic methods, i.e. elaborated with reference to liver biopsy and to a diagnostic target which is fibrosis, significant fibrosis or cirrhosis. In diagnostic methods, the relevancy of the biomarkers relative to the stage of fibrosis are assessed with comparison to the biopsies.
Diagnostic methods give efficient information on the severity of a fibrosis. However, there is a further need for help in making decisions for clinical monitoring of patients and for early identification of patients at risk of liver-related complications or death.
This further need could be fulfilled by prognostic methods, i.e. methods where the relevancy of the biomarkers is assessed on the basis of the life expectancy of the patients of the cohort of reference or on the basis of complications experienced by the patients of the cohort of reference.
In the prior art, different scores are correlated with the prognosis of cirrhosis, such as, for example the MELD score and the Child-Pugh score. The MELD score (wherein MELD stands for Model for End-Stage Liver Disease) consists in the mathematical combination of serum bilirubin and creatinine levels, International Normalized Ratio (INR) for prothrombin time and etiology of liver disease (Kamath et al, Hepatology 2001; 33:464-470). The Child-Pugh score is based on the combination of total bilirubin, serum albumin, International Normalized Ratio (INR) or prothrombin time, ascites and hepatic encephalopathy.
However, both the MELD and Child-Pugh scores are directed to the prognosis of cirrhosis only (Metavir F4 stage) and cannot help for prognostic of patients having lesser fibrosis degree like Metavir F0, F1, F2 or F3 fibrosis stages.
Other recent studies use diagnostic tests for the prediction of the survival in patients having alcoholic liver disease. For example, Naveau et al (Hepatology 2009, 49(1): 97-105) describe the use of the biomarkers of Fibrotest (a diagnostic score), in multivariate analysis, for manufacturing a prognostic score named Fibrotest prognostic.
Moreover, Vergniol et al (Gastroenterology 2011; 140: 1970-1979) studied the efficiency of liver stiffness, FibroTest, APRI or FIB-4 for predicting survival or liver-related death in hepatic disease patients. Vergniol et al concluded that liver stiffness and FibroTest may both be used for predicting survival or liver-related death, and that their mathematical combination presents a significant prognostic accuracy. However, there is no statistical difference between the accuracies of FibroTest, liver stiffness and their combination. Furthermore, Vergniol et al does not describe or suggest a prognostic test for assessing the risk of liver-related event.
Applying diagnostic tests to prognosis may not outcome in the most accurate prognosis. This invention aims at proposing scores and methods intentionally designed for prognosis, on the basis of multivariate mathematical combination of blood biomarkers, clinical markers, tests, and/or physical data. The prognostic tests of the invention allow assessing the risk of death with an increased accuracy as compared to the prognostic tests of the prior art, and further allow assessing the risk of liver-related events.
The present invention relates to an in vitro non-invasive prognostic method for assessing the risk of death and/or of liver-related event in a subject, comprising:
In one embodiment, said at least one physical data is liver stiffness measurement (LSM), preferably measured by VCTE.
In one embodiment, said at least one variable obtained in step (a) is the group of variables to be used for obtaining a blood test result selected from the group comprising FibroMeter, CirrhoMeter, ELF score, FibroSpect, APRI, FIB-4 and Hepascore, and step (d) comprises combining said group of variables in a multivariate time-related model and/or said at least one blood test obtained in step (b) is selected from the group comprising FibroMeter, CirrhoMeter, ELF score, FibroSpect, APRI, FIB-4 and Hepascore.
In one embodiment, the method comprises mathematically combining in a multivariate time-related model in step (d) the following variables, or at least the following variables:
Another object of the invention is an in vitro non-invasive prognostic method for assessing the risk of liver-related events, especially complications, in a subject, comprising obtaining and mathematically combining in a multivariate time-related model
In one embodiment, the method comprises mathematically combining in a multivariate time related model the following variables, or at least the following variables:
In one embodiment, the method further comprises mathematically combining in a multivariate time-related model at least one physical data from medical imaging or clinical measurement, preferably from elastometry, more preferably from Vibration Controlled Transient Elastography (VCTE), even more preferably said at least one physical data is LSM measured by Vibration Controlled Transient Elastography (VCTE).
In one embodiment, said multivariate time related model is a COX model, and/or wherein said multivariate COX model is time-fixed or time-dependent, preferably is time-dependent.
The present invention also relates to an in vitro non-invasive prognostic method for assessing the risk of death and/or of liver-related event in a subject, comprising:
In one embodiment, said first blood test is FibroMeter and said second blood test result is CirrhoMeter.
In one embodiment, death is all-cause death and/or liver-related death.
In one embodiment, said liver-related event, especially complication, is selected from the list comprising ascites, encephalopathy, jaundice, occurrence of large esophageal varices, variceal bleeding, gastro-intestinal hemorrhage, hepato-renal syndrome, hepatocellular carcinoma, hepatic transplantation, oesophageal varices, and portal hypertension.
In one embodiment, the subject is affected with a liver disease, preferably selected from the list comprising significant porto-septal fibrosis, severe porto-septal fibrosis, centrolobular fibrosis, cirrhosis, persinusoidal fibrosis, the fibrosis being from alcoholic or non-alcoholic origin and/or the subject is a patient affected with a chronic disease, preferably said chronic disease is selected from the group comprising chronic viral hepatitis C, chronic viral hepatitis B, chronic viral hepatitis D, chronic viral hepatitis E, non-alcoholic fatty liver disease (NAFLD), alcoholic chronic liver disease, autoimmune hepatitis, primary biliary cirrhosis, hemochromatosis and Wilson disease.
In one embodiment, the method of the invention is computer implemented.
Another object of the invention is a microprocessor comprising a computer algorithm carrying out the method of the invention.
In the present invention, the following terms have the following meanings:
“Prognosis” refers to the statistical prediction of the occurrence of death or of the occurrence of a liver-related event (or of the first liver-related event) over a specific time period for a subject, such as, for example, within a period of at least 3 months, preferably 3 months, 6 months, 9 months, or 1, 2, 3, 4, 5, 6, 7, 8, 9, or 10 years as of date of assessment. The term “prognosis” thus relates to the assessment of a specific risk of death or of liver-related event for a subject.
“Liver-related event” refers to a liver-related event or complication requiring specific therapy or care or to the progression of the liver disease or disorder in a patient, including, without limitation, jaundice, ascites, encephalopathy, variceal bleeding, hepatocellular carcinoma, severe infection (like septic shock), hepato-renal syndrome, hepato-pulmonary syndrome.
“Risk” in the context of the present invention, relates to the probability that death or a liver-related event will occur over a specific time period, and can mean a subject's “absolute” risk or “relative” risk. Absolute risk can be measured with reference to either actual observation post-measurement for the relevant time cohort, or with reference to index values developed from statistically valid historical cohorts that have been followed for the relevant time period. Relative risk refers to the ratio of absolute risks of a subject compared either to the absolute risks of low risk cohorts or an average population risk, which can vary by how clinical risk factors are assessed. Odds ratios, the proportion of positive events to negative events for a given test result, are also commonly used (odds are according to the formula p/(1-p) wherein p is the probability of event and (1-p) is the probability of no event). In one embodiment, the method of the invention may further comprise a step of comparing the prognostic score obtained from the subject with a reference score.
“Score” refers to any digit value obtained by the mathematical combination (univariate or multivariate) of at least one biomarker and/or at least one clinical data and/or at least one physical data and/or at least one binary marker and/or at least one blood test result. In one embodiment, a score is an unbound digit value. In another embodiment, a score is a bound digit value, obtained by a mathematical function. Preferably, a score ranges from 0 to 1. In one embodiment, the at least one biomarker and/or at least one clinical data and/or at least one physical data and/or at least one binary marker and/or at least one score, mathematically combined in a score are independent, i.e. give each information that is different and not linked to the information given by the others.
“About” preceding a figure means plus or less 10% of the value of said figure.
“Biomarkers” refers to a variable that may be measured in a sample from the subject, wherein the sample may be a bodily fluid sample, such as, for example, a blood, serum or urine sample, preferably a blood or serum sample.
“Clinical data” refers to a data recovered from external observation of the subject, without the use of laboratory tests and the like.
“Binary marker” refers to a marker having the value 0 or 1 (or yes or no).
“Physical data” refers to a variable obtained by a physical method.
“Blood test” corresponds to a test comprising non-invasively measuring at least one data, and, when at least two data are measured, mathematically combining said at least two data within a score. In the present invention, said data may be a biomarker, a clinical data, a physical data, a binary marker or any combination thereof (such as, for example, any mathematical combination within a score).
The present invention relates to in vitro prognostic methods, preferably to non-invasive in vitro prognostic methods, for prognosing death (preferably liver-related death) or liver-related events, especially complications, in a subject. Accordingly, the methods of the invention may also be defined as in vitro methods for assessing the risk of death or of liver-related events, especially liver-related complications, in a subject. The present invention also relates to in vitro prognostic methods, preferably to non-invasive in vitro prognostic methods, for determining an increased risk of mortality or of liver-related event in a subject. The predictive power of the prognosis method of the invention is assessed by the Harrell's C index which reflects the discriminative ability.
In one embodiment, the prognosis method of the invention has a Harrell's C index of at least about 0.70, preferably of at least about 0.71, 0.72, 0.73, 0.74, 0.75, 0.76, 0.77, 0.78, 0.79, 0.80, 0.81, 0.82, 0.83, 0.84, 0.85, 0.86, 0.87, 0.88, 0.89, 0.90, 0.91, 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, 0.98, 0.99 or more.
The present invention thus relates to an in vitro prognostic method, preferably said method being non-invasive, for assessing the risk of death and/or of liver-related event, especially complication, in a subject, comprising:
In one embodiment, said blood test obtained in step (b) is not a APRI, FIB-4 or Actitest.
In one embodiment, at least 2, preferably 2, 3, 4, 5, 6, 7, 8 or more variables are obtained in step (a).
In one embodiment, the method of the invention further comprises a step of assessing the risk of death and/or of liver-related event, especially complication, in the subject according to the value of the prognostic score obtained in step (d).
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) said at least one variable obtained in step (a) and said at least one physical data obtained in step (c).
In another embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) said at least one blood test obtained in step (b) and said at least one physical data obtained in step (c).
In one embodiment, the terms “death” and “mortality” both refer to overall death or mortality (which may also be referred as all-cause death or mortality) and/or to liver-related death or mortality. Examples of causes of liver-related deaths include, but are not limited to, deaths consecutive to a portal hypertension related hemorrhage, deaths consecutive to an oesophageal or gastric varice related hemorrhage, a hepatocellular carcinoma, ascites, encephalopathy, liver failure with sepsis, acute on chronic liver failure, hepato-renal syndrome, hepato-pulmonary syndrome or other liver decompensation.
In one embodiment, the term “liver-related event, especially complications” refers to a liver-related event or complication requiring specific therapy or care, such as, for example, ascites, encephalopathy, jaundice (which may be defined as serum bilirubin >50 μmol/l ), occurrence of large esophageal varices (preferably having a diameter ≧5 mm, and/or preferably a diameter ≧25% of esophageal circumference), variceal bleeding, gastro-intestinal hemorrhage (such as, for example, due to portal hypertension), hepatorenal syndrome, hepato-pulmonary syndrome, hepatocellular carcinoma, hepatic transplantation, oesophageal varices, portal hypertension superior or equal to a predetermined threshold (such as, for example, hepatic vein pressure gradient superior or equal to 10 mm Hg or superior or equal to 12 mm Hg), severe infection (such as, for example, septic shock).
In one embodiment, the term “liver-related event, especially complication” refers to the progression of the liver disease or disorder in a patient, such as, for example, the appearance of cirrhosis in a fibrotic non-cirrhotic patient, or the fact, for a patient, to cross a predetermined threshold (such as, for example, FibroMeter result superior or equal to 0.982, or Fibroscan result superior or equal to 14 kPa).
According to the invention, death (including all-cause death and liver-related death) does not refer to a liver-related event, especially complication.
In one embodiment, the method of the invention is for prognosing the first liver-related event, especially complication, in a patient.
In one embodiment of the invention, the prognostic method of the invention is for assessing the risk of death or of a liver-related event within a period of at least 3 months, preferably 3 months, 6 months, 9 months, 1, 2, 3, 4, 5, 6, 7, 8, 9, or 10 years as of date of assessment.
In one embodiment, the subject is a mammal, preferably a human. In one embodiment, the subject is a male or a female. In one embodiment, the subject is an adult or a child.
In one embodiment, the subject is affected, preferably is diagnosed with a liver disease or disorder, and the method of the invention may thus be defined as a method of prognosing liver disease or disorder.
Preferably, the subject is a patient, i.e. the subject is awaiting the receipt of, or is receiving medical care or is/will be the object of a medical procedure.
In one embodiment, the subject is affected with a liver disease or disorder, preferably selected from the list comprising significant porto-septal fibrosis, severe porto-septal fibrosis, centrolobular fibrosis, cirrhosis, persinusoidal fibrosis, the fibrosis being from alcoholic or non-alcoholic origin.
In one embodiment, the subject is a patient affected with a chronic disease, preferably said chronic disease is selected from the group comprising chronic viral hepatitis C, chronic viral hepatitis B, chronic viral hepatitis D, chronic viral hepatitis E, non-alcoholic fatty liver disease (NAFLD), alcoholic chronic liver disease, autoimmune hepatitis, primary biliary cirrhosis, hemochromatosis and Wilson disease.
Examples of biomarkers include, but are not limited to, glycemia, total cholesterol, HDL cholesterol (HDL), LDL cholesterol (LDL), AST (aspartate aminotransferase), ALT (alanine aminotransferase), AST/ALT (ratio), AST.ALT (product), ferritin, platelets (PLT), AST/PLT (ratio), prothrombin time (PT) or prothrombin index (PI), hyaluronic acid (HA or hyaluronate), haemoglobin, triglycerides, alpha-2 macroglobulin (A2M), gamma-glutamyl transpeptidase (GGT), urea, bilirubin, apolipoprotein A1 (ApoA1), type III procollagen N-terminal propeptide (P3NP or P3P), gamma-globulins (GBL), sodium (Na), albumin (ALB), ferritine (Fer), glucose (Glu), alkaline phosphatases (ALP), YKL-40 (human cartilage glycoprotein 39), tissue inhibitor of matrix metalloproteinase 1 (TIMP-1), TGF, cytokeratine 18 and matrix metalloproteinase 2 (MMP-2) to 9 (MMP-9), haptoglobin, ratios and mathematical combinations thereof.
Examples of clinical data include, but are not limited to, weight, body mass index, age, sex, hip perimeter, abdominal perimeter or height and mathematical combinations thereof, such as, for example, the ratio thereof, such as for example hip perimeter/abdominal perimeter.
Examples of binary non-invasive test markers include, but are not limited to, diabetes, SVR (wherein SVR stands for sustained virologic response, and is defined as aviremia 6 weeks, preferably 12 weeks, more preferably 24 weeks after completion of antiviral therapy for chronic hepatitis C virus (HCV) infection), etiology, and NAFLD. Regarding the binary marker “etiology”, the skilled artisan knows that said variable is a single or multiple binary marker, and that for liver disorders, etiology may be NAFLD, alcohol, virus or other. Thus, the binary marker might be expressed as NAFLD vs others (single binary marker) or as NAFLD vs reference etiology plus virus vs reference etiology and so on (multiple binary marker).
In one embodiment, the blood test of the invention corresponds to a blood test selected from the group comprising ELF, FibroSpect™, APRI, FIB-4, Hepascore, Fibrotest™, FibroMeter™, CirrhoMeter™, CombiMeter, Elasto-FibroMeter™, InflaMeter™, Elasto-Fibrotest. As these blood tests are diagnostic tests, they are based on multivariate mathematical combination, such as, for example, binary logistic regression.
ELF is a blood test based on hyaluronic acid, P3P, TIMP-1 and age.
FibroSpect™ is a blood test based on hyaluronic acid, TIMP-1 and A2M.
APRI is a blood test based on platelet and AST.
FIB-4 is a blood test based on platelet, AST, ALT and age.
HEPASCORE is a blood test based on hyaluronic acid, bilirubin, alpha2-macroglobulin, GGT, age and sex.
FIBROTEST™ is a blood test based on alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, age and sex.
ELASTO-FIBROTEST is a test based on the mathematical combination of variables of FIBROTEST or of the result of a FIBROTEST, with LSM measurement, measured for example by Fibroscan™.
FIBROMETER™ and CIRRHOMETER™ together form a family of blood tests, the content of which depends on the cause of chronic liver disease and the diagnostic target (such as, for example, fibrosis, significant fibrosis (moderate bridging fibrosis), severe fibrosis (extensive bridging fibrosis) or cirrhosis). This blood test family is called FM family and is detailed in the table below.
aHA may be replaced by GGT
COMBIMETER™ or Elasto-FibroMeter™ is a family of tests based on the mathematical combination of variables of the FM family (as detailed in the Table hereinabove) or of the result of a test of the FM family with FIBROSCAN™ result. In one embodiment, said mathematical combination is a binary logistic regression.
In one embodiment, CombiMeter™ or Elasto-FibroMeter™ results in a score based on the mathematical combination of physical data from liver or spleen elastometry such as dispersion index from Fibroscan™ such as IQR or IQR/median or median, preferably of Fibroscan™ median with at least 3, preferably at least 4, 5, 6, 7 or more and more preferably of 7 or 8 or 9 biomarkers and/or clinical data selected from the list comprising glycemia, total cholesterol, HDL cholesterol (HDL), LDL cholesterol (LDL), AST (aspartate aminotransferase), ALT (alanine aminotransferase), AST/ALT, AST.ALT, ferritin, platelets (PLT), AST/PLT, prothrombin time (PT) or prothrombin index (PI), hyaluronic acid (HA or hyaluronate), haemoglobin, triglycerides, alpha-2 macroglobulin (A2M), gamma-glutamyl transpeptidase (GGT), urea, bilirubin, apolipoprotein A1 (ApoA1), type III procollagen N-terminal propeptide (P3NP), gamma-globulins (GBL), sodium (Na), albumin (ALB), ferritine (Fer), glucose (Glu), alkaline phosphatases (ALP), YKL-40 (human cartilage glycoprotein 39), tissue inhibitor of matrix metalloproteinase 1 (TIMP-1), TGF, cytokeratine 18 and matrix metalloproteinase 2 (MMP-2) to 9 (MMP-9), haptoglobin, diabetes, weight, body mass index, age, sex, hip perimeter, abdominal perimeter or height and ratios and mathematical combinations thereof.
INFLAMETER™ is a companion test reflecting necro-inflammatory activity including ALT, A2M, PI, and platelets.
In one embodiment, the method of the invention does not comprise obtaining a Fibrostest™ in step (a).
In one embodiment, the method of the invention does not comprise obtaining an Actitest™ in step (a) and/or in step (b). ACTITEST™ is a blood test based on alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, ALT, age and sex.
In one embodiment, the method of the invention does not comprise obtaining an APRI in step (a) and/or in step (b).
In one embodiment, the method of the invention does not comprise obtaining a FIB-4 in step (a) and/or in step (b).
In one embodiment, the method of the invention does not comprise or consists in mathematically combining in a multivariate domain a Fibrostest™ result and a physical data, such as, for example, liver stiffness measurement (LSM) measured by VCTE (such as, for example, Fibroscan).
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of a blood test as hereinabove defined.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of ELF, i.e. hyaluronic acid, P3P, TIMP-1 and age.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FibroSpect™, i.e. hyaluronic acid, TIMP-1 and A2M.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of APRI, i.e. platelet and AST. In one embodiment, the method of the invention does not comprise or consist in obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of APRI, i.e. platelet and AST.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FIB-4, i.e. platelet, ASAT, ALT and age. In one embodiment, the method of the invention does not comprise or consist in obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FIB-4, i.e. platelet, ASAT, ALT and age.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of HEPASCORE, i.e. hyaluronic acid, bilirubin, alpha2-macroglobulin, GGT, age and sex.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FIBROTEST™, i.e. alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, age and sex. In one embodiment, the method of the invention does not comprise or consist in obtaining in step (a), and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FIBROTEST™, i.e. alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, age and sex.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of ACTITEST™, i.e. alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, ALT, age and sex. In one embodiment, the method of the invention does not comprise or consist in obtaining in step (a), and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of ACTITEST™, i.e. alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, ALT, age and sex.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of FIBROMETER™ and/or CIRRHOMETER™ as listed in the Table hereinabove.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of COMBIMETER™ or Elasto-FibroMeter™ as presented hereinabove.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of Elasto-Fibrotest as presented hereinabove.
In one embodiment, the method of the invention comprises obtaining in step (a) and mathematically combining in a multivariate time-related model with a physical data in step (d) the variables of INFLAMETER™, i.e. ALT, A2M, PI, and platelets.
Examples of physical methods include, but are not limited to, medical imaging data and clinical measurements, such as, for example, measurement of spleen, especially spleen length. According to an embodiment, the physical method is selected from the group comprising ultrasonography, especially Doppler-ultrasonography and elastometry ultrasonography and velocimetry ultrasonography (preferred tests using said data are vibration controlled transient elastography (VCTE, also known as Fibroscan™), ARFI, VTE, supersonic elastometry (supersonic imaging), MRI (Magnetic Resonance Imaging), and MNR (Magnetic Nuclear Resonance) as used in spectroscopy, especially MNR elastometry or velocimetry.
In one embodiment, the physical method is VCTE, ARFI, VTE, supersonic elastometry or MRI stiffness.
Preferably, the data is an elastometry data, preferably Liver Stiffness Evaluation (LSE) data or spleen stiffness evaluation, which may be for example obtained by VCTE. According to a preferred embodiment of the invention, the physical data is liver stiffness measurement (LSM), preferably measured by VCTE.
In a particular embodiment, the physical data is Liver stiffness measurement (LSM) by VCTE (also known as Fibroscan™, Paris, France), preferably performed with the M probe. Preferably, examination conditions are those recommended by the manufacturer, with the objective of obtaining at least 3 and preferably 10 valid measurements. Results may be expressed as the median (kilopascals) of all valid measurements, or as IQR or as the ratio (IQR/median).
In one embodiment of the invention, the method of the invention comprises carrying out a VCTE, which refers to obtaining at least 3 and preferably 10 valid measurements and recovering a physical data corresponding to the median in kilopascals of all valid measurements.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) prothrombin index, bilirubin, urea, hyaluronic acid, LSM, age and sex.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, etiology, LSM, AST and urea.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, NAFLD, LSM, AST and urea.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, LSM, AST and urea.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, NAFLD, prothrombin time, bilirubin, urea and LSM.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, NAFLD, prothrombin time, AST, urea and LSM.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, prothrombin time, bilirubin, urea, hyaluronic acid and LSM.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, prothrombin time, bilirubin, urea and LSM.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, prothrombin time, AST, urea and LSM.
In one embodiment, the method of the invention comprises mathematically combining in a multivariate time-related model in step (d) age, sex, etiology, AST, urea and LSM.
Another object of the invention is an in vitro non-invasive prognostic method for assessing the risk of liver-related events, especially complications, in a subject, comprising obtaining and mathematically combining in a multivariate time-related model
thereby obtaining a prognostic score,
wherein biomarkers, clinical data, binary markers and blood tests are as described hereinabove.
In one embodiment, the method of the invention further comprises a step of assessing the risk of liver-related events, especially complications, in the subject according to the value of the prognostic score.
In one embodiment, the method of the invention comprises or consists in obtaining and mathematically combining in a multivariate time-related model the variables of a blood test as hereinabove defined.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of ELF, i.e. hyaluronic acid, P3P, TIMP-1 and age.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of FibroSpect™, i.e. hyaluronic acid, TIMP-1 and A2M.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of APRI, i.e. platelet and AST.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of FIB-4, i.e. platelet, ASAT, ALT and age.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of HEPASCORE, i.e. hyaluronic acid, bilirubin, a1pha2-macroglobulin, GGT, age and sex.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of FIBROTEST™, i.e. alpha2-macroglobulin, haptoglobin, apolipoprotein A1, total bilirubin, GGT, age and sex.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of FIBROMETER™ and/or CIRRHOMETER™ as listed in the Table hereinabove.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of COMBIMETER™ or Elasto-FibroMeter™ as presented hereinabove.
In one embodiment, the method of the invention comprises obtaining and mathematically combining in a multivariate time-related model the variables of INFLAMETER™, i.e. ALT, A2M, PI, and platelets.
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid, AST, SVR, platelets and GGT.
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid, platelets and SVR.
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid and SVR.
In one embodiment of the invention, the in vitro non-invasive prognostic method for assessing the risk of liver-related events, especially complications, in a subject further comprises:
wherein said physical data is as defined hereinabove.
In one embodiment, the method of the invention further comprises a step of assessing the risk of liver-related events, especially complications, in the subject according to the value of the prognostic score.
In one embodiment, said physical data is a data obtained by elastometry, more preferably by Vibration Controlled Transient Elastography (VCTE), even more preferably said at least one physical data is LSM measured by Vibration Controlled Transient Elastography (VCTE).
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid, AST, SVR, platelets, GGT and LSM.
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid, platelets, SVR and LSM.
In one embodiment of the invention, the method comprises obtaining and mathematically combining in a multivariate time-related model hyaluronic acid, SVR and LSM.
In one embodiment, said multivariate time-related model is a COX proportional hazard regression model.
In one embodiment, the prognostic score of the invention is based on time-fixed Cox model, i.e. the method of the invention comprises only one measurement of variables selected from biomarkers, clinical data, binary markers and blood tests and/or physical data, which are generally recorded at baseline time.
In another embodiment, the prognostic score of the invention is based on time-dependent Cox model, i.e. the method of the invention comprises at least two successive measurements of variables selected from biomarkers, clinical data, binary markers and blood tests and/or physical data recorded at various times including baseline time.
In one embodiment, the prognosis method of the invention further comprises a step of comparing the prognostic score obtained for the subject with a reference score.
According to the invention, the prognostic score is indicative of the risk of death or of liver-related event in the subject.
In one embodiment, the prognostic method of the invention further comprises comparing the prognostic score obtained for the subject with a reference prognostic score. In one embodiment, a difference or an absence of difference between the prognostic score obtained for the subject and the reference score is indicative of an increased risk of death or of liver-related events.
In one embodiment, the reference score may be an index score or may be derived from one or more risk prediction algorithms or computed indices for the prognosis of death or liver-related events. A reference score can be relative to a number or value derived from population studies, including without limitation, such subjects having similar fibrosis level, subjects of the same or similar age range, subjects in the same or similar ethnic group.
In one embodiment of the present invention, the reference score is derived from one or more subjects (i.e. the reference population) who are substantially healthy, i.e. who are not affected with a liver disease or disorder. In another embodiment of the present invention, the reference score is derived from one or more subjects (i.e. the reference population) who are affected with a liver disease or disorder.
In one embodiment of the invention, a prognostic score obtained for the subject which is different from the reference score is indicative of a risk of death or of liver-related event. Preferably, in this embodiment, the reference score is derived from substantially healthy subjects.
In another embodiment of the invention, a prognostic score obtained for the subject which is not different from the reference score is indicative of a risk of death or of liver-related event. Preferably, in this embodiment, the reference score is derived from liver disease patients, more preferably from liver disease patients with a poor prognosis.
In one embodiment, a first value is considered different from a second value if the first one is superior or inferior by about 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100% or more compared to the second one.
In one embodiment, the prognostic scores measured in the reference population used for obtaining the reference score are classified in percentiles, wherein the prognostic scores obtained by all subjects of the reference population are ranged according to their numerical value in ascending order. In one embodiment of the invention, the percentiles are percentiles of subjects, i.e. each percentile comprises the same number of subjects. Therefore, the first percentile corresponds to subjects with the lowest prognostic scores, while the last percentile corresponds to subjects with the highest prognostic scores.
In one embodiment, the reference score corresponds to the highest prognostic score of the first percentile of the reference population.
In another embodiment, the reference score corresponds to the highest prognostic score of the second, third . . . or penultimate percentile of the reference population.
In one embodiment, when three percentiles are drawn, each percentile is named a tertile. According to this embodiment, the reference score corresponds to the highest prognostic score of the first or of the second tertile.
In another embodiment, when four percentiles are drawn, each percentile is named a quartile. According to this embodiment, the reference value corresponds to the highest prognostic score of the first, second or third quartile.
In another embodiment, when five percentiles are drawn, each percentile is named a quintile. According to this embodiment, the reference value corresponds to the highest prognostic score of the first, second, third or fourth quintile.
In one embodiment of the invention, the prognostic method of the invention thus leads to the classification of the subject in a class of a classification, wherein each class of the classification corresponds to a specific risk of death or liver-related event. Such classifications, wherein a specific prognostic score corresponds to a specific assessed risk, are also objects of the present invention.
For example, the prognostic method of the invention leads to the classification of the subject in a classification comprising 5 classes, each corresponding to a quintile: the first class (or quintile) corresponds to a risk of death or of liver-related event of 0% in 5 years, the second class (or quintile) corresponds to a risk of death or of liver-related event superior to 0% and inferior or equal to 25%, the third class (or quintile) corresponds to a risk of death or of liver-related event superior to 25% and inferior or equal to 50%, the fourth class (or quintile) corresponds to a risk of death or of liver-related event superior to 50% and inferior or equal to 75%, and the fifth and last class (or quintile) corresponds to a risk of death or of liver-related event superior to 75%.
The present invention also relates to an in vitro non-invasive prognostic method for assessing the risk of death and/or of liver-related event in a subject, comprising:
In one embodiment of the invention, the first blood test is FibroMeter and the second blood test result is CirrhoMeter.
In one embodiment of the invention, the prognostic method of the invention thus leads to the classification of the subject in a class of a combined classification (which may also be referred to as a score risk chart), wherein said combined classification is in the form of a table wherein each cell is associated with a specific assessed risk of death or liver-related event. Such score risk charts are also objects of the present invention.
In one embodiment, the score risk chart is a matrix table with 2 dimensions, wherein each dimension corresponds to a blood test and/or a physical method, for which a classification was previously constructed with the use of a reference to a reference population. For example, carrying out a Fibroscan may result in the classification of a subject in one of 3 fibrosis stages (F0-F1, F2-F3 or F4). Therefore, the first dimension of the score risk chart (lines) corresponds to stages for the first blood test or physical method, while the second dimension (columns) corresponds to stages for the second blood test or physical method.
In one embodiment, the score risk chart comprises additional dimensions, corresponding for example to clinical data, such as, for example, age and sex.
An example of 4-dimensions score risk chart of the invention is shown in
In one embodiment of the invention, the method of the invention is computer implemented.
The present invention thus also relates to a microprocessor comprising a computer algorithm carrying out the prognostic method of the invention.
The skilled artisan would easily deduce that the prognostic score of the invention being indicative of a particular risk of death or of liver-related event, it may be used by the physician willing to provide the best medical care to his/her patient. For example, a patient presenting a high risk of death or liver-related event will require more frequent follow-up than a patient with low risk.
Therefore, the present invention also relates to a method for adapting the treatment, the medical care or the follow-up of a patient, wherein said method comprises measuring the prognostic score of the invention for said patient.
The present invention also relates to a method for monitoring the treatment of a patient, wherein said method comprises measuring the prognostic score of the invention for said patient.
The present invention also relates to a method for treating a hepatic disease patient, preferably a patient with chronic liver disease, wherein said method comprises (i) measuring the prognostic score of the invention for said patient and (ii) treating the patient according to the prognostic score measured in step (i).
Moreover, the prognostic score of the invention may be used for prioritizing patients for hepatic transplantation.
The present invention is further illustrated by the following examples.
Patients and Methods
Study Population
We included all subjects over 18 years of age who consulted or were hospitalized for a chronic liver disease in the Department of Hepatology of the University Hospital of Angers from January 2005 to December 2009, and who agreed to participate in the cohort SNIFF (Score Non Invasif de Fibrose du Foie). Patients have been included in the cohort whatever severity and etiology of chronic liver disease (HVC, HVB, alcohol abuse, non-alcoholic fatty liver, other causes). They were followed until death or up to Jan. 1, 2011.
Non-Invasive Evaluation of Liver Fibrosis
All patients had a non-invasive evaluation of liver fibrosis at the time of inclusion in the cohort and at every visit during follow-up, either by a blood test or a liver stiffness measurement, or both.
Blood tests for liver fibrosis—Blood sample analyses were centralized in the laboratory of our hospital. Biological parameters measured (alpha-2-macroglobulin, platelet count, aspartate aminotransferase, alanine aminotransferase, gamma-glutamyl transpeptidase, bilirubin, hyaluronic acid, prothrombin index, urea) allowed the computation of at least one of the following blood tests for liver fibrosis: APRI (Wai et al, Hepatology 2003; 38:518-26), FIB-4 (Sterling et al, Hepatology 2006; 43:1317-25), Hepascore (Adams et al, Clin Chem 2005; 51:1867-73), FibroMeterV2G (Leroy et al, Clin Biochem 2008; 41:1368-76). The MELD score which has been specifically designed to predict short-term mortality in patients with end-stage of liver disease (Kamath et al, Hepatology 2001; 33(2):464-70) was also calculated in cirrhotic patients.
Liver stiffness measurement—Liver stiffness measurement (LSM) by Fibroscan (Echosens, Paris, France) was performed with the M probe. Examination conditions were those recommended by the manufacturer, with the objective of obtaining at least 10 valid measurements. Results were expressed as the median (kilopascals) of all valid measurements, and classified as “very reliable”, “reliable” and “poorly reliable” by using ratio of interquartile range to median value (IQR/M) according to liver stiffness level: “very reliable” (IQR/M≦0.10), “reliable” (0.10<IQR/M≦0.30, or IQR/M>0.30 with M<7.1 kPa), “poorly reliable” (IQR/M>0.30 with M>7.1 kPa).
Endpoint Assessment
The endpoints of interest for these analyses were all-cause mortality and liver-related mortality. The date and primary cause of death were recorded during follow-up and checked in May 2011 using the computerized national death registry. Data for those patients who could not be found by this way were obtained from hospital database, and confirmed by information from general practitioner. Causes of death were coded according to the 10th revision of the International Classification of Diseases (ICD-10) and pooled into categories.
Statistical Methods
Statistical analyses were performed using SAS version 9.2 (SAS Institute Inc., Cary, N.C., USA) and R version 3.0.2 (Vienna, Austria, available from: http://www.R-project.org/). All tests were two-sided, and p<0.05 was considered statistically significant.
Data are expressed as the median (interquartile range), or frequency (%), as appropriate. Group comparisons were performed using Fisher's exact test for categorical data and nonparametric Mann-Whitney U test for continuous variables. Kaplan-Meier analyses with log-rank comparisons were used to calculate the median survival times and the survival differences between quartiles of each liver fibrosis diagnostic test at baseline.
Comparison of prognostic accuracy of non-invasive fibrosis tests—The predictive performance of non-invasive liver fibrosis tests was firstly evaluated on the entire follow-up by the Harrell's C-index. This non-parametric measure corresponds to an extension of the binary outcome AUC for time-to-event outcomes (censored survival data), and evaluates the amount of agreement between predictions (that are here the baseline values of fibrosis test) and the observed survival time comparing not only subjects with events to those without events, but also subjects having events that happened at different points in time. The results of C-index vary between 0.5 (no predictive power) and 1 (perfect concordance). The 95% confidence interval of each C-index was calculated by a bootstrap method, using 1,000 random samples (with replacement) of the same size as the original data set. Paired comparisons of C-indexes of fibrosis tests were also performed using a bootstrap resampling procedure. Briefly, the difference between two C-indexes was calculated on each of the 1,000 bootstrap samples to obtain the bootstrap 95% confidence interval of the difference. Then, the estimate of the standard error from this bootstrap distribution was used to calculate the standardized difference between the C-indexes and find the corresponding p-value.
Thereafter, we examined the ability of non-invasive fibrosis tests to predict the risk of death at various time horizons (different delays of prediction) from 0 to 5 years. The reported AUC(t) were computed and compared using the timeROC R package, they are based on the cumulative/dynamic ROC curve definition at a given time t and are corrected for right-censoring. Time-dependent ROC and corresponding AUC can be defined in different ways, the cumulative cases/dynamic controls definition is usually considered the most relevant as clinicians often want to discriminate subjects who have the event before time t (and not at a specific time t, as in “incident cases” definition) and subjects free of the event at time t (and not at a later fixed pre-specified time, e.g. the end of the study, as in “static controls” definition). The interpretation of cumulative/dynamic AUC(t) is similar to that of a standard AUC: the higher the AUC(t) is, the better the performance of the test to discriminate between patients who will do the event before time t and other patients.
Multivariate prognostic models for all-cause and liver-related mortality—We randomly split the sample of 1,565 patients having all tests into a training set (n=783) and a validation set (n=782). We built prognostic models with the training set. Age, sex, etiology of chronic liver disease and all non-invasive liver fibrosis tests were evaluated as potential predictors of death in both time-fixed and time-dependent multivariate Cox models. In the time-fixed model, only the baseline values of fibrosis tests were taken into account. The time-dependent model additionally used the repeated measurements of the liver fibrosis tests until 1 year before the end of follow-up. We applied the method of Murtaugh, et al. (Hepatology 1994; 20:126-34) to update the covariate over time in the model. We have assumed that the values observed on a patient at a given time remain unchanged until the next measurement because this corresponds to the clinical situation. By dividing each patient's follow-up in successive intervals and taking into account the measurements available at the start of each, 1337 survival data points were generated for the training set and 1342 for the validation set. For each multivariate model, assessment of proportional hazards was performed. We also tested the one-way interactions between the independent predictors, and checked that continuous covariates fit the appropriate functional form in the Cox model.
The accuracy of time-fixed and time-dependent prognostic models was evaluated in the validation set both in terms of calibration and discrimination performance. The calibration was assessed by comparing observed and estimated survival rates of the patients divided into 5 groups (quintiles) according to their estimated risk of death. Observed survival rates for each quintile group were calculated according to the Kaplan-Meier method. Estimated survival rates were computed as described elsewhere (Ripoll et al, Hepatology 2005; 42:793-801): mean values of each variable (most frequent modality for categorical variable) in the multivariate Cox model were calculated in each quintile subgroup and the predicted survival rate of the subgroup was defined as the death risk score calculated for these mean values (Σbimi, where bi represents the estimated regression coefficient of each variable in model and mi represents the mean value of the variable in the subgroup). Discrimination ability of prognostic models was assessed by the Harrell's C-index.
Results
Characteristics of Patients
A total of 3,352 patients have been enrolled in the cohort at the end of December 2009, their main baseline characteristics are summarized in Table 1. More than two thirds of subjects were males, the average age at inclusion was 54.9±14.3. The chronic liver disease was due to excessive alcohol consumption in 42.7% of patients, whereas 30.4% had viral chronic hepatitis, 19.0% had non-alcoholic fatty liver disease (NAFLD) and the remaining 7.9% had other etiologies. A liver stiffness measurement was performed at baseline in 3,079 patients (91.9%), median baseline LSM result was 8.5 [1st and 3rd quartiles: 5.1-9.5]. Each liver fibrosis blood test was available in more than 2,500 patients (>75%) of the cohort.
(a)p-value from Fisher's exact test (categorical variables) or Mann-Whitney-Wilcoxon test (quantitative variables).
On 1st January 2011, the median duration of follow-up in the whole cohort was 3.3 years [interquartile range (IQR) 1.9-4.6]. During follow-up period, 588 patients (17.5%) died, 258 deaths (7.7%) were liver-related. Kaplan-Meier estimates of overall survival (without all-cause death) at 1, 3 and 5 years were respectively 91.6%, 83.0% and 75.5%. The liver-related survival (without liver-related death) at 1, 3 and 5 years were 95.7%, 92.5% and 86.6%.
According to data from the French national death registry, the other causes of death were mainly neoplasms unrelated to liver and intrahepatic bile ducts (ICD-10 COO-D48 except C22, n=104), diseases of the circulatory system (ICD-10 100-199, n=79), external causes of morbidity and mortality (ICD-10 V01-Y98, n=36), endocrine nutritional and metabolic diseases (ICD-10 E00-E90, n=17), diseases of the digestive system other than of the liver (ICD-10 K00-K93 except K70-K77, n=16). Cause of death was unknown in 12 patients.
Univariate Relationships with Mortality
Patients who died had quite different baseline characteristics, as compared to those alive at the end of follow-up (Table 2): they were older and more frequently males, they were twice to have an alcoholic liver disease and had more severe disease at baseline, as reflected of all non-invasive liver fibrosis tests at baseline. Kaplan-Meier estimates of both overall survival and liver-related survival significantly differed according to quartiles of each liver fibrosis test values at baseline (all p-values from log-rank comparisons <0.0001).
(a)All-cause death.
(b)P-value from Fisher's exact test (categorical variables) or Mann-Whitney-Wilcoxon test (quantitative variables).
(c)As classified by recently proposed reliability criteria: “very reliable” (IQR/median ≦0.10) and “reliable” (0.10 < IQR/median ≦ 0.30, or IQR/median >0.30 with LSE median <7.1 kPa).
Comparisons of the Ability of Liver Fibrosis Tests to Predict Mortality
Statistical comparisons of the predictive power of non-invasive liver fibrosis tests were performed in the subsample of patients having available data for all tests within a 3-month interval at inclusion. There was no difference in socio-demographic characteristics and causes of chronic liver disease between these 1,565 patients as compared to the other patients from the cohort (Table 1). But patients having all tests were more recently included in the cohort, and consequently had a slightly lower follow-up duration: median follow-up was 2.8 years [1st and 3rd quartiles: 1.7-3.9]. They had also a higher baseline liver stiffness median value than non-selected patients, but conversely they had lower baseline values for most liver fibrosis blood tests.
The global predictive power of APRI, as assessed by the Harrell's C-index, was much lower than other non-invasive liver fibrosis tests for both all-cause mortality and liver-related mortality (Table 3).
(a)significantly lower than LSM (p < 0.0001), FIB-4 (p < 0.0001), Hepascore (p < 0.0001), FibroMeterV2G (p < 0.0001).
(b)significantly higher than those of LSM (p < 0.040), APRI (p < 0.0001), FIB-4 (p < 0.0001), Hepascore (p = 0.0086).
(c)significantly lower than LSM (p < 0.0001), FIB-4 (p < 0.0001), Hepascore (p < 0.0001), FibroMeterV2G (p < 0.0001).
(d)significantly lower than FibroMeterV2G (p < 0.0001) and LSM very reliable/reliable results (p < 0.030).
FibroMeterV2G had the greatest performance in predicting all-cause mortality, significant difference was observed as compared to all other tests. The global prognostic performance of LSM was only marginally altered when including LSM results considered as poorly reliable for diagnostic purpose. For liver-related mortality, both FibroMeterV2G and LSM had higher prognostic value, with significant differences in comparison with APRI and FIB-4. According to time-dependent AUC, liver fibrosis tests had similar trend of predictive accuracy over time: as expected, the predictive power of tests was the greatest in the few months after the measure, then decreased to a performance level that remains relatively stable until 5 years (
We further compared the AUC(t) of liver fibrosis tests to MELD score in the subsample of 448 patients having cirrhosis, as defined by LSM≧14 KPa. Only 4 cirrhotic patients had not available data for MELD at baseline, results of the time-dependent ROC curve analysis in the remaining 444 patients are graphically presented in
Comparative analysis between liver fibrosis tests were also performed in the 3 main causes of chronic liver disease separately (alcoholic liver disease (ALD), chronic viral hepatitis, NAFLD). Results of time-dependent AUC are reported in
Multivariate Prognostic Models Based only on Baseline Liver Fibrosis Tests
Multivariate analyses were performed after randomly splitting the subsample of 1,565 patients having all tests at baseline into a training set and a validation set. The main characteristics of these two subgroups of patients are presented in Table 4.
By multivariate analysis in the training set, we observed that each blood fibrosis test and LSM were independent prognostic factors (data not shown). Age and sex were also significant in all multivariate Cox models, but not the cause of chronic liver disease, nor the interquartile range of LSM results. According to the Akaike information criteria (AIC), the best multivariate Cox model was the one that included age, sex, LSM and FibroMeterV2G, either to predict all-cause mortality or mortality due liver disease. No interactions between the independent predictors were found. Details of the coefficients of the multivariate models minimizing the AIC are given in Table 5.
The agreement between the observed and predicted survival (calibration) was less satisfying (
Multivariate Prognostic Models with Liver Fibrosis Tests Defined as Time-Varying Covariates
The non-invasive tests were carried out from 1 to 8 times from inclusion in the cohort to one year before the end of follow-up. The median interval length between measurements was 1.1 years [1st and 3rd quartiles: 0.6-1.7]. After applying the time-dependent prognostic model to the validation set, we observed a better match between predicted and observed liver-related survival (
We then compared overall prognostic performance though discriminative ability of liver-fibrosis diagnostic tests and of the prognostic model of the present invention, based on age, sex, LSM and prothrombin index, bilirubin, urea, hyaluronic acid. Prognostic performance was assessed by Harrell's index [95% CI] in the 1565 patients having all non-invasive liver fibrosis tests at inclusion in the cohort. Results are shown in Table 7.
(a)significantly lower than those of FibroMeterV2G (p <0.007), EFMV2G (p <0.002), EFMV3G (p <0.008)
(b)significantly lower than those of FibroMeterV2G (p <0.040)
(c)significantly higher than all liver fibrosis diagnostic tests (p <0.001)
(d)significantly lower than those of EFMV2G (p <0.030)
(e)significantly higher than CirrhoMeterV2G (p <0.012), FibroMeterV3G (p <0.001), CirrhoMeterV3G (p <0.001)
These results provide new insight about the prognostic value of liver fibrosis diagnostic tests, comparing their ability to identify patients with higher risk of death, at various time points, during the 5 years following the initial evaluation of liver fibrosis. Furthermore, our results indicate that combining FibroMeterV2G and LSM increased prognostic performance, and more particularly in patients with ALD where the prognostic value of LSM was lower than in other causes of chronic liver diseases.
The first aim of this study was to compare the prognostic accuracy of several blood fibrosis tests and liver biopsy for the prediction of liver-related events in a large cohort of CHC patients with long term follow-up. Moreover, the second aim was to evaluate if a combination of blood fibrosis tests can improve the liver-prognosis assessment.
Patients and Methods
Patients
We used a previously published database of CHC patients with liver biopsy and blood fibrosis tests from the University Hospital of Angers, France (tertiary referral center) (Calès et al, Liver Int 2008; 28:1352-62). Patients were included from 1994 to 2007 if they had CHC, defined as both positive anti-hepatitis C virus antibodies and hepatitis C virus RNA in serum. Exclusion criteria were other causes of chronic hepatitis (hepatitis B or HIV co-infection, alcohol consumption >30g/d in men or >20g/d in women in the 5 years before inclusion, hemochromatosis, autoimmune hepatitis), or cirrhosis complications (ascites, variceal bleeding, systemic infection, hepatocellular carcinoma). Study protocol conformed to the ethical guidelines of the current Declaration of Helsinki and all patients gave informed consent following an IRB approval.
Data Collection During Follow-Up
Follow-up started the day of liver biopsy and ended January 1st, 2011. Patients were followed and received antiviral therapy according to current EASL and AASLD guidelines. Data recorded during follow-up included: achievement of SVR, date of first significant liver-related event, and date and cause of death. SVR was defined as negative hepatitis C virus RNA 6 months after the end of antiviral treatment. The significant liver-related events (SLRE) under consideration were those requiring specific therapy or care: ascites, encephalopathy, jaundice (serum bilirubin >50 μmol/l), occurrence of large esophageal varices (diameter ≧5 mm), variceal bleeding, hepatorenal syndrome, and hepatocellular carcinoma. Large esophageal varices were considered as significant liver-related event because their occurrence significantly impacts treatment and prognosis. Follow-up data were recorded in patient files; for cases lost to hospital follow-up, the patients or their general practitioners were called for follow-up data. Date and cause of death were obtained for all patients by consulting the national death registry. Transplanted patients were censored as dead at the date of liver transplantation. Finally, 2 clinical outcomes were evaluated in the present study: SLRE as primary endpoint (patients were censored at the date of the first event) and liver-related death as secondary endpoint.
Histological Assessment
Liver fibrosis was evaluated according to Metavir fibrosis (F) staging. Initial pathological examinations were independently performed by two senior experts specialized in hepatology, with a consensus reading for discordant cases. The pathologists were blinded for patient characteristics and blood marker results.
Blood Tests
Blood collection was performed within the 7 days around the date of liver biopsy in 87.4% of the cases, and no more than 3 months thereafter in the other cases. All blood assays were performed in the laboratory of Angers hospital.
Blood fibrosis tests—The following blood fibrosis tests were calculated according to published or patented formulas: APRI (Wai et al, Hepatology 2003; 38:518-26), FIB4 (Sterling et al, Hepatology 2006; 43:1317-25), Hepascore (Adams et al, Clin Chem 2005; 51:1867-73), Fibrotest (Biopredictive, Paris, France) (Castera et al, Gastroenterology 2005; 128:343-50), FibroMeterV2G (2nd generation for virus; Echosens, Paris, France) (Leroy et al, Clin Biochem 2008; 41:1368-76). We also calculated CirrhoMeterV2G (Echosens, Paris, France), which shares the same markers as FibroMeterV2G but with specific coefficients targeted for the diagnosis of cirrhosis (Boursier et al, Eur J Gastroenterol Hepatol 2009; 21:28-38). In clinical practice, the results of some blood fibrosis tests are interpreted by using their fibrosis stage classifications that estimate Metavir F stage(s) from the test results. Fibrotest classification includes 8 classes (called here FT1 to FT8), FibroMeterV2G includes 7 classes (FM1 to FM7), and CirrhoMeterV2G 6 classes (CM1 to CM6).
Statistical Analysis
Survival curves were determined using the Kaplan-Meier method and compared with the log rank test. Discriminative ability for the prediction of clinical outcomes was evaluated by the C-index of Harrell. The C-index of Harrell is an extension of the AUROC for time-to-event (survival) data and evaluates the concordance between the predicted risk of event and the observed survival time. Its results varies from 0 to 1: 1 shows a perfect concordance (discriminative power of the risk score), 0.5 shows random prediction, and a value less than 0.5 indicates discrimination in the opposite direction to that expected. For each test evaluated, the 95% confidence interval of the Harrell C-index was calculated by a bootstrap method, using 1,000 random samples (with replacement) of the same size as the original data set. Paired comparisons of the C-indexes were performed using a bootstrap resampling procedure. To compare 2 tests, we calculated the difference between their C-indexes on each of the 1,000 random samples chosen with replacement from the original dataset. Then, the bootstrap 95% confidence interval of the difference was obtained, and we used the estimate of the standard error from this bootstrap distribution to calculate the standardized difference between the C-Indexes and find the corresponding p-value. Most statistical analyses were performed using SPSS version 18.0 software (IBM, Armonk, N.Y., USA) and SAS 9.1 (SAS Institute Inc., Cary, N.C., USA) by a professional statistician (SB).
Results
Patients
The main baseline characteristics of the 373 included patients are detailed in Table 8.
aEither no or unsuccessful antiviral therapy
bOnly the first significant liver-related event during follow-up was taken into account
Median liver biopsy length at inclusion was 20mm (1st quartile: 15mm, 3rd quartile: 25mm). 47 patients (12.6%) experienced at least one SLRE during follow-up. First SLRE were: large varices (n=20), hepatocellular carcinoma (n=13), ascites (n=4), jaundice (n=4), hepatic encephalopathy (n=1), variceal bleeding (n=1), ascites and jaundice (n=2), ascites and jaundice and encephalopathy (n=2). Median follow-up for SLRE was 7.6 years (interquartile range: 4.0-12.1; 2972 person-years). Cumulative incidence of first SLRE at median follow-up was 11.5% (
50 patients died during follow-up; cause of death was unknown in 2 cases and liver-related in 17. Also, 6 patients had liver transplantation. Death was thus considered as liver-related in 23 patients. As death occurred after the first SLRE, the median follow-up for death was longer than for SLRE: 9.5 years (interquartile range: 5.8-13.0; 3508 person-years). Survival without liver-related death and overall survival at median follow-up were, respectively: 95.1% (
37.2% of patients achieved SVR during follow-up, and 62.8% had either no or unsuccessful antiviral therapy. There was no liver-related death during the follow-up of patients who achieved SVR, and only 5 of them had a SLRE. 3 of these 5 patients (2 with Metavir F4 and 1 with F3 at inclusion) developed large oesophageal varices before achieving SVR. One patient (Metavir F2 at inclusion) had large oesophageal varices 5 years after achieving SVR in a context of metabolic syndrome, and thus attributed to an ongoing active nonalcoholic steatohepatitis. The last patient (Metavir F3 at inclusion) had jaundice and ascites before achieving SVR in a context of recovery of excessive alcohol consumption during follow-up.
Prediction of First SLRE
Liver Biopsy Versus Blood Tests
Harrell C-indexes of liver biopsy (Metavir F) and blood fibrosis tests for the prediction of the first SLRE are depicted in Table 9 (paired comparisons between tests are detailed in Table 10).
aComparison of Metavir F vs blood tests: p = 0.039 vs APRI, p = 0.002 vs FIB4, p = 0.005 vs FibroMeterV2G, p >0.340 for other comparisons
bComparison of Metavir F vs blood tests: p = 0.011 vs FIB4, p = 0.078 vs FibroMeterV2G, p >0.190 for other comparisons
0.039
0.002
0.011
0.005
0.049
0.040
0.032
0.053
0.028
Briefly, Metavir F and the 6 blood fibrosis tests had good discriminative ability for the prediction of the first SLRE (with C-indexes >0.8). 5 out 6 blood fibrosis tests had a higher C-index than liver biopsy; this increase being significant for FIB4 (p=0.002), FibroMeter (p=0.005), and APRI (p=0.039).
To evaluate each blood fibrosis test against liver biopsy, we performed successive multivariate stepwise forward Cox models including each blood test with Metavir F, adjusted on age, sex, HCV genotype (1 vs others) and SVR for the prediction of first SLRE (Table 11).
Each blood fibrosis test was an independent predictor of SLRE (at the first step, except for Hepascore). Metavir F was not selected as an independent predictor of SLRE when introduced with FibroMeter or CirrhoMeter in the multivariate analysis. This showed that FibroMeter or CirrhoMeter entirely caught the independent prediction without any significant additional prediction by liver biopsy.
Taken together, all the previous results show that blood fibrosis tests are predictors at least as accurate, and for some significantly better than Metavir F staging for the prediction of SRLE in CHC patients.
Fibrosis Classifications
As our results showed blood fibrosis tests as prognostic markers, they should individualize subgroups of patients with different prognosis. A robust a posteriori determination of these subgroups would require a large cohort (>1000 patients) with high rate of events, followed by external validation. We thus chose to evaluate the prognostic significance of blood tests fibrosis classifications. This was relevant for two reasons. First, fibrosis classifications allow for an a priori stratification of patients. Second, fibrosis classifications have been developed to estimate fibrosis stages and this original design would validate their ability to identify at-risk patients.
How to Predict SLRE in Clinical Practice?
We have already shown in cross-sectional studies that combination of fibrosis tests by binary logistic regression improves the non-invasive diagnosis of liver fibrosis (Boursier et al, Am J Gastroenterol 2011; 106:1255-63). We used the same design in the present longitudinal study to determine the most complementary combination of blood fibrosis tests for the prediction of first-SLRE.
In multivariate forward Cox model including the 6 blood fibrosis tests, age, sex, SVR, and HCV genotype (1 vs others), independent predictors of first SLRE were CirrhoMeter (p<0.001), FibroMeter (p=0.008), and SVR (p=0.010). This result was not surprising as FibroMeter and CirrhoMeter, the first dedicated to significant fibrosis and the latter to cirrhosis, are clearly complementary with CirrhoMeter “stretching” high FibroMeter values (
These results lead us to evaluate a combination of FibroMeter classification for patients included in the FM1-4 classes with CirrhoMeter classification for patients included in FM5-7 classes. This combination provided 10 intermediate classes (FM1-4 and CM1-6 if FM5-7). These 10 intermediate classes were grouped into 5 final classes owing to significantly different prognosis (FM1, FM2/3, FM≧4+CM≦4, FM≧4+CM5, FM≧4+CM6) defining the new FM/CM classification. The rates of patients included in the 5 classes of the FM/CM classification were, respectively from the 1st to the 5th: 16.7%, 41.5%, 27.2%, 10.0%, and 4.6%.
FM/CM Classification Versus Liver Biopsy
As with Metavir F staging, FM/CM classification split patients into 5 classes but provided several advantages. First, the cumulative incidence of first SLRE during follow-up was significantly different between adjacent classes of the FM/CM classification (
Liver-Related Death
Survival without liver-related death was significantly different between Metavir F4 and F3 stages (p<0.001), but not among F0 to F3 stages (p=0.456,
Evaluation of Liver-Related Death in Clinical Practice
In a multivariate forward Cox model including the 6 blood fibrosis tests, age, sex, and HCV genotype (SVR was not introduced in the model as there was no death in SVR patients), CirrhoMeter was the only independent predictor (p<0.001) of liver-related-death. No liver-related death occurred among patients included in CM1 class (
By contrast, survival without liver-related death was significantly different between the 2nd and the 3rd classes, and between the 4th and the 5th classes of the FM/CM classification (
Evaluation of Liver Prognosis in Clinical Practice
To summarize, liver-prognosis in CHC is best assessed by a combination of two complementary tests: the FibroMeter (targeted for fibrosis) in patients with mild disease and CirrhoMeter (targeted for cirrhosis) in patients with advanced disease. Both tests provide a combined FM/CM classification which is the best independent predictor of first SLRE as well as liver-related death. The incidences of 5 and 10-year SLRE or liver-related death as a function of FM/CM classification are presented in
a Either no or unsuccessful antiviral therapy
b no patient in this case
Our results demonstrate that blood fibrosis tests are globally at least as accurate as Metavir F staging for the prediction of liver prognosis in CHC patients. Indeed, none of the blood fibrosis tests had a significantly lower Harrell C-index than liver biopsy for the prediction of SLRE or liver-related death. In multivariate analyses including each blood test with Metavir F staging, almost all blood fibrosis tests were independent predictors of SLRE or liver-related death, suggesting that blood fibrosis tests provide relevant prognostic information. Moreover, FibroMeter and CirrhoMeter were independent predictors of SLRE with no role for Metavir F, as was CirrhoMeter for liver-related death.
Among the 6 blood fibrosis tests evaluated, multivariate analysis identified FibroMeter and CirrhoMeter as the best combination for the prediction of first SLRE. These two tests, the former specifically targeted for significant fibrosis and the latter for cirrhosis (Boursier et al, Eur J Gastroenterol Hepatol 2009; 21:28-38), are very complementary (
FIB4, an inexpensive and easy-to-calculate blood fibrosis test, showed very good prognostic accuracy in our study. It was not selected as an independent predictor of first SLRE by multivariate analysis, probably because it was masked by the high complementary role of FibroMeter and CirrhoMeter for this endpoint. By using its 2 recommended diagnostic cut-offs, FIB4 identified 3 subgroups of patients whose prognosis was significantly different. However, 11 patients from the best-prognosis subgroup experienced SLRE during their follow-up whereas, by comparison, no patients included in the first class of FibroMeter or CirrhoMeter classification had SLRE. This could be due to an inadequate choice for FIB4 thresholds. Indeed, no SLRE occurred during follow-up of patients included in the 3 first deciles of FIB4. This suggests that FIB4 is also able to determine a subgroup of patients with excellent prognosis, and FIB4 thresholds should be refined for this purpose.
In conclusion, blood fibrosis tests are also prognostic tests at least as accurate as liver biopsy for the prediction of liver prognosis in patients with CHC. A combination of two complementary blood tests, one targeted for fibrosis (FibroMeter) and the other for cirrhosis (CirrhoMeter), is more accurate than liver biopsy and appears as a simple and accurate method for the baseline assessment of prognosis in clinical practice.
Number | Date | Country | Kind |
---|---|---|---|
14163948.4 | Apr 2014 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2015/057644 | 4/8/2015 | WO | 00 |