Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 06 January 2026

Sec. Cardiovascular Endocrinology

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1651689

Association of estimated liver fibrosis with carotid but not femoral atherosclerotic burden: the ILERVAS cohort

  • 1Endocrinology and Nutrition Department, University Hospital Arnau de Vilanova de Lleida, Lleida, Spain
  • 2Obesity, Diabetes and Metabolism (ODIM) Research Group, Lleida Biomedicine Research Institute (IRBLleida), University of Lleida, Lleida, Spain
  • 3Department of Cell Biology, Physiology, and Immunology, Maimonides Institute for Biomedical Research of Córdoba (IMIBIC), University of Córdoba (UCO), Biomedical Research Network center on the Pathophysiology of Obesity and Nutrition (CIBERobn), Instituto de Salud Carlos III, Córdoba, Spain
  • 4Vascular and Renal Translational Research Group, Lleida Biomedicine Research Institute (IRBLleida), Renal Research Network (RICORS2040, ISCIII), Lleida, Spain
  • 5Department of Experimental Medicine , University of Lleida, Lleida, Spain
  • 6Department of Experimental Medicine, University of Lleida, Lleida Biomedicine Research Institute (IRBLleida), Lleida, Spain
  • 7Translational Research in Respiratory Medicine, University Hospital Arnau de Vilanova and Santa Maria, Lleida Biomedicine Research Institute (IRBLleida), Biomedical Research Network Center of Respiratory Diseases (CIBERes), Institute of Health Carlos III, Lleida, Spain
  • 8Department of Endocrinology and Nutrition, Hospital de la Santa Creu i Sant Pau, Biomedical Research Network Center on Diabetes and Metabolic Diseases (CIBERdem), Instituto de Salud Carlos III, Autonomous University of Barcelona, Barcelona, Spain
  • 9Endocrinology and Nutrition Department, Parc Taulí Hospital Universitari, Institute for Research and Innovation Parc Taulí-Centers of Research in Catalonia (IPT-CERCA), Sabadell, Spain
  • 10Medicine Department, Universitat Autònoma de Barcelona, Sabadell, Spain
  • 11Division of Endocrinology, Diabetes and Metabolism, University of Florida, Gainesville, FL, United States
  • 12Endocrinology and Nutrition Department, University Hospital Vall d'Hebron, Barcelona, Spain
  • 13Diabetes and Metabolism Research Group, Institut de Recerca Vall d'Hebron (VHIR), CIBER of diabetes and metabolic diseases (CIBERdem), Autonomous University of Barcelona, Barcelona, Spain

Introduction: Advanced liver fibrosis, a key complication of metabolic dysfunction-associated steatotic liver disease, has been increasingly linked to extrahepatic conditions, including type 2 diabetes, obesity, and cardiovascular disease. However, the specific association of liver fibrosis in the development and progression of subclinical atheromatous disease across vascular territories remains poorly understood. This study evaluates the utility of two non-invasive indices to predict liver fibrosis and their associations with subclinical atheromatous plaque burden and distribution.

Methods: Atheromatous plaque burden (plaque presence, number, and total area) was assessed in the carotid and femoral territories via ultrasonography in 3,981 middle-aged participants without known cardiovascular disease, diabetes, or liver disease from the ILERVAS cohort (ClinicalTrials.gov Identifier: NCT03228459). The fibrosis-4 (FIB-4) and the NAFLD Fibrosis Score (NFS) were evaluated. FIB-4 risk categories were defined as low (<1.30), intermediate (1.30–2.67), and high (>2.67).

Results: Participants in the intermediate and high-risk FIB-4 categories exhibited a higher prevalence of carotid atheromatous disease (56.8% vs. 49.5%, p<0.001), a greater number of plaques (p<0.001), and a larger total plaque area (p=0.007). Multivariable analyses confirmed FIB-4 as an independent predictor of carotid plaque burden (OR: 1.14, 95% CI 1.05-1.24, p=0.003), even adjusting for traditional cardiovascular risk factors. Moving from low to high FIB-4 cut-offs was associated with 12.6% higher odds of carotid atherosclerosis. NFS was also independently associated with carotid atheromatosis (OR 1.10, 95% CI 1.05–1.15, p<0.001). No significant associations were found in the femoral territory for either index.

Conclusions: Estimated liver fibrosis, particularly FIB-4, is a valuable marker for identifying carotid subclinical atherosclerosis in populations without known liver disease. These findings highlight the importance of vascular territory-specific evaluations and support their utility in integrated liver and cardiovascular risk assessment strategies.

Introduction

Liver fibrosis, a hallmark of progressive metabolic dysfunction-associated steatotic liver disease (MASLD), is emerging as a key determinant of cardiovascular (CV) risk (1). While MASLD remains the most prevalent chronic liver disease worldwide, affecting approximately 25% of the global population (2), it is the progression to liver fibrosis that seems to carry the most significant prognostic implications beyond hepatic complications (1, 3, 4). Fibrosis severity has been linked not only to advanced liver outcomes, such as cirrhosis and hepatocellular carcinoma, but also to extrahepatic conditions, including insulin resistance, type 2 diabetes mellitus (T2DM), hypertension, and atherogenic dyslipidemia (5). CV disease is recognized as one of the leading causes of mortality in patients with MASLD (6). Growing evidence suggests that liver fibrosis, along with liver dysfunction, contributes independently to this increased CV risk through systemic metabolic alterations, chronic inflammation, and endothelial dysfunction (7, 8). In addition, pro-inflammatory cytokines and oxidative stress associated with fibrosis may enhance vascular changes leading to a more severe atheromatous plaque formation (9). However, whether liver fibrosis represents an independent risk factor for subclinical atheromatosis distinct from other metabolic comorbidities such as T2DM, remains an open question (10). Similarly, it is still unclear whether assessing liver fibrosis itself, independently of a formal MASLD diagnosis, might offer an opportunity to identify individuals at increased CV risk.

