1. Introduction
At this time, renal transplantation is the recommended treatment for end-stage kidney disease. Unfortunately, there is a persistent lack of kidney grafts, with a mismatch between the number of new patients registered on a transplant waiting list and the number of grafts available. Even though over recent years, several strategies have been drawn up to alleviate this gap with new types of organ donations from marginal donors and controlled and non-controlled donors having died after circulatory collapse, brain death donors (DBD) remain the principal source of grafts in the world [
1,
2].
However, DBD organs are exposed to injuries induced by brain death (BD) increasing their sensitivity to ischemia-reperfusion (IR) injury and compromising organ quality. Indeed, these injuries contribute to primary non-function and delayed graft function. This phenomenon is aggravated by the growing number of “expanded criteria” donors [
3]. In this context, better understanding and improvement of DBD management protocols prior to organ explantation is of utmost importance in view of upgrading organ quality. Comprehension of the different phenomena secondary to brain death occurrence and their repercussions on potential grafts is essential to the elaboration of management protocols increasing the quantity and improving the quality the donor organs susceptible to being extracted.
In the literature the usual animal brain death models are essentially mouse models characterized by brutally induced brain death, and they present limits, particularly as regards to the kidneys [
4]. At times, the data issued from these models are contradictory to findings in health clinics for humans [
5].
A recent experimental model was focused on the impact of the kinetics of brain death on the different intra-abdominal organs (kidney and liver), and the results highlighted appreciable difference [
6]. However, this study was based on a mouse model, and its results are to be interpreted with caution, especially insofar as the grafts from the donors were finally not transplanted.
Wishing to provide clinicians with potential intervention routes in brain death donor (DBD) management, we reproduced near clinical brain death (BD) conditions in a preclinical porcine model of kidney BD donors prior to transplantation. In fact, due to anatomy and physiology similar to humans, the pig is of pronounced interest as a preclinical model [
4].
Our objective was to determine the time-dependent effects of rapid brain death induction (RgBD) as opposed to slow brain death induction (SgBD) on the kidney during 4h donor management after BD declaration. We opted for four hours of management subsequent to brain death based on the results reported by Schuurs et al. [
7], who showed that after four hours, considerable expression of the genes coding for the pro and anti-inflammatory responses secondary to brain death was to be observed.
We consequently explored the impact of RgBD and SgBD on critical pathways determining kidney graft quality, specifically targeting oxidative stress [
8] and post-transplant outcome (
Figure 1).
We highlighted a difference in early function recovery and observed more severe chronic injuries in pigs grafted with kidneys from slow generated brain death (SGBD). This observation validates previous findings in mouse models and should lead to proposal of animal models better suited to the models being tested.
3. Discussion
At present we are lacking in fundamental data that would improve the treatment of brain death donors (DBD) and achieve early characterization of future kidney grafts during the phase of donor reanimation. The relevant scientific literature presents contradictory clinical and pre-clinical data featuring sizeable differences between the data obtained in humans and those in animal models. The difference may be explained by a limited choice of usable animal models, namely mouse models with rapid, even instantaneous systemic inflammation during BD [
5]. While mouse models in renal transplantation are imperfect, pig models are anatomically and physiologically more satisfactory [
4]. Moreover, one of the other limitations in pre-clinical studies in rodents consists in systemic inflammatory response during BD, which does not correspond to usual clinical situations in which several hours elapse between a patient’s brain injury and occurrence of BD. That being said, a recent study using a mouse model compared the impact of different speeds of BD induction. While the results are of interest, showing a detrimental effect of slow-speed induction of brain death on the biological parameters of renal injuries [
6], this study did not purport to evaluate post-transplant function recovery in kidney grafts. The authors show an increased oxidative stress and apoptosis in the organs of the “slow” BD induction group as compared to the “rapid” induction group. As regards plasma creatinine concentration, once again it was higher during reanimation in the slow induction than in the rapid induction animal groups. The one hemodynamic parameter considered in this study was mean arterial pressure (MAP), and its mean values did not significantly differ between the two groups.
Here, in our preclinical study of a BD pig model, we recorded comparable results, with a deleterious impact of slowly induced BD on post-transplant function recovery in kidney grafts. The originality of our work consists in (a) the use of a BD pig model, (b) the transplantation of kidneys from said models, and (c) follow-up lasting up to three months, a time span enabling us to complement data on recovery of initial kidney function with data on chronic injuries.
Our results show that primary graft function recovery was pronouncedly less satisfactory in pigs grafted with kidneys originating in the “slow” group than with pigs grafted with kidneys originating in the “rapid” group, who displayed markedly superior post-transplantation plasma creatinine levels during the first week. This difference is also found in post-transplantation kinetics up until D14. We also investigated markers of inflammation, such as IL-18, which did not demonstrate a difference between the groups; further we attempted to measure TNFα, IL-6 and IL-10, however these were not detected. This may be due either to dilution of the cytokine in the systemic blood flow or possibly to the incompatibility of the kit with the pig, a problem that is unfortunately common in this species.
In human clinical cases, the DBD parameters liable to impact primary functional recovery are: Preexisting pathological states along with age, cardiovascular comorbidity factors (expanded criteria donors) [
3], donor hemodynamic instability [
9,
10,
11,
12] and cold ischemia duration. In the controlled setting of our pre-clinical study, all of the pigs had the same age (14 weeks) and were devoid of any pathological state prior to experimentation; moreover, reanimation duration was uniform, and the different grafts underwent cold ischemia in static conservation for the same length of time. Only the BD hemodynamics was more unstable in the RgBD group, and the variations were statistically identical.
The quality of post-BD hemodynamic reanimation was identical in the two groups of animals, and the results were similar regarding the determinants of oxygen transport (cardiac output, oxygen saturation of arterial hemoglobin and partial pressure of oxygen). Moreover, the arterially extracted blood lactate in the two groups observed the same increase kinetics at the moment of brain death and after 60 min of reanimation, a result underlining a balance between the supplies and needs of oxygen in the two groups. A possible difference in oxygen transport during reanimation consequently does not explain the observed differences between the groups regarding function recovery. On the other hand, we observed a trend towards increased LDH plasma concentration in the “progressive” group; it is explained by the peculiarities of kinetics during brain death, all other factors remaining experimentally identical. This cell lesion marker was associated with a significant ASAT increase, which may function as a marker of renal tubular lesion [
13] and be considered as deterioration of the renal function in association, notwithstanding satisfactory hemodynamics, with plasma creatinine increase. When it had been decided that reanimation would go on for four hours, the decision was based on the fact that in 100% of procedures, this time lapse allowed for stabilization of a given hemodynamic state; moreover, it allowed for presentation of maximum tissular expression of inflammatory markers [
7].
We explored oxidative stress and response to oxidative stress for two fundamental reasons: (a) Oxidative stress is associated with the occurrence and aggravation of numerous chronic renal illnesses [
14]; data collected in works by Stiegler et al. highlighted activation of oxidative stress subsequent to BD occurrence [
8]. In their study, renal nitrotyrosine through immunohistochemistry testing was shown to increase after BD. Moreover, oxidative stress in a donor tends to delay function recovery [
15].
In the study by Stiegler et al., a number of markers, parameters of response to oxidative stress (superoxide dismutase 2, gluthation synthetase, peroxisome proliferator-activated receptor alpha, oxidative stress responsive 1 and gluthation peroxidase 3) were analyzed, but only as regards to their mRNA. The authors found no increased expressions of these genes in the kidneys of DBDs; this is probably due to the fact that response to oxidative stress in the short time span of reanimation is mainly induced by the activation of anti-oxidizing proteins. In our study, we focused on three markers reflecting the intensity of oxidative stress at the level of renal tissue: 4HNE, aconitase and nitrotyrosine. In order to evaluate response to oxidative stress, we analyzed the expression of anti-oxidizing proteins by Western blot (Sirtuins, PGC1α and NRF2), which is more informative than mRNA.
Our results have shown that the kinetics of BD induction exerted little or no influence on the intensity of oxidative stress; levels of 4HNE, nitrotyrosine and aconitase were similarly overexpressed in the two BD groups. However, the kinetics of BD induction seemed to influence the response to oxidative stress mediated by Sirtuin 1 and 3 and PG1alpha. However, NRF2 was similarly overexpressed in the two groups. In addition, we studied the degree of phosphorylated mTOR expression in total mTOR. mTOR is a major, indispensable actor in numerous cellular mechanisms, particularly as regards energy regulation in renal tubular cells and fibrosis generation secondary to an aggression inducing delayed graft function loss [
16]. In the explanted kidneys of the RgBD group, activation of the mTOR protein was significantly reduced. While mTOR modulation by the ROS produced during oxidative stress is confirmed by the relevant scientific data, it can take on the form of either upregulation via AMPK [
17] or of downregulation via FoxO3 [
18]. On this topic, recent publications have highlighted direct regulation of mTOR by NRF2 [
19] as well as double downregulation involving mTOR and NRF2 [
20]. mTOR detection is also concordant with our apoptosis results [
21], which demonstrated higher apoptosis in the RgBD.
Our hypothesis was that the observed differences in primary function recovery in kidney grafts according to BD kinetics was associated not with a difference in oxidative stress intensity, but rather with modulation of anti-oxidative response intensity leading to modulation of mTOR by NRF2. This difference in secondary adaptation to an aggression corresponds to the phenomenon of allostasis, i.e., adaptation of a system to new environmental conditions. To our knowledge, this is the first time that a donor allostasis mechanism has been described. The observed difference in allostatic adaptation could originate in the protective phenomena initiated by shear stress, insofar as it has been demonstrated that sheer stress has a protective effect, particularly through eNOS activation [
21], and that the effect may entail the intervention of SIRT-1 [
22,
23]. Part of the originality of this paper consists in our having highlighted more pronounced expression of the ratio of phosphorylated eNOS to total eNOS in the rapidly generated brain death (RgBD) group. We may thus hypothesize that post-BD, the renal cells were subject to intensified shear stress involving SIRT-1 activation, which might help to explain our results. These organ level modifications of vascular homeostasis did not have observable impact on the post-BD hemodynamic parameters, possibly because such local change is unlikely to have measurable effect at the systemic level, especially measured by the typical tools of hemodynamics monitoring, which lack precision. Finally, the hemodynamic phenomena secondary to BD were more condensed in the RgBD group (5 min) than in the in the “slow” SgBD group (90 min), which strengthens the hypothesis of an involvement of shear stress.
Observation of eNOS phosphorylation and modulation of vascular homeostasis, which is of better quality in the RgBD group, is in line with the observed difference in fibrosis shown herein: As demonstrated in earlier studies from our lab [
24,
25], ischemia reperfusion targets mainly micro vessels, with deleterious consequences on the long term, including fibrosis development.
The potential interest of this work consists in its possibly explaining, for reasons other than the imperfect use of mouse models, the divergencies of the results reported in the literature. Indeed, analysis has often been based on evaluation of intensity of aggression rather than intensity of the response to the aggression of the regulation/response systems.
This work presents some limits, particularly as regards extrapolation of the results to humans and to the kidneys of donors presenting preexisting lesions, given the fact that mTOR is regulated in different ways in different renal pathologies [
16]. Even if our results possess some external validity insofar as they confirm those of Schuur et al., other experimentations remain necessary.
Were our results to be confirmed, several possible applications could be envisioned, particularly therapeutic utilization of the mTOR and/or NRF2 modulator in a donor having died following brain death. Another possible aspect would consist in early characterization of renal grafts through analysis of mTOR and NRF2 expression. While access to this type of information would necessitate an invasive approach, one possible track would consist in analysis of renal cellular material in donor urine and blood; relevant assessments have been carried out in cancerology [
26].
4. Materials and Methods
4.1. Experimental Model
All animal experiments were performed in accordance with the Care and Use of Laboratory Animals and Animal Research Reporting of In Vivo Experiments (ARRIVE) guidelines and approved by French Poitou-Charentes ethics committee of animal experimentation (protocol number CE2012-14, 15 November 2013). We used three-month-old Large White pigs weighing 40 ± 4 kg (MOPICT IBiSA Plateforme, INRA Magneraud, Surgères, France).
4.2. Brain Death Model
12 pigs (large white, 35–52 kg) were utilized for the purposes of the study, which was conducted following the agreement of the Comité d’éthique de Poitou-Charentes (C2EA-84), protocol CE2012-14 in the surgical experimentation platform MOPICT (INRA, domaine du Magneraud, Surgères, France). Pre-medication was carried out using Midazolam IM 0.05 mg/kg 20 min prior to initiation of the experimentation. The pigs were anesthetized first by inhalation of a mixture of Sevoflurane 8% and oxygen. Following loss of consciousness, a 20 G catheter was set up on an ear vein and the anesthesia was deepened by injection of propofol (3 mg/kg), fentanyl (1 μg/kg), ketamine (1 mg/kg) and rocuronium (1 mg/kg). Respiratory tract management was ensured by installation of orotracheal intubation, size 6.5 (Hi-Contour, Mallinckrodt Medical™, Athlone, Ireland). The pigs were ventilated with a volume-controlled Servo 900C™ (Siemens, Solna, Sweden) with tidal volume set at 8 mL/kg, and the respiratory rate was adjusted to obtain PaCO2 between 35 and 40 mmHg and FiO2 at 50%. Positive end-expiratory pressure was 5 cm of H20, an I/E ratio set at 1/3. No alveolar recruitment maneuver was carried out. Maintenance was achieved through iterative administration of propofol at 3 mg/kg/h, fentanyl 1.5 μg/kg/h and rocuronium 1 mg/kg according to the animal’s behavior. A femoral arterial catheter and a central venous line in the internal jugular vein put into place surgically, and monitoring was then ensured by PICCO1® Edward (Pulsion medical system, Irving, USA). A supra-pubic catheter was installed following vesical incision with hourly diuresis measurement. Transition to brain death (BD) was provoked by carving a hole at 10 cm from the coronal plane and 10 cm from the sagittal line at the level of the frontal zone, and afterwards by detachment of the dura mater and insertion of a Foley catheter Ch 18 in the epidural space; the catheter was then swollen with saline according to the Purins model. For the six animals in the SgBD (slowly generated brain death) group, brain death was provoked slowly by inflation of the Foley catheter balloon in accordance with a sequence of 5 mL at 0 min, followed by deflation of 5 mL at 10 min, inflation of 10 mL at 20 min, deflation of 5 mL at 30 min, inflation of 10 mL at 40 min, deflation of 5 mL at 50 min and inflation of 10 mL at 60 min up until brain death. For the six animals in the RgBD (rapidly generated brain death) group, brain death was rapidly provoked by inflation of the Foley catheter balloon of 20 mL in 5 min. The animals were monitored by a double-derivation electroencephalogram (EEG). Brain death diagnosis was carried out on the basis of clinical examination: A reactive bilateral mydriasis, absent corneal reflexes, abolition of spontaneous ventilation after hypercapnia response test and flat line EEG. Sedation was discontinued once brain death was validated. Reanimation consisted in maintaining mean arterial pressure above 70 mmHg, with normal hematosis under mechanical ventilation (tidal volume 8 mL/kg) and diuresis between 10 and 15 mL kg/h. Vascular filling with crystalloid solutions was guided by injection volume variation (IVV). Hemodynamic support was achieved by continuous intravenous infusion of noradrenaline. Occurrence of polyuria, which was attributed to central diabetes insipidus during brain death, was treated by low-dose desmopressin. A need to increase noradrenalin doses to a level exceeding 0.3 μg/kg/min occasioned administration of hydrocortisone at 1 mg/kg/h. Reanimation following brain death proceeded for 4 h.
Blood collections were performed at baseline time (before BD induction), at the time of BD declaration (BD time), and during donor management at 60 min after BD declaration (60 aBD), 120 min (120 aBD), 180 min (80 aBD) and 240 min (240 aBD). At 240 min after BD (4 h management) kidneys were collected. The left kidney was flushed with a cold University Wisconsin (UW) solution and preserved in cold UW in static condition during 18 h before transplantation. The right kidney was collected for biological analysis with biopsies at the end of donor management and at the end of kidney preservation in UW as previously described. The biological analyses of BD kidney were compared to control kidneys collected on three-month-old healthy pigs.
4.3. Kidney Allo-Transplantation and Renal Function Evaluation After Transplantation
At the end of static preservation, the left kidney was allotransplanted in a fully swine leukocyte antigen (SLA) compatible littermate, bilaterally nephrectomized pig (six with grafts of fast induction group and six with grafts of slow induction group). Heterotopic allotransplantation was performed in anaesthetized pig via midline incision. End-to-side aorta and inferior vena cava anastomoses were performed, just above the iliac bifurcation. Ureteroneocystostomy was carried out. After kidney transplantation, the pigs were placed in a metabolic cage for urine and blood collection during the first seven days and at days 11, 14 and 90 after transplantation (
Figure 6). Blood collections were performed at day −1 (1 day before transplantation), 60 min after kidney reperfusion and at days 1 (D1), D3, D5, D7, D11, D14 and D90 after kidney transplantation. Follow-up was performed for 90 days, animals were anesthetized at D90 post-transplantation and the transplanted kidneys were collected for histological analysis.
4.4. Soluble Protein Quantification and Renal Function Evaluation
Blood NGAL, HMGB1 and IL18 proteins were assessed in porcine plasma using ELISA kit (respectively, Eurobio, Euromedex and Invitrogen France) according to the manufacturer’s instructions. Blood creatinine, sodium, lactate dehydrogenase (LDH) and aspartate aminotransferase (ASAT) were measured with a modular bioanalyzer (Roche-Diagnostics, France). Diuresis was evaluated by urine production after transplantation. Glomerular filtration rate (GFR) was calculated as: ((Diuresis mL/h)*urine creatinine)/blood creatinine = GFR L/h, and converted in GFR mL/min.
4.5. Histological Study
Biopsy samples from corticomedullary kidney collected at day 90 after transplantation were fixed with 4% formalin and paraffin-embedded. Histological interstitial fibrosis evaluated by red Sirius staining (collagen fibers I and III staining) was quantified with ImageJ software. Renal tubulitis and leukocyte infiltration were scored on histological HES staining.
4.6. Real-Time Quantitative PCR
Total RNAs from porcine cortical renal tissue were extracted with a commercial kit (Macherey-Nagel, Hœrdt, France). Genomic DNA was removed using DNA-free kit and first-strand reverse-transcription (Applied-Biosystems, Saint-Aubin, France) was performed. Real-Time qPCR assays were performed on a RotorGene Q (Qiagen, Courtaboeuf, France). Porcine DNA primers were designed using OligoPerfect™ (Life Technologies), QuantPrim (Max-Planck-Gesellschaft, Potsdam, Germany) and OligoAnalyser (Integrated DNA Technologies, Coralville, USA); the sequences are detailed in
supplementary Table S1. RNA was expressed as “relative to expression in healthy ctl kidney” determined with the Pfaffl method (expressed as the relative fold change), using
β-actin, CYA62, SDHA, RPLPO and
L19 genes as internal controls.
4.7. Western Blot Analysis
Western Blotting was carried out using specific antibodies; Bcl-2 (D17C4, Cell Signaling, France), Bax (ab115193, Abcam, France), aconitase 2 (ab110321, Abcam, France), nitrotyrosine (05-233, Millipore, France), 4HNE (ab5605, Millipore, France), sirtuine-3 (D22A3, Cell signaling, France), sirtuine-1 (SC-9475, Cell signaling, France), phosphorylated sirtuine-1 (Sc-2314, Cell signaling, France), peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC1a, ab54481, Abcam, France), nuclear factor (erythroid-derived 2)-like 2 (NRF2, ab31163, Abcam, France), mammalian target of rapamycin (mTOR, 2983S, Cell signaling, France), phospho mTOR (5536P, Cell signaling, France), eNOS total (610296, BD Transduction, Le Pont de Claix, France) and eNOS phophorylated Ser1177 (C9C3; Cell-Signaling, Saint Quentin-Yvelines, France). Protein quantification level was expressed as ratio of total proteins loaded on precast gels with stain-free technology (BioRad, France). All blots were obtained by imaging with chemidoc and quantificated with ImageLab software (BioRad).
4.8. Statistical Methods
Kinetic results are expressed in mean ± SD, and other results are expressed as a scatter plot with the median. Statistical analysis was performed with a Kruskal–Wallis test and Dunn’s post-test for kinetic parameters, and for other parameters statistical analysis was performed as Mann–Whitney test, with R and NCSS software, statistical significance was set at p < 0.05. n = 6.