- 1Department of Morphology, Imaging, Orthopedics, Rehabilitation and Nutrition, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
- 2Department of Medical Imaging, Ghent University Hospital, Ghent, Belgium
- 3Ghent Institute for Functional and Metabolic Imaging, Ghent University, Ghent, Belgium
- 4Department of Diagnostic Sciences, Faculty of Medicine and Health Sciences, Ghent University, Ghent, Belgium
- 5Department of Veterinary and Biosciences, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
- 6Institute of Biomedical Engineering and Technology, Faculty of Engineering and Architecture, Ghent University, Ghent, Belgium
Introduction: Dynamic contrast-enhanced (DCE) MRI and arterial spin labeling (ASL) MRI enable non-invasive measurement of renal blood flow (RBF), whereas blood oxygenation level-dependent (BOLD) MRI enables non-invasive measurement of the apparent relaxation rate (R2*), an indicator of oxygenation. This study was conducted to evaluate the potential role of these MRI modalities in assessing RBF and oxygenation in dogs. The correlation between contrast-enhanced ultrasound (CEUS) and the MRI modalities was examined and also the ability of the MRI modalities to detect pharmacologically induced changes.
Methods: RBF, using CEUS, ASL- and DCE-MRI, as well as renal oxygenation, using BOLD-MRI of eight adult beagles were assessed at two time-points, 2–3 weeks apart. During each time point, the anesthetized dogs received either a control (0.9% sodium chloride) or a dopamine treatment. For each time point, measurements were carried out over 2 days. An MRI scan at 3 T was performed on day one, followed by CEUS on day two.
Results: Using the model-free model with caudal placement of the arterial input function (AIF) region of interest (ROI) in the aorta, the DCE results showed a significant correlation with ASL measured RBF and detected significant changes in blood flow during dopamine infusion. Additionally, R2* negatively correlated with ASL measured RBF at the cortex and medulla, as well as with medullary wash-in rate (WiR) and peak intensity (PI). ASL measured RBF, in its turn, showed a positive correlation with cortical WiR, PI, area under the curve (AUC) and fall time (FT), and with medullary WiR and PI, but a negative correlation with medullary rise time (RT). During dopamine infusion, BOLD-MRI observed a significant decrease in R2* at the medulla and entire kidney, while ASL-MRI demonstrated a significant increase in RBF at the cortex, medulla and the entire kidney.
Conclusion: ASL- and BOLD-MRI can measure pharmacologically induced changes in renal blood flow and renal oxygenation in dogs and might allow detection of changes that cannot be observed with CEUS. However, further research is needed to confirm the potential of ASL- and BOLD-MRI in dogs and to clarify which analysis method is most suitable for DCE-MRI in dogs.
1 Introduction
Acute kidney injury (AKI) and chronic kidney disease (CKD) are common conditions that compromise renal function and affect both humans and dogs (1–4). As the glomerular filtration rate (GFR) is directly proportional to functional renal mass, it is regarded as a highly sensitive and specific indicator of renal function impairment (5). The GFR can be assessed in a number of ways using both direct and indirect measurement techniques. Direct measurement of GFR based on the clearance of a filtration marker is considered the gold standard, but this method is cumbersome and can only assess global renal function (5). Instead, GFR is often measured indirectly using surrogate markers (e.g., blood urea nitrogen and creatinine) in routine veterinary practice because of their ease of use (5). A major drawback of these commonly employed surrogate markers, however, is their inability to identify kidney injury at an early stage (5). The imaging technique scintigraphy, on the other hand, makes it possible to determine, in addition to the total GFR, the individual contribution of both kidneys within a short time frame of 20 to 30 min (6, 7). Renal scintigraphy in dogs relies on the detection of gamma radiation from the injected radiopharmaceutical technetium-99 m-diethylenetriaminepetaacetic acid (99mTc-DTPA), which the kidney removes from the blood via glomerular filtration (6, 7). The GFR is calculated using a regression equation relating the percentage dose uptake of 99mTc-DTPA by each kidney to the GFR measured with plasma inulin clearance (6, 7). In dogs, scintigraphy has been shown to be equivalent to inulin clearance in determining GFR (8, 9). Nevertheless, scintigraphy has some disadvantages as well, such as the need for radioactive isotopes, the need for short-term patient isolation, and the fact that it is only available in veterinary specialist clinics with advanced technological capabilities (7).
In addition to GFR, renal perfusion and oxygenation are essential factors to consider in the renal function assessment (10, 11). The development and progression of various forms of AKI and CKD are significantly influenced by perfusion impairment and hypoxia of the renal tissue (10, 11). In the early stages of renal disease, hypoperfusion and hypoxia may even be present without other signs of renal dysfunction (12, 13). Assessment of tissue perfusion and oxygenation is complicated by the lack of gold standard techniques and limitations of current reference techniques (14, 15). Furthermore, this complicates the evaluation of potential new techniques to assess tissue perfusion and oxygenation (14, 15).
Contrast-enhanced ultrasound (CEUS) is an imaging technique that enables assessment of microvascular perfusion by intravenous administration of an ultrasound contrast agent (16). Ultrasound contrast agents consist of tiny, encapsulated gas bubbles that are contained within the intravascular space and reflect ultrasound waves as harmonic signals (16). Several studies have demonstrated the potential of CEUS in evaluating kidney diseases in humans and animals. Among others, CEUS was found to be a sensitive method to assess the changes in renal microvascular perfusion in canine acute kidney disease (17), as well as in chronic kidney disease in cats (18) and humans (19). A drawback of CEUS, however, is its high intrinsic variability (20).
Functional magnetic resonance imaging (fMRI), on the other hand, is an emerging imaging technique that can non-invasively assess both renal structure and function of the entire renal parenchyma. Using fMRI modalities such as dynamic contrast-enhanced (DCE) MRI, arterial spin labeling (ASL) and blood oxygenation level-dependent (BOLD) MRI, multiple functional kidney parameters can be measured in a single scan session, including perfusion, glomerular filtration and oxygenation (21).
With DCE MRI, a gadolinium contrast agent is injected intravenously and the contrast uptake is analyzed based on a pharmacokinetic (PK) model to determine RBF and GFR of each kidney (21, 22). Although there have been some concerns regarding the safety of gadolinium contrast, especially in patients with renal impairment, low doses used in clinical settings have been demonstrated to be not harmful (23). Among the advantages of DCE-MRI are its lack of radioactive radiation, its limited duration, the ability to determine single kidney GFR, and the possibility to use other MRI modalities simultaneously (23). DCE-MRI has been successfully used to measure renal function parameters (GFR and RBF) in humans (24) and animals, such as pigs (25), rabbits (26, 27), rats (28, 29), and mice (30, 31). So far, DCE-MRI has only been used in a limited number of studies on dogs. Aumann et al. (32) evaluated the feasibility of renal perfusion assessment with DCE-MRI in dogs by comparison with renal artery blood flow measurements using an implanted flow probe. Another study demonstrated that quantitative renal perfusion measurements of DCE-MRI can be used for the evaluation of renal artery stenosis in dogs (33). One canine study has described the successful quantification of the single kidney GFR using DCE-MRI in dogs (34).
ASL-MRI enables tissue perfusion quantification by using magnetically labeled blood water protons as an endogenous contrast agent (15). A labeled image is first collected in which a radiofrequency pulse alters the longitudinal magnetization of arterial blood flowing into the kidney. An image with an unchanged magnetization is obtained next as a control. By subtracting the label from the control image, a perfusion-weighted image is created where signal intensity is proportional to perfusion (15). The perfusion weighted images are analyzed by a kinetic model to quantify renal perfusion and get a quantitative perfusion map (15). The feasibility of RBF determination by ASL-MRI of kidneys with both normal and altered function was demonstrated by previous studies in humans (35, 36), pigs (37, 38), rabbits (27, 39) and rats (40, 41). There was a good correlation between ASL-MRI measurements of renal blood flow and measurements from alternative methods to evaluate renal blood flow. For example, ASL-MRI showed a good correlation with DCE-MRI (42), plasma clearance of para-aminohippuric acid (35) and technetium-99 m-mercaptoacetyltriglycine (99mTc-MAG3) scintigraphy (43) in humans, with DCE-MRI in rabbits (27), and with microspheres and ultrasound flowmeter in pigs (37, 38). To our knowledge, no other researchers have conducted any studies on dogs with renal ASL-MRI.
In BOLD-MRI the paramagnetic effect of deoxyhemoglobin is used as an endogenous contrast agent (44). Tissue deoxyhemoglobin concentration rises as tissue oxygen content decreases, causing the effective transverse relaxation time (T2*) of adjacent water protons to decrease (44). The apparent relaxation rate R2* (1/T2*), the reciprocal of T2*, rises with increasing concentrations of deoxygenated hemoglobin and is a measure of tissue oxygenation (44). The sensitivity of BOLD-MRI to assess intrarenal oxygenation under physiological and various pathophysiological conditions has been demonstrated by several studies in humans (45–47) and animal models (pig (48, 49), rabbit (39, 50), rat (41, 51) and mice (52)). However, some studies have produced contradictory findings, showing no distinction between healthy and diseased kidneys (53, 54). In dogs, there has only been one prior study that used BOLD-MRI to quantify renal oxygenation to our knowledge (55). Lee et al. (55) studied the feasibility and repeatability of BOLD-MRI derived T2* between two sessions in eight healthy beagles under anesthesia before and after furosemide administration.
A comprehensive assessment of renal function can be obtained using a combination of several fMRI techniques in addition to morphological images, i.e., multiparametric MRI (56, 57). Multiparametric MRI has the potential to elucidate (patho) physiological processes of the kidney and evaluate treatment effects. A great deal remains to be discovered about the normal physiology of renal oxygenation and perfusion, as well as the causes and consequences of their dysregulation in renal disease (56, 57). Due to the invasive nature of the present experimental and transitional research methods, client-owned dogs and human volunteers cannot be enrolled in certain experiments (48). The main aim of this study is to evaluate fMRI’s potential value for the assessment of renal perfusion and oxygenation in dogs. The objective of this study was to investigate which DCE analysis method, based on RBF measurements using ASL-MRI, is most appropriate for evaluating RBF in dogs. This study also examined the correlation between measurements of CEUS and the fMRI modalities (ASL-, BOLD- and DCE-MRI). Furthermore, the ability of fMRI modalities to detect pharmacologically induced changes in renal parameters, perfusion and oxygenation, was evaluated.
2 Materials and methods
2.1 Dogs
This study used eight purpose-bred beagles, four males (two intact, two neutered) and four females (three intact, one spayed). Their age ranged from 4 to 9.6 years [mean ± standard deviation (SD), 5.4 ± 1.6 years] and their body weight ranged from 9.4 to 15.6 kg (mean ± SD, 12.5 ± 1.9 kg). A standard dry maintenance diet (Medium mature diet, Royal Canin®) was provided to the dogs as well as ad libitum access to tap water. A physical examination, thoracic radiography, abdominal ultrasonography, Doppler blood pressure measurement, and laboratory analyses including complete blood count, serum biochemistry profile and urinalysis (urine sediment, dipstick test, specific gravity (USG), urine protein:creatinine (UPC) ratio and bacterial culture) were used to establish the baseline health status. The study protocol was approved by the Animal Ethics Committee from the Faculty of Veterinary Medicine and the Faculty of Bioscience Engineering of Ghent University, Belgium (Approval number: EC2021-085).
2.2 Protocol
The renal perfusion and renal oxygenation of the dogs were assessed at two time-points, 2–3 weeks apart, using MRI and CEUS. During each time point, the dogs received either a control treatment with 0.9% sodium chloride (NaCl) or a treatment with dopamine. A different treatment was given to each dog at each time point, and the order of the treatments was randomized. For each time point, measurements were carried out over 2 days. An fMRI scan including ASL-, BOLD- and DCE-MRI was performed on day one, followed by CEUS on day two. The study protocol is illustrated in Figure 1.
Figure 1. Flow chart of the study. Parts of the figure were drawn using modified pictures from Server Medical Art. Servier Medical Art by Servier is licensed under Creative Commons Attribution 4.0 International license (https://creativecommons.org/licenses/by/4.0/).
2.3 Anesthesia
Food and water were withheld for 12 and 2 h before anesthesia induction. In the cephalic vein, a 22-gauge IV catheter was placed and butorphanol (0.2 mg/kg) (Dolorex®; MSD Animal Health, The Netherlands) was injected intravenously to induce sedation. Propofol (4–6 mg/kg IV) (PropoVet®; Zoetis, BE) was administered in conjunction with midazolam (0.2 mg/kg IV) (Midazolam Accord Healthcare®; Accord Healthcare Limited, UK) to induce anesthesia. Anesthesia was maintained with 1.2–1.4% isoflurane (Isoflutek®; Laboratorios Karizoo, Spain) in 100% oxygen, delivered via a circle rebreathing system, after dogs were intubated. Anesthesia monitoring included measurement of inspiratory and end-expiratory isoflurane and CO2 levels by capnography, invasive and/or non-invasive blood pressure measurement, and blood oxygen saturation and heart rate measurement by pulse oximetry.
Dogs received dopamine (3 μg/kg/min) (Dopaminehydrochloride Hikma® 40 mg/mL; Hikma Pharmaceuticals, UK) or sodium chloride (0.9%) via an infusion pump at a rate of 3 mL/kg/h. The infusion started 5 min after the start of inhalation anesthesia and the evaluation of renal perfusion with the imaging techniques started 15 min after the start of the intravenous infusion of dopamine or 0.9% sodium chloride.
2.4 Contrast enhanced ultrasound
All scans were performed by a board-certified radiologist (E.S.) and an ECVDI-trained radiologist (E.V.), with 10 and 6 years of experience with CEUS, respectively. An ultrasound machine (iU22, Philips, Bothell, WA, United States) equipped with contrast specific software and a 12–5 MHz linear transducer was used for the CEUS examinations. The protocol used for the CEUS examinations was based on previous studies and was standardized (58, 59). The mechanical index was set to 0.08, in order to minimize microbubble destruction. Other standardized machine settings included focus (underneath the kidney), persistency (disabled) and dynamic range (50). The depth was adjusted to have the kidney in the field-of-view and the gain was optimized during the study, starting with a nearly anechoic image in contrast mode. The dogs were positioned in dorsal recumbency and coupling gel was applied to the area of shaved skin. The kidney’s longitudinal plane was imaged in dual screen-mode (side by side display of B- and contrast-mode images) and the transducer position was manually maintained throughout the examination. The 22-G catheter in the cephalic vein was fitted with a three-way stopcock. Using the three-way stopcock, sulfur hexafluoride-filled microbubbles (SonoVue®, Bracco Diagnostics Inc., Milan, Italy) were injected intravenously as a bolus (0.04 mL/kg) followed by 2 mL saline solution (NaCl 0.9%). Concurrent with injection of the contrast agent, a 90-s recording started. First, the left kidney was scanned, then the right kidney. Two injections were performed for the left kidney. Since the first injection typically results in lower enhancement (60), only images from the second injection were used for analysis. For the right kidney, one injection was administered. If technical difficulties were encountered during recording and images were therefore unsuitable for analysis, additional injections were performed. In all dogs, the injections were administered by the same person (A.H.) in a consistent way. Between two injections, the residual circulating microbubbles were eliminated by using high mechanical index pulses while scanning the caudal abdominal aorta and spleen. This was performed until a similar level of background echogenicity was obtained as before injection.
Quantitative analysis of the recordings was performed with dedicated software (QLAB, Philips Healthcare). Six regions of interest (ROIs) were manually positioned: 3 in the cortex, 2 in the medulla and 1 surrounding the entire kidney. Cortical and medullary ROIs were rectangular and had a size of 10 (± 0.20) mm2 and 6 (± 0.20) mm2, respectively (Figures 2A,B). Placement of the cortical and medullary ROIs was done at the same organ depth in a region with homogeneous enhancement, avoiding the inclusion of arteries. A time-intensity curve (TIC) was generated automatically for each ROI by the software. An average was calculated for the 3 ROIs in the cortex and 2 ROIs in the medulla. From the raw data, the following perfusion variables were calculated: peak intensity (PI), wash-in rate (WiR), wash-out rate (WoR), rise time (RT), fall time (FT), time to peak (TTP), total area under the curve (AUC) (Figure 2C). Blood volume is reflected by intensity-related parameters like PI, AUC, and WiR, whereas blood velocity is reflected by time-related parameters like TPP, RT, and FT (16).
Figure 2. Illustration of regions-of-interest (ROIs) placement in the kidney on a contrast-enhanced ultrasound image (A) and a B-mode ultrasound image (B), as well as of blood flow parameters in a time intensity curve (C). (A,B) Three ROIs were drawn in the cortex (white, solid line) and two ROIs were drawn in the medulla (white, striped line). (C) Peak intensity (PI), total area under the curve (AUC), wash in rate (WiR), wash-out rate (WoR), time to peak (TTP), rise time (RT) and fall time (FT).
The PI was defined as the maximum signal intensity minus the pre-rise baseline intensity. The WiR refers to the slope of the linear regression line for the data from 10% above baseline intensity to 85% of peak intensity. WoR refers to the slope of the linear regression line for data after the peak intensity, i.e., from 85 to 40% of peak intensity. TI is the time point at which the ascending linear regression line (10–85%) intersects the x-axis. TO is the time point at which the descending linear regression line (85–40%) intersects the x-axis. RT was defined as the time from the onset of enhancement to peak enhancement. FT was defined as the time from peak enhancement to TO. TTP was defined as the time between contrast injection and peak enhancement. AUC refers to the area under the time-intensity curve from timepoints TI to TO.
2.5 Magnetic resonance imaging scan protocol
All MRI examinations were performed with a 3 T MRI scanner (Siemens 3 T Magnetom Prisma Fit, Siemens AG, Healthcare Sector, Erlangen, Germany) running VE11C, using an 18 Channel Body Coil (Body 18, Siemens AG, Healthcare Sector, Erlangen, Germany). The dogs were positioned in dorsal recumbency with their head facing the MRI control room. A localizer sequence was acquired to determine the kidneys’ location and support scan planning.
2.5.1 Arterial spin labeling
ASL images were acquired by using a Flow-sensitive Alternating Inversion Recovery (FAIR) labeling scheme with Q2TIPS (Siemens Work In Progress sequence). The readout parameters of the FAIR Q2TIPS ASL sequence were: imaging matrix 64 × 64, field of view (FOV) 272 × 136 mm2, voxel size 2.1 × 2.1 × 8.0 mm3, TI1 2000 ms, bolus length 1,000 ms, repetition time (TR)/ echo time (TE) 4,500/23.58 ms, flip angle 180°, bandwidth 4,340 Hz/Px, slice thickness 8 mm, number of slices 8. The total scan time of the ASL sequence was 4 min 33 s. Background suppression pulse was used to suppress T1 values of 230 ms. At each time point two consecutive ASL scans were performed.
2.5.2 Blood oxygen level-dependent
A multi-echo gradient-recalled echo (mGRE) sequence with 12 equally spaced echo times was used to obtain BOLD images. Consensus-based technical recommendations for clinical translation of renal BOLD MRI were followed [19]. The acquisition parameters of the BOLD sequence included the following: TR 70 ms, TEs 4.20–50.40 ms (echo time intervals 4.20 ms), matrix 96 × 96, FOV 180×180, number of slices 4, slice thickness 4 mm, coronal oblique plane, flip angle 30°, voxel size 1.9×1.9×4.0 mm, bandwidth 300 Hz/Px, fat saturation was applied. The parallel imaging technique GRAPPA (generalized autocalibrating partially parallel acquisitions) was used with an acceleration factor of 2. With an acquisition time of 13 s, each kidney was scanned once. Each scan was completed during a 15 s breath-hold, which was induced by applying end-expiratory positive pressure with the breathing balloon.
2.5.3 Dynamic contrast-enhanced
For DCE MRI acquisition the RAVE (Radial Volumetric Encoding) sequence was used with the following parameters: TR/TE of 4.80/1.37 ms, flip angle of 17°, 36 slices, slice thickness of 3 mm, matrix of 256×128 mm, FOV of 380×380 mm, voxel size of 1.5×1.5×3.0 mm, bandwidth of 780 Hz/Px. A dose of 0.2 mL/kg of Gadolinium contrast agent (Clariscan®, 0.5 mmoL/mL, GE Healthcare, Nydalen, Oslo) was injected as a bolus at 1 mL/s followed by a 12 mL saline flush at the same rate using an automated contrast injector (BAYER Medrad® Spectris Solaris EP).
2.6 Post-processing
2.6.1 Arterial spin labeling
Label, control, and M0 images, as well as Perfusion-Weighted (PW) and RBF maps were generated inline on the scanner, preceded by motion correction that was applied to the acquired MR datasets using an elastic registration algorithm (Applications Guide, Arterial Spin Labeling Sequence for Kidney Perfusion Assessment, Siemens Healthineers, Erlangen DE). Using MRIcroGL (RRID:SCR_024413), all ASL DICOM images were converted to NIfTI format (61). The mid-coronal slice of each kidney was identified, and a whole kidney region of interest (ROI) was manually drawn by a fourth-year PhD student (A.H.) on the PW images with 3D Slicer (RRID:SCR_005619). RBF was determined by ASL-MRI using the average of two consecutive scans and is expressed in mL/100 g/min.
2.6.2 Blood oxygen level-dependent
Four-dimensional images were produced by the scanner, with the first three dimensions representing the acquired volume and the fourth dimension representing the voxel value at different echo times. A voxel-by-voxel calculation was performed to generate estimated R2* maps following Zhang et al. (62). On a mid-coronal slice of the T2* weighted image, which served as an anatomic template, a mask covering the renal parenchyma of the whole kidney was manually defined by a fourth-year PhD student (A.H.) using 3D Slicer (RRID:SCR_005619).
2.6.3 Regional image analysis
Analysis of ASL and BOLD images was conducted using the 12-layer concentric objects method (TLCO) and TEO method (Ten Equiangular Object) (45, 46, 63). With the TLCO-method, an automatic algorithm divides the renal parenchyma selected by the whole kidney ROI into 12 concentric layers of equal thickness (45, 46). The novel approach called the TEO method divides the kidney into ten equiangular segments from the cranial to the caudal pole, enabling evaluation along the renal cortex (63). Figures 3A,B illustrates the division of the kidney into twelve concentric layers using the TLCO method and in ten equiangular segments using the TEO method.
Figure 3. Illustration of the regional division of the kidney by the TLCO (A) and TEO method (B), as well as the RBF (C,D) and R2* plots (E,F) of all subjects using the TLCO and TEO method. (A) The Twelve Layer Concentric Objects (TLCO) technique, which divides the renal parenchyma into 12 layers of uniform thickness. (B) The Ten Equiangular Objects (TEO) technique, which divides the kidney from the cranial to the caudal pole into ten equiangular segments. The RBF plots (mean ± SD) of all subjects during 0.9% NaCl infusion (control) and dopamine infusion using the TLCO (C) and TEO (D) method. The R2* plots (mean ± SD) of all subjects during 0.9% NaCl infusion (control) and dopamine infusion using the TLCO (E) and TEO (F) method. In the TLCO method, layer 1 refers to the outer layer and layer 12 to the inner layer. In the TEO method, the segment at the caudal pole is represented by 1 and the segment at the cranial pole by 10. *p-value <0.05.
The TLCO-method involves calculating a center of mass C(x,y) for each mask (64). Each point P(x,y) within the mask is classified into 12 equal sections based on its distance from C. By implementing the TLCO method in Python [Python Programming Language (RRID:SCR_008394)], a scientific software development environment, the average perfusion and oxygenation values at different distances from the center were calculated. As a result, a label image with 12 concentric layers was obtained as well as the average RBF and R2* values for each layer. Values for the cortex, medulla, and whole kidney have been estimated based on the measurements of these 12 layers, using the following definitions: cortex (mean of the three outer layers), medulla (mean of the inner 8th to 10th layer) and whole kidney (mean of all 12 layers) (45). The TEO method, implemented in Python [Python Programming Language (RRID:SCR_008394)], defined 10 regions. A label image with 10 equiangular objects was obtained. The mean perfusion rates and oxygenation levels were calculated using TEO across the different kidney sections.
2.6.4 Dynamic contrast-enhanced
All processing was carried out offline. Image data from DCE was transferred in DICOM format and reconstructed with Golden-angle Radial Sparse Parallel imaging (GRASP) (65). The reconstructed images were analyzed with PMI 0.4, written in IDL 6.4 (RRID:SCR_025084).
The arterial input function (AIF) in this study describes the changes in contrast agent concentration over time in the abdominal aorta supplying blood to the kidneys. Four different placements for the ROI in the aorta to define the AIF were evaluated in this study (Figure 4). The AIF ROIs were rectangular and had a standardized size (10 × 4.4 mm). The first and second ROI were placed 12 and 25 mm cranial to the right renal artery branch, respectively. The third ROI was also located cranial to the right renal artery, but just caudal to the celiac and cranial mesenteric arteries to avoid this area of increased intensity. The fourth ROI was placed 12 mm caudal to the right renal artery branch. A renal ROI was drawn around the left and right kidney, respectively. The renal ROIs were defined on the middle coronal slice of the perfusion image of the respective kidney. For RBF determination following compartment models were evaluated: C2 filtration model (C2FM), C2 uptake model (C2UM), model-free (MF) and maximum slope (MS).
Figure 4. Illustration of the MRI region of interest (ROI) settings for the kidneys and arterial input function (AIF). Placement of the renal ROI on the right (A) and left (D) kidney. Placement of the AIF ROI: 12 mm cranial to the right renal artery branch (B), 25 mm cranial to the right renal artery branch (C), just caudal to the celiac and cranial mesenteric arteries (E) and 12 mm caudal to the right renal artery branch (F).
2.6.5 Kidney volume
For kidney volume determination, manual segmentation on each T1 VIBE (Volume Interpolated Breath-hold Examination) image slice was performed using 3D Slicer (RRID:SCR_005619). The software then generated a renal volume automatically. The sequence parameters of the T1 VIBE ISO sequence were: FOV 380 mm, voxel size: 1.0×1.0×1.0 mm, slice thickness: 1.0 mm, TR/TE: 3.64/ 1.45 ms, flip angle: 9.0°, number of slices: 144, bandwidth 450 Hz/Px.
2.7 Statistics
Statistical analyses were done in R version 4.3.2 (R Project for Statistical Computing (RRID:SCR_001905)). The analysis was done stepwise. First, correlations between techniques were calculated in the native phase, i.e., when the dogs did not receive treatment, using the spearman correlation coefficient. Based on the guidelines from Schober et al. (66), correlation was considered negligible if r was less than 0.10, weak if r was between 0.10 and 0.39, moderate if r was between 0.40 and 0.69, strong if r was between 0.70 and 0.89 and very strong if r was greater than 0.90. Next, the effect of receiving dopamine on various outcome parameters (ASL TLCO, ASL TEO, BOLD TLCO, BOLD TEO, DCE, CEUS, respectively) was evaluated using a linear mixed model with treatment as fixed effect and the appropriate random effect. For CEUS, the random effect was dog; for DCE, the random effect was dog and side nested within dog; for BOLD and ASL, the random effect was dog, side nested within dog and segment nested within side. In every model, significance was assessed with a likelihood ratio test. Finally, the association between kidney volume and treatment was evaluated with a linear mixed model with treatment as fixed effect and as random effects dog, and side nested within dog. Significance was set at α ≤ 0.05 and a Bonferroni-correction for multiple testing was applied when needed. Corrected p-values are reported.
3 Results
3.1 Pharmacokinetic model for dynamic contrast-enhanced-MRI
The uptake model and the model-free model were the only DCE-MRI PK models whose RBF measurements demonstrated a significant correlation with RBF measurements of ASL-MRI at baseline. With regard to the model-free model, this applied to analyses where the AIF ROI was located 25 mm cranially, 12 mm cranially, or caudally to the right renal artery. For the uptake model, this specifically concerned the analysis in which the AIF ROI was positioned 25 mm cranially to the right renal artery. For the different combinations of PK model and AIF ROI placement in DCE analysis, Table 1 shows the significance of the correlation between the measured RBF and the RBF measured with ASL-MRI.
Table 1. Overview of the significance of the correlations between RBF measured with ASL-MRI and RBF measured with DCE-MRI for different combinations of PK model and AIF ROI placement.
3.2 Correlation between imaging methods at baseline
Between ASL- and BOLD-MRI measurements a significant negative correlation was found at the level of the cortex (r = −0.42, p < 0.01) and the medulla (r = −0.46, p = 0.001). So there is an inverse relation between renal perfusion and R2*, R2* decreases with an increasing RBF. This suggests that increased renal perfusion is accompanied by increased renal oxygenation.
Table 2 summarizes the correlations of RBF and R2* with the CEUS parameters at different levels of the kidney. A significant positive correlation was observed for blood volume parameters WiR and PI with RBF, both at the level of the cortex and medulla. The RBF likewise revealed a significant positive correlation with blood volume parameter AUC and blood velocity parameter FT, however, only at the cortex level. On the contrary, a negative correlation was found between blood velocity parameter RT and RBF at the medulla. Based on these findings, elevated renal blood volume and increased rate of blood inflow are related to higher renal blood flow.
Table 2. Correlation coefficients of RBF and R2* with CEUS parameters at the cortex and the medulla.
For R2*, significant negative correlations were found with PI and WiR at the medulla. At the cortex, however, no significant correlations were observed between R2* and any of the CEUS parameters. These results suggest that increased medullary renal blood volume is related to a higher oxygenation.
3.3 Pharmacological challenge with dopamine
The mean ± SD RBF values during 0.9% NaCl infusion (control) and dopamine infusion of each layer using the TLCO method and of each segment obtained using the TEO method is illustrated in Figures 3C,D. Dopamine infusion induced a significant increase in RBF at the level of the cortex (p < 0.001; ∆RBF = +71.15 mL/100 g/min), medulla (p < 0.001; ∆RBF = +20.83 mL/100 g/min) and whole kidney (p < 0.001; ∆RBF = +47.58 mL/100 g/min), as observed with the TLCO method. A significant increase in RBF was also observed with the TEO method for the whole kidney during dopamine infusion compared to 0.9% NaCl infusion (p < 0.001; ∆RBF = +54.31 mL/100 g/min).
Figures 3E,F illustrates the mean ± SD R2* values during 0.9% NaCl infusion (control) and dopamine infusion of each layer using the TLCO method and of each segment obtained using the TEO method. With the TLCO method, a lower R2* was observed at the cortex, medulla and whole kidney during dopamine infusion compared to 0.9% NaCl infusion. However, a significant decrease in R2* was only observed at the medulla (p < 0.01; ∆R2* = −2.83 s−1) and the whole kidney (p < 0.001; ∆R2* = −2.19 s−1). The difference in R2* at the cortex was −0.72 s−1 (p = 0.50). With the TEO method, the whole kidney had a lower R2* during dopamine infusion compared to 0.9% NaCl infusion, but this difference was not significant (p = 0.08; ∆R2* = 1.36 s−1).
Table 3 summarizes the mean RBF (mean of 2 runs) and R2* during 0.9% NaCl infusion (control) and dopamine infusion of the renal cortex, medulla, and the whole kidney in both kidneys using the TLCO method.
Table 3. Regional mean and standard deviation of RBF and R2* during 0.9% NaCl infusion (control) and dopamine infusion derived with the TLCO method.
Of the DCE-MRI analysis methods that showed a significant correlation with ASL-MRI, only the model-free model with the ROI caudal to the right renal artery could detect a significant difference in RBF induced by dopamine. This DCE-MRI analysis method showed a significantly higher RBF (∆RBF = 100.65, p < 0.05) during dopamine infusion than 0.9% NaCl infusion. In Figure 5 the DCE-MRI images, as well as the Ktrans, RBF and R2* maps are presented of two subjects during each treatment.
Figure 5. Illustration of DCE-MRI images and MRI maps for Ktrans, RBF and R2* during control and dopamine treatment for two subjects. 1: subject one, 2: subject two. MRI images during control treatment: (A) Ktrans map, (B) RBF-map, (C) R2* map, and (D) DCE image. MRI images during dopamine treatment: (E) Ktrans map, (F) RBF-map, (G) R2* map, and (H) DCE image. Ktrans: volume transfer constant, RBF: renal blood flow, R2*: apparent relaxation rate, DCE: dynamic contrast enhanced.
Quantitative CEUS showed no significant differences in renal perfusion parameters between dogs receiving a 0.9% NaCl and dopamine infusion. Figure 6 shows an overview of the renal perfusion parameters during 0.9% NaCl infusion and dopamine infusion. In addition, no significant effect of dopamine was observed on kidney volume (p = 0.60).
Figure 6. Mean of CEUS renal perfusion parameters during 0.9% NaCl infusion (control) and dopamine infusion. LK, left kidney; RK, right kidney; TTP, time to peak; FT, fall time; RT, rise time; PI, peak intensity; AUC, area under the curve; WiR, wash-in rate; a.u, arbitrary units; s, seconds. There were no significant differences.
4 Discussion
This study was conducted to evaluate the potential role of the emerging imaging technique fMRI in assessing renal perfusion and oxygenation in dogs. Among the DCE-MRI analyses methods, the model-free model with caudal placement of the AIF ROI in the aorta showed most potential in dogs. Additionally, R2* negatively correlated with ASL measured RBF at the cortex and medulla, as well as with medullary WiR and PI. ASL measured RBF, on its turn, showed a positive correlation with cortical WiR, PI, AUC and FT, and with medullary WiR and PI, but a negative correlation with medullary RT. During a pharmacological challenge with dopamine, BOLD-MRI observed a significant decrease in R2* (higher oxygenation) at the medulla and entire kidney, while ASL-MRI demonstrated a significant increase in renal perfusion at the cortex, medulla and the entire kidney.
In renal DCE-MRI, a bolus of gadolinium-based contrast agent (CA) is injected intravenously and its progression through the kidney is tracked by acquiring a time series of T1-weighted images (21, 22). Applying a pharmacokinetic (PK) model to the dynamically acquired tissue concentration curve allows quantification of physiological parameters like RBF and GFR (21, 22). Several PK models have been described for determining renal perfusion and filtration based on DCE-MRI data, but there is no consensus about which model is most suitable (21, 22). Among the models to measure renal perfusion are the two-compartment filtration model (C2FM), the two-compartment uptake model (C2UM), the model-free deconvolution model (MF) and maximum slope model (MS) (22). All of these PK models require measurements of the arterial input function (AIF), which describes the changes in CA concentration over time in a major feeding artery of the kidney, usually the aorta (21, 22). The AIF is defined by the MRI signal obtained from a selected region in the aorta (22). The size and placement of this region in the aorta, however, affect renal perfusion values significantly (67, 68). Nevertheless, there is also no agreed standard for drawing the ROI in the aorta for AIF determination (67, 68). In the present study, the uptake model and the model-free model yielded results that were significantly correlated with RBF measured by ASL-MRI. As for the model-free model, this specifically pertains to analyses where the AIF ROI was positioned 25 mm cranial, 12 mm cranial, or caudal to the right renal artery. For the uptake model, this specifically concerned the analysis in which the AIF ROI was positioned 25 mm cranial to the right renal artery.
There was a significant correlation between several CEUS parameters and RBF determined by ASL-MRI at baseline. At the cortex, WiR, PI, AUC and FT were positively correlated with RBF. At the medulla, WiR and PI correlated positively with RBF, while RT showed a negative correlation. This was in line with the expectations as both RBF and CEUS parameters reflect renal perfusion, although CEUS provides a semi-quantitative assessment (16). With CEUS, a graph of contrast intensity versus time is generated following a bolus injection of contrast agent, from which parameters that semi-quantitatively describe blood perfusion are derived (16). The CEUS parameters extracted from the TIC curve can be divided into two categories: intensity-related parameters that reflect blood volume and time-related parameters that reflect blood velocity (16). The correlations in the current study therefore indicate that elevated renal blood volume and a faster rate of blood inflow are related to increased renal blood flow (16).
A moderate negative correlation between R2* and RBF determined by ASL-MRI was observed at baseline for the cortex and the medulla. As R2* reflects the absolute amount of deoxy-Hb per tissue volume, the negative correlation between R2* and RBF indicates: the lower the perfusion blood flow, the higher the deoxyhemoglobin content and therefore the lower the renal oxygenation (35). A reduction in renal blood flow leads to a decrease in oxygen delivery, which in turn lowers the oxyhemoglobin to deoxyhemoglobin ratio and results in a higher R2*. In human renal transplant patients, Peng et al. (69) also found a negative correlation between ASL and BOLD parameters. There was also a negative correlation between R2* and the CEUS parameters PI and WiR, but only at the level of the medulla. This further indicates that there is a relationship between decreased medullary perfusion and lower oxygenation.
Dopamine is a catecholamine that functions as a neurotransmitter and a hormone. Dopamine is a regulator of the circulatory, renal, and neurohumoral system, showing a dose-dependent effect. At low doses (≤3 μg/kg/min), dopamine increases renal blood flow, GFR, and sodium excretion through dilatation of renal vessels and inhibition of tubular ion reabsorption (70). At medium doses (3 to 7 μg/kg/min), dopamine additionally induces systemic hemodynamic changes such as increases in cardiac output, heart rate and reduction in peripheral resistance (70). At high doses (> 10 μg/kg/min), dopamine exerts a reverse effect through stimulation of α-adrenergic receptors. These high doses stimulate sodium and water reabsorption and induce vasoconstriction, resulting in decreased sodium excretion, urinary output, renal blood flow and GFR (70). Medium and high dopamine dose infusions can provoke arrhythmia, which requires careful subject monitoring (71). In this study, the dose of dopamine was determined based on previous studies in dogs investigating the effects of dopamine on kidney functions and hemodynamics (71–73). Dopamine at a low dose was chosen since higher doses do not induce a greater change in renal perfusion, but do affect systemic hemodynamics and cause arrhythmias, posing a greater anesthetic risk (71–73). In veterinary medicine, dopamine is used in the prevention and treatment of kidney and heart disease (72). Recently, however, researchers have questioned dopamine’s effectiveness in treating these conditions (70).
In this study, ASL-MRI was able to measure a significant increase in renal blood flow following dopamine infusion compared to control in the whole kidney, cortex and medulla. The magnitude of that increase was greater in the cortex than in the medulla. On average, renal perfusion during dopamine infusion increased by 20.1, 19.6 and 13.3%, at the entire kidney, cortex and medulla, respectively. The extent of dopamine’s effect on renal perfusion is subject to variation in the literature. In agreement with our results, a previous study investigating dopamine arrhythmogenicity using an electromagnetic flowmeter reported an increase in RBF of approximately 20% in dogs under anesthesia following a dose of 3 μg/kg/min (71). However, at the same dose, another canine study using colored microspheres observed an increase in RBF of 67.63 ± 23.08% (72). In healthy human volunteers, dopamine infusion at 3 g/kg/min significantly increased the effective renal plasma flow (ERPF) by 25–50% as measured with urine or plasma clearance of radioiodine-labelled-hippuran (74, 75).
Using BOLD-MRI, a decrease in R2* could be detected during dopamine infusion compared to 0.9% NaCl infusion, suggesting an increase in renal oxygenation. The difference in R2* reached statistical significance in the medulla and whole kidney, but not in the cortex. On average, dopamine infusion reduced the medullary R2* by 20% and the total renal R2* by 9.1%. The absence of a significant decrease in R2* for the cortex, in contrast to the medulla, may be due to the lower sensitivity of BOLD-MRI to detect changes in cortical oxygenation (47). This is because of the different relative positions of medullary and cortical oxygenation on the oxygen dissociation curve (47). Consequently, a greater change in oxygenation is required to produce a similar change in R2* for the cortex as for the medulla (47). Our results are in line with those by Redfors et al. (75) who examined the impact of low-dose dopamine on renal hemodynamics and oxygenation in humans, post-cardiac surgery. Their research has shown that dopamine increases renal oxygenation by 28–34% together with a 45–55% increase in RBF (75, 76).
The assumption that dopamine always increases renal oxygen consumption is contradicted by the findings of the current study and the study conducted by Redfors et al. (75) and Ricksten et al. (76). Dopamine administration was thought to increase renal oxygen consumption because it would increase the amount of Na that is supplied to the tubular cells for active transport-mediated reabsorption (75, 76). Dopamine was suggested to increase Na load in the tubular filtrate in two different ways. Firstly, by inhibiting Na reabsorption in the proximal tubule, more Na would flow into the distal tubule possibly increasing regional oxygen consumption (76). Secondly, a rise in renal blood flow due to dopamine-induced vasodilation was expected to coincide with a proportional increase in GFR and sodium filtration (75). However, the pronounced increase in renal blood flow found in the study by Redfors et al. (75) was not accompanied by an increase in GFR. Renal oxygen consumption remained unaffected as a result of dopamine’s lack of influence on GFR, leading to improvements in renal oxygenation (75). Most likely, balanced vasodilation of pre- and post-glomerular arterioles is responsible for the dopamine-induced increase in renal blood flow without significant effects on GFR (75, 76). In vitro experiments using isolated rabbit renal arterioles have shown that dopamine mainly acts on dopamine-1 receptors, inducing equal dilatation of afferent and efferent arterioles (77). This probably explains the findings of our study as well.
In the present study, renal perfusion parameters of CEUS did not show a significant difference between 0.9% NaCL infusion and dopamine infusion. In contrast, CEUS has demonstrated the ability to detect drug-induced renal perfusion changes in previous studies. Several human studies showed that CEUS can identify changes in renal perfusion caused by pharmaceuticals (78, 79). In dogs, contrast-enhanced power Doppler imaging detected significant changes in perfusion parameters induced by norepinephrine (80). It must be highlighted, however, that the magnitude of the significant changes induced was substantial, ranging from 47 to 92% in these human studies (78, 79) and 62 to 94% in the canine study (80). In dogs, the anticipated rise in RBF due to dopamine varies from 20 to 50% (71, 72). Possibly the dose of dopamine used did not affect renal perfusion sufficiently for a significant difference to be observed with CEUS. In fact, CEUS is known for its significant variability in the quantification of tissue perfusion, which is primarily due to patient-related factors, scanner settings, and contrast agent-related factors (81). The high inherent variability of CEUS leads to diagnostic uncertainty (81).
Among the DCE analysis methods significantly correlated with RBF measured by ASL-MRI, only the model-free model with the ROI caudal to the right renal artery could detect a significant difference in RBF with dopamine infusion compared to 0.9% NaCl infusion. These results may indicate that the model-free model with placement of the AIF ROI caudal to the right renal artery has great potential to assess renal perfusion in the dogs. However, additional research is required to determine the most suitable DCE analysis method for dogs. In future studies, it would be beneficial to use a golden reference method to calculate the true RBF.
There were some limitations to the present study. First, reference standards for quantifying renal blood flow, renal oxygenation, and GFR were not included in this study. The current gold standard methods for measuring renal blood flow and tissue oxygenation are invasive, requiring surgical intervention or animal euthanasia. These techniques were not used for animal welfare reasons. It should be noted, however, that the imaging techniques used in this study have previously been validated in animal studies against gold standard techniques. ASL-MRI has been compared with microspheres (38) and ultrasound flowmeter (37) in a porcine model, as well as with histology in rodent models (82) which all showed a good correlation. Comparisons of BOLD-MRI with oxygen-sensitive microelectrodes in multiple pig models have demonstrated that regional oxygen content of the kidneys can be reliably estimated by BOLD-MRI (48, 49). For DCE-MRI, a good correlation has been established between measurements of renal perfusion and those determined with colored microspheres in swine (83) and transit-time flow probes in dogs and rabbits (32, 84). A second limitation might be that only one mid-coronal slice was used to determine ASL-, BOLD-, and DCE-MRI biomarkers. Since the medullary volume fraction of the kidney is higher in a mid-coronal slice than in peripheral slices, this approach might not be an exact representation of the whole kidney (85). However, a study by Winter et al. (86) found that DCE-MRI perfusion parameters from a single slice yielded similar values to those obtained from whole kidneys in humans. Furthermore, in the current study, perfusion and oxygenation were also determined for the cortex and medulla separately with ALS- and BOLD-MRI. Another limitation is that the dopamine infusion and control infusion were performed two to 3 weeks apart. This is because DCE-MRI requires the use of a contrast agent; hence, repeating the MRI scans on the same day would have resulted in erroneous readings due to the residual contrast agent still present in the body. Moreover, the fMRI scan and CEUS examination were performed on two consecutive days since they were performed at different locations. As a result of the time interval between the different examinations, changes in renal blood flow and oxygenation may have been caused by spontaneous fluctuations or time-dependent effects. However, several systemic hemodynamic parameters, including blood pressure, were regularly monitored during the examinations to eliminate any possibility that these factors might have been responsible for the changes observed in renal blood flow or oxygenation.
5 Conclusion
In conclusion, ASL- and BOLD-MRI can measure pharmacologically induced changes in renal blood flow and renal oxygenation in dogs. In contrast, no significant changes in renal perfusion parameters could be observed with CEUS. This might indicate that fMRI modalities allow for the detection of changes that cannot be observed with CEUS. Our results further indicate a correlation between renal perfusion and oxygenation. Regarding DCE-MRI, the model-free model with placement of the AIF ROI caudal to the right renal artery appears to be a promising analysis method in dogs. The use of multiparametric MRI can provide insight into (patho) physiological processes in the kidney and evaluate treatment outcomes. However, further research is needed to confirm the potential of ASL- and BOLD-MRI in dogs and to clarify which analysis method is most suitable for DCE-MRI in dogs.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Animal Ethics Committee from the Faculty of Veterinary Medicine and the Faculty of Bioscience Engineering of Ghent University, Belgium (Approval number: EC2021-085). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
AH: Writing – review & editing, Writing – original draft, Visualization, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. LS: Writing – review & editing, Visualization, Validation, Software, Resources, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. SB: Writing – review & editing, Resources, Investigation, Data curation. BB: Conceptualization, Writing – review & editing, Validation, Methodology, Formal analysis. MH: Writing – review & editing, Resources, Conceptualization. EV: Writing – review & editing, Investigation. JG: Writing – review & editing, Investigation. ES: Writing – review & editing, Supervision, Methodology, Investigation, Funding acquisition, Conceptualization. PP: Software, Resources, Investigation, Writing – review & editing, Supervision, Methodology, Conceptualization. KV: Funding acquisition, Writing – review & editing, Supervision, Methodology, Conceptualization.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was funded by the Special Research Fund of Ghent University, Belgium (BOF grant 01D27919).The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
We thank Siemens Healthineers, and in particular Bernd Kühn, for the support and provision of the research software for body ASL to enable this study. The authors would also like to thank Fien Verdoodt, Sofie Salden, Camila Baptista da Silva and Yangfeng Xu, colleagues of the Department of Morphology, Imaging, Orthopedics, Rehabilitation and Nutrition; Bruna Alonso of the Department of Large Animal Surgery, Anaesthesia and Orthopaedics; the Department of Translational Physiology, Infectiology and Public Health - Comparative Physiology and the Small Animal Department for their help with sample handling and data collection. Additionally, the authors would like to express their gratitude to the animal caretakers Amy DeLuycker and Kjell Meganck who took excellent care of the Beagles.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1406343/full#supplementary-material
SUPPLEMENTARY DATA SHEET S1 | Excel file with physiological parameters of the test subjects during MRI examination.
SUPPLEMENTARY DATA SHEET S2 | Excel file with renal blood flow measurements with ASL-MRI using the TEO and TLCO method.
SUPPLEMENTARY DATA SHEET S3 | Excel file with renal blood flow measurements with DCE-MRI.
SUPPLEMENTARY DATA SHEET S4 | Excel file with oxygenation measurements with BOLD-MRI using the TEO and TLCO method.
SUPPLEMENTARY DATA SHEET S5 | Excel file with physiological parameters of the test subjects during CEUS examination.
SUPPLEMENTARY DATA SHEET S6 | Excel file with renal blood flow parameters measured with CEUS.
References
1. Kellum, JA, and Hoste, EAJ. Acute kidney injury: epidemiology and assessment. Scand J Clin Lab Invest. (2008) 68:6–11. doi: 10.1080/00365510802144813
2. Thoen, ME, and Kerl, ME. Characterization of acute kidney injury in hospitalized dogs and evaluation of a veterinary acute kidney injury staging system. J Vet Emerg Crit Care. (2011) 21:648–57. doi: 10.1111/j.1476-4431.2011.00689.x
3. O’Neill, DG, Elliott, J, Church, DB, McGreevy, PD, Thomson, PC, and Brodbelt, DC. Chronic kidney disease in dogs in UK veterinary practices: prevalence, risk factors, and survival. J Vet Intern Med. (2013) 27:814–21. doi: 10.1111/jvim.12090
4. Levin, A, Tonelli, M, Bonventre, J, Coresh, J, Donner, JA, Fogo, AB, et al. Global kidney health 2017 and beyond: a roadmap for closing gaps in care, research, and policy. Lancet. (2017) 390:1888–917. doi: 10.1016/S0140-6736(17)30788-2
5. Finch, N, and Heiene, R. Early detection of chronic kidney disease In: J Elliott, GF Grauer, and JL Westropp, editors. BSAVA manual of canine and feline nephrology and urology. Gloucester: BSAVA (2017). 130–42.
6. Twardock, AR, and Bahr, A. Renal scintigraphy In: GB Daniel and CR Berry, editors. Textbook of veterinary nuclear medicine. Raleigh, NC: North Carolina State Univ (2006). 330–50.
7. Tyson, R, and Daniel, GB. Renal scintigraphy in veterinary medicine. Semin Nucl Med. (2014) 44:35–46. doi: 10.1053/j.semnuclmed.2013.08.005
8. Krawiec, D, Badertscher, R, Twardock, A, Rubin, S, and Gelberg, H. Evaluation of 99mTc-diethylenetriaminepentaacetic acid nuclear imaging for quantitative determination of the glomerular filtration rate of dogs. Am J Vet Res. (1986) 47:2175–9.
9. Barthez, PY, Hornof, WJ, Cowgill, LD, Neal, LA, and Mickel, P. Comparison between the scintigraphic uptake and plasma clearance of 99MTc-diethylenetriaminepentacetic acid (DTPA) for the evaluation of the glomerular filtration rate in dogs. Vet Radiol Ultrasound. (1998) 39:470–4. doi: 10.1111/j.1740-8261.1998.tb01636.x
10. Mahmoud, H, Buchanan, C, Francis, ST, and Selby, NM. Imaging the kidney using magnetic resonance techniques: structure to function. Curr Opin Nephrol Hypertens. (2016) 25:487–93. doi: 10.1097/MNH.0000000000000266
11. Lankadeva, YR, Okazaki, N, Evans, RG, Bellomo, R, and May, CN. Renal medullary hypoxia: a new therapeutic target for septic acute kidney injury? Semin Nephrol. (2019) 39:543–53. doi: 10.1016/j.semnephrol.2019.10.004
12. Ross, L. Acute kidney injury in dogs and cats. Vet Clin North Am. (2011) 41:1–14. doi: 10.1016/j.cvsm.2010.09.003
13. Dong, Y, Wang, W, Cao, J, Fan, P, and Lin, X. Quantitative evaluation of contrast-enhanced ultrasonography in the diagnosis of chronic ischemic renal disease in a dog model. PLoS One. (2013) 8:e70337–7. doi: 10.1371/journal.pone.0070337
14. Niendorf, T, Pohlmann, A, Arakelyan, K, Flemming, B, Cantow, K, Hentschel, J, et al. How BOLD is blood oxygenation level-dependent (BOLD) magnetic resonance imaging of the kidney? Opportunities, challenges and future directions. Acta Physiol. (2015) 213:19–38. doi: 10.1111/apha.12393
15. Odudu, A, Nery, F, Harteveld, AA, Evans, RG, Pendse, D, Buchanan, CE, et al. Arterial spin labelling MRI to measure renal perfusion: a systematic review and statement paper. Nephrol Dial Transplant. (2018) 33:ii15–21. doi: 10.1093/ndt/gfy180
16. Sencha, AN, and Patrunov, YN. Contrast-enhanced ultrasound. Cham: Springer Nature Switzerland AG (2022).
17. Mannucci, T, Lippi, I, Rota, A, and Citi, S. Contrast enhancement ultrasound of renal perfusion in dogs with acute kidney injury. J Small Anim Pract. (2019) 60:471–6. doi: 10.1111/jsap.13001
18. Stock, E, Paepe, D, Daminet, S, Vandermeulen, E, Duchateau, L, Saunders, JH, et al. Contrast-enhanced ultrasound examination for the assessment of renal perfusion in cats with chronic kidney disease. J Vet Intern Med. (2018) 32:260–6. doi: 10.1111/jvim.14869
19. Dong, Y, Wang, WP, Cao, J, Fan, P, and Lin, X. Early assessment of chronic kidney dysfunction using contrast-enhanced ultrasound: a pilot study. Br J Radiol. (2014) 87:20140350. doi: 10.1259/bjr.20140350
20. Foster, S, Hyvelin, JM, Lakshman, M, Needles, A, Tardy, I, and Tranquart, F. Ultrasound In: F Kiessling, BJ Pichler, and P Hauff, editors. Small Animal Imaging. Cham: Springer International Publishing AG (2011). 325–55.
21. Zhang, JL, and Lee, VS. Renal perfusion imaging by MRI. J Magn Reson Imaging. (2020) 52:369–79. doi: 10.1002/jmri.26911
22. Khalifa, F, Soliman, A, el-Baz, A, Abou el-Ghar, M, el-Diasty, T, Gimel'farb, G, et al. Models and methods for analyzing DCE-MRI: a review. Med Phys. (2014) 41:124301. doi: 10.1118/1.4898202
23. Costelloe, CM, Amini, B, and Madewell, JE. Risks and benefits of gadolinium-based contrast-enhanced MRI. Semin Ultrasound CT MRI. (2020) 41:170–82. doi: 10.1053/j.sult.2019.12.005
24. Sourbron, SP, Michaely, HJ, Reiser, MF, and Schoenberg, SO. MRI-measurement of perfusion and glomerular filtration in the human kidney with a separable compartment model. Investig Radiol. (2008) 43:40–8. doi: 10.1097/RLI.0b013e31815597c5
25. Lüdemann, L, Nafz, B, Eisner, F, Groszlige-Siestrup, C, Meissler, M, Kaufels, N, et al. Absolute quantification of regional renal blood flow in swine by dynamic contrast-enhanced magnetic resonance imaging using a blood pool contrast agent. Investig Radiol. (2009) 44:125–34. doi: 10.1097/RLI.0b013e318193598c
26. Annet, L, Hermoye, L, Peeters, F, Jamar, F, Dehoux, JP, and Van Beers, BE. Glomerular filtration rate: assessment with dynamic contrast-enhanced MRI and a cortical-compartment model in the rabbit kidney. J Magn Reson Imaging. (2004) 20:843–9. doi: 10.1002/jmri.20173
27. Winter, JD, Lawrence, KS, and Cheng, HLM. Quantification of renal perfusion: comparison of arterial spin labeling and dynamic contrast-enhanced MRI. J Magn Reson Imaging. (2011) 34:608–15. doi: 10.1002/jmri.22660
28. Baumann, D, and Rudin, M. Quantitative assessment of rat kidney function by measuring the clearance of the contrast agent Gd (DOTA) using dynamic MRI. Magn Reson Imaging. (2000) 18:587–95. doi: 10.1016/S0730-725X(00)00134-X
29. Zöllner, FG, Schock-Kusch, D, Bäcker, S, Neudecker, S, Gretz, N, and Schad, LR. Simultaneous measurement of kidney function by dynamic contrast enhanced MRI and FITC-sinistrin clearance in rats at 3 tesla: initial results. PLoS One. (2013) 8:11. doi: 10.1371/journal.pone.0079992
30. Jiang, K, Tang, H, Mishra, PK, Macura, SI, and Lerman, LO. Measurement of murine single-kidney glomerular filtration rate using dynamic contrast-enhanced MRI. Magn Reson Med. (2018) 79:2935–43. doi: 10.1002/mrm.26955
31. Privratsky, JR, Wang, N, Qi, Y, Ren, J, Morris, BT, Hunting, JC, et al. Dynamic contrast-enhanced MRI promotes early detection of toxin-induced acute kidney injury. Am J Physiol Ren Physiol. (2019) 316:F351–9. doi: 10.1152/ajprenal.00416.2018
32. Aumann, S, Schoenberg, SO, Just, A, Briley-Saebo, K, Bjørnerud, A, Bock, M, et al. Quantification of renal perfusion using an intravascular contrast agent (part 1): results in a canine model. Magn Reson Med. (2003) 49:276–87. doi: 10.1002/mrm.10380
33. Schoenberg, SO, Aumann, S, Just, A, Bock, M, Knopp, MV, Johansson, LO, et al. Quantification of renal perfusion abnormalities using an intravascular contrast agent (part 2): results in animals and humans with renal artery stenosis. Magn Reson Med. (2003) 49:288–98. doi: 10.1002/mrm.10383
34. Mehl, JN, Lüpke, M, Brenner, AC, Dziallas, P, Wefstaedt, P, and Seifert, H. Measurement of single kidney glomerular filtration rate in dogs using dynamic contrast-enhanced magnetic resonance imaging and the Rutland-Patlak plot technique. Acta Vet Scand. (2018) 60:72. doi: 10.1186/s13028-018-0423-3
35. Ritt, M, Janka, R, Schneider, MP, Martirosian, P, Hornegger, J, Bautz, W, et al. Measurement of kidney perfusion by magnetic resonance imaging: comparison of MRI with arterial spin labeling to Para-aminohippuric acid plasma clearance in male subjects with metabolic syndrome. Nephrol Dial Transplant. (2010) 25:1126–33. doi: 10.1093/ndt/gfp639
36. Wu, WC, St Lawrence, KS, Licht, DJ, and Wang, DJJ. Quantification issues in arterial spin labeling perfusion magnetic resonance imaging. Top Magn Reson Imaging. (2010) 21:65–73. doi: 10.1097/RMR.0b013e31821e570a
37. Warmuth, C, Nagel, S, Hegemann, O, Wlodarczyk, W, and Lüdemann, L. Accuracy of blood flow values determined by arterial spin labeling: a validation study in isolated porcine kidneys. J Magn Reson Imaging. (2007) 26:353–8. doi: 10.1002/jmri.21011
38. Artz, NS, Wentland, AL, Sadowski, EA, Djamali, A, Grist, TM, Seo, S, et al. Comparing kidney perfusion using noncontrast arterial spin labeling MRI and microsphere methods in an interventional swine model. Investig Radiol. (2011) 46:124–31. doi: 10.1097/RLI.0b013e3181f5e101
39. Zhang, Y, Wang, J, Yang, X, Wang, X, Zhang, J, Fang, J, et al. The serial effect of iodinated contrast media on renal hemodynamics and oxygenation as evaluated by ASL and BOLD MRI. Contrast Media Mol Imaging. (2012) 7:418–25. doi: 10.1002/cmmi.1468
40. Williams, DS, Zhang, W, Koretsky, AP, and Adler, S. Perfusion imaging of the rat kidney with MR. Radiology. (1994) 190:813–8. doi: 10.1148/radiology.190.3.8115632
41. Chen, WB, Liang, L, Zhang, B, Liu, CL, Liu, HJ, Luo, HY, et al. To evaluate the damage of renal function in CIAKI rats at 3T: using ASL and BOLD MRI. Biomed Res Int. (2015) 2015:593060. doi: 10.1155/2015/593060
42. Wu, WC, Su, MY, Chang, CC, Tseng, WYI, and Liu, KL. Renal perfusion 3-T MR imaging: a comparative study of arterial spin labeling and dynamic contrast-enhanced techniques. Radiology. (2011) 261:845–53. doi: 10.1148/radiol.11110668
43. Shimizu, K, Kosaka, N, Fujiwara, Y, Matsuda, T, Yamamoto, T, Tsuchida, T, et al. Arterial transit time-corrected renal blood flow measurement with pulsed continuous arterial spin labeling MR imaging. Magn Reson Med Sci. (2017) 16:38–44. doi: 10.2463/mrms.mp.2015-0117
44. Bane, O, Mendichovszky, IA, Milani, B, Dekkers, IA, Deux, JF, Eckerbom, P, et al. Consensus-based technical recommendations for clinical translation of renal BOLD MRI. Magn Reson Mater Phys Biol Med. (2020) 33:199–215. doi: 10.1007/s10334-019-00802-x
45. Piskunowicz, M, Hofmann, L, Zuercher, E, Bassi, I, Milani, B, Stuber, M, et al. A new technique with high reproducibility to estimate renal oxygenation using BOLD-MRI in chronic kidney disease. Magn Reson Imaging. (2015) 33:253–61. doi: 10.1016/j.mri.2014.12.002
46. Milani, B, Ansaloni, A, Sousa-Guimaraes, S, Vakilzadeh, N, Piskunowicz, M, Vogt, B, 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–105. doi: 10.1093/ndt/gfw362
47. Cox, EF, Buchanan, CE, Bradley, CR, Prestwich, B, Mahmoud, H, Taal, M, et al. Multiparametric renal magnetic resonance imaging: validation, interventions, and alterations in chronic kidney disease. Front Physiol. (2017) 8:696. doi: 10.3389/fphys.2017.00696
48. Pedersen, M, Dissing, TH, Morkenborg, J, Stodkilde-Jorgensen, H, Hansen, LH, Pedersen, LB, et al. Validation of quantitative BOLD MRI measurements in kidney: application to unilateral ureteral obstruction. Kidney Int. (2005) 67:2305–12. doi: 10.1111/j.1523-1755.2005.00334.x
49. Warner, L, Glockner, JF, Woollard, J, Textor, SC, Romero, JC, and Lerman, LO. Determinations of renal cortical and medullary oxygenation using blood oxygen level-dependent magnetic resonance imaging and selective diuretics. Investig Radiol. (2011) 46:41–7. doi: 10.1097/RLI.0b013e3181f0213f
50. Wang, Y, Zhang, X, Hua, Z, Xie, L, Jiang, X, Wang, R, et al. Blood oxygen level-dependent imaging and Intravoxel incoherent motion MRI of early contrast-induced acute kidney injury in a rabbit model. Kidney Blood Press Res. (2019) 44:496–512. doi: 10.1159/000500052
51. Wang, Y, Ren, K, Liu, Y, Sun, W, Wang, J, Zhang, X, et al. Application of BOLD MRI and DTI for the evaluation of renal effect related to viscosity of iodinated contrast agent in a rat model. J Magn Reson Imaging. (2017) 46:1320–31. doi: 10.1002/jmri.25683
52. Niles, DJ, Gordon, JW, and Fain, SB. Effect of anesthesia on renal R2 * measured by blood oxygen level-dependent MRI. NMR Biomed. (2015) 28:811–7. doi: 10.1002/nbm.3314
53. Ries, M, Basseau, F, Tyndal, B, Jones, R, Deminière, C, Catargi, B, et al. Renal diffusion and BOLD MRI in experimental diabetic nephropathy. J Magn Reson Imaging. (2003) 17:104–13. doi: 10.1002/jmri.10224
54. Michaely, HJ, Metzger, L, Haneder, S, Hansmann, J, Schoenberg, SO, and Attenberger, UI. Renal BOLD-MRI does not reflect renal function in chronic kidney disease. Kidney Int. (2012) 81:684–9. doi: 10.1038/ki.2011.455
55. Lee, SK, Lee, J, Jang, S, Lee, E, Jeon, CY, Lim, KS, et al. Quantification of renal T2 relaxation rate by use of blood oxygen level–dependent magnetic resonance imaging before and after furosemide administration in healthy beagles. Am J Vet Res. (2021) 82:880–9. doi: 10.2460/AJVR.82.11.880
56. Selby, NM, Blankestijn, PJ, Boor, P, Combe, C, Eckardt, K, Eikefjord, E, et al. Magnetic resonance imaging biomarkers for chronic kidney disease: a position paper from the European Cooperation in Science and Technology action PARENCHIMA. Nephrol Dial Transplant. (2018) 33:ii4–ii14. doi: 10.1093/ndt/gfy152
57. Copur, S, Yavuz, F, Sag, AA, Tuttle, KR, and Kanbay, M. Future of kidney imaging: functional magnetic resonance imaging and kidney disease progression. Eur J Clin Investig. (2022) 52:e13765. doi: 10.1111/eci.13765
58. Haers, H, Daminet, S, Smets, PMY, Duchateau, L, Aresu, L, and Saunders, JH. Use of quantitative contrast-enhanced ultrasonography to detect diffuse renal changes in beagles with iatrogenic hypercortisolism. Am J Vet Res. (2013) 74:70–7. doi: 10.2460/ajvr.74.1.70
59. Liu, DJX, Hesta, M, Stock, E, Bogaerts, E, Broeckx, BJG, Saunders, JH, et al. Renal perfusion parameters measured by contrast-enhanced ultrasound in healthy dogs demonstrate a wide range of variability in the long-term. Vet Radiol Ultrasound. (2019) 60:201–9. doi: 10.1111/vru.12690
60. Stock, E, Vanderperren, K, Haers, H, Duchateau, L, Hesta, M, and Saunders, JH. Quantitative differences between the first and second injection of contrast agent in contrast-enhanced ultrasonography of feline kidneys and spleen. Ultrasound Med Biol. (2017) 43:500–4. doi: 10.1016/j.ultrasmedbio.2016.09.013
61. Cox, RW, Ashburner, J, Breman, H, Fissell, K, Haselgrove, C, Holmes, CJ, et al. A (Sort of) new image data format standard: NIfTI-1. In: Proceedings of the 10th annual meeting of organisation of human brain mapping; (2004) 13–17; San Diego, CA, Academic Press (2004).
62. Zhang, X, Peng, J, Wang, C, Feng, Y, Feng, Q, Li, X, et al. Improved liver R2∗ mapping by averaging decay curves. Sci Rep. (2017) 7:1. doi: 10.1038/s41598-017-05683-5
63. Sanmiguel Serpa, LC, De Visschere, P, Speeckaert, M, and Pullens, P. A new method to analyse renal perfusion: a proof of concept. In: Conference presentation abstract 2023 ISMRM & ISMRT annual meeting & exhibition; (2023) 3–8; Toronto.
64. McManus, IC, Stöver, K, and Kim, D. Arnheim’s gestalt theory of visual balance: examining the compositional structure of art photographs and abstract images. Iperception. (2011) 2:615–47. doi: 10.1068/i0445aap
65. Shahzadi, I, Aslam, I, Qazi, SA, and Omer, H. Golden-angle radial sparse parallel MR image reconstruction using SC-GROG followed by iterative soft thresholding. Appl Magn Reson. (2019) 50:977–88. doi: 10.1007/s00723-019-01127-7
66. Schober, P, Boer, C, and Schwarte, LA. Correlation coefficients: appropriate use and interpretation. Anesth Analg. (2018) 126:1763–8. doi: 10.1213/ANE.0000000000002864
67. Mendichovszky, IA, Cutajar, M, and Gordon, I. Reproducibility of the aortic input function (AIF) derived from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) of the kidneys in a volunteer study. Eur J Radiol. (2009) 71:576–81. doi: 10.1016/j.ejrad.2008.09.025
68. Cutajar, M, Mendichovszky, IA, Tofts, PS, and Gordon, I. The importance of AIF ROI selection in DCE-MRI renography: reproducibility and variability of renal perfusion and filtration. Eur J Radiol. (2010) 74:e154–60. doi: 10.1016/j.ejrad.2009.05.041
69. Peng, J, Gao, J, Hong, Y, Wu, Z, Chen, G, and Lu, G. The value of functional magnetic resonance imaging in evaluating renal allograft function. Asian J Surg. (2024) 47:1740–5. doi: 10.1016/j.asjsur.2023.12.121
70. Sigrist, NE. Use of dopamine in acute renal failure. J Vet Emerg Crit Care. (2007) 17:117–26. doi: 10.1111/j.1476-4431.2007.00225.x
71. Sato, T, Washizu, M, Kobayashi, K, Misaka, K, Hayashi, T, Orima, H, et al. Arrhythmogenicity of dopamine and its effects on hemodynamics in dogs under halothane anesthesia. Japanese J Vet Sci. (1989) 51:396–402. doi: 10.1292/jvms1939.51.396
72. Furukawa, S, Nagashima, Y, Hoshi, K, Hirao, H, Tanaka, R, Maruo, K, et al. Effects of dopamine infusion on cardiac and renal blood flows in dogs. J Vet Med Sci. (2002) 64:41–4. doi: 10.1292/jvms.64.41
73. Meyer, MB, McNay, JL, and Goldberg, LI. Effects of dopamine on renal function and hemodynamics in the dog. J Pharmacol Exp Ther. (1967) 156:186–92.
74. Olsen, NV, Hansen, JM, Kanstrup, IL, Richalet, JP, and Leyssac, PP. Renal hemodynamics, tubular function, and response to low-dose dopamine during acute hypoxia in humans. J Appl Physiol. (1993) 74:2166–73. doi: 10.1152/jappl.1993.74.5.2166
75. Redfors, B, Bragadottir, G, Sellgren, J, Swärd, K, and Ricksten, SE. Dopamine increases renal oxygenation: a clinical study in post-cardiac surgery patients. Acta Anaesthesiol Scand. (2010) 54:183–90. doi: 10.1111/j.1399-6576.2009.02121.x
76. Ricksten, SE, Bragadottir, G, and Redfors, B. Renal oxygenation in clinical acute kidney injury. Crit Care. (2013) 17:221–10. doi: 10.1186/cc12530
77. Edwards, RM. Comparison of the effects of fenoldopam, SK&F R-87516 and dopamine on renal arterioles in vitro. Eur J Pharmacol. (1986) 126:167–70. doi: 10.1016/0014-2999(86)90756-9
78. Muskiet, MHA, Emanuel, AL, Smits, MM, Tonneijck, L, Meijer, RI, Joles, JA, et al. Assessment of real-time and quantitative changes in renal hemodynamics in healthy overweight males: contrast-enhanced ultrasonography vs Para-aminohippuric acid clearance. Microcirculation. (2019) 26:e12580. doi: 10.1111/micc.12580
79. Selby, NM, Williams, JP, and Phillips, BE. Application of dynamic contrast enhanced ultrasound in the assessment of kidney diseases. Curr Opin Nephrol Hypertens. (2021) 30:138–43. doi: 10.1097/MNH.0000000000000664
80. Arger, PH, Sehgal, CM, Pugh, CR, Kirchoffer, JI, Kotlar, EY, and Bovee, KC. Evaluation of change in blood flow by contrast-enhanced power Doppler imaging during norepinephrine-induced renal vasoconstriction. J Ultrasound Med. (1999) 18:843–51. doi: 10.7863/jum.1999.18.12.843
81. Tang, MX, Mulvana, H, Gauthier, T, Lim, AKP, Cosgrove, DO, Eckersley, RJ, et al. Quantitative contrast-enhanced ultrasound imaging: a review of sources of variability. Interf Focus. (2011) 1:520–39. doi: 10.1098/rsfs.2011.0026
82. Hueper, K, Gutberlet, M, Rong, S, Hartung, D, Mengel, M, Lu, X, et al. Acute kidney injury: arterial spin labeling to monitor renal perfusion impairment in mice — comparison with histopathologic results and renal function. Radiology. (2014) 270:117–24. doi: 10.1148/radiol.13130367
83. Prasad, PV, Cannillo, J, Chavez, DR, Pinchasin, ES, Dolan, RP, Walovitch, R, et al. First-pass renal perfusion imaging using MS-325, an albumin-targeted MRI contrast agent. Investig Radiol. (1999) 34:566–71. doi: 10.1097/00004424-199909000-00003
84. Montet, X, Ivancevic, MK, Belenger, J, Jorge-Costa, M, Pochon, S, Pechère, A, et al. Noninvasive measurement of absolute renal perfusion by contrast medium-enhanced magnetic resonance imaging. Investig Radiol. (2003) 38:584–92. doi: 10.1097/01.RLI.0000077127.11949.8c
85. Bokacheva, L, Rusinek, H, Zhang, JL, and Lee, VS. Assessment of renal function with dynamic contrast-enhanced MR imaging. Magn Reson Imaging Clin N Am. (2008) 16:597–611. doi: 10.1016/j.mric.2008.07.001
Keywords: kidney, perfusion, oxygenation, ASL-MRI, BOLD-MRI, DCE-MRI, CEUS, dog
Citation: Hillaert A, Sanmiguel Serpa LC, Bogaert S, Broeckx BJG, Hesta M, Vandermeulen E, Germonpré J, Stock E, Pullens P and Vanderperren K (2024) Assessment of pharmacologically induced changes in canine kidney function by multiparametric magnetic resonance imaging and contrast enhanced ultrasound. Front. Vet. Sci. 11:1406343. doi: 10.3389/fvets.2024.1406343
Edited by:
Tommaso Banzato, University of Padua, ItalyReviewed by:
Hakyoung Yoon, Jeonbuk National University, Republic of KoreaSofia Alves-Pimenta, Universidade de Trás-os-Montes e Alto Douro, Portugal
Copyright © 2024 Hillaert, Sanmiguel Serpa, Bogaert, Broeckx, Hesta, Vandermeulen, Germonpré, Stock, Pullens and Vanderperren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katrien Vanderperren, Katrien.Vanderperren@ugent.be
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship