Tacrolimus increases the expression level of the chemokine receptor CXCR2 to promote renal fibrosis progression.

Tacrolimus is one of the most used and effective immunosuppressive agents currently available in the clinic; however, its use is limited by nephrotoxicity, which is the main secondary effect of this drug. The mechanisms underlying tacrolimus‑induced nephrotoxicity remain unknown. The present study aimed to investigate the mechanism underlying tacrolimus‑induced nephrotoxicity and to identify novel potential targets. Masson staining, Sirius red staining and periodic acid‑silver methenamine staining were used to observe kidney pathological changes. Immunohistochemical and immunofluorescent analyses were performed to examine the expression levels of vimentin, E‑cadherin and α‑smooth muscle actin (α‑SMA). Transcriptomics and bioinformatics analyses were performed to investigate the nephrotoxicity mechanism induced by tacrolimus using RNA‑sequencing, differentially expressed genes identification and annotation, and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. The present results demonstrated that compared with the normal control group, the tacrolimus nephrotoxicity group exhibited severe renal fibrosis (P<0.05), upregulated vimentin (P<0.01), downregulated E‑cadherin (P<0.05) and upregulated α‑SMA (P<0.01). Transcriptomics and bioinformatics analyses identified the pathway 'cytokine‑cytokine receptor interaction' as the most significantly enriched (P<0.05). Moreover, KEGG pathway enrichment analysis identified that tacrolimus increased the expression levels of chemokine (C‑X‑C) motif ligand (CXCL)1, CXCL2 and CXCL3 and the chemokine receptor C‑X‑C chemokine receptor type 2 (CXCR2). Collectively, the present study suggested that tacrolimus increases the level of chemokine receptor CXCR2 to promote renal fibrosis progression, which is one of the potential mechanisms underlying tacrolimus‑induced nephrotoxicity.

Introduction

Tacrolimus, a potent immunosuppressive agent, was developed in the 1990s and can be isolated from the bacteria Streptomyces tsukubaensis. FK-506-binding protein 12 (FKBP-12) is an important immunophilin targeted by tacrolimus in T cells, and tacrolimus can form a complex with FKBP-12, thus repressing the phosphatase calcineurin, an enzyme necessary to activate the nuclear factor of T cells (NF-AT) (1-4). NF-AT serves a key role in the transcription of cytokine-encoding genes in T cells (5).

Tacrolimus has primarily been used for treating patients who have received allogeneic organ transplants or patients with autoimmune diseases (4,6-11). However, clinical management of tacrolimus therapy can be challenging due to its narrow treatment range and significant variabilities within and among individuals (12,13), and these limitations can be caused by multiple factors that interfere with its metabolism. Therefore, in clinical settings, therapeutic drug monitoring is used to optimize the treatment regimen (14).

Tacrolimus has various side effects, and nephrotoxicity is the most common, occurring in ~50% of patients treated with tacrolimus (15). Renal fibrosis is commonly associated with nephrotoxicity (16-20). Fibrosis plays a crucial role in cadmium-induced nephrotoxicity (16), cyclosporine nephrotoxicity (17,19), nephrotoxicity induced by oral sodium nitrite (18) and aristolochic acid nephrotoxicity (20). In addition, renal fibrosis is an important process underlying tacrolimus nephrotoxicity (21,22). Therefore, it is important to identify the factors that lead to renal fibrosis following tacrolimus nephrotoxicity and control its development (21,22). The present study aimed to investigate the mechanism underlying renal fibrosis induced by tacrolimus and to identify novel potential targets.

Materials and methods

Animal experiments

In total, 16 specific pathogen free male Wistar rats (age, 3 weeks; weight, 64±3 g) were purchased from Shanghai SLAC Laboratory Animal Co., Ltd. The animals were housed in standard cages and maintained under standard conditions at a constant room temperature of 20-25°C, a humidity of 40-70% and a 12/12 h light/dark cycle. All rats had free access to regular chow and water. The rats were randomly divided into a normal control (NC) group (n=8) and a tacrolimus nephrotoxicity (NE) group (n=8). The NE group was intraperito neally injected with tacrolimus (Astellas Ireland Co., Ltd.) at 2 mg/kg per day (23,24). The NC group received daily intraperito neal administrations of equal volumes of tacrolimus solvent, which consisted of polyoxyethylene hydrogenated castor oil and absolute ethyl alcohol, for 2 weeks. Animal protocols and procedures were approved by The Animal Care and Use Committee of Children's Hospital of Fudan University (Shanghai, China) and complied with the appropriate institutional regulations.

Sample collection

