Research Paper Volume 11, Issue 13 pp 4382—4406

Complement component C5a induces aberrant epigenetic modifications in renal tubular epithelial cells accelerating senescence by Wnt4/βcatenin signaling after ischemia/reperfusion injury

Giuseppe Castellano1, *, , Rossana Franzin1, *, , Fabio Sallustio1,2, , Alessandra Stasi1, , Barbara Banelli3, , Massimo Romani3, , Giuseppe De Palma1, , Giuseppe Lucarelli4, , Chiara Divella1, , Michele Battaglia4, , Antonio Crovace5, , Francesco Staffieri5, , Giuseppe Grandaliano6, , Giovanni Stallone6, , Pasquale Ditonno4, , Paolo Cravedi7, , Vincenzo Cantaluppi8, , Loreto Gesualdo1, ,

  • 1 Nephrology, Dialysis and Transplantation Unit, University of Bari “Aldo Moro”, Bari, Italy
  • 2 Department of Basic Medical Sciences, Neuroscience and Sense Organs, University of Bari “Aldo Moro”, Bari, Italy
  • 3 IRCCS Ospedale Policlinico San Martino, Genova, Italy
  • 4 Urology, Andrology and Renal Transplantation Unit, University of Bari, Bari, Italy
  • 5 Veterinary Surgery Unit, Department of Emergency and Organ Transplantation, University of Bari, Bari, Italy
  • 6 Nephrology, Dialysis and Transplantation Unit, Department of Medical and Surgical Sciences, University of Foggia, Foggia, Italy
  • 7 Department of Medicine, Translational Transplant Research Center, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
  • 8 Nephrology and Kidney Transplantation Unit, Department of Translational Medicine, University of Eastern Piedmont “A. Avogadro” (UPO), Novara, Italy
* Equal contribution

Received: March 8, 2019       Accepted: June 24, 2019       Published: July 8, 2019      

https://doi.org/10.18632/aging.102059
How to Cite

Copyright: Castellano et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Epigenetic mechanisms, such as DNA methylation, affect tubular maladaptive response after Acute Kidney Injury (AKI) and accelerate renal aging. Upon ischemia/reperfusion (I/R) injury, Complement activation leads to C5a release that mediates damage; however, little is known about the effect of C5a-C5a Receptor (C5aR) interaction in Renal Tubular Epithelial Cells (RTEC).

Through a whole-genome DNA methylation analysis in cultured RTEC, we found that C5a induced aberrant methylation, particularly in regions involved in cell cycle control, DNA damage and Wnt signaling. The most represented genes were BCL9, CYP1B1 and CDK6. C5a stimulation of RTEC led to up-regulation of SA-β Gal and cell cycle arrest markers such as p53 and p21. C5a increased also IL-6, MCP-1 and CTGF gene expression, consistent with SASP development. In accordance, in a swine model of renal I/R injury, we found the increased expression of Wnt4 and βcatenin correlating with SA-β Gal, p21, p16 and IL-6 positivity. Administration of Complement Inhibitor (C1-Inh), antagonized SASP by reducing SA-β Gal, p21, p16, IL-6 and abrogating Wnt4/βcatenin activation.

Thus, C5a affects the DNA methylation of genes involved in tubular senescence. Targeting epigenetic programs and Complement may offer novels strategies to protect tubular cells from accelerated aging and to counteract progression to Chronic Kidney Disease

Introduction

Acute Kidney Injury (AKI) is a frequent condition in hospitalized patients and is mainly caused by Ischemia/Reperfusion (I/R) injury, sepsis or nephrotoxic drugs [1]. Currently, AKI is associated with unacceptably high mortality [1] or development of Chronic Kidney Disease (CKD); therefore, there is an urgent need to develop new therapeutic strategies capable to interfere with detrimental effects of AKI and inhibit AKI-to CKD transition [2, 3].

Recent studies have suggested a new hypothesis on pathogenesis of I/R induced-AKI based on epigenetic modification such as DNA methylation and histone acetylation [46]. The genome interacts with several environmental factors such as nutrients, pathogens, drugs and toxins that can modify the chromatin condensation to make specific genes accessible or not accessible to transcription factors, thereby extensively regulating gene expression. These modifications become stable and heritable upon mitosis and have been associated to the risk of AKI, allograft rejection after transplantation [7] or transition to CKD [810].

DNA methylation is catalyzed by enzyme called DNMT (DNA methyl transferase) and occurs in region highly rich in dinucleotide Cytosine-Guanine called CpG islands, localized above all in the promoters but also in intra and inter-genic regions. High methylation of promoter is associated to a closed chromatin, then to transcriptional silencing [11]. Recently, DNA methylation modifications have been shown also to accelerate renal aging [12]; interestingly, hypermethylation of Klotho promoter, the principal anti-aging and renoprotective factor, led to a reduced Klotho gene expression with a significant association with CKD severity [10].

To date, epigenomic modifications offered drug targets already clinically available in the treatment of different diseases such as cancer; moreover, they showed efficacy also in experimental model of renal fibrosis [13, 14]. Recently, a whole kidney methylome was characterized in a mouse model of I/R describing the central role of DNA regulatory enzyme (DNMTs and Tet methylcytosine dioxygenase) in promoting renal fibrosis, inflammation and apoptosis [15]. Nevertheless, untargeted DNA demethylation therapy (with DNMTs-inhibitors) has not been successful in the attenuation of I/R induced fibrosis, highlighting the need of highly specific target genes [16].

In the pathophysiology of I/R injury, Complement system plays a pivotal role and the C1 esterase inhibitor C1-INH (C1-Inhibitor) is a promising strategy to prevent I/R damage in kidney transplantation and other forms of AKI [17, 18]. During I/R injury, Complement activation after C3 convertase generation led to the release of pro-inflammatory anaphylatoxin C5a that binds to C5a receptors expressed on monocytes, neutrophils and renal resident cells [1921]. In the present study, we aimed to investigate the effect of C5a on human renal tubular epithelial cells (RTEC) by a genome-wide DNA methylation analysis.

Results

C5a induces DNA methylation changes in RTECs

To study whether Complement could modulate DNA methylation in renal cells, we stimulated human RTECs with C5a anaphylatoxin for 24 hours and performed a whole-genome DNA methylation study. Methylation analysis was performed for single CpGs and also for sets of predefined genomic regions such as CpG islands, promoters and tiling regions (5-kb spanning regions). In untreated cells (indicated as basal), the highest number of methylated tiling regions were found in chromosome 1, 12, 10, 7 and 3. They cover the 50% of the total tiling regions, whereas chromosomes 4, 18 and 21 contained the fewest methylated regions (Figure 1A). C5a stimulation decreased the DNA methylation in the overall genome when we consider the single CpG sites (Figure 1B). In particular, considering the tiling regions, most DNA regions were hypomethylated and located in specific chromosomes (Figure 1C, Supplementary Table 1).

C5a-associated changes in DNA methylation as indicated by whole-genome bisulfite assay in RTEC. (A) Chart showing the number and the frequency of methylated regions (tiling regions) identified in RTEC at basal level. The left vertical axis represents the number of methylated regions per each chromosome. The right vertical axis indicates the cumulative percentage of the total number of occurrences. The red concave curve is the cumulative function indicating that the 50% of the total methylated regions in RTEC are covered by chromosomes 1, 12, 10, 7 and 3. (B) DNA methylation levels at the single CpG sites for RTEC at basal level (blue line) and for RTEC stimulated by C5a (orange line). C5a decreased the DNA methylation in the overall genome. (C) Graph showing the mean DNA methylation levels of tiling regions, shared for chromosomes, in RTEC (blue line) and in C5a-stimulated RTEC (red line). The central axis shows the mean β methylation value. (D) Scatterplot of the CpG site methylation comparison, colored according to the combined ranks of a given site. Values are represented as mean differences (mean.diff) between stimulated and unstimulated RTEC for each CpG site. (Combined rank: difference in mean methylation levels of stimulated and non-stimulated RTEC, the quotient in mean methylation and the t test are ranked for all regions. This value aggregates them using the maximum, i.e. worst rank of a site among the three measures.)

Figure 1. C5a-associated changes in DNA methylation as indicated by whole-genome bisulfite assay in RTEC. (A) Chart showing the number and the frequency of methylated regions (tiling regions) identified in RTEC at basal level. The left vertical axis represents the number of methylated regions per each chromosome. The right vertical axis indicates the cumulative percentage of the total number of occurrences. The red concave curve is the cumulative function indicating that the 50% of the total methylated regions in RTEC are covered by chromosomes 1, 12, 10, 7 and 3. (B) DNA methylation levels at the single CpG sites for RTEC at basal level (blue line) and for RTEC stimulated by C5a (orange line). C5a decreased the DNA methylation in the overall genome. (C) Graph showing the mean DNA methylation levels of tiling regions, shared for chromosomes, in RTEC (blue line) and in C5a-stimulated RTEC (red line). The central axis shows the mean β methylation value. (D) Scatterplot of the CpG site methylation comparison, colored according to the combined ranks of a given site. Values are represented as mean differences (mean.diff) between stimulated and unstimulated RTEC for each CpG site. (Combined rank: difference in mean methylation levels of stimulated and non-stimulated RTEC, the quotient in mean methylation and the t test are ranked for all regions. This value aggregates them using the maximum, i.e. worst rank of a site among the three measures.)

We then performed differential methylation analysis according to the combined Rank parameter and found 74 hypermethylated and 320 hypomethylated CpG sites in RTEC stimulated with C5a (Figure 1D). Next, after the differential methylation analysis of the tiling regions we selected the best 150 combined ranks having a threshold of Δβ>0.05. We found 142 regions, 88 hypomethylated and 54 hypermethylated (Supplementary Table 2). In addition, when we analyzed differentially methylated CpG island, we found only 3 island hypermethylated and 8 island hypomethylated by C5a (Supplementary Figure 4F).

