Multicenter Study Evaluating Intrarenal Oxygenation and Fibrosis Using Magnetic Resonance Imaging in Individuals with Advanced CKD
Pottumarthi V. Prasad,Wei Li,Dominic S. Raj,James Carr,Molly C. Carr,Jon Thacker,Luping Li,Chi Wang,Stuart M. Sprague,Joachim H. Ix,Michel Chonchol,Geoffrey A. Block,Alfred K. Cheung,Kalani L. Raphael,Jennifer Gassman,Myles Wolf,Linda F. Fried,Tamara Isakova
DOI: https://doi.org/10.1016/j.ekir.2018.07.006
IF: 6.234
2018-01-01
Kidney International Reports
Abstract:The chronic hypoxia theory states that hypoxia and interstitial fibrosis are key contributors to progression of chronic kidney disease (CKD).1Fine L.G. Orphanides C. Norman J.T. Progressive renal disease: the chronic hypoxia hypothesis.Kidney Int Suppl. 1998; 65: S74-S78PubMed Google Scholar Presence of fibrosis may further enhance the hypoxia by limiting oxygen transport, resulting in a perpetual cycle of hypoxic injury and progressive loss of kidney function. Blood oxygenation level dependent (BOLD)2Prasad P.V. Edelman R.R. Epstein F.H. Noninvasive evaluation of intrarenal oxygenation with BOLD MRI.Circulation. 1996; 94: 3271-3275Crossref PubMed Scopus (327) Google Scholar and diffusion3Mao W. Zhou J. Zeng M. et al.Intravoxel incoherent motion diffusion-weighted imaging for the assessment of renal fibrosis of chronic kidney disease: a preliminary study.Magn Reson Imaging. 2018; 47: 118-124Crossref PubMed Scopus (27) Google Scholar magnetic resonance imaging (MRI) can provide information on renal oxygenation and fibrosis, respectively. The methods rely on endogenous contrast mechanisms that do not require exogenous contrast administration and that are widely available on major vendor platforms. We report baseline MRI data in 127 individuals with advanced CKD (mean estimated glomerular filtration rate [eGFR] = 33.4 ± 7.2 ml/min per 1.73m2) who participated in the COMBINE (CKD Optimal Management with BInders and NicotinamidE) study, a multicenter clinical trial that aimed to test whether nicotinamide and lanthanum carbonate would safely lower serum phosphate and FGF23 levels compared with placebo. Relaxation rate R2* served as the BOLD MRI index; higher values represent decreased oxygenation.2Prasad P.V. Edelman R.R. Epstein F.H. Noninvasive evaluation of intrarenal oxygenation with BOLD MRI.Circulation. 1996; 94: 3271-3275Crossref PubMed Scopus (327) Google Scholar Apparent diffusion coefficient (ADC) was the diffusion MRI index; lower values may be due to increased fibrosis.4Togao O. Doi S. Kuro-o M. et al.Assessment of renal fibrosis with diffusion-weighted MR imaging: study with murine model of unilateral ureteral obstruction.Radiology. 2010; 255: 772-780Crossref PubMed Scopus (129) Google Scholar Similar MRI data obtained in a small group (n = 13) of healthy volunteers were used for comparison. Detailed MRI methods are provided as Supplementary Methods and specific MRI parameters are listed in Supplementary Table S1. Table 1 summarizes the baseline characteristics of the study population.Table 1Baseline characteristics of study populationVariableCKD, n = 127Control, n = 13Age, yr65 ± 1259 ± 9Male, n (%)81 (64)6 (46)Race, n (%) White73 (57)7 (54) African American39 (31)3 (23) Other15 (12)3 (23)Diabetes, n (%)64 (50)0 (0)eGFR, ml/min per 1.73 m233.4 ± 7.290.7 ± 11.9UACR, mg/g161 (20–584)n/aHGB, g/dl12.9 ± 1.7n/aeGFR, estimated glomerular filtration rate; HGB, serum hemoglobin; n/a, not available; UACR, urine albumin to creatinine ratio shown as median (interquartile range). Open table in a new tab eGFR, estimated glomerular filtration rate; HGB, serum hemoglobin; n/a, not available; UACR, urine albumin to creatinine ratio shown as median (interquartile range). Figure 1 shows examples of MRI data in both a healthy volunteer and an individual with CKD. Table 2 summarizes the differences in MRI measurements between the trial participants with CKD and the healthy control group. Consistent with current literature on studies from single sites,5Inoue T. Kozawa E. Okada H. et al.Noninvasive evaluation of kidney hypoxia and fibrosis using magnetic resonance imaging.J Am Soc Nephrol. 2011; 22: 1429-1434Crossref PubMed Scopus (226) Google Scholar, 6Prasad P.V. Thacker J. Li L.P. et al.Multi-parametric evaluation of chronic kidney disease by MRI: a preliminary cross-sectional study.PLoS One. 2015; 10: e0139661Crossref PubMed Scopus (36) Google Scholar, 7Pruijm M. Milani B. Pivin E. et al.Reduced cortical oxygenation predicts a progressive decline of renal function in patients with chronic kidney disease.Kidney Int. 2018; 93: 932-940Abstract Full Text Full Text PDF PubMed Scopus (58) Google Scholar cortical R2* was modestly higher in the group with CKD versus the control group, suggesting decreased cortical oxygenation. R2*_Medulla was lower in participants with CKD compared with healthy volunteers, suggesting increased medullary oxygenation. This seemingly contradictory observation has been reported before in preclinical studies8Priyadarshi A. Periyasamy S. Burke T.J. et al.Effects of reduction of renal mass on renal oxygen tension and erythropoietin production in the rat.Kidney Int. 2002; 61: 542-546Abstract Full Text Full Text PDF PubMed Scopus (36) Google Scholar using invasive measurements and in human studies using BOLD MRI.9Milani B. Ansaloni A. Sousa-Guimaraes S. et al.Reduction of cortical oxygenation in chronic kidney disease: evidence obtained with a new analysis method of blood oxygenation level-dependent magnetic resonance imaging.Nephrol Dial Transplant. 2017; 32: 2097-2105PubMed Google Scholar, 10Wang Z.J. Kumar R. Banerjee S. et al.Blood oxygen level-dependent (BOLD) MRI of diabetic nephropathy: preliminary experience.J Magn Reson Imaging. 2011; 33: 655-660Crossref PubMed Scopus (51) Google Scholar In a remnant kidney model, direct oxygen levels measured using microelectrodes at 6 to 8 weeks showed increased tissue oxygen levels.8Priyadarshi A. Periyasamy S. Burke T.J. et al.Effects of reduction of renal mass on renal oxygen tension and erythropoietin production in the rat.Kidney Int. 2002; 61: 542-546Abstract Full Text Full Text PDF PubMed Scopus (36) Google Scholar In the same model, at an early phase (4–7 days), tissue hypoxia was observed.11Manotham K. Tanaka T. Matsumoto M. et al.Evidence of tubular hypoxia in the early phase in the remnant kidney model.J Am Soc Nephrol. 2004; 15: 1277-1288Crossref PubMed Scopus (190) Google Scholar These observations suggest the possibility that with advanced kidney disease, the decrease in renal perfusion is surpassed by decreased oxygen consumption that results from reduced delivery of glomerular filtrate and reduced tubular sodium transport.12Neugarten J. Renal BOLD-MRI and assessment for renal hypoxia.Kidney Int. 2012; 81: 613-614Abstract Full Text Full Text PDF PubMed Scopus (46) Google ScholarTable 2Comparison of magnetic resonance imaging parameters between study groupsCKD/ControlnMeanSDPaBy Student 2-tailed t test.Mean difference (confidence interval)Oxygenation R2*_Cortex (s−1)Control1318.742.370.022−1.811 (−3.322 to −0.3)CKD12320.553.10 R2*_Medulla (s−1)Control1329.033.87<0.015.278 (2.890 to 7.665)CKD12123.753.22 R2*_Kidney (s−1)Control1322.152.250.380.611 (−0.810 to 2.032)CKD12321.542.76 R2*_(Kid-Cor)(s−1)Control133.411.89<0.012.453 (1.303 to 3.604)CKD1230.960.89 R2*_MC ratioControl131.570.28<0.010.403 (0.230 to 0.577)CKD1211.170.15Response to furosemide ΔR2*_Cortex (s−1)Control130.420.860.11−0.461 (−1.035 to 0.112)CKD540.881.02 ΔR2*_Medulla (s−1)Control136.283.460.0023.747 (1.578 to 5.915)CKD542.542.47 ΔR2*_Kidney (s−1)Control132.031.130.0140.927 (0.217 to 1.638)CKD541.110.87Fibrosis ADCControl131.670.08< 0.010.219 (0.165 to 0.273) (× 10–3 mm2/s)CKD1261.450.17ADC, apparent diffusion coefficient; CKD, chronic kidney disease; R2*_(Kid-Cor) = R2*_Kidney – R2*_Cortex; R2*_MC ratio = R2*_Medulla / R2*_Cortex.a By Student 2-tailed t test. Open table in a new tab ADC, apparent diffusion coefficient; CKD, chronic kidney disease; R2*_(Kid-Cor) = R2*_Kidney – R2*_Cortex; R2*_MC ratio = R2*_Medulla / R2*_Cortex. Because cortico-medullary contrast is reduced in CKD,13Lee V.S. Kaur M. Bokacheva L. et al.What causes diminished corticomedullary differentiation in renal insufficiency?.J Magn Reson Imaging. 2007; 25: 790-795Crossref PubMed Scopus (47) Google Scholar identification of medulla may be challenging, especially when using small regions of interest (ROIs). This difficulty may result in lower interreader agreement for medullary ROIs, as has been reported previously.14Thacker J.M. Li L.P. Li W. et al.Renal blood oxygenation level-dependent magnetic resonance imaging: a sensitive and objective analysis.Invest Radiol. 2015; 50: 821-827Crossref PubMed Scopus (22) Google Scholar We had a high interreader agreement in medullary ROIs (Supplementary Table S2), and we included whole-kidney ROI specifically to mitigate the limitation of small ROIs for evaluating medulla. The difference between kidney and cortex ROIs could be used as an indirect estimate of medullary R2* with a higher degree of objectivity. Consistent with R2*_Medulla, changes in R2*_(Kid-Cor) show a net positive mean difference value when compared with controls (Table 2), whereas R2*_Cortex had a negative mean difference. This supports that the observed increase in R2*_Medulla is not an artifact of using small ROIs. Prior studies have reported increased medullary oxygenation using both small ROIs10Wang Z.J. Kumar R. Banerjee S. et al.Blood oxygen level-dependent (BOLD) MRI of diabetic nephropathy: preliminary experience.J Magn Reson Imaging. 2011; 33: 655-660Crossref PubMed Scopus (51) Google Scholar and more objective concentric ROI method.9Milani B. Ansaloni A. Sousa-Guimaraes S. et al.Reduction of cortical oxygenation in chronic kidney disease: evidence obtained with a new analysis method of blood oxygenation level-dependent magnetic resonance imaging.Nephrol Dial Transplant. 2017; 32: 2097-2105PubMed Google Scholar The response to furosemide provides an index of active tubular sodium reabsorption in the medulla.15Gomez S.I. Warner L. Haas J.A. et al.Increased hypoxia and reduced renal tubular response to furosemide detected by BOLD magnetic resonance imaging in swine renovascular hypertension.Am J Physiol Renal Physiol. 2009; 297: F981-F986Crossref PubMed Scopus (41) Google Scholar Consistent with prior reports, we observed a blunted response to furosemide in participants with CKD compared with controls.6Prasad P.V. Thacker J. Li L.P. et al.Multi-parametric evaluation of chronic kidney disease by MRI: a preliminary cross-sectional study.PLoS One. 2015; 10: e0139661Crossref PubMed Scopus (36) Google Scholar, 16Pruijm M. Hofmann L. Piskunowicz M. et al.Determinants of renal tissue oxygenation as measured with BOLD-MRI in chronic kidney disease and hypertension in humans.PLoS One. 2014; 9e95895Crossref PubMed Scopus (60) Google Scholar Notably, for the analysis of response to furosemide, we excluded participants with CKD who were chronically treated with loop diuretics, given a prior study reporting lower response in such individuals.17Hall M.E. Rocco M.V. Morgan T.M. et al.Chronic diuretic therapy attenuates renal BOLD magnetic resonance response to an acute furosemide stimulus.J Cardiovasc Magn Reson. 2014; 16: 17Crossref PubMed Scopus (14) Google Scholar Lower ADC values are associated with the presence of fibrosis.3Mao W. Zhou J. Zeng M. et al.Intravoxel incoherent motion diffusion-weighted imaging for the assessment of renal fibrosis of chronic kidney disease: a preliminary study.Magn Reson Imaging. 2018; 47: 118-124Crossref PubMed Scopus (27) Google Scholar, 4Togao O. Doi S. Kuro-o M. et al.Assessment of renal fibrosis with diffusion-weighted MR imaging: study with murine model of unilateral ureteral obstruction.Radiology. 2010; 255: 772-780Crossref PubMed Scopus (129) Google Scholar The cortical ADC estimates from this study in the healthy volunteers are comparable with a recent report,3Mao W. Zhou J. Zeng M. et al.Intravoxel incoherent motion diffusion-weighted imaging for the assessment of renal fibrosis of chronic kidney disease: a preliminary study.Magn Reson Imaging. 2018; 47: 118-124Crossref PubMed Scopus (27) Google Scholar and as also previously reported, we observed lower cortical ADC in participants with CKD compared with controls (Table 2). The cortical ADC values in our participants with CKD were lower than previously reported3Mao W. Zhou J. Zeng M. et al.Intravoxel incoherent motion diffusion-weighted imaging for the assessment of renal fibrosis of chronic kidney disease: a preliminary study.Magn Reson Imaging. 2018; 47: 118-124Crossref PubMed Scopus (27) Google Scholar (1.45 ± 0.17 vs. 1.63 ± 0.14 × 10–3mm2/s), probably due to the lower mean eGFR in our study in comparison with the mean eGFR in the prior study (33.4 vs. 71.2 ml/min per 1.73 m2). Interestingly, 50% of the group with CKD had diabetes, and participants with CKD with diabetes had significantly lower ADC compared with those without diabetes (Table 3).Table 3Comparison of magnetic resonance imaging measures by diabetes statusCKD/ControlnMeanSDPaP < 0.05 by Student 2-tailed t test.Mean difference (confidence interval)Oxygenation R2*_Cortex (s−1)Diabetes6120.553.000.990.007 (−1.106 to 1.118)No diabetes6220.543.23 R2*_Medulla (s−1)Diabetes6024.133.320.210.745 (−0.412 to 1.901)No diabetes6123.383.10 R2*_Kidney (s−1)Diabetes6121.612.620.780.143 (−0.846 to 1.132)No diabetes6221.462.92 R2*_MC ratioDiabetes601.190.170.110.043 (−0.010 to 0.097)No diabetes611.150.13Response to furosemide ΔR2*_Cortex (s−1)Diabetes170.781.1780.6340.156 (−0.510 to 0.822)No diabetes370.930.958 ΔR2*_Medulla (s−1)Diabetes171.761.470.0541.137 (−0.022 to 2.296)No diabetes372.902.76 ΔR2*_Kidney (s−1)Diabetes170.960.930.4390.208 (−0.335 to 0.751)No diabetes371.170.84Fibrosis ADCDiabetes631.40aP < 0.05 by Student 2-tailed t test.0.15< 0.020.097 (0.038 to 0.156) (× 10–3 mm2/s)No diabetes631.500.19ADC, apparent diffusion coefficient; CKD, chronic kidney disease.Bold values indicate P < 0.05.a P < 0.05 by Student 2-tailed t test. Open table in a new tab ADC, apparent diffusion coefficient; CKD, chronic kidney disease. Bold values indicate P < 0.05. To characterize MRI measurements across all clinical sites, we examined box-whisker plots for R2*_MC Ratio, ADC_Cortex, eGFR, and urine albumin to creatinine ratio (UACR) by individual sites (Figure 2). Analysis of variance showed no significant differences in any of these parameters between sites. Finally, we explored associations of MRI measurements with clinical parameters in participants with CKD (Table 4). R2*_Medulla and R2*_MC Ratio showed significant association with eGFR and UACR (Table 5). These relationships remained significant after adjusting for age, race, gender, and diabetes status. The association of medullary R2* with UACR may be interesting because higher UACR is considered to be a predictor of fast progression.18Lorenzo V. Saracho R. Zamora J. et al.Similar renal decline in diabetic and non-diabetic patients with comparable levels of albuminuria.Nephrol Dial Transplant. 2010; 25: 835-841Crossref PubMed Scopus (44) Google Scholar R2*_Medulla was associated with ADC and remained significant even after adjusting for eGFR and UACR (Table 6). ADC was also associated with diabetes status, and remained significant even after adjusting for eGFR and UACR (Table 7). Further studies are necessary to confirm these findings. Race and gender did not show differences in any of the MRI parameters (data not shown).Table 4Spearman correlations between magnetic resonance imaging parameters and renal function in the participants with CKDADCdiffusioneGFRLog(UACR)Oxygenation:ρ0.011−0.027−0.111 R2*_CortexP0.9030.7690.225n122123122ρ0.184aCorrelation is significant at the 0.05 level (2-tailed).0.195aCorrelation is significant at the 0.05 level (2-tailed).−0.311bCorrelation is significant at the 0.01 level (2-tailed). R2*_MedullaP0.0440.0320.001n120121120ρ0.0620.026−0.209aCorrelation is significant at the 0.05 level (2-tailed). R2*_KidneyP0.5000.7750.021n122123122ρ0.223aCorrelation is significant at the 0.05 level (2-tailed).0.150−0.225aCorrelation is significant at the 0.05 level (2-tailed). R2*_MC ratioP0.0140.1000.013n120121120Response to furosemide:ρ−0.0090.045−0.080 ΔR2*_CortexP0.9520.7480.571n535453ρ0.257−0.018−0.139 ΔR2*_MedullaP0.0630.8990.321n535453ρ0.0530.086−0.239 ΔR2*_KidneyP0.7050.5390.085n535453Fibrosis:ρ1.000−0.053−0.167 ADCP0.5570.063n126125125Conventional parameters:ρ−0.0531.000−0.105 eGFRP0.5570.244n125126125ρ−0.167−0.1051.000 Log(UACR)P0.0630.244n125125126ADC, apparent diffusion coefficient; CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; UACR, urinary albumin to creatinine ratio.Bold values indicate P < 0.05.a Correlation is significant at the 0.05 level (2-tailed).b Correlation is significant at the 0.01 level (2-tailed). Open table in a new tab Table 5Linear regression of MRI indices with UACR and eGFRDependent variablePredictorsβSEPAdjustedaAdjusted for age, race, gender, and diabetes status.βSEPR2*_MedullaUACRbLog transformed.−1.1250.3420.001−0.7980.3640.021eGFR0.0810.0400.0480.0800.0390.044R2*_MC ratioUACRbLog transformed.−0.0390.0160.017−0.0350.0170.048eGFR0.0040.0020.0490.0040.0020.038R2*_KidneyUACRbLog transformed.−0.5690.3000.061−0.3530.3160.267eGFR0.0160.0350.6490.0130.0340.703eGFR, estimated glomerular filtration rate; MRI, magnetic resonance imaging; UACR, urine albumin to creatinine ratio.Bold values indicate P < 0.05.a Adjusted for age, race, gender, and diabetes status.b Log transformed. Open table in a new tab Table 6Linear regression of ADC with R2* indicesDependent variablePredictorsβSEPAdjustedaAdjusted for estimated glomerular filtration rate and log-transformed urine albumin to creatinine ratio.βSEPADCR2*_Medulla0.0100.0050.0360.0110.0050.043R2*_MC ratio0.1570.1070.1440.1560.1120.165ADC, apparent diffusion coefficient.Bold values indicate P < 0.05.a Adjusted for estimated glomerular filtration rate and log-transformed urine albumin to creatinine ratio. Open table in a new tab Table 7Linear regression of ADC with diabetes status (0 = no; 1 = yes)Dependent variablePredictorsβSEPAdjustedaAdjusted for estimated glomerular filtration rate and log-transformed urine albumin to creatinine ratio.βSEPADCDiabetes status−0.0970.0300.002−0.0960.0310.003ADC, apparent diffusion coefficient.Bold values indicate P < 0.05.a Adjusted for estimated glomerular filtration rate and log-transformed urine albumin to creatinine ratio. Open table in a new tab ADC, apparent diffusion coefficient; CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; UACR, urinary albumin to creatinine ratio. Bold values indicate P < 0.05. eGFR, estimated glomerular filtration rate; MRI, magnetic resonance imaging; UACR, urine albumin to creatinine ratio. Bold values indicate P < 0.05. ADC, apparent diffusion coefficient. Bold values indicate P < 0.05. ADC, apparent diffusion coefficient. Bold values indicate P < 0.05. Our study was performed in patients with a 20 < GFR < 45 ml/min, so conclusions based on these results cannot be generalized to patients with all stages of CKD. Sodium intake was not controlled, which has been shown to affect renal medullary oxygenation.19Pruijm M. Hofmann L. Maillard M. et al.Effect of sodium loading/depletion on renal oxygenation in young normotensive and hypertensive men.Hypertension. 2010; 55: 1116-1122Crossref PubMed Scopus (53) Google Scholar Use of hand-drawn ROIs for the medulla may not be objective.14Thacker J.M. Li L.P. Li W. et al.Renal blood oxygenation level-dependent magnetic resonance imaging: a sensitive and objective analysis.Invest Radiol. 2015; 50: 821-827Crossref PubMed Scopus (22) Google Scholar Ideally, future studies should use fully automated segmentation of the kidneys performed on high-contrast anatomic images that are coregistered to the BOLD MRIs.20Cox E.F. Buchanan C.E. Bradley C.R. et al.Multiparametric renal magnetic resonance imaging: validation, interventions, and alterations in chronic kidney disease.Front Physiol. 2017; 8: 696Crossref PubMed Scopus (49) Google Scholar It is not clear whether stopping angiotensin-converting enzyme inhibitors/angiotensin receptor blockers for 1 day before MRI and oral loop diuretics only on the day of the MRI is sufficient in this group of participants with advanced CKD. Use of the same dose of furosemide in individuals with lower renal function may not be optimal. Participants in this study did not undergo kidney biopsy, so we cannot directly determine correlations of the MRI findings with biopsy-proven fibrosis measures. The 6 centers enrolled participants from varying clinic settings. The control group had a limited number (n = 13) and were all from a single center. All participating sites used MRI scanners from a single vendor. In conclusion, our data support the feasibility of using renal BOLD and diffusion MRI in multicenter trials and the data are consistent with prior reports based on single-site studies. These data combined with other results from a recent report21Li L.P. Thacker J. Li W. et al.Consistency of multiple renal functional MRI measurements over 18 months.J Magn Reson Imaging. 2018; 48: 514-521Crossref PubMed Scopus (3) Google Scholar support planned international initiatives, such as BEAt-DKD,22Biomarker Enterprise to Attack Diabetic Kidney Disease (BEAt-DKD). Available at: www.beat-dkd.eu. Accessed on August 6, 2018.Google Scholar that involve longitudinal multicenter trials using multiparametric MRI. Overall, our observations in advanced CKD further confirm the reduced renal cortical oxygenation and presence of renal fibrosis consistent with the chronic hypoxia hypothesis. The medullary oxygenation was significantly increased compared with controls and is also consistent with prior reports. The significantly lower ADC in participants with diabetes is novel and may have clinical relevance. Future studies to monitor progressive changes in eGFR are needed to verify if and which of these MRI parameters are specific to progressive CKD. In participants with advanced CKD, evaluating response to furosemide may be limited by the low baseline medullary R2* and their potential chronic use of loop diuretics. All the authors declared no competing interests. The COMBINE study is supported by grants U01DK099877 , U01DK097093 , U01DK099930 , U01DK099933 , U01DK099924 , and R01DK102438 (TI) from the National Institute of Diabetes and Digestive and Kidney Diseases . The work also was supported in part by R01DK093793 (PVP) and F31HL123360 (JT). Download .docx (.03 MB) Help with docx files Supplementary Methods Download .pdf (.01 MB) Help with pdf files Table S1MRI acquisition parameters. Download .pdf (.02 MB) Help with pdf files Table S2Intra- and interreader agreement.