# Identification of endophenotypes supporting outcome prediction in hemodialysis patients based on mechanistic markers of statin treatment

Johannes Leierer, Madonna Salib, Michail Evgeniou, Patrick Rossignol, Ziad A Massy, Klaus Kratochwill, Gert Mayer, Bengt Fellström, Nicolas Girerd, Faiez Zannad, et al.

► To cite this version:

Johannes Leierer, Madonna Salib, Michail Evgeniou, Patrick Rossignol, Ziad A Massy, et al. Identification of endophenotypes supporting outcome prediction in hemodialysis patients based on mechanistic markers of statin treatment. Heliyon, 2024, 10 (9), pp.e30709. 10.1016/j.heliyon.2024.e30709. hal-04582153

HAL Id: hal-04582153

https://hal.univ-lorraine.fr/hal-04582153

Submitted on 21 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Identification of endophenotypes supporting outcome prediction in hemodialysis patients based on mechanistic markers of statin treatment

Johannes Leierer, Madonna Salib, Michail Evgeniou, Patrick Rossignol, Ziad A Massy, Klaus Kratochwill, Gert Mayer, Bengt Fellström, Nicolas Girerd, Faiez Zannad, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Johannes Leierer, Madonna Salib, Michail Evgeniou, Patrick Rossignol, Ziad A Massy, et al. Identification of endophenotypes supporting outcome prediction in hemodialysis patients based on mechanistic markers of statin treatment. Heliyon, 2024, 10 (9), pp.e30709. 10.1016/j.heliyon.2024.e30709. hal-04582153

HAL Id: hal-04582153

https://hal.univ-lorraine.fr/hal-04582153v1

Submitted on 21 May 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research article

<figure><img src="image_2.png" /></figure>

Identification of endophenotypes supporting outcome prediction
in hemodialysis patients based on mechanistic markers of
statin treatment

Johannes Leierer a,<fnref n="1" />, Madonna Salib b,<fnref n="1" />, Michail Evgeniou c, Patrick Rossignol b, d, Ziad A. Massy e, Klaus Kratochwill c, Gert Mayer a, Bengt Fellström f, Nicolas Girerd b, Faiez Zannad b,<fnref n="2" />, Paul Perco a,*,<fnref n="2" />

ARTICLE INFO

Keywords:

Statins

Gene expression profiling

Predictive biomarkers

Patient stratification

Cox proportional hazards regression models

Cardiovascular risk

Hemodialysis

AURORA cohort

A B S T R A C T

**Background:** Statins are widely used to reduce the risk of cardiovascular disease (CVD). Patients with end-stage renal disease (ESRD) on hemodialysis have significantly increased risk of developing CVD. Statin treatment in these patients however did not show a statistically significant benefit in large trials on a patient cohort level.

**Methods:** We generated gene expression profiles for statins to investigate the impact on cellular programs in human renal proximal tubular cells and mesangial cells in-vitro. We subsequently selected biomarkers from key statin-affected molecular pathways and assessed these biomarkers in plasma samples from the AURORA cohort, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration that have been treated with rosuvastatin. Patient clusters (phenotypes) were created based on the identified biomarkers using Latent Class Model clustering and the associations with outcome for the generated phenotypes were assessed using Cox proportional hazards regression models. The multivariable models were adjusted for clinical and biological covariates based on previously published data in AURORA.

**Results:** The impact of statin treatment on mesangial cells was larger as compared with tubular cells with a large overlap of differentially expressed genes identified for atorvastatin and rosuvastatin indicating a predominant drug class effect. Affected molecular pathways included TGFB-, TNF-, and MAPK-signaling and focal adhesion among others. Four patient clusters were identified based on the baseline plasma concentrations of the eight biomarkers. Phenotype 1 was

<footnote>ª Medical University of Innsbruck, Department of Internal Medicine IV, Innsbruck, Austria</footnote>

<footnote>ᵇ Université de Lorraine, Inserm, Centre d'Investigations Cliniques- 1433, and Inserm U1116, CHRU Nancy, F-CRIN INI-CRCT, Nancy, France</footnote>

<footnote>Medical University of Vienna, Comprehensive Center for Pediatrics, Department of Pediatrics and Adolescent Medicine, Division of Pediatric Nephrology and Gastroenterology, Vienna, Austria</footnote>

<footnote n="5">⁵ Medical Specialties and Nephrology departments, Princess Grace Hospital, Monaco, Monaco</footnote>

<footnote>e Association pour l'Utilisation du Rein Artificiel (AURA) Paris and Department of Nephrology, CHU Ambroise Paré, APHP, 92104, Boulogne Billancourt, and Centre for Research in Epidemiology and Population Health (CESP), University Paris-Saclay, University Versailles-Saint Quentin, Inserm UMRS, 1018, Clinical Epidemiology Team, Villejuif, France</footnote>

<footnote n="1">¹ Uppsala University, Department of Medical Sciences, Uppsala, Sweden</footnote>

<footnote>* Corresponding author. Department of Internal Medicine IV (Nephrology and Hypertension) Medical University Innsbruck Anichstrasse 35 6020 Innsbruck, Austria.

E-mail address: paul.perco@i-med.ac.at (P. Perco).</footnote>

<footnote n="1">¹ J.L. and M.S. contributed equally.</footnote>

<footnote n="2">² F.Z. and P.P contributed equally.</footnote>
Research article

<figure><img src="image_3.png" /></figure>

# Identification of endophenotypes supporting outcome prediction in hemodialysis patients based on mechanistic markers of statin treatment

Johannes Leierer a,<fnref n="1" />, Madonna Salib b,<fnref n="1" />, Michail Evgeniou c, Patrick Rossignol b, d, Ziad A. Massy e, Klaus Kratochwill c, Gert Mayer a, Bengt Fellström f, Nicolas Girerd b, Faiez Zannad b,<fnref n="2" />, Paul Perco a,*,<fnref n="2" />

ARTICLE INFO

Keywords:

Statins

Gene expression profiling

Predictive biomarkers

Patient stratification

Cox proportional hazards regression models

Cardiovascular risk

Hemodialysis

AURORA cohort

A B S T R A C T

**Background:** Statins are widely used to reduce the risk of cardiovascular disease (CVD). Patients with end-stage renal disease (ESRD) on hemodialysis have significantly increased risk of developing CVD. Statin treatment in these patients however did not show a statistically significant benefit in large trials on a patient cohort level.

**Methods:** We generated gene expression profiles for statins to investigate the impact on cellular programs in human renal proximal tubular cells and mesangial cells in-vitro. We subsequently selected biomarkers from key statin-affected molecular pathways and assessed these biomarkers in plasma samples from the AURORA cohort, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration that have been treated with rosvastatin. Patient clusters (phenotypes) were created based on the identified biomarkers using Latent Class Model clustering and the associations with outcome for the generated phenotypes were assessed using Cox proportional hazards regression models. The multivariable models were adjusted for clinical and biological covariates based on previously published data in AURORA.

**Results:** The impact of statin treatment on mesangial cells was larger as compared with tubular cells with a large overlap of differentially expressed genes identified for atorvastatin and rosvastatin indicating a predominant drug class effect. Affected molecular pathways included TGFB-, TNF-, and MAPK-signaling and focal adhesion among others. Four patient clusters were identified based on the baseline plasma concentrations of the eight biomarkers. Phenotype 1 was

<footnote>ª Medical University of Innsbruck, Department of Internal Medicine IV, Innsbruck, Austria</footnote>

<footnote>ᵇ Université de Lorraine, Inserm, Centre d'Investigations Cliniques- 1433, and Inserm U1116, CHRU Nancy, F-CRIN INI-CRCT, Nancy, France</footnote>

<footnote>Medical University of Vienna, Comprehensive Center for Pediatrics, Department of Pediatrics and Adolescent Medicine, Division of Pediatric Nephrology and Gastroenterology, Vienna, Austria</footnote>

<footnote n="5">⁵ Medical Specialties and Nephrology departments, Princess Grace Hospital, Monaco, Monaco</footnote>

<footnote>e Association pour l'Utilisation du Rein Artificiel (AURA) Paris and Department of Nephrology, CHU Ambroise Paré, APHP, 92104, Boulogne Billancourt, and Centre for Research in Epidemiology and Population Health (CESP), University Paris-Saclay, University Versailles-Saint Quentin, Inserm UMRS, 1018, Clinical Epidemiology Team, Villejuif, France</footnote>

<footnote n="1">¹ Uppsala University, Department of Medical Sciences, Uppsala, Sweden</footnote>

<footnote>* Corresponding author. Department of Internal Medicine IV (Nephrology and Hypertension) Medical University Innsbruck Anichstrasse 35 6020 Innsbruck, Austria.

E-mail address: paul.perco@i-med.ac.at (P. Perco).</footnote>

<footnote n="1">¹ J.L. and M.S. contributed equally.</footnote>

<footnote n="2">² F.Z. and P.P contributed equally.</footnote>
characterized by low to medium levels of the hepatocyte growth factor (HGF) and high levels of interleukin 6 (IL6) or matrix metalloproteinase 2 (MMP2) and it was significantly associated with outcome showing increased risk of developing major adverse cardiovascular events (MACE) or cardiovascular death. Phenotype 2 had high HGF but low Fas cell surface death receptor (FAS) levels and it was associated with significantly better outcome at 1 year.

**Conclusions:** In this translational study, we identified patient subgroups based on mechanistic markers of statin therapy that are associated with disease outcome in patients on hemodialysis.

# 1. Background

Statins are among the most prescribed drugs in the world. They exert their lipid-lowering effect by targeting the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), better known as HMG-CoA reductase, one of the key enzymes in the cholesterol synthesis pathway. Statin use is recommended for adults at risk of developing cardiovascular disease (CVD) [1,2]. Variability in response to statin therapy is however substantial and there are currently no good parameters for defining statin responders. Next to this variability in response, considerable side effects seem to be associated with statin therapy such as an elevated risk of diabetes mellitus, elevated transaminase levels, myalgia, or rhabdomyolysis [3-5]. In addition, some interfere with the metabolism of other drugs making combination therapy complex.

Patients with end-stage renal disease (ESRD) on hemodialysis have a significantly increased risk of developing CVD. Large trials such as AURORA or the German 4D study investigating the impact of rosuvastatin and atorvastatin on clinical outcomes in hemodialysis patients however failed to demonstrate a beneficial impact of statin treatment in this patient population on the cohort level [6, 7].

Researchers have therefore started to search for predictive markers for statin therapy to identify responders and optimize therapeutic use of statins. These approaches included pharmacogenomic markers [8] as also reviewed by Arrigoni and colleagues [9], pharmacometabolomic markers [10], and recently also in linking changes in the human gut microbiome with response to statin therapy [11]. Associations between gene polymorphisms in HMGCR and concentrations of inflammatory serum markers that might influence the response to statin therapy regarding the anti-inflammatory potential have been recently described [12]. Other groups focused on clinical parameters in calculating risk scores to guide therapeutic intervention in hemodialysis patients [13].

In this study we investigated the effect of atorvastatin and rosuvastatin on the transcriptional level in renal mesangial and proximal tubular cells using whole genome microarrays. We analyzed changes in gene expression on the level of molecular pathways and selected a panel of biomarkers for measurements in plasma samples from the AURORA study, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration. The marker levels were used to stratify patients, and resulting phenotypes were evaluated regarding their association with disease outcome and response to statin therapy.

# 2. Methods

## 2.1. Renal cell culture models

Normal human mesangial cells (MCs) were purchased from Lonza (CC-2559, Basel, Switzerland) and cultured according to manufacturer's protocol in MsGM medium (Lonza, CC-3147). Proximal tubular human kidney cells (HK2) were purchased from American Type Culture Collection (CRL-2190, Wesel, Germany) and cultured in Keratinocyte-Serum Free Medium (KSFM) containing 10 % fetal bovine serum (FBS), 5 ng/ml recombinant epidermal growth factor (rEGF), 0.05 mg/ml bovine pituitary extract (BPE), 100 U/ml penicillin and 100 µg/ml streptomycin. All cells were grown at 37 °C in a humidified atmosphere with 5 % CO₂. After growth to confluency, cells were serum deprived and stimulated with high glucose (30 mM) for 24 h. Thereafter cells were treated with 20 µM atorvastatin or rosuvastatin for 24 h. Total RNA was isolated with RNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's protocol. RNA yield and quality was determined using a DS-11 FX + spectrophotometer (DeNovix, Wilmington, DEL, USA).

## 2.2. Transcriptomics profiling and gene expression analysis

All materials for microarray experiments were purchased from Agilent (Agilent Technologies, Inc, USA) unless otherwise specified. For each sample 200 ng total RNA was used. Cyanine 3 (Cy3) labelled cRNA was generated with the LowInput QuickAmp Labeling Kit (5190-2305) for hybridization to oligonucleotide microarrays (Human GE 4 × 44K v2, G4845A). Arrays were scanned at 5 µm resolution using an Axon Gene Pix 4000B scanner (Molecular Devices, Sunnyvale, CA) and signal intensities were extracted using Agilent Feature Extraction software (v.9.5.3.1).

Microarray preprocessing and statistical analysis was conducted with the limma package of the statistical software framework R [14]. Arrays of HK2 and MCs were preprocessed separately. Preprocessing consisted of background correction, quantile normalization as well as summarization of probes with the same sequence spotted multiple times on the array. The moderated t-statistics with adjustment for multiple testing setting the false discovery rate (FDR) to < 5 % was used to identify differentially expressed transcripts between cells treated with either atorvastatin or rosuvastatin and untreated cells respectively. The overlap between sets of identified
characterized by low to medium levels of the hepatocyte growth factor (HGF) and high levels of interleukin 6 (IL6) or matrix metalloproteinase 2 (MMP2) and it was significantly associated with outcome showing increased risk of developing major adverse cardiovascular events (MACE) or cardiovascular death. Phenotype 2 had high HGF but low Fas cell surface death receptor (FAS) levels and it was associated with significantly better outcome at 1 year.

**Conclusions:** In this translational study, we identified patient subgroups based on mechanistic markers of statin therapy that are associated with disease outcome in patients on hemodialysis.

# 1. Background

Statins are among the most prescribed drugs in the world. They exert their lipid-lowering effect by targeting the 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), better known as HMG-CoA reductase, one of the key enzymes in the cholesterol synthesis pathway. Statin use is recommended for adults at risk of developing cardiovascular disease (CVD) [1,2]. Variability in response to statin therapy is however substantial and there are currently no good parameters for defining statin responders. Next to this variability in response, considerable side effects seem to be associated with statin therapy such as an elevated risk of diabetes mellitus, elevated transaminase levels, myalgia, or rhabdomyolysis [3-5]. In addition, some interfere with the metabolism of other drugs making combination therapy complex.

Patients with end-stage renal disease (ESRD) on hemodialysis have a significantly increased risk of developing CVD. Large trials such as AURORA or the German 4D study investigating the impact of rosuvastatin and atorvastatin on clinical outcomes in hemodialysis patients however failed to demonstrate a beneficial impact of statin treatment in this patient population on the cohort level [6, 7].

Researchers have therefore started to search for predictive markers for statin therapy to identify responders and optimize therapeutic use of statins. These approaches included pharmacogenomic markers [8] as also reviewed by Arrigoni and colleagues [9], pharmacometabolomic markers [10], and recently also in linking changes in the human gut microbiome with response to statin therapy [11]. Associations between gene polymorphisms in HMGCR and concentrations of inflammatory serum markers that might influence the response to statin therapy regarding the anti-inflammatory potential have been recently described [12]. Other groups focused on clinical parameters in calculating risk scores to guide therapeutic intervention in hemodialysis patients [13].

In this study we investigated the effect of atorvastatin and rosuvastatin on the transcriptional level in renal mesangial and proximal tubular cells using whole genome microarrays. We analyzed changes in gene expression on the level of molecular pathways and selected a panel of biomarkers for measurements in plasma samples from the AURORA study, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration. The marker levels were used to stratify patients, and resulting phenotypes were evaluated regarding their association with disease outcome and response to statin therapy.

# 2. Methods

## 2.1. Renal cell culture models

