Impaired ATG16L-Dependent Autophagy Promotes Renal Interstitial Fibrosis in Chronic Renal Graft Dysfunction Through Inducing EndMT By NF-kB Signal Pathway

Zeping Gui ^1,2†, Chuanjian Suo ^1†, Zijie Wang ^1†, Ming Zheng ¹, Shuang Fei ¹, Hao Chen ¹, Li Sun ¹, Zhijian Han ¹, Jun Tao ¹, Xiaobin Ju ¹, Haiwei Yang ¹, Min Gu ^2* and Ruoyun Tan ^1*

Chronic renal graft dysfunction (CAD) is caused by multiple factors, including glomerular sclerosis, inflammation, interstitial fibrosis, and tubular atrophy (IF/TA). However, the most prominent elements of CAD are IF/TA. Our studies have confirmed that endothelial- mesenchymal transition (EndMT) is an important source to allograft IF/TA. The characteristic of EndMT is the loss of endothelial marker and the acquisition of mesenchymal or fibroblastic phenotypes. Autophagy is an intracellular degradation pathway that is regulated by autophagy-related proteins and plays a vital role in many fibrotic conditions. However, whether or not autophagy contributes to fibrosis of renal allograft and how such mechanism occurs still remains unclear. Autophagy related 16 like gene (ATG16L) is a critical autophagy-related gene (ARG) necessary for autophagosome formation. Here, we first analyzed kidney transplant patient tissues from Gene Expression Omnibus (GEO) datasets and 60 transplant patients from our center. Recipients with stable kidney function were defined as a non-CAD group and all patients in the CAD group were histopathologically diagnosed with CAD. Results showed that ATG16L, as one significant differential ARG, was less expressed in the CAD group compared to the non-CAD group. Furthermore, we found there were fewer autophagosomes and autolysosomes in transplanted kidneys of CAD patients, and downregulation of autophagy is a poor prognostic factor. In vitro, we found out that the knockdown of ATG16L enhanced the process of EndMT in human renal glomerular endothelial cells (HRGECs). In vivo, the changes of EndMT and autophagic flux were then detected in rat renal transplant models of CAD. We demonstrated the occurrence of EndMT and indicated that the abundance of ATG16L was accompanied by the dynamic autophagic flux change along different stages of kidney transplantation. Mechanistically, knockdown of ATG16L, specifically in endothelial cells, reduced NF-kB degradation and excreted inflammatory cytokines (IL-1b, IL-6, and TNF-a), which could facilitate EndMT. In conclusion, ATG16L-dependent autophagic flux caused by transplant showed progressive loss increase over time.

Inflammatory cytokines from this process promoted EndMT, thereby leading to the progression of CAD. ATG16L served as a negative regulator of EndMT and development of renal graft fibrosis, and autophagy can be explored as a potential therapeutic target for chronic renal graft dysfunction.

Keywords: chronic renal graft dysfunction, renal interstitial fibrosis, ATG16L, autophagy, EndMT, inflammatory cytokines

For more information:ali.ma@wecistanche.com

INTRODUCTION

Kidney transplantation is one of the optimal treatments for patients with uremia because it significantly improves their quality of life (1). However, there is a relatively high number of renal allograft dysfunction after transplantation due to the chronic progressive deterioration of renal function, which is a key factor affecting the long-term survival of transplanted kidneys (2). Chronic renal graft dysfunction (CAD), formerly known as chronic allograft nephropathy, is a multifactorial condition associated with progressive renal interstitial fibrosis. Various factors are known to contribute to the loss of renal allograft function, including, but not limited to, acute and chronic rejection, ischemia and reperfusion inflammatory and tissue rebuilding process, and drug-related nephrotoxicity (3).

CAD is morphologically characterized by inflammation, progressive interstitial fibrosis and tubular atrophy(IF/TA), and glomerular sclerosis (4). Among them, IF/TA is a key factor determining renal allograft function, but the mechanisms by which it occurs are still unknown. Several factors have been identified to contribute to the high proportion of renal allograft loss. In our previous study, we also confirmed that the main pathological process of CAD was renal allograft IF/TA, which is characterized by excessive deposition of extracellular matrix in transplanted renal tubular and interstitial tissue (5). Studies have verified that collagens are mainly secreted by myofibroblasts, and collagens secretion lead to extracellular matrix sedimentation and subsequently transplant kidney IF/TA. There are four principal cells involved in the formation of myofibroblasts: epithelial cells, endothelial cells, bone marrow-derived fibroblasts, and microvascular pericytes (6, 7). They play a crucial role in repairing and protecting the integrity of kidney tissue. Among them, external stimuli enable the differentiation of intrinsic kidney cells such as epithelial cells to myofibroblasts, which can produce an extracellular matrix. This process is called epithelial-to-mesenchymal transition (EMT). EMT usually involved three types: (I) embryogenesis, (II) tissue repair and fibrosis, (III) metastasis (8). Several studies have been done to elucidate the molecular and cellular mechanisms of type II EMT in organ fibrosis (9, 10), however, roles of EMT in different primary fibroblast-generating processes during the CAD progression still remain inconclusive.

