Next Article in Journal
Health Impact of Silver Nanoparticles: A Review of the Biodistribution and Toxicity Following Various Routes of Exposure
Next Article in Special Issue
Antibodies, Nanobodies, or Aptamers—Which Is Best for Deciphering the Proteomes of Non-Model Species?
Previous Article in Journal
Genetic and Clinical Heterogeneity in Thirteen New Cases with Aceruloplasminemia. Atypical Anemia as a Clue for an Early Diagnosis
Previous Article in Special Issue
Gel-Based Proteomics of Clinical Samples Identifies Potential Serological Biomarkers for Early Detection of Colorectal Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomic Analysis of Renal Biomarkers of Kidney Allograft Fibrosis—A Study in Renal Transplant Patients

1
Department of Nephrology, Odense University Hospital, DK-5000 Odense, Denmark
2
Department of Epidemiology, Biostatistics and Biodemography, University of Southern Denmark, DK-5000 Odense, Denmark
3
Department of Clinical Genetics, Odense University Hospital, DK-5000 Odense, Denmark
4
Department of Pathology, Odense University Hospital, DK-5000 Odense, Denmark
5
Department of Clinical Biochemistry and Pharmacology, Centre for Clinical Proteomics, Odense University Hospital, DK-5000 Odense, Denmark
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(7), 2371; https://doi.org/10.3390/ijms21072371
Submission received: 18 February 2020 / Revised: 12 March 2020 / Accepted: 25 March 2020 / Published: 30 March 2020

Abstract

:
Renal transplantation is the preferred treatment of end stage renal disease, but allograft survival is limited by the development of interstitial fibrosis and tubular atrophy in response to various stimuli. Much effort has been put into identifying new protein markers of fibrosis to support the diagnosis. In the present work, we performed an in-depth quantitative proteomics analysis of allograft biopsies from 31 prevalent renal transplant patients and correlated the quantified proteins with the volume fraction of fibrosis as determined by a morphometric method. Linear regression analysis identified four proteins that were highly associated with the degree of interstitial fibrosis, namely Coagulation Factor XIII A chain (estimate 18.7, adjusted p < 0.03), Uridine Phosphorylase 1 (estimate 19.4, adjusted p < 0.001), Actin-related protein 2/3 subunit 2 (estimate 34.2, adjusted p < 0.05) and Cytochrome C Oxidase Assembly Factor 6 homolog (estimate −44.9, adjusted p < 0.002), even after multiple testing. Proteins that were negatively associated with fibrosis (p < 0.005) were primarily related to normal metabolic processes and respiration, whereas proteins that were positively associated with fibrosis (p < 0.005) were involved in catabolic processes, cytoskeleton organization and the immune response. The identified proteins may be candidates for further validation with regards to renal fibrosis. The results support the notion that cytoskeleton organization and immune responses are prevalent processes in renal allograft fibrosis.

1. Introduction

Kidney transplantation is the preferred treatment of end stage renal disease, reducing mortality when compared to any type of dialysis [1], and increasing quality of life [2]. The introduction of modern induction therapy, as well as highly effective immunosuppressive regimens, has reduced graft loss due to acute rejection [3]. Over time, however, kidney graft function inevitably deteriorates. Improving long-term graft survival remains a key issue in renal transplantation.
The histological finding of chronic allograft lesions with no known etiology is referred to as IFTA (interstitial fibrosis and tubular atrophy) [4]. IFTA is often accompanied by deteriorating renal function, and the presence of IFTA predicts an adverse renal outcome [5]. The development of kidney fibrosis is a multifactorial process including inflammation and ischemia, which ultimately leads to the deposition of extracellular matrix proteins [6].
The causes of chronic allograft failure are diverse and might include acute or chronic rejection, recurrent disease, drug toxicity or infection [7]. Irrespective of the primary insult, IFTA is characterized by largely irreversible damage. Hence, IFTA is important to distinguish from reversible conditions that need specific interventions. Also, insight into the pathogenesis of IFTA might provide the opportunity for early, preventative therapy and guide the development of pharmacological interventions.
Currently, the diagnosis of IFTA is made by kidney biopsy. Although a minor procedure, a biopsy introduces a risk of a gross hematuria of 3%–4%, and a risk of a perirenal hematoma of 2%–4% [8], which necessitates post-procedural observation. Furthermore, there is a risk of sampling variability [9]. Much effort has been devoted to identifying biomarkers to guide diagnostics in acute or chronic renal failure, but so far, no biomarkers are in routine clinical use. In tissue, proteins are prevalent in both the intracellular and extracellular compartments, and have diverse functions, serving as structural components and as mediators of a wide range of biological functions (i.e., as transcription factors, hormones, antibodies and enzymes) [10]. The emergence of large-scale proteomic approaches provides a unique opportunity to gain insight into proteins involved in specific diseases. The use of proteomics is, however, critically dependent on interpreting results in context to make biological sense of the large amounts of information provided [11].
Only a few studies have investigated the proteome in chronic allograft dysfunction. Early discovery-driven proteogenomic analyses of biopsies with varying degrees of IFTA revealed several differences in the proteome between fibrotic and non-fibrotic kidneys [12]. The objective of the current study was to identify proteins that correlate with the degree of fibrosis in renal biopsies from stable kidney transplant patients. The purpose was to gain insight into the mechanisms of fibrogenesis and to identify biomarker candidates to guide the diagnosis of IFTA. To address this, a set of biopsies from 31 prevalent renal transplant patients were evaluated by explorative in-depth proteomics based on nano-liquid chromatography combined with orbitrap mass spectrometry analysis and 10-plex tandem mass tags.

2. Results

Biopsies from 31 individuals were examined and their clinical characteristics are summarized in Table 1. Twenty two (71%) participants were male and the median age was 60 (ranging from 24–73) years. The median time since transplantation was 3 years, but ranged from 0.3 to 12.8 years. All patients were treated with a calcineurin inhibitor. The most prevalent immunosuppressive regimen was tacrolimus in combination with mycophenolate mofetil, but four subjects (13%) were additionally treated with a small dosage of prednisolone. The majority of the included subjects had hypertension (94%). Two patients (6%) had diabetes at the time of inclusion. Three patients had a previous rejection occurring between 2 and 3.5 years prior to inclusion in the study. In all cases, the rejection episode was treated with prednisolone, and renal function returned to pre-rejection levels. The fibrosis quantifications in these subjects at inclusion were below average, at 28%, 29% and 32%, respectively.

2.1. Characterization of Fibrosis in Kidney Allograft Biopsies