Normal human mesangial cells (MCs) were purchased from Lonza (CC-2559, Basel, Switzerland) and cultured according to manufacturer's protocol in MsGM medium (Lonza, CC-3147). Proximal tubular human kidney cells (HK2) were purchased from American Type Culture Collection (CRL-2190, Wesel, Germany) and cultured in Keratinocyte-Serum Free Medium (KSFM) containing 10% fetal bovine serum (FBS), 5 ng/ml recombinant epidermal growth factor (rEGF), 0.05 mg/ml bovine pituitary extract (BPE), 100 U/ml penicillin and 100 µg/ml streptomycin. All cells were grown at 37 °C in a humidified atmosphere with 5% CO₂. After growth to confluency, cells were serum deprived and stimulated with high glucose (30 mM) for 24 h. Thereafter cells were treated with 20 µM atorvastatin or rosuvastatin for 24 h. Total RNA was isolated with RNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's protocol. RNA yield and quality was determined using a DS-11 FX + spectrophotometer (DeNovix, Wilmington, DEL, USA).

## 2.2. Transcriptomics profiling and gene expression analysis

All materials for microarray experiments were purchased from Agilent (Agilent Technologies, Inc, USA) unless otherwise specified. For each sample 200 ng total RNA was used. Cyanine 3 (Cy3) labelled cRNA was generated with the LowInput QuickAmp Labeling Kit (5190-2305) for hybridization to oligonucleotide microarrays (Human GE 4 × 44K v2, G4845A). Arrays were scanned at 5 µm resolution using an Axon Gene Pix 4000B scanner (Molecular Devices, Sunnyvale, CA) and signal intensities were extracted using Agilent Feature Extraction software (v.9.5.3.1).

Microarray preprocessing and statistical analysis was conducted with the limma package of the statistical software framework R [14]. Arrays of HK2 and MCs were preprocessed separately. Preprocessing consisted of background correction, quantile normalization as well as summarization of probes with the same sequence spotted multiple times on the array. The moderated t-statistics with adjustment for multiple testing setting the false discovery rate (FDR) to < 5% was used to identify differentially expressed transcripts between cells treated with either atorvastatin or rosuvastatin and untreated cells respectively. The overlap between sets of identified
differentially expressed transcripts was visualized with a Venn diagram. Differentially expressed transcripts were further annotated with Ensembl GeneIDs via Ensembl Biomart and the number of unique protein coding differentially expressed genes (DEGs) was determined.

Microarray data are deposited in the Gene Expression Omnibus with accession number GSE208579.

## 2.3. Pathway enrichment analysis and biomarker selection

The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 was used for functional enrichment analysis [15]. The sets of DEGs for the two drugs and the two cell lines were separately used as input for identification of enriched molecular pathways making use of the KEGG pathway database [16]. Pathways with enrichment p-values smaller than 0.1 were considered as relevant. Relevant pathways were ranked in ascending order based on the sum of the four individual p-values.

Biomarkers for measurements in the AURORA study were selected based on their regulation by statin treatment in-vitro, their involvement in mechanistic processes linked to statin treatment, their associations to CVD and CKD, as well as based on the availability of analytes on the respective OLink panels.

## 2.4. Biomarker measurements in the AURORA study

Biomarkers were measured in 398 baseline samples of the AURORA (A study to evaluate the Use of Rosuvastatin in subjects On Regular hemodialysis: an Assessment of survival and cardiovascular events) trial, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration for at least three months. After providing written informed consent, eligible patients were randomly assigned in a 1:1 ratio to the treatment arm (10 mg daily dose of rosuvastatin) or the placebo group. Description, baseline data and main results of the AURORA study have been published previously [6,17,18]. In brief, rosuvastatin did not show a significant beneficial impact in the overall cohort with respect to the composite primary endpoint of MACE nor the secondary endpoint of all-cause mortality.