With kidney transplantation, there is always an increased likelihood of damage to the allograft macro- and microvasculature, due to the physical location of the allograft endothelium that makes it an initial target of choice for allograft injury (11). Allografts act as potential targets for immune response mediated by quite a number of serum inflammatory cytokines. Studies indicated that many inflammatory cytokines involved in the CAD process (12) and endothelial cells at sites of renal allografts were not only participants in inflammation but also regulators of inflammation (11). Our previous research revealed that TNF‐a is more highly expressed in the CAD group than the non-CAD group. TNF‐a also facilitated the EMT process in human proximal tubular cells (HK2) (13). Given the intimate contact between endothelial cells and blood, circulating inflammatory cells and cytokines will first attack the allograft endothelium. Chronic stimulation from inflammatory cytokines leads to endothelial cells undergoing disorders of cell structure and internal environment, which results in endothelial injury. Emerging studies suggest that endothelial injury contributes to extracellular matrix deposition and plays a key role in organ fibrosis diseases. Recently, Endothelial-to-mesenchymal transition (EndMT) is considered as the principal cause of endothelial injury and promotes the progress of IF/TA during renal fibrosis disease (5, 9). EndMT is a distinctive type of EMT. It is characterized by cells that gradually lose endothelial markers, such as CD31 and CD34, and gain mesenchymal or myofibroblastic phenotype, such as a-smooth muscular actin (a-SMA), collagen I, and fibronectin (FN). The regulatory mechanism for EndMT remains a complex issue. The occurrence of EndMT could be affected by many factors such as oxidative stress, hypoxic, and various injuries (14, 15). But almost all of these regulatory factors eventually come together in one direct action, which is inflammatory cytokines. In our previous studies, we showed that the EndMT was an important factor in the pathogenesis of IF/TA and CAD through the TGF-b/Smad signaling pathways (5). In addition, still other scholars suggested that TGF-b expression was up-regulated by TNF‐a (16). So, inflammatory cytokines such as TNF‐a may also be an important pathogenic factor in vascular endothelial injury that characterizes CAD progression. Although the initial factors of EndMT were different in various microenvironments, the final outcome was the same. Hence, it sufficed to hypothesize that there must be a central regulatory mechanism in the production of inflammatory cytokines and progression of EndMT.

Autophagy is a cellular pathway responsible for protein and organelle degradation (17). The process of autophagy is also known as autophagic flux, and the strength of this process often represents the degree of activation of autophagy. Autophagy is involved in different renal pathophysiological processes, including glomerulosclerosis, diabetic nephropathy, and cystic kidney disease (18). Some studies reported that autophagy was a cytoprotective process. A classic defense for this presupposition is as seen in mice with proximal tubule ATG5 deletion where the deletion promoted more severe renal fibrosis due to impaired autophagy (19). Other studies have also found that proximal tubule epithelial cells conditional knockdown ATG5 aggravated acute kidney injury in the early stage, but decelerated the progression of kidney fibrosis in the recovery or repair process (20). However, the specific interaction between autophagy and EMT or EndMT on fibrosis disease still remains controversial. Several studies suggested that autophagy inhibition could induce EMT and fibrosis by affecting the aberrant epithelial–fibroblast crosstalk in idiopathic pulmonary fibrosis (21). Others showed that Rapamycin promotes EndMT through the activation of autophagy during premature senescence (22). ATG16L gene mainly includes ATG16L1 and ATG16L2. ATG16L2 homo- and hetero-oligomerizes with ATG16L1. But ATG16L2 has less associated with autophagy (23). ATG16L1 has been identified as an important autophagy-related gene, promotes autophagosome formation at the plasma membrane, and plays a critical role in the lipidated form of LC3 (23). Here, ATG16L1 is collectively referred to as ATG16L. As a key factor in autophagy, ATG16L was once reported that be related to the pathogenesis of several inflammatory diseases. For instance, the mutation of ATG16L conferred a strong predisposition to Crohn’s disease development (24). However, the mechanism of ATG16L action and autophagy has not yet been well-studied in a targeted approach such as in the management and prognostication of kidney transplantation. So the dynamic role of autophagic flux and ATG16L-dependent autophagy must be the focus of research for potential therapeutics in CAD progression.

In this study, we explored the potential role of ATG16L- dependent autophagy in CAD progression. We revealed that the loss of ATG16L and autophagy were associated with the occurrence of EndMT in clinical CAD patients and kidney transplanted rats. We then found that the levels of EndMT and renal allograft interstitial fibrosis were enhanced after altering autophagy activity by knocking down the ATG16L gene. In addition, potential underlying mechanisms were also explored by RNA sequencing, we found that impairment of ATG16L contributed to the nuclear factor-kB (NF-kB) pathway activation and the inflammatory cytokines secretion, then promoted EndMT progression.

MATERIALS AND METHODS

Ethics Statement

The study’s protocols complied with the Declaration of Helsinki and Istanbul. Human studies were examined and authorized by the local ethics committee of the First Affiliated Hospital of Nanjing Medical University (ID: IACUC-2010020). Informed consent was acquired from all transplant recipients as well as nephrectomy patients who were included in this research.

The studies involving rats were examined and authorized by the local ethics committee of the First Affiliated Hospital of Nanjing Medical University (ID: IACUC-2010020). Written informed consent is given by the owners for the participation of their animals in this research.

Sample Collection

Sixty of adults who had received living or deceased donor kidney transplants started to be followed up from January 2010 to December 2017 at First Affiliated Hospital of Nanjing Medical University. Patients with serum creatinine level consistently < 141.46 mmol/L (1.6 mg/dl) for at least 12 months after kidney transplantation and no other complications such as episodes of significant rejection, drug toxicity injury, and infection were assigned to the non-CAD group. Patients in the CAD group were defined as elevated serum creatinine greater than or equal to 141.46 mmol/L (1.6 mg/dl) for at least three months, and they were diagnosed by two independent pathology experts combined with biopsy results, laboratory indexes, and imaging features. Adjacent normal kidney tissues were obtained from regions outside the tumor margin (>5cm) in patients with radical nephrectomy operations. Blood samples were obtained after patients had given informed consent. The specific step can be referred to in our previous studies (5). The baseline characteristics of patients in the CAD group and the non-CAD group were as shown in Table 1.

Rats and Animals Models

Adult male F344 and Lewis rats (Weight 250 ± 10.3 g) were obtained from Charles River Laboratories (Beijing, China) and abided by the guidelines of the Institutional Animal Care and Use Committee at Nanjing Medical University. The animals were handled in accordance with the norms of Nanjing Medical University and the guidelines published by the US National Institutes of Health.

All these animals experienced orthotopic left kidney transplantation. Lewis rats were used as recipients and syngeneic donors (Syn group), F344 rats as allogeneic donors (Allo group). The right kidney was excised simultaneously. The average time of cold ischemia was less than 20 minutes and warm ischemia was less than 35 minutes. Cyclosporine A (5 mg/kg, QD, IP; Neoral, Novartis, Switzerland) was used for 14 days to avoid acute rejection,