The distribution of interstitial fibrosis was evaluated by point counting (Figure 1a) using the Banff Lesion Score method, which relies on the assessment of the presence and degree of different histopathological changes in different compartments of the renal biopsy, such as fibrosis, interstitial inflammation, and mild tubulitis, important for the diagnosis of graft rejection. The median number of points counted per biopsy was 284 (IQR 252–329). The median extent of fibrosis, estimated as the volume fraction, was 33% (IQR 30%–41%). Interstitial inflammation was present in 16 biopsies (n = 12 (29%) for Banff I1 score, n = 3 (10%) for Banff I2 score and n = 1 (3%) for Banff I3 score). Mild tubulitis (Banff T1 score) was present in three biopsies (10%) (see Supplementary Table S1 for a complete summary of the Banff scores). No biopsies gave the suspicion of acute or chronic rejection.
The association between renal function and the extent of fibrosis is shown in Figure 1b. The variation in the extent of fibrosis increased with decreasing renal function. There was a weak, but not significant, correlation between renal function and the volume fraction of interstitial fibrosis (Spearman’s ρ = −0.23, p = 0.21).

2.2. Kidney Allograft Proteomics

Nano-LC-MSMS analysis analyzed 4717 proteins across the analyzed biopsies, whereof 1973 proteins were quantified in all patient biopsies (Supplementary Table S2). Having missing values is a well-known phenomenon in quantitative mass spectrometry-based proteomics when using data-dependent acquisition, and we found that the number of missing values per quantified protein acceptable for downstream statistical analyses was 11, leaving a dataset of 2687 proteins analyzed in 31 kidney biopsies for further evaluation.

2.3. Proteins Correlated with the Degree of Fibrosis

The multiple linear regression identified 26 proteins with an adjusted p-value of less than 0.5 (Table 2). Interestingly, four proteins were highly significantly associated with the degree of fibrosis after correction for multiple testing. Of these, three were positively associated: Coagulation Factor XIII A chain (estimate = 18.7, adjusted p < 0.03, Figure 2a), Uridine Phosphorylase 1 (estimate = 19.4, adjusted p < 0.001, Figure 2d) and Actin-related protein 2/3 subunit 2 (estimate = 34.2, adjusted p < 0.05, Figure 2b). Cytochrome C Oxidase Assembly Factor 6 homolog was negatively associated (estimate = −44.9, adjusted p < 0.002, Figure 2c).
Then we performed a LASSO analysis in order to identify the most predictive proteins for fibrosis (Supplementary Table S3). The nine proteins selected by the LASSO in at least 10% of the subsamples are listed in Table 3.
Interestingly, Cytochrome C Oxidase Assembly Factor 6 homolog (selection probability 56.7%), Coagulation Factor XIII A chain (selection probability 21.7%) and Actin-related protein 2/3 subunit 2 (selection probability 10.7%), which all strongly correlated with the degree of fibroses, were also highly predictive of fibrosis. The fourth protein that strongly correlated with fibrosis, Uridine Phosphorylase 1, was not identified in the analysis, most probably because this protein was not detected in all biopsies, and only proteins with no missing values across patients were included in the LASSO analysis.
In order to identify molecular pathways and biological processes, we evaluated 133 proteins positively correlated with fibrosis (unadjusted p < 0.005) by Gene Ontology enrichment analysis. This analysis revealed an abundance of proteins mainly related to catabolic processes, cytoskeleton organization and immune responses (Table 4). By contrast, the 76 proteins negatively associated with fibrosis (unadjusted p < 0.005) were mainly involved in normal metabolic processes and respiration (Table 4 and Supplementary Table S4).

3. Discussion