After 2 weeks of intervention, Wistar rats were anesthetized with 10% chloral hydrate (300 mg/kg) by intraperitoneal injection. Animals did not present obvious signs of peritonitis in the present study. Blood samples were collected from the abdominal aorta. The rats were euthanized by cervical dislocation under deep anesthesia with 10% chloral hydrate. The death of rats was verified by the heartbeat, breathing and neural reflex. After Wistar rats were sacrificed, the kidneys were removed and weighed. The kidney index was calculated using the following formula: Kidney index = (average kidney weight/body weight) (25). Paraformaldehyde (4%) was used to fix the fresh kidneys, and liquid nitrogen was used to rapidly freeze the residual kidney tissues, which were stored at -80°C until further analysis.

Observation of kidney histological

Kidney tissues were fixed in 4% paraformaldehyde at 4°C for 24 h, embedded in paraffin and cut into 4-µm-thick sections (26). The sections were then stained by Masson staining, Sirius red staining and periodic acid-silver methenamine (PASM) staining (all from Wuhan Servicebio Technology Co., Ltd.). Masson staining was performed at room temperature (25°C), with iron hematoxylin for 3 min, Ponceau S solution for 5 min, phosphomolybdic acid for 1 min and aniline blue for 3 min. Masson staining was used to dye the collagen fibers blue. Sirius red staining was performed at room temperature (25°C) for 1 h. Following Sirius red staining, the collagen fibers were stained red. For PASM staining, the sections were stained with periodic acid for 15 min at room temperature (25°C) and hexamine silver working solution for 40 min at 58°C. Following PASM staining, elastic and mesh fibers were stained in black. The aforementioned staining was observed using a light microscope (Nikon Eclipse E100) at a magnification of ×400. The sum of integral optical density (IOD) was used to quantify the results and calculated using Image-Pro Plus 6.0 software (Media Cybernetics, Inc.).

Immunohistochemistry assays

The 4-µm-thick sections were baked at 60°C for 2 h and then deparaffinized in xylene, rehydrated in 100, 85 and 75% alcohol, and then washed in water. Subsequently, antigen retrieval was performed. The samples were treated with 3% hydrogen peroxide to block endogenous peroxidase activity and then blocked with BSA for 30 min at room temperature (25°C). Sections were incubated with vimentin antibody (1:500; Wuhan Servicebio Technology Co., Ltd.; cat. no. GB11192) at 4°C overnight. Subsequently, the appropriate secondary antibodies (1:200; HRP-goat anti-rabbit IgG; Wuhan Servicebio Technology Co., Ltd.; cat. no. GB23303) were incubated with the sections in the dark at room temperature (25°C) for 50 min. Subsequently, 3,3'-diami-nobenzidine dye solution (Dako; Agilent Technologies, Inc.) was incubated with the sections at room temperature in the dark. Sections were observed and images were acquired using a light microscope (magnification, ×400). The sum of IOD was used to quantify the results and calculated using Image-Pro Plus 6.0 software (Media Cybernetics, Inc.). The color of vimentin staining was claybank.

Immunofluorescence method

The 4-µm-thick sections were baked at 60°C for 2 h and then deparaffinized in xylene, rehydrated in 100, 85 and 75% alcohol, and then washed in water. Subsequently, antigen retrieval was performed. Then, to reduce the spontaneous auto-fluorescence, tissue auto-fluorescence quencher (Wuhan Servicebio Technology Co., Ltd.; cat. no. G1221) was added for 5 min and samples were rinsed in water for 10 min. Sections were blocked using BSA for 30 min at room temperature (25°C). Samples were incubated with anti-E-cadherin (1:5,000; Wuhan Servicebio Technology Co., Ltd.; cat. no. GB12082) and anti-α-smooth muscle actin (α-SMA) antibodies (1:10,000; Wuhan Servicebio Technology Co., Ltd.; cat. no. GB13044) at 4°C overnight. Following incubation with the primary antibodies, the sections were washed and incubated with the corresponding secondary antibodies (1:500; HRP-goat anti-mouse IgG; Wuhan Servicebio Technology Co., Ltd.; cat. no. GB23301) in the dark at room temperature (25°C) for 50 min. DAPI solution (Wuhan Servicebio Technology Co., Ltd.) was then added and the samples were incubated in the dark at room temperature for 10 min. An anti-fluorescence quenching solution was used to seal the samples. Samples were observed and images were acquired using a confocal microscope (Nikon Corporation; magnification, ×400). The sum of IOD was used to quantify the results and calculated using Image-Pro Plus 6.0 software (Media Cybernetics, Inc.). E-cadherin was stained in red and α-SMA was stained in green.

RNA library construction and sequencing