Pharmaceutical Treatment and Tissue Harvest

The transplant kidneys from rats were harvested at weeks 4, 8, 12, and 16 after surgery. Paraffin-embedded formalin-fixed (10% neutral formalin) renal allograft tissue specimens were obtained for histological and IHC/IF staining. Remaining kidney stored at − 80°C refrigerators for detection of RNA and protein.

HE and Masson Trichrome Staining Assay

Protocols of HE and Masson trichrome staining assay can be referred to in our previous studies (5). To identify the severity of CAD and the extent of the fibrotic area, the renal allograft fibrosis area was quantified with Masson trichrome staining. For each slice, five random visual fields under × 400 microscopes were selected. Area positive for Masson trichrome was measured by two pathologists blinded to the experimental design using the Image‐Pro Plus (Media Cybernetics, Rockville, MD).

Immunohistochemistry Staining Assay

Kidney tissues fixed with formalin were cut into 3 mm thick paraffin sections. The process of deparaffinized, hydrated, and antigen-retrieved was consistent with our previous studies (5). The antibodies [anti‐ATG16L (1:100; Abcam, USA), anti‐a‐ SMA (1:200; Abcam, USA), anti‐Fibronectin (1:100; Abcam, USA) and anti‐CD31 (1:100; CST, USA)] were incubated overnight at 4°C after sections blocked with 10% normal donkey serum. Next steps were also performed as described earlier with biotinylated goat anti‐mouse/rabbit IgG (0.5 mg/mL;

Abcam) and substrate 3-amino-9-ethyl carbazole or 3,3′

diaminobenzidine (Vector Laboratories, Burlingame, CA). The stained slides were photographed using a Nikon Eclipse 80i microscope equipped with a digital camera (DS-Ri1, Nikon, Shanghai, China).

Indirect Immunofluorescence Staining Assay

After fixing with 4% formaldehyde solution, renal allograft sections and cell climbing slices were penetrated with 0.1% Triton X-100 for 1 h, then blocked with 5% goat serum for 1 h. Afterward, sections and cell climbing slices were incubated with the anti‐ NF-kB p65 (1:200; Abcam, USA) at 4°C overnight. Incubation conditions of secondary antibody Dapi can refer to our previous studies (5). Slides were viewed with a Nikon Eclipse 80i fluorescence microscope (DS-Ri1, Nikon, Shanghai, China).

Autophagic Flux Detection

AAV-mRFP-GFP-LC3 (Hanbio, Shanghai, China) was stereotactically injected into the left kidney of Lewis rats (3 mL) 14 days before transplantation. The kidney was removed at 4, 8, 12, 16 weeks after transplantation and fixed with 4% paraformaldehyde for 24 h. Slides (30 mm) were viewed with a Nikon Eclipse 80i fluorescence microscope (DS-Ri1, Nikon, Shanghai, China). Yellow spots indicate autophagosomes and red spots indicate autolysosomes. When the red signal is stronger than the yellow signal represent autophagic flux is activated. When more yellow signals than red signals represent autophagic flux is impaired.

Electron Microscopy

The samples were fixed with ice-cold glutaraldehyde (3% in 0.1 M cacodylate buffer, pH 7.4) and further processed by the Core Facility (Service, Wuhan, China). The observation was performed on a JEOL JEM-2100 transmission electron microscope.

Renal Function Detection

We used a rat QuantiChrom Creatinine Assay Kit and QuantiChrom Urea Assay Kit (Jiancheng, Beijing, China) to detect the concentrations of rat blood creatinine and urea nitrogen according to the manufacturer’s instructions.

Uria Protein Detection

The 24 h urine samples were collected using metabolic cages, and urinary protein excretion was tested with the commercial kit (Mlbio, Shanghai, China) according to the instructions of the manufacturer.

Real-Time PCR Assessment

Briefly, total RNA was purified from HRGECs using the RNA extraction kits (TIANGEN, Beijing, China). cDNA (cDNA) was synthesized as described previously (5). Gene expression was measured by real-time PCR assay (Vazyme) and a DNA Engine Opticon 2 System (BioRad Laboratories, Hercules, CA). The primer sequences were described as follows:

IL-1 b: 5′‐ TTCCTGTTGTCTACACCAATGC‐3′ (F) 5′‐CGGGCTTTAAGTGAGTAGGAGA‐3′ (R);

IL-6: 5′‐TCTCTCCGCAAGAGACTTCCA‐3′ (F)

5′‐ ATACTGGTCTGTTGTGGGTGG‐3′ (R);

TNF-a: 5′‐ CCTCTCTCTAATCAGCCCTCTG‐3′ (F) 5′‐ GAGGACCTGGGAGTAGATGAG ‐3′ (R);

Actin: 5′‐ TGACGTGGACATCCGCAAAG‐3′ (F) 5′‐ CTGGAAGGTGGACAGCGAGG‐3′ (R);

Western Blot and Elisa Assay

Briefly, proteins from cells and tissues were extracted in RIPA buffer (Thermo ScientificTM, Chelmsford, MA, USA) containing phosphatases and proteases inhibitor cocktails (Sigma, St Louis, MO, USA). Proteins (20 mg) were transferred to a PVDF membrane (Millipore, IPVH00010, Massachusetts, USA) following SDS-PAGE. The slides were incubated in blocking solution (5% non-fat milk) for 60 min, and then the PVDF membrane was incubated with primary antibody without washing in a cold room overnight. After that, the process of PVDF membrane wash by tris buffered saline-tween (TBST) buffer, incubated by secondary antibody, ECL tableting, and exposed were consistent with our previous studies (5). The primary antibodies were listed as follows: anti‐GAPDH (1:1000; CST, USA), anti‐CD31 (1:1000; CST, USA), anti‐a‐SMA (1:1000;

Abcam, USA), anti‐fibronectin (1:1000; BD Biosciences, USA),

anti‐LC3 (1:1000; CST, USA), anti‐ATG16L (1:1000; CST, USA),

