Abstract
Accurate vascular input function (VIF) derivation is essential in brain dynamic contrast-enhanced (DCE) MRI. The optimum site for VIF estimation is, however, debated. This study sought to compare VIFs extracted from the internal carotid artery (ICA) and its branches with an arrival-corrected vascular output function (VOF) derived from the superior sagittal sinus (VOFSSS). DCE-MRI datasets from sixty-six patients with different brain tumours were retrospectively analysed and plasma gadolinium-based contrast agent (GBCA) concentration-time curves used to extract VOF/VIFs from the SSS, the ICA, and the middle cerebral artery. Semi-quantitative parameters across each first-pass VOF/VIF were compared and the relationship between these parameters and GBCA dose was evaluated. Through a test–retest study in 12 patients, the repeatability of each semiquantitative VOF/VIF parameter was evaluated; and through comparison with histopathological data the accuracy of kinetic parameter estimates derived using each VOF/VIF and the extended Tofts model was also assessed. VOFSSS provided a superior surrogate global input function compared to arteries, with greater contrast-to-noise (p < 0.001), higher peak (p < 0.001, repeated-measures ANOVA), and a greater sensitivity to interindividual plasma GBCA concentration. The repeatability of VOFSSS derived semi-quantitative parameters was good to excellent (ICC = 0.717–0.888) outperforming arterial based approaches. In contrast to arterial VIFs, kinetic parameters obtained using a SSS derived VOF permitted detection of intertumoural differences in both microvessel surface area and cell density within resected tissue specimens. These results support the usage of an arrival-corrected VOFSSS as a surrogate vascular input function for kinetic parameter mapping in brain DCE-MRI.
Similar content being viewed by others
Introduction
Dynamic contrast-enhanced (DCE) MRI has a developing role as an imaging tool for quantifying brain tumour microvasculature and tumour response to anti-angiogenic therapy1,2,3,4,5,6. In human clinical studies, non-invasive measurement of a suitable vascular input function (VIF) is essential for deriving microvascular kinetic parameters from brain DCE-MRI. Non-invasive measurement of a VIF with high temporal resolution is especially important when first-pass bolus tracking is necessary for quantitative kinetic analysis, such as when using the extended Tofts model7,8,9,10. Use of a fixed, experimentally derived population-averaged VIF for all subjects has been previously proposed11,12,13, thereby simplifying data acquisition, but large variations can occur in the actual VIF between subjects and scan visits due to both technical (e.g. differences in injection timing and dose) and patient specific factors (e.g. cardiac output, haematocrit, caffeine intake and atherosclerosis related vessel narrowing)11,14,15. For this reason methods have also been developed that enable simultaneous measurement of plasma gadolinium-based contrast agent (GBCA) concentration changes in both the blood and tissue under study, permitting VIF to be measured on an individual patient basis11,14,16.
In DCE-MRI studies of brain tumours the ideal choice for defining an individual VIF is the feeding artery of the tumour but due to either data acquisition constraints, the small size of the feeding vessel or lack of the feeding artery within the imaging field of view (FOV) this is often not possible11,16,17,18. Large intracranial arteries such as the internal carotid artery (ICA) and middle cerebral artery (MCA) are therefore often used as a surrogate global VIF measurement11,16,19. Use of the superior sagittal sinus (SSS), a venous structure, as a surrogate global input function has also been adopted, however, in numerous cross-sectional and longitudinal DCE-MRI studies1,2,5,11,16,19,20,21,22,23,24,25,26. The presented rationales for the use of the SSS vascular output function (VOF) as a surrogate global input function were said to be an assumption that the venous and arterial GBCA concentration is equal; that the degree of dispersion between arterial and venous structures in the brain during the first-pass circulation of the GBCA bolus is minimal; and that due to its larger comparative size the SSS is less susceptible to partial volume errors (PVE) than smaller intracranial arteries such as the ICA8.
The most commonly used 3D DCE-MRI acquisition is the 3D spoiled gradient echo method, and a VIF (plasma GBCA concentration-time curve) is typically determined from GBCA related magnitude changes in signal intensity27. Using a whole-brain, low-dose high temporal resolution (LDHT) axial 3D spoiled gradient-recalled echo sequence and a magnitude-based determination method we simultaneously measured a VOF/VIF from the SSS and large intracranial arteries (ICA, MCA) in a brain tumour patient cohort. Through this we sought to compare and understand features of each VOF/VIF and evaluate the respective ability of each VOF/VIF to accurately capture interindividual changes in patient dosing and plasma GBCA concentration. Through an included test–retest study we sought to establish the respective reproducibility of parameters derived from these arterial and venous VOF/VIFs; and through comparison with resected tumour specimens in a patient cohort, evaluate the ability of kinetic parameters derived using each VOF/VIF to detect intertumoural differences in histopathological data.
Methods
Study population
Previously acquired dual temporal resolution (DTR), dual injection DCE-MRI data in three groups of patients were analysed for this study: twenty-five patients with newly diagnosed WHO grade IV glioma synonym glioblastoma (GBM); twenty-nine patients with sporadic vestibular schwannoma (VS) listed for either radiological surveillance or treatment with surgery or stereotactic radiosurgery (SRS); and twelve patients with neurofibromatosis type 2 (NF2) related VS undergoing treatment with the anti-vascular endothelial growth factor (anti-VEGF) antibody, bevacizumab (Avastin ©). Ethical approvals were in place from the National Research Ethics Service Greater Manchester North-West research ethics committee (REC references: 13/NW/0131, 13/NW/0247 and 15/NW/0429). All patients had provided informed consent for study participation and later analysis of their MRI data and all research was performed in accordance with the Declaration of Helsinki and with local guidelines and policies.
MR imaging
The 25 patients with GBM and 29 patients with sporadic VS were all imaged once on a 1.5 T scanner (Philips Achieva, Best, Netherlands). Thirteen of the included patients with GBM and twenty-two of the included patients with sporadic VS had been recruited and scanned as part of previous published studies at our institution1,2,6. The twelve patients with NF2-related VS had similarly been recruited as part of an earlier published study investigating bevacizumab (Avastin ©) related changes in DCE-MRI derived kinetic parameters in VS and these patients had been imaged twice at 1.5 T: pre-treatment (day 0) and 3 months (day 90) following bevacizumab (Avastin ©) treatment.
DCE-MRI data was acquired using a previously described DTR, dual injection technique6,16,28. Single dose macrocyclic GBCA (gadoterate meglumine; Dotarem, Guerbet S.A.) was used at a dose of 0.2 ml/kg. For VOF/VIF estimation and as the first part of this DTR technique, a low-dose fixed volume pre-bolus (either 2 or 3mls) of GBCA was administered over 1 s during acquisition of a high temporal resolution (LDHT) DCE-MRI dataset. All intravenous injections were performed using a two-cylinder power injector (MEDRAD® Spectris Solaris EP, Bayer, PA, US). The GBCA and 0.9% saline are contained within separate cylinders and the pre-bolus injection was followed by a chaser of 20 ml of 0.9% saline administered at the same rate (2 or 3 ml/s). A 3D spoiled gradient recalled echo (GRE) sequence with axial slab orientation and anterior–posterior frequency encoding was used for data acquisition and acquisition parameters for this LDHT acquisition were as follows: flip angle of 20°, TR/TE of 2.5 ms/0.696 ms, SENSE acceleration factor of 1.8, reconstructed matrix size of 96 × 96 × 22, voxel size of 2.5 × 2.5 × 6.35 mm3, pixel bandwidth of 700 Hz, frame duration (Δt) 1.0 s (n = 300)6,16,28. The minimum TE and fixed volume low GBCA dose used for the LDHT DCE series and VOF/VIF estimation was designed to avoid signal magnitude saturation and expected to produce minimal T2* and water exchange effects29,30,31.
As the second part of this DTR DCE-MRI technique, a full-dose of GBCA (dose = 0.2 ml/kg ·weight – dose of pre-bolus) was administered at the same rate (2 to 3 ml/s) as the pre-bolus (followed by a chaser of 20 ml of 0.9% saline administered at the same rate) during acquisition of a high-spatial resolution (FDHS) sequence6,28. FFT (Fast Fourier Transform) reconstruction in the z-direction was used for both the low-dose high temporal resolution (LDHT) and full-dose high spatial resolution (FDHS) acquisitions, doubling the number of slices6. Variable flip-angle (VFA; α = 2°, 8°, 15° and 20°) acquisitions were undertaken prior to both the LDHT and FDHS DCE-MRI series for baseline longitudinal relaxation rate (R10) mapping, and the spatial resolution of each VFA acquisition series was chosen to match the LDHT and FDHS DCE series respectively6.
To eliminate unsaturated flowing spins entering the imaging slab and improve the accuracy of VIF estimation27,32 a large 3D acquisition volume covering the top of the brain, the circle of Willis and the terminations of the internal carotid arteries bilaterally was used16,33. Throughout the FOV, the number of radiofrequency (RF) pulses and the gradient spoiling that spins in flowing blood received was also maximized through the use of a fast spoiled GRE sequence with short TR and phase cycling. Gradient spoilers were applied along both the read and slice/slab selection directions and phase cycling with a phase increment angle of 117° was used. This allowed for more complete dephasing of residual transverse magnetization, minimizing blood inflow-induced errors within each imaging slice27,30,34. A pulse sequence diagram of the 3D spoiled gradient recalled echo (GRE) sequence used for both the LDHT and FDHS acquisition is shown in Fig. 1A.
Vascular input function extraction
For VOF/VIF extraction, acquired 4D LDHT DCE-MRI datasets (voxel size 2.5 × 2.5 × 6.35 mm3) were spatially aligned with and resliced to the FDHS data (voxel size 1 × 1 × 2 mm3) using SPM35. Although native LDHT datasets can also be used for input function extraction without prior co-registration and re-slicing, use of this co-registration step improved delineation of the blood vessel from surrounding tissues and permitted greater flexibility in ROI delineation and voxel selection. Plasma GBCA concentration-time curves, Cp(t), were derived from signal intensity (SI)—time curves measured at three different sites: 1) vertical/posterior third of the SSS, VIFSSS; 2) the ICA just distal to the carotid syphon, VIFICA: and 3) horizontal segment of the middle cerebral artery (MCA), VIFMCA. The first-pass data from these GBCA concentration-time curves were then fitted using a gamma variate function with a recirculation cut-off36, i.e.
For each extracted VOF/VIF the following semi-quantitative parameters were derived and compared: the contrast-to-noise ratio (CNR), the bolus arrival time (BAT), the bolus peak-amplitude, the bolus peak width (FWHM, full-width at half-maximum) and the bolus peak∙FWHM product (PWP), a parameter associated with the area under the bolus curve.
For measurement of VOF in the SSS, blood inflow-induced errors were reduced by including the anterior and middle third of the SSS in the field of view thereby maximizing the number of RF pulses that the flowing blood experiences before it reaches the posterior third of the SSS. Following manual delineation of a small rectangle ROI within the vertical/posterior SSS (Fig. 1B), an automatic extraction method was used to search and identify voxels within this segment of the SSS that display maximum enhancement area under the SI curve within 30 s of the bolus arrival time (AUC30)37. Through this semi-automatic extraction method voxels with maximal AUC30 and thereby less inflow and PVE were chosen for inclusion in the VOF. The above semi‐automatic extraction method was also applied for VIF extraction from the ICA and MCA. A pair of rectangle ROIs was manually drawn on an axial image section over the ICA (just distal to the carotid syphon) or MCAs bilaterally respectively (Fig. 1). Similar to VOFSSS an automatic method was then used to search and identify voxels within neighbouring contiguous axial slices that displayed maximum enhancement (AUC30).
For each VOF/VIF a mean SI-time curve was calculated from 20 voxels with the highest AUC30 and this mean SI-time curve was then converted to a plasma GBCA concentration-time curve Cp(t) using previously described methods33. Due to the difficulties of accurately measuring the pre-contrast T1 of flowing blood using standard DCE-MRI sequences and the bias introduced through in vivo experiments38,39,40, a literature value of blood R10 of 0.694 s-1 was used for the conversion41,42.
Kinetic parameter analysis
To evaluate the effect of VIF approach on kinetic parameter estimates, high spatial resolution (1 × 1 × 2 mm3) voxelwise maps of the microvascular kinetic parameters Ktrans (transfer constant), vp (fractional plasma volume) and ve (the fractional volume of extravascular extracellular space or EES) were derived using the extended Tofts model (ETM)43 and the previously described LEGATOS (LEvel and rescale the Gadolinium contrast concentrations curves of high-temporal TO high Spatial DCE-MRI) method6. Pre-surgery DCE-MRI datasets from 15 sporadic VS with available comparative tissue histology were chosen as test group for this analysis and for each VOF/VIF (SSS, ICA, MCA) separate kinetic parameter estimates were derived43. In addition to derivation of Ktrans, vp and ve within each tumour voxel through the ETM, the voxelwise intracellular fraction (vi) was also estimated through the relationship vi = 1−ve−vp.
The LEGATOS method for deriving high-spatial resolution kinetic parameter maps from DTR, dual-injection DCE-MRI data has been previously described6. In key step 1 of this method, errors through temporal jitter uncertainty are reduced through construction of a merged DTR 4D GBCA concentration volume containing a high temporal (HT) resolution ‘arterial’ phase followed by a later low temporal but high spatial (HS) resolution ‘parenchymal’ phase6,44. In key step II the high temporal but low spatial resolution arterial phase of each pixel concentration curve is then re-scaled using the LEGATOS method and a derived pixelwise calibration ratio, to increase the spatial resolution of derived kinetic parameter maps. For the LEGATOS method, a combined VOF/VIF is adopted. The whole Cp (t) from the LDHT-derived VOF/VIF is concatenated with the dose-calibrated late part of the Cp (t) measured from the FDHS-derived VOF/VIF, and is used for kinetic analysis of the LEGATOS-generated 4D high spatiotemporal resolution GBCA concentration volume6.
The BAT for each tissue voxel is calculated as part of each fitting procedure and the Cp(t) measured from each VOF/VIF time-shifted and aligned with the BAT of each tissue voxel contrast agent concentration-time curve6. As part of the fitting procedure and to assess the discrepancy between the original data and the derived curve a map of scaled fitting error (SFE) was also generated, with voxels displaying an SFE value > 50% being excluded from the tumour statistics6,45. For all patients, the SFE and derived kinetic parameter maps, both before and after exclusion of voxels with SFE > 50%, were visually inspected to confirm the acceptance of using SFE > 50% for outlier tumour voxel exclusion6.
Tissue analysis
For the 15 resected sporadic VS, previously obtained tissue metrics were compared against derived kinetic parameter estimates using each VOF/VIF6. Collected paraffin blocks from each case were cut into serial 5-µm tissue sections and assessed for cell density (haematoxylin and eosin, H&E), vascular permeability (fibrinogen) and microvessel surface area (CD31) using immunoperoxidase immunohistochemistry and established protocols1,2,6. Ethical approval was obtained for tissue analyses (REC reference 15/NW/0429 and 19/NS/0167) and detailed protocols are described in prior publications1,2,6.
Statistical analysis
The SPSS statistical software package (version 25, IBM Corp.) and Stata version 11 were used for all statistical tests. Extracted semi-quantitative parameters were compared across each vessel (SSS, ICA and MCA) using a repeated-measures ANOVA with Greenhouse–Geisser correction for non-sphericity. Post hoc analysis of pairwise comparisons between different VOF/VIF locations was performed using the Bonferroni method. Due to the design of a fixed volume injection approach, the dose of pre-bolus slightly varied with the patient body mass, whilst keeping the same length of the bolus (1 s) for all subjects. This allowed the sensitivity of different VOF/VIFs to small variations in GBCA dose to be assessed. For each vessel region (SSS, ICA, MCA) the relationship between extracted features and administered GBCA dose (mmol/kg) was assessed using scatterplots and correlation analysis. Correlation analysis was also used to evaluate the relationship of the GBCA bolus arrival time delay between each arterial (ICA, MCA) VIF and VOFSSS, and differences in bolus peak-amplitude, bolus peak width (FWHM) and the bolus PWP between the SSS and either ICA or MCA. In particular, the correlation of the BAT delay with either the absolute difference in each semiquantitative parameter (bolus peak-amplitude, bolus peak FWHM, bolus PWP) or the ratio of each parameter between the SSS and ICA/MCA (e.g., PeakSSS/PeakICA/MCA ratio) was assessed.
The intra-subject repeatability of each VOF/VIF semiquantitative parameter was assessed across the twelve patients with NF2 related VS who were imaged twice using the test–retest coefficient of variation (CoV). The CoV is the standard deviation, σ, across all measurements for each subject, divided by the mean, μ, for that subject. For a group of N subjects the global test–retest CoV is defined as \(\sqrt{\sum {(\sigma /\mu )}^{2}/N}\)46,47. As a supporting measure of repeatability the average measures intraclass correlation coefficient (ICC) of each VOF/VIF semiquantitative parameter across the two visits was also calculated using an absolute-agreement, 2-way mixed-effects model48. The inter-tumour correlation between DCE-MRI derived parameter estimates (Ktrans, vp, ve and vi) and tissue-derived metrics (H&E cell density, CD31% microvessel surface area, fibrinogen optical density) for the 15 resected sporadic VS are reported as Pearson’s product moment correlation coefficient (r).
Results
Compared to large intracranial arteries (ICA, MCA) VOFSSS provides a global surrogate input function with higher CNR, higher peak and higher peak∙FWHM product (PWP)
In Fig. 2, representative plasma GBCA concentration curves Cp(t) derived from the SSS, the ICA, and the horizontal segment of the middle cerebral artery MCA following a bolus injection of 0.016 mmol/kg of GBCA are shown. As shown in Table 1, across all patient datasets the CNR of VOFSSS (mean CNR = 197.2 ± 95.7) was significantly higher than either VIFICA (mean CNR = 51.3 ± 25.9, p < 0.001) or VIFMCA (52.9 ± 22.7, p < 0.001, repeated measures ANOVA). There was no difference in CNR between VIFICA and VIFMCA (p > 0.05, repeated measures ANOVA, Table 1). Compared to either the VIFICA or VIFMCA, VOFSSS displayed significantly longer BAT (p < 0.001), higher peak (p < 0.001), larger bolus PWP (p < 0.001) and a non-significantly narrower FWHM (p > 0.05, repeated measures ANOVA, Table 1). There was no significant difference in either BAT (p > 0.05) or FWHM (p > 0.05) between VIFICA and VIFMCA, but VIFMCA displayed a higher peak (p < 0.001) and higher PWP (p < 0.001) than VIFICA.
In Fig. 3 a comparison is shown between an individual patient derived VOFSSS and the population VIF measured in the descending aorta by Parker et al13. In keeping with the observed lack of bolus widening (bolus FWHM) between the large intracranial arteries and the SSS, after converting the patient derived VOFSSS to a full-dose input function by summing several time-shifted low-dose GBCA concentration–time curves49, there was a close resemblance between the summed SSS defined VOF and the Parker population VIF13.
VOFSSS demonstrates a greater sensitivity to interindividual changes in plasma GBCA concentration compared to arterial approaches
Across all sixty-six patients, body weight adjusted GBCA dose for the pre-bolus injection and LDHT acquisition varied from 0.0091 to 0.027 mmol/kg with a mean injected dose of 0.016 mmol/kg. As shown in Fig. 4 both bolus peak (mM) and bolus PWP (mM⋅s) correlated significantly (p ≤ 0.002) with body weight adjusted GBCA dose with the strongest correlation observed for VOFSSS peak (r = 0.70, p < 0.001) and VOFSSS bolus PWP (r = 0.65, p < 0.001). No significant correlation was observed between body weight adjusted GBCA dose and bolus peak width (FWHM) for any derived VOF/VIF (p > 0.05, Pearson’s correlation co-efficient). As shown in Table 2, no significant correlation (p > 0.05) was observed between GBCA bolus PWP ratio (PWPSSS/PWPICA or PWPSSS/PWPMCA) and the BAT difference between SSS and ICA or between SSS and MCA respectively. A longer time interval between BATICA and BATSSS correlated with an increase in the FWHMSSS/FWHMICA ratio (r = 0.35, p = 0.004) and was non-significantly associated with a lower PeakSSS/PeakICA ratio (r = − 0.18, p = 0.15).
VOFSSS demonstrates greater repeatability compared to arterial VIF approaches
Table 3 shows the intra-subject variability and repeatability of each semiquantitative VOF/VIF parameter across the twelve patients with NF2-related VS. In the case of VOFSSS, global CoV values for BAT, bolus peak, bolus FWHM and bolus PWP were 3.98%, 17.0%, 16.8% and 12.4% respectively. Except for BAT, global CoV values for VIFSSS were lower than the corresponding global CoV values for VIFICA or VIFMCA. Across all semi-quantitative parameters extracted from the VOFSSS repeatability was good to excellent (ICC = 0.717–0.888)50,51.
Kinetic parameters obtained using a SSS derived VOF permitted detection of intertumoural differences in histopathological data
In Fig. 5 the inter-tumour correlation between LEGATOS derived kinetic parameter estimates and tissue metrics for each VOF/VIF are shown. There was a significant correlation between cell density and mean tumour vi when using VOFSSS (r = 0.54, p = 0.04, Fig. 5A). No such correlation was seen, however, when using either VIFICA or VIFMCA. In many tumours there was overestimation of ve when using arterial VIF (6/15 and 4/15 VS had ve > 0.7 when using VIFICA and VIFMCA respectively), and such high EES fractions were not evident on collected tissue from these tumours (Fig. 6). For VOFSSS and VIFICA a significant positive correlation was seen between CD31% microvessel surface area and mean tumour vp (p < 0.05), but this correlation was strongest for vp maps derived using VOFSSS (r = 0.85, p < 0.001, Fig. 5B). Estimates of both vp and Ktrans when using an arterial VIF were higher than VOFSSS derived estimates. Across all VOF/VIF approaches, there was a correlation of Ktrans with both CD31% microvessel surface area (p < 0.05) and perivascular leak as measured through fibrinogen optical density (p < 0.05).
Discussion
Compared to large intracranial arteries (ICA, MCA), an arrival-corrected VOF derived from the SSS (VOFSSS) can provide a superior surrogate global input function for kinetic parameter analysis in brain DCE-MRI, demonstrating higher CNR, higher peak and a greater sensitivity to interindividual changes in plasma GBCA concentration. Through an included test–retest study we demonstrate that semi-quantitative parameters derived from VOFSSS display greater repeatability than arterial based approaches. Furthermore, through comparison with matched tissue datasets in patients with resected sporadic VS we demonstrate that microvascular kinetic parameters obtained using an arrival-corrected VOFSSS permit evaluation of intertumoural differences in microvessel surface area and cell density, a feature not seen with large artery based VIFs.
Previous studies comparing arterial and venous based VOF/VIFs have reported similar results to our findings, and several factors can be hypothesized to contribute to the observed higher CNR and higher peak measured within VOFSSS compared to large arteries19,32. Keil et al. demonstrated that arterial VIF showed lower peak Cp values and more often led to implausibly high ve values in subsequent pharmacokinetic parameter fitting compared to a SSS derived VOF19. The larger cross-sectional area of the SSS relative to both the ICA and MCA may serve to reduce PVE and contribute to the higher CNR, higher peak GBCA bolus measured32,39. Reduced ‘in-flow’ effects secondary to lower flow velocities within the SSS may also contribute to the higher CNR observed25,52,53,54. To maximise T1 weighting, a 3D T1W GRE sequence with the shortest TE was used for the DCE-MRI acquisition. Such GRE sequences are, however, prone to ‘in-flow’ related enhancement, which can lead to a significant attenuation of contrast enhancement by GBCA and subsequent reduction of measured VIF CNR27,30. Within our study, measures were taken to prevent these ‘in-flow’ artefacts such as RF phase cycling, use of a large acquisition volume covering the whole brain, and orientation of the read gradient of the 3D slab parallel to the vessels. It is, however, possible that such ‘in-flow’ related artefacts still occurred39 and because the contrast between flowing and static tissue depends directly upon the flow velocity53,54, the higher velocity arteries are much more prone to increased ‘in-flow’ effects and resulting lower CNR than lower velocity venous structures such as the SSS25,39,52.
An additional factor that may have contributed to the observed decreased CNR and peak within arterial VIFs relative to the SSS is differences in inter- and intravoxel velocity dispersion within these vessels39. Compared to large intracranial arteries, blood flow within the SSS is thought to be less turbulent and more laminar, with a parabolic velocity profile55,56. Assuming laminar flow and a parabolic velocity profile, then the velocity v of a given proton spin across a circular vessel, is a function of distance to the isocentre of the flow, i.e., v(r) = vmax (1−r2/a2), where a is the radius of the vessel, r is the distance between the position of the flowing spin and the centre of the vessel, v(r) is the velocity at r, and vmax is the maximum velocity at the centre of the vessel. Whereas stationary spins are completely rephased during data acquisition, moving spins retain a degree of phase shift that is a function of both the first moment of the gradient M1 and the flow velocity v. The phase increases linearly with the flow velocity but there is no effect on the spatial location of the signal by 2D Fourier transform reconstruction because the velocity is set constant from view to view. The location would, however, be affected if acceleration in the vessel is heterogenous and in vessels where there is significant intra- and intervoxel velocity heterogeneity, phase dispersion and loss of CNR can occur. Moreover, changes of the parabolic velocity profile can influence the width of first-pass of the contrast bolus. In vivo blood flow measurement is considerably more complex, blood is non-Newtonian for example and the flow profile depends on the geometry of the vessels and the sampling point of the pulsed flow, i.e., the time points in the cardiac cycle. The field gradients for spatial encoding and the spoiler gradients used in our acquisition were not specifically designed for ‘flow compensation’ and extra spin dephasing caused by spins moving during the actual signal encoding may have occurred. Such flow-induced loss of phase coherence increases with higher velocity blood flow and may in addition to PVE and residual ‘in-flow’ effects explain the lower CNR and lower peak measured within arterial VIFs compared to VOFSSS.
Within this study a rapid bolus injection protocol was used, with a compact low-dose bolus of GBCA injected over 1 s. Following injection into an upper limb vein, the GBCA bolus is expected to undergo a degree of dispersion and widening as it passes through first the pulmonary and then cerebral circulation. This dispersion process is commonly modelled as a mathematical convolution with the vascular system of the brain and a simple yet useful model for the dispersion process might contain a single artery that feeds a set of parallel pathways (an arteriole, a blood-tissue exchange unit, and a venule) that drain into a single vein57,58,59,60,61. Data from our study, however, suggests that when first-pass low-dose GBCA bolus uptake curves are measured within the distal SSS such dispersion is minimal, with no significant difference in peak width (FWHM) between VOFSSS and the VIF measured from feeding larger arteries such as the ICA and MCA. This absence of significance bolus dispersion between the ICA and SSS has also been reported in other studies and several factors may underlie this unexpected observation8,39. A large contributing factor is likely inaccuracies in the measurement of the VIF within the smaller arterial vessels due to increased PVE, increased vessel inflow effects and increased velocity dispersion of flowing blood relative to the SSS. Such effects likely contribute to the increased CNR seen within the VOFSSS but also counterbalance the broadening effect by dispersion8.
Previous works have demonstrated that incorrect estimation of VIF peak and the bolus arrival time can have significant impact on pharmacokinetic parameter fitting and accuracy62,63,64,65. In particular the bolus arrival time in the SSS is 4–5 s later than the true input arrival time and a correction must be made for the delayed arrival of the contrast agent within the venous system by incorporating voxel-wise BAT estimation in DCE-MRI analysis66,67. Within this study use of a DTR acquisition incorporating high temporal sampling during pre-bolus injection was adopted, improving voxel-by-voxel estimation of BAT delay, and reducing fitting errors induced by uncertainty in time alignment of the VOF/VIF and tissue uptake curves. Alongside measurement effects, physiological factors may also play a role in increasing CNR and reducing the observed dispersion seen within the SSS relative to arterial VIF sources. The SSS serves as the principal draining vein for the cerebral cortex and it’s comparatively straight course, absence of valves, and wide diameter relative to large intracranial arteries means that laminar blood flow characteristics within it resemble that of similar sized arteries22,24,25,68. Indeed a similarity between the measured VOF within the SSS and Parker’s population averaged VIF has been reported in both this and other studies13,19. Human histological studies and microsphere studies in rhesus monkeys have demonstrated that within normal brain under physiological conditions, a small percentage of blood from the arteries passes through small pial arteriovenous shunts (12.5 µm or more in diameter) directly into cerebral veins and the SSS69,70,71 Direct shunting of blood into cerebral veins has also been demonstrated angiographically within supratentorial GBM and VS, although in the latter case the source of principal venous drainage is the transverse and sigmoid sinuses rather than the SSS itself72,73,74,75,76,77,78. Although in theory the presence of such shunts may serve to reduce dispersion of the GBCA bolus within the venous system, the predicted low volume of blood passing into the SSS via these direct channels makes them unlikely to be a significant factor in driving the observed VOFSSS characteristics.
The presented study is one of the largest comparisons of different arterial and venous based input functions for brain DCE-MRI and to our knowledge one of the first studies to evaluate the sensitivity of different VOF/VIF approaches to interindividual differences in plasma GBCA concentration and histopathological data. A limitation of the present study though is that characterization of plasma GBCA concentration curves, Cp(t), was limited to semi-quantitative parameters such as the bolus peak-amplitude, bolus width and GBCA bolus PWP. Future studies should seek to undertake more detailed and sophisticated shape analysis of the whole length of the VOF/VIF. In this study, a semi-automatic extraction method was employed to measure each VOF/VIF, in which the twenty vessel voxels with the highest AUC30 are selected and averaged together to create the final input curve 18,37,79. This is under the assumption that these voxels demonstrate less inflow effects and PVE. Although it is possible that such an approach may overestimate the true value of VOFSSS, our demonstration that the VOFSSS showed an initial peak height and first-pass bolus shape very close to the ‘gold standard’ Parker VIF model13,80,81, and that the variation in the VOFSSS bolus peak was strongly correlated with dose variation, suggests that such overestimation was minimal and that peak estimates are not dominated by bolus shape distortion. Our demonstration that kinetic parameter estimates, obtained using VOFSSS correlated well with tissue derived measures of microvessel surface area and cell density, in contrast to arterial VIFs that showed lower peak values and overestimation of ve, further supports the robustness of input function surrogate measurements from the SSS and the absence of significant peak overestimation through this semi-automatic voxel selection method. Larger studies incorporating matched imaging-tissue cohorts in a range of different tumours should, however, be undertaken to better evaluate the effect of VOF/VIF location and voxel extraction method on kinetic parameter accuracy.
Conclusion
Accurate derivation of a vascular input function (VIF) is essential for quantitative kinetic analysis of brain DCE-MRI data. In this in vivo patient study, we compared VIFs extracted from either the internal carotid artery and its branches with an arrival-corrected vascular output function derived from the superior sagittal sinus (VOFSSS). We demonstrated that compared to large intracranial arteries VOFSSS can provide a superior surrogate global VIF, with lower noise, higher repeatability, and greater sensitivity to interindividual changes in plasma GBCA concentration. Through comparison with matched histopathological data, we furthermore demonstrate that microvascular parameters obtained using a SSS derived VOF permitted evaluation of intertumoural differences in microvessel surface area and cell density. These results support the use of venous sinus-based approaches for input function extraction and pharmacokinetic parameter mapping in brain DCE-MRI.
Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- ANOVA:
-
Analysis of variance
- AUC30:
-
Area under the enhancing curve within 30 s of the bolus arrival time, mM·s
- BAT:
-
Bolus arrival time, seconds (s)
- Cp(t):
-
Plasma GBCA concentration time course, mM
- CD31:
-
Cluster of differentiation 31 synonym platelet endothelial cell adhesion molecule (PECAM-1), an endothelial marker
- CNR:
-
Contrast to noise ratio
- COV:
-
Test–retest coefficient of variation, %
- Δt:
-
Update time/frame rate of DCE-MRI acquisition, seconds (s)
- DCE-MRI:
-
Dynamic contrast enhanced MRI
- DTR:
-
Dual temporal resolution
- EES:
-
Extravascular-extracellular space
- ETM:
-
Extended Tofts model
- FDHS DCE:
-
Full dose high spatial DCE
- FOV:
-
Field of view
- FWHM:
-
Full-width at half maximum of bolus curve or bolus width, seconds (s)
- GBCA:
-
Gadolinium-based contrast agent/s
- GBM:
-
Glioblastoma multiforme synonym WHO grade IV glioma
- GRE:
-
Gradient recalled echo
- ICA:
-
Internal carotid artery
- ICC:
-
Average measures intraclass correlation coefficient
- Ktrans :
-
Volume transfer constant, min-1
- LDHT DCE:
-
Low-dose high temporal DCE
- LEGATOS:
-
LEvel and rescale the GAdolinium contrast concentrations curves of high temporal TO high spatial
- MCA:
-
Middle cerebral artery
- NF2:
-
Neurofibromatosis type II
- PWP:
-
Bolus peak∙FWHM product, a parameter associated with the area under the bolus curve, mM·s
- PVE:
-
Partial volume errors
- R10:
-
Native longitudinal relaxation rate (1/T10), s-1
- ROI:
-
Region of interest
- SI:
-
Signal intensity
- SFE:
-
Scaled fitting error, %
- SSS:
-
Superior sagittal sinus
- T1W:
-
Longitudinal relaxation (T1) weighted (MR imaging)
- TE/TR:
-
Echo time/ Repetition time, ms
- ve :
-
Volume of extravascular-extracellular space per unit volume of tissue, %
- vi :
-
Volume of intracellular space per unit volume of tissue (estimated), %
- vp :
-
Blood plasma volume per unit volume of tissue, vp = CBV(1—Haematocrit), where CBV is cerebral blood volume, %
- VIFICA :
-
Vascular input function (plasma concentration time course) derived from internal carotid artery just distal to the carotid syphon, mM
- VIFMCA :
-
Vascular input function (plasma concentration time course) derived from horizontal segment of the middle cerebral artery, mM
- VOFSSS :
-
Vascular output function (plasma concentration time course) derived from vertical/posterior third of the superior sagittal sinus, mM
- VFA:
-
Variable flip angle
- VS:
-
Vestibular schwannoma
References
Lewis, D. et al. The microenvironment in sporadic and neurofibromatosis type II–related vestibular schwannoma: the same tumor or different? A comparative imaging and neuropathology study. J. Neurosurg. https://doi.org/10.3171/2020.3.jns193230 (2020).
Lewis, D. et al. Inflammation and vascular permeability correlate with growth in sporadic vestibular schwannoma. Neuro. Oncol. 21, 314–325 (2019).
O’Connor, J. P. B., Jackson, A., Parker, G. J. M., Roberts, C. & Jayson, G. C. Dynamic contrast-enhanced MRI in clinical trials of antivascular therapies. Nat. Rev. Clin. Oncol. 9, 167–177 (2012).
Jain, R. Measurements of tumor vascular leakiness using DCE in brain tumors: clinical applications. NMR Biomed. 26, 1042–1049 (2013).
Li, K. L. et al. Vascular biomarkers derived from dynamic contrast-enhanced MRI predict response of vestibular schwannoma to antiangiogenic therapy in type 2 neurofibromatosis. Neuro. Oncol. 18, 275–282 (2016).
Li, K.-L. et al. The LEGATOS technique: a new tissue-validated dynamic contrast-enhanced MRI method for whole-brain, high-spatial resolution parametric mapping. Magn. Reson. Med. https://doi.org/10.1002/mrm.28842 (2021).
Larsson, H. B. W., Courivaud, F., Rostrup, E. & Hansen, A. E. Measurement of brain perfusion, blood volume, and blood-brain barrier permeability, using dynamic contrast-enhanced T1-weighted MRI at 3 tesla. Magn. Reson. Med. 62, 1270–1281 (2009).
Jelescu, I. O. et al. Dual-temporal resolution dynamic contrast-enhanced MRI protocol for blood-brain barrier permeability measurement in enhancing multiple sclerosis lesions. J. Magn. Reson. Imaging https://doi.org/10.1002/jmri.22565 (2011).
Van De Haar, H. J. et al. Blood-brain barrier leakage in patients with early Alzheimer disease. Radiology https://doi.org/10.1148/radiol.2016152244 (2016).
Tofts, P. S. Modeling tracer kinetics in dynamic Gd-DTPA MR imaging. J. Magn. Reson. Imaging 7, 91–101 (1997).
Yankeelov, T. & Gore, J. Dynamic contrast enhanced magnetic resonance imaging in oncology:theory, data acquisition, analysis, and examples. Curr. Med. Imaging Rev. 3, 91–107 (2007).
Weinmann, H. J., Laniado, M. & Mützel, W. Pharmacokinetics of GdDTPA/dimeglumine after intravenous injection into healthy volunteers. Physiol. Chem. Phys. Med. NMR 16, 167–172 (1984).
Parker, G. J. M. et al. Experimentally-derived functional form for a population-averaged high-temporal-resolution arterial input function for dynamic contrast-enhanced MRI. Magn. Reson. Med. 56, 993–1000 (2006).
Parker, G. J. M. & Buckley, D. L. Tracer kinetic modelling for T1-weighted DCE-MRI. In Dynamic Contrast-Enhanced Magnetic Resonance Imaging in Oncology (eds Jackson, A. et al.) 81–92 (Springer, Berlin, 2005). https://doi.org/10.1007/3-540-26420-5_6.
Buckley, D. L. Uncertainty in the analysis of tracer kinetics using dynamic contrast-enhanced T1-weighted MRI. Magn. Reson. Med. 47, 601–606 (2002).
Li, K. L., Lewis, D., Jackson, A., Zhao, S. & Zhu, X. Low-dose T1W DCE-MRI for early time points perfusion measurement in patients with intracranial tumors: a pilot study applying the microsphere model to measure absolute cerebral blood flow. J. Magn. Reson. Imaging 48, 543–557 (2018).
Yang, C. et al. Comparison of quantitative parameters in cervix cancer measured by dynamic contrast-enhanced MRI and CT. Magn. Reson. Med. 63, 1601–1609 (2010).
Sourbron, S., Ingrisch, M., Siefert, A., Reiser, M. & Herrmann, K. Quantification of cerebral blood flow, cerebral blood volume, and blood-brain-barrier leakage with DCE-MRI. Magn. Reson. Med. 62, 205–217 (2009).
Keil, V. C. et al. Effects of arterial input function selection on kinetic parameters in brain dynamic contrast-enhanced MRI. Magn. Reson. Imaging 40, 83–90 (2017).
Li, K. L., Zhu, X. P., Waterton, J. & Jackson, A. Improved 3D quantitative mapping of blood volume and endothelial permeability in brain tumors. J. Magn. Reson. Imaging 12, 347–357 (2000).
Lavini, C. & Verhoeff, J. J. C. Reproducibility of the gadolinium concentration measurements and of the fitting parameters of the vascular input function in the superior sagittal sinus in a patient population. Magn. Reson. Imaging 28, 1420–1430 (2010).
Hirabuki, N. et al. Quantitation of flow in the superior sagittal sinus performed with cine phase-contrast MR imaging of healthy and achondroplastic children. AJNR. Am. J. Neuroradiol. 21, 1497–1501 (2000).
Haroon, H. A. et al. A comparison of Ktrans measurements obtained with conventional and first pass pharmacokinetic models in human gliomas. J. Magn. Reson. Imaging 19, 527–536 (2004).
Mattle, H., Edelman, R. R., Reis, M. A. & Atkinson, D. J. Flow quantification in the superior sagittal sinus using magnetic resonance. Neurology 40, 813–815 (1990).
Inao, S., Kuchiwaki, H., Yoshida, J. & Furuse, M. Magnetic resonance imaging quantitation of superior sagittal sinus flow: correlation to cerebral blood flow measured by xenon-enhanced computed tomography. Neurol. Res. 19, 35–40 (1997).
Gowland, P. et al. Dynamic studies of gadolinium uptake in brain tumors using inversion-recovery echo-planar imaging. Magn. Reson. Med. 26, 241–258 (1992).
Roberts, C. et al. The effect of blood inflow and B(1)-field inhomogeneity on measurement of the arterial input function in axial 3D spoiled gradient echo dynamic contrast-enhanced MRI. Magn. Reson. Med. 65, 108–119 (2011).
Li, K.-L. et al. An improved coverage and spatial resolution-using dual injection dynamic contrast-enhanced (ICE-DICE) MRI: a novel dynamic contrast-enhanced technique for cerebral tumors. Magn. Reson. Med. 68, 452–462 (2012).
Wang, S. et al. Arterial input functions (AIFs) measured directly from arteries with low and standard doses of contrast agent, and AIFs derived from reference tissues. Magn. Reson. Imaging 34, 197–203 (2016).
Garpebring, A., Wirestam, R., Ostlund, N. & Karlsson, M. Effects of inflow and radiofrequency spoiling on the arterial input function in dynamic contrast-enhanced MRI: a combined phantom and simulation study. Magn. Reson. Med. 65, 1670–1679 (2011).
Barth, M. & Moser, E. Proton NMR relaxation times of human blood samples at 1.5 T and implications for functional MRI. Cell. Mol. Biol. (Noisy-le-grand) 43, 783–791 (1997).
Foottit, C., Cron, G. O., Hogan, M. J., Nguyen, T. B. & Cameron, I. Determination of the venous output function from MR signal phase: feasibility for quantitative DCE-MRI in human brain. Magn. Reson. Med. 63, 772–781 (2010).
Zhu, X. P. et al. Quantification of endothelial permeability, leakage space, and blood volume in brain tumors using combined T1 and T2* contrast-enhanced dynamic MR imaging. J. Magn. Reson. Imaging 11, 575–585 (2000).
van Schie, J. J. N., Lavini, C., van Vliet, L. J. & Vos, F. M. Estimating the arterial input function from dynamic contrast-enhanced MRI data with compensation for flow enhancement (I): theory, method, and phantom experiments. J. Magn. Reson. Imaging 47, 1190–1196 (2018).
Ashburner, J. & Friston, K. Multimodal image coregistration and partitioning - a unified framework. Neuroimage 6, 209–217 (1997).
Rosen, B. R., Belliveau, J. W., Vevea, J. M. & Brady, T. J. Perfusion imaging with NMR contrast agents. Magn. Reson. Med. 14, 249–265 (1990).
Li, K.-L., Zhu, X., Zhao, S. & Jackson, A. Blood–brain barrier permeability of normal-appearing white matter in patients with vestibular schwannoma: a new hybrid approach for analysis of T1 W DCE-MRI. J. Magn. Reson. Imaging https://doi.org/10.1002/jmri.25573 (2017).
Gwilliam, M. N., Collins, D. J., Leach, M. O. & Orton, M. R. Quantifying MRI T(1) relaxation in flowing blood: implications for arterial input function measurement in DCE-MRI. Br. J. Radiol. 94, 20191004 (2021).
Bourassa-Moreau, B., Lebel, R., Gilbert, G., Mathieu, D. & Lepage, M. Robust arterial input function surrogate measurement from the superior sagittal sinus complex signal for fast dynamic contrast-enhanced MRI in the brain. Magn. Reson. Med. https://doi.org/10.1002/mrm.28922 (2021).
Duan, C. et al. Are complex DCE-MRI models supported by clinical data?. Magn. Reson. Med. 77, 1329–1339 (2017).
Stanisz, G. J. et al. T1, T2 relaxation and magnetization transfer in tissue at 3T. Magn. Reson. Med. 54, 507–512 (2005).
Shen, Y. et al. T1 relaxivities of gadolinium-based magnetic resonance contrast agents in human whole blood at 1.5, 3, and 7 T. Invest. Radiol. 50, 330–338 (2015).
Tofts, P. S. et al. Estimating kinetic parameters from dynamic contrast-enhanced T1- weighted MRI of a diffusable tracer: standardized quantities and symbols. J. Magn. Reson. Imaging 10, 223–232 (1999).
Henderson, E., Rutt, B. K. & Lee, T. Y. Temporal sampling requirements for the tracer kinetics modeling of breast disease. Magn. Reson. Imaging 16, 1057–1073 (1998).
Li, K.-L. et al. Heterogeneity in the angiogenic response of a BT474 human breast cancer to a novel vascular endothelial growth factor-receptor tyrosine kinase inhibitor: assessment by voxel analysis of dynamic contrast-enhanced MRI. J. Magn. Reson. Imaging 22, 511–519 (2005).
Jahng, G. H. et al. Human brain: reliability and reproducibility of pulsed arterial spin-labeling perfusion MR imaging. Radiology https://doi.org/10.1148/radiol.2343031499 (2005).
Waterton, J. C. et al. Diurnal variation in the femoral articular cartilage of the knee in young adult humans. Magn. Reson. Med. https://doi.org/10.1002/(sici)1522-2594(200001)43:1%3c126::aid-mrm15%3e3.0.co;2-%23 (2000).
Koo, T. K. & Li, M. Y. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med. 15, 155–163 (2016).
Köstler, H. et al. Prebolus quantitative MR heart perfusion imaging. Magn. Reson. Med. 52, 296–299 (2004).
Erickson, N. J. et al. Koos classification of vestibular schwannomas: a reliability study. Neurosurgery 85, 409–414 (2019).
Lenth, R. V. & Fleiss, J. L. Design and analysis of clinical experiments (Wiley Classics Library). J. Am. Stat. Assoc. 94, 1384 (1999).
Blackshear, W. M. et al. Carotid artery velocity patterns in normal and stenotic vessels. Stroke 11, 67–71 (1980).
Singer, J. R. & Crooks, L. E. Nuclear magnetic resonance blood flow measurements in the human brain. Science 221, 654–656 (1983).
Feinberg, D. A., Crooks, L., Hoenninger, J. 3rd., Arakawa, M. & Watts, J. Pulsatile blood velocity in human arteries displayed by magnetic resonance imaging. Radiology 153, 177–180 (1984).
Schuchardt, F. et al. Acute cerebral venous thrombosis: three-dimensional visualization and quantification of hemodynamic alterations using 4-dimensional flow magnetic resonance imaging. Stroke 48, 671–677 (2017).
Schuchardt, F. et al. In vivo analysis of physiological 3D blood flow of cerebral veins. Eur. Radiol. 25, 2371–2380 (2015).
King, R. B., Deussen, A., Raymond, G. M. & Bassingthwaighte, J. B. A vascular transport operator. Am. J. Physiol. 265, H2196–H2208 (1993).
Calamante, F. Bolus dispersion issues related to the quantification of perfusion MRI data. J. Magn. Reson. Imaging 22, 718–722 (2005).
Calamante, F., Willats, L., Gadian, D. G. & Connelly, A. Bolus delay and dispersion in perfusion MRI: implications for tissue predictor models in stroke. Magn. Reson. Med. 55, 1180–1185 (2006).
Ewing, J. R. et al. Patlak plots of Gd-DTPA MRI data yield blood-brain transfer constants concordant with those of 14C-sucrose in areas of blood-brain opening. Magn. Reson. Med. 50, 283–292 (2003).
Zierler, K. L. Theory of the use of arteriovenous concentration differences for measuring metabolism in steady and non-steady states. J. Clin. Invest. 40, 2111–2125 (1961).
Port, R. E., Knopp, M. V. & Brix, G. Dynamic contrast-enhanced MRI using Gd-DTPA: interindividual variability of the arterial input function and consequences for the assessment of kinetics in tumors. Magn. Reson. Med. 45, 1030–1038 (2001).
O’Connor, J. P. B. et al. Dynamic contrast-enhanced imaging techniques: CT and MRI. Br. J. Radiol. 84, S112-20 (2011).
Lavini, C. Simulating the effect of input errors on the accuracy of Tofts’ pharmacokinetic model parameters. Magn. Reson. Imaging 33, 222–235 (2015).
Li, X. et al. Relative sensitivities of DCE-MRI pharmacokinetic parameters to arterial input function (AIF) scaling. J. Magn. Reson. 269, 104–112 (2016).
Singh, A. et al. Improved bolus arrival time and arterial input function estimation for tracer kinetic analysis in DCE-MRI. J. Magn. Reson. Imaging 29, 166–176 (2009).
Nadav, G., Liberman, G., Artzi, M., Kiryati, N. & Bashat, D. B. Optimization of two-compartment-exchange-model analysis for dynamic contrast-enhanced mri incorporating bolus arrival time. J. Magn. Reson. Imaging 45, 237–249 (2017).
Beards, S. C., Yule, S., Kassner, A. & Jackson, A. Anatomical variation of cerebral venous drainage: the theoretical effect on jugular bulb blood samples. Anaesthesia 53, 627–633 (1998).
Rowbotham, G. F. & Little, E. A new concept of the circulation and the circulations of the brain. The discovery of surface arteriovenous shunts. Br. J. Surg. 52, 539–542 (1965).
Rowed, D. W., Stark, V. J., Hoffer, P. B. & Mullan, S. Cerebral arteriovenous shunts re-examined. Stroke 3, 592–600 (1972).
Elands, S. et al. Early venous filling following thrombectomy: association with hemorrhagic transformation and functional outcome. Front. Neurol. 12, 649079 (2021).
Teranishi, Y., Kohno, M., Sora, S., Sato, H. & Nagata, O. Hypervascular vestibular schwannomas: clinical characteristics, angiographical classification, and surgical considerations. Oper. Neurosurg. (Hagerstown, Md) 15, 251–261 (2018).
Perneczky, A. Blood supply of acoustic neurinomas. Acta Neurochir. (Wien) 52, 209–218 (1980).
Maruki, C. et al. Acoustic neurinoma with early venous drainage and caput medusae-like vasculature on angiogram. Case report. Neurol. Med. Chir. (Tokyo) 26, 989–992 (1986).
Yamada, S., Aiba, T. & Hara, M. Early venous filling in acoustic schwannoma. Radiat. Med. 5, 10–13 (1987).
Takahashi, M., Okudera, T., Tomanaga, M. & Kitamura, K. Angiographic diagnosis of acoustic neurinomas: analysis of 30 lesions. Neuroradiology 2, 191–200 (1971).
Mariani, L., Schroth, G., Wielepp, J. P., Haldemann, A. & Seiler, R. W. Intratumoral arteriovenous shunting in malignant gliomas. Neurosurgery 48, 353–358 (2001).
Yoshikawa, A. et al. Visualization of angiographical arteriovenous shunting in perisylvian glioblastomas. Acta Neurochir. (Wien) 155, 715–719 (2013).
Larsson, C. et al. Sampling requirements in DCE-MRI based analysis of high grade gliomas: simulations and clinical results. J. Magn. Reson. Imaging 37, 818–829 (2013).
He, D., Xu, L., Qian, W., Clarke, J. & Fan, X. A simulation study comparing nine mathematical models of arterial input function for dynamic contrast enhanced MRI to the Parker model. Australas. Phys. Eng. Sci. Med. 41, 507–518 (2018).
Aerts, H. J. W. L., Jaspers, K. & Backes, W. H. The precision of pharmacokinetic parameters in dynamic contrast-enhanced magnetic resonance imaging: the effect of sampling frequency and duration. Phys. Med. Biol. 56, 5665–5678 (2011).
Acknowledgements
This study was funded by the CRUK and EPSRC Cancer Imaging Centre in Cambridge and Manchester (Award Number: C8742/A18097). Thanks go to Dr Erjon Agushi, Mr Omar Pathmanaban and all the staff at the adult neuro-oncology service and Manchester Skull Base unit at Salford Royal NHS Hospital for their help with patient recruitment and tissue collection. Our thanks also go to Dr Josephine Naish and Dr David Higgins for the many useful discussions. A portion of this work has been presented as an electronic poster abstract at the 2021 ISMRM & SMRT Annual Meeting & Exhibition.
Author information
Authors and Affiliations
Contributions
D.L., X.P., D.J.C., T.C., A.J. and K.L. contributed to conception and design of this study. D.L., D.J.C. and A.T.K. contributed to patient recruitment. D.L., X.P., D.J.C., S.Z., T.C., A.J., and K.L. contributed to data analysis. All authors contributed to the writing and review of this manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lewis, D., Zhu, X., Coope, D.J. et al. Surrogate vascular input function measurements from the superior sagittal sinus are repeatable and provide tissue-validated kinetic parameters in brain DCE-MRI. Sci Rep 12, 8737 (2022). https://doi.org/10.1038/s41598-022-12582-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-12582-x
This article is cited by
-
Blood–brain barrier integrity is linked to cognitive function, but not to cerebral arterial pulsatility, among elderly
Scientific Reports (2024)
-
Low-dose GBCA administration for brain tumour dynamic contrast enhanced MRI: a feasibility study
Scientific Reports (2024)
-
Model-free dynamic contrast-enhanced MRI analysis: differentiation between active tumor and necrotic tissue in patients with glioblastoma
Magnetic Resonance Materials in Physics, Biology and Medicine (2022)