Early stages of atherogenesis are characterized by the hyperplasia and expansion of adventitial vasa vasorum (VV) into the avascular intima (11). Although MASLD severity, as assessed by liver biopsy, has shown no clear association with VV expansion in individuals with severe obesity (12), the specific role of fibrosis in promoting atheromatosis has not been fully explored (13). Some studies suggest that liver fibrosis, more than steatosis itself, could correlate with increased atherosclerotic plaque burden in various vascular territories, including carotid, coronary, and femoral arteries (1314). In fact, it has been shown that liver fibrosis, assessed by the fibrosis-4 score (FIB-4), identifies acute coronary syndrome patients not only at higher risk of in-hospital mortality but also at increased risk of heart failure after discharge (15). In addition, longitudinal studies report that patients with steatosis are more likely to develop coronary arterial calcification (CAC) or incident carotid plaques over time (1618). However, other studies have found no significant association between MASLD and surrogate markers of CV disease, including carotid intima-media thickness, aortic stiffness, brachial artery vasodilatory function, and CAC (811).

Based on the hypothesis that advanced liver fibrosis may be associated with a greater burden of subclinical atheromatous disease, we aimed to evaluate the relationship between two non-invasive liver fibrosis indices, the FIB-4 index and the NAFLD fibrosis (NFS) score, and subclinical atheromatosis (plaque presence, affected territories, number, and total plaque area) in a large population of individuals with low to moderate CV risk from the province of Lleida, Spain (1921).

Methods

Study population, metabolic status, and selection of patients

The ongoing prospective ILERVAS study (ClinicalTrials.gov Identifier: NCT03228459) examines the progression of subclinical atherosclerotic disease (SAD) in individuals with low to moderate cardiovascular risk (22, 23). Of 8,330 participants recruited between January 2015 and December 2017, 3,940 had complete FIB-4 data and 3,981 NFS.

The inclusion criteria were: age 45–70 years, absence of previous cardiovascular disease (angina, myocardial infarction, stroke, peripheral artery disease, heart failure, or any vascular surgery/procedure), and the presence of at least one cardiovascular risk factor [dyslipidemia, arterial hypertension, obesity, smoking, or a first-degree relative with premature cardiovascular disease (myocardial infarction, stroke, or peripheral arterial disease] before the age of 55 years in men or 65 years in women)]. Exclusion criteria included any form of diabetes mellitus, chronic hepatitis or kidney disease, active neoplasia, hematologic disorders, life expectancy less than 18 months, or pregnancy. Regarding liver disease, participants with known chronic hepatopathies (including hepatitis B or C infection, cirrhosis, liver cancer, or other chronic liver disorders), identified through clinical history and diagnostic codes (ICD-9/ICD-10, as applicable at the time of recruitment) in electronic medical records, were excluded from the cohort. Prediabetes was defined according to the ADA criteria: fasting plasma glucose of 100–125 mg/dL or HbA1c 5.7–6.4% (24).

The prescribed treatments for hypertension (ACE inhibitors, diuretics, angiotensin receptor blockers, beta-blockers, and calcium channel blockers) and lipid management (statins, fibrates, ezetimibe, and omega-3 fatty acids) were extracted from prescription data provided by the Catalan Health Service.

Evaluation of clinical variables

Height and body weight were measured with participants wearing light clothing without shoes. Total cholesterol levels (mg/dL) were determined from non-fasting dried capillary blood samples obtained via fingertip puncture and analyzed with the REFLOTRON® Plus system (Roche Diagnostics, GmbH, Germany) (22, 23). A complete lipid profile was measured in participants with total cholesterol levels ≥ 200 mg/dL after fasting for at least 6 hours, or ≥ 250 mg/dL regardless of fasting status. Pulse pressure was calculated as the difference between systolic blood pressure (SBP) and diastolic blood pressure (DBP). Smoking status was categorized as non-smoker, current smoker, or former smoker, with the latter defined as participants who had quit smoking one year prior.

Non-invasive methods for predicting liver fibrosis

