- 1Neuroimaging and Biotechnology Laboratory (NOBEL), Clinical Neurosciences Research Laboratory (LINC), Health Research Institute of Santiago de Compostela (IDIS), Santiago de Compostela, Spain
- 2Centro de Supercomputación de Galicia (CESGA), Santiago de Compostela, Spain
- 3Nasasbiotech S. L., A Coruña, Spain
- 4Translational Stroke Laboratory (TREAT), Clinical Neurosciences Research Laboratory (LINC), Health Research Institute of Santiago de Compostela (IDIS), Santiago de Compostela, Spain
- 5Neurological Sciences and Cerebrovascular Research Laboratory, Department of Neurology and Stroke Center, La Paz University Hospital, Neuroscience Area of IdiPAZ Health Research Institute, Universidad Autónoma de Madrid, Madrid, Spain
- 6NeuroAging Group (NEURAL), Clinical Neurosciences Research Laboratory (LINC), Health Research Institute of Santiago de Compostela (IDIS), Santiago de Compostela, Spain
- 7Centro de Investigación Biomédica en Red en Enfermedades Neurodegenerativas (CIBERNED), Instituto de Salud Carlos III, Madrid, Spain
- 8Department of Neuroradiology, Hospital Clínico Universitario, Health Research Institute of Santiago de Compostela (IDIS), Santiago de Compostela, Spain
Introduction: Dynamic susceptibility-weighted contrast-enhanced (DSC) perfusion studies in magnetic resonance imaging (MRI) provide valuable data for studying vascular cerebral pathophysiology in different rodent models of brain diseases (stroke, tumor grading, and neurodegenerative models). The extraction of these hemodynamic parameters via DSC-MRI is based on tracer kinetic modeling, which can be solved using deconvolution-based methods, among others. Most of the post-processing software used in preclinical studies is home-built and custom-designed. Its use being, in most cases, limited to the institution responsible for the development. In this study, we designed a tool that performs the hemodynamic quantification process quickly and in a reliable way for research purposes.
Methods: The DSC-MRI quantification tool, developed as a Python project, performs the basic mathematical steps to generate the parametric maps: cerebral blood flow (CBF), cerebral blood volume (CBV), mean transit time (MTT), signal recovery (SR), and percentage signal recovery (PSR). For the validation process, a data set composed of MRI rat brain scans was evaluated: i) healthy animals, ii) temporal blood–brain barrier (BBB) dysfunction, iii) cerebral chronic hypoperfusion (CCH), iv) ischemic stroke, and v) glioblastoma multiforme (GBM) models. The resulting perfusion parameters were then compared with data retrieved from the literature.
Results: A total of 30 animals were evaluated with our DSC-MRI quantification tool. In all the models, the hemodynamic parameters reported from the literature are reproduced and they are in the same range as our results. The Bland–Altman plot used to describe the agreement between our perfusion quantitative analyses and literature data regarding healthy rats, stroke, and GBM models, determined that the agreement for CBV and MTT is higher than for CBF.
Conclusion: An open-source, Python-based DSC post-processing software package that performs key quantitative perfusion parameters has been developed. Regarding the different animal models used, the results obtained are consistent and in good agreement with the physiological patterns and values reported in the literature. Our development has been built in a modular framework to allow code customization or the addition of alternative algorithms not yet implemented.
1. Introduction
Dynamic susceptibility-weighted contrast-enhanced (DSC) perfusion studies in magnetic resonance imaging (MRI) provide valuable data for studying vascular cerebral pathophysiology. It is possible to acquire functional information about perfusion-related parameters such as cerebral blood flow (CBF), mean transit time (MTT), and cerebral blood volume (CBV) (Keston et al., 2003; Østergaard, 2004). The accurate quantification of these parameters has several clinical applications, including the identification and evaluation of ischemic stroke prior to treatment (Calamante et al., 2002; Callewaert et al., 2021), the description of lesions associated with multiple sclerosis (Haselhorst et al., 2000; D'haeseleer et al., 2015), the diagnosis of tumors (Stadlbauer et al., 2015; Choi et al., 2016), or as trackers of Alzheimer's disease progression (Lacalle-Aurioles et al., 2014; Warpechowski et al., 2023).
This imaging technique is based on monitoring MRI signal strength variations after the injection of a bolus of a paramagnetic contrast agent, such as gadolinium salt-diethylenetriaminepentaacetic acid (Gd DTPA). On T2*-weighted images, the circulation of the bolus results in a decrease in signal intensity due to small variations in the local magnetic field, displaying, as a result, the time course of this tracer through the tissue. In clinical practice, there are different available software to process these images, contrary to preclinical studies. Nowadays, rodent models of human brain disorders represent more than 80% of the animals in research. In this sense, in vivo perfusion by MRI in rodents is applied to monitor disease progressions, such as stroke, tumor size, neurodegeneration, and evaluation of therapeutic response in longitudinal studies, as well as to develop new animal models of different pathologies (Boisserand et al., 2017; Park et al., 2020; Qi et al., 2021). Therefore, it would be appealing to have a tool capable of performing the whole hemodynamic quantification process in a fast and reliable way for research purposes.
On the technical side, the extraction of these hemodynamic parameters via DSC-MRI is based on tracer kinetic modeling, which can be solved using deconvolution-based methods, among others. Most of the post-processing software used in preclinical studies are home-built and custom-designed, its use being, in most cases, limited to the institution responsible for the development. While there do exist some commercial software tools, they are expensive and specifically designed for clinical use (Gordaliza et al., 2015; López-Larrubia, 2018; Hartmann et al., 2020; Tsai et al., 2021). These tools often only provide relative values (no absolute tissue hemodynamic parameters), do not compute all of the main perfusion parameters (CBF, CBV, and MTT), and overlook direct parameters such as signal recovery (SR) and percentage signal recovery (PSR), which may add valuable diagnostic information without requiring additional measurements (Huhndorf et al., 2016). These restrictions do not allow direct comparison with other software, highlighting the need for an open-source implementation of a DSC-MRI perfusion software application for use in preclinical studies. Furthermore, making it free and open-access could help toward the standardization of DSC-MRI methodology, a pressing issue limiting this method's potential (Boxerman et al., 2020).
In this study, we present the implementation of an open-source DSC quantification tool, named Perfusion-NOBEL, developed as a Python project. Our semi-automatic approach requires no pre-processing aside from the manual delineation of masks, it provides absolute perfusion maps (CBF, CBV, MTT, including SR, and PSR), and it was validated on a large and diverse data set composed of 30 MRI rat brain scans of different models of brain diseases. These data set included i) healthy animals, ii) temporal blood–brain barrier (BBB) dysfunction animals, iii) cerebral chronic hypoperfusion (CCH) model, iv) ischemic stroke model, and v) glioblastoma multiforme (GBM) model. In addition, the resulting hemodynamic parameters were then compared with the literature data. Our tool has been built in a modular framework to allow code customization and the addition of alternative algorithms not yet implemented.
2. Materials and methods
2.1. Theoretical basis and mathematical description
Quantitative analysis of tissue perfusion and blood volume was performed using established tracer kinetic models, which have been extensively reviewed in the literature (Ostergaard et al., 1996; Wu et al., 2003). In brief, the mathematical approach to this process starts from the linear relation between the T2* relaxation time variation [transverse relaxation rate change ()], signal intensity change, and the concentration of the contrast agent. These relations go as follows:
Cm(t) is the concentration of contrast agent measured in the tissue, TE is the echo time, S(t) is the signal intensity at a given time, S0 is the signal intensity before the contrast agent injection, and k is a proportionality constant which, to a first approximation, does not depend on the size and geometry of the vessel, or in the density of the capillary bed.
The tissue response to the arrival of the contrast agent is very sensitive to many factors such as the bolus arrival time and the vascular system, causing the delay and dispersion of the bolus injected. In order to obtain an accurate quantification taking into account these effects, an arterial input function (AIF) is used, measuring the concentration of contrast agent arriving at the vessel. The use of arteries as reference vessels is widely documented in the literature, but it is also possible to use veins [venous output function (VOF)]. The AIF is semi-manually obtained from the DSC-MRI images by delineating and masking the desired region.
In terms of system analysis, the measured concentration curve, Cm(t), the output (response of the tissue to the injection of contrast agent), is related to the ideal concentration curve with no delay or dispersion, k(t), the input, through the convolution of the former with the AIF, the impulse response function, as shown in the following equation:
The * represents the convolution and it can be interpreted as the CAIF(t) modifying the shape of the ideal concentration curve k(t), to give the measured concentration curve Cm(t). Through the inverse operation, known as deconvolution, it is possible to obtain k(t), needed for the following analysis, given that both Cm(t) and CAIF(t) are known. This ideal concentration curve, k(t), it can be also referred as flow-scaled residue function (Fieselmann et al., 2011).
An optional step before the deconvolution process can be taken by fitting the curve CAIF(t) to a gamma variate function (Calamante et al., 2000):
where t0 is the bolus arrival time (BAT) and determines the arrival of the bolus to any given region. This step helps to reduce noise and avoid other effects such as the recirculation of the contrast agent.
The deconvolution step can be computed using different methods. Truncated single value decomposition (TSVD) and Tikhonov regularization are two of these methods, both being appropriate for the resolution of ill-posed problems like this one (Calamante et al., 2000; Fieselmann et al., 2011). In our program, we have implemented both methods of resolution. Then, the quantitative perfusion parameters, CBF, CBV, and MTT, can be obtained from the concentration curves resulting.
CBF and CBV are computed using deconvolution (Fieselmann et al., 2011):
kH being a constant that groups together all the parameters needed to obtain absolute measurements of the perfusion parameters. is the tissue-to-artery concentration scale factor ratio, HA(= 0.45) is the assumed hematocrit in large arterial vessels, HT(= 0.25) is the assumed hematocrit in the capillary bed in the tissue, and ρVoi(= 1.04 g/ml) is the apparent brain density. It should be noted that, CBV can be obtained via an alternative non-deconvolution method, relating the contrast agent in the tissue, Cm(t), with the contrast agent in the AIF, CAIF(t) (Calamante et al., 1999, 2000; Konstas et al., 2009; Fieselmann et al., 2011),
Finally, the MTT values are related to the CBF and CBV through the central volume theorem,
This equation holds true for the volume of interest, and it can be applied independently of the method used to calculate the CBV.
Due to being supported by a reliable mathematical model, CBF, CBV, and MTT are considered the three main perfusion parameters. However, there are other parameters that also provide important information outside those three main quantities. SR and PSR are two of those and they play important roles in the clinical diagnosis of GBM. These variables provide insightful information even when no underpinning mathematical model is present. Contrary to previous hemodynamic parameters, whose calculation requires time-consuming post-processing of the DSC-MRI images, obtaining SR and PSR is much faster and more straightforward. Following (Huhndorf et al., 2016) they can be calculated as,
These two parameters only depend on the signal intensities at different times of the bolus passage through the tissue. Spost is the signal intensity at a time after the bolus arrival, usually 60 s after, Spre is the signal intensity at baseline, before the bolus arrival to the tissue, and Smin is the minimum in signal, corresponding to the peak of the bolus. Restoring baseline signal intensity corresponds to a value of 100% on the PSR and SR maps.
2.2. Software implementation
This study presents an open DSC quantification tool for preclinical studies (Perfusion-Nobel). The implementation was carried out with Python (V 3.8.1), and the following libraries were used: NumPy (V 1.24.1) (Harris et al., 2020), SciPy (V 1.10.0) (Virtanen et al., 2020), Re (V Python 3.9) (Van Rossum, 2020), Matplotlib (V 3.6.2) (Hunter, 2007), Pillow (V 9.4) (Murray et al., 2023), OpenCV (V 4.7.0.68) (Bradski, 2000), Pydicom (V 2.3.1) (Mason, 2011), and IPython (V 8.8.0) (Perez, 2007).
Modular construction has been used to implement the DSC quantification method. As a result, it is possible to easily replace existing stages or add new ones to the quantification workflow (for instance, new pre-processing algorithms or additional fitting models). The following is the processing workflow from the images acquired in the MR system (Figure 1):
i) The inputs requested by the software are: DSC-MRI in DICOM format, the brain mask, the AIF-VOF (PNG or NPY format) mask, and the frame number corresponding to the minimum signal intensity. The masking process should be done beforehand, either semi-manually by drawing the ROIs or using any automatic segmentation tool.
ii) DSC-MRI images are automatically pre-processed, reducing the noise by applying a 2D Gaussian Low Pass Filter (X and Y directions with a 5-pixel diameter and a standard deviation of 0.5 pixels). Then, the PSR and SR maps are computed, shown on screen, and automatically saved (TIFF format). Moreover, the mean image signal intensity curve (S(t)) and mean contrast concentration curve (Cm(t)) are presented.
iii) Next, the AIF or VOF plot is shown, calculated as the mean contrast concentration curve in the AIF or VOF masked region. In brain perfusion, the mean cerebral or anterior carotid arteries are frequently utilized as reference vessels, corresponding to the highest signal intensity regions in T2*-weighted imaging. This function will be the global reference vessel for each pixel in the image.
iv) If we continue the analysis, it will display again the concentration curve corresponding to the reference vessel. Then, two values must be selected: (1) the arrival time of the contrast agent to the tissue and (2) the time after the first pass of the contrast agent through the tissue. These values are needed for the gamma variate function fitting to the concentration curve, and it will be shown on the screen. This process can be repeated for an accurate adjustment.
v) Finally, the deconvolution process will start, and the perfusion maps CBF, CBV, and MTT will be calculated, shown, and automatically saved (TIFF format). The deconvolution step can be done using TSVD or Tikhonov methods.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g001.jpg)
Figure 1. Diagram of the main steps implemented in our software to extract the perfusion-weighted images.
2.3. Validation against a data set
All experimental animal procedures were conducted under procedure numbers: 15011/2021/002, 15011/2021/003, 15011/2023/002, and 15011/2021/001 approved by the Animal Care Committee, according to European Union Rules and the Spanish regulation (2010/63/EU and RD53/2013). Animals were kept in a controlled environment at 22 ± 1°C and 60 ± 5% humidity, with 12:12 h light: darkness cycles, and were fed ad libitum with standard diet pellets and tap water. All surgical procedures and MRI studies were conducted under sevoflurane (Abbott Laboratories, IL, USA) anesthesia (3–4%) using a carrier 70:30 gas mixture of N2O:O2.
Thirty Sprague-Dawley rats with a weight between 250 and 350 g were used to test the developed program in animal models of different pathologies. Five groups of experiments are shown in Figure 2.
2.4. Healthy animals
We used healthy Sprague-Dawley (SD) rats with no surgery or treatment as the control healthy group (n = 6).
2.5. BBB disfunction model
Healthy rats (n = 4) were anesthetized, and X gr/Kg of mannitol [25% mannitol dissolved in isotonic 0.9% sodium chloride solution (B. Braun Medical SA, Barcelona, Spain)] was then intravenously injected (i.v.) in the lateral tail vein. MRI-DSC measurements were made starting 10 min after the mannitol injection (Duong et al., 2000).
2.6. Cerebral chronic hypoperfusion model
A group of rats (n = 4) were induced cerebral chronic hypoperfusion (CCH) by two-vessel occlusion (2VO), also known as permanent bilateral common carotid artery occlusion, as previously described by Cao et al., 2018. The animals were anesthetized, and the right common carotid artery was isolated and ligated with a 4-0, not absorbent suture through a ventral median incision in the cervical area. After 1 week, the left common carotid artery was ligated following the same protocol. MRI-DSC measurements were made 3 weeks after the surgery.
2.7. Glioblastoma multiforme model
GBM was induced in a group of rats (n = 7) by the implementation of the F98 cell line (Bulin et al., 2022). After anesthesia, animals were then mounted onto a stereotaxic frame. A midline incision was performed, and a burr hole was punctured using a 16-gauge needle. The burr hole coordinates were, using the bregma as a reference point, 1 mm anterior and 3 mm lateral to the right. Then 1 × 105 cells suspended in a volume of 5 μl of non-supplemented DMEM were injected into the brain with a 26-gauge needle at a depth of 6 mm from the surface of the skull (injection rate of 0.2 ul/min). MRI-DSC measurements were evaluated 10 and 20 days after the surgery.
2.8. Ischemic stroke animal model
Regarding the stroke group (n = 9), transient focal ischemia was induced in rats by transient middle cerebral artery occlusion (MCAo) following surgical procedures previously described in Vieites-Prado et al. (2016). In brief, using 6-0 silk sutures, the right external carotid artery as well as the pterygopalatine artery of the internal carotid was ligated. A silicon rubber-coated size 4–0 monofilament (diameter 0.19 mm, length 23 mm; diameter with coating 0.37 ± 0.02 mm; coating length 3–4 mm) (Doccol Corporation, Sharon, MA) was inserted into the stump of the right common carotid artery and advanced into the internal carotid artery to 20 mm from the bifurcation to occlude the origin of the MCA. A laser Doppler flow probe (tip diameter 1 mm) attached to a PeriFlux 5000 Laser Doppler Flowmeter (Perimed AB, Stockholm, Sweden) was placed over the thinned skull in the MCA territory (4 mm lateral to bregma) to obtain a continuous measure of relative cerebral blood flow during the experiment Diffusion-weighted imaging (DWI), magnetic resonance angiography (MRA), and MRI-DSC measurements were made starting 20–30 min after the onset of MCA occlusion. The suture was removed 75 min after the occlusion.
2.9. DSC-MRI
All studies were conducted on a Bruker BioSpec 9.4 T MR scanner (horizontal bore magnet with 12 cm wide Bruker BioSpin) equipped with actively shielded gradients (440 mT m−1). Animals were imaged with a combination of a linear birdcage resonator (7 cm in diameter) for signal transmission and a 2 × 2 surface coil array for signal detection, positioned over the head of the animal, which was fixed with a teeth bar, earplugs, and adhesive tape. Animals were physiologically monitored throughout the MR imaging experiments. Transmission and reception coils were actively decoupled from each other.
DSC-MR images were acquired using ultra-fast gradient-echo methods (EPI-T2*) with the following parameters: 6.2 ms of echo time (ET), 1 s of repetition time (RT), number of repetitions (NRs) = 180, 1 average, 50 kHz spectral bandwidth (SW), flip angle (FA) of 90°, and 5 slices of 1.5 mm. Field of view (FOV) of 2.2 × 2.2 cm2, and a matrix size of 128 × 128, giving an in-plane resolution of 172 μm/pixel, implemented without fat suppression. The contrast agent was quickly administrated as a bolus 20 s after starting acquisition [intravenous (i.v.) bolus injection of gadolinium contrast agent (0.3 mmol/kg)]. The total DSC scan acquisition time was 3 min.
Brain volumetry study was also evaluated from T2-wi using a rapid acquisition with relaxation enhancement (RARE) sequence (axial and coronal orientations): with an ET = 11 ms, RT = 2.5 s, Rare Factor (RF) = 8, FA = 180°, NA = 3, SW = 37 KHz, 14 slices of 1 mm, 25.6 × 25.6 mm2 FOV, and a matrix size of 256 × 256 (isotropic in-plane resolution of 100 μm2/pixel).
To evaluate the status of MCAo in a non-invasive manner during occlusion in the stroke animal model, time-of-flight magnetic resonance angiography (TOF-MRA) was performed. The TOF-MRA scan was performed with a 3D FLASH sequence with TE = 2.5 ms, TR = 15 ms, FA = 20°, NA = 2, SW = 98 KHz, 1 slice of 14 mm, FOV = 30.72 × 30.72 × 14 mm3, and a matrix size of 256 × 256 × 58 (resolution of 120 × 120 × 241 μm3/pixel). Moreover, apparent diffusion coefficient (ADC) maps were obtained from diffusion-weighted image (DWI) using a spin echo-planar imaging sequence (DTI-EPI): ET = 26.91 ms; RT = 4 s; SW = 200 KHz; seven b-values of 0, 300, 600, 900, 1,200, 1,600, and 2,000 s/mm2; FA = 90°; NA = 4; 14 consecutive slices of 1 mm; FOV = 24 × 16 mm2, and a matrix size of 96 × 64 (isotropic in-plane resolution of 250 μm2/pixel).
2.10. Data analysis and statistics
We used region of interest (ROI) analysis to quantify the absolute maps CBV, CBF, MTT, and SR, PSR obtained in our new program. Next, through ImageJ software (Rasband W, NIH, Bethesda, MD, USA), we placed a circular ROI over each hemisphere cortex at the plane and a whole-brain ROI. For the ischemic stroke animal model, ROIs were situated in the core, penumbra region of the ipsilateral (IL) hemisphere, and cortex of the contralateral (CL) hemisphere. In the last group, the peripheral tumor area (tumor rim) and core ROIs were placed. Adjacent slices were measured to permit discrimination between intra- and inter-subject variance. To test whether the two slices can be considered samples of the same mean, a paired t-test was performed with the values from the two slices (p < 0.05). We generated Bland–Altman (BA) plots to compare our results (CBV, CBF, and MTT) with the literature data, where the horizontal axis represents the mean value [(our data + literature data)/2], and the vertical axis represents the difference; this method was used to compare two different measurement techniques (Giavarina, 2015). Data are presented as the mean ± 1.96*SD.
3. Results
A total of 30 animals were evaluated with our DSC-MRI quantification tool, and these results were compared in good agreement with values reported in the literature for rats. Table 1 reports our results, as well as a literature overview of absolute perfusion values obtained with different procedures, experimental conditions, rat strains, and post-processing methods.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-t001.jpg)
Table 1. Literature review: absolute perfusion parameters CBV, CBF, and MTT in healthy rats, rats with a mannitol i.v. injection, rat brain hypoperfusion model, ischemic stroke rat model, and glioblastoma rat model.
Our analysis has been developed as a Python script and has been published under a free software license (GNU GPL). Source code and binary downloads are available at https://github.com/MRI-NOBEL/Perfusion-NOBEL. Our tool successfully generated and displayed, in common animal models in neuroscience, all the perfusion parametric images: CBF, CBV, MTT, including PSR, and SR (Figure 3). Moreover, the anatomic details of the rat brain model used have been represented by means of a T2-wi or DWI (ADC) during the MCAo in the axial orientation (Figure 3). Regions of impaired flow under ischemic conditions and tumor influence areas can be delineated. Two distinct regions can be detected during the acute phase of stroke: an ischemic core that is severely and irreversibly damaged, and a penumbra region defined as ischemic tissue that is functionally impaired and at risk of infarction but has the potential to be salvaged (Boisserand et al., 2017). By means of hemodynamic variables, it is possible to accurately differentiate the core and peripheral tumor areas (Aprile et al., 2015). However, changes associated with chronic cerebral hypoperfusion or temporary BBB dysfunction models are only appreciable due to the quantitative values obtained.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g003.jpg)
Figure 3. Examples of mean cerebral l blood flow (CBF), cerebral blood volume (CBV), mean transit time (MTT), signal recovery (SR), and percentage signal recovery (PSR) images generated with our tool for the different animal models. Anatomic MRI of the rat brain model used is detailed with T2-wi or DWI (ADC) (red arrows indicate the ischemic lesion and presence of glioma).
3.1. Quantification of absolute perfusion parameters
Changes in hemodynamic parameters were compared using left, right, and whole-brain ROIs between healthy animals and animals subjected to chronic hypoperfusion or temporal BBB dysfunction (Figure 4A). For all regions, CBF decreases by ~49–35.6% in both models compared to healthy animals. The same is for the CBV, although the reduction is only 32–13.7%, respectively. Mean MTT increased in all ROIs by ~40.7%, and SR or PSR shows a lower capacity to discriminate changes in these models.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g004.jpg)
Figure 4. (A) Comparison of mean cerebral blood flow (CBF), cerebral blood volume (CBV), mean transit time (MTT), signal recovery (SR), and percentage signal recovery (PSR) in selected regions of interest (ROIs) of healthy SD rats, during hypoperfusion model, and mannitol injection. (B) Mean CBF, CBV, MTT, SR, and PSR in contralateral hemisphere, penumbra, and core ROIs during occlusion in an ischemic stroke animal model. (C) Mean CBF, CBV, MTT, SR, and PSR in the peripherical and core ROIs for glioblastoma model at 10 and 20 days after the surgery.
Figure 4B illustrates the hemodynamic parameter evolution in a stroke animal model during MCAo. As can be seen, we found a CBF reduction of 36.4% and 77.2% in the penumbra and core regions in accordance with previous studies (Robertson et al., 2011; Reid et al., 2012; Cipolla et al., 2017), respectively [75.1 ± 21.1 vs. 47.8 ± 13.4 vs. 17.16 ± 9.0 ml/(100g*min)]. Moreover, CBF decreases by ~22.5 and 47.5% in both regions (12.0 ± 5.6 vs. 9.3 ± 4.2 vs. 6.3 ± 3.2 ml/100 g). Due to the absence of vascularization in these regions, the MTT values show an increase of 26.1 and 111.5% in the penumbra and core (9.6 ± 2.7 vs. 12.1 ± 3.9 vs. 20.3 ± 6.9 s). We found that SR and PSR have also the capacity to accurately differentiate the core from the other brain regions: CL (−3.3 ± 1.1%, 64.3 ± 11.1%), penumbra region (−2.3 ± 2.6%, 81.7 ± 18.1%), and lesion core (1.1 ± 3.8%, 113.5 ± 46.1%).
Parameters of tumor perfusion were validated in a preclinical GBM model at 10 and 20 days after the surgery known to produce different levels of vascularization (Figure 4C). We can appreciate that it is possible to accurately identify core and peripheral tumor regions at both time points. The CBF values in the tumor core were reduced at 10 and 20 days to 74.5 and 48.8%, respectively. At 10 days, animals showed increased CBV and MTT values in the core compared to the tumor rim (3.4 ± 2.1 vs. 1.3 ± 1.2 ml/100 g, and 16.3 ± 8.5 vs. 2.2 ± 0.9 s). However, at 20 days of CBV and MTT, a change in behavior is detected (1.2 ± 0.4 vs. 4.1 ± 3 ml/100 g, and 0.67 ± 0.1 vs. 1 ± 0.43 s). These results may reflect the heterogeneous structure of the tumor regions, and the time evolution probably also contributes to this heterogeneity. As previously described, SR and PSR maps provided a GBM spatial distribution and add valuable diagnostic information. At 20 days, we found SR core 59.1 ± 6.7% vs. tumor rim 20.4 ± 3.4%, and PSR core −249.4 ± 181.1% vs. tumor rim 253.5 ± 27.3%.
3.2. Bland–altman analysis
Figure 5 shows the BA plots for our data and healthy SD rats, corresponding to the left and right cortex. We observed low mean differences for CBV (2.5 ml/100 g) and MTT (2.7 s) variables and large differences for CBF [49.9 ml/(100 g*min)]. Regarding the ischemic stroke model, we compared the lesion core and CL cortex, for our results and literature data (Figure 6). The BA plot highlighted low mean values of bias for CBF [2.6 ml/(100 g*min)], CBV (1.17 ml/100 g), and MTT (9.7 s) for the core lesion. We found the same for CBV (2.8 ml/100 g) and MTT (3 s) in the CL cortex, but not for CBF [65.2 ml/(100 g*min)]. Regional BA analysis of the tumor core and CL cortex is illustrated in Figure 7. We observed low mean discrepancies for CBV (−3.45 and 6.63 ml/100 g) and MTT (−4.71 and 0.73 s) variables, and large discrepancies for CBF [30.9 and 32.2 ml/(100 g*min)] both in the tumor core and in the CL cortex.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g005.jpg)
Figure 5. Bland–Altman plots for CBV, CBF, and MTT for our results vs. literature data regarding SD healthy rats (Adam et al., 2003; Perles-Barbacaru and Lahrech, 2007; Rouine et al., 2013; Lee et al., 2021). Dashed lines represent the bias, +95% (upper line) and −95% (lower line) of the limits of agreements.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g006.jpg)
Figure 6. Bland–Altman plots for CBV, CBF, and MTT for our results vs. literature data regarding the ischemic stroke animal model (Thomas et al., 2006; Zhang et al., 2014; Boisserand et al., 2017; Livingston et al., 2020; Tsai et al., 2021). Dashed lines represent the bias, +95% (upper line) and −95% (lower line) of the limits of agreements.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1202156/fninf-17-1202156-HTML-r1/image_m/fninf-17-1202156-g007.jpg)
Figure 7. Bland–Altman plots for CBV, CBF, and MTT for our results vs. literature data regarding the glioblastoma animal model (García-Palmero et al., 2013; Stokes et al., 2014; Gonawala et al., 2018; Clément et al., 2021). Dashed lines represent the bias, +95% (upper line) and −95% (lower line) of the limits of agreements.
4. Discussion
The absolute quantification of perfusion-weighted parameters is useful for preclinical research studies of different neurological diseases. The integrity of the BBB is critical to normal brain function, and the perfusion variables allow the detection of brain flow dysregulation areas as the origin of future brain events. Recently, cerebrovascular dysfunction has emerged as an early manifestation of Alzheimer's disease with a direct impact on the evolution and clinical expression of dementia. It has been shown that mice expressing mutated tau exhibit a selective suppression of neural activity-induced CBF increases that precedes tau pathology and cognitive impairment (Warpechowski et al., 2023). Regarding stroke, these parametric maps help to diagnose and predict the final cerebral ischemic stroke volume identifying core and penumbra regions (Demeestere et al., 2020). In addition, it is also a valuable diagnostic tool for tumor grading that may contribute to a better understanding of tumor evolution and characterize the vascularization of the tumor rim and core (Aprile et al., 2015; Huhndorf et al., 2016).
Regarding preclinical MRI data processing, there is variability in the programs for computing the perfusion variables that use different methods for the analysis, both regarding the mathematical modeling and the protocols of image acquisition. Moreover, most of them have closed access, limited to the institution responsible for the software development (Gordaliza et al., 2015; Huhndorf et al., 2016; López-Larrubia, 2018; Hartmann et al., 2020; Tsai et al., 2021). This issue is present even in human health, causing concerns about the accuracy of software quantitative perfusion parameters. Therefore, the obtained hemodynamics parameters cannot be easily compared between different studies, and perfusion-weighted methods remain an active area of research.
In this study, the DSC quantification tool designed performs the basic mathematical steps to generate the main hemodynamic parametric maps CBV, CBF, MTT, including PSR, and SR solving the indicator dilution model through deconvolution-based methods, being implemented in both TSVD and Tikhonov regularization. Furthermore, it contains an additional processing step, fitting the data to a gamma function in order to remove the effects of both noise and contrast recirculation. The tool works as Python script and has been published under a free software license with a modular architecture to allow own improvements or new pre-processing algorithms (filters, alignments, or image corrections). The parametric maps were obtained with an acceptable signal-to-noise ratio (SNR) to identify the small blood signal with a spatial resolution of 0.172 × 0.172 × 1.5 mm3/pixel and can be stored using standard image formats to carry out the posterior analysis. In general, SR and PSR maps appear noisier. However, these parameters can be rapidly acquired, and they do not require complex mathematical post-processing steps and help to characterize, e.g., tumor vasculature between the core and periphery over time.
It is important to note that the absolute quantification of hemodynamic parameters requires the capability of detecting signal variations in or adjacent to a large vessel, usually the paraclinoid or middle cerebral arteries. However, due to the low spatial resolution, image-based determination of AIF in rodents is challenging and very difficult to achieve. This limitation usually overcome using different techniques, though it may be useful to perform a semiquantitative quantification using laboratory arbitrary units (a.u.) or normalization of the obtained data. There are also more complex techniques that automatically determine AIF based on, among other things, the detection of pixels with the greatest contrast enhancement in various regions of interest, the use of clustering techniques, or the creation of physical models for the echo-planar signal intensity from an artery (Ostergaard et al., 1996; Wu et al., 2003; López-Larrubia, 2018). In this study, we provide the possibility of semi-automatic AIF or VOF adjustment with real curve monitoring for an accurate AIF-pixel selection.
To the extent of our knowledge, there is no standard recommendation software for preclinical perfusion with a universally accepted rating. Thus, for the validation and examination of our tool, we used diverse data sets from different rat models of brain diseases: i) temporal BBB dysregulation animals, ii) cerebral CCH model, iii) ischemic stroke, and vi) GBM. In general, our software yielded excellent results in terms of agreements for the expected brain regional analysis based on the animal model addressed.
It has been reported that in the cerebral chronic hypoperfusion model 2VO, 2 weeks after the surgery, the CBF recovers to 55–65% of the control level (Cao et al., 2018). Our results showed 51% of CBF getting back 3 weeks after the surgery. Regarding the mannitol group, a previous study described that CBF increased predominantly in the hemisphere in which mannitol was injected. In our study, an alternative experimental model was employed, mannitol was i.v. injected in the tail as a recognized and reliable procedure to temporal modulate the BBB permeability. We found that CBF decreases by ~35.6% related to healthy animals, while CBV and MTT values remained similar or increased at 10 min after mannitol injection (Duong et al., 2000; Tanaka et al., 2011). As has been previously described (Robertson et al., 2011; Reid et al., 2012; Cipolla et al., 2017), ischemic core and penumbra areas were calculated based on a CBF reduction compared with the equivalent CL hemisphere ROI (36.4 and 77.2%). From these regional definitions, the CBV and MTT maps obtained are in agreement with the expected metabolic values for these regions due to the limited vascularization. The preclinical tumor model at 10 and 20 days after the surgery allowed us to compare the time level vascularization. Previous studies determined that 4–9 days after tumor detection, the capillary permeability significantly increases in the core (Stadlbauer et al., 2015; Choi et al., 2016; Huhndorf et al., 2016). This increase in capillary permeability most likely reflects the beginning of necrosis within the tumor core that will continue to evolve over time. The regional analysis obtained in our study reveals that at 10 days, areas of increased CBV are mostly found at the core (3.4 ± 2.1 vs. 1.3 ± 1.2 ml/100 g). However, we can appreciate that this behavior is reversed after 20 days (1.2 ± 0.4 vs. 4.1 ± 3 ml/100 g), which reflects that the nucleus has necrotized and there is a vascular proliferation at the tumor rim.
The BA plot was also used to describe the agreement between our perfusion quantitative analyses and literature data regarding healthy rats, stroke, and GBM models. In general, the agreement for CBV and MTT is higher than for CBF. Possible explanations are the high CBF sensitivity errors due to the localization and delineation of ROIs, and the microvasculature mechanism effects of the mathematical tracer kinetic models. The higher sensitivity of DSC-MRI to susceptibility artifacts, an effect known and already studied (Maral et al., 2020), could also be another reason, being supported by the higher sensitivity to artifacts observed in the CBV and MTT maps, compared to the corresponding CBF map.
Finally, Table 1 reports a literature overview of absolute perfusion values obtained with different procedures, experimental conditions (including several rat strains), and different post-processing methods. Although it is difficult to compare, all hemodynamic parameters reported from the literature are in the same range as our results. In general, few studies assess CBF, CBV, and MTT simultaneously, and many studies show high deviations probably as a result of the regional analysis developed. As it can be appreciated, arterial spin labeling (ASL) MR perfusion is another technique commonly used that does not require intravenous administration of contrast. By using arterial blood water protons that have been magnetically labeled as endogenous tracers, this non-invasive and non-ionizing MRI technique assesses tissue perfusion (blood flow). Different techniques have been described to achieve ASL perfusion: i) pulsed (PASL), ii) continuous (CASL), iii) pseudo-continuous (PCASL), and iv) velocity-selective ASL (VS-ASL). However, the parameter most commonly derived in these protocols is CBF because CBV and MTT are difficult to reliably obtain, and this method usually presents low SNR and long acquisition time.
5. Conclusion
In order to facilitate the use and comparison of perfusion-weighted imaging in preclinical studies, we provide an open-source DSC post-processing software package. This software allows the calculation of several key quantitative perfusion parameters, such as CBF, CBV, and MTT, including SR and PSR maps from semi-automatic AIF or VOF adjustment. The open-source computational steps allow for possible improvements or pre-processing algorithms for other situations. The results obtained in diverse data sets of brain disease models in rats are consistent and in good agreement with values and behavior reported in the literature.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by 15011/2021/002, 15011/2021/003, 15011/2023/002, and 15011/2021/001.
Author contributions
JC, RI-R, PH, and MA-A: organization and design of the study. JC, RI-R, PH, MA-A, SF-R, and GF-C: manuscript drafting. AO, AS-V, MB-B, MP-M, AF, JP, and AM: supervision, review, and critique. MP-M, EL-A, SF-R, and AS-V: data acquisition. All authors read, reviewed, and agreed upon the manuscript version.
Funding
This research was funded by the Spanish Ministry of Science and Innovation (SAF2017-84267-R), PDC2021-121455-I00, Xunta de Galicia (Axencia Galega de Innovación: IN607A2022-03), Instituto de Salud Carlos III (ISCIII) (PI17/01103, ISCIII/PI21/01256/Co-financed by the European Union), Spanish Research Network on Cerebrovascular Diseases RETICS-INVICTUS PLUS (RD16/0019/0001), and RICORS-ICTUS (Cerebrovascular diseases) D21/0006/0003. MB-B is a PFIS Researcher (FI22/00200) of Instituto de Salud Carlos III. MP-M is a Sara Borrell Researcher (CD19/00033) of Instituto de Salud Carlos III. RI-R (CP22/00061) from the Miguel Servet Program of Instituto de Salud Carlos III and Co-financed by the EU. Sponsors did not participate in the study design, collection, analysis, or interpretation of the data, or in writing the report.
Conflict of interest
AF was employed by company Nasasbiotech.
The remaining 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.
References
Adam, J. F., Elleaume, H., Le Duc, G., Corde, S., Charvet, A. M., Troprès, I., et al. (2003). Absolute cerebral blood volume and blood flow measurements based on synchrotron radiation quantitative computed tomography. J. Cereb. Blood Flow Metab. 23, 499–512. doi: 10.1097/01.WCB.0000050063.57184.3C
Aprile, I., Giovannelli, G., Fiaschini, P., Muti, M., Kouleridou, A., Caputo, N., et al. (2015). High- and low-grade glioma differentiation: the role of percentage signal recovery evaluation in MR dynamic susceptibility contrast imaging. Radiol. Med. 120, 967–74. doi: 10.1007/s11547-015-0511-7
Baskerville, T. A., McCabe, C., Weir, C. J., Macrae, I. M., and Holmes, W. M. (2012). Noninvasive MRI measurement of CBF: evaluating an arterial spin labelling sequence with 99mTc-HMPAO CBF autoradiography in a rat stroke model. J. Cereb. Blood Flow Metab. 32, 973–7. doi: 10.1038/jcbfm.2012.19
Boisserand, L. S. B., Lemasson, B., Hirschler, L., Moisan, A., Hubert, V., Barbier, E. L., et al. (2017). Multiparametric magnetic resonance imaging including oxygenation mapping of experimental ischaemic stroke. J. Cereb. Blood Flow Metab. 37, 2196–2207. doi: 10.1177/0271678X16662044
Boxerman, J. L., Quarles, C. C., Hu, L. S., Erickson, B. J., Gerstner, E. R., Smits, M., et al. (2020). Jumpstarting brain tumor drug development coalition imaging standardization steering committee. Consensus recommendations for a dynamic susceptibility contrast MRI protocol for use in high-grade gliomas. Neuro. Oncol. 22, 1262–1275. doi: 10.1093/neuonc/noaa141
Bulin, A. L., Adam, J. F., and Elleaume, H. (2022). Stereotaxic implantation of F98 cells in fischer rats: a syngeneic model to investigate photodynamic therapy response in glioma. Methods Mol. Biol. 2451, 203–210. doi: 10.1007/978-1-0716-2099-1_15
Calamante, F., Gadian, D. G., and Connelly, A. (2000). Delay and dispersion effects in dynamic susceptibility contrast MRI: simulations using singular value decomposition. Magn. Reson. Med. 44, 466–473. doi: 10.1002/1522-2594(200009)44:3<466::aid-mrm18>3.0.co;2-m
Calamante, F., Gadian, D. G., and Connelly, A. (2002). Quantification of perfusion using bolus tracking magnetic resonance imaging in stroke: assumptions, limitations, and potential implications for clinical use. Stroke. 33, 1146–1151. doi: 10.1161/01.STR.0000014208.05597.33
Calamante, F., Thomas, D. L., Pell, G. S., Wiersma, J., and Turner, R. (1999). Measuring cerebral blood flow using magnetic resonance imaging techniques. J. Cereb. Blood Flow Metab. 19, 701–735. doi: 10.1097/00004647-199907000-00001
Callewaert, B., Jones, E. A. V., Himmelreich, U., and Gsell, W. (2021). Non-invasive evaluation of cerebral microvasculature using pre-clinical MRI: principles, advantages and limitations. Diagnostics. 11, 926. doi: 10.3390/diagnostics11060926
Cao, D., Bai, Y., and Li, L. (2018). Common carotid arteries occlusion surgery in adult rats as a model of chronic cerebral hypoperfusion. Bio. Protoc. 8, e2704. doi: 10.21769/BioProtoc.2704
Choi, S. H., Jung, S. C., Kim, K. W., Lee, J. Y., Choi, Y., Park, S. H., et al. (2016). Perfusion MRI as the predictive/prognostic and pharmacodynamic biomarkers in recurrent malignant glioma treated with bevacizumab: a systematic review and a time-to-event meta-analysis. J. Neurooncol. 128, 185–194. doi: 10.1007/s11060-016-2102-4
Cipolla, M. J., Sweet, J. G., and Chan, S. L. (2017). Effect of hypertension and peroxynitrite decomposition with FeTMPyP on CBF and stroke outcome. J. Cereb. Blood Flow Metab. 37, 1276–1285. doi: 10.1177/0271678X16654158
Clément, A., Doyen, M., Fauvelle, F., Hossu, G., Chen, B., Barberi-Heyob, M., et al. (2021). In vivo characterization of physiological and metabolic changes related to isocitrate dehydrogenase 1 mutation expcression by multiparametric MRI and MRS in a rat model with orthotopically grafted human-derived glioblastoma cell lines. NMR Biomed. 34, e4490. doi: 10.1002/nbm.4490
Demeestere, J., Wouters, A., Christensen, S., Lemmens, R., and Lansberg, M. G. (2020). Review of perfusion imaging in acute ischemic stroke: from time to tissue. Stroke. 51, 1017–1024. doi: 10.1161/STROKEAHA.119.028337
D'haeseleer, M., Hostenbach, S., Peeters, I., Sankari, S. E., Nagels, G., De Keyser, J., et al. (2015). Cerebral hypoperfusion: a new pathophysiologic concept in multiple sclerosis? J. Cereb. Blood Flow Metab. 35, 1406–10. doi: 10.1038/jcbfm.2015.131
Duong, T. Q., Silva, A. C., Lee, S. P., and Kim, S. G. (2000). Functional MRI of calcium-dependent synaptic activity: cross correlation with CBF and BOLD measurements. Magn. Reson. Med. 43, 383–392. doi: 10.1002/(sici)1522-2594(200003)43:3<383::aid-mrm10>3.0.co;2-q
Fieselmann, A., Kowarschik, M., Ganguly, A., Hornegger, J., and Fahrig, R. (2011). Deconvolution-based CT and MR brain perfusion measurement: theoretical model revisited and practical implementation details. Int. J. Biomed. Imaging. 2011, 467563. doi: 10.1155/2011/467563
García-Palmero, I., López-Larrubia, P., Cerdán, S., and Villalobo, A. (2013). Nuclear magnetic resonance imaging of tumour growth and neovasculature performance in vivo reveals Grb7 as a novel antiangiogenic target. NMR Biomed. 26, 1059–1069. doi: 10.1002/nbm.2918
Giavarina, D. (2015). Understanding bland altman analysis. Biochem. Med. 25, 141–151. doi: 10.11613/BM.2015.015
Gonawala, S., Aryal, M., Ewing, J. R., deCarvalho, A. C., Kalkanis, S., Ali, M. M. M. R. I., et al. (2018). Monitoring of cerebral blood flow after the delivery of nanocombretastatin across the blood brain tumor barrier. J. Nanomed. Nanotechnol. 9, 516. doi: 10.4172/2157-7439.1000516
Gordaliza, P. M., Mateos-Pérez, J. M., Montesinos, P., Guzmán-de-Villoria, J. A., Desco, M., and Vaquero, J. J. (2015). Development and validation of an open source quantification tool for DSC-MRI studies. Comput. Biol. Med. 58, 56–62. doi: 10.1016/j.compbiomed.2015.01.002
Harris, C. R., Millman, K. J., van der Walt, S. J., Gommers, R., Virtanen, P., Cournapeau, D., et al. (2020). Array programming with NumPy. Nature. 585, 357–362. doi: 10.1038/s41586-020-2649-2
Hartmann, A., Scholz, V. B., and Lautenschläger, I. E. (2020). MRI perfusion analysis using freeware, standard imaging software. BMC Vet. Res. 16, 141. doi: 10.1186/s12917-020-02352-0
Haselhorst, R., Kappos, L., Bilecen, D., Scheffler, K., Möri, D., Rad,ü, E. W., et al. (2000). Dynamic susceptibility contrast MR imaging of plaque development in multiple sclerosis: application of an extended blood-brain barrier leakage correction. J. Magn. Reson. Imaging. 11, 495–505. doi: 10.1002/(sici)1522-2586(200005)11:5<495::aid-jmri5>3.0.co;2-s
Huhndorf, M., Moussavi, A., Kramann, N., Will, O., Hattermann, K., Stadelmann, C., et al. (2016). Alterations of the blood-brain barrier and regional perfusion in tumor development: MRI insights from a rat C6 glioma model. PLoS ONE. 11, e0168174. doi: 10.1371/journal.pone.0168174
Hunter, J. D. (2007). Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90–95. doi: 10.1109/MCSE.2007.55
Keston, P., Murray, A. D., and Jackson, A. (2003). Cerebral perfusion imaging using contrast-enhanced MRI. Clin. Radiol. 58, 505–513. doi: 10.1016/S0009-9260(03)00130-2
Konstas, A. A., Goldmakher, G. V., Lee, T. Y., and Lev, M. H. (2009). Theoretic basis and technical implementations of CT perfusion in acute ischemic stroke, part 1: Theoretic basis. AJNR Am. J. Neuroradiol. 30, 662–8. doi: 10.3174/ajnr.A1487
Lacalle-Aurioles, M., Mateos-Pérez, J. M., Guzmán-De-Villoria, J. A., Olazarán, J., Cruz-Orduña, I., Alemán-Gómez, Y., et al. (2014). Cerebral blood flow is an earlier indicator of perfusion abnormalities than cerebral blood volume in Alzheimer's disease. J. Cereb. Blood Flow Metab. 34, 654–9. doi: 10.1038/jcbfm.2013.241
Larkin, J. R., Simard, M. A., Khrapitchev, A. A., Meakin, J. A., Okell, T. W., Craig, M., et al. (2019). Quantitative blood flow measurement in rat brain with multiphase arterial spin labelling magnetic resonance imaging. J. Cereb. Blood Flow Metab. 39, 1557–1569. doi: 10.1177/0271678X18756218
Lee, K., Bohnert, S., Vair, C., Mikler, J., and Dunn, J. F. (2021). Cerebral blood flow and oxygenation in rat brain after soman exposure. Toxicol. Lett. 336, 50–56. doi: 10.1016/j.toxlet.2020.10.009
Livingston, J. M., McDonald, M. W., Gagnon, T., Jeffers, M. S., Gomez-Smith, M., Antonescu, S., et al. (2020). Influence of metabolic syndrome on cerebral perfusion and cognition. Neurobiol. Dis. 137, 104756. doi: 10.1016/j.nbd.2020.104756
López-Larrubia, P. (2018). Dynamic susceptibility contrast MRI in small animals. Methods Mol. Biol. 1718, 41–57. doi: 10.1007/978-1-4939-7531-0_3
Maral, H., Ertekin, E., Tunçyürek, Ö., and Özsunar, Y. (2020). Effects of susceptibility artifacts on perfusion MRI in patients with primary brain tumor: a comparison of arterial spin-labeling versus DSC. AJNR Am. J. Neuroradiol. 41, 255–261. doi: 10.3174/ajnr.A6384
Mason, D. (2011). SU-E-T-33: pydicom: an open source DICOM library. Med. Phys. 38, 3493–3493. doi: 10.1118/1.3611983
Murray, A., Kemenade, H., van wiredfool, Clark (Alex), J. A., Karpinsky, A., Baranovič, O., et al. (2023). python-pillow/Pillow: 9, 4.0.
Østergaard, L. (2004). Cerebral perfusion imaging by bolus tracking. Top Magn. Reson. Imaging. 15, 3–9. doi: 10.1097/00002142-200402000-00002
Ostergaard, L., Weisskoff, R. M., Chesler, D. A., Gyldensted, C., and Rosen, B. R. (1996). High resolution measurement of cerebral blood flow using intravascular tracer bolus passages. Part I: mathematical approach and statistical analysis. Magn. Reson. Med. 36, 715–725. doi: 10.1002/mrm.1910360510
Park, L., Hochrainer, K., Hattori, Y., Ahn, S. J., Anfray, A., Wang, G., et al. (2020). Tau induces PSD95-neuronal NOS uncoupling and neurovascular dysfunction independent of neurodegeneration. Nat. Neurosci. 23, 1079–1089. doi: 10.1038/s41593-020-0686-7
Perez, F. (2007). IPython: a system for interactive scientific computing. Comput. Sci. Eng. 9, 21–29. doi: 10.1109/MCSE.2007.53
Perles-Barbacaru, A. T., and Lahrech, H. (2007). A new magnetic resonance imaging method for mapping the cerebral blood volume fraction: the rapid steady-state T1 method. J. Cereb. Blood Flow Metab. 27, 618–31. doi: 10.1038/sj.jcbfm.9600366
Qi, Q., Fox, M. S., Lim, H., Bartha, R., Scholl, T. J., Hoffman, L., et al. (2021). Multimodality In Vivo imaging of perfusion and glycolysis in a rat model of C6 glioma. Mol. Imaging Biol. 23, 516–526. doi: 10.1007/s11307-021-01585-1
Reid, E., Graham, D., Lopez-Gonzalez, M. R., Holmes, W. M., Macrae, I. M., McCabe, C., et al. (2012). Penumbra detection using PWI/DWI mismatch MRI in a rat stroke model with and without comorbidity: comparison of methods. J. Cereb. Blood Flow Metab. 32, 1765–77. doi: 10.1038/jcbfm.2012.69
Robertson, C. A., McCabe, C., Gallagher, L., Lopez-Gonzalez Mdel, R., Holmes, W. M., Condon, B., et al. (2011). Stroke penumbra defined by an MRI-based oxygen challenge technique: 2. Validation based on the consequences of reperfusion. J. Cereb. Blood Flow Metab. 31, 1788–98. doi: 10.1038/jcbfm.2011.67
Rouine, J., Gobbo, O. L., Campbell, M., Gigliucci, V., Ogden, I., McHugh Smith, K., et al. (2013). 'ecstasy' increases cerebral cortical perfusion determined by bolus-tracking arterial spin labelling (btASL) MRI. Br. J. Pharmacol. 169, 974–87. doi: 10.1111/bph.12178
Stadlbauer, A., Pichler, P., Karl, M., Brandner, S., Lerch, C., Renner, B., et al. (2015). Quantification of serial changes in cerebral blood volume and metabolism in patients with recurrent glioblastoma undergoing antiangiogenic therapy. Eur. J. Radiol. 84, 1128–36. doi: 10.1016/j.ejrad.2015.02.025
Stokes, A. M., Skinner, J. T., and Quarles, C. C. (2014). Assessment of a combined spin- and gradient-echo (SAGE) DSC-MRI method for preclinical neuroimaging. Magn. Reson. Imaging. 32, 1181–90. doi: 10.1016/j.mri.2014.08.027
Tanaka, Y., Nagaoka, T., Nair, G., Ohno, K., and Duong, T. Q. (2011). Arterial spin labeling and dynamic susceptibility contrast CBF MRI in postischemic hyperperfusion, hypercapnia, and after mannitol injection. J. Cereb. Blood Flow Metab. 31, 1403–11. doi: 10.1038/jcbfm.2010.228
Thomas, D. L., Lythgoe, M. F., van der Weerd, L., Ordidge, R. J., and Gadian, D. G. (2006). Regional variation of cerebral blood flow and arterial transit time in the normal and hypoperfused rat brain measured using continuous arterial spin labeling MRI. J. Cereb. Blood Flow Metab. 26, 274–82. doi: 10.1038/sj.jcbfm.9600185
Tsai, Y. H., Su, C. H., Lee, I. N., Yang, J. T., Lin, L. C., Huang, Y. C., et al. (2021). Effects of early rehydration on brain perfusion and infarct core after middle cerebral artery occlusion in rats. Brain Sci. 11, 439. doi: 10.3390/brainsci11040439
Vieites-Prado, A., Iglesias-Rey, R., Fernández-Susavila, H., da Silva-Candal, A., Rodríguez-Castro, E., Gröhn, O. H., et al. (2016). Protective Effects and magnetic resonance imaging temperature mapping of systemic and focal hypothermia in cerebral ischemia. Stroke. 47, 2386–96. doi: 10.1161/STROKEAHA.116.014067
Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., et al. (2020). SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods. 17, 261–272. doi: 10.1038/s41592-019-0686-2
Warpechowski, M., Warpechowski, J., Kulczyńska-Przybik, A., and Mroczko, B. (2023). Biomarkers of activity-dependent plasticity and persistent enhancement of synaptic transmission in alzheimer disease: a review of the current status. Med. Sci. Monit. 29, e938826. doi: 10.12659/MSM.938826
Wu, O., Østergaard, L., Weisskoff, R. M., Benner, T., Rosen, B. R., Sorensen, A. G., et al. (2003). Tracer arrival timing-insensitive technique for estimating flow in MR perfusion-weighted imaging using singular value decomposition with a block-circulant deconvolution matrix. Magn. Reson. Med. 50, 164–174. doi: 10.1002/mrm.10522
Keywords: DSC-MRI imaging, glioblastoma (GBM), neuroimaging, perfusion analysis, Python, stroke
Citation: Fernández-Rodicio S, Ferro-Costas G, Sampedro-Viana A, Bazarra-Barreiros M, Ferreirós A, López-Arias E, Pérez-Mato M, Ouro A, Pumar JM, Mosqueira AJ, Alonso-Alonso ML, Castillo J, Hervella P and Iglesias-Rey R (2023) Perfusion-weighted software written in Python for DSC-MRI analysis. Front. Neuroinform. 17:1202156. doi: 10.3389/fninf.2023.1202156
Received: 07 April 2023; Accepted: 27 June 2023;
Published: 01 August 2023.
Edited by:
Xiaohao Cai, University of Southampton, United KingdomReviewed by:
Ivana Galinovic, Charité University Medicine Berlin, GermanyMalgorzata Burek, Julius Maximilian University of Würzburg, Germany
Copyright © 2023 Fernández-Rodicio, Ferro-Costas, Sampedro-Viana, Bazarra-Barreiros, Ferreirós, López-Arias, Pérez-Mato, Ouro, Pumar, Mosqueira, Alonso-Alonso, Castillo, Hervella and Iglesias-Rey. 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: Ramón Iglesias-Rey, ramon.iglesias.rey@sergas.es; Pablo Hervella, pablo.hervella.lorenzo@sergas.es
†These authors have contributed equally to this work