The selected biomarkers were measured with the Olink Proseek Multiplex CVD panels II and III as well as the Inflammation panel in plasma samples (Olink Proteomics AB, Uppsala, Sweden). The platform is based on a proximity extension assay (PEA) technology where 92 oligonucleotide-labelled antibody probe pairs per panel are allowed to bind to their respective targets in 1 µL plasma in a 96-well plate format. When binding to their correct targets, they give rise to new DNA amplicons each ID-barcoding their respective antigens. The amplicons are subsequently quantified using a Fluidigm BioMark HD realtime polymerase chain reaction (PCR) platform. The platform provides log2 Normalized Protein eXpression (NPX) data wherein high protein values correspond to high protein concentrations, but not an absolute quantification (https://olink.com/faq/what-is-npx/).

## 2.5. Patient stratification based on molecular biomarkers in the AURORA study

We computed residuals from fitted linear regression models, where each protein biomarker was used as dependent variable and sex, age, and calculated Kt/V were used as independent variables. Latent Class Model clustering was performed on the residuals to identify patient patterns using the VarSelLCM R package.

To assess the biomarker relevance for each phenotype, we generated a decision tree to predict the phenotypes based on biomarker concentrations using the rpart R package. The minimum number of observations to be present in a node in order to perform additional splits was set to be >10 % of the total patient population.

## 2.6. Statistical outcome analysis

Comparisons of patient characteristics across the identified phenotypes were analyzed using analysis of variance (ANOVA), Kruskal-Wallis-, and chi-squared tests, as appropriate. Categorical variables are expressed as frequencies (percentages), whereas continuous variables are expressed as means ± standard deviation (SD) for normally distributed data and median along with interquartile ranges (25th and 75th percentiles) for non-normally distributed data.

The associations between phenotypes and outcome (CV death, all-cause mortality, and MACE) were assessed using Cox proportional hazards regression models. The multivariable models were adjusted for clinical and biological covariates based on previously published data in AURORA (i.e. age, history of CV disease, diabetes mellitus, serum albumin, high sensitivity C-reactive protein (hs-CRP), dialysis vintage, and B-type natriuretic peptide (BNP)) [19]. History of CV disease is a composite variable of history of coronary heart disease (i.e. prior myocardial infarction, prior coronary angioplasty or stent, and coronary artery bypass graft), history of vascular disease (i.e. peripheral artery disease, abdominal aortic aneurysm, carotid artery disease, carotid stenosis ≥ 50 %, and carotid endarterectomy), and history of neurovascular disease (i.e. prior ischemic vascular accident and transient ischemic attack). BNP and hs-CRP were best modelled by using their natural logarithm.

The Kaplan Meier method was used to estimate risk for each outcome according to phenotypes and are presented as event curves. The risk differences at one year, two years, and three years with confidence intervals (CI) at 95 % between rosuvastatin and placebo groups are provided.

All statistical analyses were performed using R version 4.0.2 software (R Development Core Team).
differentially expressed transcripts was visualized with a Venn diagram. Differentially expressed transcripts were further annotated with Ensembl GeneIDs via Ensembl Biomart and the number of unique protein coding differentially expressed genes (DEGs) was determined.

Microarray data are deposited in the Gene Expression Omnibus with accession number GSE208579.

## 2.3. Pathway enrichment analysis and biomarker selection

The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 was used for functional enrichment analysis [15]. The sets of DEGs for the two drugs and the two cell lines were separately used as input for identification of enriched molecular pathways making use of the KEGG pathway database [16]. Pathways with enrichment p-values smaller than 0.1 were considered as relevant. Relevant pathways were ranked in ascending order based on the sum of the four individual p-values.

Biomarkers for measurements in the AURORA study were selected based on their regulation by statin treatment in-vitro, their involvement in mechanistic processes linked to statin treatment, their associations to CVD and CKD, as well as based on the availability of analytes on the respective OLink panels.

## 2.4. Biomarker measurements in the AURORA study

Biomarkers were measured in 398 baseline samples of the AURORA (A study to evaluate the Use of Rosuvastatin in subjects On Regular hemodialysis: an Assessment of survival and cardiovascular events) trial, a double-blind, randomized, multi-center study of patients on hemodialysis or hemofiltration for at least three months. After providing written informed consent, eligible patients were randomly assigned in a 1:1 ratio to the treatment arm (10 mg daily dose of rosuvastatin) or the placebo group. Description, baseline data and main results of the AURORA study have been published previously [6,17,18]. In brief, rosuvastatin did not show a significant beneficial impact in the overall cohort with respect to the composite primary endpoint of MACE nor the secondary endpoint of all-cause mortality.

The selected biomarkers were measured with the Olink Proseek Multiplex CVD panels II and III as well as the Inflammation panel in plasma samples (Olink Proteomics AB, Uppsala, Sweden). The platform is based on a proximity extension assay (PEA) technology where 92 oligonucleotide-labelled antibody probe pairs per panel are allowed to bind to their respective targets in 1 µL plasma in a 96-well plate format. When binding to their correct targets, they give rise to new DNA amplicons each ID-barcoding their respective antigens. The amplicons are subsequently quantified using a Fluidigm BioMark HD realtime polymerase chain reaction (PCR) platform. The platform provides log2 Normalized Protein eXpression (NPX) data wherein high protein values correspond to high protein concentrations, but not an absolute quantification (https://olink.com/faq/what-is-npx/).

## 2.5. Patient stratification based on molecular biomarkers in the AURORA study

We computed residuals from fitted linear regression models, where each protein biomarker was used as dependent variable and sex, age, and calculated Kt/V were used as independent variables. Latent Class Model clustering was performed on the residuals to identify patient patterns using the VarSelLCM R package.

To assess the biomarker relevance for each phenotype, we generated a decision tree to predict the phenotypes based on biomarker concentrations using the rpart R package. The minimum number of observations to be present in a node in order to perform additional splits was set to be >10 % of the total patient population.

## 2.6. Statistical outcome analysis

Comparisons of patient characteristics across the identified phenotypes were analyzed using analysis of variance (ANOVA), Kruskal-Wallis-, and chi-squared tests, as appropriate. Categorical variables are expressed as frequencies (percentages), whereas continuous variables are expressed as means ± standard deviation (SD) for normally distributed data and median along with interquartile ranges (25th and 75th percentiles) for non-normally distributed data.

The associations between phenotypes and outcome (CV death, all-cause mortality, and MACE) were assessed using Cox proportional hazards regression models. The multivariable models were adjusted for clinical and biological covariates based on previously published data in AURORA (i.e. age, history of CV disease, diabetes mellitus, serum albumin, high sensitivity C-reactive protein (hs-CRP), dialysis vintage, and B-type natriuretic peptide (BNP)) [19]. History of CV disease is a composite variable of history of coronary heart disease (i.e. prior myocardial infarction, prior coronary angioplasty or stent, and coronary artery bypass graft), history of vascular disease (i.e. peripheral artery disease, abdominal aortic aneurysm, carotid artery disease, carotid stenosis ≥ 50 %, and carotid endarterectomy), and history of neurovascular disease (i.e. prior ischemic vascular accident and transient ischemic attack). BNP and hs-CRP were best modelled by using their natural logarithm.

The Kaplan Meier method was used to estimate risk for each outcome according to phenotypes and are presented as event curves. The risk differences at one year, two years, and three years with confidence intervals (CI) at 95 % between rosuvastatin and placebo groups are provided.

All statistical analyses were performed using R version 4.0.2 software (R Development Core Team).
### 3. Results

Statin treatment significantly impacts gene expression in proximal tubular cells and mesangial cells in-vitro.

34,127 unique transcripts remained after background correction, quantile normalization and probe summarization for statistical analysis. Rosuvastatin treatment significantly affected 188 transcripts in HK2 cells and 340 transcripts in MCs which could be mapped to 167 and 304 unique protein coding genes respectively (see Fig. 1A). Atorvastatin treatment had a larger impact on gene expression changes than rosuvastatin, significantly affecting gene expression of 423 transcripts in HK2 cells and 646 transcripts in MCs. These transcripts could be mapped to 388 and 570 unique protein coding genes respectively.

Significant overlaps between DEGs affected by both drugs could be detected in HK2 cells (p-value <0.0001) as well as in MCs (p-value <0.0001) based on chi-squared tests. This set of molecules affected by both drugs represents the drug class effect, and of note that not a single feature in this set was upregulated by one drug and downregulated by the other drug thus further strengthening the quality of the generated gene expression dataset.

The overall effect on gene expression changes was larger in MCs as compared to HK2 cells with in total 684 differentially expressed transcripts in MCs as compared to 462 differentially expressed transcripts in HK2 cells.

A set of 17 genes was affected by both drugs in both cell lines. 13 were upregulated after statin treatment including for example the pro-fibrotic connective tissue growth factor (CTGF). Three genes were downregulated after statin treatment including for example the Kruppel like factor 2 (KLF2). One gene, namely the zinc finger BED-type containing 2 (ZBED2), was upregulated in MCs but down-regulated in HK2 cells.

Statin treatment affects gene expression of its direct target an associated genes linked to cholesterol biosynthesis.

We observed a significant upregulation of HMGCR on the mRNA level in MCs but not in HK2 cells (see Fig. 1B). In addition to HMGCR, also the receptor for low density lipoprotein (LDLR) was upregulated in MCs but not in HK2 cells in our study. The proprotein convertase subtilisin/kexin type 9 (PCSK9) and squalene epoxidase (SQLE), one of the rate-limiting enzymes in the sterol biosynthesis pathway was also significantly upregulated in MCs after statin treatment in our experiments.

Next to the direct target of statins and additional members of the cholesterol biosynthesis we focused on selected key processes being enriched in DEGs in order to investigate the pleiotropic effects of statins on renal cells.

Statin treatment significantly impacts molecular pathways and molecular biomarkers linked to CKD and CVD progression.

29 molecular KEGG pathways were identified as being enriched based on at least one of the four DEG lists.

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Statin treatment effect on gene expression in renal cells</figcaption></figure>

*Panel A:* The number of differentially expressed transcripts between treated and untreated cells for the two cell lines studied is depicted in the Venn diagram. A set of 17 transcripts was deregulated in both cell lines by both drugs studied.

*Panel B:* Graphical representation of log fold-change (logFC) values of the direct statin drug target HMGCR and associated genes LDLR, PCSK9, and SQLE. Significant fold-change values are indicated with ** (FDR <5 %) and * (FDR <10 %) respectively.
### 3. Results

Statin treatment significantly impacts gene expression in proximal tubular cells and mesangial cells in-vitro.

34,127 unique transcripts remained after background correction, quantile normalization and probe summarization for statistical analysis. Rosuvastatin treatment significantly affected 188 transcripts in HK2 cells and 340 transcripts in MCs which could be mapped to 167 and 304 unique protein coding genes respectively (see Fig. 1A). Atorvastatin treatment had a larger impact on gene expression changes than rosuvastatin, significantly affecting gene expression of 423 transcripts in HK2 cells and 646 transcripts in MCs. These transcripts could be mapped to 388 and 570 unique protein coding genes respectively.

Significant overlaps between DEGs affected by both drugs could be detected in HK2 cells (p-value <0.0001) as well as in MCs (p-value <0.0001) based on chi-squared tests. This set of molecules affected by both drugs represents the drug class effect, and of note that not a single feature in this set was upregulated by one drug and downregulated by the other drug thus further strengthening the quality of the generated gene expression dataset.

The overall effect on gene expression changes was larger in MCs as compared to HK2 cells with in total 684 differentially expressed transcripts in MCs as compared to 462 differentially expressed transcripts in HK2 cells.

A set of 17 genes was affected by both drugs in both cell lines. 13 were upregulated after statin treatment including for example the pro-fibrotic connective tissue growth factor (CTGF). Three genes were downregulated after statin treatment including for example the Kruppel like factor 2 (KLF2). One gene, namely the zinc finger BED-type containing 2 (ZBED2), was upregulated in MCs but down-regulated in HK2 cells.

Statin treatment affects gene expression of its direct target an associated genes linked to cholesterol biosynthesis.

We observed a significant upregulation of HMGCR on the mRNA level in MCs but not in HK2 cells (see Fig. 1B). In addition to HMGCR, also the receptor for low density lipoprotein (LDLR) was upregulated in MCs but not in HK2 cells in our study. The proprotein convertase subtilisin/kexin type 9 (PCSK9) and squalene epoxidase (SQLE), one of the rate-limiting enzymes in the sterol biosynthesis pathway was also significantly upregulated in MCs after statin treatment in our experiments.

Next to the direct target of statins and additional members of the cholesterol biosynthesis we focused on selected key processes being enriched in DEGs in order to investigate the pleiotropic effects of statins on renal cells.

Statin treatment significantly impacts molecular pathways and molecular biomarkers linked to CKD and CVD progression.

29 molecular KEGG pathways were identified as being enriched based on at least one of the four DEG lists.

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Statin treatment effect on gene expression in renal cells</figcaption></figure>

*Panel A:* The number of differentially expressed transcripts between treated and untreated cells for the two cell lines studied is depicted in the Venn diagram. A set of 17 transcripts was deregulated in both cell lines by both drugs studied.

*Panel B:* Graphical representation of log fold-change (logFC) values of the direct statin drug target HMGCR and associated genes LDLR, PCSK9, and SQLE. Significant fold-change values are indicated with ** (FDR <5 %) and * (FDR <10 %) respectively.
The top ranked enriched molecular KEGG pathway based on the sum of p-values of the four comparisons was the MAPK signaling pathway. Among the top enriched pathways were also the fibrosis related TGF-beta signaling pathway, the inflammation related pathways TNF-signaling and leukocyte transendothelial migration but also pathways on focal adhesion and mineral absorption (see Table 1).

Eleven genes linked to the fibrosis related TGF-beta signaling pathway were identified as differentially expressed after statin treatment, mainly in MCs. The set of downregulated genes contained the two TGF-beta isoforms 2 (TGFB2) and 3 (TGFB3) as well as the growth differentiation factor 6 (GDF6), bone morphogenetic protein 4 (BMP4) as well as follistatin (FST). One of the endogenous inhibitors of the TGF-beta signaling pathway, SMAD6, on the other hand, was induced by statin treatment in-vitro. The impact of atorvastatin and rosuvastatin on gene expression of TGF-beta associated transcripts was very consistent indicating a drug class effect on this molecular process.

A clear difference in the regulation of MAPK components is apparent in the two cell types under study. Whereas a large number of molecules from the MAPK signaling pathway is upregulated in HK2 cells with almost none of the MAPK signaling members being downregulated, there are twice as many downregulated than upregulated MAPK signaling molecules in MCs.

Factors being associated with leukocyte transendothelial migration that were upregulated by statins in our cell lines included

<table><caption>**Table 1**

Affected molecular pathways after statin therapy.</caption><thead><tr><th rowspan="2">Pathway</th><th colspan="2">HK2</th><th colspan="2">MC</th></tr><tr><th>ATO</th><th>ROSU</th><th>ATO</th><th>ROSU</th></tr></thead><tbody><tr><td>MAPK signaling pathway</td><td>0.06451</td><td>-</td><td>0.00008</td><td>0.00057</td></tr><tr><td>Focal adhesion</td><td>0.04461</td><td>-</td><td>0.05932</td><td>0.03017</td></tr><tr><td>Leukocyte transendothelial migration</td><td>-</td><td>-</td><td>0.00280</td><td>0.00114</td></tr><tr><td>Mineral absorption</td><td>0.00367</td><td>-</td><td>0.00101</td><td>-</td></tr><tr><td>TNF signaling pathway</td><td>-</td><td>-</td><td>0.01199</td><td>0.00056</td></tr><tr><td>TGF-beta signaling pathway</td><td>-</td><td>-</td><td>0.00907</td><td>0.00369</td></tr><tr><td>Renin secretion</td><td>-</td><td>-</td><td>0.00687</td><td>0.02676</td></tr><tr><td>Biosynthesis of antibiotics</td><td>-</td><td>-</td><td>0.00123</td><td>0.03493</td></tr><tr><td>cGMP-PKG signaling pathway</td><td>0.01308</td><td>-</td><td>0.03178</td><td>-</td></tr><tr><td>Estrogen signaling pathway</td><td>-</td><td>-</td><td>0.02279</td><td>0.03127</td></tr><tr><td>Hippo signaling pathway</td><td>-</td><td>-</td><td>0.00674</td><td>0.05210</td></tr><tr><td>Vascular smooth muscle contraction</td><td>-</td><td>-</td><td>0.02390</td><td>0.06060</td></tr><tr><td>Gastric acid secretion</td><td>0.08381</td><td>0.00504</td><td>-</td><td>-</td></tr><tr><td>Oxytocin signaling pathway</td><td>0.02724</td><td>0.06369</td><td>-</td><td>-</td></tr><tr><td>Neurotrophin signaling pathway</td><td>-</td><td>-</td><td>0.06021</td><td>0.06238</td></tr><tr><td>Calcium signaling pathway</td><td>0.05076</td><td>0.09113</td><td>-</td><td>-</td></tr><tr><td>Inositol phosphate metabolism</td><td>0.00107</td><td>-</td><td>-</td><td>-</td></tr><tr><td>FoxO signaling pathway</td><td>0.00333</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Phosphatidylinositol signaling system</td><td>0.02371</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Pentose and glucuronate interconversions</td><td>-</td><td>-</td><td>-</td><td>0.02564</td></tr><tr><td>Dorso-ventral axis formation</td><td>-</td><td>0.02705</td><td>-</td><td>-</td></tr><tr><td>Cysteine and methionine metabolism</td><td>-</td><td>-</td><td>-</td><td>0.02953</td></tr><tr><td>Inflammatory mediator regulation of TRP channels</td><td>-</td><td>-</td><td>-</td><td>0.03010</td></tr><tr><td>Pyruvate metabolism</td><td>-</td><td>-</td><td>-</td><td>0.03372</td></tr><tr><td>Adherens junction</td><td>-</td><td>-</td><td>0.03801</td><td>-</td></tr><tr><td>Fc gamma R-mediated phagocytosis</td><td>-</td><td>-</td><td>0.07454</td><td>-</td></tr><tr><td>Regulation of actin cytoskeleton</td><td>-</td><td>-</td><td>-</td><td>0.08180</td></tr><tr><td>Metabolic pathways</td><td>0.08439</td><td>-</td><td>-</td><td>-</td></tr><tr><td>PPAR signaling pathway</td><td>-</td><td>-</td><td>0.08586</td><td>-</td></tr></tbody></table>
The top ranked enriched molecular KEGG pathway based on the sum of p-values of the four comparisons was the MAPK signaling pathway. Among the top enriched pathways were also the fibrosis related TGF-beta signaling pathway, the inflammation related pathways TNF-signaling and leukocyte transendothelial migration but also pathways on focal adhesion and mineral absorption (see Table 1).

Eleven genes linked to the fibrosis related TGF-beta signaling pathway were identified as differentially expressed after statin treatment, mainly in MCs. The set of downregulated genes contained the two TGF-beta isoforms 2 (TGFB2) and 3 (TGFB3) as well as the growth differentiation factor 6 (GDF6), bone morphogenetic protein 4 (BMP4) as well as follistatin (FST). One of the endogenous inhibitors of the TGF-beta signaling pathway, SMAD6, on the other hand, was induced by statin treatment in-vitro. The impact of atorvastatin and rosuvastatin on gene expression of TGF-beta associated transcripts was very consistent indicating a drug class effect on this molecular process.

A clear difference in the regulation of MAPK components is apparent in the two cell types under study. Whereas a large number of molecules from the MAPK signaling pathway is upregulated in HK2 cells with almost none of the MAPK signaling members being downregulated, there are twice as many downregulated than upregulated MAPK signaling molecules in MCs.

Factors being associated with leukocyte transendothelial migration that were upregulated by statins in our cell lines included

<table><caption>**Table 1**

Affected molecular pathways after statin therapy.</caption><thead><tr><th rowspan="2">Pathway</th><th colspan="2">HK2</th><th colspan="2">MC</th></tr><tr><th>ATO</th><th>ROSU</th><th>ATO</th><th>ROSU</th></tr></thead><tbody><tr><td>MAPK signaling pathway</td><td>0.06451</td><td>-</td><td>0.00008</td><td>0.00057</td></tr><tr><td>Focal adhesion</td><td>0.04461</td><td>-</td><td>0.05932</td><td>0.03017</td></tr><tr><td>Leukocyte transendothelial migration</td><td>-</td><td>-</td><td>0.00280</td><td>0.00114</td></tr><tr><td>Mineral absorption</td><td>0.00367</td><td>-</td><td>0.00101</td><td>-</td></tr><tr><td>TNF signaling pathway</td><td>-</td><td>-</td><td>0.01199</td><td>0.00056</td></tr><tr><td>TGF-beta signaling pathway</td><td>-</td><td>-</td><td>0.00907</td><td>0.00369</td></tr><tr><td>Renin secretion</td><td>-</td><td>-</td><td>0.00687</td><td>0.02676</td></tr><tr><td>Biosynthesis of antibiotics</td><td>-</td><td>-</td><td>0.00123</td><td>0.03493</td></tr><tr><td>cGMP-PKG signaling pathway</td><td>0.01308</td><td>-</td><td>0.03178</td><td>-</td></tr><tr><td>Estrogen signaling pathway</td><td>-</td><td>-</td><td>0.02279</td><td>0.03127</td></tr><tr><td>Hippo signaling pathway</td><td>-</td><td>-</td><td>0.00674</td><td>0.05210</td></tr><tr><td>Vascular smooth muscle contraction</td><td>-</td><td>-</td><td>0.02390</td><td>0.06060</td></tr><tr><td>Gastric acid secretion</td><td>0.08381</td><td>0.00504</td><td>-</td><td>-</td></tr><tr><td>Oxytocin signaling pathway</td><td>0.02724</td><td>0.06369</td><td>-</td><td>-</td></tr><tr><td>Neurotrophin signaling pathway</td><td>-</td><td>-</td><td>0.06021</td><td>0.06238</td></tr><tr><td>Calcium signaling pathway</td><td>0.05076</td><td>0.09113</td><td>-</td><td>-</td></tr><tr><td>Inositol phosphate metabolism</td><td>0.00107</td><td>-</td><td>-</td><td>-</td></tr><tr><td>FoxO signaling pathway</td><td>0.00333</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Phosphatidylinositol signaling system</td><td>0.02371</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Pentose and glucuronate interconversions</td><td>-</td><td>-</td><td>-</td><td>0.02564</td></tr><tr><td>Dorso-ventral axis formation</td><td>-</td><td>0.02705</td><td>-</td><td>-</td></tr><tr><td>Cysteine and methionine metabolism</td><td>-</td><td>-</td><td>-</td><td>0.02953</td></tr><tr><td>Inflammatory mediator regulation of TRP channels</td><td>-</td><td>-</td><td>-</td><td>0.03010</td></tr><tr><td>Pyruvate metabolism</td><td>-</td><td>-</td><td>-</td><td>0.03372</td></tr><tr><td>Adherens junction</td><td>-</td><td>-</td><td>0.03801</td><td>-</td></tr><tr><td>Fc gamma R-mediated phagocytosis</td><td>-</td><td>-</td><td>0.07454</td><td>-</td></tr><tr><td>Regulation of actin cytoskeleton</td><td>-</td><td>-</td><td>-</td><td>0.08180</td></tr><tr><td>Metabolic pathways</td><td>0.08439</td><td>-</td><td>-</td><td>-</td></tr><tr><td>PPAR signaling pathway</td><td>-</td><td>-</td><td>0.08586</td><td>-</td></tr></tbody></table>
among others CXCL12, CLDN4, Platelet F11 receptor (F11R), or MMP2. Components of the TNF signaling pathway were mainly downregulated by statins in MCs in our experiment. This downregulation could be observed on the level of intracellular signaling molecules such as MAP2K3 but also more importantly on the effector level of this signaling pathway. The inflammatory molecules IL6, CCL20, CSF2, and ICAM1 for example were significantly downregulated after statin treatment. Also endothelin 1 (EDN1) was markedly downregulated by both statins in MCs.

Another sign of evidence for the impact on cell proliferation and cell differentiation of statins is the fact that a number of components of cell-matrix adhesion proteins is dysregulated as compared to unstimulated cells. Various actinin isoforms (ACTN1, ACTN3, ACTN4) were found to be downregulated after statin treatment next to downregulation of caveolin 1 (CAV1) and 3 (CAV3), linking integrins to intracellular signaling pathways thus promoting cell cycle progression. In addition, filamin (FLNA), vinculin (VCL) and zyxin (ZYX) were found to be downregulated after statin treatment. Extracellular molecules like integrin beta 4 (ITGB4), integrin alpha 7 (ITGA7) or the hepatocyte growth factor (HGF) on the other hand were significantly upregulated by statins.

The expression of another major renoprotective factor, namely heme oxygenase 1 (HMOX1), was significantly induced by atorvastatin but also upregulated by rosuvastatin in both cell types under study.

The following eight biomarkers that are linked to the affected molecular processes after statin treatment in renal cells were available in patients at baseline in the AURORA study (HGF, IL6, MMP2, FAS, NTF3, HMOX1, FST, and F11R) which were subsequently evaluated regarding the effect of rosuvastatin in hemodialysis patients (see Fig. 2).

## 3.1. Biomarker-based patient phenotypes are associated with outcome in AURORA

We identified four patient phenotype clusters in the AURORA study based on the residual values of the eight measured biomarkers. The phenotype clusters ranged from a minimum of 35 to a maximum 185 patients (Table 2).

The largest cluster (phenotype 4, n = 185), labelled reference phenotype, had the lowest concentration of hs-CRP (p < 0.01) and showed the highest levels of both aspartate aminotransferase (p < 0.01) as well as creatine kinase (p < 0.01). Phenotype 1 exhibited

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Heatmap of pathway-related genes and biomarkers</figcaption></figure>

Log fold-change values of DEGs assigned to enriched molecular pathways are depicted as heatmaps with purple indicating significant upregulation and yellow indicating significant downregulation after statin treatment. Selected molecular biomarkers that have been evaluated in the AURORA study are depicted in bold font with a border. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)
among others CXCL12, CLDN4, Platelet F11 receptor (F11R), or MMP2. Components of the TNF signaling pathway were mainly downregulated by statins in MCs in our experiment. This downregulation could be observed on the level of intracellular signaling molecules such as MAP2K3 but also more importantly on the effector level of this signaling pathway. The inflammatory molecules IL6, CCL20, CSF2, and ICAM1 for example were significantly downregulated after statin treatment. Also endothelin 1 (EDN1) was markedly downregulated by both statins in MCs.

Another sign of evidence for the impact on cell proliferation and cell differentiation of statins is the fact that a number of components of cell-matrix adhesion proteins is dysregulated as compared to unstimulated cells. Various actinin isoforms (ACTN1, ACTN3, ACTN4) were found to be downregulated after statin treatment next to downregulation of caveolin 1 (CAV1) and 3 (CAV3), linking integrins to intracellular signaling pathways thus promoting cell cycle progression. In addition, filamin (FLNA), vinculin (VCL) and zyxin (ZYX) were found to be downregulated after statin treatment. Extracellular molecules like integrin beta 4 (ITGB4), integrin alpha 7 (ITGA7) or the hepatocyte growth factor (HGF) on the other hand were significantly upregulated by statins.

The expression of another major renoprotective factor, namely heme oxygenase 1 (HMOX1), was significantly induced by atorvastatin but also upregulated by rosuvastatin in both cell types under study.

The following eight biomarkers that are linked to the affected molecular processes after statin treatment in renal cells were available in patients at baseline in the AURORA study (HGF, IL6, MMP2, FAS, NTF3, HMOX1, FST, and F11R) which were subsequently evaluated regarding the effect of rosuvastatin in hemodialysis patients (see Fig. 2).

## 3.1. Biomarker-based patient phenotypes are associated with outcome in AURORA

We identified four patient phenotype clusters in the AURORA study based on the residual values of the eight measured biomarkers. The phenotype clusters ranged from a minimum of 35 to a maximum 185 patients (Table 2).

The largest cluster (phenotype 4, n = 185), labelled reference phenotype, had the lowest concentration of hs-CRP (p < 0.01) and showed the highest levels of both aspartate aminotransferase (p < 0.01) as well as creatine kinase (p < 0.01). Phenotype 1 exhibited

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Heatmap of pathway-related genes and biomarkers</figcaption></figure>

Log fold-change values of DEGs assigned to enriched molecular pathways are depicted as heatmaps with purple indicating significant upregulation and yellow indicating significant downregulation after statin treatment. Selected molecular biomarkers that have been evaluated in the AURORA study are depicted in bold font with a border. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)
<table><caption>**Table 2**
Baseline characteristics of the analyzed AURORA samples and the four phenotypes.</caption><thead><tr><th></th><th>Overall population (n = 398)</th><th>Phenotype 1 (n = 125)</th><th>Phenotype 2 (n = 35)</th><th>Phenotype 3 (n = 53)</th><th>Phenotype 4 (n = 185)</th><th>p-value</th></tr></thead><tbody><tr><td>Baseline characteristics</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Female gender (%)</td><td>83 (20.9)</td><td>28 (22.4)</td><td>10 (28.6)</td><td>14 (26.4)</td><td>31 (16.8)</td><td>0.23</td></tr><tr><td>Age (years)</td><td>67 [60, 74]</td><td>68 [62, 74]</td><td>69 [62, 74]</td><td>68 [59, 74]</td><td>66 [59, 73]</td><td>0.29</td></tr><tr><td>Years on RRT</td><td>2.7 [1.1, 6.0]</td><td>3.7 [1.7, 6.3]</td><td>1.9 [1.1, 3.8]</td><td>7.0 [1.9, 13.6]</td><td>2.1 [0.7, 4.3]</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Calculated K<sub>f</sub>/V</td><td>1.1 [1.0, 1.3]</td><td>1.1 [1.0, 1.3]</td><td>1.1 [0.9, 1.3]</td><td>1.1 [1.0, 1.2]</td><td>1.1 [1.0, 1.2]</td><td>0.75</td></tr><tr><td>Albumin (g/L)</td><td>39.3 (3.1)</td><td>39.6 (2.9)</td><td>37.8 (2.8)</td><td>39.0 (3.3)</td><td>39.5 (3.1)</td><td><strong>0.01</strong></td></tr><tr><td>Hemoglobin (g/dL)</td><td>11.8 (1.5)</td><td>11.7 (1.6)</td><td>12.2 (1.4)</td><td>11.9 (1.5)</td><td>11.7 (1.5)</td><td>0.33</td></tr><tr><td>hs-CRP (mg/L)</td><td>1.2 (1.2)</td><td>1.6 (1.3)</td><td>1.8 (1.4)</td><td>1.2 (1.2)</td><td>0.8 (0.9)</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Fasting Glucose level (mmol/L)</td><td>4.9 [4.6, 5.6]</td><td>5.0 [4.6, 5.7]</td><td>4.9 [4.3, 6.0]</td><td>4.9 [4.7, 5.4]</td><td>4.9 [4.6, 5.6]</td><td>0.83</td></tr><tr><td>Creatine Kinase (IU/L)</td><td>56 [36, 84]</td><td>50 [33, 82]</td><td>47 [32, 62]</td><td>49 [36, 75]</td><td>65 [41, 91]</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Aspartate Aminotransferase (IU/L)</td><td>16.0 (8.7)</td><td>15.6 (8.3)</td><td>15.9 (5.4)</td><td>20.2 (12.0)</td><td>15.2 (7.9)</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Calculated LDL-Cholesterol (mmol/L)</td><td>2.6 (0.9)</td><td>2.5 (0.8)</td><td>2.8 (1.1)</td><td>2.7 (0.9)</td><td>2.6 (0.8)</td><td>0.12</td></tr><tr><td>BMI (kg/m2)</td><td>25.0 (3.4)</td><td>24.8 (3.4)</td><td>24.4 (3.3)</td><td>24.5 (3.3)</td><td>25.4 (3.5)</td><td>0.2</td></tr><tr><td>Systolic Blood Pressure (mmHg)</td><td>138.7 (24.4)</td><td>142.8 (24.2)</td><td>129.5 (17.7)</td><td>136.2 (25.8)</td><td>138.5 (24.7)</td><td><strong>0.03</strong></td></tr><tr><td>Diastolic Blood Pressure (mmHg)</td><td>76.4 (12.7)</td><td>77.3 (12.9)</td><td>73.8 (12.4)</td><td>75.7 (13.3)</td><td>76.6 (12.6)</td><td>0.52</td></tr><tr><td>Pulse pressure (mmHg)</td><td>60 [50, 74]</td><td>69 [55, 77]</td><td>50 [43, 67]</td><td>57 [49, 75]</td><td>60 [50, 73]</td><td><strong>0.01</strong></td></tr><tr><td>Medical history, n (%)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Current smoker</td><td>66 (16.6)</td><td>25 (20.0)</td><td>9 (25.7)</td><td>3 (5.7)</td><td>29 (15.7)</td><td><strong>0.05</strong></td></tr><tr><td>Diabetes</td><td>84 (21.1)</td><td>27 (21.6)</td><td>9 (25.7)</td><td>9 (17.0)</td><td>39 (21.1)</td><td>0.8</td></tr><tr><td>History of CVD</td><td>172 (43.2)</td><td>58 (46.4)</td><td>15 (42.9)</td><td>21 (39.6)</td><td>78 (42.2)</td><td>0.83</td></tr><tr><td>Peripheral artery disease</td><td>69 (17.3)</td><td>27 (21.6)</td><td>6 (17.1)</td><td>7 (13.2)</td><td>29 (15.7)</td><td>0.46</td></tr><tr><td>History of coronary heart disease</td><td>74 (18.6)</td><td>27 (21.6)</td><td>9 (25.7)</td><td>8 (15.1)</td><td>30 (16.2)</td><td>0.39</td></tr><tr><td>History of vascular disease</td><td>98 (24.6)</td><td>31 (24.8)</td><td>7 (20.0)</td><td>14 (26.4)</td><td>46 (24.9)</td><td>0.92</td></tr><tr><td>History of neurovascular disease</td><td>68 (17.1)</td><td>23 (18.4)</td><td>6 (17.1)</td><td>4 (7.5)</td><td>35 (18.9)</td><td>0.26</td></tr><tr><td>Cause of ESRD, n (%)</td><td></td><td></td><td></td><td></td><td></td><td>0.21</td></tr><tr><td>Diabetes</td><td>54 (13.6)</td><td>18 (14.4)</td><td>6 (17.1)</td><td>5 (9.4)</td><td>25 (13.5)</td><td></td></tr><tr><td>Genetic conditions</td><td>51 (12.8)</td><td>16 (12.8)</td><td>5 (14.3)</td><td>2 (3.8)</td><td>28 (15.1)</td><td></td></tr><tr><td>Glomerulonephritis or Vasculitis</td><td>69 (17.3)</td><td>24 (19.2)</td><td>3 (8.6)</td><td>12 (22.6)</td><td>30 (16.2)</td><td></td></tr><tr><td>Nephropathy or Nephrosclerosis</td><td>94 (23.6)</td><td>32 (25.6)</td><td>5 (14.3)</td><td>9 (17.0)</td><td>48 (25.9)</td><td></td></tr><tr><td>Pyelonephritis or Interstitial</td><td>63 (15.8)</td><td>19 (15.2)</td><td>9 (25.7)</td><td>13 (24.5)</td><td>22 (11.9)</td><td></td></tr><tr><td>Unknown/unspecified</td><td>44 (11.1)</td><td>10 (8.0)</td><td>3 (8.6)</td><td>8 (15.1)</td><td>23 (12.4)</td><td></td></tr><tr><td>Other</td><td>23 (5.8)</td><td>6 (4.8)</td><td>4 (11.4)</td><td>4 (7.5)</td><td>9 (4.9)</td><td></td></tr></tbody></table>