Suzhou Basepair Biotechnology Company (http://www.basepair.cn/) performed the RNA library construction and sequencing in biological triplicates using an Illumina HiSeq X Ten Sequencing system (Illumina, Inc.). The main sequencing processes included sequencing data quality preprocessing, reference genome alignment, gene expression analysis and differential expression analysis (25,27).

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment

Sequencing data were analyzed as raw reads, and were saved in a FASTQ format document. To obtain clean reads, adaptor-contaminated and low-quality sequences were removed using filtering methods, as previously described (28). FastQC (version 0.11.4; http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to examine the quality of the clean reads (25). Subsequently, the reads were aligned to the reference genome Ensemble-Rnor6.0 using Hisat2 (version 2.1.0). Gene coverage was analyzed using the percentage of genes covered by the reads. Functional annotation was performed using ANNOVAR (29). Differentially expressed genes (DEGs) were identified using the DEGseq panormalage method (30), according to negative binomial distributions. The gene expression levels were analyzed according to the number of fragments per kilobase of transcript per million reads and counts values. Genes with an adjusted |log2(fold change)| >1 and P<0.05 were considered as DEGs. KEGG pathway analysis was performed for the DEGs using the KEGG Orthology-Based Annotation System (31). The significantly enriched KEGG pathways exhibited a P<0.05 (25,27,28,32).

Statistical analysis

Data are presented as the mean ± standard error of the mean. The statistical analysis was performed using SPSS (version 13.0; SPSS, Inc. and GraphPad Prism (version 5.0; GraphPad Software Inc.). Each experiment was performed three times. An unpaired Student's t-test was used to compare the protein expression level between the two groups. P<0.05 was considered to indicate a statistically significant difference.

Results

Pathological changes in kidneys induced by tacrolimus

Kidneys were examined using pathological and immunofluorescent staining (Figs. 1A and 2A, respectively). The NE group exhibited a significantly higher of kidney index compared with the NC group (Fig. 1B; P<0.05). Compared with the NC group, the NE group exhibited severe renal fibrosis. This was indicated by a significantly higher sum of IOD in the NE group following Masson staining (Fig. 1C; P<0.05), Sirius red staining (Fig. 1D; P<0.01) and PASM staining (Fig. 1E; P<0.01). In addition, the NE group exhibited significantly upregulated vimentin (Fig. 1F; P<0.01), significantly down-regulated E-cadherin (Fig. 2B; P<0.05) and significantly upregulated α-SMA (Fig. 2C; P<0.01).

Figure 1

Renal histology analysis. (A) Pathological staining. Scale bar, 100 µm. Magnification, ×400. (B) Kidney index, kidney index = average kidney weight/body weight. (C) Sum of IOD of Masson staining. (D) Sum of IOD of Sirius red staining. (E) Sum of IOD of PASM staining. (F) Sum of IOD of Vimentin expression. *P<0.05, **P<0.01 vs. NC. NC, normal control group; NE, nephrotoxicity group; IOD, integral optical density; PASM, periodic acid-silver metheramine.

Figure 2

Expression of E-cadherin and α-SMA. (A) Immunofluorescence staining. Scale bar, 100 µm, Magnification, ×400. (B) Sum of IOD of E-cadherin expression. (C) Sum of IOD of α-SMA expression. *P<0.05, **P<0.01 vs. NC. NC, normal control group; NE, nephrotoxicity group; α-SMA, α-smooth muscle actin; IOD, integral optical density.

Transcriptomic and bioinformatics analyses reveal the nephrotoxicity mechanism underlying tacrolimus

A total of six next-generation sequencing libraries, including three from the NC group and three from the NE group, were analyzed. As presented in Fig. 3, the transcriptomic analysis identified 453 DEGs, including 173 upregulated and 280 downregulated genes. The pathway demonstrating the highest enrichment following KEGG analysis was 'cytokine-cytokine receptor interaction' (Fig. 4A; P<0.05). Further analysis of this signaling pathway using KEGG identified that tacrolimus increased the expression levels of CXCL1, CXCL2, CXCL3 and the chemokine receptor CXCR2 (Fig. 4B; P<0.05). The potential mechanism underlying tacrolimus-induced nephrotoxicity is presented in Fig. 5.

Figure 3

Differentially expressed genes of the NE and NC groups identified by |log2(fold change)| >1 and P<0.05. NC, normal control group. NE, nephrotoxicity group.

Figure 4

Differential gene expression screening. (A) KEGG pathway enrichment of differentially expressed genes. (B) The pathway map of cytokine-cytokine receptor interaction. Yellow indicates upregulated genes, red indicates downregulated genes and green indicates that the gene is both upregulated and down-regulated. KEGG, Kyoto Encyclopedia of Genes and Genomes.

Figure 5

Potential mechanism of nephrotoxicity induced by tacrolimus. Tacrolimus increases the expression levels of CXCL1, CXCL2 and CXCL3 and the chemokine receptor CXCR2, which further promotes renal fibrosis progression. CXCL, chemokine (C-X-C) motif ligand; CXCR2, C-X-C chemokine receptor type 2; α-SMA, α-smooth muscle actin.

Discussion

Tacrolimus is one of the most used and effective clinical immunosuppressive agents currently available in the clinic, and it has been widely used for treating patients receiving renal (33-38), liver (39-44) and lung transplants (45), and for patients with idiopathic membranous nephropathy (46), nephritic syndrome (4) and systemic-onset juvenile idiopathic arthritis (7). However, although tacrolimus has been demonstrated to exhibit evident benefits, immunosuppressive agents are associated with the occurrence of acute or chronic renal toxicity, limiting their clinical use (15).

In addition, it has previously been reported that the epithetlial-mesenchymal transition (EMT) is significantly correlated with renal fibrosis, which is associated with calcineurin inhibitor-mediated nephrotoxicity (21). EMT is one of the basic mechanisms of renal fibrosis and involves various processes in which epithelial cells stop exhibiting epithelial characteristics, including the expression of E-cadherin, and obtain traits specific of mesenchymal cells, including the upregulation of α-SMA (25,47-48). In addition, it has previously been reported that vimentin is a potential novel biomarker in renal fibrosis (25,49).

In the present study, Masson staining, Sirius red staining and PASM staining were used to examine the pathological alterations occurring in kidneys. Compared with the NC group, the tacrolimus nephrotoxicity group exhibited severe renal fibrosis. Further analysis confirmed that vimentin was upregulated, E-cadherin was downregulated and α-SMA was upregulated in the tacrolimus-induced nephrotoxicity group. Subsequently, the exact mechanism underlying tacrolimus-induced nephrotoxicity was examined.

Transcriptome analysis can identify both coding and non-coding RNA, quantifying gene expression heterogeneity in cells, tissues, organs and even in the whole organism (50). Transcriptome analysis is important to investigate various processes (51) and it has been widely used to identify key factors in the progression of multiple diseases (27,52-53). Kim et al (52), using transcriptome analysis, identified Tensin 4 as a key effector of cetuximab and a regulator of the oncogenic activity of KRAS-mutant colorectal cancer cell lines. Yang et al (27) reported the transcriptome profiling of brain myeloid cells, and identified an upregulation of integrin subunit α L, triggering receptor expressed on myeloid cells 1 and secreted phosphoprotein 1 in Western diet-induced obesity. Siena et al (53) performed a whole transcriptome analysis in melanoma and identified a correlation between the expression level of the long non-coding RNA ZEB1-AS1 with invasive ability of melanoma cells. These previous studies suggested that the transcriptome has become a reliable tool to identify key factors in the development and progression of various diseases (27,52,53). The present study aimed to investigate the mechanism underlying tacrolimus nephrotoxicity and to identify novel potential targets via transcriptomic and bioinformatics analyses.

The KEGG enrichment analysis identified 'cytokine-cytokine receptor interaction' as the pathway most significantly enriched following tacrolimus-mediated nephrotoxicity induction. By analyzing components of the 'cytokine-cytokine receptor interaction' signaling pathway, tacrolimus was identified to increase the expression levels of CXCL1, CXCL2, CXCL3 and the chemokine receptor CXCR2.

CXCR2 is a seven-transmembrane G-protein-coupled receptor that medi ates chemotaxis during immune response, and is expressed in renal parenchymal cells and neutrophils (54-55). Dornelles et al (56) firstly reported the association between the increase in CXCR2 expression and nephrotoxicity following cyclophosphamide treatment. In addition, upregulation of CXCR2 has been reported in inflammatory diseases, including psoriasis, atherosclerosis and rheumatoid arthritis (57,58). CXCR2-knockout mice were identified to be protected against dextran sodium sulfate-mediated colitis and acute kidney injury. In addition, the expression of cytokines and chemokines and the level of neutrophil infiltration were reduced in the colon and kidney of CXCR2-knockout mice (59). Collectively, these previous studies suggest that CXCR2 may be a promoter of kidney damage.

In conclusion, the mechanism underlying tacrolimus-induced nephrotoxicity may involve the increase of the chemokine receptor CXCR2 to promote the upregulation of vimentin and α-SMA, and the downregulation of E-cadherin, thus accelerating the renal fibrosis progression. However, the present analysis was performed in animal models, and validation of the present results is required in the future by analyzing blood and kidney biopsies from patients with nephrotoxicity caused by tacrolimus.