The current study investigated kidney graft biopsies from patients with creatinine clearance ≥ 30 mL/min and varying degrees of interstitial fibrosis by in-depth quantitative proteomic analysis using nano-LC-MSMS combined with 10-plex tandem mass tags (TMT), and the current study constitutes a first step towards the identification of non-invasive biomarkers of renal allograft fibrosis. Our proteomic analysis of renal tissue provides a detailed characterization of the protein composition of renal allografts; however, proteins of very low abundance may not be detected. Despite this, we identified four proteins—Coagulation Factor XIII A chain, Actin-related protein 2/3 subunit 2, Cytochrome C Oxidase Assembly Factor 6 homolog and Uridine Phosphorylase 1—that, even after multiple testing, showed a strong correlation with renal allograft fibrosis. Moreover, three of these proteins were also shown to be strongly predictive for fibrosis.
Coagulation factor XIII (FXIII) A chain was positively correlated with the degree of fibrosis and was also chosen by the LASSO in 21.7% of cases. FXIII is synthesized in cells of bone-marrow origin (i.e., thrombocytes and monocytes) [13]. Circulating FXIII is a tetramer consisting of two A chains and two B chains. The A chain can be activated by thrombin and calcium to become a transglutaminase that cross-links and stabilizes fibrin, thereby catalyzing the final steps in the coagulation cascade [14]. FXIII also cross-links proteins of the extracellular matrix [15], increases the proliferation of monocytes, and further increases the migration and decreases the apoptosis of both monocytes and fibroblasts [16]. While the recruitment and longevity of fibroblasts and the cross-linking of ECM proteins are beneficial in the context of wound healing, they constitute core processes in the progression of tissue fibrosis. Thus, it is feasible that the overexpression of FXIII might serve to promote fibrogenesis through increased fibroblast activity and inflammation, and reduced ECM degradation.
Actin-related protein 2/3 (ARP2/3) subunit 2 was similarly strongly associated with fibrosis. ARP2/3 is a protein complex involved in actin cytoskeleton organization and has been associated with a range of cell functions, primarily relating to cell motility [17]. Of possible interest in relation to fibrosis is the observation that disrupting actin cytoskeleton assembly via ARP2/3 in fibroblasts reduces fibroblast motility in an in vitro model of wound healing, presumably by preventing the polarization of the Golgi apparatus [18]. A high renal abundance of ARP2/3 might be an indicator of ongoing inflammation or increased fibroblast motility.
Cytochrome C oxidase is located in the mitochondrial membrane and plays a major role in cell respiration and ATP synthesis in eukaryotic cells [19]. Cytochrome C oxidase assembly factor, also located in the mitochondrion, contributes to the correct assembly and function of Cytochrome C oxidase [20]. Cytochrome C oxidase assembly factor was negatively correlated to fibrosis and chosen by the LASSO in 56.7% of cases, which indicates that it is a strong predictor of fibrosis. This observation may reflect a lack of healthy cells in the fibrotic kidneys. A proteomic study of CsA toxicity to human proximal tubular cells in vitro found a reduction of proteins from the inner mitochondrial membrane associated with CsA treatment [21]. One obvious difference to the current study, however, was the use of cells rather than renal biopsies. The proteomic analysis of biopsy specimens includes all parts of the nephron, as well as invading cells and fibrotic areas, which constitute a more heterogeneous sample than a single cell line.
Uridine phosphorylase catalyzes the reversible conversion of uridine to uracil, which is an important step in the pyrimidine salvage pathway [22]. We found a positive correlation between this protein and the extent of fibrosis. The previously mentioned in vitro study of human proximal tubule cells treated with CsA found an upregulation of proteins associated with purine salvage in relation to CsA treatment. The authors hypothesized that respiratory chain dysfunction led the cells to use this pathway rather than the more energy-consuming de novo synthesis [21]. Uridine phosphorylase is present in most human cells and can be induced by inflammatory cytokines in a number of tumor cell lines [23]. It has previously been shown to co-localize with intermediate filament vimentin in fibroblasts and colon cells; however, the significance of this co-localization is unknown [22].
Functional annotation of the proteins most significantly associated with fibrosis and eGFR revealed a marked difference between those positively and those negatively associated, i.e., proteins positively associated with fibrosis were negatively associated with eGFR and vice versa, indicating that kidney function and fibrosis are linked processes at the proteome microenvironment level (Table 4 and Supplementary Table S5). The observation that proteins involved in cytoskeleton organization and immune responses were abundant in fibrotic tissue has been reported in previous studies.
Nakorchevsky et al. performed a proteomic analysis of renal allograft biopsies divided into three groups: no IFTA, mild IFTA or moderate to severe IFTA, and subsequently compared protein abundances between the groups. The authors identified 492 proteins uniquely expressed across the groups, and a further 904 proteins differentially expressed between the groups. Functional annotation revealed that proteins related to acute phase responses, actin cytoskeleton signaling and complement activation were particularly related to severe IFTA, which points towards immunologic factors playing a role in the progression of IFTA [12]. Two of the proteins identified in the current study, Coagulation factor XIII and Actin-related protein 2/3, belong to the same functional annotations identified by Nakorchevsky to be of importance in the development of IFTA.
In a proteomic analysis of renal allografts from rats, Reuter et al. identified ten proteins that were differentially regulated in allogeneic transplants compared to in syngeneic transplants. The authors pointed towards imbalances of energy homeostasis and oxidative stress as causal explanations for the identified proteins [24]. None of the identified proteins were found to correlate with the degree of fibrosis in the current study. These differences may be due to differences between the species, but additionally, the renal changes associated with the allogeneic rat transplantation model resemble chronic rejection [25] and thus may reflect more immune-driven processes than the present biopsies from stable renal transplant patients. Moreover, a study by Späth et al. confirmed our findings on the involvement of proteins related to inflammatory and fibrotic responses in kidney function and fibrosis [26].
Traditional evaluation of fibrosis is performed by semi-quantitative Banff-scoring in which the pathologist grades the extent of fibrosis as CI1 (<25%), CI2 (26%–50%) or CI3 (>50%) [27]. Since this scale is non-linear and allows for large variations within each category, a morphometric method (point counting) was chosen to quantify the extent of fibrosis in the current study. A previous study of point counting in renal allograft biopsies found a coefficient of variability of 7% with repeated measurements and established point counting as a reproducible way of quantifying interstitial fibrosis [28]. In renal allograft biopsies obtained 6 months post-transplantation, interstitial fibrosis evaluated by point counting was inversely correlated with renal function and predicted allograft survival [29]; however, no significant correlation was found in the current study. This is likely due to the limited sample size. Furthermore, sampling error may have influenced the results in some cases. Previous studies have similarly noted that there is not a strictly monotonous relationship between chronic allograft damage and renal function [5].
The proteome analysis does not discriminate between the renal compartments nor distinguish resident cells from invading cells. The morphometric evaluation of interstitial fibrosis aimed to provide an objective and continuous measure of the extent of fibrosis. The decision to include perivascular fibrosis served to minimize the bias introduced by subjective evaluation, but may have caused the overestimation of the extent of fibrosis, and thus weakened the correlation to renal function.
We are aware that the plethora of information yielded from the proteomic analysis poses the risk of giving false positive results. Moreover, it is well-known that the quantitative proteomic approach based on relative quantification using 10-plex TMT isobaric tags applied in the present study suffers from dynamic range compression that affects the linearity of the signal, which may also affect biological conclusions [30]. This underlines the need to interpret results in a biological context. In the current study, we have, however, performed corrections for multiple testing to keep the risk of type I error low. In summary, we have provided an extensive characterization of the renal proteome in stable renal transplant patients, and the suggestion of four novel biomarkers of renal allograft fibrosis. Our results support the notion that cytoskeleton organization and immune responses are prevalent processes in renal allograft fibrosis. The prognostic value of the identified biomarkers of fibrosis remains to be validated in independent cohorts.

4. Methods

4.1. Participants

The study population included 31 patients from the SPIREN trial [31]. In brief, kidney transplant patients, at any time after the transplantation, were included from the outpatient clinic at Odense University Hospital. The eligibility criteria are listed in Table 5. At baseline, we performed a clinical examination, draw blood and urine samples, measured chrome-EDTA clearance, performed ambulatory blood pressure measurements, collected 24-hour urine samples and performed a kidney graft biopsy.
The collection of kidney graft biopsies and laboratory variables, fibrosis quantification and ambulatory blood pressure measurements were done as previously described [28]. GFR was determined by chrome-EDTA clearance following the standard procedure, in the Department of Nuclear Medicine, Odense University Hospital, in which a single injection of 51CrEDTA was given, followed by the taking of a venous blood sample after 240 min to determine residual radioactivity. An additional blood sample 24 h after injection was taken from male patients with p-creatinine ≥200 µmol/L and female patients with p-creatinine ≥150 µmol/L. Ultrasound-guided kidney graft biopsies were performed by trained senior physicians in the Department of Nephrology, kept in formalin and subsequently embedded in paraffin. All biopsies were scored according to the Banff classification by an experienced renal pathologist. The Banff classification is a semi-quantitative score, where each biopsy is given a score between 0 and 3 for acute (tubulitis (T), interstitial mononuclear infiltration (I), glomerulitis (G), vascular (V) and arteriolar hyalinosis (AH)) or chronic changes (chronic glomerulopathy (CG), interstitial fibrosis (CI), tubular atrophy (CT) and chronic vascular changes (CV)) [32]. The morphometric quantification of fibrosis was determined by point counting performed on Masson Trichrome-stained sections using the Cast 2.0 software. In brief, this method estimates the volume fraction of fibrosis by superimposing a grid with 12 intersection points on a computerized image of the biopsy. After manually delineating the renal cortex, the program randomly selects sections of the biopsy for quantification. The volume fraction of fibrosis is calculated by determining the number of intersection points that overlie fibrotic areas relative to the number of intersection points overlying normal renal tissue. In the current study, no exclusion of perivascular fibrosis was performed.