Key clinical parameters of the investigated sample set of the AURORA trial are given for the whole cohort as well as for the four identified phenotypes. Mean values plus standard deviation are provided for normally distributed continuous variables with median values plus interquartile ranges for non-normally distributed continuous data. Categorical variables are expressed as frequencies (percentages). Years on RRT = years on renal replacement therapy; hs-CRP = high sensitivity C-reactive protein; BMI = body mass index; ESRD = end-stage renal disease.

higher pulse pressure (p = 0.01), whereas phenotype 2 showed higher proportion of smokers (p = 0.05). Patients in phenotype 3 had significantly longer dialysis vintage (7.0 years as compared to 2.7 years of the overall population) (p < 0.01).

The decision tree that was built with the eight biomarkers identified distinctive biological profiles (Fig. 3). The global accuracy of the decision tree was 81.9 %. Phenotype 1 had low HGF (a protein upregulated by statins) and high IL6/NTF3 (proteins down-regulated by statins) levels or moderately low HGF and high MMP2 levels respectively. In contrast, phenotype 2 had high HGF (a protein upregulated by statins) and low FAS (a protein down-regulated by statins) levels.

Phenotype 1 was significantly associated with increased rates of CV death, all-cause mortality, and MACE in crude models [HR = 1.78; 95 % CI (1.26–2.52), HR = 1.72; 95 % CI (1.28–2.30), and HR = 1.62; 95 % CI (1.19–2.22), respectively] as compared with the reference phenotype (Table 3). After adjustment for clinical and biological factors, phenotype 1 remained significantly associated with CV death and MACE [adjusted HR = 1.55; 95 % CI (1.03–2.34) and HR = 1.48; 95 % CI (1.03–2.14), respectively] and showed borderline association with all-cause mortality [adjusted HR = 1.35; 95 % CI (0.95–1.91) (Table 3).

## 3.2. Biomarker-based patient phenotypes are associated with response to statin therapy

Rosuvastatin was associated with significantly lower risk of all-cause mortality at 3 years [risk difference = −21.0; 95 % CI (−38.3 to −3.8); p = 0.017] in phenotype 1. In contrast, rosuvastatin was significantly associated with increased risk of all-cause mortality at 1 year and tended to be associated with high all-cause mortality at 2 and 3 years in phenotype 2 (Fig. 4). No significant differences regarding treatment response as compared with placebo were identified for phenotypes 3 and 4 [all p > 0.5 at risk difference 1, 2, and 3 years] (Fig. 4). A similar pattern was observed when considering MACE and CV death.
<table><caption>**Table 2**
Baseline characteristics of the analyzed AURORA samples and the four phenotypes.</caption><thead><tr><th></th><th>Overall population (n = 398)</th><th>Phenotype 1 (n = 125)</th><th>Phenotype 2 (n = 35)</th><th>Phenotype 3 (n = 53)</th><th>Phenotype 4 (n = 185)</th><th>p-value</th></tr></thead><tbody><tr><td>Baseline characteristics</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Female gender (%)</td><td>83 (20.9)</td><td>28 (22.4)</td><td>10 (28.6)</td><td>14 (26.4)</td><td>31 (16.8)</td><td>0.23</td></tr><tr><td>Age (years)</td><td>67 [60, 74]</td><td>68 [62, 74]</td><td>69 [62, 74]</td><td>68 [59, 74]</td><td>66 [59, 73]</td><td>0.29</td></tr><tr><td>Years on RRT</td><td>2.7 [1.1, 6.0]</td><td>3.7 [1.7, 6.3]</td><td>1.9 [1.1, 3.8]</td><td>7.0 [1.9, 13.6]</td><td>2.1 [0.7, 4.3]</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Calculated K<sub>f</sub>/V</td><td>1.1 [1.0, 1.3]</td><td>1.1 [1.0, 1.3]</td><td>1.1 [0.9, 1.3]</td><td>1.1 [1.0, 1.2]</td><td>1.1 [1.0, 1.2]</td><td>0.75</td></tr><tr><td>Albumin (g/L)</td><td>39.3 (3.1)</td><td>39.6 (2.9)</td><td>37.8 (2.8)</td><td>39.0 (3.3)</td><td>39.5 (3.1)</td><td><strong>0.01</strong></td></tr><tr><td>Hemoglobin (g/dL)</td><td>11.8 (1.5)</td><td>11.7 (1.6)</td><td>12.2 (1.4)</td><td>11.9 (1.5)</td><td>11.7 (1.5)</td><td>0.33</td></tr><tr><td>hs-CRP (mg/L)</td><td>1.2 (1.2)</td><td>1.6 (1.3)</td><td>1.8 (1.4)</td><td>1.2 (1.2)</td><td>0.8 (0.9)</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Fasting Glucose level (mmol/L)</td><td>4.9 [4.6, 5.6]</td><td>5.0 [4.6, 5.7]</td><td>4.9 [4.3, 6.0]</td><td>4.9 [4.7, 5.4]</td><td>4.9 [4.6, 5.6]</td><td>0.83</td></tr><tr><td>Creatine Kinase (IU/L)</td><td>56 [36, 84]</td><td>50 [33, 82]</td><td>47 [32, 62]</td><td>49 [36, 75]</td><td>65 [41, 91]</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Aspartate Aminotransferase (IU/L)</td><td>16.0 (8.7)</td><td>15.6 (8.3)</td><td>15.9 (5.4)</td><td>20.2 (12.0)</td><td>15.2 (7.9)</td><td><strong>&lt;0.01</strong></td></tr><tr><td>Calculated LDL-Cholesterol (mmol/L)</td><td>2.6 (0.9)</td><td>2.5 (0.8)</td><td>2.8 (1.1)</td><td>2.7 (0.9)</td><td>2.6 (0.8)</td><td>0.12</td></tr><tr><td>BMI (kg/m2)</td><td>25.0 (3.4)</td><td>24.8 (3.4)</td><td>24.4 (3.3)</td><td>24.5 (3.3)</td><td>25.4 (3.5)</td><td>0.2</td></tr><tr><td>Systolic Blood Pressure (mmHg)</td><td>138.7 (24.4)</td><td>142.8 (24.2)</td><td>129.5 (17.7)</td><td>136.2 (25.8)</td><td>138.5 (24.7)</td><td><strong>0.03</strong></td></tr><tr><td>Diastolic Blood Pressure (mmHg)</td><td>76.4 (12.7)</td><td>77.3 (12.9)</td><td>73.8 (12.4)</td><td>75.7 (13.3)</td><td>76.6 (12.6)</td><td>0.52</td></tr><tr><td>Pulse pressure (mmHg)</td><td>60 [50, 74]</td><td>69 [55, 77]</td><td>50 [43, 67]</td><td>57 [49, 75]</td><td>60 [50, 73]</td><td><strong>0.01</strong></td></tr><tr><td>Medical history, n (%)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Current smoker</td><td>66 (16.6)</td><td>25 (20.0)</td><td>9 (25.7)</td><td>3 (5.7)</td><td>29 (15.7)</td><td><strong>0.05</strong></td></tr><tr><td>Diabetes</td><td>84 (21.1)</td><td>27 (21.6)</td><td>9 (25.7)</td><td>9 (17.0)</td><td>39 (21.1)</td><td>0.8</td></tr><tr><td>History of CVD</td><td>172 (43.2)</td><td>58 (46.4)</td><td>15 (42.9)</td><td>21 (39.6)</td><td>78 (42.2)</td><td>0.83</td></tr><tr><td>Peripheral artery disease</td><td>69 (17.3)</td><td>27 (21.6)</td><td>6 (17.1)</td><td>7 (13.2)</td><td>29 (15.7)</td><td>0.46</td></tr><tr><td>History of coronary heart disease</td><td>74 (18.6)</td><td>27 (21.6)</td><td>9 (25.7)</td><td>8 (15.1)</td><td>30 (16.2)</td><td>0.39</td></tr><tr><td>History of vascular disease</td><td>98 (24.6)</td><td>31 (24.8)</td><td>7 (20.0)</td><td>14 (26.4)</td><td>46 (24.9)</td><td>0.92</td></tr><tr><td>History of neurovascular disease</td><td>68 (17.1)</td><td>23 (18.4)</td><td>6 (17.1)</td><td>4 (7.5)</td><td>35 (18.9)</td><td>0.26</td></tr><tr><td>Cause of ESRD, n (%)</td><td></td><td></td><td></td><td></td><td></td><td>0.21</td></tr><tr><td>Diabetes</td><td>54 (13.6)</td><td>18 (14.4)</td><td>6 (17.1)</td><td>5 (9.4)</td><td>25 (13.5)</td><td></td></tr><tr><td>Genetic conditions</td><td>51 (12.8)</td><td>16 (12.8)</td><td>5 (14.3)</td><td>2 (3.8)</td><td>28 (15.1)</td><td></td></tr><tr><td>Glomerulonephritis or Vasculitis</td><td>69 (17.3)</td><td>24 (19.2)</td><td>3 (8.6)</td><td>12 (22.6)</td><td>30 (16.2)</td><td></td></tr><tr><td>Nephropathy or Nephrosclerosis</td><td>94 (23.6)</td><td>32 (25.6)</td><td>5 (14.3)</td><td>9 (17.0)</td><td>48 (25.9)</td><td></td></tr><tr><td>Pyelonephritis or Interstitial</td><td>63 (15.8)</td><td>19 (15.2)</td><td>9 (25.7)</td><td>13 (24.5)</td><td>22 (11.9)</td><td></td></tr><tr><td>Unknown/unspecified</td><td>44 (11.1)</td><td>10 (8.0)</td><td>3 (8.6)</td><td>8 (15.1)</td><td>23 (12.4)</td><td></td></tr><tr><td>Other</td><td>23 (5.8)</td><td>6 (4.8)</td><td>4 (11.4)</td><td>4 (7.5)</td><td>9 (4.9)</td><td></td></tr></tbody></table>

Key clinical parameters of the investigated sample set of the AURORA trial are given for the whole cohort as well as for the four identified phenotypes. Mean values plus standard deviation are provided for normally distributed continuous variables with median values plus interquartile ranges for non-normally distributed continuous data. Categorical variables are expressed as frequencies (percentages). Years on RRT = years on renal replacement therapy; hs-CRP = high sensitivity C-reactive protein; BMI = body mass index; ESRD = end-stage renal disease.

higher pulse pressure (p = 0.01), whereas phenotype 2 showed higher proportion of smokers (p = 0.05). Patients in phenotype 3 had significantly longer dialysis vintage (7.0 years as compared to 2.7 years of the overall population) (p < 0.01).

The decision tree that was built with the eight biomarkers identified distinctive biological profiles (Fig. 3). The global accuracy of the decision tree was 81.9 %. Phenotype 1 had low HGF (a protein upregulated by statins) and high IL6/NTF3 (proteins down-regulated by statins) levels or moderately low HGF and high MMP2 levels respectively. In contrast, phenotype 2 had high HGF (a protein upregulated by statins) and low FAS (a protein down-regulated by statins) levels.

Phenotype 1 was significantly associated with increased rates of CV death, all-cause mortality, and MACE in crude models [HR = 1.78; 95 % CI (1.26–2.52), HR = 1.72; 95 % CI (1.28–2.30), and HR = 1.62; 95 % CI (1.19–2.22), respectively] as compared with the reference phenotype (Table 3). After adjustment for clinical and biological factors, phenotype 1 remained significantly associated with CV death and MACE [adjusted HR = 1.55; 95 % CI (1.03–2.34) and HR = 1.48; 95 % CI (1.03–2.14), respectively] and showed borderline association with all-cause mortality [adjusted HR = 1.35; 95 % CI (0.95–1.91) (Table 3).

## 3.2. Biomarker-based patient phenotypes are associated with response to statin therapy

Rosuvastatin was associated with significantly lower risk of all-cause mortality at 3 years [risk difference = −21.0; 95 % CI (−38.3 to −3.8); p = 0.017] in phenotype 1. In contrast, rosuvastatin was significantly associated with increased risk of all-cause mortality at 1 year and tended to be associated with high all-cause mortality at 2 and 3 years in phenotype 2 (Fig. 4). No significant differences regarding treatment response as compared with placebo were identified for phenotypes 3 and 4 [all p > 0.5 at risk difference 1, 2, and 3 years] (Fig. 4). A similar pattern was observed when considering MACE and CV death.
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Biomarker relevance per patient phenotype</figcaption></figure>

A decision tree was constructed to categorize phenotypes based on biomarker levels. Color-coding of nodes in the decision tree reflect the dominant phenotype (red = phenotype 1; green = phenotype 2; yellow = phenotype 3; blue = phenotype 4). Simplified rulesets for all terminal nodes holding more than 10 patients of the dominant phenotype are given on the right. Phenotype 3 for example is characterized by high levels of HGF as well as FAS. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

<table><caption>**Table 3**

Association of patient phenotypes with outcome.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Univariable</th><th colspan="2">Multivariable</th></tr><tr><th>HR (95 % CI)</th><th>p-value</th><th>HR (95 % CI)</th><th>p-value</th></tr></thead><tbody><tr><td><strong>CV death (n = 155)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.78 (1.26-2.52)</td><td><strong>0.001</strong></td><td>1.55 (1.03-2.34)</td><td><strong>0.036</strong></td></tr><tr><td>Phenotype 2</td><td>0.71 (0.36-1.43)</td><td>0.34</td><td>0.74 (0.36-1.52)</td><td>0.42</td></tr><tr><td>Phenotype 3</td><td>0.94 (0.55-1.61)</td><td>0.82</td><td>0.83 (0.46-1.50)</td><td>0.54</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>All-cause mortality (n = 230)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.72 (1.28-2.30)</td><td><strong>0.0003</strong></td><td>1.35 (0.95-1.91)</td><td>0.09</td></tr><tr><td>Phenotype 2</td><td>0.86 (0.51-1.47)</td><td>0.59</td><td>0.80 (0.46-1.39)</td><td>0.44</td></tr><tr><td>Phenotype 3</td><td>1.18 (0.78-1.77)</td><td>0.43</td><td>1.04 (0.67-1.62)</td><td>0.86</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>MACE (n = 199)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.62 (1.19-2.22)</td><td><strong>0.002</strong></td><td>1.48 (1.03-2.14)</td><td><strong>0.035</strong></td></tr><tr><td>Phenotype 2</td><td>0.73 (0.40-1.33)</td><td>0.30</td><td>0.78 (0.42-1.46)</td><td>0.43</td></tr><tr><td>Phenotype 3</td><td>1.18 (0.77-1.83)</td><td>0.45</td><td>1.12 (0.69-1.79)</td><td>0.65</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Hazard ratios of univariate and multivariate Cox proportional hazards regression models are provided for the three phenotypes with respect to the reference phenotype 4. Parameters included in the multivariate analysis for adjustment included age, albumin, log hs-CRP, history of CV disease, diabetes mellitus, dialysis vintage and log BNP.

## 4. Discussion

In this translational study we investigated the impact of statin treatment on renal cells in-vitro on the transcriptional levels and evaluated a set of biomarkers being linked to statin-affected molecular pathways in plasma samples of patients from the AURORA trial regarding their association with disease outcome and response to statin treatment. We show that statins had a larger impact on gene expression changes in MCs as compared to HK2 cells. Atorvastatin had a larger effect on expression changes than rosuvastatin with a number of transcripts, however, being affected in a consistent way by both statins being an indication for a strong drug class effect. Key affected molecular pathways after statin treatment in renal cells included the fibrosis-associated TGFB signaling pathway, the inflammatory processes of TNF signaling and leukocyte trans-endothelial migration, focal adhesion, mineral absorption as well as MAPK signaling. These findings are in-line with recently reported effects of statin treatment in the context of contrast-induced nephropathy [20]. When building phenotypes based on the eight selected plasma protein biomarkers in the AURORA trial, phenotype 1 was associated with worse outcome overall, but appeared to have improved outcome when treated with statins. Interestingly, these patients had low HGF levels at baseline, a marker that is up-regulated by statins, and high IL6 and NTF3 levels, two proteins that are
<figure><img src="image_6.png" /><figcaption>**Fig. 3.** Biomarker relevance per patient phenotype</figcaption></figure>

A decision tree was constructed to categorize phenotypes based on biomarker levels. Color-coding of nodes in the decision tree reflect the dominant phenotype (red = phenotype 1; green = phenotype 2; yellow = phenotype 3; blue = phenotype 4). Simplified rulesets for all terminal nodes holding more than 10 patients of the dominant phenotype are given on the right. Phenotype 3 for example is characterized by high levels of HGF as well as FAS. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

<table><caption>**Table 3**

Association of patient phenotypes with outcome.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Univariable</th><th colspan="2">Multivariable</th></tr><tr><th>HR (95 % CI)</th><th>p-value</th><th>HR (95 % CI)</th><th>p-value</th></tr></thead><tbody><tr><td><strong>CV death (n = 155)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.78 (1.26-2.52)</td><td><strong>0.001</strong></td><td>1.55 (1.03-2.34)</td><td><strong>0.036</strong></td></tr><tr><td>Phenotype 2</td><td>0.71 (0.36-1.43)</td><td>0.34</td><td>0.74 (0.36-1.52)</td><td>0.42</td></tr><tr><td>Phenotype 3</td><td>0.94 (0.55-1.61)</td><td>0.82</td><td>0.83 (0.46-1.50)</td><td>0.54</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>All-cause mortality (n = 230)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.72 (1.28-2.30)</td><td><strong>0.0003</strong></td><td>1.35 (0.95-1.91)</td><td>0.09</td></tr><tr><td>Phenotype 2</td><td>0.86 (0.51-1.47)</td><td>0.59</td><td>0.80 (0.46-1.39)</td><td>0.44</td></tr><tr><td>Phenotype 3</td><td>1.18 (0.78-1.77)</td><td>0.43</td><td>1.04 (0.67-1.62)</td><td>0.86</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>MACE (n = 199)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Per phenotype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenotype 1</td><td>1.62 (1.19-2.22)</td><td><strong>0.002</strong></td><td>1.48 (1.03-2.14)</td><td><strong>0.035</strong></td></tr><tr><td>Phenotype 2</td><td>0.73 (0.40-1.33)</td><td>0.30</td><td>0.78 (0.42-1.46)</td><td>0.43</td></tr><tr><td>Phenotype 3</td><td>1.18 (0.77-1.83)</td><td>0.45</td><td>1.12 (0.69-1.79)</td><td>0.65</td></tr><tr><td>Phenotype 4</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Hazard ratios of univariate and multivariate Cox proportional hazards regression models are provided for the three phenotypes with respect to the reference phenotype 4. Parameters included in the multivariate analysis for adjustment included age, albumin, log hs-CRP, history of CV disease, diabetes mellitus, dialysis vintage and log BNP.

## 4. Discussion

In this translational study we investigated the impact of statin treatment on renal cells in-vitro on the transcriptional levels and evaluated a set of biomarkers being linked to statin-affected molecular pathways in plasma samples of patients from the AURORA trial regarding their association with disease outcome and response to statin treatment. We show that statins had a larger impact on gene expression changes in MCs as compared to HK2 cells. Atorvastatin had a larger effect on expression changes than rosuvastatin with a number of transcripts, however, being affected in a consistent way by both statins being an indication for a strong drug class effect. Key affected molecular pathways after statin treatment in renal cells included the fibrosis-associated TGFB signaling pathway, the inflammatory processes of TNF signaling and leukocyte trans-endothelial migration, focal adhesion, mineral absorption as well as MAPK signaling. These findings are in-line with recently reported effects of statin treatment in the context of contrast-induced nephropathy [20]. When building phenotypes based on the eight selected plasma protein biomarkers in the AURORA trial, phenotype 1 was associated with worse outcome overall, but appeared to have improved outcome when treated with statins. Interestingly, these patients had low HGF levels at baseline, a marker that is up-regulated by statins, and high IL6 and NTF3 levels, two proteins that are
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Association of treatment response with outcome
Kaplan-Meier plots for the four phenotypes showing significant associations with the outcome all-cause mortality when comparing rosuvastatin treated vs untreated patients per phenotype.

\* p-value <0.1 \*\* p-value <0.05.</figcaption></figure>

down-regulated by statins. In contrast, phenotype 2 was characterized by high HGF and low FAS levels respectively. These patients experienced worse outcome when treated with statins as compared to untreated patients with this molecular phenotype.

Traversing finding from the in-vitro setting to the human situation is in most cases challenging and requires quality control at various steps. We first investigated the changes in gene expression of the direct statin drug target and associated proteins. Based on data from the Human Protein Atlas, there is low protein expression of the direct drug target of statins, HMGCR, in glomeruli and medium protein expression in tubular cells [21]. HMGCR was significantly upregulated in MCs in our analysis on the transcriptional level after statin treatment. This upregulation of HMGCR is most likely a cellular counter mechanism due to reduced intracellular cholesterol concentrations being in-line with reports from other groups [22,23]. In addition to HMGCR, also the receptor for low density lipoprotein (LDLR) was reported to be upregulated after statin treatment, a phenomenon that we also detected in MCs but not in HK2 cells [23]. Lower intracellular cholesterol levels not only increase expression of LDLR but also expression of the proprotein convertase subtilisin/kexin type 9 (PCSK9), another molecule with induced gene expression in MCs in our study. Squalene epoxidase (SQLE), one of the rate-limiting enzymes in the sterol biosynthesis pathway was also significantly upregulated in MCs after statin treatment in our experiments [24]. These findings solidified the conduct of the performed in-vitro gene expression experiments. Further evidence regarding the quality of the in-vitro data was found when inspecting the set of genes that were significantly regulated
<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Association of treatment response with outcome
Kaplan-Meier plots for the four phenotypes showing significant associations with the outcome all-cause mortality when comparing rosuvastatin treated vs untreated patients per phenotype.

\* p-value <0.1 \*\* p-value <0.05.</figcaption></figure>

down-regulated by statins. In contrast, phenotype 2 was characterized by high HGF and low FAS levels respectively. These patients experienced worse outcome when treated with statins as compared to untreated patients with this molecular phenotype.

Traversing finding from the in-vitro setting to the human situation is in most cases challenging and requires quality control at various steps. We first investigated the changes in gene expression of the direct statin drug target and associated proteins. Based on data from the Human Protein Atlas, there is low protein expression of the direct drug target of statins, HMGCR, in glomeruli and medium protein expression in tubular cells [21]. HMGCR was significantly upregulated in MCs in our analysis on the transcriptional level after statin treatment. This upregulation of HMGCR is most likely a cellular counter mechanism due to reduced intracellular cholesterol concentrations being in-line with reports from other groups [22,23]. In addition to HMGCR, also the receptor for low density lipoprotein (LDLR) was reported to be upregulated after statin treatment, a phenomenon that we also detected in MCs but not in HK2 cells [23]. Lower intracellular cholesterol levels not only increase expression of LDLR but also expression of the proprotein convertase subtilisin/kexin type 9 (PCSK9), another molecule with induced gene expression in MCs in our study. Squalene epoxidase (SQLE), one of the rate-limiting enzymes in the sterol biosynthesis pathway was also significantly upregulated in MCs after statin treatment in our experiments [24]. These findings solidified the conduct of the performed in-vitro gene expression experiments. Further evidence regarding the quality of the in-vitro data was found when inspecting the set of genes that were significantly regulated
by both statins in both cell lines. Among the significantly downregulated molecules was the connective tissue growth factor (CTGF), a pro-fibrotic stimulus, which has been shown to be affected by statin treatment, both in renal as well as endothelial cells [25–27]. The Kruppel like factor 2 (KLF2) was among the significantly upregulated molecules by both drugs in both cell lines. KLF2 was reported to be downregulated after high glucose treatment in human endothelial cells and reduced KLF2 expression was associated with aggravated endothelial injury in DN [28,29].

The affected genes in the cell culture setup that were responsible for patient stratification into the four endophenotypes in the AURORA were HGF, IL6, MMP2, FAS, and NTF3. HGF is a cytokine with a number of regulatory functions in the context of immunomodulation, angiogenesis, morphogenesis, as well as tissue repair [30]. HGF was shown to be elevated in hemodialysis patients and a link to cardiovascular risk was postulated although it is still under debate whether HGF increases CVD risk or whether it is actually involved in counterbalancing endothelial damage [30–32]. Burgazli and colleagues reported on the role of statins to modulate HGF-induced proliferation of human umbilical vein endothelial cells [33]. Our data rather indicate that HGF has a protective role with phenotype clusters showing higher values of HGF being linked to improved outcome whereas phenotype cluster 1 with the worst prognosis showed the lowest HGF values. This cluster was also characterized by high MMP2 levels, a molecule involved in fibrotic processes. HGF was previously reported to induce MMP2 levels in human glomeruli next to reducing levels of TGFB1 and type IV collagens thus exerting its anti-fibrotic effects [34]. MMP2 levels are elevated in ESRD with hemodialysis reducing levels of MMP2 [35]. Serum levels of MMP2 are associated with mortality in hemodialysis patients [36]. An inhibitory impact of statin treatment on MMP2 levels was also reported by de Araujo Junior and colleagues as well as by Janardhanan et al. respectively [37,38]. Statins were also shown to reduce levels of pro-inflammatory IL6 in human vascular smooth muscle cells [39] and monocytes from hypercholesterolemic patients [40]. In addition statins reduced IL6 levels in patients on hemodialysis or the pre-dialysis state [41–43]. High IL6 levels in our analysis where characteristic for phenotype 1 holding those patients with the worst prognosis but on the other hand showing the most benefit of statin treatment.

The analysis performed in the AURORA trial strengthens the validity of our results. First, phenotype 1 (n = 125) had the most opposing biological profile (low HGF, high IL6, and NTF3) to the expected biological modifications triggered by statins. Patients of this phenotype had the worst prognosis with respect to overall survival, CV death and also risk of MACE. In addition, this profile appeared to have some benefit from statin treatment when considering all-cause mortality. This suggests that patients in whom the biomarkers are most opposed at baseline to statin actions will experience more benefit from statins. These patients have the most severely impaired biological profile on which statins can have an impact on. The only biological feature that is not the “opposite” of statin-induced biological modification is high MMP2. This may be related to the increased inflammation and oxidative stress in hemodialysis patients, as serum MMP-2 is known as a marker of inflammation which is a key mechanistic pathway in the accelerated atherosclerosis in dialysis patients. Although MMP2 concentration decreases with the initiation of hemodialysis, it does not fall to concentrations seen in patients without ESRD [36].

In contrast, phenotype 2 (n = 35), which had spontaneously high levels of HGF and low levels of FAS appeared to experience worse outcome when treated with statins (especially when considering short-term outcomes). As these patients have already high levels of HGF and low levels of FAS, it is likely that statins cannot further modify these biological systems, or could even modify these already extreme concentrations to levels that are deleterious to patients. In other words, statins cannot do good in patients that already have a very favorable profile, but statins might induce adverse events.

The two other clusters (3 and 4) had baseline biomarker profile that were not the inverse of the expected modification induced by statins: cluster 3 had already high levels of HGF, and cluster 4 had already low levels of IL6. In these profiles, no significant response to statin therapy was detected.

This study has its limitations. First, we only had biomarker data available before initiation of rosuvastatin treatment in dialysis patients of the AURORA trial. We therefore were not able to assess the direct impact of statin treatment on concentrations of those mechanistic biomarkers that we have identified in the in-vitro setting. We were nevertheless able to identify four phenotypes based on the baseline biomarker levels that showed differences in association with outcome. Although we were able to identify these phenotypes, we cannot rule out that the patient segmentation might be enhanced by adding further markers that have not been measured in the current study. In our set of DEGs there have been a few other markers of interest that might be worth investigating in future studies like e.g. the connective tissue growth factor (CTGF), the growth differentiation factor 6 (GDF6), or the C-X-C motif chemokine ligand 12 (CXCL12). We are also aware that the cohort of hemodialysis patients is a very special one and we can only speculate at that point in time how well the identified markers will perform in non-dialysis patients with increased risk to develop CVD prone to statin therapy. Most importantly, validating this biological approach in other trials would be key to ascertain external validity.

Logical next steps for future studies are therefore to (i) investigate the potential of the proposed biomarkers in identifying statin responders in patients with advanced CKD not treated with hemodialysis as well as to (ii) measure additional biomarkers being linked to statin mechanism of action in order to evaluate whether the patient clustering can be further enhanced to properly predict statin effect.

## 5. Conclusions

In this translational study we investigated the impact of statin therapy on gene expression changes in renal cells and identified dysregulated mechanisms and markers. A panel of eight markers was able to stratify hemodialysis patients into four phenotypes that showed significant differences in patient outcome. Patients with low HGF and high IL6 appeared to experience better response to statins, while patients with high HGF and low FAS appeared to have worse outcome when treated by statins. Intermediate biological phenotype (with either low IL6 or high FAS) seemed not to be affected by statin treatment.
by both statins in both cell lines. Among the significantly downregulated molecules was the connective tissue growth factor (CTGF), a pro-fibrotic stimulus, which has been shown to be affected by statin treatment, both in renal as well as endothelial cells [25–27]. The Kruppel like factor 2 (KLF2) was among the significantly upregulated molecules by both drugs in both cell lines. KLF2 was reported to be downregulated after high glucose treatment in human endothelial cells and reduced KLF2 expression was associated with aggravated endothelial injury in DN [28,29].

The affected genes in the cell culture setup that were responsible for patient stratification into the four endophenotypes in the AURORA were HGF, IL6, MMP2, FAS, and NTF3. HGF is a cytokine with a number of regulatory functions in the context of immunomodulation, angiogenesis, morphogenesis, as well as tissue repair [30]. HGF was shown to be elevated in hemodialysis patients and a link to cardiovascular risk was postulated although it is still under debate whether HGF increases CVD risk or whether it is actually involved in counterbalancing endothelial damage [30–32]. Burgazli and colleagues reported on the role of statins to modulate HGF-induced proliferation of human umbilical vein endothelial cells [33]. Our data rather indicate that HGF has a protective role with phenotype clusters showing higher values of HGF being linked to improved outcome whereas phenotype cluster 1 with the worst prognosis showed the lowest HGF values. This cluster was also characterized by high MMP2 levels, a molecule involved in fibrotic processes. HGF was previously reported to induce MMP2 levels in human glomeruli next to reducing levels of TGFB1 and type IV collagens thus exerting its anti-fibrotic effects [34]. MMP2 levels are elevated in ESRD with hemodialysis reducing levels of MMP2 [35]. Serum levels of MMP2 are associated with mortality in hemodialysis patients [36]. An inhibitory impact of statin treatment on MMP2 levels was also reported by de Araujo Junior and colleagues as well as by Janardhanan et al. respectively [37,38]. Statins were also shown to reduce levels of pro-inflammatory IL6 in human vascular smooth muscle cells [39] and monocytes from hypercholesterolemic patients [40]. In addition statins reduced IL6 levels in patients on hemodialysis or the pre-dialysis state [41–43]. High IL6 levels in our analysis where characteristic for phenotype 1 holding those patients with the worst prognosis but on the other hand showing the most benefit of statin treatment.

The analysis performed in the AURORA trial strengthens the validity of our results. First, phenotype 1 (n = 125) had the most opposing biological profile (low HGF, high IL6, and NTF3) to the expected biological modifications triggered by statins. Patients of this phenotype had the worst prognosis with respect to overall survival, CV death and also risk of MACE. In addition, this profile appeared to have some benefit from statin treatment when considering all-cause mortality. This suggests that patients in whom the biomarkers are most opposed at baseline to statin actions will experience more benefit from statins. These patients have the most severely impaired biological profile on which statins can have an impact on. The only biological feature that is not the “opposite” of statin-induced biological modification is high MMP2. This may be related to the increased inflammation and oxidative stress in hemodialysis patients, as serum MMP-2 is known as a marker of inflammation which is a key mechanistic pathway in the accelerated atherosclerosis in dialysis patients. Although MMP2 concentration decreases with the initiation of hemodialysis, it does not fall to concentrations seen in patients without ESRD [36].

In contrast, phenotype 2 (n = 35), which had spontaneously high levels of HGF and low levels of FAS appeared to experience worse outcome when treated with statins (especially when considering short-term outcomes). As these patients have already high levels of HGF and low levels of FAS, it is likely that statins cannot further modify these biological systems, or could even modify these already extreme concentrations to levels that are deleterious to patients. In other words, statins cannot do good in patients that already have a very favorable profile, but statins might induce adverse events.

The two other clusters (3 and 4) had baseline biomarker profile that were not the inverse of the expected modification induced by statins: cluster 3 had already high levels of HGF, and cluster 4 had already low levels of IL6. In these profiles, no significant response to statin therapy was detected.

This study has its limitations. First, we only had biomarker data available before initiation of rosuvastatin treatment in dialysis patients of the AURORA trial. We therefore were not able to assess the direct impact of statin treatment on concentrations of those mechanistic biomarkers that we have identified in the in-vitro setting. We were nevertheless able to identify four phenotypes based on the baseline biomarker levels that showed differences in association with outcome. Although we were able to identify these phenotypes, we cannot rule out that the patient segmentation might be enhanced by adding further markers that have not been measured in the current study. In our set of DEGs there have been a few other markers of interest that might be worth investigating in future studies like e.g. the connective tissue growth factor (CTGF), the growth differentiation factor 6 (GDF6), or the C-X-C motif chemokine ligand 12 (CXCL12). We are also aware that the cohort of hemodialysis patients is a very special one and we can only speculate at that point in time how well the identified markers will perform in non-dialysis patients with increased risk to develop CVD prone to statin therapy. Most importantly, validating this biological approach in other trials would be key to ascertain external validity.

Logical next steps for future studies are therefore to (i) investigate the potential of the proposed biomarkers in identifying statin responders in patients with advanced CKD not treated with hemodialysis as well as to (ii) measure additional biomarkers being linked to statin mechanism of action in order to evaluate whether the patient clustering can be further enhanced to properly predict statin effect.

## 5. Conclusions

In this translational study we investigated the impact of statin therapy on gene expression changes in renal cells and identified dysregulated mechanisms and markers. A panel of eight markers was able to stratify hemodialysis patients into four phenotypes that showed significant differences in patient outcome. Patients with low HGF and high IL6 appeared to experience better response to statins, while patients with high HGF and low FAS appeared to have worse outcome when treated by statins. Intermediate biological phenotype (with either low IL6 or high FAS) seemed not to be affected by statin treatment.
## Ethics approval and consent to participate

The study is being conducted in accordance with the ethical principles of the Declaration of Helsinki, the International Conference of Harmonisation Good Clinical Practice guidelines and local regulatory requirements. All randomized patients provided a written informed consent.

## Consent for publication

Not applicable.

## Availability of data and materials

Microarray data are deposited in the Gene Expression Omnibus with accession number GSE208579. All other data from this study are included within the manuscript.

## Funding

This project has received funding from the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement No 115974. The JU receives support from the European Union's Horizon 2020 research and innovation programme and EFPIA and JDRF. Any dissemination of results reflects only the author's view; the JU is not responsible for any use that may be made of the information it contains. This project has also received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 812699 (IMPROVE-PD).

## CRediT authorship contribution statement

**Johannes Leierer:** Data curation, Formal analysis, Visualization, Writing – original draft, Writing – review & editing. **Madonna Salib:** Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. **Michail Evgeniou:** Data curation, Formal analysis, Writing – review & editing. **Patrick Rossignol:** Formal analysis, Investigation, Writing – review & editing. **Ziad A. Massy:** Formal analysis, Investigation, Writing – review & editing. **Klaus Kratochwill:** Formal analysis, Investigation, Writing – review & editing. **Gert Mayer:** Conceptualization, Supervision, Writing – review & editing. **Bengt Fellström:** Formal analysis, Investigation, Writing – review & editing. **Nicolas Girerd:** Conceptualization, Formal analysis, Investigation, Supervision, Writing – review & editing. **Faiez Zannad:** Conceptualization, Formal analysis, Investigation, Supervision, Writing – review & editing. **Paul Perco:** Conceptualization, Formal analysis, Investigation, Supervision, Visualization, Writing – original draft, Writing – review & editing.

## Declaration of competing interest

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Paul Perco reports a relationship with Delta4 GmbH that includes: employment.Klaus Kratochwill reports a relationship with Delta4 GmbH that includes: co-founder.

## Acknowledgements

The authors want to thank all working group members of the EU COST Action Program PerMediK (Personalized Medicine in Chronic Kidney Disease: improved outcome based on Big Data) for fruitful discussions during the regular COST Action project meetings.

## List of abbreviations

<table><tr><td>CV</td><td>cardiovascular</td></tr><tr><td>CVD</td><td>cardiovascular disease</td></tr><tr><td>DEG</td><td>differentially expressed gene</td></tr><tr><td>ESRD</td><td>end-stage renal disease</td></tr><tr><td>FDR</td><td>false discovery rate</td></tr><tr><td>HK2</td><td>proximal tubular human kidney cell</td></tr><tr><td>MACE</td><td>major adverse cardiovascular event</td></tr><tr><td>MC</td><td>mesangial cell</td></tr></table>
## Ethics approval and consent to participate

The study is being conducted in accordance with the ethical principles of the Declaration of Helsinki, the International Conference of Harmonisation Good Clinical Practice guidelines and local regulatory requirements. All randomized patients provided a written informed consent.

## Consent for publication

Not applicable.

## Availability of data and materials

Microarray data are deposited in the Gene Expression Omnibus with accession number GSE208579. All other data from this study are included within the manuscript.

## Funding

This project has received funding from the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement No 115974. The JU receives support from the European Union's Horizon 2020 research and innovation programme and EFPIA and JDRF. Any dissemination of results reflects only the author's view; the JU is not responsible for any use that may be made of the information it contains. This project has also received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 812699 (IMPROVE-PD).

## CRediT authorship contribution statement

**Johannes Leierer:** Data curation, Formal analysis, Visualization, Writing – original draft, Writing – review & editing. **Madonna Salib:** Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. **Michail Evgeniou:** Data curation, Formal analysis, Writing – review & editing. **Patrick Rossignol:** Formal analysis, Investigation, Writing – review & editing. **Ziad A. Massy:** Formal analysis, Investigation, Writing – review & editing. **Klaus Kratochwill:** Formal analysis, Investigation, Writing – review & editing. **Gert Mayer:** Conceptualization, Supervision, Writing – review & editing. **Bengt Fellström:** Formal analysis, Investigation, Writing – review & editing. **Nicolas Girerd:** Conceptualization, Formal analysis, Investigation, Supervision, Writing – review & editing. **Faiez Zannad:** Conceptualization, Formal analysis, Investigation, Supervision, Writing – review & editing. **Paul Perco:** Conceptualization, Formal analysis, Investigation, Supervision, Visualization, Writing – original draft, Writing – review & editing.

## Declaration of competing interest

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Paul Perco reports a relationship with Delta4 GmbH that includes: employment.Klaus Kratochwill reports a relationship with Delta4 GmbH that includes: co-founder.

## Acknowledgements

The authors want to thank all working group members of the EU COST Action Program PerMediK (Personalized Medicine in Chronic Kidney Disease: improved outcome based on Big Data) for fruitful discussions during the regular COST Action project meetings.

## List of abbreviations

<table><tr><td>CV</td><td>cardiovascular</td></tr><tr><td>CVD</td><td>cardiovascular disease</td></tr><tr><td>DEG</td><td>differentially expressed gene</td></tr><tr><td>ESRD</td><td>end-stage renal disease</td></tr><tr><td>FDR</td><td>false discovery rate</td></tr><tr><td>HK2</td><td>proximal tubular human kidney cell</td></tr><tr><td>MACE</td><td>major adverse cardiovascular event</td></tr><tr><td>MC</td><td>mesangial cell</td></tr></table>
References

[1] S.M. Grundy, N.J. Stone, A.L. Bailey, C. Beam, K.K. Birtcher, R.S. Blumenthal, L.T. Braun, S. de Ferranti, J. Faiella-Tommasino, D.E. Forman, R. Goldberg, P. A. Heidenreich, M.A. Hlatky, D.W. Jones, D. Lloyd-Jones, N. Lopez-Pajares, C.E. Ndumele, C.E. Orringer, C.A. Peralta, J.J. Saseen, S.C. Smith, L. Sperling, S. S. Virani, J. Yeboah, AHA/ACC/AACVPR/AAPA/ABC/APCP/ADA/AGS/AphA/ASPC/NLA/PCNA guideline on the management of blood cholesterol, J. Am. Coll. Cardiol. 73 (2018) e285-e350, https://doi.org/10.1016/j.jacc.2018.11.003, 2019.

[2] J. Yoo, J. Jeon, M. Baek, S.O. Song, J. Kim, Impact of statin treatment on cardiovascular risk in patients with type 1 diabetes: a population-based cohort study, J. Transl. Med. 21 (2023) 806, https://doi.org/10.1186/s12967-023-04691-6.

[3] H. Naci, J. Brugs, T. Ades, Comparative tolerability and harms of individual statins: a study-level network meta-analysis of 246 955 participants from 135 randomized, controlled trials, Circ. Cardiovasc. Qual. Outcomes 6 (2013) 390-399, https://doi.org/10.1161/CIRCOUTCOMES.111.000071.

[4] I.A. Mansi, M. Chansard, I. Lingvay, S. Zhang, E.A. Halm, C.A. Alvarez, Association of statin therapy initiation with diabetes progression: a retrospective matched-cohort study, JAMA Intern. Med. 181 (2021) 1562-1574, https://doi.org/10.1001/jamainternmed.2021.5714.

[5] M. Ruscica, N. Ferri, M. Banach, C.R. Sirtori, A. Corsini, Side effects of statins-from pathophysiology and epidemiology to diagnostic and therapeutic implications, Cardiovasc. Res. (2022), https://doi.org/10.1093/cvr/cvac020.

[6] B.C. Fellström, A.G. Jardine, R.E. Schmieder, H. Holdaas, K. Bannister, J. Beutler, D.-W. Chae, A. Chevaile, S.M. Cobbe, C. Grönhagen-Riska, J.J. De Lima, R. Lins, G. Mayer, A.W. McMahon, H.-H. Parving, G. Remuzzi, O. Samuelsson, S. Sonkodi, G. Süleymanlar, D. Tsakiris, V. Tesar, V. Todorov, A. Wiecek, R. P. Wüthrich, M. Gottlow, E. Johnsson, F. Zannad, Rosuvastatin and cardiovascular events in patients undergoing hemodialysis, N. Engl. J. Med. 360 (2009) 1395-1407, https://doi.org/10.1056/NEJMoa0810177.

[7] C. Wanner, V. Krane, W. März, M. Olschewski, J.F.E. Mann, G. Ruf, E. Ritz, German Diabetes and Dialysis Study Investigators, Atorvastatin in patients with type 2 diabetes mellitus undergoing hemodialysis, N. Engl. J. Med. 353 (2005) 238-248, https://doi.org/10.1056/NEJMoa043545.

[8] R.B.R. León-Cachón, J.A. Ascacio-Martínez, M.E. Gamino-Pena, R.M. Cerda-Flores, I. Meester, H.L. Gallardo-Blanco, M. Gómez-Silva, E. Piñeyro-Garza, H. A. Barrera-Saldana, A pharmacogenetic pilot study reveals MTHFR, DRD3, and MDR1 polymorphisms as biomarker candidates for slow atorvastatin metabolizers, BMC Cancer 16 (2016) 74, https://doi.org/10.1186/s12885-016-2062-2.

[9] E. Arrigoni, M. Del Re, L. Fidilio, S. Fogli, R. Danesi, A. Di Paolo, Pharmacogenetic foundations of therapeutic efficacy and adverse events of statins, Int. J. Mol. Sci. 18 (2017), https://doi.org/10.3390/ijms1801043.

[10] M. Elbadawi-Sidhu, R.A. Baillie, H. Zhu, Y.-D.I. Chen, M.O. Goodarzi, J.I. Rotter, R.M. Krauss, O. Fiehn, R. Kaddurah-Daouk, Pharmacometabolomic signature links simvastatin therapy and insulin resistance, Metabolomics 13 (2017), https://doi.org/10.1007/s11306-016-1141-3.

[11] T. Wilmanski, S.A. Kornilov, C. Diener, M. Conomos, J.C. Lovejoy, P. Sebastiani, E.S. Orwoll, L. Hood, N.D. Price, N. Rappaport, A.T. Magis, S.M. Gibbons, Heterogeneity in statin responses explained by variation in the human gut microbiome, medRxiv (2021), https://doi.org/10.1101/2021.12.02.21267193, 2021.12.02.21267193.

[12] B. Perrone, P. Ruffo, G. Augimeri, D. Sisci, M.S. Sinicropi, G. Tripepi, C. Mammi, D. Bonofiglio, F.L. Conforti, Association between HMGCR, CRP, and CETP gene polymorphisms and metabolic/inflammatory serum profile in healthy adolescents, J. Transl. Med. 21 (2023) 718, https://doi.org/10.1186/s12967-023-04571-z.

[13] B. Genser, C. Wanner, W. März, A scoring system for predicting individual treatment effects of statins in type 2 diabetes patients on haemodialysis, Eur. J. Prev. Cardiol. 28 (2021) 838-851, https://doi.org/10.1177/2047487320905721.

[14] M.E. Ritchie, B. Phipson, D. Wu, Y. Hu, C.W. Law, W. Shi, G.K. Smyth, Limma powers differential expression analyses for RNA-sequencing and microarray studies, Nucleic Acids Res. 43 (2015) e47, https://doi.org/10.1093/nar/gkv007.

[15] D.W. Huang, B.T. Sherman, R.A. Lempicki, Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources, Nat. Protoc. 4 (2009) 44-57, https://doi.org/10.1038/nprot.2008.211.

[16] M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, K. Morishima, KEGG: new perspectives on genomes, pathways, diseases and drugs, Nucleic Acids Res. 45 (2017) D353-D361, https://doi.org/10.1093/nar/gkw1092.

[17] B. Fellström, H. Holdaas, A.G. Jardine, H. Rose, R. Schmieder, W. Wilpshaar, F. Zannad, Effect of rosuvastatin on outcomes in chronic haemodialysis patients: baseline data from the AURORA study, kidney blood press, Res. 30 (2007) 314-322, https://doi.org/10.1159/000106803.

[18] B. Fellström, F. Zannad, R. Schmieder, H. Holdaas, A. Jardine, H. Rose, W. Wilpshaar, Effect of rosuvastatin on outcomes in chronic haemodialysis patients – design and rationale of the AURORA study, Curr. Control Trials Cardiovasc. Med. 6 (2005) 9, https://doi.org/10.1186/1468-6708-6-9.

[19] K. Quach, L. Lvtvyn, C. Baigent, J. Bueti, A.X. Garg, C. Hawley, R. Haynes, B. Manns, V. Perkovic, C.G. Rabbat, R. Wald, M. Walsh, The safety and efficacy of mineralocorticoid receptor antagonists in patients who require dialysis: a systematic review and meta-analysis, Am. J. Kidney Dis. 68 (2016) 591-598, https://doi.org/10.1053/j.ajkd.2016.04.011.

[20] J. Deng, G. Wu, C. Yang, Y. Li, Q. Jing, Y. Han, Rosuvastatin attenuates contrast-induced nephropathy through modulation of nitric oxide, inflammatory responses, oxidative stress and apoptosis in diabetic male rats, J. Transl. Med. 13 (2015) 53, https://doi.org/10.1186/s12967-015-0416-1.

[21] M. Uhlén, L. Fagerberg, B.M. Hallström, C. Lindskog, P. Oksvold, A. Mardinoglu, Å. Sivertsson, C. Kampf, E. Sjöstedt, A. Asplund, I. Olsson, K. Edlund, E. Lundberg, S. Navani, C.A.-K. Szigyarto, J. Odeber, D. Djureinovic, J.O. Takanen, S. Hober, T. Alm, P.-H. Edqvist, H. Berling, H. Tegel, J. Mulder, J. Rockberg, P. Nilsson, J.M. Schwenk, M. Hamsten, K. von Feilitzen, M. Forsberg, L. Persson, F. Johansson, M. Zwahlen, G. von Heijne, J. Nielsen, F. Pontén, Proteomics. Tissue-based map of the human proteome, Science 347 (2015) 1260419, https://doi.org/10.1126/science.1260419.

[22] A.J. Tremblay, B. Lamarche, V. Lemelin, L. Hoos, S. Benjannet, N.G. Seidah, H.R. Davis, P. Couture, Atorvastatin increases intestinal expression of NPC1L1 in hyperlipidemic men, J. Lipid Res. 52 (2011) 558-565, https://doi.org/10.1194/jlr.M011080.

[23] C. Stormo, M.K. Kringen, R.M. Grimholt, J.P. Berg, A.P. Piehler, A novel 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) splice variant with an alternative exon 1 potentially encoding an extended N-terminus, BMC Mol. Biol. 13 (2012) 29, https://doi.org/10.1186/1471-2199-13-29.

[24] S. Gill, J. Stevenson, I. Kristiana, A.J. Brown, Cholesterol-dependent degradation of squalene monooxygenase, a control point in cholesterol synthesis beyond HMG-CoA reductase, Cell Metabol. 13 (2011) 260-273, https://doi.org/10.1016/j.cmet.2011.01.015.

[25] E.H. Bae, I.J. Kim, J.W. Park, S.K. Ma, J.U. Lee, S.W. Kim, Renoprotective effect of rosuvastatin in DOCA-salt hypertensive rats, Nephrol. Dial. Transplant. 25 (2010) 1051-1059, https://doi.org/10.1093/ndt/gfp604.

[26] B. Zhu, A.-Q. Ma, L. Yang, X.-M. Dang, Atorvastatin attenuates bleomycin-induced pulmonary fibrosis via suppressing iNOS expression and the CTGF (CCN2)/ERK signaling pathway, Int. J. Mol. Sci. 14 (2013) 24476-24491, https://doi.org/10.3390/ijms141224476.

[27] S. Muehlich, N. Schneider, F. Hinkmann, C.D. Garlich, M. Goppelt-Struebe, Induction of connective tissue growth factor (CTGF) in human endothelial cells by lysophosphatidic acid, sphingosine-1-phosphate, and platelets, Atherosclerosis 175 (2004) 261-268, https://doi.org/10.1016/j.atherosclerosis.2004.04.011.

[28] F. Zhong, H. Chen, C. Wei, W. Zhang, Z. Li, M.K. Jain, P.Y. Chuang, H. Chen, Y. Wang, S.K. Mallipattu, J.C. He, Reduced Krüppel-like factor 2 expression may aggravate the endothelial injury of diabetic nephropathy, Kidney Int. 87 (2015) 382-395, https://doi.org/10.1038/ki.2014.286.

[29] Y.-S. Liu, D.-L. Xu, Z.-W. Huang, L. Hao, X. Wang, Q.-H. Lu, Atorvastatin counteracts high glucose-induced Krüppel-like factor 2 suppression in human umbilical vein endothelial cells, Postgrad. Med. 127 (2015) 446-454, https://doi.org/10.1080/00325481.2015.1039451.

[30] C. Libetta, P. Esposito, C. Martinelli, F. Grosjean, M. Gregorini, T. Rampino, A. Dal Canton, Hepatocyte growth factor (HGF) and hemodialysis: physiopathology and clinical implications, Clin. Exp. Nephrol. 20 (2016) 371-378, https://doi.org/10.1007/s10157-015-1211-2.

[31] T. Rampino, C. Libetta, W. De Simone, A. Ranghino, M. Soccio, M. Gregorini, P. Guallini, L. Tamagnone, A. Dal Canton, Hemodialysis stimulates hepatocyte growth factor release, Kidney Int. 53 (1998) 1382-1388, https://doi.org/10.1046/j.1523-1755.1998.00870.x.

[32] K. Janda, M. Krzanowski, P. Dumnicka, B. Kuśnierz-Cabala, D. Sorysz, W. Sulowicz, Hepatocyte growth factor as a long-term predictor for total and cardiovascular mortality in patients on peritoneal dialysis, Pol. Arch. Med. Wlew. 123 (2013) 453-459, https://doi.org/10.20452/pawm.1883.

[33] K.M. Burgazli, K.L. Bui, M. Mericililer, A.T. Albayrak, M. Parahuleva, A. Erdogan, The effects of different types of statins on proliferation and migration of HGF-induced Human Umbilical Vein Endothelial Cells (HUVECs), Eur. Rev. Med. Pharmacol. Sci. 17 (2013) 2874-2883.
References

[1] S.M. Grundy, N.J. Stone, A.L. Bailey, C. Beam, K.K. Birtcher, R.S. Blumenthal, L.T. Braun, S. de Ferranti, J. Faiella-Tommasino, D.E. Forman, R. Goldberg, P. A. Heidenreich, M.A. Hlatky, D.W. Jones, D. Lloyd-Jones, N. Lopez-Pajares, C.E. Ndumele, C.E. Orringer, C.A. Peralta, J.J. Saseen, S.C. Smith, L. Sperling, S. S. Virani, J. Yeboah, AHA/ACC/AACVPRA/AAPC/ABC/APCM/ADA/AGS/AphA/ASPC/NLA/PCNA guideline on the management of blood cholesterol, J. Am. Coll. Cardiol. 73 (2018) e285-e350, https://doi.org/10.1016/j.jacc.2018.11.003, 2019.

[2] J. Yoo, J. Jeon, M. Baek, S.O. Song, J. Kim, Impact of statin treatment on cardiovascular risk in patients with type 1 diabetes: a population-based cohort study, J. Transl. Med. 21 (2023) 806, https://doi.org/10.1186/s12967-023-04691-6.

[3] H. Naci, J. Brugs, T. Ades, Comparative tolerability and harms of individual statins: a study-level network meta-analysis of 246 955 participants from 135 randomized, controlled trials, Circ. Cardiovasc. Qual. Outcomes 6 (2013) 390-399, https://doi.org/10.1161/CIRCOUTCOMES.111.000071.

[4] I.A. Mansi, M. Chansard, I. Lingvay, S. Zhang, E.A. Halm, C.A. Alvarez, Association of statin therapy initiation with diabetes progression: a retrospective matched-cohort study, JAMA Intern. Med. 181 (2021) 1562-1574, https://doi.org/10.1001/jamainternmed.2021.5714.

[5] M. Ruscica, N. Ferri, M. Banach, C.R. Sirtori, A. Corsini, Side effects of statins-from pathophysiology and epidemiology to diagnostic and therapeutic implications, Cardiovasc. Res. (2022), https://doi.org/10.1093/cvr/cvac020.

[6] B.C. Fellström, A.G. Jardine, R.E. Schmieder, H. Holdaas, K. Bannister, J. Beutler, D.-W. Chae, A. Chevaile, S.M. Cobbe, C. Grönhagen-Riska, J.J. De Lima, R. Lins, G. Mayer, A.W. McMahon, H.-H. Parving, G. Remuzzi, O. Samuelsson, S. Sonkodi, G. Süleymanlar, D. Tsakiris, V. Tesar, V. Todorov, A. Wiecek, R. P. Wüthrich, M. Gottlow, E. Johnsson, F. Zannad, Rosuvastatin and cardiovascular events in patients undergoing hemodialysis, N. Engl. J. Med. 360 (2009) 1395-1407, https://doi.org/10.1056/NEJMoa0810177.

[7] C. Wanner, V. Krane, W. März, M. Olschewski, J.F.E. Mann, G. Ruf, E. Ritz, German Diabetes and Dialysis Study Investigators, Atorvastatin in patients with type 2 diabetes mellitus undergoing hemodialysis, N. Engl. J. Med. 353 (2005) 238-248, https://doi.org/10.1056/NEJMoa043545.

[8] R.B.R. León-Cachón, J.A. Ascacio-Martínez, M.E. Gamino-Pena, R.M. Cerda-Flores, I. Meester, H.L. Gallardo-Blanco, M. Gómez-Silva, E. Piñeyro-Garza, H. A. Barrera-Saldana, A pharmacogenetic pilot study reveals MTHFR, DRD3, and MDR1 polymorphisms as biomarker candidates for slow atorvastatin metabolizers, BMC Cancer 16 (2016) 74, https://doi.org/10.1186/s12885-016-2062-2.

[9] E. Arrigoni, M. Del Re, L. Fidilio, S. Fogli, R. Danesi, A. Di Paolo, Pharmacogenetic foundations of therapeutic efficacy and adverse events of statins, Int. J. Mol. Sci. 18 (2017), https://doi.org/10.3390/ijms1801043.

[10] M. Elbadawi-Sidhu, R.A. Baillie, H. Zhu, Y.-D.I. Chen, M.O. Goodarzi, J.I. Rotter, R.M. Krauss, O. Fiehn, R. Kaddurah-Daouk, Pharmacometabolomic signature links simvastatin therapy and insulin resistance, Metabolomics 13 (2017), https://doi.org/10.1007/s11306-016-1141-3.

[11] T. Wilmanski, S.A. Kornilov, C. Diener, M. Conomos, J.C. Lovejoy, P. Sebastiani, E.S. Orwoll, L. Hood, N.D. Price, N. Rappaport, A.T. Magis, S.M. Gibbons, Heterogeneity in statin responses explained by variation in the human gut microbiome, medRxiv (2021), https://doi.org/10.1101/2021.12.02.21267193, 2021.12.02.21267193.

[12] B. Perrone, P. Ruffo, G. Augimeri, D. Sisci, M.S. Sinicropi, G. Tripepi, C. Mammi, D. Bonofiglio, F.L. Conforti, Association between HMGCR, CRP, and CETP gene polymorphisms and metabolic/inflammatory serum profile in healthy adolescents, J. Transl. Med. 21 (2023) 718, https://doi.org/10.1186/s12967-023-04571-z.

[13] B. Genser, C. Wanner, W. März, A scoring system for predicting individual treatment effects of statins in type 2 diabetes patients on haemodialysis, Eur. J. Prev. Cardiol. 28 (2021) 838-851, https://doi.org/10.1177/2047487320905721.

[14] M.E. Ritchie, B. Phipson, D. Wu, Y. Hu, C.W. Law, W. Shi, G.K. Smyth, Limma powers differential expression analyses for RNA-sequencing and microarray studies, Nucleic Acids Res. 43 (2015) e47, https://doi.org/10.1093/nar/gkv007.

[15] D.W. Huang, B.T. Sherman, R.A. Lempicki, Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources, Nat. Protoc. 4 (2009) 44-57, https://doi.org/10.1038/nprot.2008.211.

[16] M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, K. Morishima, KEGG: new perspectives on genomes, pathways, diseases and drugs, Nucleic Acids Res. 45 (2017) D353-D361, https://doi.org/10.1093/nar/gkw1092.

[17] B. Fellström, H. Holdaas, A.G. Jardine, H. Rose, R. Schmieder, W. Wilpshaar, F. Zannad, Effect of rosuvastatin on outcomes in chronic haemodialysis patients: baseline data from the AURORA study, kidney blood press, Res. 30 (2007) 314-322, https://doi.org/10.1159/000106803.

[18] B. Fellström, F. Zannad, R. Schmieder, H. Holdaas, A. Jardine, H. Rose, W. Wilpshaar, Effect of rosuvastatin on outcomes in chronic haemodialysis patients – design and rationale of the AURORA study, Curr. Control Trials Cardiovasc. Med. 6 (2005) 9, https://doi.org/10.1186/1468-6708-6-9.

[19] K. Quach, L. Lvtvyn, C. Baigent, J. Bueti, A.X. Garg, C. Hawley, R. Haynes, B. Manns, V. Perkovic, C.G. Rabbat, R. Wald, M. Walsh, The safety and efficacy of mineralocorticoid receptor antagonists in patients who require dialysis: a systematic review and meta-analysis, Am. J. Kidney Dis. 68 (2016) 591-598, https://doi.org/10.1053/j.ajkd.2016.04.011.

[20] J. Deng, G. Wu, C. Yang, Y. Li, Q. Jing, Y. Han, Rosuvastatin attenuates contrast-induced nephropathy through modulation of nitric oxide, inflammatory responses, oxidative stress and apoptosis in diabetic male rats, J. Transl. Med. 13 (2015) 53, https://doi.org/10.1186/s12967-015-0416-1.

[21] M. Uhlén, L. Fagerberg, B.M. Hallström, C. Lindskog, P. Oksvold, A. Mardinoglu, Å. Sivertsson, C. Kampf, E. Sjöstedt, A. Asplund, I. Olsson, K. Edlund, E. Lundberg, S. Navani, C.A.-K. Szigyarto, J. Odeber, D. Djureinovic, J.O. Takanen, S. Hober, T. Alm, P.-H. Edqvist, H. Berling, H. Tegel, J. Mulder, J. Rockberg, P. Nilsson, J.M. Schwenk, M. Hamsten, K. von Feilitzen, M. Forsberg, L. Persson, F. Johansson, M. Zwahlen, G. von Heijne, J. Nielsen, F. Pontén, Proteomics. Tissue-based map of the human proteome, Science 347 (2015) 1260419, https://doi.org/10.1126/science.1260419.

[22] A.J. Tremblay, B. Lamarche, V. Lemelin, L. Hoos, S. Benjannet, N.G. Seidah, H.R. Davis, P. Couture, Atorvastatin increases intestinal expression of NPC1L1 in hyperlipidemic men, J. Lipid Res. 52 (2011) 558-565, https://doi.org/10.1194/jlr.M011080.

[23] C. Stormo, M.K. Kringen, R.M. Grimholt, J.P. Berg, A.P. Piehler, A novel 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) splice variant with an alternative exon 1 potentially encoding an extended N-terminus, BMC Mol. Biol. 13 (2012) 29, https://doi.org/10.1186/1471-2199-13-29.