anti‐NF-kB p65 (1:1000; CST, USA), anti‐ Phospho-NF-kB p65 (1:1000; CST, USA). Quantification was performed by measuring the intensity of the signals with the aid of NIH image analysis software. MAP1LC3 was detected using an ELISA kit (ml000829, Mlbio, China) according to the manufacturer’s instructions.

Cell Culture, Treatment, and shRNA Transfection

Human renal glomerular endothelial cells (HRGECs) were cultured in an Endothelial cell medium (ECM, ScienCell Research Laboratories Carlsbad, CA, USA) containing 5% fetal bovine serum. Equipment Setup CO2 Incubator Culture HRGECs in a humidified atmosphere containing 5% CO2 at 37°C. Endothelial cells were preserved in a serum-free medium overnight and then treated with IL-1b, IL-6, and TNF‐a for different hours or different concentrations.

HRGECs stable ATG16L-knockdown and ATG16L- overexpression were constructed using the lentivirus (Jikai, Shanghai, China). HRGECs were first transfected with the lipofectamine 2000, then transfected with 1.5 mL ATG16L shRNA viruses or 2.0 mL ATG16L overexpressing viruses according to the manufacturer’s protocol. The endothelial cells were then switched into ECM containing 2.0 mL polybrene and incubated for an additional 6 h. We detected the infection rate of viruses with a Nikon Eclipse 80i fluorescence microscope (DS-Ri1, Nikon, Shanghai, China). Puromycin (2 µg/mL; Gibco, Thermo Fisher Scientific) was used to select stable virus-infected cells. Thereafter, HRGECs were collected for further experiments.

Statistical Analysis

GraphPad Prism 5.0 (GraphPad Software, Inc., La Jolla, CA, USA) was used for statistical analysis. Results were expressed as mean ± SD (mean ± SD) from at least three independent experiments. Comparison between and within multiple groups was performed using one-way analysis of variance followed by the Student-Newman-Keuls test. P values of <0.05 were considered significant.

RESULTS

Changes of Histomorphology and Pathology in Kidney Tissues From the Non-CAD and CAD Patients

In this study, histological examination with hematoxylin-eosin (HE) and Masson’s trichrome staining showed significant renal allograft IF/TA in CAD patients tissues (n = 4) compared to the kidney transplanted patients with stable renal function from our center (non-CAD group, n=4) (Figures 1A, B). Sample collection of materials and methods and Table 1 describes the detailed grouping information for the original features. Immunohistochemistry staining (IHC) of human renal biopsy samples demonstrated a higher expression of a‐SMA, and fibronectin (FN), and remarkably less expression of CD31 in the CAD group (Figures 1C–H). CD31 is a marker protein of endothelial cells, and endothelial cells lost their feature CD31 but acquired mesenchymal features such as a‐SMA, FN when EndMT occurs. These results indicated that the expressions of EndMT markers were notably higher in the CAD groups compared to the non-CAD group.

ATG16L Elevated After Kidney Transplantation, but ATG16L Was Downregulated in CAD Patients Comparing to the Non-CAD Patients

RNA-seq and clinical data from 25 CAD samples and 21 non- CAD patients’ renal biopsies tissue samples were downloaded from GSE9493. Patients’ information in the database was enrolled at Novartis Institutes for BioMedical Research, Switzerland, and processed for histopathology. Histological diagnosis and biopsy classification based on Banff ‘05 criteria. Demographic and clinical data are provided in the GEO database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc= GSE9493). The schematic diagram of the flow of analyses was as shown in Figure 2A. There were 1463 significant differential genes were extracted with P-value <0.05. Using the criteria for [log2 fold change (log2FC)]> 1, we selected 92 up-regulated and 115 down-regulated differential genes from the CAD group compared to the non-CAD group (Figure 2B). Then we searched genes associated with autophagy in the GeneCards database. A total of 149 ARGs with a relevance score >7 were chosen. There were 14 ARGs included in CAD and non-CAD patients' ARGs. Among them, we found that the ATG16L which got the highest relevance score gene was downregulated in CAD samples (Figure 2C).

In order to avoid the batch effect and population differences, we performed a series of experiments using human specimens from our center to confirm public database results. In the first set, we stained sections of normal, non-CAD, and CAD renal samples with an anti-ATG16L antibody, and demonstrated that the expression of ATG16L was significantly reduced in CAD patients compared with the non-CAD group (Figure 2E). The same conclusion was reached by the western blot(WB) assay. The increase of FN represented aggravated fibrosis level of transplant kidneys and one of the indicators of CAD. The WB tendency of ATG16L was congruity with the results of our IHC stain verification and consistent with the results in GEO (Figure 2D). These suggested that ATG16L expression was instead decreased with the progression of CAD. When performing Pearson’s correlation analysis, we found that ATG16L had a strong negative correlation in the expression pattern of a‐SMA, FN proteins and positive correlation with CD31 expression. ‘r’ represents the Pearson correlation value and P values denote significance of correlation (Figures 2F–I).

Autophagy Was Downregulated in CAD Patients, and Downregulation of MAP1LC3 Predicted Poor Prognosis in CAD Patients

In the previous experiment, a number of ARGs that regulate autophagy were identified from the GEO database. Among them, ATG16L is involved in the early step of autophagy, we speculate about the possibility that autophagy could be involved in CAD. To examine the levels of autophagy in human renal allograft tissues, we first found fewer autophagic vacuoles in the CAD group through a transmission electron microscope (TEM), compared to the non-CAD group (Figures 3A, B). Besides the gold standard TEM for monitoring autophagy, decreased levels of SQSTM1 and increase of LC3-II expression reflect the increase of autophagy. Similar trends of protein expression levels of LC3-I (non-lipidation of LC3), LC3-II (lipidation of LC3), and SQSTM1 in the CAD group were verified by WB assay as shown in Figures 3C, D. P62/SQSTM1 (SQSTM1) is a selective autophagy receptor and is degraded by autophagy. Autophagosomes are associated with the lipidation of the cytosolic form of LC3 and forming of LC3-II. Thus, WB results also suggested that autophagy occurs.

Some studies have shown that there are differences of ARGs in the serum of kidney transplant patients; serum LC3 level probably was the prognostic indicator correlated with the severity of the disease (25, 26). To further explore the role of autophagy in CAD procession, we examined the serum LC3 concentration of 30 non-CAD and 30 CAD patients by enzyme-linked immunosorbent assay (ELISA). The demographic data presenting the basic features of the patients was in Table 1. Distributions of the case number, patient sex, patient age, panel reactive antibodies (PRA), immunosuppressive regimen, renal function of the non-CAD and CAD groups were as shown in Table 1. Sex, age, and PRA showed no major differences. There was a significant difference in serum creatinine and blood urea nitrogen (BUN) between the two groups. Then, we tried to determine the prognostic value of autophagy in CAD patients. Through serum ELISA results, we found that the LC3 expression was lower in the CAD group. There were also a progressive CAD and a simultaneous decrease in LC3 level (Figure 3E). In addition, we divided LC3 into high and low expression groups based on the median, survival analysis, the results revealed LC3 low expression group had a shorter allograft dysfunction time (Figure 3F). Overall, we speculated that the deficiency of autophagy had a significant association with the progression of CAD.

Downregulation of ATG16L Reduced Autophagy and Promoted Progression of EndMT

To determine whether ATG16L is involved in the change of autophagy, we first transfected HRGECs with ATG16L short hairpin RNA (shRNA). TEM was used to directly demonstrate autophagic vacuoles formation and we found that there are fewer autophagic vacuoles in cells transfected ATG16L shRNA compared with those transfected with scramble shNC (Figures 4A, B). Both ATG16L protein expression and autophagy activity were also downregulated in knocking down the ATG16L group compared with shNC group (Figure 4C).

The role of ATG16L-dependent autophagy in the progression of EndMT was then investigated. HRGECs were transfected with ATG16L shRNA or shNC. Figure 4D showed that ATG16L shRNA transfection could markedly facilitate FN and a-SMA expressions and decrease CD31 expression compared with scramble shNC-transfected HRGECs. In contrast, when the ATG16L was overexpressed in HRGECs, the progression of EndMT did not significantly promote (Figure 4E). IHC staining results also confirmed WB trends for CD31 and a-SMA expressions in HRGECs (Figure 4F).

Changes of Renal Interstitial Fibrosis, Urine Protein, Renal Function and Time-Dependent Increases of EndMT in Chronic Rejection Rat Model

Next, we investigated the changes in renal interstitial fibrosis, urine protein, renal function, and EndMT in vivo. Rat renal transplanted models of chronic rejection were established by kidney transplantation from F344 rats to Lewis rats. Renal allografts were taken at weeks 8, 12, and 16 after kidney transplantation. HE staining and Masson staining assays indicated varying degrees of inflammatory cell infiltration, glomerulosclerosis, tubular atrophy, and interstitial fibrosis after kidney transplantation. The results of the Masson staining also showed a large amount of blue-stained collagen fibers in the glomerular and renal interstitium at weeks 8, 12, 16 after kidney transplantation, compared with the syn group (Figures 5A–C). In order to test the renal function of rats with an allogeneic transplant, the native right kidneys were extirpated during the kidney transplantation. Four weeks after the resection of the right kidney, 24 h urine protein and renal functions including the serum creatinine and BUN deteriorated gradually in all groups until their death (Figure 5D). The outcomes of the IHC assay revealed high expressions of FN and a-SMA while CD31 expression reduced in all groups at 16 weeks (Figure 5E). Western blot assay confirmed the outcomes of immunohistochemistry assay and showed that progression of EndMT was intensified with the extension of the time after kidney transplantation (Figure 5F).

Kidney Transplantation Activated Autophagy but Impaired Autophagic Flux in Late Stage of Chronic Rejection Rat Model

LC3, which is considered a defined marker for autophagy, is critical for autophagosome formation and the activation of autophagic flux. Adeno-associated virus (AAV)-mRFP-GFP- LC3 was injected into the rat renal allograft in vivo stereotactically. We create an mRFP-GFP-LC3 model in order to monitor autophagic activity in kidney transplantation. AAV- mRFP-GFP-LC3 is usually used to monitor autophagic flux in vivo. GFP signal has higher sensitivity to the acidic conditions of lysosomal lumen than mRFP. Thus, autophagosomes exhibited yellow fluorescence because of the co-location of GFP and mRFP fluorescence indicate and autolysosomes showed red fluorescence. An increase of red fluorescence indicated an accumulation of autolysosomes and activation of autophagy in the kidney. The positive areas of red fluorescence significantly increased at weeks 4, 8,12, and 16 in the allo group compared with the syn group. However, the red puncta weakened rapidly after reaching the peak at 8 and 12 weeks. Changes in the level of autophagy were mainly observed in the glomerular and peripheral small vessels, rather than in the distal tubules or collecting duct (Figures 6A, B). These results revealed that autophagy was activated in the early stages of kidney transplantation but not the late stages, and autophagosome clearance was impaired after about 8 weeks. Immunohistochemistry assay demonstrated that the expression of ATG16L was higher in 8 weeks than 16 weeks (Figure 6C). Western blot results also showed LC3-II and ATG16L levels peaked at 8 weeks, and then gradually decreased until 16 weeks in rat allograft kidneys after transplantation (Figure 6D).

NF-kB Pathway Was Activated and Promoted EndMT After Knockdown of ATG16L Expression

In light of the association between ATG16L downregulation and EndMT in a chronic rejection rat model and CAD disease, we hypothesized that loss of ATG16L could directly lead to the progression of EndMT. To further explore the potential mechanism by which this progression occurs, we first performed total transcriptome sequencing (RNA-seq) in HRGECs with ShNC or ShATG16L. A total of 453 differential genes (DEGs)were obtained between the two groups. Among them, 358 DEGs were up-regulated and 95 were down-regulated (Figure 7A). Then 20 pathways were significantly selected by KEGG enrichment. We found that the NF-kB pathway was one of the pathways that were significantly activated (Figure 7B). Upregulation of the NF-kB pathway genes was also confirmed by Gene set enrichment analysis (GSEA) (Figure 7C).