4.2. Proteomics

4.2.1. Sample Preparation

Sections of formalin-fixed, paraffin-embedded kidney biopsies were washed three times in chloroform. The proteins were extracted by dissolving the deparaffinized tissue sections in extraction buffer (1M dithiothretiol (DTT), 0.2 M tetraethylammonium bicarbonate (TEAB) and 10% sodium dodecyl sulphate (SDS)) followed by two rounds of ultrasonification (15 min of ultrasonification/15 min of cooling on ice), and incubated at 99 °C for 20 min then at 80 °C for 120 min. Protein alkylation was done by adding a 200 mM iodoacetamide (IAA) solution to a final DTT/IAA concentration ratio of 1:3. The acetone-precipitated proteins were re-dissolved in 5 µL 8 M urea with 1 µg LysC and incubated at 30 °C for 4 h, followed by a further dilution to 1 M urea, the addition of 2 µg trypsin, and an overnight incubation at 30 ˚C. The resulting tryptic peptides were isotopically labelled using the 10-plex tandem mass tag (TMT, ThermoFisher Scientific, Waltham, MA, USA). Peptide samples were randomly labelled with the 127N, 127C, 128N, 128C, 129N, 129C and 130N mass tags, whereas a pool of all samples were labelled with the mass tag 126 that served as an internal standard. Tagged peptides were mixed into 7 mixed peptide samples that were fractionated using hydrophilic interaction chromatography (HILIC), as described below.

4.2.2. HILIC Fractionation

Briefly, the labelled peptide mixtures were re-dissolved in 90% ACN/0.1% TFA, and 15 µL aliquots corresponding to approximately 25 µg peptides were injected onto an in-house-packed TSK gel Amide-80 HILIC 300 μm × 300 mm capillary high-performance liquid chromatography (HPLC) column and fractionated into 42 fractions by a Dionex UltiMate 3000 nanoHPLC using a linear 59 min gradient (85.5% B to 54% B; solvent A: 0.1%TFA; solvent B: 90% ACN/0.1% TFA) at a flow rate of 6 µL per minute. The fractions were automatically collected in micro-well plates at 1-minute intervals after UV detection at 210 nm. The fractions were dried by vacuum centrifugation, re-dissolved in 10 µL 0.1% TFA and analyzed by nano-LC–MS/MS, as described below.

4.2.3. NanoLC-MS/MS

A Q-Exactive mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) equipped with a nano-HPLC interface (Dionex UltiMate 3000 nano HPLC) was applied for nano-LC-MSMS analyses. The samples (5 µL) were loaded onto a custom-made, fused capillary pre-column (2 cm length, 360 µm OD, 75 µm ID) and separated on a custom-made fused capillary column (20 cm length, 360 µm OD, 100 µm ID—both columns were packed with ReporSil Pur C13 3 µm resin) using a linear gradient from 95% solution A (0.1% formic acid) to 30% B (100% acetonitrile in 0.1% formic acid) over 51 min, followed by 5 min at 90% B and 5 min at 98% A, at a flow rate of 300 nL per minute. The acquisition of mass spectra was done in positive ion mode using an automatic data-dependent switch between an Orbitrap survey MS scan in the mass range from 400 to 1200 m/z, and high-energy collisional dissociation fragmentation (HCD) and Orbitrap detection of the 15 most intense ions observed in the MS scan. The Orbitrap target values for the MS and MSMS scans were 1,000,000 and 50,000 ions, at resolutions of 70,000 and 35,000 at m/z 200. Fragmentation in the HCD cell was performed at the normalized collision energy of peptides and 31 eV for TMT-labelled peptides. The ion selection threshold was set to 17,000 counts. Selected sequenced ions were dynamically excluded for 60 s. The data are available via ProteomeXchange with the identifier PXD017867. Reviewer account details; Username: [email protected]; password: AOuDknvS.

4.3. Data Analysis