[24] S. Gill, J. Stevenson, I. Kristiana, A.J. Brown, Cholesterol-dependent degradation of squalene monooxygenase, a control point in cholesterol synthesis beyond HMG-CoA reductase, Cell Metabol. 13 (2011) 260-273, https://doi.org/10.1016/j.cmet.2011.01.015.

[25] E.H. Bae, I.J. Kim, J.W. Park, S.K. Ma, J.U. Lee, S.W. Kim, Renoprotective effect of rosuvastatin in DOCA-salt hypertensive rats, Nephrol. Dial. Transplant. 25 (2010) 1051-1059, https://doi.org/10.1093/ndt/gfp604.

[26] B. Zhu, A.-Q. Ma, L. Yang, X.-M. Dang, Atorvastatin attenuates bleomycin-induced pulmonary fibrosis via suppressing iNOS expression and the CTGF (CCN2)/ERK signaling pathway, Int. J. Mol. Sci. 14 (2013) 24476-24491, https://doi.org/10.3390/ijms141224476.

[27] S. Muehlich, N. Schneider, F. Hinkmann, C.D. Garlich, M. Goppelt-Struebe, Induction of connective tissue growth factor (CTGF) in human endothelial cells by lysophosphatidic acid, sphingosine-1-phosphate, and platelets, Atherosclerosis 175 (2004) 261-268, https://doi.org/10.1016/j.atherosclerosis.2004.04.011.

[28] F. Zhong, H. Chen, C. Wei, W. Zhang, Z. Li, M.K. Jain, P.Y. Chuang, H. Chen, Y. Wang, S.K. Mallipattu, J.C. He, Reduced Krüppel-like factor 2 expression may aggravate the endothelial injury of diabetic nephropathy, Kidney Int. 87 (2015) 382-395, https://doi.org/10.1038/ki.2014.286.

[29] Y.-S. Liu, D.-L. Xu, Z.-W. Huang, L. Hao, X. Wang, Q.-H. Lu, Atorvastatin counteracts high glucose-induced Krüppel-like factor 2 suppression in human umbilical vein endothelial cells, Postgrad. Med. 127 (2015) 446-454, https://doi.org/10.1080/00325481.2015.1039451.

[30] C. Libetta, P. Esposito, C. Martinelli, F. Grosjean, M. Gregorini, T. Rampino, A. Dal Canton, Hepatocyte growth factor (HGF) and hemodialysis: physiopathology and clinical implications, Clin. Exp. Nephrol. 20 (2016) 371-378, https://doi.org/10.1007/s10157-015-1211-2.

[31] T. Rampino, C. Libetta, W. De Simone, A. Ranghino, M. Soccio, M. Gregorini, P. Guallini, L. Tamagnone, A. Dal Canton, Hemodialysis stimulates hepatocyte growth factor release, Kidney Int. 53 (1998) 1382-1388, https://doi.org/10.1046/j.1523-1755.1998.00870.x.