In order to gain insight into the biological mechanisms involved in the response to C5a, we performed a pathway and a gene ontology (GO) analysis of genes spanned in the identified DNA regions (Supplementary Figure 1). The pathway analysis of regions hypo- and hypermethylated by C5a showed the involvement of biological processes as “Aryl Hydrocarbon Receptor Signaling”, “cell cycle: G2/M DNA damage checkpoint regulation” and “Wnt/βcatenin pathway” (p=0.019, 0.027 and p=0.033, respectively). In particular, we found three hypomethylated genes as central node of the networks: BCL9, CYP1B1 and CDK6 (Supplementary Figure 2). These genes are involved in the cell cycle progression and apoptosis. In addition, also the GO analysis showed an enrichment of DNA regions containing genes involved in the cell cycle checkpoint (Supplementary Figure 3).

We then validated the methylation status of the three genes by pyrosequencing (Figure 2A, 2B). We confirmed that BCL9, CYP1B1 and CDK6 genes were hypomethylated by C5a (p<0.05) with mean differences in methylation levels of about 8-10% (Figure 2A, 2B). In particular, the CYP1B1 region contained multiple CpG sites most of which were concordantly hypomethylated.

Analysis of methylation levels for BCL9, CYP1B1 and CDK6 in different lots of C5a-stimulated RTEC compared to the respective basal conditions. Pyrosequencing assays were performed on the same regions that we found to be methylated in the whole-genome assay. The pyrosequencing assays confirmed that the DNA of these regions were differentially methylated in C5a-stimulated RTEC (B panel) compared to basal condition (A panel). Representative pyrograms show methylation levels of the indicated CpG sites in the gene promoter.

Figure 2. Analysis of methylation levels for BCL9, CYP1B1 and CDK6 in different lots of C5a-stimulated RTEC compared to the respective basal conditions. Pyrosequencing assays were performed on the same regions that we found to be methylated in the whole-genome assay. The pyrosequencing assays confirmed that the DNA of these regions were differentially methylated in C5a-stimulated RTEC (B panel) compared to basal condition (A panel). Representative pyrograms show methylation levels of the indicated CpG sites in the gene promoter.

Overall, these data demonstrated the functional effect of the C5a anaphylatoxin to re-program in-depth the tubular cells epigenetic state, promoting a global hypomethylation in different chromosomes. In addition, the DNA methylation changes were consistent with cell cycle pathway.

C5a induced DNA methylation regulating BCL9, CYP1B1 and CDK6 gene expression

To confirm the critical role of DNA methylation in regulating the BCL9, CYP1B1 and CDK6, we assessed the gene expression from C5a-stimulated RTEC compared with untreated cells (Figure 3). The genes were differentially expressed in C5a-stimulated RTEC; their expression was inversely associated with the relative DNA methylation (Figure 3D3F): BCL9, CYP1B1 and CDK6 genes were hypomethylated with a significant up-regulation of their gene expression. (Figure 3A3C: mean fold changes of BLC9 2.19±0.5, CYP1B1 2.38±0.59 and CDK6 2.49±0.62, *p<0.05, **p<0.01).

BCL9, CYP1B1 and CDK6 gene expression is regulated by the DNA methylation. (A–C) Gene expression of BCL9, CYP1B1 and CDK6 evaluated by qRT-PCR in the C5a stimulated-RTEC and compared to normal RTEC cultured for 24h. The gene expression validated the genes analyzed for the methylation using the pyrosequencing assay. Expression levels were significantly different in C5a-stimulated RTEC compared with normal cells. Gene expression levels were normalized to the housekeeping gene GAPDH. (D–F) BCL9, CYP1B1 and CDK6 methylation levels of C5a-stimulated RTEC compared to basal condition. Results are means±SD, n=3. (G–I) Gene expression of BCL9, CYP1B1 and CDK6 of C5a stimulated-RTEC not treated or treated with 1 μM 5-aza-2’-deoxycytidine. The DNA demethylation agent reduced the methylation levels in the promoters of the three genes. The DNA methylation status of these three DNA regions was determined by qMSP real-time analysis. The degree of fully methylated molecules at a specific locus was expressed as a PMR index. The percentage PMR was calculated as described in the Materials and methods. Qiagen methylation control DNA was used as full methylated reference. (J–L) Gene expression of BCL9, CYP1B1 and CDK6 of C5a stimulated-RTEC not treated or treated with 1 μM aza. Results are means±SD, n = 3.*p

Figure 3. BCL9, CYP1B1 and CDK6 gene expression is regulated by the DNA methylation. (AC) Gene expression of BCL9, CYP1B1 and CDK6 evaluated by qRT-PCR in the C5a stimulated-RTEC and compared to normal RTEC cultured for 24h. The gene expression validated the genes analyzed for the methylation using the pyrosequencing assay. Expression levels were significantly different in C5a-stimulated RTEC compared with normal cells. Gene expression levels were normalized to the housekeeping gene GAPDH. (DF) BCL9, CYP1B1 and CDK6 methylation levels of C5a-stimulated RTEC compared to basal condition. Results are means±SD, n=3. (GI) Gene expression of BCL9, CYP1B1 and CDK6 of C5a stimulated-RTEC not treated or treated with 1 μM 5-aza-2’-deoxycytidine. The DNA demethylation agent reduced the methylation levels in the promoters of the three genes. The DNA methylation status of these three DNA regions was determined by qMSP real-time analysis. The degree of fully methylated molecules at a specific locus was expressed as a PMR index. The percentage PMR was calculated as described in the Materials and methods. Qiagen methylation control DNA was used as full methylated reference. (JL) Gene expression of BCL9, CYP1B1 and CDK6 of C5a stimulated-RTEC not treated or treated with 1 μM aza. Results are means±SD, n = 3.*p<0.05.

Next, to further validate the regulatory involvement of DNA methylation in BCL9, CYP1B1 and CDK6 gene expression, we treated RTEC for 24h with the DNA methyltransferase inhibitor 5-aza-2’-deoxycytidine (5-aza, 1μM), a DNA hypomethylating agent. Following 5-aza treatment, we assessed BCL9, CYP1B1 and CDK6 hypomethylation (Figure 3G, 3H, 3I). Effectively, after 5-aza induced-DNA demethylation, we found that the mRNA expression of BCL9, CYP1B1 and CDK6 by qPCR significantly increased (Figure 3J, 3K, 3L). These results indicated that gene expression of BCL9, CYP1B1 and CDK6 is directly regulated by DNA methylation and that C5a is an epigenetic mediator that caused hypomethylation of genes that are involved, at least in part, in the Wnt/βcatenin signaling.

Characterization of C5a-induced cellular senescence in RTEC

Because the three validated genes were clustered in Wnt/βcatenin pathway [22, 23] that is involved in aging processes [24], we tested the hypothesis that C5a might directly promote RTEC senescence in vitro. We then stimulated RTEC by C5a for 3h and 24h (Figure 4B and 4F, 4J) followed by culture in fresh medium for additional 48h (4B) or 24h (4J); interestingly, C5a stimulated-RTEC acquired a senescent phenotype as determined by significant increase in SA-βGAL positivity (Figure 4B, 4F, 4J). C5a exposed RTEC become larger (Figure 4F, arrow) and polynucleated (Figure 4B, arrow). The inhibition of receptor C5aR1 appeared to significantly protect tubular cells from C5a-induced senescence (Figure 4C, 4G and 4K), indicating that C5a effect was mediated by this receptor. RTEC exposed to H2O2 (Figure 4D, 4H and 4L) were used as positive control of cellular senescence. Quantification of senescent cells in Figure 4M showed that the induction of C5a mediated RTEC senescence was statistically significant in all the condition tested and also in presence of different time of activation (3h and 24h; *p<0.05, images at 3h+24 not showed).

