1. Introduction
Chronic liver disease (CLD) is a common clinical condition consisting of a progressive deterioration of liver function related to a continuous process of inflammation, destruction and regeneration of the liver parenchyma [
1]. This process can be due to a wide spectrum of etiologies including prolonged excessive alcohol consumption (alcoholic liver disease), metabolic alterations (non-alcoholic fatty liver disease), infections (chronic viral hepatitis), autoimmune diseases, genetic disorders, and hepatotoxic drugs, among others [
2,
3]. This severe damage causes the liver to lose its ability to repair itself, beginning with hepatosteatosis that leads to fibrosis, cirrhosis and, ultimately, the appearance of hepatocellular carcinoma [
1,
2].
Representing the fifth cause of death in the European Union, CLD is also responsible for high rates of disability and intensive use of healthcare services [
3,
4]. Its prevalence has experienced an increasing trend in recent decades. Approximately, 1.5 billion people suffer from CLD in the world [
4,
5], and its average prevalence in Europe is 0.83% [
6]. According to data from the Global Burden of Disease Study, the age-standardized incidence rate of this disease was 20.7 cases per 100,000 inhabitants in 2015, which represents an increase of 13% since 2000 [
4]. The epidemiology of CLD has undergone a change in recent years, reflecting, on the one hand, the implementation of large-scale hepatitis B vaccination and hepatitis C treatment programs and, on the other hand, the increasing prevalence of metabolic syndrome and excessive alcohol consumption [
3].
The initial symptoms of CLD may initially be nonspecific, including fatigue, anorexia, and weight loss, or may manifest clinically with the development of complications such as esophageal varices, ascites, jaundice, hepatic encephalopathy, hepatorenal syndrome, hepatopulmonary syndrome, and coagulopathies [
2,
7]. Once diagnosed, the goal of CLD treatment is to stop disease progression and complications, which requires a multidisciplinary approach to correct the underlying cause, control portal hypertension, and treat individual signs of the disease. This requires a holistic and integrated approach by an interprofessional team including primary care physicians, specialists in gastroenterology, hepatology and nephrology, liver transplant teams if necessary, nutritionists to provide advice on diet, and social workers and community nurses [
2].
The spectrum of symptoms and conditions of the patient with CLD is not limited to those caused by the disease itself; there may be many other coexisting diseases related to other pathophysiological mechanisms, both concordant and discordant. The coexistence of two or more chronic diseases in the same individual is known as multimorbidity [
8]. Multimorbidity is currently the norm rather than the exception in clinical practices around the world, and it can result in interactions between diseases and between the drugs used to treat them, complicating the clinical management of chronic patients [
9].
Comorbidity of CLD has already been studied in the literature; however, most studies focused on describing the specific comorbidities of some of the etiological diseases of CLD, such as alcoholic liver disease [
10] and nonalcoholic fatty liver [
11], or focus on its complications, such as ascites [
12] and hepatocellular carcinoma [
13]. Predictors of CLD have also been reported, including hypertension, insulin resistance, diabetes mellitus, and obesity. Most previous studies on the comorbidity of CLD were based on the calculation of disease prevalence rates or on regression models to calculate the likelihood of appearance of specific comorbidities based on the presence of CLD. No studies except for one has analyzed the existence of multimorbidity patterns in patients with CLD; however, it only focused on patients with hepatocellular carcinoma [
13]. To the best of our knowledge, there is no study that has analyzed the complete multimorbidity spectrum of CLD as a whole.
The comprehensive analysis of multimorbidity in patients with CLD could help identify the most frequent and associated coexisting diseases, providing us with useful information for their early identification and diagnosis and even prevention, whereas the analysis of their multimorbidity patterns could help identify the profiles of patients susceptible to differentiated clinical management based on their specific pattern of coexisting diseases.
Our population study based on real-life data aimed to comprehensively analyze the comorbidity of CLD by means of the characterization of the most prevalent chronic comorbidities in these patients, the identification of those comorbidities systematically associated with CLD regardless of their prevalence, and the identification and clinical description of the existing multimorbidity patterns in patients with CLD.
2. Materials and Methods
2.1. Study Design and Population
We conducted a retrospective analytical observational study in the EpiChron Cohort, which integrates pseudonymized demographic and clinical information of all the users of the public health system of Aragón (Spain), who represent approximately 95% of its reference population (1.3 million inhabitants). The baseline characteristics of the EpiChron Cohort can be found elsewhere [
14]. The conformation of this cohort for research on the epidemiology of multimorbidity and chronic diseases was favorably evaluated by the Research Ethics Committee of the Community of Aragón (CEICA; approval code PI17/0024). Given the epidemiological nature of the study, which used anonymized data that was presented only at an aggregate level, the obligation to obtain informed consent from the patients was waived by CEICA.
For this study, we selected as reference population all the people aged 18 years and older from the cohort who were registered as users of the public health system in 2015 and who were also registered during the previous year to have at least one year of complete information for their demographic and clinical characterization. Of these, we selected all the patients with a diagnosis of CLD recorded in their electronic health records for the study of their multimorbidity.
2.2. Variables and Data Sources
For all the patients included in the study, we analyzed their sex and age as of January 1, 2015, and all their chronic comorbidities registered in both primary and hospital care.
The diagnoses were initially coded using the International Classification of Primary Care, 1st edition (ICPC-1), in the case of primary care, and the International Classification of Diseases, 9th edition, Clinical Modification (ICD-9-CM), in the case of hospital care. We grouped these codes into broader mutually exclusive diagnostic categories called Expanded Diagnostic Clusters (EDCs) based on their clinical and diagnostic similarities using the Johns Hopkins Adjusted Clinical Groups (ACG
®) [
15]. The purpose of using this software was to transform the thousands of potential different codes into a more manageable list of 264 diagnoses that avoid duplicate diagnoses and facilitate the counting of diseases. From this list, only the 114 EDCs proposed as chronic in the study by Salisbury et al. [
16] were taken into account. Chronic diseases were defined as those diseases that typically last six months or more, including past conditions that require current care, conditions with a high risk of recurrence, or past conditions that have ongoing implications for patient management. We grouped some diagnoses so that a total of 110 conditions were finally analyzed; the list of diseases used is presented as supplementary material (
Table S1). Subjects with CLD were identified as those having the EDC code “GAS05: Chronic liver disease”.
All the information was gathered from patients’ electronic medical history of primary care, hospital discharge reports (through the Minimum Basic Data Set; CMBD), and clinical-administrative databases (user database, BDU).
3.3. Statistical Analysis
For the comprehensive analysis of the multimorbidity of patients with CLD, we conducted a three-phase analysis. First, we carried out a descriptive analysis of the demographic and clinical characteristics of the study population. For the identification of the most frequent comorbidities, we analyzed the prevalence of each comorbidity in the total population and stratified by sex and age (i.e., 18-44, 45-65, and >65 years), which was presented as frequencies and percentages.
Secondly, for the identification of those comorbidities systematically associated with CLD, we performed logistic regression models to calculate the risk of occurrence of each comorbidity (dependent variable) based on the presence or absence of CLD (independent variable). As a result, we obtained odds ratios (ORs) accompanied by their 95% confidence intervals (95% CI), which were calculated as unadjusted crude ORs and as age- and sex-adjusted ORs. For the comparison of adjusted ORs, we used the Bonferroni correction method for multiple comparisons (for a total of 85 disease comparisons with at least five cases in both men and women), establishing statistical significance at p <0.00059.
Finally, for the analysis of the existence of associations among diseases in the form of multimorbidity patterns in patients with CLD, we carried out an exploratory factor analysis. This technique allowed us to reduce the number of explanatory variables in the data set to a smaller number of latent variables, which can be interpreted as groups of variables (in our case diseases) that share a common underlying causal factor. This analysis was carried out in each sex and age group. To increase the clinical relevance of the results and facilitate their interpretation, we only included in the analysis those comorbidities with a prevalence greater than 5%. Each comorbidity was represented in the analysis as a binary variable indicating present (1) or absence (0), which allowed for the construction of a tetrachoric correlation matrix. To extract the factors, we used the main factor method, and to facilitate its interpretation, we applied an oblique rotation of the factors (Oblimin). To choose the number of factors (i.e., patterns) to extract, we used scree plots of the eigenvalues from the correlation matrix ordered in descending order. This visualization allowed us to identify the inflection point, which indicates the optimal number of factors to extract. This statistical criterion was accompanied by the clinical assessment of the different results.
In each pattern, we included those chronic diseases with a loading factor or factor score (a value that ranges between -1 and 1, representing the strength of association of each diagnosis within each pattern) equal to or greater than 0.30, allowing each disease to be in more than one pattern. To determine the degree of suitability of the sample to use factor analysis, we performed the Kaiser–Meyer–Olkin (KMO) test; the value of this parameter varies between 0 and 1, with values closer to 1 representing greater goodness of fit.
Once the patterns in each subpopulation were obtained, and to facilitate their clinical interpretation, we named the patterns based on the composition of diseases and the names already described in the bibliography with the help of two clinical care specialists in family and community medicine from the research group.
All statistical analyzes were carried out in Stata (Version 12.0, StataCorp LLC, College Station, TX, US) and R (Version 3.6.3, The R Foundation for Statistical Computing, Vienna, Austria).
4. Discussion
In this study, we observed that multimorbidity is the norm rather than the exception in patients with CLD, which is frequently accompanied by other diseases, like arterial hypertension, lipid metabolism disorders, diabetes, obesity, and musculoskeletal diseases, with some relevant differences between sexes, such as the higher prevalence in men of disorders related to substance use from an early age, as well as of depression and cancer. This multimorbidity seems to cluster in three different multimorbidity patterns (Cardiovascular, Metabolic-Geriatric, and Mental-Substance use) in both men and women, although with slight differences in their composition depending on sex.
The prevalence rates of CLD obtained in our study are similar and consistent with the average prevalence of CLD described in the European Union, which is around 0.83% [
6], and with the fact that excessive alcohol consumption, one of the most common etiologies of CLD, is more common in men than in women [
19].
The prevalence rate of multimorbidity is higher than those observed in previous articles reporting multimorbidity rates of 80% and a lower number of comorbidities [
11,
20]; this is probably due to the fact that our study exhaustively analyzed all possible chronic conditions, and it was not based on a limited number of them. It has been demonstrated that the risk and severity of non-alcoholic liver disease increase with the number of components of metabolic syndrome present, and obesity is considered as the biggest risk factor for it [
11]. Similar results were observed in people with non-alcoholic liver disease in Russia; almost 80% had at least two metabolic comorbidities, with the most common ones being overweight/ obesity (81%), hypercholesterolemia (75%), and type 2 diabetes (17%) [
11]. In a population-based study in Sweden, high blood pressure was more prevalent (33%) than type 2 diabetes (29%) and obesity (24%). In the Manitoba Follow-up Study, hypertension, insulin resistance/ diabetes mellitus, and obesity were reported as greatest predictors of the appearance of CLD [
21]. However, the majority of studies focused on describing the specific comorbidities of some of the etiological diseases of CLD, such as alcoholic liver disease [
10] and nonalcoholic fatty liver [
11], or on its complications, such as ascites [
12] and hepatocellular carcinoma [
13].
Surprisingly, 75 conditions of the original list of 110 were associated with CLD in our study, and their prevalence was higher than expected. This could be due to the multifactorial etiology of CLD and to the fact that the liver represents an important organ with many functions [
2], whose alteration and damage could in turn be the cause of many other chronic disorders. In fact, among the conditions that showed a higher degree of association with CLD, we observed conditions that are both possible causes (e.g., inherited metabolic disorders, chronic viral infection, substance use, autoimmune diseases, diabetes, obesity [
2,
3]) and potential consequences (e.g., coagulation disorders, thrombosis, gastroesophageal reflux, chronic renal/respiratory failure) of CLD [
2,
22].
Although the relationship between CLD and comorbidities has been studied, there are not many studies that analyze the association of these comorbidities in the form of patterns. As well as we have observed the presence of three patterns, Mu et al. identified three main comorbidity patterns in patients with hepatocellular carcinoma, which included: 1) cirrhosis, hepatitis B, portal hypertension, and ascites; 2) hypertension, diabetes mellitus, coronary heart disease and cerebral infarction; and 3) hypoproteinemia, electrolyte disorders, gastrointestinal bleeding, and hemorrhagic anemia [
13].
Of the three patterns analyzed, the Cardiovascular pattern has been widely described in the literature and is one of the most consistent in the general population [
23,
24]. Although it had not been specifically described in patients with CLD, it is one of the main cause of hospital admissions [
12] and death in patients with non-alcoholic fatty liver disease [
25,
26]. In Germany, it was observed that among hospital admissions of CLD patients, diagnosis of respiratory diseases with infection had the highest mortality rate of 21.6% followed by cerebrovascular disease with a rate of 15.5% [
12]. The development of these cardiovascular complications are related to obesity [
12,
25]. It is thought that the hyperdynamic circulation in cirrhosis provides some protection against overt heart failure, atherosclerosis, and ischemic events, but peripheral arterial disease, acute myocardial infarction, and heart failure were predictors of mortality in the CirCom study [
22].
The Metabolic-Geriatric pattern is also one of the most consistently described in the general population [
23], and in the specific case of CLD it makes sense that there is a profile of patients with this type of multimorbidity, since the metabolic syndrome and the disease of the associated fatty liver disease is one of the most common causes of CLD [
11,
27]. We observed that this pattern was mainly composed of obesity, diabetes and lipid metabolism disorder (due to low HDL and high triglycerides), that are considered metabolic risk factors of CLD, and diabetes is a predictor of severe outcomes [
27]. Non-alcoholic fatty liver disease is recognized as the hepatic manifestation of the metabolic syndrome that represents a cluster of metabolic abnormalities, such as hyperlipidemia, glucose intolerance, obesity, and systemic hypertension [
11], comorbidities that we found associated to CLD. The fact that this pattern is associated with geriatric or aging diseases could be due to the fact that older CLD patients are precisely those who have had this metabolic etiology, and not the one related to excessive consumption of alcohol and related factors that may cause these patients to have a lower life expectancy [
28]. However, due to the technique used (factor analysis), which groups diseases and not people, it is not possible to know the average age of the patients of this multimorbidity pattern to corroborate this hypothesis.
The Mental-substance use pattern was composed by substance use, depression, sleep disorders, and neuritis has been previously described in the literature, but mainly in young and adult men, while in women a mental pattern with anxiety and depression is usually observed but not associated with substance use [
23,
24]. As in the previous case, this pattern of multimorbidity makes clinical sense since it would represent another of the most frequent etiologies of CLD: alcoholic liver disease due to excessive alcohol consumption [
10,
22]. Knowing that this pattern is also associated with problems such as depression and sleep disorders is important in order to characterize the specific needs of these patients [
29] and proactively seek their appearance and diagnosis to prevent worsening health outcomes. In many cases, depression is associated to stigmatization in patients with liver disease, and it is associated to a lack of social support, and a decrease in the tendency to seek health care [
29]. In the development of CirCom score, Jepsen et al. considered that mental disease could be a predictor of mortality due to its association with substance abuse and suicide risk, but they observed that schizophrenia was indeed an adverse prognostic factor, but depression and bipolar disorder were not associated with mortality [
22]. They did not find other studies that have examined the prognostic impact of psychiatric diseases in cirrhosis [
22].
The negatives outcomes of CLD patients are not only due to their comorbidities and clinical complexity, but also to biological (i.e., ageing, frailty, multimorbidity, mental disease, dependency, malnutrition) and non-biological (i.e., socioeconomic, behavioral, environmental, cultural) variables [
30]. Lifestyle factors, including smoking, alcohol, physical inactivity, adiposity [
26], and poor diet [
25], played a key role in the incidence of CLD as well as to their complications. That is the reason why promotion of lifestyle interventions among CLD patients in early stage of disease course is necessary, in order to prevent cardiovascular risk factors [
25], transitions to metabolic complications and death [
26]. Also, social factors like limitations in daily living due to their disease, loneliness, low income, stigmatisation, and isolation, play an important role in the evolution of these disease decreasing the quality of life of these patients [
29]. Patients highlighted the need for information to understand and manage CLD, and awareness and support from healthcare professionals to better cope with the disease [
29]. An interprofessional team that provides a holistic and integrated approach is needed for a CLD patient to achieve the best possible outcomes. The earlier signs and symptoms of a complication are identified, the better is the prognosis and outcome [
2].
One of the most important limitations of our study was related to its cross-sectional design, which did not allow the establishment of causal relationships between the existence of CLD and the comorbidities analyzed, nor to study the order of appearance or time of progression of the diseases within each multimorbidity pattern. Furthermore, it was not possible to analyze variables of interest for the interpretation of the results, such as the socioeconomic and educational level of the patients, genetic variables, or information on lifestyles (e.g., physical activity, consumption of alcohol, tobacco), which were not available in the EpiChron Cohort. On the other hand, the factor analysis presented limitations for the study of the different sex and age strata, and the analysis had to be done stratified only by sex. Lastly, given the multifactorial etiology of CLD, it would be interesting to carry out future studies that include various subanalyses according to the origin of the disease (e.g., alcoholism, fatty liver, viral infection) to know whether this influences the formation and composition of multimorbidity patterns.
One of the main strengths of this study lies in its population-based nature, since all cases of CLD in the reference population were analyzed. Furthermore, we exhaustively analysed the comorbidity of CLD based on the analysis of virtually all chronic diseases (and not only the most relevant or frequent ones) contained in the clinical history and recorded by a professional (and not self-referred by the patients or from surveys).