[32] K. Janda, M. Krzanowski, P. Dumnicka, B. Kuśnierz-Cabala, D. Sorysz, W. Sulowicz, Hepatocyte growth factor as a long-term predictor for total and cardiovascular mortality in patients on peritoneal dialysis, Pol. Arch. Med. Wlew. 123 (2013) 453-459, https://doi.org/10.20452/pawm.1883.

[33] K.M. Burgazli, K.L. Bui, M. Mericililer, A.T. Albayrak, M. Parahuleva, A. Erdogan, The effects of different types of statins on proliferation and migration of HGF-induced Human Umbilical Vein Endothelial Cells (HUVECs), Eur. Rev. Med. Pharmacol. Sci. 17 (2013) 2874-2883.
[34] C. Esposito, B. Parrilla, A. De Mauri, F. Cornacchia, G. Fasoli, A. Foschi, T. Mazzullo, A. Plati, R. Scudellaro, A. Dal Canton, Hepatocyte growth factor (HGF) modulates matrix turnover in human glomeruli, Kidney Int. 67 (2005) 2143–2150, https://doi.org/10.1111/j.1523-1755.2005.00319.x.

[35] B.P. Marson, R. Lacchini, V. Belo, S. Dickel, B.P. da Costa, C.E. Poli de Figueiredo, J.E. Tanus-Santos, Matrix metalloproteinase (MMP)-2 genetic variants modify the circulating MMP-2 levels in end-stage kidney disease, Am. J. Nephrol. 35 (2012) 209–215, https://doi.org/10.1159/000336108.