NF-kB is one of the crucial transcription factors for inflammatory signaling pathways, and NF-kB will translocate from cytoplasm to nucleus when it is activated. In order to confirm the NF-kB activation and nuclear translocation, we performed the immunofluorescence staining assay using anti- p65 antibody and found p65 and the nuclear dye were colocalized which indicated nuclear translocation (Figure 7D). In addition, we found that QNZ, an NF-kB signaling pathway inhibitor, significantly inhibited the expression of NF-kB in the nucleus, and reversed the expressions of EndMT related proteins induced by knockdown of ATG16L (Figure 7E).

WB assay was also carried out to verify the specific mechanism of NF-kB pathway activation. We detected quantification of p65 (total NF-kB) and p-p65 (phosphorylated NF-kB) in HRGECs. WB results showed that the p65 level was increased in ShATG16Lgroup compared to the ShNC group, while p-p65 followed a similar increasing trend (Figure 7F). The increase in p-p65 level in HRGECs after knockdown ATG16L was accompanied by an increase in total p65 protein expression, which likely reflected an additional effect on p65 expression in HRGECs. To detect whether stabilization of NF-kB is due to decreased degradation or increased synthesis, protein synthesis inhibitor cycloheximide (CHX) was added to the cells to inhibit protein synthesis and the NF-kB protein level was analyzed. Western blot results showed that there is a continued expression of p65 protein in the presence of CHX (Figure 7G). Meanwhile, NF-kB mRNA levels were not yet significantly altered by ShATG16L. The quantitative data for the levels of NF-kB mRNA were indicated graphically in Figure 7H. This suggested that NF-kB nuclear translocation was due to decreased degradation, not increased synthesis after knockdown of ATG16L in HRGECs. To further confirm the expression of NF-kB in human tissue, we performed an immunofluorescence staining assay on achieved human renal allograft tissues and found an expression level of p65 significant increase in renal allografts of CAD patient (Figures 7I, J).

NF-kB Pathway Activation Induced Cytokines Secretion, Among Them IL-1b, IL-6, and TNF‐a Could Stimulate Progression of EndMT in HRGECs

NF-kB signaling is a master regulator of the inflammatory response. We then explored the soluble factors responsible for shATG16L-induced EndMT by RT-PCR. Among the several cytokines that are relevant to the autophagy defects and inflammatory response, IL-1b, IL-6, and TNF-a, the typical downstream inflammatory factors induced by the NF-kB pathway, were significantly elevated when compared with shNC- transfected HRGECs during impairment of autophagy activity (Figure 8A). This is in line with our previous findings of inflammatory factors such as TNF-a an increase in the rat serum with chronic rejection (27).

To investigate the pathogenesis of EndMT stimulated by IL-1b, IL-6, and TNF-a, we treated the HRGECs with different concentrations of IL-1b (0~20 ng/mL), IL-6 (0~200 ng/mL) and

TNF‐a (0~20 ng/mL) for 24 hours. All of IL-1b, IL-6, and TNF-a could upregulate the protein level of FN and a‐SMA with a dose-dependent decline of the endothelial cell markers CD31 simultaneously. These effects were peaked when HRGECs were treated with 10 ng/mL IL-1b, 200 ng/mL IL-6, or 10 ng/mL TNF‐a, respectively (Figure 8B). Time dependencies experiments were also performed to validate the related expressions of EndMT. As shown in Figure 8B, IL-1b, IL-6 or TNF-a treatment for 0~48 h could remarkably induce FN and a- SMA expressions and led to the loss of CD31 in HRGECs as the stimulus time increases. To test the effects of IL-1b, IL-6, and TNF-a on EndMT, Anakinra, Tocilizumab, and Infliximab treatment were initiated when knocking down ATG16L. Anakinra is an antagonist of interleukin-1 receptor (IL-1R). Tocilizumab is an IL-6R neutralizing antibody that blocks IL-6 from binding to IL-6R. Infliximab is a monoclonal IgG1 antibody that binds specifically to TNF-a and is a TNF-a inhibitor that has been successfully used in the management of inflammatory bowel disease. WB results showed that Anakinra and Infliximab could reduce the EndMT more pronouncedly (Figure 8C). These results suggested that the loss of ATG16L was not only reduced autophagy activity but also promoted the secretion of inflammatory cytokines IL-1b, IL-6, and TNF-a by suppressing p65 degradation to activate the NF-kB pathway (Figure 8D).

DISCUSSION

In this present study, we showed that kidney transplantation impaired the autophagic flux, which was associated with an increase in the inflammatory response, EndMT, and the severity of renal allografts interstitial fibrosis. Although autophagic activity was temporarily increased at the early stages of allograft transplantation, it gradually decreased due to loss of ATG16L. ATG16L-dependent autophagy, as a cytoprotective process, could attenuate the transplanted kidney interstitial fibrosis via the regulation of the EndMT induced by IL-1b, IL-6, and TNF-a. To the best of our knowledge, this is the first study that elucidates the role and dynamic change of autophagic flux in kidney transplantation. In this study, we also investigated the molecular pathophysiologic details underlying the way in which impairment of autophagic flux results in EndMT and renal allograft interstitial fibrosis.

Although many types of research have reported that autophagy is involved in the formation and progression of renal interstitial fibrosis, ARGs have not been comprehensively analyzed to explore their clinical value in the progression of chronic renal graft rejection. Hundreds of proteins are associated with the process of autophagy. Given the importance of autophagy in kidney transplantation, it is reasonable to speculate that ARGs hold great promise in prognostic prediction and therapeutic targets. In order to explore the original pathogenetic factors of CAD, we profiled the mRNA expressions of 149 ARGs in the GSE9493 cohort and emphatically analyzed 14 ARGs, which were most associated with autophagy between CAD and non-CAD patients. Among the downregulated ARGs in the CAD group, autophagy negative regulatory genes, such as Raptor and MAPK3 (28, 29), were downregulated in the CAD group. On the other hand, ATG16L as an up-regulated ARG was a scaffold for LC3 lipidation by dynamically localizing to the putative source membranes and promoting autophagic vacuole formation (30). However, some studies have shown that ATG16L was not required for the elongation of the isolation membrane (23). Hence, the association between the autophagy-independent or autophagy-dependent role of ATG16L and EndMT was studied.

To investigate the association between ATG16L, autophagy, and EndMT in the transplanted kidneys, we knocked down ATG16L by shRNA in HRGECs. Our findings revealed that shATG16L could reduce autophagy levels and induce EndMT. Similarly, in vivo research and clinical samples also revealed that activation of ATG16L and autophagic flux in the early stage of kidney transplantation, autophagic flux declined in a dynamic range as the fibrotic changes increased. It may be concluded that the reason for EndMT progression is probably due to autophagic flux reduction, and ATG16L-dependent autophagy plays a protective role in the progression of EndMT and renal allograft interstitial fibrosis. However, some findings showed that rapamycin-induced autophagy led to activation of EndMT in HCAECs under anoxic conditions (22). These seem to be distinct from our results. We consider that role of autophagy might be different in diverse cell types. It could be cellular context- and upstream regulatory factors- dependent. Given that ATG16L-dependent autophagy is a means of self-protection when EndMT occurs in CAD, it is meaningful to determine the specific mechanism of shATG16L-induced EndMT.

We performed next-generation sequencing from the group of ATG16L knockdown and control group in HRGECs. KEGG pathways and GSEA analysis revealed that the NF-kB pathway was one of the most significantly enriched pathways. NF-kB normally binds to IkB, which stabilizes NF-kB in the cytoplasm. When the NF-kB pathway is activated, NF-kB heterodimer p65 translocates to the nucleus and binds to its specific promoter (31). One of the canonical (classical) pathways for NF-kB nuclear translocation is p65 phosphorylation. However, an increase in the number of non-canonical pathways has been reported in different diseases. It has also been reported that the suppressed p62/SQSTM1-mediated selective autophagy could reduce NF-kB pathway lysosomal degradation (32), autophagy defect also leads to NF-kB activation because of sustained SQSTM1 expression, which in turn promotes tumorigenesis in mouse models (33). In our study, we observed an increase in p65 and p-p65 protein levels while the mRNA levels were unchanged. In the presence of a protein synthesis inhibitor CHX, accumulation of the p65 protein still continues compared with the absence of inhibitor. We, therefore, speculated that NF-kB activation induced by knockdown of ATG16L resulted from reduced protein degradation and this process might be autophagic dependent degradation suppression.

Inflammations are strictly interconnected with important consequences of kidney transplants at clinical and therapeutic levels (34). The NF-kB p65 is the most important transcription factor in the regulation of cellular inflammatory factors. In this study, we also found knocking down ATG16L increased cytokine IL-1b, IL-6, and TNF-a through the NF-kB pathway. This and other findings agree with the many other studies such as the ones that concluded that loss of ATG16L enhanced endotoxin-induced IL- 1b production (35) and TNFAIP3/A20 binds ATG16L1 to control the autophagic response, NF-kB activation (36). Furthermore, the role of inflammatory cytokines on EndMT has been reported in some earlier studies. Nevertheless, it is probable these inflammatory cytokines play a diverse role depending on the underlying pathology and its ambient levels. In a diet-induced obesity mouse model, endothelial autophagy deficiency induced IL6-dependent EndMT (37). In vascular calcification disease, IL-1b and TNF-a induced EndMT in human primary aortic endothelial cells (38). Our results screened and identified IL-1b, IL-6, and TNF-a as the promoting EndMT inflammatory cytokines in pathological environments of CAD. Herewith, the reciprocal relationships between loss of ATG16L and inflammatory cytokines were researched to underline the possible therapeutic targets to control loss of ATG16L in CAD progression.

Our research has several limitations. The current validation of the effects of ATG16L on EndMT and fibrosis progression was only conducted in a knockdown cell line model, the complexity of autophagy alterations in vivo and in the development of CAD can hardly be mimicked and underlines the need for endothelial ATG16L conditional knockout mice. We will establish the ATG16L knockout mouse renal transplant model for further mechanism study. Our results showed that the loss of ATG16L and autophagy promoted EndMT, QNZ inhibited NF-kB activation, and IL-1b, TNF-a monoclonal antibody antagonized specific receptors thereby slowing down the allograft fibrotic process, but it could not completely reverse this outcome. Whether other mechanisms were playing a similar role remains unknown. Additionally, besides endothelial cells, the source of myofibroblasts consists of a variety of cells, including epithelial, fibroblasts, and pericytes (39). Whether the loss of ATG16L and autophagy deficiency affect these cells is still inconclusive.

In conclusion, the current studies proved that knockdown of ATG16L was vital for impairment of autophagic flux and EndMT formation induced by IL-1b, IL-6, and TNF‐a in the progression of transplanted renal interstitial fibrosis. This function was mediated by NF-kB reducing degradation and pathway activation. In summary, the results of our study provide novel insight into the dynamic relationship of autophagic flux and EndMT. Preventing ATG16L loss and autophagy flux inhibition could be a new option for the treatment and prevention of progression of renal interstitial fibrosis and CAD in kidney transplanted recipients.

Cistanche stem prevents kidney disease, click here to get the sample

DATA AVAILABILITY STATEMENT

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

ETHICS STATEMENT

The studies involving human participants were reviewed and approved by the local ethics committee of the First Affiliated Hospital of Nanjing Medical University. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by the local ethics committee of the First Affiliated Hospital of Nanjing Medical University. Written informed consent was obtained from the owners for the participation of their animals in this study.

AUTHOR CONTRIBUTIONS

RT and MG supervised and conceived the project. ZG, CS, and ZW designed and carried out most of the experiments. MZ and SF collected the samples of rats. LS, HC, ZH, JT, XJ, and HY analyzed the data. ZG and ZW made the figures. ZG and RT drafted and revised the paper. All authors contributed to the article and approved the submitted version.

FUNDING

This work has been supported by the National Natural Science Foundation of China [grant numbers 82070769, 81870512, 81770751, 81570676, 81470981, 81100532].

REFERENCES

1. Carney EF. Transplantation: Survival benefit of accepting a diabetic deceased donor kidney. Nat Rev Nephrol (2017) 13(8):444. DOI: 10.1038/nrneph.2017.88

2. El-Husseini A, Aghil A, Ramirez J, Sawaya B, Rajagopalan N, Baz M, et al. The outcome of kidney transplant in primary, repeat, and kidney-after-nonrenal solid-organ transplantation: 15-year analysis of recent UNOS database. Clin Transplant (2017) 31(11). DOI: 10.1111/ctr.13108

3. Schinstock CA, Stegall M, Cosio F. New insights regarding chronic antibody-mediated rejection and its progression to transplant glomerulopathy. Curr Opin Nephrol Hypertens (2014) 23(6):611–8. doi: 10.1097/MNH. 0000000000000070

4. Goldberg RJ, Weng FL, Kandula P. Acute and Chronic Allograft Dysfunction in Kidney Transplant Recipients. Med Clin North Am (2016) 100(3):487–503. doi: 10.1016/j.mcna.2016.01.002

5. Wang Z, Han Z, Tao J, Wang J, Liu X, Zhou W, et al. Role of endothelial-to- mesenchymal transition induced by TGF-beta1 in transplant kidney interstitial fibrosis. J Cell Mol Med (2017) 21(10):2359–69. doi: 10.1111/ jump.13157

6. An C, Jia L, Wen J, Wang Y. Targeting Bone Marrow-Derived Fibroblasts for Renal Fibrosis. Adv Exp Med Biol (2019) 1165:305–22. doi: 10.1007/978-981- 13-8871-2_14

7. Cruz-Solbes AS, Youker K. Epithelial to Mesenchymal Transition (EMT) and Endothelial to Mesenchymal Transition (EndMT): Role and Implications in Kidney Fibrosis. Results Probl Cell Differ (2017) 60:345–72. doi: 10.1007/978- 3-319-51436-9_13

8. Yang J, Antin P, Berx G, Blanpain C, Brabletz T, Bronner M, et al. Guidelines and definitions for research on epithelial-mesenchymal transition. Nat Rev Mol Cell Biol (2020) 21(6):341–52. doi: 10.1038/s41580-020-0237-9

9. Zeisberg EM, Potenta SE, Sugimoto H, Zeisberg M, Kalluri R. Fibroblasts in kidney fibrosis emerge via endothelial-to-mesenchymal transition. J Am Soc Nephrol (2008) 19(12):2282–7. doi: 10.1681/ASN.2008050513

10. Di Gregorio J, Robuffo I, Spalletta S, Giambuzzi G, De Iuliis V, Toniato E, et al. The Epithelial-to-Mesenchymal Transition as a Possible Therapeutic Target in Fibrotic Disorders. Front Cell Dev Biol (2020) 8:607483. doi: 10.3389/ fell.2020.607483

11.Cardinal H, Dieude M, Hebert MJ. Endothelial Dysfunction in Kidney Transplantation. Front Immunol (2018) 9:1130. DOI: 10.3389/fimmu.2018. 01130

12. Zhao Y, Cooper DKC, Wang H, Chen P, He C, Cai Z, et al. Potential pathological role of pro-inflammatory cytokines (IL-6, TNF-alpha, and IL-17) in xenotransplantation. Xenotransplantation (2019) 26(3):e12502. doi: 10.1111/Xen.12502

13. Zhao C, Xu Z, Wang Z, Suo C, Tao J, Han Z, et al. Role of tumor necrosis factor-alpha in epithelial-to-mesenchymal transition in transplanted kidney cells in recipients with chronic allograft dysfunction. Gene (2018) 642:483–90. DOI: 10.1016/j.gene.2017.11.059

14. Zhang B, Niu W, Dong HY, Liu ML, Luo Y, Li ZC. Hypoxia induces endothelial mesenchymal transition in pulmonary vascular remodeling. Int J Mol Med (2018) 42(1):270–8. doi: 10.3892/ijmm.2018.3584

15. Li J, Zhang Q, Ren C, Wu X, Zhang Y, Bai X, et al. Low-Intensity Pulsed Ultrasound Prevents the Oxidative Stress-Induced Endothelial-Mesenchymal Transition in Human Aortic Endothelial Cells. Cell Physiol Biochem (2018) 45 (4):1350–65. doi: 10.1159/000487561

16. Guan Q, Nguan CY, Du C. Expression of transforming growth factor-beta1 limits renal ischemia-reperfusion injury. Transplantation (2010) 89(11):1320– 7. doi: 10.1097/TP.0b013e3181d8e9dc

17. Levine B, Kroemer G. Biological Functions of Autophagy Genes: A Disease Perspective. Cell (2019) 176(1-2):11–42. DOI: 10.1016/j.cell.2018.09.048

18.De Rechter S, Decuypere JP, Ivanova E, van den Heuvel LP, De Smedt H,

Levtchenko E, et al. Autophagy in renal diseases. Pediatr Nephrol (2016) 31 (5):737–52. doi: 10.1007/s00467-015-3134-2

19. Li H, Peng X, Wang Y, Cao S, Xiong L, Fan J, et al. Atg5-mediated autophagy deficiency in proximal tubules promotes cell cycle G2/M arrest and renal fibrosis. Autophagy (2016) 12(9):1472–86. doi: 10.1080/ 15548627.2016.1190071

Baisantry A, Bhayana S, Rong S, Ermeling E, Wrede C, Hermann J, et al. Autophagy Induces Prosenescent Changes in Proximal Tubular S3 Segments. J Am Soc Nephrol (2016) 27(6):1609–16. doi: 10.1681/ASN.2014111059