SA-β Gal staining in RTEC treated with C5a with or without C5aR1 blocking. SA-β-Gal activity in early passage RTEC exposed to 100 nM C5a for 3h (B, F) or 24h (J). For C5aR1 inhibition, mouse monoclonal anti-C5aR1 was pre-incubated for 1h before the C5a exposure and then maintained in fresh medium for 24h or 48h (C, G, K). More SA-β-gal+ cells were observed after C5a exposure, senescent RTEC appeared enlarged and morphologically distinct from the normal cells at the same passage with formation of larger and polynucleated cells (B and F, arrows). Untreated cells are also named as Basal (A, E, I). H2O2 exposed cells were used as positive control of senescence (D, H, L). Representative images were acquired by phase contrast microscopy. (M) Quantification of SA-β-Gal+ cells cultures. The ratio of cells positive for SA-β-gal activity was calculated by examining five not overlapping fields per condition (6-well plate). The results are presented as the mean ± SD of three independent experiments (*p

Figure 4. SA-β Gal staining in RTEC treated with C5a with or without C5aR1 blocking. SA-β-Gal activity in early passage RTEC exposed to 100 nM C5a for 3h (B, F) or 24h (J). For C5aR1 inhibition, mouse monoclonal anti-C5aR1 was pre-incubated for 1h before the C5a exposure and then maintained in fresh medium for 24h or 48h (C, G, K). More SA-β-gal+ cells were observed after C5a exposure, senescent RTEC appeared enlarged and morphologically distinct from the normal cells at the same passage with formation of larger and polynucleated cells (B and F, arrows). Untreated cells are also named as Basal (A, E, I). H2O2 exposed cells were used as positive control of senescence (D, H, L). Representative images were acquired by phase contrast microscopy. (M) Quantification of SA-β-Gal+ cells cultures. The ratio of cells positive for SA-β-gal activity was calculated by examining five not overlapping fields per condition (6-well plate). The results are presented as the mean ± SD of three independent experiments (*p< 0.05), Magnification 40X.

Since cellular senescence is characterized by cell cycle arrest, we performed the MTT assay to assess the number of viable RTEC (Figure 5A). C5a led to a reduction in cellular proliferation compared to untreated cells (MTT C5a 24h: 0.134±0.03 vs basal 24h: 0.219±0.02); the inhibition of proliferation was stable even in presence of short C5a stimulation with 48h of culture in fresh medium (C5a 3h+48h: 0.013±0.04 vs basal). RTEC exposed to H2O2 were used as internal control showing a 54% inhibition of proliferation. In addition, cells counting performed by Trypan blue from the supernatants of cell culture did not indicate significant variations in the viable/dead cells number ratio (data not shown).

Altered cell proliferation and p53 expression after C5a exposure in RTEC. (A) MTT assay of RTEC in logarithmic growth phase treated with the same dosage of C5a showed an inhibition of cell proliferation. (B) TP53 gene expression level in C5a stimulated RTEC. Data are expressed as the mean ± SD (*p

Figure 5. Altered cell proliferation and p53 expression after C5a exposure in RTEC. (A) MTT assay of RTEC in logarithmic growth phase treated with the same dosage of C5a showed an inhibition of cell proliferation. (B) TP53 gene expression level in C5a stimulated RTEC. Data are expressed as the mean ± SD (*p< 0.05 versus basal condition).

Next, we tested the activation of p53 /TP53 gene pathway that is crucial in regulating premature senescence and cell- cycle arrest [25]; therefore, we performed a cell cycle assay and assessed the expression of TP53. Three hours after C5a exposure, we observed a significant increase in p53 transcript that further increased after 24h and 48h of culture (mean fold change: C5a 3h+48h: 5.6±1.64, p=0.008; C5a 3h+24h: 3.17±1.32, p=0.004, versus basal 1±0.0). The p53 expression increase was higher even after a short C5a stimulation (3h) followed by 48h of normal culture in fresh medium (Figure 5B); effectively most of cells were halted in G0/G1 phase of the cell cycle (83% versus 58% of non-stimulated cells), indicating that C5a affected the arrest of the cell cycle (Supplementary Figure 4A).

C5a induced the development of Secretory Associated Senescence Phenotype (SASP) in RTEC

After p53 activation, cells transiently express p21Cip1/Waf1 [26] a Cyclin-Dependent Kinase Inhibitor (CKI) that can lead to a chronic state of senescence, triggering G1 cell cycle arrest or leading to apoptosis. Therefore, we measured the p21 protein level in C5a stimulated RTEC and we found a significant protein increase after both short time (3h) and longer time (24h) of C5a exposure. (Figure 6A, 6B, p21 C5a 3h+24: 1.45±0.17; p<0.01, C5a 24h+24: 1.28±0.37, p<0.05). In particular, this effect was more evident after 24h of C5a exposure followed by 72h of normal culture. Since p21 could transiently increase after an injury, we assessed also the p16INK4a protein level. Stimulation with C5a significantly induced a constant augment in protein expression of p16INK4a compared to untreated condition (basal) (Figure 6A and 6B).

Cell cycle negative regulators expression levels and SASP pro-inflammatory cytokines transcripts are up regulated by C5a. (A–B) Representative p21 Waf1/Cip1 and p16INK4a western blot of C5a-stimulated RTEC and quantification. Protein expression was normalized to βactin. (C–I) Time course of IL-6, MCP1, CTGF, SERPIN1 (PAI-1), BCL9, CYP1B1 and CDK6 gene expression in RTEC stimulated for 3h or 24h by C5a followed or not by 24h or 48h of normal culture was assessed by qPCR. Data were normalized to GAPDH. *p2O2 (100-300 μm) was used as positive control of senescence.

Figure 6. Cell cycle negative regulators expression levels and SASP pro-inflammatory cytokines transcripts are up regulated by C5a. (AB) Representative p21 Waf1/Cip1 and p16INK4a western blot of C5a-stimulated RTEC and quantification. Protein expression was normalized to βactin. (CI) Time course of IL-6, MCP1, CTGF, SERPIN1 (PAI-1), BCL9, CYP1B1 and CDK6 gene expression in RTEC stimulated for 3h or 24h by C5a followed or not by 24h or 48h of normal culture was assessed by qPCR. Data were normalized to GAPDH. *p<0.05, **p<0,01, n=3. H2O2 (100-300 μm) was used as positive control of senescence.

Subsequently, we investigated the gene expression of IL-6, MCP-1, CTGF and PAI-1, that are chemokines typical associated with the Senescence-Associated Secretory Phenotype (SASP) [27]. As observed for p21 protein, we found a significant increase in IL-6, SERPINE1 (PAI-1), CCL2 (MCP-1) and CCN2 (CTGF) after three hours of C5a exposure followed by 24h of normal culture (Figure 6C6F; mean fold change IL-6: 8.05±0.24; PAI-1: 6.58±1.72 versus baseline, p<0.01). Interestingly, IL-6, PAI-1 and CTGF gene expression levels remained high even after 48h of culture from C5a stimulation. Significant expression increases for all these genes, except MCP1, was found also after 24h of C5a exposure followed by 24h of normal culture. Increase in the gene expression of each pro-inflammatory cytokine was found also by 300 μM H2O2, that is a stressor inducer of SASP (Figure 6C6F).

Finally, in order to gain insight in BCL9, CYP1B1 and CDK6 genes expression we performed qPCR at different time points. All genes showed a significant up-regulation both after a short C5a stimulation (3h) followed by 24 or 48h of culture in normal medium and after 24h of C5a stimulation followed by 24h of culture. A similar increase was detected after 300 μM H2O2 exposure (Figure 6G6I).

Taken together these results demonstrate that C5a can induce tubular senescence in vitro, characterized by morphological and functional changes such as growth arrest and acquirement of a pro-inflammatory phenotype named SASP.

C1-INH modulated Wnt pathway activation during I/R leads counteracting renal senescence

We further investigated whether I/R might promote the development of tubular senescence in vivo. To this aim, we used a swine model of renal I/R injury in which pigs were treated with C1-Inhibitor to modulate Complement activation [28, 29]. We first investigated Wnt4 expression on renal sections obtained at 24 hours after I/R injury since the in vitro data showed an aberrant methylation and increased expression of FZD6 (Supplementary Figure 4B, 4D), a well-known Wnt4 receptor [30].

During I/R injury, the protein level of Wnt4 increased compared to T0 (Figure 7A). The expression mainly localized to the apical side (Figure 7A, T24 CTRL arrow) of epithelial tubular cells and was more evident in area of tubular damage characterized by vacuolization and interstitial edema. Moreover, Wnt4 expression was significantly reduced by treatment with C1 inhibitor (Figure 7A, T24 C1-INH). Subsequently, we assessed the Wnt/β-catenin pathway activation by measuring the level of β-catenin activation. One day after I/R injury, an increased number of β-catenin+ cells were detected (Figure 7B). Compared with the T0, β- catenin protein was up regulated predominantly in renal tubules of the I/R injured kidney. β-catenin was localized in the cytoplasm and the nuclei of tubular epithelial cells. (Figure 7B, T24 CTRL). After 24h from the reperfusion, relative β-catenin levels were increased compared to the T0, suggesting that induction of Wnt expression would result in an accumulation of β-catenin in injured and senescent kidney. Less β-catenin+ cells were present in C1 inhibitor-treated pigs (Figure 7B, T24 C1-INH), again showing the involvement of Complement in the activation of this pathway.

Wnt4/βcatenin pathway and inflammaging markers are activated in tubular cells after I/R and modulated by C1-INH treatment. (A–B) Immunoistochemical stainings showing the tubular Wnt4 and βcatenin increase after 1 day of I/R injury and the C1-INH-mediated modulation. IHC was performed on paraffin kidney sections. Arrows indicate positive tubular staining. (in the right) Graphical representation of Wnt4 and βcatenin protein expression level in the different groups. (n=5, p value as indicated, scale bar as indicated).(C) Representative SA-β Gal stained kidney tissues revealed an higher number of senescent cells after 24 h of I/R injury compared to T0. Treatment with C1-INH restored SA-β Gal at basal expression. Arrows indicate positive tubular staining on cryo tissues. (in the right) Graphical representation %SA-βGal area fraction. (n=5, p value as indicated, scale bar as indicated). (D) Representative micrographs indicating p16INK4a protein expression in different groups of swine, as indicated. Boxed areas are enlarged in the bottom of each micrographs. In the T0, p16INK4a had constitutive level and was localized in tubular nuclei (black arrow). Biopsies after 24h from reperfusion (T24 CTRL) showed increased nuclear (black arrow) and cytoplasmic staining (white dotted arrows). C1-INH restored p16INK4a at basal expression; limiting the cytoplasmic p16INK4a expression. (in the right) Graphical representation of p16INK4a area fraction in the different groups. (n=5, p value as indicated, scale bar as indicated).

Figure 7. Wnt4/βcatenin pathway and inflammaging markers are activated in tubular cells after I/R and modulated by C1-INH treatment. (AB) Immunoistochemical stainings showing the tubular Wnt4 and βcatenin increase after 1 day of I/R injury and the C1-INH-mediated modulation. IHC was performed on paraffin kidney sections. Arrows indicate positive tubular staining. (in the right) Graphical representation of Wnt4 and βcatenin protein expression level in the different groups. (n=5, p value as indicated, scale bar as indicated).(C) Representative SA-β Gal stained kidney tissues revealed an higher number of senescent cells after 24 h of I/R injury compared to T0. Treatment with C1-INH restored SA-β Gal at basal expression. Arrows indicate positive tubular staining on cryo tissues. (in the right) Graphical representation %SA-βGal area fraction. (n=5, p value as indicated, scale bar as indicated). (D) Representative micrographs indicating p16INK4a protein expression in different groups of swine, as indicated. Boxed areas are enlarged in the bottom of each micrographs. In the T0, p16INK4a had constitutive level and was localized in tubular nuclei (black arrow). Biopsies after 24h from reperfusion (T24 CTRL) showed increased nuclear (black arrow) and cytoplasmic staining (white dotted arrows). C1-INH restored p16INK4a at basal expression; limiting the cytoplasmic p16INK4a expression. (in the right) Graphical representation of p16INK4a area fraction in the different groups. (n=5, p value as indicated, scale bar as indicated).

To further clarify the final effect of the Wnt activation, we analyzed the tubular senescence by SA β-gal, p16 and p21 staining. As shown in Figures 7C and 7D and 8A (T24 CTRL), the Wnt/β-catenin activation was associated with SA β-gal positivity and increase in tubular p16 and p21 24h after I/R injury, Regards to p16 and p21, we detected beyond the constitutive nuclear staining (black arrow) (Figures 7D, 8A, T24 CTRL], a strong and diffuse cytoplasmic staining (white dotted arrow) that was described to have a negative clinical significance [3133]. Furthermore, to assess the acquirement of SASP in vivo, we stained for IL-6 (Figure 8B) and found a statistically significant increase, as shown in Figure 8B and 8D.

p21 and IL-6 are induced in tubular cells after I/R and modulated by C1-INH treatment. (A) Representative micrographs indicating p21/CIP1 protein expression in different groups of swine after I/R injury, as indicated. Boxed areas are enlarged at the bottom of each micrographs. Compared to T0, biopsies after 24h from reperfusion (T24 CTRL) showed increased nuclear (black arrow) and cytoplasmic staining (white dotted arrow). C1-INH treatment restored p21 at basal expression; limiting the cytoplasmic p21 expression. (B) Representative micrographs show the expression and localization of IL-6, a marker of SASP in T0 (right), T24 CTRL (middle) and T24 C1-INH (right) groups. Arrows indicate positive tubular cells. Compared to T0, biopsies after 24h from reperfusion (T24 CTRL) showed increased IL-6 expression, predominantly at tubular cells. Treatment with C1-INH counteracts IL-6 increase. (C, D) Graphical representation of p21 and IL-6 area fraction in the different groups. (n=5, p value as indicated).

Figure 8. p21 and IL-6 are induced in tubular cells after I/R and modulated by C1-INH treatment. (A) Representative micrographs indicating p21/CIP1 protein expression in different groups of swine after I/R injury, as indicated. Boxed areas are enlarged at the bottom of each micrographs. Compared to T0, biopsies after 24h from reperfusion (T24 CTRL) showed increased nuclear (black arrow) and cytoplasmic staining (white dotted arrow). C1-INH treatment restored p21 at basal expression; limiting the cytoplasmic p21 expression. (B) Representative micrographs show the expression and localization of IL-6, a marker of SASP in T0 (right), T24 CTRL (middle) and T24 C1-INH (right) groups. Arrows indicate positive tubular cells. Compared to T0, biopsies after 24h from reperfusion (T24 CTRL) showed increased IL-6 expression, predominantly at tubular cells. Treatment with C1-INH counteracts IL-6 increase. (C, D) Graphical representation of p21 and IL-6 area fraction in the different groups. (n=5, p value as indicated).

As shown in Figure 7A and 7B, Complement inhibition significantly abolished the Wnt/β-catenin activation induced by I/R. Consistently with the decreased activation of Wnt/β-catenin pathway by C1-inhibitor, the reduction of Wnt4 was associated to a block of SA-β gal enzymatic tubular activity (Figure 7C), to a reduction of p16 and p21 cytoplasmic expression (Figures 7D, 8A and 8C) and to a downregulation of IL-6 production. (Figure 8B and 8D). These results provide in vivo proof that Complement modulation can reduce tubular senescence by regulating Wnt/β-catenin pathway and inducing SASP in vivo.

Discussion

In this study, we provide evidences that Complement anaphylatoxin C5a significantly triggers DNA methylation changes in renal tubular epithelial cells. C5a epigenetically induced the expression of genes involved in the Wnt/βcatenin pathway promoting the acquirement of a premature senescence-related phenotype called SASP. We confirmed the occurrence of this pathogenic process in vivo in a swine model of renal I/R injury where the treatment by C1-INH could inhibit the Wnt/βcatenin signaling and the occurrence of tubular Inflammaging.

It is known that Complement is the driving force in the pathophysiology of I/R injury-induced AKI [34, 35] and C5a, the most powerful anaphylatoxin, is able to bind two kind of receptors (C5aR1 and C5aR2), both expressed on renal tubular epithelial cells [36, 37]. Recently, a growing body of evidence suggested that epigenetic mechanisms are involved during acute injury in the transition from AKI to CKD, [13, 3840]; therefore we reasoned that C5a might have an effect on the DNA methylation profile in tubular cells by regulating the expression of inflammatory or injury factors after I/R. In our study, in order to extensively cover the whole human genome, we analyzed all CpGs regions into promoters, 5’UTR, 3’UTR gene exons and introns. Our observations that C5a induced the hypomethylation of 88 DNA regions and the hypermethylation of 54 DNA sequences are in line with data from literature showing a significant correlation between DNA methylation changes in kidney donors after I/R and the risk to progression to CKD [39]. Interestingly, we found that C5a exposure induced a global genome-wide hypo-methylation (Figure 1B), a DNA pattern that has already been observed in human tissues during aging by several groups [4144]. This conclusion could be supported by recent evidences showing that Complement has a central role in several aged related disease [45, 46].

Between all the pathways involved in C5a-stimulated RTEC, we showed that Complement changes the methylation level of genes involved in Wnt/βcatenin pathway, the DNA damage response, the cell cycle checkpoint regulation and in particular the Aryl Hydrocarbon Receptor Signaling. Concerning Wnt/βcatenin pathway, [47] it has been demonstrated its significant role in the aging kidney [48, 49] in several experimental models of chronic fibrosis [5052]. Interestingly, DiRocco Derek P et al, demonstrated that Wnt/βcatenin signaling constitutive activation was sufficient to drive spontaneous myofibroblast differentiation in absence of injury [53], with Wnt4 playing a pivotal role in chronic fibrosis [54].

By network analysis (Supplementary Figure 2), we determined the biological and functional connectivity of all differentially methylated genes and we were able to identify and validate (Figures 2; 3) three main genes as central nodes in the network: BCL9, CDK6 and CYP1B1. These hypomethylated genes functionally belonged to Wnt/βcatenin pathway. BCL9 is an essential co-factor in the Wnt/β-catenin [55]; in addition, the protein kinase CDK6 that regulates the cell cycle is under the control of Wnt [56]. Our findings are in line with previous data that showed firstly a pivotal role of Wnt signaling during AKI [5760] and secondly a new pathogenic and pro-inflammatory effect of CDK6 during renal I/R injury [61]. Finally, also CYP1B1, a cytochrome P450 monooxygenase involved in estrogen metabolism, can be regulated by β-catenin [62] and can activate Wnt signaling [6365]. Nevertheless, the CYP1B1 finding, that acts as a well-known oncogenic protein in renal cell carcinoma by the Aryl Hydrocarbon Receptor Signaling suggested an unexplored link between aberrant Complement activation, DNA damage and tumorigenesis [66]. Even if it is known that the C5a/C5aR1 pathway can mediate tumorigenesis in several animal models [6669], more studies are required to elucidate how Complement activation and C5aR signaling could differentially lead to premature senescence by Wnt pathway or increase susceptibility to cancer after transplantation.

The double presence of Wnt pathway (Supplementary Figure 1) in the most significant biological processes and its pivotal involvement in renal aging [49] led us to investigate the tubular senescence. Cellular senescence is defined as a state of irreversible cell cycle arrest, DNA damage, apoptosis inhibition and a persistent secretion of pro-inflammatory cytokines, which is referred as the Senescence-Associated Secretory Phenotype (SASP) [70]. The growth arrest is maintained by an increase of cell cycle regulators as p16INK4A, p21CIP1 and p53 [49, 71], the SASP by the cytokines IL-6, MCP-1, CTGF and PAI-1 [72, 73]. Our results demonstrated that Complement is involved in premature senescence in tubular cells and that this event is induced after a short period of stimulation and then became permanent and stable even after the anaphylatoxin removal. This is consistent with results showing the predominant role of C5a in the early phases of I/R injury [74] and with the chronic effect of Complement activation after renal transplantation. In particular, by SA-β-gal assay we found the C5aR1 inhibition was able to attenuate the acquirement of a senescence phenotype, modulating also the C5a-induced enlargement of cell size. These results are in line with several findings showing that C5aR1 inhibition protects kidneys from I/R injury [75, 76].

Recent studies have revealed that Wnt signaling leading to mammalian aging can be induced by Complement component C1q [24] and antagonized by the protein Klotho [77, 78]. Furthermore, Complement severely reduced Klotho level in tubular cells [29]. Thus, we hypothesized that Complement modulation could regulate this vicious cycle that lead to tubular senescence and graft failure. C1-INH is a Complement inhibitor safely used as therapy for hereditary angioedema and antibody-mediated rejection [79] (NCT 01134510, NCT 01035593, NCT 01147302). Our in vivo studies confirmed that renal I/R activated Wnt/β-catenin signaling [52, 80] and that accelerates renal aging as showed by p16 and p21 higher level (Figures 7D, 8A T24 CTRL versus T0) [32, 81, 82]. The interesting observation is that C1-INH treatment abrogated Wnt4 and reduced the β-catenin activation. This inhibition led to a reduced tubular senescence as observed by decreased SA-βgal, p16, p21 and IL-6 positivity. These data suggest a consequential correlation between Complement activation, the epigenetically induced Wnt signaling and tubular senescence after I/R injury. Considering the strong link between senescence and fibrosis as long-term consequences of renal graft deterioration, these findings are supported by several studies showing anti-fibrotic role of C1-INH [28, 83, 84].

We recognize that overall our study has some limitations. The methylation data used included only 3 different lots of human RTEC and further investigation is required to increase the number of samples. However, by analyzing one single renal population we obtained results that are not affected by cellular heterogeneity since DNA methylation is cell-type specific. In addition, we evaluated the Wnt pathway modulation in a swine model of renal I/R injury. Even if pigs represents one of the best large animal models for kidney transplantation [85], in a next future human renal biopsies should be investigated.

In conclusion, with the present study we demonstrate for the first time, that C5a can modulate the DNA methylation of specific genes in RTEC and that following Complement exposure, the Wnt/βcatenin pathway was significantly activated, leading to a senescent phenotype. Therefore, understanding of the epigenetic mechanisms underlying the AKI-accelerated renal senescence could pave the way for the use of Complement inhibitors to counteract the progression from AKI to CKD.

Materials and Methods

RTEC culture and DNA/RNA extraction

RTEC (human Renal Proximal Tubular Epithelial Cells) between 3-5 passages were grown to confluence in RGEM medium (Renal epithelial cell growth medium, Lonza) in which a basal medium RBEM was supported with rhEGF (recombinant human EGF), transferrin, insulin, hydrocortisone, epinephrine, triidothyronine, and Fetal Bovine Serum (FBS), preserved at +4°C and changed every two days in the culture. Cells with medium were incubated at 37°C with CO2 at 5%. Once they became confluent, cells were exposed to C5a anaphylotoxin (10−7 M) at different times. After C5a stimulation, cells were washed 3 times in PBS1X and fresh medium was replaced. In addition, cells were exposed to H2O2 (100-300 μm) as positive control of senescence. For C5aR1 inhibition, mouse monoclonal anti-C5aR (Abcam) was pre-incubated for 1h before the C5a exposure. DNA and RNA were simultaneously extracted from RTEC using the AllPrep DNA/RNA Mini Kit (Qiagen) according to the manufacturer’s protocol. Nucleic acid concentration and quality were assessed using a NanoDrop Spectrophotometer (NanoDrop Technologies).

DNA methylation analysis

DNA methylation analysis was performed utilizing the Illumina methylation profiling platform consisting of HiScanSQ system and Infinium HumanMethylation450 BeadChips. Bisulfite conversion of DNA extracted from RTEC was performed with the EZ-96 DNA Methylation Kit (ZymoResearch), and subsequent hybridization of this DNA was carried out on the Infinium Human Methylation 450 Bead Chips. Microarray data of the C5a stimulated RTEC and normal cultured RTEC are available under accession number GSE115227 at the Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/). All statistical analysis was performed using R and the RnBeads R package. Differential methylation analysis on site and region level was computed considering the difference in mean methylation levels C5a stimulated RTEC and normal cultured RTEC, the quotient in mean methylation and the t test assessing whether the methylation values in the two groups originate from distinct distributions. Additionally each site was assigned a rank based on each of these three criteria. A combined rank was computed as the maximum (i.e. worst) rank among the three ranks. Coverage of differential methylation analysis on region level was targeted across gene regions with sites in the promoter region, 5’-UTR, first exon, gene body and 3’-UTR in order to provide the broadest most comprehensive view of methylation state possible.

Pyrosequencing assay

We validated the methylation status of differentially methylated CpG sites of CDK6, BAX, CYP1B1-AS1, BCL9, FZD6, CCDC6 genes using the pyrosequencing technology. The pyrosequencing assays (two PCR primers to amplify a specific target region and a sequencing primer) were designed with the Pyrosequencing Assay Design software (Biotage). The assays were set up sequencing control DNA samples chemically methylated and unmethylated respectively (EpiTect PCR Control DNA Set, Qiagen). The PCRs were performed using ImmolaseTM DNA Polymerase (Bioline, Aurogene) according to the manufacturer’s instructions. The amplicons of the specific targets were analyzed with a PSQ 96MA instrument (Qiagen). Sequencing reactions were performed with the Pyro Gold reagent kit PSQ 96MA (Qiagen) according to the manufacturer’s instructions and the sequencing analysis was conducted with the PSQTM 96MA software (version 2.02). The sequence of all primers, the PCR amplification conditions, and the sequence of the pyrosequenced fragments are indicated in Supplementary Table 3.

Quantitative methylation-specific PCR assay (qMSP)

Bisulfite converted DNA was used to perform the qMSP by Methylamp MS-qPCR Fast Kit (Epigentek Group) according the manufacturer instructions. For each reaction 20 ng of bisulfite-treated DNA was used as template. Primer design and sequences and reaction methods are detailed in the Supplementary Methods section. Primers for the genes of interest were designed using MethPrimer (http://www.urogene.org/methprimer/index.html) (Supplementary Table 1). ACTB (β-actin) gene was used as a reference gene. No-template controls were included in each run as negative controls. An EpiTect Control DNA, a 100% methylated DNA (Qiagen), was used as a positive control for all genes studied. The PMR (percentage of methylated reference) (i.e. degree of methylation) was used to define the percentage of fully methylated molecules at a specific locus and was calculated as reported previously. Briefly, the PMR value was calculated by dividing the gene/ACTB ratio in a sample by the gene/ACTB ratio in SssI-treated leucocyte DNA (Qiagen) and multiplied by 100. Parallel PCRs were carried out for the genes of interest and reference. PMR values were detected using the comparative CT method. The relationship between the percentages of methylated DNA molecules and CT is described as PMR = 2−DDCT×100%.

qPCR

qPCR was carried out with SsoAdvanced™ Universal SYBR® Green Supermix (Biorad) and the Light Cycler@96 (Roche). Cycling conditions and primer list sequence in Supplementary Table 4.

SA-β-gal test

Senescence-associated SA-β-gal staining was performed as described in manufacturer's protocol (Cell Signaling Technology). After washing with PBS (pH 6.0), the cells were fixed in 4% paraformaldehyde for 15 minutes and stained with freshly prepared SA- β –gal solution (1 mg/mL X-gal, 40 mM citric acid/sodium phosphate (pH 6.0), 5 mM potassium ferrocyanide, 5 mM potassium ferricyanide, 150 mM NaCl, and 2 mM MgCl). Next, cells are incubated at 37 °C overnight in a dry incubator. The, we removed the staining solution, immersed the samples in 70% glycerol, and assessed the development of blue color. The number of SA-β-gal-positive cells was determined by counting at least 500 cells per dish (60mm) in five not overlapping fields. The % of the number of SA-β-gal strong positive cells that appeared with increased dimension is calculated as ratio considering the total numbers of cells counted in each dish.

MTT assay

Primary epithelial tubular cells were incubated with and without human recombinant C5a for 3 or 24h followed by 24, 48 h of normal culture after medium culture change. Cultured RTEC proliferation was measured by MTT Cell Proliferation Assay Kit, according to the manufacturer instructions (Sigma Aldrich). Briefly, 3×104cells/well were seeded in a 96-well plate, and then cells were treated with C5a as indicated. Absorbance at 570 nm was then measured by a spectrophotometer.

Western blot

Protein lysates were homogenized by RIPA buffer with phosphatase and protease inhibitors. Proteins (30 μg) were separated in 4–15% polyacrylamide gel and then transferred to PVDF membrane (0.2mM) by Trans-Blot Turbo (BioRad, Hercules, CA). After blocking in BSA at 5%, the membranes were incubated overnight with the following primary antibodies: p21 (Abcam), p16INK4a (Abcam) and then with secondary antibody (hrp-conjugated, Santa Cruz). The same membrane was probed with mouse monoclonal anti-βactin antibody (1:20000; Sigma). The Electrochemiluminescence (ECL) system was used to detect the antibody binding, (Amersham, UK). The chemiluminescent signal was acquired by Chemidoc and quantified using Image J software.

Statistical analysis

Graphs were displayed using GraphPad Prism Software 5. Data were expressed as median ± interquartile range (IQR) and compared with a Mann–Whitney test for tissue immunostainings. For qPCR, MTT and WB data were expressed as the mean± SD. Statistical analysis was assessed using unpaired Student’t-test. A p value of <0.05 was considered significant.

Swine model of I/R injury

Animal studies were carried out under protocol approved by Ethical Committee of the Italian Ministry of Health. Briefly, I/R was induced by clamping the renal artery for 30min followed by reperfusion, as described previously. A biopsy was performed before ischemia (T0). Pigs were divided into two groups: control (CTRL, n=5, vehicle infused) and C1 Inhibitor treated group (C1-INH, n=5). Five minutes before the beginning of the reperfusion, rC1-INH was injected in the ear vein (500 U/kg). Biopsies were performed at 15, 30, 60 min and 24h after reperfusion. All animals were sacrificed 24 h after the procedure. Controlateral kidney was not removed.

Immunohistochemistry

Renal sections underwent deparaffination and heat-mediated antigen retrieval (citrate buffer, pH=6.00) as previously described. For p16, p21 and βcatenin detection, sections were permeabilized with Triton 0.25% for 5 min, then blocked by Protein Block Solution (DakoCytomation, USA) for 10 min. Incubation was performed with antibodies against: p16, p21 (Abcam, Cambridge UK), IL-6 (Novus Biologicals), wnt4 and βcatenin. Then, the positive staining was detected by the Peroxidase/DAB Dako Real EnVision Detection System (Dako, Glostrup, Denmark). The peroxidase reaction was shown by a brown precipitate, counterstained with Mayers hematoxylin (blue). Negative controls were prepared by incubation with a control irrelevant antibody. Images were scanned by Aperio ScanScope CS2 device and signals were analyzed with the ImageScope V12.1.0.5029 (Aperio Technologies, Vista, CA).

Cell cycle analysis

After treatment with C5a for 3h, 24h, 48h, RTEC cells were fixed in 70% ethanol then treated with 50 mg of RNase A/mL at 37 °C for 30 minutes. Aliquots of 300.000 cells were added to 0.2 mL of propidium iodide solution (from 50 μg/ml stock solution) with 0.1% sodium citrate and 0.1% Nonidet P-40) and incubated for 30 minutes in the dark at 4°C. The cellular DNA content was analyzed by flow cytometry FC500 flow cytometer (Beckmann Coulter) and analyzed by Kaluza software.

Abbreviations

AhR: Aryl Hydrocarbon Receptor Signaling; AKI: Acute Kidney Injury; Bas: Basal condition; BCL9: B-Cell CLL/Lymphoma 9; C1-INH: C1-Inhibitor; CCDC6: Coiled-Coil Domain Containing 6; CKD: Chronic Kidney Disease; CTRL: Control group; CYP1B1: Cytochrome P450 Family 1 Subfamily B Member 1; DNMT: DNA methyl transferase; FDZ6: Frizzled Class Receptor 6; GAPDH: Glyceraldehyde 3-phosphate dehydrogenase; I/R: Ischemia/Reperfusion; RTEC: Renal Tubular Epithelial Cells; SASP: Senescence-Associated Secretory Phenotype; SA-βGAL: Senescence-Associated-βgalactosidase; SD: Standard Deviation; SERPINE1 (PAI-1): Serpin Family E Member 1 (Plasminogen Activator Inhibitor, Type I); TP53: Tumor Protein P53; qMSP: Quantitative Methylation Specific PCR; 5-aza: 5-aza-2’-deoxycytidine or Decitabine.

Author Contributions

Giuseppe Castellano, Rossana Franzin, Fabio Sallustio and Alessandra Stasi performed all the experiments and statistical analyses, Giuseppe De Palma performed MSP-qPCR validation analysis, Barbara Banelli and Massimo Romani performed the pyrosequencing assay and analysis, Giuseppe Lucarelli, Michele Battaglia, Antonio Crovace, Francesco Staffieri, Pasqual Ditonno performed the swine model of IRI, Chiara Divella collected and provided the pig samples, Giuseppe Grandaliano, Paolo Cravedi, Vincenzo Cantaluppi and Loreto Gesualdo supervised all the project. All authors were involved in data interpretation. Giuseppe Castellano and Rossana Franzin wrote the paper, all authors gave final approval of the submitted version.

Acknowledgments

This work was supported by University of Bari ‘‘Aldo Moro’’ and Ministero della Salute (Ricerca Finalizzata 2009, GR-2009-1608662, granted to G. C. and G. G.). The authors thank Angelica Intini and Angela Picerno for excellent technical assistance. The authors thank Beatris Oortwijn and Edwin van Amersfoort from Pharming Group NV, Leiden, the Netherlands for the unrestricted research grant supporting the study.

Conflicts of Interest

The authors of this manuscript have no conflicts of interest to disclose.

References

  • 1. Kellum JA, Prowle JR. Paradigms of acute kidney injury in the intensive care setting. Nat Rev Nephrol. 2018; 14:217–30. https://doi.org/10.1038/nrneph.2017.184 [PubMed]
  • 2. Kellum JA, Bellomo R, Ronco C. Kidney attack. JAMA. 2012; 307:2265–66. https://doi.org/10.1001/jama.2012.4315 [PubMed]
  • 3. Barasch J, Zager R, Bonventre JV. Acute kidney injury: a problem of definition. Lancet. 2017; 389:779–81. https://doi.org/10.1016/S0140-6736(17)30543-3 [PubMed]
  • 4. Tang J, Zhuang S. Epigenetics in acute kidney injury. Curr Opin Nephrol Hypertens. 2015; 24:351–58. [PubMed]
  • 5. Wanner N, Bechtel-Walz W. Epigenetics of kidney disease. Cell Tissue Res. 2017; 369:75–92. https://doi.org/10.1007/s00441-017-2588-x [PubMed]
  • 6. Parker MD, Chambers PA, Lodge JP, Pratt JR. Ischemia- reperfusion injury and its influence on the epigenetic modification of the donor kidney genome. Transplantation. 2008; 86:1818–23. https://doi.org/10.1097/TP.0b013e31818fe8f9 [PubMed]
  • 7. Boer K, de Wit LE, Peters FS, Hesselink DA, Hofland LJ, Betjes MG, Looman CW, Baan CC. Variations in DNA methylation of interferon gamma and programmed death 1 in allograft rejection after kidney transplantation. Clin Epigenetics. 2016; 8:116. https://doi.org/10.1186/s13148-016-0288-0 [PubMed]
  • 8. Smyth LJ, McKay GJ, Maxwell AP, McKnight AJ. DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease. Epigenetics. 2014; 9:366–76. https://doi.org/10.4161/epi.27161 [PubMed]
  • 9. Chu AY, Tin A, Schlosser P, Ko YA, Qiu C, Yao C, Joehanes R, Grams ME, Liang L, Gluck CA, Liu C, Coresh J, Hwang SJ, et al. Epigenome-wide association studies identify DNA methylation associated with kidney function. Nat Commun. 2017; 8:1286. https://doi.org/10.1038/s41467-017-01297-7 [PubMed]
  • 10. Chen J, Zhang X, Zhang H, Lin J, Zhang C, Wu Q, Ding X. Elevated Klotho promoter methylation is associated with severity of chronic kidney disease. PLoS One. 2013; 8:e79856. https://doi.org/10.1371/journal.pone.0079856 [PubMed]
  • 11. Peters FS, Manintveld OC, Betjes MG, Baan CC, Boer K. Clinical potential of DNA methylation in organ transplantation. J Heart Lung Transplant. 2016; 35:843–50. https://doi.org/10.1016/j.healun.2016.02.007 [PubMed]
  • 12. Shiels PG, McGuinness D, Eriksson M, Kooman JP, Stenvinkel P. The role of epigenetics in renal ageing. Nat Rev Nephrol. 2017; 13:471–82. https://doi.org/10.1038/nrneph.2017.78 [PubMed]
  • 13. Bechtel W, McGoohan S, Zeisberg EM, Müller GA, Kalbacher H, Salant DJ, Müller CA, Kalluri R, Zeisberg M. Methylation determines fibroblast activation and fibrogenesis in the kidney. Nat Med. 2010; 16:544–50. https://doi.org/10.1038/nm.2135 [PubMed]
  • 14. Zhang Q, Yin S, Liu L, Liu Z, Cao W. Rhein reversal of DNA hypermethylation-associated Klotho suppression ameliorates renal fibrosis in mice. Sci Rep. 2016; 6:34597. https://doi.org/10.1038/srep34597 [PubMed]
  • 15. Zhao Y, Ding C, Xue W, Ding X, Zheng J, Gao Y, Xia X, Li S, Liu J, Han F, Zhu F, Tian P. Genome-wide DNA methylation analysis in renal ischemia reperfusion injury. Gene. 2017; 610:32–43. https://doi.org/10.1016/j.gene.2017.02.005 [PubMed]
  • 16. Vervaet BA, Moonen L, Godderis L, Poels K, D’Haese PC. Untargeted DNA-Demethylation Therapy Neither Prevents Nor Attenuates Ischemia-Reperfusion-Induced Renal Fibrosis. Nephron. 2017; 137:124–36. https://doi.org/10.1159/000477507 [PubMed]
  • 17. Tatapudi VS, Montgomery RA. Pharmacologic Complement Inhibition in Clinical Transplantation. Curr Transplant Rep. 2017; 4:91–100. https://doi.org/10.1007/s40472-017-0148-7 [PubMed]
  • 18. Legendre C, Sberro-Soussan R, Zuber J, Frémeaux-Bacchi V. The role of complement inhibition in kidney transplantation. Br Med Bull. 2017; 124:5–17. https://doi.org/10.1093/bmb/ldx037 [PubMed]
  • 19. Fayyazi A, Scheel O, Werfel T, Schweyer S, Oppermann M, Götze O, Radzun HJ, Zwirner J. The C5a receptor is expressed in normal renal proximal tubular but not in normal pulmonary or hepatic epithelial cells. Immunology. 2000; 99:38–45. https://doi.org/10.1046/j.1365-2567.2000.00911.x [PubMed]
  • 20. Zhang K, Li GQ, He QH, Li Y, Tang M, Zheng QY, Xu GL, Zhang KQ. C5a/C5aR pathway accelerates renal ischemia-reperfusion injury by downregulating PGRN expression. Int Immunopharmacol. 2017; 53:17–23. https://doi.org/10.1016/j.intimp.2017.10.006 [PubMed]
  • 21. Liu F, Gou R, Huang J, Fu P. [The effect of C3a, C5a and their receptor antagonists on the expression of beta-catenin in renal tubular cells]. Sichuan Da Xue Xue Bao Yi Xue Ban. 2011; 42:74–7. [PubMed]
  • 22. Thompson B, Townsley F, Rosin-Arbesfeld R, Musisi H, Bienz M. A new nuclear component of the Wnt signalling pathway. Nat Cell Biol. 2002; 4:367–73. https://doi.org/10.1038/ncb786 [PubMed]
  • 23. Kawasaki Y, Komiya M, Matsumura K, Negishi L, Suda S, Okuno M, Yokota N, Osada T, Nagashima T, Hiyoshi M, Okada-Hatakeyama M, Kitayama J, Shirahige K, Akiyama T. MYU, a Target lncRNA for Wnt/c-Myc Signaling, Mediates Induction of CDK6 to Promote Cell Cycle Progression. Cell Rep. 2016; 16:2554–64. https://doi.org/10.1016/j.celrep.2016.08.015 [PubMed]
  • 24. Naito AT, Sumida T, Nomura S, Liu ML, Higo T, Nakagawa A, Okada K, Sakai T, Hashimoto A, Hara Y, Shimizu I, Zhu W, Toko H, et al. Complement C1q activates canonical Wnt signaling and promotes aging-related phenotypes. Cell. 2012; 149:1298–313. https://doi.org/10.1016/j.cell.2012.03.047 [PubMed]
  • 25. Lessel D, Wu D, Trujillo C, Ramezani T, Lessel I, Alwasiyah MK, Saha B, Hisama FM, Rading K, Goebel I, Schütz P, Speit G, Högel J, et al. Dysfunction of the MDM2/p53 axis is linked to premature aging. J Clin Invest. 2017; 127:3598–608. https://doi.org/10.1172/JCI92171 [PubMed]
  • 26. Georgakilas AG, Martin OA, Bonner WM. p21: A Two-Faced Genome Guardian. Trends Mol Med. 2017; 23:310–19. https://doi.org/10.1016/j.molmed.2017.02.001 [PubMed]
  • 27. Tchkonia T, Zhu Y, van Deursen J, Campisi J, Kirkland JL. Cellular senescence and the senescent secretory phenotype: therapeutic opportunities. J Clin Invest. 2013; 123:966–72. https://doi.org/10.1172/JCI64098 [PubMed]
  • 28. Castellano G, Franzin R, Stasi A, Divella C, Sallustio F, Pontrelli P, Lucarelli G, Battaglia M, Staffieri F, Crovace A, Stallone G, Seelen M, Daha MR, et al. Complement Activation During Ischemia/Reperfusion Injury Induces Pericyte-to-Myofibroblast Transdifferentiation Regulating Peritubular Capillary Lumen Reduction Through pERK Signaling. Front Immunol. 2018; 9:1002. https://doi.org/10.3389/fimmu.2018.01002 [PubMed]
  • 29. Castellano G, Intini A, Stasi A, Divella C, Gigante M, Pontrelli P, Franzin R, Accetturo M, Zito A, Fiorentino M, Montinaro V, Lucarelli G, Ditonno P, et al. Complement Modulation of Anti-Aging Factor Klotho in Ischemia/Reperfusion Injury and Delayed Graft Function. Am J Transplant. 2016; 16:325–33. https://doi.org/10.1111/ajt.13415 [PubMed]
  • 30. Heinonen KM, Vanegas JR, Lew D, Krosl J, Perreault C. Wnt4 enhances murine hematopoietic progenitor cell expansion through a planar cell polarity-like pathway. PLoS One. 2011; 6:e19279. https://doi.org/10.1371/journal.pone.0019279 [PubMed]
  • 31. McGlynn LM, Stevenson K, Lamb K, Zino S, Brown M, Prina A, Kingsmore D, Shiels PG. Cellular senescence in pretransplant renal biopsies predicts postoperative organ function. Aging Cell. 2009; 8:45–51. https://doi.org/10.1111/j.1474-9726.2008.00447.x [PubMed]
  • 32. Melk A, Schmidt BM, Vongwiwatana A, Rayner DC, Halloran PF. Increased expression of senescence-associated cell cycle inhibitor p16INK4a in deteriorating renal transplants and diseased native kidney. Am J Transplant. 2005; 5:1375–82. https://doi.org/10.1111/j.1600-6143.2005.00846.x [PubMed]
  • 33. Sis B, Tasanarong A, Khoshjou F, Dadras F, Solez K, Halloran PF. Accelerated expression of senescence associated cell cycle inhibitor p16INK4A in kidneys with glomerular disease. Kidney Int. 2007; 71:218–26. https://doi.org/10.1038/sj.ki.5002039 [PubMed]
  • 34. Sacks SH, Zhou W. The role of complement in the early immune response to transplantation. Nat Rev Immunol. 2012; 12:431–42. https://doi.org/10.1038/nri3225 [PubMed]
  • 35. Ricklin D, Reis ES, Lambris JD. Complement in disease: a defence system turning offensive. Nat Rev Nephrol. 2016; 12:383–401. https://doi.org/10.1038/nrneph.2016.70 [PubMed]
  • 36. Manthey HD, Woodruff TM, Taylor SM, Monk PN. Complement component 5a (C5a). Int J Biochem Cell Biol. 2009; 41:2114–17. https://doi.org/10.1016/j.biocel.2009.04.005 [PubMed]
  • 37. Ricklin D, Mastellos DC, Reis ES, Lambris JD. The renaissance of complement therapeutics. Nat Rev Nephrol. 2018; 14:26–47. https://doi.org/10.1038/nrneph.2017.156 [PubMed]
  • 38. Pratt JR, Parker MD, Affleck LJ, Corps C, Hostert L, Michalak E, Lodge JP. Ischemic epigenetics and the transplanted kidney. Transplant Proc. 2006; 38:3344–46. https://doi.org/10.1016/j.transproceed.2006.10.112 [PubMed]
  • 39. Heylen L, Thienpont B, Naesens M, Busschaert P, Depreeuw J, Smeets D, Jochmans I, Monbaliu D, Pirenne J, Lerut E, Ghesquiere B, Kuypers D, Lambrechts D, Sprangers B. Ischemia-Induced DNA Hypermethylation during Kidney Transplant Predicts Chronic Allograft Injury. J Am Soc Nephrol. 2018; 29:1566–76. https://doi.org/10.1681/ASN.2017091027 [PubMed]
  • 40. Rodríguez-Romo R, Berman N, Gómez A, Bobadilla NA. Epigenetic regulation in the acute kidney injury to chronic kidney disease transition. Nephrology (Carlton). 2015; 20:736–43. https://doi.org/10.1111/nep.12521 [PubMed]
  • 41. Day K, Waite LL, Thalacker-Mercer A, West A, Bamman MM, Brooks JD, Myers RM, Absher D. Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape. Genome Biol. 2013; 14:R102. https://doi.org/10.1186/gb-2013-14-9-r102 [PubMed]
  • 42. Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, Beyan H, Whittaker P, McCann OT, Finer S, Valdes AM, Leslie RD, Deloukas P, Spector TD. Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010; 20:434–39. https://doi.org/10.1101/gr.103101.109 [PubMed]
  • 43. Johansson A, Enroth S, Gyllensten U. Continuous Aging of the Human DNA Methylome Throughout the Human Lifespan. PLoS One. 2013; 8:e67378. https://doi.org/10.1371/journal.pone.0067378 [PubMed]
  • 44. Sun D, Yi SV. Impacts of Chromatin States and Long-Range Genomic Segments on Aging and DNA Methylation. PLoS One. 2015; 10:e0128517. https://doi.org/10.1371/journal.pone.0128517 [PubMed]
  • 45. Maugeri A, Barchitta M, Mazzone MG, Giuliano F, Agodi A. Complement System and Age-Related Macular Degeneration: Implications of Gene-Environment Interaction for Preventive and Personalized Medicine. BioMed Res Int. 2018; 2018:7532507. https://doi.org/10.1155/2018/7532507 [PubMed]
  • 46. Morgan BP. Complement in the pathogenesis of Alzheimer’s disease. Semin Immunopathol. 2018; 40:113–24. https://doi.org/10.1007/s00281-017-0662-9 [PubMed]
  • 47. Gewin LS. Renal Tubule Repair: Is Wnt/β-Catenin a Friend or Foe? Genes (Basel). 2018; 9:E58. https://doi.org/10.3390/genes9020058 [PubMed]
  • 48. Liu H, Fergusson MM, Castilho RM, Liu J, Cao L, Chen J, Malide D, Rovira II, Schimel D, Kuo CJ, Gutkind JS, Hwang PM, Finkel T. Augmented Wnt signaling in a mammalian model of accelerated aging. Science. 2007; 317:803–06. https://doi.org/10.1126/science.1143578 [PubMed]
  • 49. O’Sullivan ED, Hughes J, Ferenbach DA. Renal Aging: causes and Consequences. J Am Soc Nephrol. 2017; 28:407–20. https://doi.org/10.1681/ASN.2015121308 [PubMed]
  • 50. Maarouf OH, Aravamudhan A, Rangarajan D, Kusaba T, Zhang V, Welborn J, Gauvin D, Hou X, Kramann R, Humphreys BD. Paracrine Wnt1 Drives Interstitial Fibrosis without Inflammation by Tubulointerstitial Cross-Talk. J Am Soc Nephrol. 2016; 27:781–90. https://doi.org/10.1681/ASN.2014121188 [PubMed]
  • 51. Tan RJ, Zhou D, Zhou L, Liu Y. Wnt/β-catenin signaling and kidney fibrosis. Kidney Int Suppl (2011). 2014; 4:84–90. https://doi.org/10.1038/kisup.2014.16 [PubMed]
  • 52. Rong S, Zhao X, Jin X, Zhang Z, Chen L, Zhu Y, Yuan W. Vascular calcification in chronic kidney disease is induced by bone morphogenetic protein-2 via a mechanism involving the Wnt/beta-catenin pathway. Cell Physiol Biochem. 2014; 34:2049–60. https://doi.org/10.1159/000366400 [PubMed]
  • 53. DiRocco DP, Kobayashi A, Taketo MM, McMahon AP, Humphreys BD. Wnt4/β-catenin signaling in medullary kidney myofibroblasts. J Am Soc Nephrol. 2013; 24:1399–412. https://doi.org/10.1681/ASN.2012050512 [PubMed]
  • 54. Surendran K, McCaul SP, Simon TC. A role for Wnt-4 in renal fibrosis. Am J Physiol Renal Physiol. 2002; 282:F431–41. https://doi.org/10.1152/ajprenal.0009.2001 [PubMed]
  • 55. Xu W, Zhou W, Cheng M, Wang J, Liu Z, He S, Luo X, Huang W, Chen T, Yan W, Xiao J. Hypoxia activates Wnt/β-catenin signaling by regulating the expression of BCL9 in human hepatocellular carcinoma. Sci Rep. 2017; 7:40446. https://doi.org/10.1038/srep40446 [PubMed]
  • 56. Witkiewicz AK, Knudsen ES. Retinoblastoma tumor suppressor pathway in breast cancer: prognosis, precision medicine, and therapeutic interventions. Breast Cancer Res. 2014; 16:207. https://doi.org/10.1186/bcr3652 [PubMed]
  • 57. Wang Z, Havasi A, Gall JM, Mao H, Schwartz JH, Borkan SC. Beta-catenin promotes survival of renal epithelial cells by inhibiting Bax. J Am Soc Nephrol. 2009; 20:1919–28. https://doi.org/10.1681/ASN.2009030253 [PubMed]
  • 58. Zhou D, Tan RJ, Fu H, Liu Y. Wnt/β-catenin signaling in kidney injury and repair: a double-edged sword. Lab Invest. 2016; 96:156–67. https://doi.org/10.1038/labinvest.2015.153 [PubMed]
  • 59. He W, Kang YS, Dai C, Liu Y. Blockade of Wnt/β-catenin signaling by paricalcitol ameliorates proteinuria and kidney injury. J Am Soc Nephrol. 2011; 22:90–103. https://doi.org/10.1681/ASN.2009121236 [PubMed]
  • 60. Wang SY, Yang SX, Zhao XX, Chen F, Shi J. Expression of the Wnt/β-catenin signal pathway in patients with acute renal injury. Eur Rev Med Pharmacol Sci. 2017; 21:4661–67. [PubMed]
  • 61. DiRocco DP, Bisi J, Roberts P, Strum J, Wong KK, Sharpless N, Humphreys BD. CDK4/6 inhibition induces epithelial cell cycle arrest and ameliorates acute kidney injury. Am J Physiol Renal Physiol. 2014; 306:F379–88. https://doi.org/10.1152/ajprenal.00475.2013 [PubMed]
  • 62. Ziegler N, Awwad K, Fisslthaler B, Reis M, Devraj K, Corada M, Minardi SP, Dejana E, Plate KH, Fleming I, Liebner S. β-Catenin Is Required for Endothelial Cyp1b1 Regulation Influencing Metabolic Barrier Function. J Neurosci. 2016; 36:8921–35. https://doi.org/10.1523/JNEUROSCI.0148-16.2016 [PubMed]
  • 63. Kwon YJ, Baek HS, Ye DJ, Shin S, Kim D, Chun YJ. CYP1B1 Enhances Cell Proliferation and Metastasis through Induction of EMT and Activation of Wnt/β-Catenin Signaling via Sp1 Upregulation. PLoS One. 2016; 11:e0151598. https://doi.org/10.1371/journal.pone.0151598 [PubMed]
  • 64. Mitsui Y, Chang I, Fukuhara S, Hiraki M, Arichi N, Yasumoto H, Hirata H, Yamamura S, Shahryari V, Deng G, Wong DK, Majid S, Shiina H, et al. CYP1B1 promotes tumorigenesis via altered expression of CDC20 and DAPK1 genes in renal cell carcinoma. BMC Cancer. 2015; 15:942. https://doi.org/10.1186/s12885-015-1951-0 [PubMed]
  • 65. Park YS, Kwon YJ, Chun YJ. CYP1B1 Activates Wnt/β-Catenin Signaling through Suppression of Herc5-Mediated ISGylation for Protein Degradation on β-Catenin in HeLa Cells. Toxicol Res. 2017; 33:211–18. https://doi.org/10.5487/TR.2017.33.3.211 [PubMed]
  • 66. Reis ES, Mastellos DC, Ricklin D, Mantovani A, Lambris JD. Complement in cancer: untangling an intricate relationship. Nat Rev Immunol. 2018; 18:5–18. https://doi.org/10.1038/nri.2017.97 [PubMed]
  • 67. Corrales L, Ajona D, Rafail S, Lasarte JJ, Riezu-Boj JI, Lambris JD, Rouzaut A, Pajares MJ, Montuenga LM, Pio R. Anaphylatoxin C5a creates a favorable microenvironment for lung cancer progression. J Immunol. 2012; 189:4674–83. https://doi.org/10.4049/jimmunol.1201654 [PubMed]
  • 68. Gunn L, Ding C, Liu M, Ma Y, Qi C, Cai Y, Hu X, Aggarwal D, Zhang HG, Yan J. Opposing roles for complement component C5a in tumor progression and the tumor microenvironment. J Immunol. 2012; 189:2985–94. https://doi.org/10.4049/jimmunol.1200846 [PubMed]
  • 69. Vadrevu SK, Chintala NK, Sharma SK, Sharma P, Cleveland C, Riediger L, Manne S, Fairlie DP, Gorczyca W, Almanza O, Karbowniczek M, Markiewski MM. Complement c5a receptor facilitates cancer metastasis by altering T-cell responses in the metastatic niche. Cancer Res. 2014; 74:3454–65. https://doi.org/10.1158/0008-5472.CAN-14-0157 [PubMed]
  • 70. Sturmlechner I, Durik M, Sieben CJ, Baker DJ, van Deursen JM. Cellular senescence in renal ageing and disease. Nat Rev Nephrol. 2017; 13:77–89. https://doi.org/10.1038/nrneph.2016.183 [PubMed]
  • 71. Baker DJ, Childs BG, Durik M, Wijers ME, Sieben CJ, Zhong J, Saltness RA, Jeganathan KB, Verzosa GC, Pezeshki A, Khazaie K, Miller JD, van Deursen JM. Naturally occurring p16(Ink4a)-positive cells shorten healthy lifespan. Nature. 2016; 530:184–89. https://doi.org/10.1038/nature16932 [PubMed]
  • 72. Eren M, Boe AE, Murphy SB, Place AT, Nagpal V, Morales-Nebreda L, Urich D, Quaggin SE, Budinger GR, Mutlu GM, Miyata T, Vaughan DE. PAI-1-regulated extracellular proteolysis governs senescence and survival in Klotho mice. Proc Natl Acad Sci USA. 2014; 111:7090–95. https://doi.org/10.1073/pnas.1321942111 [PubMed]
  • 73. Wang WJ, Cai GY, Chen XM. Cellular senescence, senescence-associated secretory phenotype, and chronic kidney disease. Oncotarget. 2017; 8:64520–33. https://doi.org/10.18632/oncotarget.17327 [PubMed]
  • 74. Peng Q, Li K, Smyth LA, Xing G, Wang N, Meader L, Lu B, Sacks SH, Zhou W. C3a and C5a promote renal ischemia-reperfusion injury. J Am Soc Nephrol. 2012; 23:1474–85. https://doi.org/10.1681/ASN.2011111072 [PubMed]
  • 75. Martin IV, Bohner A, Boor P, Shagdarsuren E, Raffetseder U, Lammert F, Floege J, Ostendorf T, Weber SN. Complement C5a receptors C5L2 and C5aR in renal fibrosis. Am J Physiol Renal Physiol. 2018; 314:F35–46. https://doi.org/10.1152/ajprenal.00060.2017 [PubMed]
  • 76. van Werkhoven MB, Damman J, Daha MR, Krikke C, van Goor H, van Son WJ, Hillebrands JL, van Dijk MC, Seelen MA. Novel insights in localization and expression levels of C5aR and C5L2 under native and post-transplant conditions in the kidney. Mol Immunol. 2013; 53:237–45. https://doi.org/10.1016/j.molimm.2012.08.013 [PubMed]
  • 77. Mencke R, Olauson H, Hillebrands JL. Effects of Klotho on fibrosis and cancer: A renal focus on mechanisms and therapeutic strategies. Adv Drug Deliv Rev. 2017; 121:85–100. https://doi.org/10.1016/j.addr.2017.07.009 [PubMed]
  • 78. Zhou L, Li Y, Zhou D, Tan RJ, Liu Y. Loss of Klotho contributes to kidney injury by derepression of Wnt/β-catenin signaling. J Am Soc Nephrol. 2013; 24:771–85. https://doi.org/10.1681/ASN.2012080865 [PubMed]
  • 79. Montgomery RA, Orandi BJ, Racusen L, Jackson AM, Garonzik-Wang JM, Shah T, Woodle ES, Sommerer C, Fitts D, Rockich K, Zhang P, Uknis ME. Plasma-Derived C1 Esterase Inhibitor for Acute Antibody-Mediated Rejection Following Kidney Transplantation: Results of a Randomized Double-Blind Placebo-Controlled Pilot Study. Am J Transplant. 2016; 16:3468–78. https://doi.org/10.1111/ajt.13871 [PubMed]
  • 80. Luo C, Zhou S, Zhou Z, Liu Y, Yang L, Liu J, Zhang Y, Li H, Liu Y, Hou FF, Zhou L. Wnt9a Promotes Renal Fibrosis by Accelerating Cellular Senescence in Tubular Epithelial Cells. J Am Soc Nephrol. 2018; 29:1238–56. https://doi.org/10.1681/ASN.2017050574 [PubMed]
  • 81. Chkhotua A, Shohat M, Tobar A, Magal N, Kaganovski E, Shapira Z, Yussim A. Replicative senescence in organ transplantation-mechanisms and significance. Transpl Immunol. 2002; 9:165–71. https://doi.org/10.1016/S0966-3274(02)00003-5 [PubMed]
  • 82. Braun H, Schmidt BM, Raiss M, Baisantry A, Mircea-Constantin D, Wang S, Gross ML, Serrano M, Schmitt R, Melk A. Cellular senescence limits regenerative capacity and allograft survival. J Am Soc Nephrol. 2012; 23:1467–73. https://doi.org/10.1681/ASN.2011100967 [PubMed]
  • 83. Delpech PO, Thuillier R, SaintYves T, Danion J, Le Pape S, van Amersfoort ES, Oortwijn B, Blancho G, Hauet T. Inhibition of complement improves graft outcome in a pig model of kidney autotransplantation. J Transl Med. 2016; 14:277. https://doi.org/10.1186/s12967-016-1013-7 [PubMed]
  • 84. Curci C, Castellano G, Stasi A, Divella C, Loverre A, Gigante M, Simone S, Cariello M, Montinaro V, Lucarelli G, Ditonno P, Battaglia M, Crovace A, et al. Endothelial-to-mesenchymal transition and renal fibrosis in ischaemia/reperfusion injury are mediated by complement anaphylatoxins and Akt pathway. Nephrol Dial Transplant. 2014; 29:799–808. https://doi.org/10.1093/ndt/gft516 [PubMed]
  • 85. Golriz M, Fonouni H, Nickkholgh A, Hafezi M, Garoussi C, Mehrabi A. Pig kidney transplantation: an up-to-date guideline. Eur Surg Res. 2012; 49:121–9. https://doi.org/10.1159/000343132 [PubMed]