Following a literature review to identify predictive models of liver fibrosis validated in MASLD and based on routinely available clinical data and serum markers, two indices were included in our study: the FIB-4 index and the NFS (1921). The formulas and risk categories used to calculate these scores were: (i) FIB-4 index = (Age × AST [U/l]/(platelet count [109/L] × √ALT (U/l), with the risk categories defined as <1.30 for low risk, 1.30-2.67 for intermediate risk, and >2.67 for high risk; and (ii) NFS = −1.675 + (0.037 × age) + (0.094 × BMI [kg/m²]) + (1.13 × diabetes [yes = 1, no = 0] + (0.99 × AST/ALT ratio) – (0.013 × platelet count [× 109/l]) – (0.66 × albumin [g/dl]), with the risk categories defined as < -1.455 for no fibrosis, -1.455-0.675 for possible fibrosis, and >0.675 for suggestive of fibrosis.

Assessment of atheromatous plaque burden

The bilateral carotid arteries (common, bifurcation, internal, and external) and femoral arteries (common and superficial) were examined using a Vivid-I Doppler Ultrasound system (General Electric Healthcare, Waukesha, WI, USA) equipped with a 12L-RS broadband linear probe operating at frequencies between 5 and 13 MHz. Images were acquired by trained sonographers following standardized and validated protocols (25). The sonographers were blinded to the participants’ clinical histories to avoid potential bias.

Subclinical atheromatosis was defined as the presence of any plaque in the twelve evaluated vascular regions (26). A plaque was well-defined as a focal intima-media thickness ≥ 1.5 mm projecting into the arterial lumen (27). All detected plaques were measured, and the total plaque area (cm²) was calculated (28).

Statistical analysis

Non-normal variables were expressed as the median [interquartile range]. Categorical variables were reported as absolute frequencies. Group comparisons were performed using the Pearson’s Chi-squared test for categorical variables, whilst Mann-Whitney U test and Kruskal Wallis test were applied for quantitative variables. The association between continuous variables was evaluated using the Spearman correlation test. In line with the guideline-endorsed two-step interpretation of FIB-4 (<1.3 ‘rule out’, ≥1.3 ‘rule in’) and because only a small number of participants had high FIB-4 values (n=61), the moderate and high fibrosis categories were combined for analysis (29, 30).

A multivariable logistic regression model was developed to assess the presence of SAD, incorporating the following potential confounding factors: sex, prediabetes status, LDL cholesterol, triglycerides, pulse pressure, smoking status, and each non-invasive liver fibrosis predictor (FIB-4 or NFS) separately. The constituent components of each score were not included as separate covariates in the models to minimize multicollinearity. A multivariable lineal regression model was also performed to evaluate the number of plaques in the different territories and the total plaque area using the same covariates. Furthermore, stratified analyses were conducted by sex to explore potential sex-specific associations between FIB-4 and carotid plaque burden, applying the same set of confounding variables as in the main models. Multicollinearity between FIB-4 and covariates, including age, was assessed using the variance inflation factor (VIF), and no significant collinearity was detected (all VIF values < 5). Linearity of the association between FIB-4 and plaque was examined by fitting logistic regression models including both a quadratic term and restricted cubic splines (4 df). Model fit was assessed using likelihood ratio tests (quadratic vs. linear: LR χ²(1)=0.0004, p=0.98; spline vs. linear: LR χ²(3)=2.42, p=0.49) and the Akaike Information Criterion (AIC; 1327.4 vs. 1323.8, respectively). As no departures from linearity were detected, FIB-4 was modeled as a continuous linear predictor. Missing data were handled using a complete-case approach, whereby participants with missing values in variables included in each analysis were excluded. The proportion of missing data was low (<5% in all key variables), and therefore the potential impact on the results is expected to be minimal. All statistical analyses were conducted using STATA v.16 software, with a significant threshold set at 0.05.

Ethical considerations

The ILERVAS study protocol received approval from the Ethics Committee of the University Hospital Arnau de Vilanova (Initial approval: CEIC-1410, 19/12/2014). All participants provided written informed consent prior to inclusion in the study. The research was conducted in accordance with the Declaration of Helsinki.

Results

The main clinical and metabolic characteristics of the participants, stratified by FIB-4 index thresholds, are summarized in Table 1 and Figure 1. According to the FIB-4 classification, most participants were categorized as low-risk for fibrosis (78.3%). As expected from the formula used to calculate FIB-4, these individuals also exhibited lower platelet counts and impaired liver function test results. In addition, participants with higher FIB-4 values showed increased pulse pressure and a greater use of antihypertensive medications. Interestingly, despite these unfavorable cardiovascular risk factors, they displayed a less atherogenic lipid profile, with no significant differences in the use of lipid-lowering medications.

Table 1
www.frontiersin.org

Table 1. Main clinical characteristics, metabolic data and medical treatment in the study population according to the FIB-4 index thresholds.

Figure 1
Bar graph comparing plaque prevalence in two FIB-4 risk groups: Intermediate-High and Low. For any plaque, 74.9% in Indeterminate-High and 71.7% in Low (p=0.063). For carotid plaque, 56.8% in Indeterminate-High and 49.5% in Low (p<0.001). For femoral plaque, 56.4% in Indeterminate-High and 53.5% in Low (p=0.127).

Figure 1. Prevalence of atheromatous disease according to the FIB-4 index thresholds. Affected territories include bilateral carotid (common, bifurcation, internal, and external) and femoral (common and superficial) arteries. FIB-4 risk categories are defined as <1.30 for low risk, 1.30-2.67 for intermediate risk, and >2.67 for high risk.

Regarding atherosclerotic disease, individuals with the higher FIB-4 values demonstrated a significantly greater number of affected vascular territories [2.3 (0 to 3) vs. 1.9 (0 to 3); p < 0.001] and a higher total plaque area [0.6 (0.2 to 1.1) cm² vs. 0.4 (0.2 to 0.9) cm²; p = 0.003] (Table 1). These differences were primarily driven by carotid artery involvement. Specifically, participants in the high-risk group exhibited an approximately 10% higher prevalence of carotid artery disease (56.8% vs. 49.5%, p <0.001), a greater number of carotid plaques [1.2 (0 to 2) vs. 1.0 (0 to 2), p <0.001], and a larger carotid plaque area (p = 0.007). A similar pattern of increased plaque in the carotid arteries was observed in participants with higher FIB-4 values compared to those at low risk, affecting both the right (42.8 vs 37.1%, p=0.003) and left (40.8 vs. 35.5%, p=0.005) sides. Conversely, no significant differences between groups were observed in the femoral territories.

Univariate logistic regression analysis identified FIB-4 as an independent predictor of carotid plaque presence [OR 1.30 (95% CI 1.13 to 1.49), p < 0.001] but not for femoral plaques (Table 2). Further analysis using multivariate logistic regression, confirmed that FIB-4, along with established cardiovascular risk factors such as smoking, prediabetes, elevated pulse pressure, increased LDL cholesterol and male gender, independently predicted the presence of carotid plaques [OR 1.14 (1.05 to 1.24), p =0.003] (Table 3). Using clinical cut-offs, moving from the low-risk threshold (1.30) to the high-risk threshold (2.67) was associated with approximately 12.6% higher odds of carotid atherosclerosis. When analyses were stratified by sex, the association between FIB-4 and the number of carotid plaques remained significant in men [OR 1.18 (1.03 to 1.36), p =0.018] but not in women [OR 1.11 (0.99 to 1.24), p =0.067]. Full results of the sex-stratified multivariable models are provided in Supplementary Tables 1 and 2.

Table 2
www.frontiersin.org

Table 2. Univariable logistic regression model for the analysis between the presence (yes/no) of plaque and the non-invasive methods.

Table 3
www.frontiersin.org

Table 3. Logistic regression model for the analysis between the presence (yes/no) of carotid plaque and FIB4-index.

When using the NFS, participants with fibrosis-suggestive values exhibited a higher overall prevalence of atheromatous plaques (75.3% vs. 70.5%, p = 0.002), mainly driven by greater carotid involvement. Specifically, they showed an increased prevalence of affected carotid territories (55.4% vs. 47.9%, p < 0.001) and a higher number of plaques in this region (1.15 vs. 0.95, p < 0.001) (Table 4). In the unadjusted logistic regression analysis the NFS was identified as an independent predictor for the presence of any plaque, and specifically within the carotid territory (Table 2). In the multivariable logistic model the NFS persisted as an independent predictor for the presence of carotid [OR 1.10 (1.05 to 1.15), p <0.001] but no femoral plaque (Table 5).

Table 4
www.frontiersin.org

Table 4. Atheromatous disease characteristics in the study population according to the NAFLD fibrosis score thresholds.

Table 5
www.frontiersin.org

Table 5. Logistic regression model for the analysis between the presence (yes/no) of carotid plaque and NAFLD fibrosis score.

Discussion

Our study demonstrates that individuals in the intermediate and high-risk FIB-4 category exhibit a significantly higher overall plaque burden compared to those in the low-risk group, a finding that is largely driven by a greater prevalence of carotid, but not femoral, SAD. This finding aligns with prior research suggesting a link between MASLD and CV disease but provides a novel perspective by focusing on fibrosis-specific associations across vascular territories. Although several studies have explored the impact of hepatic steatosis, often estimated using the fatty liver index, on CV outcomes, the role of advanced liver fibrosis as a contributor to atherosclerotic disease across different vascular sites has remained largely underexplored (3136).

Our findings highlight that liver-related fibrosis, as estimated using FIB-4 thresholds, is associated with differential patterns of atherosclerosis. Specifically, participants with intermediate and high-risk FIB-4 values showed higher prevalence, number and area of carotid plaques, while no significant associations were observed in the femoral territory. Although the observed 12.6% increased odds of carotid atherosclerosis associated with intermediate-high FIB-4 values reached statistically significance, the magnitude of this effect is moderate and should be interpreted in the context of overall clinical risk and multifactorial disease etiology. At the same time, this divergence underscores the need for tailored risk assessment strategies based on vascular site-specific mechanisms. Factors associated with chronic liver disease and leading to advanced liver fibrosis may contribute to carotid plaque formation, potentially mediated by insulin resistance, cardiometabolic risk factors linked to metabolic dysfunction and chronic subclinical inflammation (3739). Notably, metabolic factors such as prediabetes and pulse pressure emerged as significant predictors in the multivariate model, reinforcing their role in carotid atherogenesis. Additionally, our findings raise the possibility that systemic fibro-inflammatory mediators released by the fibrotic liver may contribute to endothelial dysfunction in specific vascular territories. Experimental studies have shown that fibrotic livers release cytokines, chemokines, adhesion molecules, and extracellular vesicles, as well as lipids (mainly triglycerides), that promote systemic endothelial dysfunction (40). Moreover, the crosstalk between liver sinusoidal endothelial cells and the hepatic microenvironment is increasingly recognized as a key driver of fibrogenesis and systemic pro-inflammatory signaling (41). Whether such mechanisms preferentially affect certain vascular beds remains an open question that requires dedicated mechanistic studies.

The differential impact of fibrosis on carotid versus femoral atheromatosis may be explained by anatomical differences between these vascular territories. Femoral arteries have a distinct microscopic structure and flow dynamics compared to carotid arteries, which may make them less susceptible to fibrosis-driven mechanisms (42, 43). These findings are consistent with preclinical studies demonstrating variability in plaque composition, microscopic structure, and susceptibility across vascular territories (4244).

Understanding the factors that influence atherosclerosis at different anatomical sites is clinically important, as these differences directly impact the rate and nature of clinical events. For instance, the localization of SAD plays a pivotal role, with carotid arteries being more prone to foam cell lesions and lipid-rich plaques than femoral arteries. Accordingly, lipid-rich plaques are significantly more common in the carotid territory, particularly among patients who die from coronary atherosclerosis (42). Non-invasive magnetic resonance imaging studies have also demonstrated notable differences in plaque composition between vascular territories, revealing that carotid plaques tend to exhibit larger necrotic cores and more frequent hemorrhagic areas compared to those in the femoral arteries (43). Therefore, our findings align with prior research, such as the PDAY (Pathological Determination of Atherosclerosis in Youth) study, which highlighted that cardiovascular risk factors exert variable effects across different vascular sites (44).

Our findings also add to the growing evidence suggesting variability in the predictive power of liver fibrosis indices. FIB-4 index emerged as a strong independent predictor of carotid plaque presence and burden, while the NFS was also independently associated with carotid atheromatosis, albeit with somewhat weaker effect estimates. This discrepancy likely arises from differences in the indices’ sensitivity to advanced fibrosis versus metabolic factors, such as BMI and diabetes, which are integral to the NFS formula (1921). Beyond liver- and site-specific mechanisms, our findings also support the interpretation of FIB-4 and NFS as integrative markers of cumulative, age-related cardiometabolic and fibrotic burden. In our cohort, the intermediate–high FIB-4 category was significantly older and exhibited a more adverse cardiometabolic profile. This observation is consistent with histology-based data in MASLD, where the prevalence and severity of fibrosis increase in a stepwise fashion with the number of metabolic comorbidities, and the count of these comorbidities independently predicts significant fibrosis (45). Similarly, hypertension [1.92 (1.17 to 3.16), p <0.001] and diabetes mellitus [2.00 (1.22 to 3.28), p <0.001] significantly contributed to advanced fibrosis (≥ F3) on multivariate analysis among 458 patients diagnosed with MASLD via liver biopsy (46). Taken together, these observations reinforce the utility of fibrosis scores, particularly FIB-4 and NFS, as reliable, non-invasive markers for identifying individuals at risk of subclinical carotid atherosclerosis, even in populations without overt liver or cardiovascular disease.

Atherosclerosis is a multifaceted process influenced by genetic, metabolic, and environmental factors. The interplay between systemic inflammation, oxidative stress, and metabolic dysfunction likely drives the progression of carotid atheromatosis. Emerging genetic evidence supports this notion in MASLD, showing that variants like PNPLA3 and TM6SF2, while strongly associated with hepatic fat and fibrosis, do not independently increase CVD risk in the absence of metabolic syndrome (4749). This suggests that MASLD’s cardiovascular impact is predominantly mediated by systemic metabolic and fibrotic abnormalities rather than direct hepatic effects. Gender-specific differences also warrant attention. In our study, liver fibrosis (as estimated by the FIB-4 index) was independently associated with carotid plaque burden in men but not in women. This finding is consistent with the recognized sexual dimorphism in cardiovascular risk, whereby premenopausal women are relatively protected from atherosclerosis due to estrogen’s vasoprotective and anti-inflammatory effects. However, postmenopausal hormonal changes, together with an increase in visceral fat distribution and loss of muscle mass, tend to align women’s cardiovascular risk with that of men over time (5052). The absence of a significant association in women in our study may reflect both a lower overall fibrosis burden and the complex interplay between hormonal status, metabolic factors, and vascular remodeling. Interestingly, in women, prediabetes and the presence of obesity remained significant independent predictors of carotid plaque, all of which are closely related to overall and central adiposity. Future studies including hormonal data or menopausal status could help clarify these sex-specific differences. Taken together, these findings underscore the need for sex-specific approaches when using liver fibrosis indices for cardiovascular risk stratification.

Our study has several limitations. First, its cross-sectional design does not allow for causal inference. Second, liver fibrosis was assessed using non-invasive indices rather than biopsy, which remains the diagnostic gold standard. Other more accurate non-invasive alternatives, such as transient elastography or the Enhanced Liver Fibrosis (ELF) score were also not used. Third, potential confounding by age should be acknowledged, given that age is a component of FIB-4 formula. Fourth, the number of participants with a high probability of advanced liver fibrosis based on FIB-4 scores was relatively small, which precluded treating this group as a separate analytical category. Fifth, atherosclerosis evaluation was confined to carotid and femoral territories, excluding other regions such as coronary arteries. Finally, potential survivor or selection bias cannot be entirely ruled out. Nonetheless, the study’s strength lies in its focus on a large, well-characterized cohort of individuals at low-to-moderate cardiovascular risk, without previous liver disease, providing valuable insights into the relationship between liver fibrosis and vascular health in the primary prevention setting. A central illustration summarizing the conceptual framework and main findings of the study is provided in Figure 2.

Figure 2
Flowchart illustrating the progression from metabolic dysfunction to vascular injury, highlighting FIB-4 as a predictor of fibrosis and carotid plaque presence. The graphic mentions higher prevalence and total plaque area associated with FIB-4. There's a significant association in men but not in women. Femoral atherosclerosis has no predictive value. The FIB-4 risk stratification uses a traffic light system: low, intermediate/high risk. Clinical implications support identifying increased plaque burden in low-to-moderate CV risk populations and emphasize integrated assessments and territory-specific vulnerability.

Figure 2. Central illustration: Association between FIB-4 and vascular atherosclerotic burden across carotid and femoral territories.

In summary, this study highlights the site-specific relationship between liver fibrosis and atherosclerosis, with advanced fibrosis preferentially associated with carotid plaque formation. These findings underscore the utility of non-invasive indices of liver fibrosis, particularly FIB-4, as a potential tool for identifying individuals at higher risk of subclinical atherosclerosis, even in populations without known liver disease. Future research should prioritize longitudinal studies to elucidate the temporal and causal links between liver fibrosis and vascular disease progression. Additionally, exploring the underlying mechanisms driving site-specific atherogenesis will be crucial for developing tailored therapeutic interventions.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the University Hospital Arnau de Vilanova. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

JL-M: Writing – original draft, Writing – review & editing. MM: Writing – original draft, Writing – review & editing. MB-L: Writing – original draft, Writing – review & editing. JV: Writing – original draft, Writing – review & editing. RP: Writing – original draft, Writing – review & editing. GT: Writing – original draft, Writing – review & editing. DM: Writing – original draft, Writing – review & editing. EC: Writing – original draft, Writing – review & editing. EF: Writing – original draft, Writing – review & editing. AC: Writing – original draft, Writing – review & editing. MH: Writing – original draft, Writing – review & editing. CL-C: Writing – original draft, Writing – review & editing. AG: Writing – original draft, Writing – review & editing. RG-R: Writing – original draft, Writing – review & editing. KC: Writing – original draft, Writing – review & editing. AL: Writing – original draft, Writing – review & editing.

Group member of ILERVAS project collaborators

Eva Miquel, Marta Ortega, Ferran Barbé, Jordi de Batlle, Silvia Barril, Manuel Portero-Otín, Mariona Jové, Josep Franch-Nadal, Esmeralda Castelblanco, Pere Godoy, Montse Martinez- Alonso, and Cristina Farràs.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research was supported by grant PI21/00462 (Institute of Health Carlos III, European Regional Development Fund, “A way to build Europe). Authors in the study also participate in this other grants: PI21/01099, PMP21/00109, PMP22/00073, PI23/00237 and RD24/0004/0015 (Instituto de Salud Carlos III, European Regional Development Fund, “A way to build Europe), PID2022-141964OB-I00 (Ministerio de Ciencia, Innovación y Universidades, co-funded by the European Regional Development Fund, “A way to build Europe”).

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor CN declared a past co-authorship with the author DM.

Generative AI statement

The author(s) declared that generative AI was used in the creation of this manuscript. Generative AI tools (ChatGPT, OpenAI) were used exclusively to refine the grammar and language of the manuscript. No part of the scientific content, data analysis, or interpretation was generated by AI. All authors have carefully reviewed the text and take full responsibility for its content.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1651689/full#supplementary-material

References

1. Simon TG, Roelstraete B, Hagström H, Sundström J, and Ludvigsson JF. Non-alcoholic fatty liver disease and incident major adverse cardiovascular events: results from a nationwide histology cohort. Gut. (2022) 71:1867–75. doi: 10.1136/gutjnl-2021-325724

PubMed Abstract | Crossref Full Text | Google Scholar

2. Rinella ME, Lazarus JV, Ratziu V, Francque SM, Sanyal AJ, Kanwal F, et al. A multisociety Delphi consensus statement on new fatty liver disease nomenclature. J Hepatol. (2023) 79:1542–56. doi: 10.1016/j.jhep.2023.06.003

PubMed Abstract | Crossref Full Text | Google Scholar

3. Diehl AM and Day C. Cause, pathogenesis, and treatment of nonalcoholic steatohepatitis. N Engl J Med. (2017) 377:2063–72. doi: 10.1056/NEJMra1503519

PubMed Abstract | Crossref Full Text | Google Scholar

4. Mundi MS, Velapati S, Patel J, Kellogg TA, Abu Dayyeh BK, and Hurt RT. Evolution of NAFLD and its management. Nutr Clin Pract. (2020) 35:72–84. doi: 10.1002/ncp.10449

PubMed Abstract | Crossref Full Text | Google Scholar

5. Byrne CD and Targher G. NAFLD: a multisystem disease. J Hepatol. (2015) 62:S47–64. doi: 10.1016/j.jhep.2014.12.012

PubMed Abstract | Crossref Full Text | Google Scholar

6. Targher G, Byrne CD, Lonardo A, Zoppini G, and Barbui C. Non-alcoholic fatty liver disease and risk of incident cardiovascular disease: A meta-analysis. J Hepatol. (2016) 65:589–600. doi: 10.1016/j.jhep.2016.05.013

PubMed Abstract | Crossref Full Text | Google Scholar

7. Abdallah LR, de Matos RC, Souza YPDM, Vieira-Soares D, Muller-MaChado G, and Pollo-Flores P. Non-alcoholic fatty liver disease and its links with inflammation and atherosclerosis. Curr Atheroscler Rep. (2020) 22:7. doi: 10.1007/s11883-020-0820-8

PubMed Abstract | Crossref Full Text | Google Scholar

8. Josloff K, Beiriger J, Khan A, Gawel RJ, Kirby RS, Kendrick AD, et al. Comprehensive review of cardiovascular disease risk in nonalcoholic fatty liver disease. J Cardiovasc Dev Dis. (2022) 26:419. doi: 10.3390/jcdd9120419

PubMed Abstract | Crossref Full Text | Google Scholar

9. Santos RD, Valenti L, and Romeo S. Does nonalcoholic fatty liver disease cause cardiovascular disease? Current knowledge and gaps. Atherosclerosis. (2019) 282:110–20. doi: 10.1016/j.atherosclerosis.2019.01.029

PubMed Abstract | Crossref Full Text | Google Scholar

10. Buckley AJ, Louise E, Lessan N, Trovato FM, Trovato GM, and Taylor-Robinson SD. Nonalcoholic fatty liver disease: relationship with cardiovascular risk markers and clinical endpoints. Diabetes Res Clin Pract. (2018) 144:144–52. doi: 10.1016/j.diabres.2018.08.011

PubMed Abstract | Crossref Full Text | Google Scholar

11. Subbotin VM. Excessive intimal hyperplasia in human coronary arteries before intimal lipid depositions is the initiation of coronary atherosclerosis and constitutes a therapeutic target. Drug Discov Today. (2016) 21:1578–95. doi: 10.1016/j.drudis.2016.05.017

PubMed Abstract | Crossref Full Text | Google Scholar

12. León-Mengíbar J, Sánchez E, Herrerías F, de la Fuente MC, Santamaría M, Valdivielso JM, et al. Influence of nonalcoholic fatty liver disease severity on carotid adventitial vasa vasorum. Front Endocrinol (Lausanne). (2024) 15:1366015. doi: 10.3389/fendo.2024.1366015

PubMed Abstract | Crossref Full Text | Google Scholar

13. Pais R, Redheuil A, Cluzel P, Ratziu V, and Giral P. Relationship among fatty liver, specific and multiple-site atherosclerosis, and 10-year Framingham score. Hepatology. (2019) 69:1453–63. doi: 10.1002/hep.30223

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ampuero J and Romero-Gómez M. Influencia de la enfermedad por hígado graso no alcohólico en la enfermedad cardiovascular. Gastroenterol Hepatol. (2012) 35:585–93. doi: 10.1016/j.gastrohep.2012.02.005

PubMed Abstract | Crossref Full Text | Google Scholar

15. Cordero A, Escribano D, Quintanilla MA, López-Ayala JM, Masiá MD, Cazorla D, et al. Prognostic value of liver fibrosis assessed by the FIB-4 index in patients with acute coronary syndrome. Rev Esp Cardiol (Engl Ed). (2023) 76:618–25. doi: 10.1016/j.rec.2022.12.013

PubMed Abstract | Crossref Full Text | Google Scholar

16. Sinn DH, Kang D, Chang Y, Ryu S, Gu S, Kim H, et al. Nonalcoholic fatty liver disease and progression of coronary artery calcium score: a retrospective cohort study. Gut. (2017) 66:323–9. doi: 10.1136/gutjnl-2016-311854

PubMed Abstract | Crossref Full Text | Google Scholar

17. Sinn DH, Cho SJ, Gu S, Seong D, Kang D, Kim H, et al. Persistent nonalcoholic fatty liver disease increases risk for carotid atherosclerosis. Gastroenterology. (2016) 151:481–8. doi: 10.1053/j.gastro.2016.06.001

PubMed Abstract | Crossref Full Text | Google Scholar

18. Pais R, Giral P, Khan JF, Rosenbaum D, Housset C, Poynard T, et al. Fatty liver is an independent predictor of early carotid atherosclerosis. J Hepatol. (2016) 65:95–102. doi: 10.1016/j.jhep.2016.02.023

PubMed Abstract | Crossref Full Text | Google Scholar

19. Siddiqui MS, Yamada G, Vuppalanchi R, Van Natta M, Loomba R, Guy C, et al. Diagnostic accuracy of noninvasive fibrosis models to detect change in fibrosis stage. Clin Gastroenterol Hepatol. (2019) 17:1877–85.e5. doi: 10.1016/j.cgh.2018.12.031

PubMed Abstract | Crossref Full Text | Google Scholar

20. Shah AG, Lydecker A, Murray K, Tetri BN, Contos MJ, Sanyal AJ, et al. Comparison of noninvasive markers of fibrosis in patients with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. (2009) 7:1104–12. doi: 10.1016/j.cgh.2009.05.033

PubMed Abstract | Crossref Full Text | Google Scholar

21. Siddiqui MS, Patidar KR, Boyett S, Luketic VA, Puri P, and Sanyal AJ. Performance of noninvasive models of fibrosis in predicting mild to moderate fibrosis in patients with nonalcoholic fatty liver disease (NAFLD). Liver Int. (2015) 36:572–9. doi: 10.1111/liv.13054

PubMed Abstract | Crossref Full Text | Google Scholar

22. Betriu À, Farràs C, Abajo M, Martinez-Alonso M, Arroyo D, Barbé F, et al. Randomised intervention study to assess the prevalence of subclinical vascular disease and hidden kidney disease and its impact on morbidity and mortality: the ILERVAS project. Nefrologia. (2016) 36:389–96. doi: 10.1016/j.nefro.2016.02.008

PubMed Abstract | Crossref Full Text | Google Scholar

23. Bermúdez-López M, Martínez-Alonso M, Castro-Boqué E, Betriu À, Cambray S, Farràs C, et al. Subclinical atheromatosis localization and burden in a low-to-moderate cardiovascular risk population: the ILERVAS study. Rev Esp Cardiol (Engl Ed). (2021) 74:1042–53. doi: 10.1016/j.rec.2020.09.015

PubMed Abstract | Crossref Full Text | Google Scholar

24. American Diabetes Association Professional Practice Committee. Diagnosis and classification of diabetes: standards of care in diabetes-2025. Diabetes Care. (2025) 48:S27–49. doi: 10.2337/dc25-S002

PubMed Abstract | Crossref Full Text | Google Scholar

25. Sabetai MM, Tegos TJ, Nicolaides AN, Dhanjil S, Pare GJ, and Stevens JM. Reproducibility of computer-quantified carotid plaque echogenicity: can we overcome the subjectivity? Stroke. (2000) 31:2189–96. doi: 10.1161/01.str.31.9.2189

PubMed Abstract | Crossref Full Text | Google Scholar

26. Touboul PJ, Hennerici MG, Meairs S, Adams H, Amarenco P, Desvarieux M, et al. Mannheim intima-media thickness consensus. Cerebrovasc Dis. (2004) 18:346–9. doi: 10.1159/000081812

PubMed Abstract | Crossref Full Text | Google Scholar

27. Stein JH, Korcarz CE, Hurst RT, Lonn E, Kendall CB, Mohler ER, et al. Use of carotid ultrasound to identify subclinical vascular disease and evaluate cardiovascular disease risk: a consensus statement from the American Society of Echocardiography Carotid Intima-Media Thickness Task Force endorsed by the Society for Vascular Medicine. J Am Soc Echocardiogr. (2008) 21:93–111. doi: 10.1016/j.echo.2007.11.011

PubMed Abstract | Crossref Full Text | Google Scholar

28. Tiozzo E, Gardener H, Hudson BI, Dong C, Della-Morte D, Crisby M, et al. High-density lipoprotein subfractions and carotid plaque: the Northern Manhattan Study. Atherosclerosis. (2014) 237:163–8. doi: 10.1016/j.atherosclerosis.2014.09.002

PubMed Abstract | Crossref Full Text | Google Scholar

29. European Association for the Study of the Liver (EASL), European Association for the Study of Diabetes (EASD), and European Association for the Study of Obesity (EASO). EASL-EASD-EASO Clinical Practice Guidelines on the management of metabolic dysfunction-associated steatotic liver disease (MASLD). J Hepatol. (2024) 81:492–542. doi: 10.1016/j.jhep.2024.04.031

PubMed Abstract | Crossref Full Text | Google Scholar

30. Cusi K, Abdelmalek MF, Apovian CM, Balapattabi K, Bannuru RR, Barb D, et al. Metabolic dysfunction-associated steatotic liver disease (MASLD) in people with diabetes: the need for screening and early intervention. A Consensus Report of the American Diabetes Association. Diabetes Care. (2025) 48:1057–82. doi: 10.2337/dci24-0094

PubMed Abstract | Crossref Full Text | Google Scholar

31. Kim D, Choi SY, Park EH, Lee W, Kang JH, Kim W, et al. Nonalcoholic fatty liver disease is associated with coronary artery calcification. Hepatology. (2012) 56:605–13. doi: 10.1002/hep.25593

PubMed Abstract | Crossref Full Text | Google Scholar

32. Gastaldelli A, Kozakova M, Hojlund K, Flyvbjerg A, Favuzzi A, Mitrakou A, et al. Fatty liver is associated with insulin resistance, risk of coronary heart disease, and early atherosclerosis in a large European population. Hepatology. (2009) 49:1537–44. doi: 10.1002/hep.22845

PubMed Abstract | Crossref Full Text | Google Scholar

33. Kozakova M, Palombo C, Eng MP, Dekker J, Flyvbjerg A, Mitrakou A, et al. Fatty liver index, gamma-glutamyltransferase, and early carotid plaques. Hepatology. (2012) 55:1406–15. doi: 10.1002/hep.25555

PubMed Abstract | Crossref Full Text | Google Scholar

34. Yi M, Peng W, Teng F, Kong Q, and Chen Z. The role of noninvasive scoring systems for predicting cardiovascular disease risk in patients with nonalcoholic fatty liver disease: a systematic review and meta-analysis. Eur J Gastroenterol Hepatol. (2022) 34:1277–84. doi: 10.1097/MEG.0000000000002462

PubMed Abstract | Crossref Full Text | Google Scholar

35. Yan Z, Liu Y, Li W, Zhao X, Lin W, Zhang J, et al. Liver fibrosis scores and prognosis in patients with cardiovascular diseases: A systematic review and meta-analysis. Eur J Clin Invest. (2022) 52:e13855. doi: 10.1111/eci.13855

PubMed Abstract | Crossref Full Text | Google Scholar

36. Vieira Barbosa J, Milligan S, Frick A, Broestl J, Younossi Z, Afdhal N, et al. Fibrosis-4 index can independently predict major adverse cardiovascular events in nonalcoholic fatty liver disease. Am J Gastroenterol. (2022) 117:453–61. doi: 10.14309/ajg.0000000000001606

PubMed Abstract | Crossref Full Text | Google Scholar

37. Cusi K and Godinez-Leiva E. Cardiovascular risk in patients with nonalcoholic fatty liver disease: looking at the liver to shield the heart. Curr Opin Lipidol. (2020) 31:364–6. doi: 10.1097/MOL.0000000000000717

PubMed Abstract | Crossref Full Text | Google Scholar

38. DeBoer MD, Lin B, Filipp S, Cusi K, and Gurka MJ. Severity of metabolic syndrome is greater among nonalcoholic adults with elevated ALT and advanced fibrosis. Nutr Res. (2021) 88:34–43. doi: 10.1016/j.nutres.2020.12.023

PubMed Abstract | Crossref Full Text | Google Scholar

39. Bril F, Sanyal A, and Cusi K. Metabolic syndrome and its association with nonalcoholic steatohepatitis. Clin Liver Dis. (2023) 27:187–210. doi: 10.1016/j.cld.2023.01.002

PubMed Abstract | Crossref Full Text | Google Scholar

40. Nasiri-Ansari N, Androutsakos T, Flessa CM, Kyrou I, Siasos G, Randeva HS, et al. Endothelial cell dysfunction and nonalcoholic fatty liver disease (NAFLD): A concise review. Cells. (2022) 11:2511. doi: 10.3390/cells11162511

PubMed Abstract | Crossref Full Text | Google Scholar

41. Du W and Wang L. The crosstalk between liver sinusoidal endothelial cells and hepatic microenvironment in NASH related liver fibrosis. Front Immunol. (2022) 13:936196. doi: 10.3389/fimmu.2022.936196

PubMed Abstract | Crossref Full Text | Google Scholar

42. Dalager S, Paaske WP, Kristensen IB, Laurberg JM, and Falk E. Artery-related differences in atherosclerosis expression: implications for atherogenesis and dynamics in intima-media thickness. Stroke. (2007) 38:2698–705. doi: 10.1161/STROKEAHA.107.486480

PubMed Abstract | Crossref Full Text | Google Scholar

43. Helck A, Bianda N, Canton G, Yuan C, Hippe DS, Reiser MF, et al. Intra-individual comparison of carotid and femoral atherosclerotic plaque features with in vivo MR plaque imaging. Int J Cardiovasc Imaging. (2015) 31:1611–8. doi: 10.1007/s10554-015-0737-4

PubMed Abstract | Crossref Full Text | Google Scholar

44. McGill HC Jr, McMahan CA, Herderick EE, Tracy RE, Malcom GT, Zieske AW, et al. Effects of coronary heart disease risk factors on atherosclerosis of selected regions of the aorta and right coronary artery. PDAY Research Group. Pathobiological Determinants of Atherosclerosis in Youth. Arterioscler Thromb Vasc Biol. (2000) 20:836–45. doi: 10.1161/01.atv.20.3.836

PubMed Abstract | Crossref Full Text | Google Scholar

45. Hong S, Sun L, Hao Y, Li P, Zhou Y, Liang X, et al. From NAFLD to MASLD: When metabolic comorbidity matters. Ann Hepatol. (2024) 29:101281. doi: 10.1016/j.aohep.2023.101281

PubMed Abstract | Crossref Full Text | Google Scholar

46. Inukai Y, Ito T, Yokoyama S, Yamamoto K, Imai N, Ishizu Y, et al. Type 2 diabetes and hypertension as risk factors for advanced fibrosis in biopsy proven metabolic dysfunction-associated steatotic liver disease. J Dig Dis. (2024) 25:685–93. doi: 10.1111/1751-2980.13325

PubMed Abstract | Crossref Full Text | Google Scholar

47. Dongiovanni P, Petta S, Maglio C, Fracanzani AL, Pipitone R, Mozzi E, et al. Transmembrane 6 superfamily member 2 gene variant disentangles nonalcoholic steatohepatitis from cardiovascular disease. Hepatology. (2015) 61:506–14. doi: 10.1002/hep.27490

PubMed Abstract | Crossref Full Text | Google Scholar

48. Simons N, Isaacs A, Koek GH, Kuč S, Schaper NC, and Brouwers MCGJ. PNPLA3, TM6SF2, and MBOAT7 genotypes and coronary artery disease. Gastroenterology. (2017) 152:912–3. doi: 10.1053/j.gastro.2016.12.020

PubMed Abstract | Crossref Full Text | Google Scholar

49. Lauridsen BK, Stender S, Kristensen TS, Kofoed KF, Køber L, Nordestgaard BG, et al. Liver fat content, non-alcoholic fatty liver disease, and ischaemic heart disease: Mendelian randomization and meta-analysis of 279–013 individuals. Eur Heart J. (2018) 39:385–93. doi: 10.1093/eurheartj/ehx662

PubMed Abstract | Crossref Full Text | Google Scholar

50. Connelly PJ, Azizi Z, Alipour P, Delles C, Pilote L, and Raparelli V. The importance of gender to understand sex differences in cardiovascular disease. Can J Cardiol. (2021) 37:699–710. doi: 10.1016/j.cjca.2021.02.005

PubMed Abstract | Crossref Full Text | Google Scholar

51. Vu H, Khanh Tuong TT, Hoang Lan N, Quoc Thang T, Bilgin K, Hoa T, et al. Correlation between nonalcoholic fatty liver disease and coronary atherosclerosis. Clin Ter. (2022) 173:565–71. doi: 10.7417/CT.2022.2483

PubMed Abstract | Crossref Full Text | Google Scholar

52. O’Neil A, Scovelle AJ, Milner AJ, and Kavanagh A. Gender/sex as a social determinant of cardiovascular risk. Circulation. (2018) 137:854–64. doi: 10.1161/CIRCULATIONAHA.117.028595

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: cardiovascular risk factors (CVRF), FIB-4 index, liver fibrosis, MASLD, metabolic dysfunction associated steatotic liver disease, NAFLD fibrosis score (NFS), subclinical atheromatous disease

Citation: León-Mengíbar J, Malagón MM, Bermúdez-López M, Valdivielso JM, Pamplona R, Torres G, Mauricio D, Castro E, Fernández E, Caixàs A, Hernández M, Lopez-Cano C, Gordon A, Guzman-Ruiz R, Cusi K and Lecube A (2026) Association of estimated liver fibrosis with carotid but not femoral atherosclerotic burden: the ILERVAS cohort. Front. Endocrinol. 16:1651689. doi: 10.3389/fendo.2025.1651689

Received: 22 June 2025; Accepted: 08 December 2025; Revised: 01 December 2025;
Published: 06 January 2026.

Edited by:

Cristiane Nogueira, Federal University of Rio de Janeiro, Brazil

Reviewed by:

Shihua Shi, Friedrich Miescher Institute for Biomedical Research (FMI), Switzerland
Shan Hong, Capital Medical University, China

Copyright © 2026 León-Mengíbar, Malagón, Bermúdez-López, Valdivielso, Pamplona, Torres, Mauricio, Castro, Fernández, Caixàs, Hernández, Lopez-Cano, Gordon, Guzman-Ruiz, Cusi and Lecube. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Albert Lecube, YWxlY3ViZUBnbWFpbC5jb20=; Kenneth Cusi, S2VubmV0aC5DdXNpQG1lZGljaW5lLnVmbC5lZHU=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.