The Sequest search engine and Mascot search Engine (v. 2.2.3) integrated with the Proteome Discoverer (PD) version 2.1 software (Thermo Scientific) were used to search raw data files with the following criteria; Protein database: Uniprot/Swissprot (downloaded 7th November 2012, 452,768 entries) and restricted to humans. Trypsin, one missed cleavage allowed, carbamidomethylation at cysteines, and 10-plex TMT labelling at lysine and N-terminal amines were set as fixed parameters, while methionine oxidation and deamidation were set as dynamic. The precursor mass tolerance was set to 8 PPM and the MSMS tolerance was set to 0.05 Da. The peptide data were extracted using a Mascot significance threshold of 0.05 and a minimum peptide length of 6. The false discovery rate (FDR) was calculated using a decoy database search, and only high-confidence peptide identifications (false discovery rate <1%) were included. The data were normalized by using the “Total peptide amount” setting in the Reporter Ions Quantifier node in PD, i.e., the sum of abundance values for each channel over all peptides identified within a file was calculated, and the all channels were normalized against the channel with the highest total abundance. Search files were further processed using the Proteome Discoverer software. Gene Ontology analysis was performed using the GOrila tool for GO enrichment analysis (http://cbl-gorilla.cs.technion.ac.il/GOrilla/kh37217t/GOResults.html) [33].

4.4. Statistical Analysis

Analyses of baseline values were performed using Stata15 (StataCorp, College Station, TX, USA). The data are described by the median (interquartile range). Due to the non-Normal distribution of the fibrosis data, a non-parametric correlation test was used (Spearman’s rank correlation). The association of specific proteins with fibrosis was evaluated using the R software. For each protein, we applied a multiple linear regression model for the volume fraction of fibrosis, with sex, age and donor type (deceased/living related) as covariates. The resulting p-values were corrected for multiple testing by the Benjamini-Hochberg procedure. To find a subset of proteins predicting the degree of fibrosis, we used LASSO with stability selection [34]. This method identifies the proteins with the highest predictive value with regards to fibrosis. The LASSO was applied to a linear model where all proteins with no missing values (1973 proteins), sex, age, donor type, and time since transplantation were used as covariates. This was repeated 1000 times on a randomly selected subsample consisting of 16 samples. The probability of selecting each protein was computed.

4.5. Ethics

All subjects provided written consent to participate in the study. The project was performed in accordance with the International Conference of Harmonization—Good Clinical Practice (ICH—GCP) guidelines and has been approved by the Danish Research Ethics Committee and the Danish Data Protection Agency. The SPIREN trial was registered in ClinicalTrials.gov (05/17/2012): NCT01602861, EudraCT (05/31/2011): 2011-002243-98.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/7/2371/s1.

Author Contributions

L.A.M., H.C.T. and C.B. are the local investigators/sponsor. L.A.M. performed morphometric evaluation of biopsies, N.M. performed Banff scoring and supervised morphometric evaluation of biopsies. H.C.B. directed the proteomic analysis, performed gene ontology analysis, and wrote the proteomic chapter of the method section. A.M.S. and M.B. performed the statistical analyses. L.A.M. drafted the manuscript. All co-authors critically revised and approved the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for the proteomic study was obtained from Odense University Hospital Free Research Fund. The SPIREN trial was supported by the following funds: the Region of Southern Denmark Research- and PhD-funds, the Danish Kidney Association Research Fund, Helen and Einar Bjørnows Fund, Danish Society of Nephrology Travel Fund, King Christian X Fund, the Danish Medical Association Research Fund, the Medicine Fund of the Danish Regions and Odense University Hospital Board of Consultants Research Fund. The funding bodies had no influence on the design of the study, the collection, analysis and interpretation of data or in writing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wolfe, R.A.; Ashby, V.B.; Milford, E.L.; Ojo, A.; Ettenger, R.E.; Agodoa, L.Y.; Held, P.J.; Port, F.K. Comparison of Mortality in All Patients on Dialysis, Patients on Dialysis Awaiting Transplantation, and Recipients of a First Cadaveric Transplant. N. Engl. J. Med. 1999, 341, 1725–1730. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Khan, I.H.; Garratt, A.M.; Cody, D.J.; Catto, G.R.D.; Kumar, A.; Edward, N.; MacLeod, A.M. Patients’ perception of health on renal replacement therapy: Evaluation using a new instrument. Nephrol. Dial. Transplant. 1995, 10, 684–689. [Google Scholar] [PubMed]
  3. Hariharan, S.; Bresnahan, B.A.; Taranto, S.E.; McIntosh, M.J.; Stablein, N.; Johnson, C.P. Improved Graft Survival after Renal Transplantation in the United States, 1988 to 1996. N. Engl. J. Med. 2000, 342, 605–612. [Google Scholar] [CrossRef] [PubMed]
  4. Solez, K.; Colvin, R.B.; Racusen, L.C.; Sis, B.; Halloran, P.F.; Birk, P.E.; Campbell, P.M.; Cascalho, M.; Collins, A.B.; Demetris, A.J.; et al. Banff ’05 Meeting Report: Differential Diagnosis of Chronic Allograft Injury and Elimination of Chronic Allograft Nephropathy (‘CAN’). Am. J. Transplant. 2007, 7, 518–526. [Google Scholar] [CrossRef]
  5. Serón, D.; Moreso, F.; Ser, D. Protocol biopsies in renal transplantation: Prognostic value of structural monitoring. Kidney Int. 2007, 72, 690–697. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, Y. Cellular and molecular mechanisms of renal fibrosis. Nat. Rev. Nephrol. 2011, 7, 684–696. [Google Scholar] [CrossRef]
  7. Fletcher, J.T.; Nankivell, B.J.; Alexander, S.I. Chronic allograft nephropathy. Pediatr. Nephrol. 2009, 24, 1465–1471. [Google Scholar] [CrossRef] [Green Version]
  8. Schwarz, A.; Gwinner, W.; Hiss, M.; Radermacher, J.; Mengel, M.; Haller, H. Safety and Adequacy of Renal Transplant Protocol Biopsies. Am. J. Transplant. 2005, 5, 1992–1996. [Google Scholar] [CrossRef]
  9. Farris, A.B.; Colvin, R.B. Renal interstitial fibrosis: Mechanisms and evaluation. Curr. Opin. Nephrol. Hypertens. 2012, 21, 289–300. [Google Scholar] [CrossRef] [Green Version]
  10. Rhee, E.P. How Omics Data Can Be Used in Nephrology. Am. J. Kidney Dis. 2018, 72, 129–135. [Google Scholar] [CrossRef]
  11. Knepper, M.A. Common sense approaches to urinary biomarker study design. J. Am. Soc. Nephrol. 2009, 20, 1175–1178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Nakorchevsky, A.; Hewel, J.A.; Kurian, S.; Mondala, T.S.; Campbell, D.; Head, S.R.; Marsh, C.; Yates, J.R.; Salomon, D.R. Molecular mechanisms of chronic kidney transplant rejection via large-scale proteogenomic analysis of tissue biopsies. J. Am. Soc. Nephrol. 2010, 21, 362–373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Bagoly, Z.; Katona, É.; Muszbek, L. Factor XIII and inflammatory cells. Thromb. Res. 2012, 129, S77–S81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bereczky, Z.; Katona, Ë.; Muszbek, L. Fibrin Stabilization (Factor XIII), Fibrin Structure and Thrombosis. Pathophysiol. Haemost. Thromb. 2003, 33, 430–437. [Google Scholar] [CrossRef] [PubMed]
  15. Kohler, H.P.; Schroeder, V. Factor XIII: Structure and Function. Semin. Thromb. Hemost. 2016, 42, 422–428. [Google Scholar] [CrossRef]
  16. Dardik, R.; Krapp, T.; Rosenthal, E.; Loscalzo, J.; Inbal, A. Effect of FXIII on Monocyte and Fibroblast Function. Cell. Physiol. Biochem. 2007, 19, 113–120. [Google Scholar] [CrossRef]
  17. Welch, M.D.; DePace, A.H.; Verma, S.; Iwamatsu, A.; Mitchison, T.J. The human Arp2/3 complex is composed of evolutionarily conserved subunits and is localized to cellular regions of dynamic actin filament assembly. J. Cell Biol. 1997, 138, 375–384. [Google Scholar] [CrossRef] [Green Version]
  18. Magdalena, J.; Millard, T.H.; Etienne-Manneville, S.; Launay, S.; Warwick, H.K.; Machesky, L.M. Involvement of the Arp2/3 complex and Scar2 in Golgi polarity in scratch wound models. Mol. Biol. Cell. 2003, 14, 670–684. [Google Scholar] [CrossRef] [Green Version]
  19. Barrientos, A.; De Barros, M.H.; Valnot, I.; Rötig, A.; Rustin, P.; Tzagoloff, A. Cytochrome oxidase in health and disease. Gene 2002, 286, 53–63. [Google Scholar] [CrossRef]
  20. Pacheu-Grau, D.; Bareth, B.; Dudek, J.; Juris, L.; Vögtle, F.-N.; Wissel, M.; Leary, S.C.; Dennerlein, S.; Rehling, P.; Deckers, M. Cooperation between COA6 and SCO2 in COX2 Maturation during Cytochrome c Oxidase Assembly Links Two Mitochondrial Cardiomyopathies. Cell Metab. 2015, 21, 823–833. [Google Scholar] [CrossRef] [Green Version]
  21. Puigmulé, M.; Lopez-Hellin, J.; Suñè, G.; Tornavaca, O.; Camano, S.; Tejedor, A.; Meseguer, A. Differential proteomic analysis of cyclosporine A-induced toxicity in renal proximal tubule cells. Nephrol. Dial. Transplant. 2009, 24, 2672–2686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Russell, R.L.; Cao, D.; Zhang, D.; Handschumacher, R.E.; Pizzorno, G. Uridine Phosphorylase Association with Vimentin. Intracellular distribution and localization. J. Boil. Chem. 2001, 276, 13302–13307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pizzorno, G.; Cao, D.; Leffert, J.J.; Russell, R.L.; Zhang, D.; E Handschumacher, R. Homeostatic control of uridine and the role of uridine phosphorylase: A biological and clinical update. Biochim. Biophys. Acta (BBA)—Mol. Basis Dis. 2002, 1587, 133–144. [Google Scholar] [CrossRef] [Green Version]
  24. Reuter, S.; Reiermann, S.; Wörner, R.; Schröter, R.; Edemir, B.; Buck, F.; Henning, S.; Peter-Katalinić, J.; Vollenbröker, B.; Amann, K.; et al. IF/TA-related metabolic changes—Proteome analysis of rat renal allografts. Nephrol. Dial. Transplant. 2010, 25, 2492–2501. [Google Scholar] [CrossRef] [Green Version]
  25. Kunter, U.; Floege, J.; Von Jurgensonn, A.S.; Stojanovic, T.; Merkel, S.; Grone, H.J.; Ferran., C. Expression of A20 in the vessel wall of rat-kidney allografts correlates with protection from transplant arteriosclerosis. Transplant. 2003, 75, 3–9. [Google Scholar] [CrossRef]
  26. Späth, M.R.; Bartram, M.P.; Palacio-Escat, N.; Hoyer, K.J.R.; Debes, C.; Demir, F.; Schroeter, C.B.; Mandel, A.M.; Grundmann, F.; Ciarimboli, G.; et al. The proteome microenvironment determines the protective effect of preconditioning in cisplatin-induced acute kidney injury. Kidney Int. 2019, 95, 333–349. [Google Scholar] [CrossRef]
  27. Racusen, L.C.; Solez, K.; Colvin, R.B.; Bonsib, S.M.; Castro, M.C.; Cavallo, T.; Croker, B.P.; Demetris, A.J.; Drachenberg, C.B.; Fogo, A.B.; et al. The Banff 97 working classification of renal allograft pathology. Kidney Int. 1999, 55, 713–723. [Google Scholar] [CrossRef] [Green Version]
  28. Ellingsen, A.R.; Nyengaard, J.R.; Østerby, R.; Jørgensen, K.A.; Petersen, S.E.; Marcussen, N. Measurements of cortical interstitium in biopsies from human kidney grafts: How representative and how reproducible? Nephrol. Dial. Transplant. 2002, 17, 788–792. [Google Scholar] [CrossRef] [Green Version]
  29. Nicholson, M.L.; McCulloch, T.A.; Harper, S.J.; Wheatley, T.J.; Edwards, C.M.; Feehally, J.; Furness, P.N. Early measurement of interstitial fibrosis predicts long-term renal function and graft survival in renal transplantation. BJS 1996, 83, 1082–1085. [Google Scholar] [CrossRef]
  30. Williamson, J.; Edwards, A.V.G.; Verano-Braga, T.; Schwämmle, V.; Kjeldsen, F.; Jensen, O.N.; Larsen, M.R. High-performance hybrid Orbitrap mass spectrometers for quantitative proteome analysis: Observations and implications. Proteomics 2016, 16, 907–914. [Google Scholar] [CrossRef]
  31. Mortensen, L.A.; Thiesson, H.C.; Tougaard, B.; Egfjord, M.; Fischer, A.S.L.; Bistrup, C. The effect of spironolactone on calcineurin inhibitor induced nephrotoxicity: A multicenter randomized, double-blind, clinical trial (the SPIREN trial). BMC Nephrol. 2018, 19, 105. [Google Scholar] [CrossRef] [Green Version]
  32. Solez, K.; Axelsen, R.A.; Benediktsson, H.; Burdick, J.F.; Cohen, A.H.; Colvin, R.B.; Croker, B.P.; Droz, D.; Dunnill, M.S.; Halloran, P.F.; et al. International standardization of criteria for the histologic diagnosis of renal allograft rejection: The Banff working classification of kidney transplant pathology. Kidney Int. 1993, 44, 411–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Eden, E.; Navon, R.; Steinfeld, I.; Lipson, D.; Yakhini, Z. GOrilla: A tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinform. 2009, 10, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Meinshausen, N.; Buhlmann, P. Stability selection. J. R. Stat. Soc. B 2010, 72, 417–473. [Google Scholar] [CrossRef]
Figure 1. (a) The distribution of interstitial fibrosis by point counting, whiskers indicate upper and lower adjacant values, and dot represent an outlier above the upper adjacant value; (b) The association between renal function and the degree of fibrosis (Spearmans ρ = −0.23, p = 0.21).
Figure 1. (a) The distribution of interstitial fibrosis by point counting, whiskers indicate upper and lower adjacant values, and dot represent an outlier above the upper adjacant value; (b) The association between renal function and the degree of fibrosis (Spearmans ρ = −0.23, p = 0.21).
Ijms 21 02371 g001
Figure 2. The associations of specific proteins with the degree of fibrosis. The x-axis shows the relative protein abundances (AD).
Figure 2. The associations of specific proteins with the degree of fibrosis. The x-axis shows the relative protein abundances (AD).
Ijms 21 02371 g002
Table 1. The clinical characteristics of the study population (n = 31). The data are presented as medians (interquartile range).
Table 1. The clinical characteristics of the study population (n = 31). The data are presented as medians (interquartile range).
ItemValue
Male sex22 (71%)
Body Mass Index (kg/m2)26 (24–28)
Age (years)60 (50–66)
Smoking (never/previous/current)16/9/6 (52%/29%/19%)
Age of current transplant (years)3.0 (1.6–9.5)
Donor type
Living related12
Deceased19
Previous dialysis24 (77%)
Previous rejection3 (10%)
Borderline1 (3%)
Grade 1A2 (6%)
Cause of renal failure
Hypertension6 (19%)
Diabetes2 (6%)
Polycystic kidney disease9 (29%)
IgA nephritis2 (6%)
Glomerulonepritis3 (10%)
Congenital malformation3 (10%)
Reflux nefropathy1 (3%)
Unknown5 (16%)
Comorbidity
Hypertension29 (94%)
Diabetes2 (6%)
Apoplexia1 (3%)
Myocardial infarction3 (10%)
Heart failure0 (0%)
Current medication
ACE-inhibitors14 (45%)
Angiotensin II receptor blockers2 (6%)
Calcium channel antagonists18 (58%)
β-adrenergic antagonists17 (55%)
α-adrenergic antagonists5 (16%)
Loop diuretics4 (13%)
Thiazide diuretics2 (6%)
Cholesterol lowering drugs4 (13%)
Calcineurin inhibition:
 -Tacrolimus25 (81%)
 -Cyclosporine A6 (19%)
Antimetabolite:
 -Mycophenolate mofetil29 (94%)
 -Azathioprine2 (6%)
Prednisolone4 (13%)
Systolic blood pressure (mmHg)131 (118–141)
Diastolic blood pressure (mmHg)79 (75–83)
Chrome-EDTA clearance (mL/min)49 (39–60)
Fibrosis (%)33 (30–41)
Table 2. Proteins with FDR < 0.5 multiple linear regressions after correction for multiple testing.
Table 2. Proteins with FDR < 0.5 multiple linear regressions after correction for multiple testing.
ProteinAccesion Number AEstimate BSE Cp-ValueAdjusted p-Value DMissing values E
Uridine phosphorylase 1Q1683119.42.53.06 × 10-78.21 × 10-48
Cytochrome c oxidase assembly factor 6 homologQ5JTJ3.2−44.96.96,98 × 10-71.88 × 10-30
Coagulation factor XIII A chainP0048818.73.49.09 × 10-62.44 × 10-20
Actin-related protein 2/3 complex subunit 2O1514434.36.51.78 × 10-54.77 × 10-20
WD repeat-containing protein 82Q6UXN934.86.03.38 × 10-59.08 × 10-211
Actin-related protein 2/3 complex subunit 1BO1514321.64.55.17 × 10-51.39 × 10-10
Heterogeneous nuclear ribonucleoprotein UQ0083948.310.05.38 × 10-51.44 × 10-10
Mitochondrial aconitate hydrataseA2A274−47.210.18.20 × 10-52.20 × 10-10
Hematopoietic lineage cell-specific proteinP1431711.32.59.49 × 10-52.54 × 10-10
Lamina-associated polypeptide 2, isoforms beta/gammaP4216725.75.61.06 × 10-42.84 × 10-10
NADH dehydrogenase 1 alpha subcomplex subunit 2O43678−38.48.51.13 × 10-43.01 × 10-10
60S ribosomal protein L7aP6242456.212.41.17 × 10-43.14 × 10-10
Cathepsin GP083119.92.21.18 × 10-43.16 × 10-10
Perilipin-3O6066433.67.31.21 × 10-43.24 × 10-13
Dimethylglycine dehydrogenase, mitochondrialQ9UI17−26.65.91.22 × 10-43.26 × 10-10
60S ribosomal protein L5P4677736.38.11.28 × 10-43.43 × 10-10
Proteasome subunit beta type-8P2806226.86.01.35 × 10-43.61 × 10-10
HLA class I histocompatibility antigen, B-42 alpha chainP3048017.63.91.37 × 10-43.65 × 10-10
Xaa-Pro aminopeptidase 2O43895−27.56.21.41 × 10-43.78 × 10-10
Coronin-1AP311465.91.31.80 × 10-43.94 × 10-10
Tubulin beta chainP0743736.38.21.54 × 10-44.10 × 10-10
Splicing factor U2AF 65 kDa subunitP2636825.75.81.54 × 10-44.11 × 10-10
Short/branched chain specific acyl-CoA dehydrogenase, mitochondrialP45954−26.96.11.58 × 10-44.20 × 10-10
X-ray repair cross-complementing protein 6P1295640.29.11.63 × 10-44.33 × 10-10
Delta(3,5)-Delta(2,4)-dienoyl-CoA isomerase, mitochondrialQ13011−37.28.51.72 × 10-44.58 × 10-10
Sodium-dependent neutral amino acid transporter B(0)AT1Q695T7−21.44.91.79 × 10-44.77 × 10-10
A: Swissprot accession number; B: Estimated regression coefficient; C: Standard Error; D: Values are corrected according to the Benjamin-Hochberg procedure; E: The number of missing values for the specific protein out of the study population of 31 patients.
Table 3. Proteins chosen by the LASSO in more than 10% of all subsamples.
Table 3. Proteins chosen by the LASSO in more than 10% of all subsamples.
ProteinAccession NumberCorrelationBiological ProcessSelection Probability
Cytochrome c oxidase assembly factor 6 homologQ5JTJ3.2NegativeCell respiration0.567
Apoptosis inhibitor 5Q9BZZ5.2PositiveInhibition of (fibroblast) apoptosis0.234
Coagulation factor XIII A chainP00488PositiveCoagulation/cross-linking ECM proteins/inhibition of fibroblast apoptosis0.217
Microtubule-associated protein 1BP46821PositiveNeuronal cell structure maintainance0.205
Splicing factor U2AF 65 kDa subunitP26368PositivemRNA processing0.199
SymplekinQ92797PositiveCell adhesion/mRNA processing0.156
Threonine synthase-like 1Q8IYQ7NegativeUnknown0.131
AP-2 complex subunit sigmaP53680PositiveProtein transport0.118
Actin-related protein 2/3 complex subunit 2O15144PositiveActin cytoskeleton organization0.107
Table 4. Significantly enriched Gene Ontology categories. Gene ontology (GO) enrichment analysis of proteins that are significantly (p < 0.005) positively and negatively associated with fibrosis as determined by multiple regression analysis (Supplementary Table S4). GO categories were identified using two unranked list options of the GOrilla Gene Ontology (GO) program (http://cbl-gorilla.cs.technion.ac.il/). N is the total number of recognized proteins, B is the total number of recognized proteins associated with a specific GO term, n is the number of recognized proteins in the target set, and b is the number of proteins in the intersection. Enrichment = (b/n)/(B/N). The p-Value is the enrichment p-value computed according to the mHG or HG model. This p-value is not corrected for the multiple testing of 2893 GO terms, and the FDR q-value is the correction of the above p-value for multiple testing using the Benjamini and Hochberg method.
Table 4. Significantly enriched Gene Ontology categories. Gene ontology (GO) enrichment analysis of proteins that are significantly (p < 0.005) positively and negatively associated with fibrosis as determined by multiple regression analysis (Supplementary Table S4). GO categories were identified using two unranked list options of the GOrilla Gene Ontology (GO) program (http://cbl-gorilla.cs.technion.ac.il/). N is the total number of recognized proteins, B is the total number of recognized proteins associated with a specific GO term, n is the number of recognized proteins in the target set, and b is the number of proteins in the intersection. Enrichment = (b/n)/(B/N). The p-Value is the enrichment p-value computed according to the mHG or HG model. This p-value is not corrected for the multiple testing of 2893 GO terms, and the FDR q-value is the correction of the above p-value for multiple testing using the Benjamini and Hochberg method.
GO TermDescriptionp-ValueFDR q-ValueEnrichment (N, B, n, b)
Categories positively associated with fibrosis
GO:0010468regulation of gene expression2.09 × 10-71.97 × 10-31.89 (2240,546,113,52)
GO:0051704multi-organism process2.37 × 10-71.11 × 10-32.22 (2240, 348, 113, 39)
GO:0050789regulation of biological process5.94 × 10-71.86 × 10-31.34 (2240, 1371, 113, 93)
GO:0010629negative regulation of gene expression6.29 × 10-71.48 × 10-32.47 (2240, 249, 113, 31)
GO:0044419interspecies interaction between organisms8.74 × 10-71.64 × 10-32.38 (2240, 266, 113, 32)
GO:0002376immune system process1.27 × 10-61.99 × 10-31.96 (2240, 446, 113, 44)
GO:0002682regulation of immune system process1.41 × 10-61.90 × 10-32.38 (2240, 258, 113, 31)
GO:0045321leukocyte activation1.46 × 10-61.72 × 10-32.43 (2240, 245, 113, 30)
GO:0006952defense response1.83 × 10-61.91 × 10-32.87 (2240, 159, 113, 23)
GO:0002252immune effector process3.34 × 10-63.14 × 10-32.21 (2240, 296, 113, 33)
GO:0048518positive regulation of biological process3.65 × 10-63.12 × 10-31.55 (2240, 844, 113, 66)
GO:0006955immune response5.25 × 10-64.11 × 10-32.98 (2240, 133, 113, 20)
GO:0050794regulation of cellular process5.83 × 10-64.21 × 10-31.36 (2240, 1239, 113, 85)
GO:0007166cell surface receptor signaling pathway8.47 × 10-65.68 × 10-32.16 (2240, 294, 113, 32)
GO:0001775cell activation9.28 × 10-65.81 × 10-32.23 (2240, 267, 113, 30)
GO:0002684positive regulation of immune system process9.39 × 10-65.52 × 10-32.54 (2240, 187, 113, 24)
Categories negatively associated with fibrosis
GO:0055114oxidation-reduction process1.94 × 10-151.79 × 10-113.75 (2198, 318, 70, 38)
GO:0019752carboxylic acid metabolic process1.16 × 10-115.37 × 10-83.16 (2198, 348, 70, 35)
GO:0043436oxoacid metabolic process4.53 × 10-111.39 × 10-73.02 (2198, 364, 70, 38)
GO:0044281small molecule metabolic process6.03 × 10-111.39 × 10-72.40 (2198, 576, 70, 44)
GO:0006082organic acid metabolic process8.68 × 10-111.60 × 10-72.95 (2198, 372, 70, 35)
Table 5. The eligibility criteria.
Table 5. The eligibility criteria.
Inclusion CriteriaExclusion Criteria
Age > 18 yearsFormer intolerance of spironolactone
Tacrolimus/cyclosporine treatment Potassium binder or digoxin treatment
Proteinuria < 3 g/dayPregnancy or planned pregnancy
Creatinine clearance ≥ 30 mL/minClinically relevant organic, systemic or psychological disorder
Plasma potassium < 5.5 mmol/LExpectation of non-compliance
Negative pregnancy test at inclusion for women of childbearing potential and adequate contraception throughout the trial

Share and Cite

MDPI and ACS Style

Mortensen, L.A.; Svane, A.M.; Burton, M.; Bistrup, C.; Thiesson, H.C.; Marcussen, N.; Beck, H.C. Proteomic Analysis of Renal Biomarkers of Kidney Allograft Fibrosis—A Study in Renal Transplant Patients. Int. J. Mol. Sci. 2020, 21, 2371. https://doi.org/10.3390/ijms21072371

AMA Style

Mortensen LA, Svane AM, Burton M, Bistrup C, Thiesson HC, Marcussen N, Beck HC. Proteomic Analysis of Renal Biomarkers of Kidney Allograft Fibrosis—A Study in Renal Transplant Patients. International Journal of Molecular Sciences. 2020; 21(7):2371. https://doi.org/10.3390/ijms21072371

Chicago/Turabian Style

Mortensen, Line Aas, Anne Marie Svane, Mark Burton, Claus Bistrup, Helle Charlotte Thiesson, Niels Marcussen, and Hans Christian Beck. 2020. "Proteomic Analysis of Renal Biomarkers of Kidney Allograft Fibrosis—A Study in Renal Transplant Patients" International Journal of Molecular Sciences 21, no. 7: 2371. https://doi.org/10.3390/ijms21072371

APA Style

Mortensen, L. A., Svane, A. M., Burton, M., Bistrup, C., Thiesson, H. C., Marcussen, N., & Beck, H. C. (2020). Proteomic Analysis of Renal Biomarkers of Kidney Allograft Fibrosis—A Study in Renal Transplant Patients. International Journal of Molecular Sciences, 21(7), 2371. https://doi.org/10.3390/ijms21072371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop