Figures
Abstract
Experimentally renal tissue hypoxia appears to play an important role in the pathogenesis of chronic kidney disease (CKD) and arterial hypertension (AHT). In this study we measured renal tissue oxygenation and its determinants in humans using blood oxygenation level-dependent magnetic resonance imaging (BOLD-MRI) under standardized hydration conditions. Four coronal slices were selected, and a multi gradient echo sequence was used to acquire T2* weighted images. The mean cortical and medullary R2* values ( = 1/T2*) were calculated before and after administration of IV furosemide, a low R2* indicating a high tissue oxygenation. We studied 195 subjects (95 CKD, 58 treated AHT, and 42 healthy controls). Mean cortical R2 and medullary R2* were not significantly different between the groups at baseline. In stimulated conditions (furosemide injection), the decrease in R2* was significantly blunted in patients with CKD and AHT. In multivariate linear regression analyses, neither cortical nor medullary R2* were associated with eGFR or blood pressure, but cortical R2* correlated positively with male gender, blood glucose and uric acid levels. In conclusion, our data show that kidney oxygenation is tightly regulated in CKD and hypertensive patients at rest. However, the metabolic response to acute changes in sodium transport is altered in CKD and in AHT, despite preserved renal function in the latter group. This suggests the presence of early renal metabolic alterations in hypertension. The correlations between cortical R2* values, male gender, glycemia and uric acid levels suggest that these factors interfere with the regulation of renal tissue oxygenation.
Citation: Pruijm M, Hofmann L, Piskunowicz M, Muller M-E, Zweiacker C, Bassi I, et al. (2014) Determinants of Renal Tissue Oxygenation as Measured with BOLD-MRI in Chronic Kidney Disease and Hypertension in Humans. PLoS ONE 9(4): e95895. https://doi.org/10.1371/journal.pone.0095895
Editor: Jaap A. Joles, University Medical Center Utrecht, Netherlands
Received: December 14, 2013; Accepted: April 1, 2014; Published: April 23, 2014
Copyright: © 2014 Pruijm et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by research grants from the Swiss National Science Foundation (FN 32003B-132913 and 149309) and by the Centre d'Imagerie BioMédicale (CIBM) of the University of Lausanne (UNIL). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Numerous experimental studies have suggested that disturbed oxygenation plays a role in the development and progression of kidney disease including hypertensive nephropathy[1]–[3].Thus, using micro-electrodes for direct pO2 measurements, low cortical O2 levels have been found in spontaneous hypertensive rats and in rats with streptozotocine-induced diabetes or in the subtotal nephrectomy model [4], [5].
Until recently, human data were largely lacking mainly due to the lack of non-invasive methods to estimate renal tissue oxygenation. Over the last decade, blood oxygenation level- dependent magnetic resonance imaging (BOLD-MRI) has become a powerful tool to estimate renal tissue oxygenation non-invasively in humans. The basic principle of BOLD-MRI is that changes in renal tissue desoxyhemoglobin concentrations involve generation of phase incoherence of magnetic spins, leading to an increase in apparent relaxation rate R2* (expressed in sec−1). Under the assumption that blood pO2 is in equilibrium with tissue pO2, R2* values provide estimates of tissue oxygenation, a low R2* indicating a high tissue oxygenation [6]. Using post-processing programs, several circles - called regions of interest (ROI's) - are placed manually per slice in the cortex and in the medulla. This allows the assessment of the average cortical and medullary R2* value, per kidney or for both kidneys together, without the need to administer contrast product [7].
Several studies have used BOLD-MRI in humans to investigate renal oxygenation in different forms of chronic kidney disease (CKD) [8], [9]. The largest study which assessed renal oxygenation at different degrees of kidney dysfunction was published by Michaely et al and included 400 patients [10]. Interestingly, no correlation was found between cortical and medullary R2* values and the estimated glomerular filtration rate (eGFR, according to the MDRD-formula) in this study. However this evaluation had several limitations: renal function was not measured in all patients (280 with measured creatinine levels), no information was collected on medication intake or baseline characteristics such as blood pressure or underlying renal disease, and BOLD-MRI was not performed under standardized conditions of fluid and sodium intake which are known to have a strong influence on renal oxygenation measured by BOLD-MRI [11], [12].
Changes in renal tissue oxygenation may also contribute to the development of ischemic and hypertensive nephropathies [5]. For these reasons, most studies performed in hypertensive patients so far have been conducted in patients with renal artery stenosis [13], [14]. Other studies in patients with essential hypertension have focused on the role of sodium, angiotensin II or race on cortical and medullary oxygenation [12], [15], [16]. These studies were performed in a relatively small number of participants, or did not include a control group. Whether or not essential hypertension is characterized by chronic hypoxia in humans remains therefore an open question.
The aim of the present study was therefore to assess renal tissue oxygenation at baseline and after an acute administration of furosemide in patients with various levels of renal function as well as in patients with essential hypertension and a normal renal function and in control subjects. Moreover, we analyzed the potential determinants of renal cortical and medullary tissue oxygenation in these patient groups.
Methods
This research project was approved by the local institutional review committee (Ethical Committee of the Canton de Vaud, Switzerland) and conducted according to the principles expressed in the Declaration of Helsinki. Written informed consent was obtained from each participant.
Subjects
Patients with CKD stage 1–5, or with hypertension without CKD were eligible for this study. CKD was defined as an estimated glomerular filtration rate (eGFR) ≤60 ml/min/1.73 m2, or the presence of structural or functional abnormalities for at least three months [17]. Arterial hypertension (AHT) was defined as mean office blood pressure (BP) ≥140/90 mmHg measured at more than one occasion, or an office BP <140/90 mmHg while taking one or more antihypertensive drugs. Other inclusion criteria were: age ≥18 years and the ability to understand the study protocol and to sign an informed consent. Controls were normotensive, untreated healthy individuals without a history of kidney disease or hypertension or any other concomitant disease. Exclusion criteria for all participants were: a contra-indication to MR-imaging such as claustrophobia or the presence of a pacemaker or other implanted metallic device.
Study protocol
Patients were recruited at the outpatient nephrology and hypertension clinic of the university hospital in Lausanne (CHUV). Controls were recruited by local advertisement. Participants were maintained on their regular diet. Salt intake, proteinuria, and creatinine clearance were measured before BOLD-MRI by a 24 h urine collection. On the day of each BOLD-MRI measurement, an identical oral hydration protocol was followed by each participant at home (loading dose of 5 ml/kg of water at 8am, followed by 3 ml/kg every hour till 12am), see extended methods for further details and justification of this hydration protocol. Subjects joined our research unit at 11.30 am. BP was measured three times by an experienced research nurse using an automated Omron 705IT oscillometric device according to the recommendations of the European Society of Hypertension [18]. Peripherical oxygen saturation was measured using a fingertip pulse oxymeter (Oxy, Medair AB, Delsbo Sweden), and blood was drawn by a catheter inserted into an antecubital vein.
BOLD-MRI
BOLD-MRI was performed between 1 and 2 pm in the radiology department. Magnetic resonance measurements were carried out on a 3T whole-body MR system (MAGNETOM Trio, Siemens Medical Systems, Erlangen, Germany), as described previously [12], [19], [20]. Briefly, four coronal slices with good cortico-medullary differentiation obtained in expiration were selected from morphological images for functional evaluation with BOLD-MRI. Twelve T2*-weighted images were recorded within a single breath-hold of 12.4 seconds with a modified Multi Echo Data Image Combination sequence (MEDIC) with the following parameters: repetition time (TR) 68 ms, echo time (TE) 6–52.2 ms (equidistant echo time spacing 4.2 ms), flip angle 20°, field of view (FOV) 400×400 mm2, voxel size 1.6×1.6×5 mm3, bandwidth 700 Hz/pixel, matrix 256×256 (interpolation to 512×512). All images were exported for analysis with a home-built IDL program (Interactive Data Language, Boulder, CO, USA). R2* maps were calculated voxel by voxel using a Levenberg-Marquardt least-squares algorithm to fit an exponential function to the signal intensities measured for each echo time. ROIs were traced in the form of circles of equal size (containing approximately 20 voxels each) in the medulla and the cortex (two in the cortex and two in the medulla in each kidney), as illustrated in Figure 1. The reported R2* value was the mean value of four slices (16 ROIs) for the medulla and the cortex. This technique has been shown to have a good reproducibility [7], [20]. Medullary-cortical R2*ratio (MCR2*), as a marker of metabolic workload, was calculated for each participant by dividing the mean medullary R2* by cortical R2* levels [21]. The procedure was repeated for all four coronal series obtained fifteen minutes after the administration of 20 mg furosemide intravenously.
Anatomical templates are shown on the left, R2* maps in the middle, and color maps on the right (low R2* levels corresponding to higher tissue oxygenation in red, high R2* levels in yellow). Two regions of interest (ROIs) are traced in the form of circles (20 voxels each) in the cortex and medulla of each kidney; this procedure is repeated on four different slices.
Statistics
Clinical data were analyzed using STATA 11.0 (StataCorp, College Station, Texas, USA). Quantitative variables were expressed as mean ± standard deviation, or as median (25th–75th percentile range), as appropriate. Qualitative variables were expressed as number of patients and percentage. Comparisons between baseline characteristics of study groups were analyzed with ANOVA. Distribution of variables was also expressed using the probability density function according to Kernel. In case of non-normal distribution, variables were log-transformed. Multivariable logistic regression, adjusting for the predefined variables age, sex, current smoking, body mass index, diabetes, hemoglobin, and 24 h urinary sodium excretion, was used to determine the independent association of cortical and medullary R2* values with CKD status. These variables were selected based on their known association with CKD (body mass index, diabetes, age), or because of a theoretical association with renal tissue oxygenation (hemoglobin, 24 h urinary sodium excretion, smoking).
Similarly, in order to determine the independent association of each independent variable of interest with the dependent variables cortical and medullary R2* values, multivariate linear regression was performed, adjusted for the same predefined covariates.
Associations between the eGFR slope and R2* values were examined with Spearman's rank correlation and multivariable linear regression analysis. Results of the logistic multivariate analysis are presented as Odds-ratio (OR) and 95% confidence interval (95% CI). Results of all linear multivariate analyses are presented as beta-coefficients (β) and their 95% confidence intervals. P values were derived from maximum likelihood ratio tests. Statistical significance was considered for a two-sided P<0.05.
Results
Clinical characteristics of the patients and controls
Details of the screening procedure are provided in the 'extended Methods and Results S1' section. Baseline clinical characteristics of patients with CKD and hypertension and of controls are shown in Table 1. Patients in the CKD-and hypertensive group were well-matched for age, sex and other baseline characteristics. Healthy volunteers were significantly younger, and more often female. Medical treatments of subjects according to group are shown in Table 2 and across different stages of CKD in Table S1.
R2* levels in hypertension and CKD as compared with healthy controls
There were no differences in mean or median cortical and medullary R2* levels between the right and left kidney in all groups (cortex (median (range)) respectively 17.4 (16.1; 18.6) and 17.6 (16.3; 18.8) sec−1, p = 0.29; and medulla (mean±SD): 28.7±4.2 vs 28.7±2.5 sec−1, p = 0.94). Therefore, only the mean cortical and medullary R2* values of both kidneys are shown and used for statistical analysis. The distribution of medullary R2* was comparable in all groups, but that of cortical R2* values differed markedly between groups (Figure 2A and B). A bimodal cortical distribution was observed in controls and AHT patients, yet not in CKD patients (Figure 2B). The bimodality of the probability density function of cortical R2* values was partly explained by gender (Figure 3); no gender-differences were seen in the distribution of medullary R2* values (Figure S1). The mean medullary and cortical R2* values and medullary-cortical ratio's of each group were comparable between the CKD, hypertensive and control group (Table 3). No differences were seen in cortical and medullary R2* levels between stage I–V of CKD (Figure 4 and Figure S2 for results by gender).
A/medullary R2* values and B/cortical R2* values.
The number of subjects was n = 10 (for the eGFR>125 ml/min/1.73 m2 category), 64 (eGFR 90–125), 65 (60–89), 36 (30–59), 15 (15–29) and 5 (<15), respectively.
In addition, within the CKD group no significant difference was observed according to the renal diagnosis (Figure 5).
Values shown as boxplots; No (no CKD, n = 100), DM (diabetic nephropathy, n = 20), AHT (hypertensive nephropathy, n = 31), GN(glomerulonephritis, n = 17), reflux (reflux nephropathy, n = 6), one (solitary kidney, n = 7), IF (interstitial nephritis, n = 6), other (other cause of kidney disease, n = 8). There were no differences between cortical and medullary R2* levels between the groups (ANOVA p = 0.10 and p = 0.99, respectively)
Changes in R2* after furosemide in hypertension or CKD patients as compared with controls
The situation was different in acute conditions, since the response to furosemide differed significantly between groups. Medullary R2* but not cortical values decreased significantly (suggesting an increase in medulla oxygenation) in the three groups, both in men and women (Table 3) but the reduction was significantly smaller in the CKD and AHT groups (Figure 6). Although there was a trend towards larger furosemide-induced decreases in medullary R2* in men, this difference was not statistically significant. In age-and gender adjusted multivariable linear regression analysis the furosemide-induced change in medullary R2* (f-R2*) correlated significantly with eGFRmdrd (adjusted β per ml/min/1.73 m2: −0.03 (95% CI 0.01; 0.05), p = 0.003); the latter association persisted when intake of loop diuretics was introduced in the model (see extended Methods and Results S1).
The difference in furosemide-induced change between the AHT and control group was not explained by a difference in eGFR. In multivariate logistic regression adjusted for age, sex, and eGFRmdrd, the furosemide-induced change in medullary change in R2* level remained significantly smaller in the presence of AHT (OR 0.81; 0.64–0.98, p = 0.047).
Associations between renal tissue R2* levels, CKD status, hypertension, and kidney function
In multivariate logistic regression, adjusted for age, sex, body mass index, hemoglobin, current smoking and urinary sodium excretion, with CKD status as dichotomized dependent variable (CKD compared with healthy controls and hypertensive subjects pooled together), medullary R2* levels were not associated with CKD status (OR(95% CI): 1.01 (0.86;1.18), p = 0.92). Medullary R2* levels were also not associated with CKD status when comparing the CKD group separately with healthy controls (CKD vs controls, OR: 0.95 (0.79; 1.15), p = 0.59), or with the hypertensive group (CKD vs. AHT, OR: 1.1 (0.91;1.35), p = 0.3).
Although cortical R2* levels were higher in CKD patients, this difference was not statistically significant in adjusted models (CKD versus non-CKD, OR: 1.20 (0.95; 1.51, p = 0.14)). The same trend was seen when comparing the CKD group separately with healthy controls (OR: 1.10 (0.85; 1.41), p = 0.49).
The determinants of renal cortical and medullary oxygenation were also analyzed in the entire population using linear models. There was no association between medullary R2* levels and eGFR or any of the other predefined variables in the adjusted multivariable linear regression analysis (Table 4A). No significant difference was observed when antihypertensive treatment (renin-angiotensin system (RAAS) blockers or antihypertensive treatment in general) was added to the model. Similar results were found for the association between medullary R2* values and creatinine clearance as measured by 24 h urinary collection (adjusted β per ml/min: −0.01, p = 0.96).
Results for cortical R2* levels and its associations with predefined variables were also largely negative (Table 4A). There was a positive association between gender and R2* levels (adjusted β −1.8, p = 0.01). Hence, male participants had higher cortical R2* levels, corresponding to lower cortical oxygenation.
Since only a small proportion of the variability in R2* levels was explained by the predefined variables, some additional analyses were performed on the entire population. Thus, a negative correlation was found between medullary R2* levels and mean arterial blood pressure (Table 4B). However, the correlation between mean arterial pressure (MAP) and medullary R2* values was weak and disappeared when including antihypertensive treatment in the model (β −0.04 per mmHg, p = 0.043 before, and β −0.03, p = 0.13 after adjustment for antihypertensive treatment). Results of multivariable linear regression analysis are shown in more detail in Table S2, S3, S4. Diabetic subjects (n = 32) had higher cortical R2* values than non-diabetic subjects (mean R2* 18.7±5.9 vs 17.5±1.6 sec−1, p = 0.05), yet diabetes was not associated with cortical oxygenation in adjusted logistic regression models (Table S3). In contrast, cortical R2* levels correlated positively and significantly with circulating glucose levels as measured just before BOLD-MRI (Table 4B and Table S4). Including antidiabetic treatments in the multivariate regression model did not alter this result.
The positive correlation found between uric acid and cortical R2* levels did not change after the inclusion of allopurinol and diuretic treatments, or mean arterial blood pressure in the model. Results for cortical R2* levels stratified by sex are shown in Table 5. Similar positive associations were found between cortical R2* levels and uric acid; the association between cortical R2* levels and glycemia were only present in men, whereas there was a positive association between cortical R2* and age in women yet not in men.
In order to assess whether some of the results were linked to an over-representation of young women in the control group, sensitivity analyses were performed in a subgroup excluding healthy women <40 years (n = 13) and men aged >75 years suffering from CKD (n = 5). The main results of the study remained unchanged (see Methods and Results S1).
Discussion
The main findings of this study are that: 1) mean cortical and medullary R2* values as a proxy for renal tissue oxygenation are similar in hypertensive patients, CKD patients and healthy controls; however, the distribution of cortical R2* values differs markedly between groups, 2) the medullary R2* response to furosemide is blunted in hypertensive patients and markedly reduced in CKD patients, 3) baseline renal tissue oxygenation appears to be remarkably stable over different degrees of kidney dysfunction, independently of the cause of kidney disease and 4) cortical R2* levels are positively associated with male gender, glycemia and uric acid levels.
The first interesting observation of this paper is that although mean cortical and medullary R2* look identical in hypertensive and CKD patients and controls, the distribution of cortical R2* differs between groups with a clear bimodal distribution in controls and AHT and a unimodal shape of the distribution in CKD. As found in our population, part of the bimodal distribution is linked to the male/female ratio. However, this ratio is similar in hypertensives and CKD patients suggesting other mechanisms explaining the differences in distribution, which should be the subject of further study.
The second finding of our study is that renal cortical and medullary oxygenation as measured by BOLD-MRI appears to be extremely well maintained in patients with CKD. Indeed, in contrast to what has been observed experimentally we did not find any decrease in cortical or medullary oxygenation even in advanced CKD; the nature of the underlying renal disease does not appear to play a major role. In this respect, our data are in accordance with the study by Michaely et al [10], although the latter study could be criticized for some methodological issues. Despite the fact that our study has been performed under very standardized hydration conditions, and detailed information was available on possible confounders such as medication, sodium intake and the type of underlying kidney disease we were not able to demonstrate a reduction of tissue oxygenation as renal disease progressed. In addition, this study also found no alterations in renal R2* values in persons with essential hypertension in comparison with healthy controls, despite the fact that numerous previous animal studies have reported renal tissue hypoxia in AHT [5].
The discrepancy between animal studies and BOLD-studies in humans regarding oxygenation can be interpreted in different ways. First of all, it might be that BOLD-MRI is not sensitive enough or simply not as good a tool to assess renal tissue oxygenation in CKD-patients as is direct invasive measurements using microelectrodes. Nonetheless, early animal studies performed to validate the BOLD-MRI technique have reported linear relationships between directly measured renal pO2 values and the BOLD signal [22]. Alternatively, the ROI-technique used to analyze BOLD-images might be less applicable in CKD patients, where loss of cortico-medullary differentiation hampers manual placement of ROIs. This might play a role in the interpretation of medullary R2* values, but cortical R2* values are expected to be reliable, thanks to its anatomic proximity with the kidney capsula. Besides, the majority of subjects in the present study had less advanced stage I–III kidney disease, and largely maintained cortico-medullary differentiation; in patients with preserved differentiation, BOLD-MRI is a reliable, reproducible method. Nevertheless, it remains an open question if studies using different methods of analysis such as the recently described compartmental model [23] would obtain similar results. Assuming that BOLD-MRI is able to correctly estimate tissue oxygenation in humans, our findings together with those of Michaely et al put into question whether ‘chronic renal hypoxia’ truly exists in humans. In conditions of acute hypoxia, hypoxia-inducible factors (HIF) are stabilized, which stimulates erythropoietin production, angiogenesis and metabolic reprogramming, offering protection and increasing survival of renal tubular cells under hypoxic conditions [2], [24]. In contrast, chronic HIF stabilization stimulates interstitial fibrosis, a process that is generally believed to compress peritubular capillaries, which further decreases local tissue oxygenation and leads to a vicious circle of HIF-stabilization, increased formation of interstitial fibrosis, glomerulosclerosis and renal function decline [3], [25]. However, hypoxia-induced fibrosis may not be a vicious circle that further worsens oxygenation, but merely a way to maintain renal oxygenation by adapting oxygen consumption to the demand. Identical mechanisms might operate in animals. Animal models reporting hypoxia in CKD such as the 5/6 nephrectomy model may not have been carried out long enough to simulate the long-term adaptive changes in the kidney that might occur after several decades of chronic kidney disease or hypertension in humans. In line with this hypothesis, Priardarshy et al [26] evaluated renal oxygenation six to eight weeks after remnant kidney creation (instead of the 2 weeks applied by Manotham et al.)[1], and found that renal oxygenation was not decreased, but rather increased in the remnant kidney. A similar observation was made in renal artery stenosis with a reduction of renal tissue oxygenation acutely but a good maintenance of tissue oxygenation in the chronically stenotic kidney [14].
A role of renal handling of sodium in mediating oxygen consumption is supported by our observation that the medullary R2* response to furosemide differs between controls and hypertensives with a blunted response in hypertension and an even more marked reduction in CKD patients. In CKD, the markedly reduced response to furosemide can be explained by the reduced renal function leading to lower concentrations of furosemide within the kidney. However, this cannot be the explanation for hypertensive patients who had a comparable renal function as controls. Pratt et al have previously demonstrated ethnic differences in the response to furosemide [27]; to the best of our knowledge, differences in response to furosemide between normo-and hypertensive white subjects have not been reported previously. The blunted response to furosemide observed in hypertensive patients may be an indirect marker of the alterations in renal sodium handling in hypertension. Persons with essential hypertension have an increased proximal tubular reabsorption and a reduced distal delivery of sodium, which might blunt the effect of furosemide [28]. Alternatively, there might be differences in mitochondrial metabolism and oxygen consumption in the thick ascending limb of Henle between hypertensive and normotensive subjects, in analogy with recently described differences between Dahl salt-sensitive rats and salt-resistant control strains [29].
Our multivariate analysis enabled us to identify several new factors associated with renal tissue oxygenation. Thus, cortical R2* levels was positively and strongly associated with male gender. The relationship with male gender was robust and persisted in sensitivity analyses, and suggests that cortical oxygenation might be regulated differently in men and women. It may also provide some clues why renal function declines faster in men. However, our data do not offer any explanation for the higher R2* values in men, and studies measuring simultaneously renal perfusion and tubular sodium handling are necessary to clarify this issue.
Renal tissue hypoxia has been a consistent finding in mouse and rat models of diabetic nephropathy [4], [30]. The reason why this occurs is less certain, but oxydative stress and glomerular hyperfiltration leading to increased tubular sodium load and oxygen consumption have been proposed as the main mechanisms [31], [32]. BOLD-MRI studies in humans have shown mixed results, some studies reporting higher, others lower R2* values in patients with diabetes [8], [9], [33]. In this study, DM2 patients had slightly higher cortical R2* levels than controls, suggesting lower tissue oxygenation, yet this association was not maintained in adjusted models. Interestingly, cortical oxygenation correlated linearly and positively with glucose levels as measured just before the MRI exam, suggesting that glycemia directly or indirectly (via sodium transport) influences cortical oxygenation. In a study assessing the influence of renin-angiotensin blockers on renal oxygenation in patients with type 2 DM, we have recently reported similar results [19]. These data suggest that episodes of hyperglycemia might induce transient increases in oxygen consumption due to hyperfiltration which increases sodium transport in the proximal tubule [34]. The pre-MRI glycemia has not been systematically reported in previous studies, which might partly explain previous mixed results. The herein described association offers an alternative explanation in humans for to the relationship between badly controlled diabetes mellitus and adverse renal outcome [35]. However, the described relationship is based on a cross-sectional analysis which limits causal inferences, and needs confirmation in interventional studies.
Our observation that circulating uric acid levels correlate positively with cortical R2* levels adds information to the ongoing debate about the role of uric acid in cardiovascular disease. Hyperuricemia-induced cortical vasoconstriction could be partly responsible for the herein described correlation between uric acid and cortical hypoxia [36], but more studies are necessary to clarify this.
This study has some limitations. Firstly, no information was acquired on renal perfusion. We therefore cannot determine whether changes in oxygenation were the result of alterations in oxygen delivery. Nevertheless, Textor and colleagues have shown previously that renal tissue oxygenation is largely independent of renal blood flow, and that cortical (but not medullary) oxygenation only falls in case of severe renal artery stenosis occluding more than 90% of the lumen [13]. Another limitation of our study is its cross-sectional design. Finally, all measurements were performed without stopping concomitant medication, which might have weakened eventual underlying differences in R2* levels between the groups
Conclusion
Our data suggest that renal tissue oxygenation at rest is comparable in controls, treated hypertensives and CKD patients. However, the response to furosemide differs between groups providing some insights on the mechanisms linking renal tubular function, oxygen consumption and renal function in hypertension and/or chronic kidney diseases. Our data do not support the concept that chronic kidney disease is associated with decreased renal tissue oxygenation in humans, as repeated observed in animal models. However, acute changes in oxygenation most likely occur in humans, as illustrated by the changes in R2* in response to furosemide. Future interventional studies should clarify the role of renal sodium handling, blood glucose and serum uric acid in the regulation of renal tissue oxygenation as newly described in this study. The differences in cortical R2* levels between men and women suggest that renal oxygenation is possibly regulated differently in men and women; this point also needs further study. Although no correlation was found between R2* values and (previous) kidney function, it is still possible that renal tissue R2* values predict kidney function decline but this is to be demonstrated in prospective studies.
Supporting Information
Figure S1.
Probability density function (Kernel density) of medullary R2* values by gender, according to group.
https://doi.org/10.1371/journal.pone.0095895.s001
(TIF)
Figure S2.
Medullary and cortical R2* values over decreasing eGFRmdrd values, by gender. The number of subjects was n = 44 (for the eGFR>90 ml/min/1.73 m2 category), 41 (eGFR 60–90), 25 (eGFR 30–60), 15 (eGFR <30) in men and respectively 30,24,11, and 5 in women.
https://doi.org/10.1371/journal.pone.0095895.s002
(TIF)
Table S1.
Drug treatment of patients across different stages of chronic kidney disease (CKD).
https://doi.org/10.1371/journal.pone.0095895.s003
(DOCX)
Table S2.
Multivariate linear regression analysis examining associations between baseline characteristics and medullary R2* levels. Associations between medullary R2* levels and baseline characteristics are expressed as regression coefficient β (95% CI).
https://doi.org/10.1371/journal.pone.0095895.s004
(DOCX)
Table S3.
Multivariate linear regression analysis examining associations between baseline characteristics and cortical R2* levels. Associations between cortical R2* levels and baseline characteristics are expressed as regression coefficient β (95% CI).
https://doi.org/10.1371/journal.pone.0095895.s005
(DOCX)
Table S4.
Multivariate linear regression analysis examining correlations between baseline characteristics and medullary and cortical R2* levels. This analysis includes the additional variables glycemia, serum uric acid level and 24 h urinary proteinuria.
https://doi.org/10.1371/journal.pone.0095895.s006
(DOCX)
Methods and Results S1.
Extended methods and results. - Screening. -Methods: Justification hydration protocol. -Results: * Furosemide-induced change in R2* in patients already on diuretics. * Sensitivity analysis in subgroup excluding healthy women <40 years and male CKD patients >75 years.
https://doi.org/10.1371/journal.pone.0095895.s007
(DOCX)
Author Contributions
Conceived and designed the experiments: M. Pruijm BV LH MS MB. Performed the experiments: M. Pruijm M. Piskunowicz CZ MM IB. Analyzed the data: M. Pruijm M. Piskunowicz LH IB MM BV MB. Wrote the paper: M. Pruijm MB.
References
- 1. Manotham K, Tanaka T, Matsumoto M, Ohse T, Miyata T, et al. (2004) Evidence of tubular hypoxia in the early phase in the remnant kidney model. J Am Soc Nephrol 15: 1277–1288.
- 2.
Haase VH (2013) Mechanisms of Hypoxia Responses in Renal Tissue. J Am Soc Nephrol.
- 3. Higgins DF, Kimura K, Bernhardt WM, Shrimanker N, Akai Y, et al. (2007) Hypoxia promotes fibrogenesis in vivo via HIF-1 stimulation of epithelial-to-mesenchymal transition. J Clin Invest 117: 3810–3820.
- 4. dos Santos EA, Li LP, Ji L, Prasad PV (2007) Early changes with diabetes in renal medullary hemodynamics as evaluated by fiberoptic probes and BOLD magnetic resonance imaging. Invest Radiol 42: 157–162.
- 5. Welch WJ, Baumgartl H, Lubbers D, Wilcox CS (2001) Nephron pO2 and renal oxygen usage in the hypertensive rat kidney. Kidney Int 59: 230–237.
- 6. Prasad PV (2006) Evaluation of intra-renal oxygenation by BOLD MRI. Nephron Clin Pract 103: c58–c65.
- 7. Simon-Zoula SC, Hofmann L, Giger A, Vogt B, Vock P, et al. (2006) Non-invasive monitoring of renal oxygenation using BOLD-MRI: a reproducibility study. NMR Biomed 19: 84–89.
- 8. Inoue T, Kozawa E, Okada H, Inukai K, Watanabe S, et al. (2011) Noninvasive evaluation of kidney hypoxia and fibrosis using magnetic resonance imaging. J Am Soc Nephrol 22: 1429–1434.
- 9.
Yin WJ, Liu F, Li XM, Yang L, Zhao S, et al. (2011) Noninvasive evaluation of renal oxygenation in diabetic nephropathy by BOLD-MRI. Eur J Radiol.
- 10. Michaely HJ, Metzger L, Haneder S, Hansmann J, Schoenberg SO, et al. (2012) Renal BOLD-MRI does not reflect renal function in chronic kidney disease. Kidney Int 81: 684–689.
- 11. Prasad PV, Epstein FH (1999) Changes in renal medullary pO2 during water diuresis as evaluated by blood oxygenation level-dependent magnetic resonance imaging: effects of aging and cyclooxygenase inhibition. Kidney Int 55: 294–298.
- 12. Pruijm M, Hofmann L, Maillard M, Tremblay S, Glatz N, et al. (2010) Effect of sodium loading/depletion on renal oxygenation in young normotensive and hypertensive men. Hypertension 55: 1116–1122.
- 13. Gloviczki ML, Glockner JF, Crane JA, McKusick MA, Misra S, et al. (2011) Blood oxygen level-dependent magnetic resonance imaging identifies cortical hypoxia in severe renovascular disease. Hypertension 58: 1066–1072.
- 14. Rognant N, Guebre-Egziabher F, Bacchetta J, Janier M, Hiba B, et al. (2011) Evolution of renal oxygen content measured by BOLD MRI downstream a chronic renal artery stenosis. Nephrol Dial Transplant 26: 1205–1210.
- 15. Schachinger H, Klarhofer M, Linder L, Drewe J, Scheffler K (2006) Angiotensin II decreases the renal MRI blood oxygenation level-dependent signal. Hypertension 47: 1062–1066.
- 16. Textor SC, Gloviczki ML, Flessner MF, Calhoun DA, Glockner J, et al. (2012) Association of filtered sodium load with medullary volumes and medullary hypoxia in hypertensive African Americans as compared with whites. Am J Kidney Dis 59: 229–237.
- 17. Levey AS, Eckardt KU, Tsukamoto Y, Levin A, Coresh J, et al. (2005) Definition and classification of chronic kidney disease: a position statement from Kidney Disease: Improving Global Outcomes (KDIGO). Kidney Int 67: 2089–2100.
- 18. O'Brien E, Asmar R, Beilin L, Imai Y, Mallion JM, et al. (2003) European Society of Hypertension recommendations for conventional, ambulatory and home blood pressure measurement. J Hypertens 21: 821–848.
- 19.
Pruijm M, Hofmann L, Zanchi A, Maillard M, Forni V, et al. (2012) Blockade of the renin-angiotensin system and renal tissue oxygenation as measured with BOLD-MRI in patients with type 2 diabetes. Diabetes Res Clin Pract.
- 20. Hofmann L, Simon-Zoula S, Nowak A, Giger A, Vock P, et al. (2006) BOLD-MRI for the assessment of renal oxygenation in humans: acute effect of nephrotoxic xenobiotics. Kidney Int 70: 144–150.
- 21. Djamali A, Sadowski EA, Muehrer RJ, Reese S, Smavatkul C, et al. (2007) BOLD-MRI assessment of intrarenal oxygenation and oxidative stress in patients with chronic kidney allograft dysfunction. Am J Physiol Renal Physiol 292: F513–F522.
- 22. Pedersen M, Dissing TH, Morkenborg J, Stodkilde-Jorgensen H, Hansen LH, et al. (2005) Validation of quantitative BOLD MRI measurements in kidney: application to unilateral ureteral obstruction. Kidney Int 67: 2305–2312.
- 23. Ebrahimi B, Gloviczki M, Woollard JR, Crane JA, Textor SC, et al. (2012) Compartmental analysis of renal BOLD MRI data: introduction and validation. Invest Radiol 47: 175–182.
- 24. Kojima I, Tanaka T, Inagi R, Kato H, Yamashita T, et al. (2007) Protective role of hypoxia-inducible factor-2alpha against ischemic damage and oxidative stress in the kidney. J Am Soc Nephrol 18: 1218–1226.
- 25. Fine LG, Norman JT (2008) Chronic hypoxia as a mechanism of progression of chronic kidney diseases: from hypothesis to novel therapeutics. Kidney Int 74: 867–872.
- 26. Priyadarshi A, Periyasamy S, Burke TJ, Britton SL, Malhotra D, et al. (2002) Effects of reduction of renal mass on renal oxygen tension and erythropoietin production in the rat. Kidney Int 61: 542–546.
- 27. Chun TY, Bankir L, Eckert GJ, Bichet DG, Saha C, et al. (2008) Ethnic differences in renal responses to furosemide. Hypertension 52: 241–248.
- 28. Burnier M, Bochud M, Maillard M (2006) Proximal tubular function and salt sensitivity. Curr Hypertens Rep 8: 8–15.
- 29. Zheleznova NN, Yang C, Ryan RP, Halligan BD, Liang M, et al. (2012) Mitochondrial proteomic analysis reveals deficiencies in oxygen utilization in medullary thick ascending limb of Henle in the Dahl salt-sensitive rat. Physiol Genomics 44: 829–842.
- 30. Prasad P, Li LP, Halter S, Cabray J, Ye M, et al. (2010) Evaluation of renal hypoxia in diabetic mice by BOLD MRI. Invest Radiol 45: 819–822.
- 31. Baines A, Ho P (2002) Glucose stimulates O2 consumption, NOS, and Na/H exchange in diabetic rat proximal tubules. Am J Physiol Renal Physiol 283: F286–293.
- 32. Hansell P, Welch WJ, Blantz RC, Palm F (2013) Determinants of kidney oxygen consumption and their relationship to tissue oxygen tension in diabetes and hypertension. Clin Exp Pharmacol Physiol 40: 123–137.
- 33. Wang ZJ, Kumar R, Banerjee S, Hsu CY (2011) Blood oxygen level-dependent (BOLD) MRI of diabetic nephropathy: preliminary experience. J Magn Reson Imaging 33: 655–660.
- 34. Vallon V, Schroth J, Satriano J, Blantz RC, Thomson SC, et al. (2009) Adenosine A(1) receptors determine glomerular hyperfiltration and the salt paradox in early streptozotocin diabetes mellitus. Nephron Physiol 111: 30–38.
- 35. Holman RR, Paul SK, Bethel MA, Matthews DR, Neil HA (2008) 10-year follow-up of intensive glucose control in type 2 diabetes. N Engl J Med 359: 1577–1589.
- 36. Messerli FH, Frohlich ED, Dreslinski GR, Suarez DH, Aristimuno GG (1980) Serum uric acid in essential hypertension: an indicator of renal vascular involvement. Ann Intern Med 93: 817–821.