Abstract
Objective
To investigate changes in renal oxygenation levels by blood-oxygenation-level dependent (BOLD)-magnetic resonance imaging (MRI), and to evaluate BOLD-MRI for detecting early diabetic renal injury.
Methods
Seventy-five rats, with unilateral nephrectomy, were randomly divided into streptozotocin-induced diabetes mellitus (DM, n = 65) and normal control (NC, n = 10) groups. BOLD-MRI scans were performed at baseline (both groups) and at 3, 7, 14, 21, 28, 35, 42, 49, 56, 63 and 70 days (DM only). Renal cortical (C) and medullary (M) R2* signals were measured and R2* medulla/cortex ratio (MCR) was calculated.
Results
DM-group CR2* and MR2* values were significantly higher than NC values following diabetes induction. R2* values increased gradually and peaked at day 35 (CR2*, 33.95 ± 0.34 s–1; MR2*, 43.79 ± 1.46 s–1), then dropped gradually (CR2*, 33.17 ± 0.69 s–1; MR2*, 41.61 ± 0.95 s–1 at day 70). DM-group MCR rose gradually from 1.12 to 1.32 at day 42, then decreased to 1.25 by day 70.
Conclusions
BOLD-MRI can be used to non-invasively evaluate renal hypoxia and early diabetic renal injury in diabetic rats. MCR may be adopted to reflect dynamic changes in renal hypoxia.
Introduction
Diabetic nephropathy is a microvascular complication of diabetes mellitus and the leading cause of end-stage renal failure. 1 Renal hypoxia caused by microvascular disease, is closely related to continuous renal injury and eventually leads to diabetic nephropathy. 2 Functional (f) magnetic resonance imaging (MRI) is a rapidly developing imaging technology that enables the evaluation of organ function and pathology.3,4 Blood-oxygenation-level dependent (BOLD)-MRI, a type of fMRI, can be used to noninvasively monitor oxygen content in tissue, and to elucidate metabolic status.5,6 This technique is performed with the aid of an endogenous contrast agent, namely paramagnetic deoxyhaemoglobin, interfering with magnetic spin dephasing of blood water protons, which results in a BOLD enhancement effect. 7 It has been shown that the renal R2* value is positively related to the level of deoxyhaemoglobin in tissue, while the latter is correlated with oxygen partial pressure.5,8
The application of BOLD-MRI in diabetic nephropathy has been studied in animal models and in patients with diabetes.9–11 Use of BOLD-MRI with an animal model of diabetic nephropathy revealed that diabetic renal oxygen content was significantly decreased compared with controls. 9 In addition, a decrease in oxyhaemoglobin was detected as early as 2 days following administration of streptozotocin (STZ) in a rat model of diabetes, and was suggestive of renal hypoxia. 10 The above animal studies generally focused on selected time-points during early diabetic renal injury, and lacked dynamic and continuous observation, while clinical research has mainly concentrated on patients at later stages, as early stages are difficult to diagnose due to lack of symptoms.11,12 In addition, results from clinical studies may be influenced by certain drugs (such as antihypertensive drugs or diuretics), and limited availability of renal biopsy tissue remains a problem.
Streptozotocin can accelerate pathological kidney changes in rats following unilateral nephrectomy,13,14 and is a preferred choice for creating an animal model of diabetes, as it selectively damages pancreatic β cells, leading to increased blood glucose, while having little effect on other organs. 15
Based on the above findings, the aim of the present study was to use BOLD-MRI to track changes in renal oxygenation over time, in an STZ-induced rat model of diabetes compared with normal control rats.
Materials and methods
Animal model and study design
A total of 80 healthy male Sprague-Dawley rats (aged 4–6 weeks; weight, 150–170 g) were provided by the experimental animal centre of the First Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, China. All procedures were approved by the ethics committee on the use of animals in experimentation, namely, the Animal Experimental Ethical Inspection committee of the First Affiliated Hospital, School of Medicine, Zhejiang University, and were performed in accordance with the standards of the Chinese Association for Laboratory Animal Sciences (http://english.calas.org.cn/). The rats were provided with free access to water and normal diet, and housed in a well-ventilated environment maintained at 18–25°C, with a relative humidity of 40–70% and 12 h of light daily.
As the left kidney is easily affected by bowel gas, all rats in the present study underwent left kidney resection, followed by one week of enhanced care (daily or more regular bedding change, increased food and water changes). A total of 80 rats then received the normal care regime for one week, and those without weight loss during the week were further tested for blood glucose, serum creatinine and urea nitrogen levels via tail vein samples. A total of 75 rats with normal blood glucose and renal function ranges (weight, 200–230g) were randomly divided into a nephrectomised diabetes mellitus group (DM, n = 65) and normal control group (NC, n = 10). The rat model of diabetes (DM group) was established by intraperitoneal injection of 65 mg/kg STZ 16 (Sigma-Aldrich, St Louis, MO, USA) dissolved in 0.1 mol/l citric acid-sodium citrate buffer (pH 4.6), following a 12 h fast, and the NC group was administered an equivalent volume of 0.1 mol/l citric acid-sodium citrate buffer (pH 4.6) by intraperitoneal injection. Blood glucose from tail vein samples was measured using Contour® TS blood glucose test strips (Bayer Health Care, Leverkusen, Germany) at 3 days following injection. Fasting blood glucose levels ≥ 16.7 mmol/l were considered to indicate successful modelling. The upper limit for the test strip assay is 33.3 mmol/l, so blood glucose measurements above 33.3 mmol/l were recorded as ≥ 33.3 mmol/l.
Rats allocated to the DM group were divided into five subgroups, giving five DM groups and one NC group. During the study, one group of rats was randomly selected for MRI scanning and for 24 h urinary albumin excretion at each-time point (3, 7, 14, 21, 28, 35, 42, 49, 56, 63 and 70 days) following induction of diabetes. At day 7, 21, 35, 56 and 70, respectively, blood creatinine and urea nitrogen were also tested from tail vein samples, following the BOLD-MRI scan. These rats were then euthanized by intraperitoneal injection of 150 mg/kg pentobarbital, and the kidney was excised for histopathological examination.
MRI and data analysis
Magnetic resonance imaging measurements were acquired using a Signa HDxt 3.0T MRI scanner (GE Healthcare, Chicago, IL, USA). The rats were anaesthetized by intraperitoneal injection with 0.007 ml/g of 4% chloral hydrate (supplied by the pharmaceutical preparation centre of the First Affiliated Hospital, School of Medicine, Zhejiang University). During imaging, animals were placed in a supine and head advanced position within an HD T/R Knee Array 8-channel phase array coil (Invivo, Gainesville, FL, USA). To reduce image artefacts caused by breathing movements, limbs were fixed by adhesive tape and the abdomen was wrapped in gauze.
The BOLD MRI T2* weighted images were acquired in coronal plane using a fast-spoiled gradient recovery sequence. MRI parameters included: repetition time, 100 ms; echo time (TE) series, 1.7 ms, 3.8 ms, 5.9 ms, 8 ms, 10.1 ms, 12.2 ms, 14.3 ms, and 16.4 ms; echo numbers, 8; flip angle, 30°; receiver bandwidth, 31.5 kHz; field of view, 60 mm × 60 mm; slice thickness, 2.0 mm; intersection gap, 0.2 mm; acquisition matrix, 64 × 64; and number of excitations, 2. The T2* and R2* images were obtained using Advantage workstation 4.4 (GE Healthcare), which supports FuncTool software (GE Healthcare) for automated images post-processing, and automatically generates curves for BOLD signal intensities against different TEs (Figure 1). R2* medulla to cortex ratio (MCR) was also calculated.

Representative images obtained using blood-oxygenation-level dependent (BOLD)- magnetic resonance imaging (MRI) with Advantage workstation 4.4 and FuncTool software to assess the right kidney in a rat model of diabetes: (a) T2* weighted image used to show structures of the cortex and medulla and to provide a reference for measuring R2* values; (b) and (c) R2* and T2* map, respectively; and (d) time signal curve of region of interest (ROI) in the cortex and medulla
All Bold measurements were obtained from the processed R2* images using the region of interest (ROI) method. To avoiding the influence of artefacts, ROIs should be placed as far as possible from each other and avoiding vessels visible to the naked eye. Besides artefacts, the setting of ROIs must also be considered, because a slight variation in placement, size or shape of an ROI in the rat kidney may lead to different results. Thus, the placement of ROIs on an R2* image was determined by two experienced radiologists (QW and CG) reaching a consensus through discussion. As a result, five ROIs each in the cortex and outer medulla were chosen, respectively, with each ROI comprising 15–30 pixels. The R2* values (mean ± SDs) for cortex and medulla were then calculated. The outer medulla in the renal medullary region was chosen for measurement because Na+-K+-ATPases are mainly distributed in this area, which consumes almost half the oxygen of the whole kidney,17–19 and renal injury has been confirmed by histopathology to occur more easily in this region.9,19 In addition, the inner medulla is influenced by the structure of the renal sinus in R2* maps, which may cause inaccurate measurements.
24-h urinary albumin excretion and blood biochemical examinations
Before BOLD-MRI, rats from the randomly chosen DM group were placed in metabolic cages for 24 h to collect 24-h urine samples. The urine albumin concentration was determined using a Beckman Coulter IMMAGE® 800 Immunochemistry system and associated reagents (Beckman Coulter, Brea, CA, USA) according to the manufacturer’s instructions, and then multiplied by 24 h urine volume to get 24 h urinary albumin excretion. Blood creatinine and urea nitrogen were tested from tail vein blood samples using a Beckman Coulter AU5800 Clinical Chemistry System (Beckman Coulter) according to the manufacturer’s instructions, and this group was euthanized for histopathological examination, as previously stated.
Histopathological analysis
Rat kidneys were perfused via the abdominal aorta with 1.25 % phosphate buffered glutaraldehyde, pH 7.4, at a constant pressure (100 mmHg). The fixed rat kidney fragments were then cut into 3–4 μm sections and stained with haematoxylin and eosin, and periodic acid-Schiff. The glomerulosclerosis index, stereology of mesangial volume density and capillary surface density were estimated at day 7, 21, 35, 56 and 70, respectively. The glomerulosclerosis index was obtained by measuring the frequency of sclerosis in 100 glomeruli from each animal, as reported previously, 20 and the mean ± SD was calculated. Scores ranged from 0 to 4, depending on the extent of sclerosis: 0, normal glomerulus; 1+, ≤25% sclerosis; 2+, > 25% to ≤50% sclerosis; 3+, >50% to ≤75% sclerosis; and 4+, >75% glomerulosclerosis. The glomerulosclerosis index score for every animal was the sum of individual glomerular scores multiplied by the percent of glomeruli with that same score. Stereology of mesangial volume density and capillary surface density were estimated using a test grid (M42) as described previously. 21 Electron micrographs of each glomeruli were analysed at a magnification of x 10 500. The capillary surface density (Sv) was estimated by counting the number of intersections between the capillary surface and the test line (IL), where Sv: = 2 × IL. The reference volume was determined by counting the test points that hit the glomeruli (PT). The volume density of this structure was estimated by counting the number of points that hit the mesangium (PP), where Vv: = PP/PT.
Statistical analyses
Data are presented as mean ± SD. Statistical analyses were performed using IBM SPSS Statistics package, version 19.0 (IBM, Chicago, IL, USA). Comparison of R2* values and histopathological results at each time-point between the DM and NC groups were performed using Student’s t-test. A P value <0.05 was considered statistically significant.
Results
Experimental model characteristics
A total of 65 rats were successfully induced by STZ in the DM group, however, three rats from the DM group died during the study, one died of anaesthesia and two died of diabetic ketoacidosis. Thus, 62 rats in the DM group (12 or 13 rats in each of the five subgroups), and 10 rats in the NC group were included in the final analyses (Table 1). Rats in the DM group displayed polydipsia, polyphagia, polyuria, dispirited spirit, weight loss and poor fur colour. Some rats appeared with genital infection and gradually developed cataract after 6 weeks. Body weights of rats in the DM group were lower than in the NC group after day 3 (P < 0.001), and continued to decrease gradually. Kidney size (evaluated by measuring the largest diameter point of the right kidney) increased gradually following day 7, but was not significantly different from the NC group (P > 0.05). There were no statistically significant differences in blood sugar levels between the control and DM groups at baseline, however, blood sugar levels ≥ 33.3 mmol/l were reached in the DM group by day 21 (Table 1).
Physiological changes in a streptozotocin-induced rat model of diabetes (n = 62) and in normal rat controls (n = 10)
Data presented as mean ± SD.
NC, normal controls; DM, diabetes mellitus group.
aDM subgroups were randomly selected for testing at each time-point.
*P < 0.01 or **P < 0.001, compared with controls (Student’s t-test).
Urinary albumin excretion, blood biochemistry and pathology findings
Urinary albumin excretion in the DM group increased gradually, and was significantly higher than that of the NC group from day 7 (DM, 1.20 ± 0.08 mg/24 h; NC, 1.11 ± 0.17 mg/24 h; P < 0.01; Table 2). Serum creatinine and blood urea nitrogen values increased in the DM group over time (serum creatinine levels increased from 12 μmol/l at day 28 to 38 μmol/l at day 70, and blood urea nitrogen levels increased from 7.7 mmol/l at day 28 to 16.6 mmol/l at day 70). Both parameters were higher in the DM than NC group (creatinine range, 10–11 μmol/l and blood urea nitrogen range, 6–7 mmol/l).
Urinary albumin excretion over time in a streptozotocin-induced rat model of diabetes (n = 62) and normal controls (n = 10)
Data presented as mean ± SD.
NC, normal controls; DM, diabetes mellitus group.
aDM subgroups were randomly selected for testing at each time-point.
*P < 0.01 or **P < 0.001, compared with controls (Student’s t-test).
The glomerulosclerosis index was higher in the DM group than the NC group from day 56 (P < 0.01). Mesangial volume density and capillary surface density was increased in the DM group compared with the NC group (P < 0.01 on day 70, and P < 0.01 from day 35, respectively; Table 3).
Histopathological examination of the right kidney in a streptozotocin-induced rat model of diabetes (n = 62) and normal controls (n = 10)
Data presented as mean ± SD.
NC, normal controls; DM, diabetes mellitus group.
*P < 0.01 compared with controls (Student’s t-test).
BOLD-MRI results
Renal R2* values in the DM group rats were significantly higher than in the NC group rats at all time-points following induction of diabetes (P < 0.01; Table 4 and Table 5). CR2* and MR2* values in the DM group were 28.30 ± 0.59 s−1 and 31.87 ± 0.74 s−1 at baseline, then increased gradually to peak at day 35 (CR2*, 33.95 ± 0.34 s−1, MR2*, 43.79 ± 1.46 s−1). Following day 35, values gradually decreased to 33.17 ± 0.69 s−1 (CR2*) and 41.61 ± 0.95 s−1 (MR2*) at day 70 (Table 4 and Table 5, and Figure 2). Baseline MCR in the DM group was 1.12, which rose to 1.32 at day 42, then dropped down to 1.25 at day 70 (Figure 2).
Changes in CR2* values over time, measured using blood oxygenation level dependent-magnetic resonance imaging, in a streptozotocin-induced rat model of diabetes (n = 62) and normal controls (n = 10)
Data presented as mean ± SD.
NC, normal controls; DM, diabetes mellitus group.
aDM subgroups were randomly selected for testing at each time-point.
**P < 0.001 compared with controls (Student’s t-test).
Changes in MR2* values over time, measured using blood oxygenation level dependent-magnetic resonance imaging, in a streptozotocin-induced rat model of diabetes (n = 62) and normal controls (n = 10)
Data presented as mean ± SD.
NC, normal controls; DM, diabetes mellitus group.
aDM subgroups were randomly selected for testing at each time-point.
**P < 0.001 compared with controls (Student’s t-test).

Box-whisker plots showing: (a) Changes in cortex (C)R2* values in a streptozotocin-induced rat model of diabetes mellitus (DM group) over 70 days (s–1); (b) Changes in medulla (M)R2* values in the DM group over 70 days (s–1); and (c) Changes in the medulla to cortex ratio (MCR) in the DM group over 70 days; the heavy black horizontal lines are medians, the extremities of the box are the 25th and 75th percentiles, the error bars represent minimum and maximum outliers, the circles above the day 14 and 21 bars represent extreme outliers
Discussion
According to published diabetic nephropathy diagnostic criteria, 22 data from the present study regarding urinary albumin excretion and pathological changes indicate that the present diabetic rats were in nephropathy stage I within 7 days, nephropathy stage II within 7–35 days, and nephropathy stage III within days 35–70 days. The present BOLD-MRI results showed that CR2* and MR2* values within the experimental period in the DM group were increased compared with controls, suggesting decreased oxygenation in both the renal cortex and medulla at different time-points in the presence of diabetes, consistent with previous reports.9–11,23 Hypoxic changes in diabetic kidneys may be due to the hyperfiltration phenomenon seen early after induction of diabetes, whereby the high filtration and high metabolic rates activate renal tubular reabsorption and Na+/K+-ATPase activity, resulting in increased oxygen consumption.9,24,25 A previous study documented that renal tubules in rats with diabetes consumed > 40% more oxygen than those in non-diabetic rats. 9 In addition, an evaluation of the influence of STZ-induced diabetes on renal pO2 and blood flow, using invasive microprobes and BOLD-MRI, 10 suggested that hypoxic changes are not dependent on blood flow until 28 days after induction of diabetes, and these changes can be detected as early as 2 days in rat kidneys by BOLD-MRI. The present results showed that microalbuminuria differed between normal controls and a rat model of diabetes from day 7 following induction of diabetes, but cortex and medulla R2* values were significantly higher than controls from day 3, and showed an increasing trend, in the DM group. These results may indicate that BOLD-MRI can detect renal oxygen changes earlier and more sensitively than microalbuminuria.
The present study continued for up to 70 days, and the remaining diabetic rats developed hyperglycaemia, eye cataracts and raised urinary microproteins, suggesting progressive kidney damage. Hypoxia peaked on day 35, and after this time, kidney R2* values in diabetic rats decreased. The present results showed a rising then falling trend in CR2* and MR2* values during diabetic renal injury. That is, when the glomerulosclerosis index and capillary surface density increased significantly from day 35, the diabetic rats in the present study were at nephropathy stage II to stage III. Interestingly, decreasing CR2* and MR2* values suggested that renal hypoxia status was improved during this period, however, this phenomenon was not found in similar published reports, and BOLD-MRI cannot differentiate between changes oxygen supply versus demand.
The mechanism of hypoxia change cannot be explained clearly by the present study, however, an invasive investigation of renal function in diabetic rats in nephropathy stages II and III, 14 demonstrated that glomerular filtration rate (GFR) was always much higher in rats with diabetes than controls. Also, renal plasma flow was gradually increased and renal vascular resistance reduced accordingly, from stages II to stages III. Thus, the kidney oxygen supply should be increased because of increasing blood flow. This may explain the descending part of the dynamic R2* curve in the present results, which may have been due to haemodynamic changes leading to increased blood supply in the renal cortex and medulla, and subsequently, improved renal oxygenation.
In support of the present results, a similar BOLD-MRI R2* pattern has been shown in patients with type 2 diabetes, 11 and hypoxia was documented to be predominantly caused by hyper metabolism during the early clinical stage of diabetic nephropathy, with oxygen consumption decreasing gradually in later stages, as a result of lower workload. That is, the deterioration of tubular function and lower GFR probably raised pO2 in the medullary region, although MR2* in diabetic nephropathy was markedly higher compared with controls. 11 The present results are also supported by other studies that investigated rats with renal artery stenosis.26,27 For example, after 4 weeks, no renal hypoxia was detected by BOLD-MRI in the rat kidney downstream of a renal artery stenosis, 26 and in another study, remnant kidneys had 73% increased oxygen tensions (measured using needle electrodes) in both the cortex and outer medulla for 6–8 weeks following ligation of branches of the renal artery, 27 suggested to be due to higher oxygen tensions caused by a decrease in filtration fraction, or due to a possible adaptive mechanism. In addition to the above explanations, other factors may be involved in reduced oxygen consumption, such as renal tubule cells maintaining concentration and dilution function of the urine, and elevated levels of oxidative stress and inflammatory responses of the renal tissues causing dysfunction of renal tubule cells.22,28,29 Combinations of all these factors have been reported to be associated with changes in BOLD-MRI signals compared with controls.
The CR2* and MR2* value can be influenced by magnetic field inhomogeneity, coil position and the physiological condition of subject,8,29,30 thus, the medulla to cortex R2* ratio (MCR) was also calculated in the present study, to decrease the effect of these factors. In the present study, dynamic changes in the diabetic MCR resembled a downward parabola. This complicated MCR reverse change is similar to that seen in patients with diabetic nephropathy from stage III into stage IV, 11 and is important for patients, as it means the disease deteriorates from early to later stages. Mesangial proliferation is significant and nodular sclerosis appears pathologically at this point, however, patients take a long time to develop from stage I to stage IV and are more difficult to diagnose in the early stages, due to absence of symptoms. 31 Patients may, therefore, lose treatment opportunities to improve renal oxygenation. BOLD-MRI can detect oxygenation changes earlier than pathological changes, so the reversion of MCR in the present results may indicate that the kidney hypoxia was a result of increasing oxygen consumption during the early clinical stage of diabetic nephropathy leading to haemodynamic and architectural modification. Thus, the reversion of MCR may indicate kidney change from functional compensatory changes to organic lesions. Observation of MCR may, therefore, have great significance for physicians in prevention and control of diabetic renal injury.
The results of the present study may be limited by several factors. First, the sample number in the experimental group was relatively small at the later stages, due to the limited number of remaining rats. Secondly, the change in renal oxygenation levels in diabetic rats treated with insulin, particularly in stage II diabetic rats, was not investigated, but would be the present authors’ future direction of research. Thirdly, changes in serum glucose and glycosylated haemoglobin (HbA1c) levels were not measured at time-points during which R2* values were measured. HbAlc is the standard measure of glycaemic control, and the HbAlc value can be used to estimate average glucose. 32 The impact of HbAlc on the R2* value has been recognised, and the most obvious effect of hyperglycaemia in the human body is on brain tissue. Acute changes in glycaemia can lead to hemiparesis and confusion. 33 A study to determine the effect of plasma glucose on BOLD-MRI R2* values in six healthy humans, 34 found that severe acute hyperglycaemia did not appear to substantially effect BOLD-MRI values. In addition, there are currently no published reports regarding the effects of hyperglycaemia and HbA1C on R2* values in the kidneys of diabetic rats. Hyperglycaemia also affects volume and tissue water content in the kidney. In the present study, the rats had free access to water and food, and differences in renal R2* values between diabetic rats and normal controls has been suggested to be caused by water diuresis (however, no statistical significance was shown). 9 The present study did not control for diuresis, so the effect may not be excluded.
In conclusion, the dynamic changes in R2* values in a rat model of diabetes showed a rising then falling trend in early diabetic renal injury. At first, R2* values appear to reflect renal hypoxia, and are more sensitive in detecting early renal changes compared with microalbuminuria. Then, renal oxygenation appears to be elevated and renal hypoxia is improved. The reversal of MCR seen in the present study may indicate renal injury from increased oxygen consumption to increased blood supply, and might be a reversible transition stage for diabetes, and thus, may be expected to become a sensitive indicator for evaluating the degree of renal injury.
Footnotes
Declaration of conflicting interests
The authors declare that there is no conflict of interest.
Funding
This study was supported by grants from the Natural Science Foundation of Zhejiang Province (Grant No. 2015C33138), and from the Science and Technology Department of Public Welfare Projects of Zhejiang Province (Grant No. LZ16H180001).