[36] K.-C. Hsiao, J.-P. Tsai, S.-F. Yang, W.-C. Lee, J.-Y. Huang, S.-C. Chang, C.-S. Hso, H.-R. Chang, MMP-2 serum concentrations predict mortality in hemodialysis patients: a 5-year cohort study, Clin. Chim. Acta 452 (2016) 161–166, https://doi.org/10.1016/j.cca.2015.11.019.

[37] R.F. de Araújo Júnior, T.O. Souza, L.M. de Moura, K.P. Torres, L.B. de Souza, M. do S.C.F. Alves, H.O. Rocha, A.A. de Araújo, Atorvastatin decreases bone loss, inflammation and oxidative stress in experimental periodontitis, PLoS One 8 (2013) e75322, https://doi.org/10.1371/journal.pone.0075322.

[38] R. Janardhanan, B. Yang, P. Vohra, B. Roy, S. Withers, S. Bhattacharya, J. Mandrekar, H. Kong, E.B. Leof, D. Mukhopadhyay, S. Misra, Simvastatin reduces venous stenosis formation in a murine hemodialysis vascular access model, Kidney Int. 84 (2013) 338–352, https://doi.org/10.1038/ki.2013.112.

[39] H. Loppnow, L. Zhang, M. Buerke, M. Lautenschläger, L. Chen, A. Frister, A. Schlitt, T. Luther, N. Song, B. Hofmann, S. Rose-John, R.-E. Silber, U. Müller-Werdan, K. Werdan, Statins potently reduce the cytokine-mediated IL-6 release in SMC/MNC cocultures, J. Cell Mol. Med. 15 (2011) 994–1004, https://doi.org/10.1111/j.1582-4934.2010.01036.x.

[40] A. Rezaie-Majd, T. Maca, R.A. Bucek, P. Valent, M.R. Müller, P. Husslein, A. Kashanipour, E. Minar, M. Baghestanian, Simvastatin reduces expression of cytokines interleukin-6, interleukin-8, and monocyte chemoattractant protein-1 in circulating monocytes from hypercholesterolemic patients, Arterioscler. Thromb. Vasc. Biol. 22 (2002) 1194–1199, https://doi.org/10.1161/01.at.0000022694.16328.cc.

[41] J. Tian, X. Hou, L. Hu, T. Chen, K. Wu, C. Cai, X. Bai, Efficacy comparison of atorvastatin versus rosuvastatin on blood lipid and microinflammatory state in maintenance hemodialysis patients, Ren. Fail. 39 (2017) 153–158, https://doi.org/10.1080/0886022X.2016.1256309.

[42] V. Panichi, S. Paoletti, E. Mantuano, G. Manca-Rizza, C. Filippi, S. Santi, D. Taccola, C. Donadio, G. Tramonti, M. Innocenti, G. Casto, C. Consani, G. Sbragia, F. Franzoni, F. Galetta, E. Panicucci, G. Barsotti, In vivo and in vitro effects of simvastatin on inflammatory markers in pre-dialysis patients, Nephrol. Dial. Transplant. 21 (2006) 337–344, https://doi.org/10.1093/ndt/gfi224.

[43] G. Tsirpanlis, F. Boufidou, S. Manganas, K. Chantzis, A. Bleta, K. Stamatelou, E. Psimenou, C. Nicolaou, Treatment with fluvastatin rapidly modulates, via different pathways, and in dependence on the baseline level, inflammation in hemodialysis patients, Blood Purif. 22 (2004) 518–524, https://doi.org/10.1159/000082166.
[34] C. Esposito, B. Parrilla, A. De Mauri, F. Cornacchia, G. Fasoli, A. Foschi, T. Mazzullo, A. Plati, R. Scudellaro, A. Dal Canton, Hepatocyte growth factor (HGF) modulates matrix turnover in human glomeruli, Kidney Int. 67 (2005) 2143–2150, https://doi.org/10.1111/j.1523-1755.2005.00319.x.

[35] B.P. Marson, R. Lacchini, V. Belo, S. Dickel, B.P. da Costa, C.E. Poli de Figueiredo, J.E. Tanus-Santos, Matrix metalloproteinase (MMP)-2 genetic variants modify the circulating MMP-2 levels in end-stage kidney disease, Am. J. Nephrol. 35 (2012) 209–215, https://doi.org/10.1159/000336108.

[36] K.-C. Hsiao, J.-P. Tsai, S.-F. Yang, W.-C. Lee, J.-Y. Huang, S.-C. Chang, C.-S. Hso, H.-R. Chang, MMP-2 serum concentrations predict mortality in hemodialysis patients: a 5-year cohort study, Clin. Chim. Acta 452 (2016) 161–166, https://doi.org/10.1016/j.cca.2015.11.019.

[37] R.F. de Araújo Júnior, T.O. Souza, L.M. de Moura, K.P. Torres, L.B. de Souza, M. do S.C.F. Alves, H.O. Rocha, A.A. de Araújo, Atorvastatin decreases bone loss, inflammation and oxidative stress in experimental periodontitis, PLoS One 8 (2013) e75322, https://doi.org/10.1371/journal.pone.0075322.

[38] R. Janardhanan, B. Yang, P. Vohra, B. Roy, S. Withers, S. Bhattacharya, J. Mandrekar, H. Kong, E.B. Leof, D. Mukhopadhyay, S. Misra, Simvastatin reduces venous stenosis formation in a murine hemodialysis vascular access model, Kidney Int. 84 (2013) 338–352, https://doi.org/10.1038/ki.2013.112.

[39] H. Loppnow, L. Zhang, M. Buerke, M. Lautenschläger, L. Chen, A. Frister, A. Schlitt, T. Luther, N. Song, B. Hofmann, S. Rose-John, R.-E. Silber, U. Müller-Werdan, K. Werdan, Statins potently reduce the cytokine-mediated IL-6 release in SMC/MNC cocultures, J. Cell Mol. Med. 15 (2011) 994–1004, https://doi.org/10.1111/j.1582-4934.2010.01036.x.

[40] A. Rezaie-Majd, T. Maca, R.A. Bucek, P. Valent, M.R. Müller, P. Husslein, A. Kashanipour, E. Minar, M. Baghestanian, Simvastatin reduces expression of cytokines interleukin-6, interleukin-8, and monocyte chemoattractant protein-1 in circulating monocytes from hypercholesterolemic patients, Arterioscler. Thromb. Vasc. Biol. 22 (2002) 1194–1199, https://doi.org/10.1161/01.at.0000022694.16328.cc.

[41] J. Tian, X. Hou, L. Hu, T. Chen, K. Wu, C. Cai, X. Bai, Efficacy comparison of atorvastatin versus rosuvastatin on blood lipid and microinflammatory state in maintenance hemodialysis patients, Ren. Fail. 39 (2017) 153–158, https://doi.org/10.1080/0886022X.2016.1256309.

[42] V. Panichi, S. Paoletti, E. Mantuano, G. Manca-Rizza, C. Filippi, S. Santi, D. Taccola, C. Donadio, G. Tramonti, M. Innocenti, G. Casto, C. Consani, G. Sbragia, F. Franzoni, F. Galetta, E. Panicucci, G. Barsotti, In vivo and in vitro effects of simvastatin on inflammatory markers in pre-dialysis patients, Nephrol. Dial. Transplant. 21 (2006) 337–344, https://doi.org/10.1093/ndt/gfi224.

[43] G. Tsirpanlis, F. Boufidou, S. Manganas, K. Chantzis, A. Bleta, K. Stamatelou, E. Psimenou, C. Nicolaou, Treatment with fluvastatin rapidly modulates, via different pathways, and in dependence on the baseline level, inflammation in hemodialysis patients, Blood Purif. 22 (2004) 518–524, https://doi.org/10.1159/000082166.