- Open Access
Population pharmacokinetic modeling of CSF to blood clearance: prospective tracer study of 161 patients under work-up for CSF disorders
Fluids and Barriers of the CNS volume 19, Article number: 55 (2022)
Quantitative measurements of cerebrospinal fluid to blood clearance has previously not been established for neurological diseases. Possibly, variability in cerebrospinal fluid clearance may affect the underlying disease process and may possibly be a source of under- or over-dosage of intrathecally administered drugs. The aim of this study was to characterize the cerebrospinal fluid to blood clearance of the intrathecally administered magnetic resonance imaging contrast agent gadobutrol (Gadovist, Bayer Pharma AG, GE). For this, we established a population pharmacokinetic model, hypothesizing that cerebrospinal fluid to blood clearance differs between cerebrospinal fluid diseases.
Gadobutrol served as a surrogate tracer for extra-vascular pathways taken by several brain metabolites and drugs in cerebrospinal fluid. We estimated cerebrospinal fluid to blood clearance in patients with different cerebrospinal fluid disorders, i.e. symptomatic pineal and arachnoid cysts, as well as tentative spontaneous intracranial hypotension due to cerebrospinal fluid leakage, idiopathic intracranial hypertension, or different types of hydrocephalus (idiopathic normal pressure hydrocephalus, communicating- and non-communicating hydrocephalus). Individuals with no verified cerebrospinal fluid disturbance at clinical work-up were denoted references.
Population pharmacokinetic modelling based on 1,140 blood samples from 161 individuals revealed marked inter-individual variability in pharmacokinetic profiles, including differences in absorption half-life (time to 50% of tracer absorbed from cerebrospinal fluid to blood), time to maximum concentration in blood and the maximum concentration in blood as well as the area under the plasma concentration time curve from zero to infinity. In addition, the different disease categories of cerebrospinal fluid diseases demonstrated different profiles.
The present observations of considerable variation in cerebrospinal fluid to blood clearance between individuals in general and across neurological diseases, may suggest that defining cerebrospinal fluid to blood clearance can become a useful diagnostic adjunct for work-up of cerebrospinal fluid disorders. We also suggest that it may become useful for assessing clearance capacity of endogenous brain metabolites from cerebrospinal fluid, as well as measuring individual cerebrospinal fluid to blood clearance of intrathecal drugs.
While the renal glomerular filtration rate (GFR) is clinically used as marker of clearance of drugs and solutes from blood , the cerebrospinal fluid (CSF) to blood clearance has not previously been defined in either healthy individuals nor in individuals with neurological diseases. Possibly, direct measurement of CSF to blood clearance might be useful for understanding diseases of the brain, and consequently lay ground for personalized intrathecal drug administration to the central nervous system (CNS).
Since the dual discoveries of the glymphatic system in 2012  and the meningeal lymphatic system in 2015 , there have been renewed interest in how various waste substances are cleared from the brain , and not at least the role of meningeal lymphatic vessels . Impaired glymphatic clearance of toxic by-products from brain metabolism to CSF causing deposition of toxic substances in the brain, e.g. deposition of amyloid-β and tau in Alzheimer’s disease and α-synuclein in Parkinson’s disease, has been proposed as a common pathogenic pathway behind several neurodegenerative disorders . Meningeal lymphatic function seems to be affected in a wide range of diseases, as indicated in animal models of traumatic brain injury , malignant brain disease [7,8,9], stroke [10, 11] and Alzheimer’s disease , and in patients with Parkinson’s disease. Given that impaired molecular clearance from CSF to blood may have a pivotal role in the development of neurological disease; it might be desirable to obtain quantitative data about CSF to blood clearance on an individual basis. For years, levels of brain metabolites from single time points have been measured in CSF, as well as in blood, aiming at identifying the pre-symptomatic phase of dementia disease [14, 15]. On the other hand, direct assessment of clearance dynamics from CSF to blood has not been possible.
Assessment of CSF to blood clearance might as well be useful to tailor dosage of intrathecal drugs. Today, intrathecal drug administration seems promising in order to treat a wide range of diseases within the CNS, such as neuro-inflammatory, neuro-degenerative, neuro-oncologic, and neuro-vascular diseases [16,17,18,19,20]. Many systemically administered drugs, which are supposed to function in CNS, remain to a considerable degree within the systemic circulation due to their inability to cross the blood–brain-barrier (BBB) . Given previous observations of brain wide distribution of CSF tracer administered to the lumbar subarachnoid space in humans , intrathecally administered drugs have potential to better target brain disease directly by their by-passing of the BBB, and assumedly in much lower doses than applied systemically, thereby reducing side effects.
Our group has used intrathecal administration of the magnetic resonance imaging (MRI) contrast agent gadobutrol (serving as a CSF tracer) to explore molecular passage from CSF to the brain [22, 23], meninges , calvarial bone , extra-cranial lymph nodes , as well as to the blood . From this, we suggest that measurements of CSF to blood clearance of gadobutrol may provide an overall estimation of the ability of CSF to remove macromolecules. Since tracer levels in blood are highly correlated with levels of tracer in CSF at MRI , resource-demanding imaging may be omitted as part of CSF clearance assessment. Gadobutrol is a hydrophilic substance unable to cross the BBB, which after administration to CSF is excreted along the same pathways as other endogenous substances within CSF, such as the paravascular  and meningeal lymphatic pathways  suggested from animal studies. In the present work, we investigated the CSF to blood clearance of gadobutrol in patients under clinical work-up of various neurological diseases and CSF disturbances, employing a population pharmacokinetic model based on a large patient material spanning multiple disease categories. The hypothesis was that different CSF diseases present a characteristic profile of CSF to blood clearance.
A prospective and observational study design was utilized; randomization or a priori sample size calculation was not relevant.
The study included patients referred to the Department of neurosurgery, Oslo University Hospital—Rikshospitalet, Oslo, Norway, who were examined for tentative CSF disorders, and in whom intrathecal contrast enhanced MRI was considered indicated for clinical reasons. Individuals who were not eligible for inclusion included subjects with a history of hypersensitivity reactions to contrast media agents, severe allergic reactions in general, evidence of renal dysfunction, i.e. glomerular filtration rate (GFR) < 30, age < 18 or > 80 years, or pregnant or breastfeeding women.
Patients were categorized according to tentative diagnosis prior to MRI, and underwent work-up, including blood sampling, prior to any treatment. The category reference subjects (REF) includes individuals in whom we found no apparent evidence of CSF disturbance and no indication for surgery. The group with spontaneous intracranial hypotension (SIH) had an identified CSF leakage that required surgery to close the leakage. The present subjects in the category idiopathic intracranial hypertension (IIH) were shunted and demonstrated clinical improvement thereafter. Patients with pineal cysts (PC) or arachnoid cysts (AC) underwent surgery with cyst removal and demonstrated post-operative clinical improvement. The category idiopathic normal pressure hydrocephalus (iNPH) included patients who based on clinical workup, imaging findings and results of intracranial pressure (ICP) monitoring [29, 30], underwent shunting with a demonstrated clinical improvement thereby qualifying for the diagnosis Definite iNPH according to the Japanese guidelines .
Intrathecal administration of gadobutrol
The MRI contrast agent gadobutrol (Gadovist™, Bayer Pharma AG, Berlin, Germany) was administered intrathecally in volumes of 0.10, 0.25 or 0.5 mL (1.0 mmol/mL) at a speed of a few seconds. The intrathecal injection procedure was done at the lumbar level. Correct entrance to the subarachnoid space was verified by CSF backflow from the spinal needle.
The first 80 patients received intrathecal gadobutrol in a dose of 0.50 mmol only, and the latter patients received intrathecal gadobutrol in alternatig doses of 0.10 mmol, 0.25 mmol or 0.5 mmol.
Quantification of gadolinium in blood
Venous blood samples were obtained at empirically determined regular time points up to about 48 h after intrathecal administration of gadobutrol, and were stored at -80 °C. Quantification of gadolinium to estimate concentrations of gadobutrol in blood and plasma was performed as previously described . In short, the whole blood samples were homogenized using an Ultra-Turrax homogenizer (IKA T18). Both plasma and the homogenized whole blood samples were subjected to digestion with ultrapure distilled nitric acid and deionized Milli-Q water in a closed-vessel microwave technique system (UltraCLAVE, Milestone, Italy). The samples were digested according to a 60-min stepwise heating program, with a maximum temperature of 250 °C held for 15 min. Following dilution, samples were analyzed for gadolinium by inductively coupled plasma mass spectrometry (Agilent 7700x, Agilent Technologies), employing indium at 0.1 μg/L as an internal standard. A 5-point standard curve (0.01–10 μg/L) was used. All analytical results were corrected for procedural blank values. Measured gadolinium concentrations were recalculated to gadobutrol concentrations.
In this work, both plasma and whole-blood gadobutrol were used. Linear regression through the origin was used to determine the plasma to whole blood ratio, and whole blood concentration of gadobutrol was interpolated to plasma concentrations for the purpose of pharmacokinetic modelling.
Gadobutrol population pharmacokinetic modelling
A population pharmacokinetic model was developed in order to determine individual pharmacokinetic parameters of intrathecally administered gadobutrol. A non-parametric adaptive grid approach implemented in Pmetrics (version 1.9.7) for R (version 4.0.0) was used . Based on available literature [33, 34] and previous work , both one- and two-compartment structural models were initially considered. The structural models provide the hypothesized framework for which transfer of gadobutrol occurs between compartments. Both the one- and two-compartmental models estimate the transfer of gadobutrol from CSF to blood, and elimination from blood. However, in the two-compartmental model, a peripheral tissue compartment was implemented, allowing gadobutrol to distribute into and from tissue. For the purpose of internal model validation, the dataset was split into a development- (80%) and validation-set (20%). Patients with more than six samples were allocated to the development set, with additional random allocation of profiles until 80% of total profiles. Model selection was primarily based on comparison of the relative root mean squared predictive error (RMSE, %) calculated from the relative predictive error of all gadobutrol concentrations in the development dataset. Additionally, the linear regression slope, R2-values of the observed versus predicted plots, Akaike’s information criteria (AIC) and the Bayesian information criteria (BIC) also guided model development to some extent. Covariates were not included, due to sole interest in individual predictions.
Posterior individual parameter values, as well as posterior individually predicted concentrations obtained from the final population pharmacokinetic model run with the complete dataset, were used for all pharmacokinetic calculations. Predictions were made in one-minute intervals from time of administration and up to 72 h. The following pharmacokinetic variables were evaluated:
The absorption half-life (T1/2, abs) is defined as the time for half the amount of gadobutrol in the CSF to be cleared to blood. This parameter was used as a surrogate marker for CSF to blood clearance of gadobutrol. T1/2, abs was calculated by dividing the natural logarithm of 2 over the model-estimated coefficient of absorption (Ka) from CSF to blood.
Time to maximum concentration (Tmax) in plasma and maximum concentration (Cmax) in plasma were obtained directly from the individual predictions.
Lag-time of absorption to blood (Tlag) is the model-estimated time for the tracer to reach the site of clearance in CSF. Longer Tlag thus implies that the molecule stays longer within CSF or that it takes longer before the clearance process to blood starts.
Area under the concentration–time curve from zero to infinity (AUC0-∞) was calculated with the trapezoidal approximation from individual posterior predicted concentrations using the ‘makeAUC’-function in the Pmetrics package for R. The AUC is a measure of systemic exposure of gadobutrol.
In order to compare parameters across multiple doses, Cmax and AUC0-∞ were normalized with respect to dose.
Comparisons between groups were performed using two-tailed individual samples t-test for continuous variables, and Fishers exact test for categorical variables. Values were visually assessed for normality prior to testing. Differences in parameters and normalized parameters between different doses were assessed using an analysis of variance. For the predefined α = 0.05, we considered 95% confidence intervals not including zero and P-values lower than 0.05 to be statistically significant.
The study included 161 patients, with a mean ± SD age of 54 ± 19 years (range 19 to 82 years), and with a mean body mass index (BMI) of 28 ± 5 kg/m2 (range 18 to 41 kg/m2). Patients were under clinical work-up for possible CSF disorders, with diagnosis categories as indicated in Table 1.
Several groups were statistically significantly different from the reference (REF) cohort, with respect to gender, age, height, body mass index (BMI), and kidney function. A total of 1,140 samples were analyzed for gadobutrol in plasma or whole blood; the mean number of samples was 8 ± 2 in each subject (range 1 to 11 samples).
Gadobutrol blood-to-plasma ratio
In 24 patients, 204 samples were concomitantly analyzed for gadobutrol in both plasma and whole blood. Concentration of gadobutrol in plasma was linearly associated (β = 1.795, R2adjusted = 0.997; P < 0.001) with whole blood concentration of gadobutrol (Fig. 1), and whole blood concentration of gadobutrol was interpolated to plasma concentrations for the purpose of pharmacokinetic modelling, using linear regression through the origin.
Population pharmacokinetic modeling
Both one- and two-compartment models were initially evaluated. Compared to the one-compartment structural model, a two-compartment model improved the goodness of fit. Furthermore, addition of an absorption lag-time improved the individual predictions, especially in the absorption phase.
The final population pharmacokinetic model consisted of two compartments with first-order transfer from CSF to blood, and first-order elimination from the central compartment (blood) and absorption lag-time, and the model described the data well (Fig. 2A). The final model ran on the complete dataset achieved a mean prediction error of -0.032, a root mean squared error of 0.283, and a percentage root mean squared error of 18.5%. Akaike’s Information Criteria (AIC) and Bayesian Information Criteria (BIC) were 667 and 702, respectively. When assessing residual error for different times, a trend for underprediction was shown during times between 5 and 15 h (Fig. 2B). No systematic trends were found when comparing residual error to the observed concentration of gadobutrol (Fig. 2C). Individual predictions for a random subset of patients are shown in Fig. 2D, demonstrating goodness of fit.
Mean pharmacokinetic profiles across intrathecally-administered doses of gadobutrol are shown in Fig. 3. No differences in neither absorption half-life, time to maximum concentration, nor dose-normalized maximum concentration were found across the administered doses of gadobutrol. However, a statistically significant difference in dose-normalized AUC0-∞ between the dose levels of 0.1 mmol and 0.5 mmol was found (Δ = − 5.22 [95% CI: − 9.68, − 0.77] μM h). Additionally, mean predictive error of the population pharmacokinetic model was not different between dose levels.
Inter-individual variability in gadobutrol CSF to blood clearance
Irrespective of diagnosis category, a large degree of inter-individual variability was observed with respect to the pharmacokinetic parameters of intrathecally administered gadobutrol. For the complete dataset, mean absorption half-life was 3.83 ± 2.50 h, with a coefficient of variation (CV) of 65%, which did not vary with dose. Time to maximum concentration (Tmax) and dose-normalized maximum concentration (Cmax) were 8.60 ± 4.58 h (CV 53%) and 0.69 ± 0.42 μM (CV 61%) respectively. The large inter-individual variability of pharmacokinetic parameters irrespective of diagnosis is shown in Fig. 4.
Disease categories show different profiles
A notable degree of variability in pharmacokinetic parameters was observed both within and between disease categories. Individual predicted profiles with group-wise mean predicted profiles are shown in Fig. 5, and pharmacokinetic parameters at group level with comparisons are presented in Table 2. Variability in mean concentration profiles of gadobutrol for the different patient groups is further presented in Fig. 6, illustrating the between group differences.
When compared with the reference cohort, patients with pineal cysts demonstrated a 0.46 [95% CI: 0.03, 0.88] hours longer absorption lag time (Table 2). In this group, several demographic factors were associated with the pharmacokinetic parameters (Fig. 7); Tmax and dose-normalized AUC0-Inf were positively associated with age, while dose-normalized Cmax was negatively associated with age, height and weight.
Neither patients with arachnoid cysts nor patients with spontaneous intracranial hypotension demonstrated any difference in pharmacokinetic parameters of intrathecally administered gadobutrol, compared with the reference cohort (Table 2).
In contrast, patients with idiopathic intracranial hypertension showed a 2.25 [95% CI: 0.74, 3.77] hours shorter absorption half-life when compared with the reference group, indicating a greater CSF to blood clearance of gadobutrol. Additionally, a 3.09 [95% CI: − 5.74, − 0.43] μM h greater dose-normalized AUC0-∞ was found compared with the reference cohort (Table 2).
In iNPH patients, compared with the reference cohort, time to maximum concentration was 2.36 [95% CI: 0.30, 4.41] hours longer, and showed a 5.91 [95% CI: 8.18, 3.63] μM h greater mean AUC0-∞. Additionally, the lag-time was 0.42 [95% CI: 0.09, 0.74] hours longer compared with reference (Table 2); hence, in iNPH the CSF tracer stays longer within the CSF compartment prior to clearance to blood.
Patients with communicating hydrocephalus demonstrated a 0.19 [95% CI: 0.02, 0.37] μM lower dose-normalized maximum concentration of gadobutrol when compared with the reference group, which was the lowest concentration measured in the included disease categories (Table 2). Even though no statistically significant differences between patients with non-communicating hydrocephalus and the reference cohort were found, most likely due to a low number of subjects in the aforementioned group, time to maximum concentration was numerically higher (12.33 ± 7.17 h), compared to the reference cohort (7.49 ± 4.09 h), as shown in Table 2.
In this work, we present a population pharmacokinetic model applied to intrathecally administered gadobutrol that precisely estimates the clearance from CSF to blood in patients with various diseases. The included patients showed a high degree of inter-individual variability in pharmacokinetic parameters both within and between different disease categories of CSF disturbances.
Up to now, the literature on CSF to blood clearance has been scarce. The presently described model is derived from 1,140 blood samples in 161 patients, referring to plasma levels of gadobutrol measured subsequently to intrathecal injections of predefined quantities. Utilizing positron emission tomography (PET), others  previously examined clearance of intrathecal 99mTc-DPTA (technetium-99-diethylene-triamine-pentaacetate) to urine. It also has been demonstrated reduced clearance of a PET ligand from cerebral ventricles to the nasal turbinate in Alzheimer patients . Furthermore, another recent PET study  showed reduced clearance of two PET tracers (18F-THK5351 and 11C-PiB) from ventricular CSF in patients with Alzheimer’s disease, providing support to the hypothesis that impaired clearance of amyloid-β from CSF underlies the amyloid cerebral deposition characterizing Alzheimer’s disease. However, with regard to PET, a drawback is that radioactive ligands provide a radiation dose to the individual , have short half-life (about 6 h for 99mTc-DPTA), and the diagnostic process is both expensive and time-consuming.
The most significant observation of the present study is the large inter-individual variation in CSF to blood clearance, as well as the differences between CSF disease. Compared to the reference cohort, patients diagnosed with pineal or arachnoid cysts, and to some degree patients with spontaneous intracranial hypotension, did not present any differences in pharmacokinetics of intrathecally administered gadobutrol. On the other hand, a statistically significant longer lag-time was found in patients with pineal cysts, but no difference in CSF to blood clearance was found. We conclude that on the group level, these categories may possibly reflect the normal variation.
Patients with idiopathic intracranial hypertension, on the other hand, demonstrated a significantly reduced absorption-half life, possibly indicating faster egress of molecules from CSF to blood due to increased ICP. Furthermore, in patients with iNPH, the time to maximum concentration was significantly longer compared to the reference group, and lag-time of absorption was significantly increased. Therefore, in patients with iNPH, the CSF tracer stays longer in the CSF compartment and it requires longer time to reach maximum concentration. The senior authors previously found evidence of reduced CSF turnover in iNPH [23, 39]. In iNPH patients, high grade ventricular reflux of tracer  may as well contribute to the increased lag time in these individuals. The CSF to blood clearance of gadobutrol per se was not affected at group level since the absorption half-life or maximum concentration was not different.
We may not from the present data decipher which transport routes gadobutrol follow from the CSF to blood. Emerging evidence points at the role of meningeal lymphatic vessels for molecular egress from CSF to blood, which is supported by findings of reduced clearance of neurotoxic metabolites from CSF when meningeal lymphatic clearance routes are impaired . In humans, the parasagittal dura may be a direct passage route to the meningeal lymphatic structures , though molecular efflux from CSF via the cribriform plate seems to be minor . Other possible efflux routes are the cranial and spinal nerve roots , and spinal lymphatic pathways . The arachnoid membrane itself has traditionally been considered impermeable to larger molecules . Hence, a CSF tracer study of mice found no signs of tracer propagation beyond the arachnoid layer . Traditionally, it has been thought that CSF egresses via arachnoid granulations to veins, but this view is up to debate . A microscopy study showed endothelial lined gaps and fissures in parasagittal dura of pigs, which might serve as a CSF drainage pathway . In humans, a subset of arachnoid granulations might drain CSF via lymphatic vessels to the venous circulation . Our group showed that the presently used intrathecal tracer gadobutrol enriched in parasagittal dura , bone marrow at the skull vertex adjacent or remote to intradiploic dural extensions , and in extracranial lymph nodes , and demonstrated the feasibility of measuring CSF to blood clearance . The time course of CSF clearance with peak in plasma after 8.60 ± 4.58 h may indicate a major role of the spinal canal given that tracer clearance from CSF peaks to blood occurred far earlier than peak enhancement at the skull vertex . Differences in lag time (Tlag) might be related to passage capacity within the thecal sac, but we have previously not found differences between groups for time between lumbar injection and first appearance at the cranio-cervical junction, i.e. spinal transit time . We suggest that the meningeal lymphatic vessels are the main route for egress of molecules from CSF, and that meningeal lymphatic impairment may hamper CSF clearance. In this regard, it is of particular interest that evidence from animal and human studies suggest the meningeal lymphatic function deteriorates with increasing age [45, 51], and that impaired meningeal lymphatic function aggravates pathology seen in animal models of Parkinson’s  and Alzheimer’s  diseases. Experimentally, it was shown that impaired meningeal lymphatic function reduced paravascular influx of macromolecules into the brain, and reduced efflux from the interstitial space . In comparison, we previously found in humans that peak CSF tracer enhancement in human brain and cervical lymph nodes concurred in time, supporting a role of meningeal lymphatic vessels in molecular drainage from CSF .
While plasma levels of gadobutrol primarily reflect clearance from CSF along extra-vascular pathways, a minor leakage of tracer through the BBB may to a limited extent contribute to the clearance as ageing as well as neurodegenerative disease may be accompanied with impaired BBB integrity [55, 56]. Evidence of BBB disruption has also been reported for CSF disease such as IIH  and iNPH . After entering to the blood, the plasma half-life of gadobutrol in blood is 1.5 h .
The present observations may have several clinical implications; we would like to highlight three areas. First, the present observations suggest that assessing CSF to blood clearance adds to characterization of CSF diseases on the individual level. One example is the identification of CSF leakage in individuals with spontaneous intracranial hypotension; it is well established that it may be very difficult to identify the site of CSF leakage . Currently, the visualization of CSF leakage utilizes MRI [60, 61], contrast enhanced computer tomography (CT) myelography  as well as intrathecal 99mTc-DPTA nuclear imaging , though the risk of not identifying any leakage site is high. A strategy to measure CSF to blood clearance might be expected to aid in identifying individuals with the most pronounced CSF leakage, even though signs of hyper-accelerated clearance could not be shown at group level for the CSF leakage sub-cohort in this study.
Second, direct measurements of CSF to blood clearance might prove useful in preclinical stages of neurodegenerative and dementia disease. Measurements of circulating substances such as amyloid-β and tau may be used for screening purposes, providing an indicative risk of disease . However, direct measures of CSF to blood clearance may be useful in a subset of individuals at risk. In this regard, it should be noted that about ¼ of amyloid-β is cleared via CSF in rodents [63, 64], and a significant amount of tau is excreted via CSF as the majority does not pass across the BBB. For example, mice without dural lymphatic drainage showed significantly reduced excretion of tau , and demonstrated a significant association between blood and CSF levels of tau . We here found that the dementia subtype iNPH was characterized with altered pharmacokinetic variables, including longer time to maximum concentration (Tmax), longer lag time (Tlag) and higher AUC, as compared with reference subjects. However, the difference in AUC may be attributed to the difference in age and renal function compared to the reference cohort.
Third, estimation of CSF to blood clearance may be useful preceding intrathecal drug administration for treatment of neurological disease. Even though it was traditionally thought that a substance within the CSF only passed a few millimeters into the cortical substance , intrathecally injected gadobutrol showed brain-wide distribution in humans . Therefore, intrathecal drugs may directly access the entire extra-vascular part of the CNS in contrast to systemically administered substances that are restricted by the BBB . This, however, may depend on the physiochemical properties of drugs. Examples of intrathecal drugs are antisense oligonucleotides [20, 66], such as Spinraza used for spinal muscular atrophy [17, 67], intrathecal chemotherapy, e.g. methotrexate, used for cancer [68, 69], and adeno-associated viral vector-mediated gene-delivery to CNS in amyotrophic lateral sclerosis, dementia disease and spinocerebellar ataxia [16, 70,71,72,73]. However, given the high degree of variation in CSF to blood clearance, there is risk of both over- and under-dosage.
Some limitations should be noted. Gadobutrol is administered off-label as it is not approved for intrathecal use. However, here we used gadobutrol in intrathecal doses of 0.10, 0.25 and 0.50 mmol, which have all been proven safe [50, 74]. Toxic effects have previously not been reported for intrathecal gadobutrol in doses of 1.0 mmol or below . We established dose linearity for the range of 0.10–0.50 mmol, and found no difference in the predictive performance of the population pharmacokinetic model between dose levels. As such, for estimating CSF to blood clearance with population pharmacokinetic modelling, an intrathecal dose of 0.10 mmol appears sufficient. Intrathecal gadobutrol is detected in blood with high sensitivity and accuracy; the present detection threshold was about 1.35 nM, well below the observed concentrations, rendering for use of even lower doses. As gadobutrol shares many of the same molecular properties with radiopaque contrast agents, where many are approved for intrathecal use, utility of on-label contrast agents intrathecally for CSF clearance assessment could be explored in future studies.
In this work, the less tangible absorption half-life was used as a surrogate marker for CSF to blood clearance, instead of actual clearance, due to the lack of accurate determinations of individual CSF volume. However, this does not affect the interpretation or accuracy of the results. With regard to the possible normal CSF to blood clearance in healthy people, it may as well be considered a limitation that we included a range of patients spanning multiple defined CSF disturbances. It was, however, beyond the scope of this work to discuss in detail the underlying diagnoses and the clinical significance of each disease category. Additional work on the subject would benefit from the inclusion of individuals without evident neurological disorders, in order to establish a reference value and level of variability in a healthy population. Furthermore, it remains to be determined whether gadobutrol is a valid marker for clearance of other intrathecally administered drugs and endogenous metabolites of interest in disease such as amyloid-β, tau, and α-synuclein.
In conclusion, this work provides a population pharmacokinetic model of CSF to blood clearance based on 1,140 blood samples from 161 subjects. Our data demonstrates a large degree of inter-individual variability in CSF to blood clearance as well as different clearance profiles across disease categories. CSF clearance function might both be a secondary feature of various neurological diseases, and a primary driver behind disease. As such, extensive clearance may characterize CSF leakage and spontaneous intracranial hypotension, while protracted clearance may be a contributing factor in neurodegenerative diseases. In the therapeutic setting, CSF to blood clearance may prove useful for tailoring dosage of intrathecal drugs, an administration route with prospects of increased utility in the near future.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author, upon reasonable request.
Area under the curve
Akaike’s Information Criteria
Blood brain barrier
Bayesian Information Criteria
Body mass index
- Cmax :
Central nervous system
- Comm HC:
Coefficient of variations
Glomerular filtration rate
Idiopathic intracranial hypertension
Idiopathic normal pressure hydrocephalus
Magnetic resonance imaging
- Non-comm HC:
Positron emission tomography
Root mean squared predictive error
Spontaneous intracranial hypotension
- T1/2, abs :
Absorption half life
Lag-time (of absorption)
- Tmax :
Time to maximum concentration
Levey AS, Coresh J, Balk E, Kausz AT, Levin A, Steffes MW, Hogg RJ, Perrone RD, Lau J, Eknoyan G. National Kidney Foundation practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Ann Intern Med. 2003;139:137–47.
Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, Benveniste H, Vates GE, Deane R, Goldman SA, et al. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid beta. Sci Transl Med. 2012;4:147ra111.
Louveau A, Smirnov I, Keyes TJ, Eccles JD, Rouhani SJ, Peske JD, Derecki NC, Castle D, Mandell JW, Lee KS, et al. Structural and functional features of central nervous system lymphatic vessels. Nature. 2015;523:337–41.
Nedergaard M, Goldman SA. Glymphatic failure as a final common pathway to dementia. Science. 2020;370:50–6.
Louveau A, Plog BA, Antila S, Alitalo K, Nedergaard M, Kipnis J. Understanding the functions and relationships of the glymphatic system and meningeal lymphatics. J Clin Invest. 2017;127:3210–9.
Bolte AC, Dutta AB, Hurt ME, Smirnov I, Kovacs MA, McKee CA, Ennerfelt HE, Shapiro D, Nguyen BH, Frost EL, et al. Meningeal lymphatic dysfunction exacerbates traumatic brain injury pathogenesis. Nat Commun. 2020;11:4524.
Song E, Mao T, Dong H, Boisserand LSB, Antila S, Bosenberg M, Alitalo K, Thomas JL, Iwasaki A. VEGF-C-driven lymphatic drainage enables immunosurveillance of brain tumours. Nature. 2020;577:689–94.
Hu X, Deng Q, Ma L, Li Q, Chen Y, Liao Y, Zhou F, Zhang C, Shao L, Feng J, et al. Meningeal lymphatic vessels regulate brain tumor drainage and immunity. Cell Res. 2020;30:229–43.
Ma Q, Schlegel F, Bachmann SB, Schneider H, Decker Y, Rudin M, Weller M, Proulx ST, Detmar M. Lymphatic outflow of cerebrospinal fluid is reduced in glioma. Sci Rep. 2019;9:14815.
Yanev P, Poinsatte K, Hominick D, Khurana N, Zuurbier KR, Berndt M, Plautz EJ, Dellinger MT, Stowe AM. Impaired meningeal lymphatic vessel development worsens stroke outcome. J Cereb Blood Flow Metab. 2020;40:263–75.
Chen J, Wang L, Xu H, Xing L, Zhuang Z, Zheng Y, Li X, Wang C, Chen S, Guo Z, et al. Meningeal lymphatics clear erythrocytes that arise from subarachnoid hemorrhage. Nat Commun. 2020;11:3159.
Da Mesquita S, Papadopoulos Z, Dykstra T, Brase L, Farias FG, Wall M, Jiang H, Kodira CD, de Lima KA, Herz J, et al. Meningeal lymphatics affect microglia responses and anti-Aβ immunotherapy. Nature. 2021;593:255–60.
Ding XB, Wang XX, Xia DH, Liu H, Tian HY, Fu Y, Chen YK, Qin C, Wang JQ, Xiang Z, et al. Impaired meningeal lymphatic drainage in patients with idiopathic Parkinson’s disease. Nat Med. 2021;27:411–8.
Rissman RA, Trojanowski JQ, Shaw LM, Aisen PS. Longitudinal plasma amyloid beta as a biomarker of Alzheimer’s disease. J Neural Transm (Vienna). 2012;119:843–50.
Blennow K, Hampel H, Weiner M, Zetterberg H. Cerebrospinal fluid and plasma biomarkers in Alzheimer disease. Nat Rev Neurol. 2010;6:131–44.
Chen W, Hu Y, Ju D. Gene therapy for neurodegenerative disorders: advances, insights and prospects. Acta Pharm Sin B. 2020;10:1347–59.
Wurster CD, Winter B, Wollinsky K, Ludolph AC, Uzelac Z, Witzel S, Schocke M, Schneider R, Kocak T. Intrathecal administration of nusinersen in adolescent and adult SMA type 2 and 3 patients. J Neurol. 2019;266:183–94.
Petrescu GED, Sabo AA, Torsin LI, Calin GA, Dragomir MP. MicroRNA based theranostics for brain cancer: basic principles. J Exp Clin Cancer Res. 2019;38:231.
Jadhav S, Avila J, Schöll M, Kovacs GG, Kövari E, Skrabana R, Evans LD, Kontsekova E, Malawska B, de Silva R, et al. A walk through tau therapeutic strategies. Acta Neuropathol Commun. 2019;7:22.
Chen Y, Mazur C, Luo Y, Sun L, Zhang M, McCampbell A, Tomassy GS. Intrathecal delivery of antisense oligonucleotides in the rat central nervous system. J Vis Exp. 2019. https://doi.org/10.3791/60274.
Pardridge WM. The blood-brain barrier: bottleneck in brain drug development. NeuroRx J Am Soc Exp NeuroTher. 2005;2:3–14.
Ringstad G, Valnes LM, Dale AM, Pripp AH, Vatnehol SS, Emblem KE, Mardal KA, Eide PK. Brain-wide glymphatic enhancement and clearance in humans assessed with MRI. JCI Insight. 2018;3:1–16.
Ringstad G, Vatnehol SAS, Eide PK. Glymphatic MRI in idiopathic normal pressure hydrocephalus. Brain. 2017;140:2691–705.
Ringstad G, Eide PK. Cerebrospinal fluid tracer efflux to parasagittal dura in humans. Nat Commun. 2020;11:1–9.
Ringstad G, Eide PK. Molecular trans-dural efflux to skull bone marrow in humans with cerebrospinal fluid disorders. Brain. 2021;145(4):1464–72.
Eide PK, Vatnehol SAS, Emblem KE, Ringstad G. Magnetic resonance imaging provides evidence of glymphatic drainage from human brain to cervical lymph nodes. Sci Rep. 2018;8:7194.
Eide PK, Mariussen E, Uggerud H, Pripp AH, Lashkarivand A, Hassel B, Christensen H, Hovd MH, Ringstad G. Clinical application of intrathecal gadobutrol for assessment of cerebrospinal fluid tracer clearance to blood. JCI insight. 2021;6:1–13.
Patel TK, Habimana-Griffin L, Gao X, Xu B, Achilefu S, Alitalo K, mckee CA, Sheehan PW, Musiek ES, Xiong C, et al. Dural lymphatics regulate clearance of extracellular tau from the CNS. Mol Neurodegener. 2019;14:11.
Eide PK, Sorteberg W. Diagnostic intracranial pressure monitoring and surgical management in idiopathic normal pressure hydrocephalus: a 6-year review of 214 patients. Neurosurgery. 2010;66:80–91.
Eide PK, Sorteberg W. Outcome of surgery for idiopathic normal pressure hydrocephalus: role of preoperative static and pulsatile intracranial pressure. World Neurosurg. 2016;86:186–93.
Mori E, Ishikawa M, Kato T, Kazui H, Miyake H, Miyajima M, Nakajima M, Hashimoto M, Kuriyama N, Tokuda T, et al. Guidelines for management of idiopathic normal pressure hydrocephalus: second edition. Neurol Med Chir (Tokyo). 2012;52:775–809.
Neely MN, van Guilder MG, Yamada WM, Schumitzky A, Jelliffe RW. Accurate detection of outliers and subpopulations with Pmetrics, a nonparametric and parametric pharmacometric modeling and simulation package for R. Ther Drug Monit. 2012;34:467–76.
Hahn G, Sorge I, Gruhn B, Glutig K, Hirsch W, Bhargava R, Furtner J, Born M, Schröder C, Ahlström H, et al. Pharmacokinetics and safety of gadobutrol-enhanced magnetic resonance imaging in pediatric patients. Invest Radiol. 2009;44:776–83.
Kunze C, Mentzel HJ, Krishnamurthy R, Fleck R, Stenzel M, Bhargava R, Burrowes D, Sutter G, Schultze-Mosgau M, Santiuste M, Hahn G. Pharmacokinetics and safety of macrocyclic gadobutrol in children aged younger than 2 years including term newborns in comparison to older populations. Invest Radiol. 2016;51:50–7.
Verma A, Hesterman JY, Chazen JL, Holt R, Connolly P, Horky L, Vallabhajosula S, Mozley PD. Intrathecal (99m)Tc-DTPA imaging of molecular passage from lumbar cerebrospinal fluid to brain and periphery in humans. Alzheimers Dement (Amst). 2020;12: e12030.
de Leon MJ, Li Y, Okamura N, Tsui WH, Saint-Louis LA, Glodzik L, Osorio RS, Fortea J, Butler T, Pirraglia E, et al. cerebrospinal fluid clearance in Alzheimer disease measured with dynamic PET. J Nucl Med. 2017;58:1471–6.
Li Y, Rusinek H, Butler T, Glodzik L, Pirraglia E, Babich J, Mozley PD, Nehmeh S, Pahlajani S, Wang X, et al. Decreased CSF clearance and increased brain amyloid in Alzheimer’s disease. Fluids Barriers CNS. 2022;19:21.
Hesterman JY, Kost SD, Holt RW, Dobson H, Verma A, Mozley PD. Three-dimensional dosimetry for radiation safety estimates from intrathecal administration. J Nucl Med. 2017;58:1672–8.
Eide PK, Ringstad G. Delayed clearance of cerebrospinal fluid tracer from entorhinal cortex in idiopathic normal pressure hydrocephalus: a glymphatic magnetic resonance imaging study. J Cereb Blood Flow Metab. 2019;39:1355–68.
Eide PK, Pripp AH, Ringstad G. Magnetic resonance imaging biomarkers of cerebrospinal fluid tracer dynamics in idiopathic normal pressure hydrocephalus. Brain Commun. 2020;2:1–16.
Melin E, Eide PK, Ringstad G. In vivo assessment of cerebrospinal fluid efflux to nasal mucosa in humans. Sci Rep. 2020;10:1–10.
Pollay M. The function and structure of the cerebrospinal fluid outflow system. Cerebrospinal Fluid Res. 2010;7:9.
Jacob L, Boisserand LSB, Geraldo LHM, de Brito NJ, Mathivet T, Antila S, Barka B, Xu Y, Thomas JM, Pestel J, et al. Anatomy and function of the vertebral column lymphatic network in mice. Nat Commun. 2019;10:4594.
Weller RO. Microscopic morphology and histology of the human meninges. Morphologie. 2005;89:22–34.
Ma Q, Ineichen BV, Detmar M, Proulx ST. Outflow of cerebrospinal fluid is predominantly through lymphatic vessels and is reduced in aged mice. Nat Commun. 2017;8:1434.
Brinker T, Stopa E, Morrison J, Klinge P. A new look at cerebrospinal fluid circulation. Fluids Barriers CNS. 2014;11:10.
Kutomi O, Takeda S. Identification of lymphatic endothelium in cranial arachnoid granulation-like dural gap. Microscopy (Oxf). 2020;69(6):391–400.
Yağmurlu K, Sokolowski J, Soldozy S, Norat P, Çırak M, Tvrdik P, Shaffrey ME, Kalani MYS: A subset of arachnoid granulations in humans drain to the venous circulation via intradural lymphatic vascular channels. J Neurosurg 2021:1–10.
Eide PK, Vatnehol SAS, Emblem KE, Ringstad G. Magnetic resonance imaging provides evidence of glymphatic drainage from human brain to cervical lymph nodes. Sci Rep. 2018;8:1–10.
Halvorsen M, Edeklev CS, Fraser-Green J, Lovland G, Vatnehol SAS, Gjertsen O, Nedregaard B, Sletteberg R, Ringstad G, Eide PK. Off-label intrathecal use of gadobutrol: safety study and comparison of administration protocols. Neuroradiology. 2021;63:51–61.
Zhou Y, Cai J, Zhang W, Gong X, Yan S, Zhang K, Luo Z, Sun J, Jiang Q, Lou M. Impairment of the Glymphatic pathway and putative meningeal lymphatic vessels in the aging human. Ann Neurol. 2020;87:357–69.
Zou W, Pu T, Feng W, Lu M, Zheng Y, Du R, Xiao M, Hu G. Blocking meningeal lymphatic drainage aggravates Parkinson’s disease-like pathology in mice overexpressing mutated alpha-synuclein. Transl Neurodegeneration. 2019;8:7.
Wang L, Zhang Y, Zhao Y, Marshall C, Wu T, Xiao M. Deep cervical lymph node ligation aggravates AD-like pathology of APP/PS1 mice. Brain Pathol. 2019;29:176–92.
Da Mesquita S, Louveau A, Vaccari A, Smirnov I, Cornelison RC, Kingsmore KM, Contarino C, Onengut-Gumuscu S, Farber E, Raper D, et al. Functional aspects of meningeal lymphatics in ageing and Alzheimer’s disease. Nature. 2018;560:185–91.
Montagne A, Barnes SR, Sweeney MD, Halliday MR, Sagare AP, Zhao Z, Toga AW, Jacobs RE, Liu CY, Amezcua L, et al. Blood-brain barrier breakdown in the aging human hippocampus. Neuron. 2015;85:296–302.
Montagne A, Nation DA, Sagare AP, Barisano G, Sweeney MD, Chakhoyan A, Pachicano M, Joe E, Nelson AR, D’Orazio LM, et al. APOE4 leads to blood-brain barrier dysfunction predicting cognitive decline. Nature. 2020;581:71–6.
Hasan-Olive MM, Hansson HA, Enger R, Nagelhus EA, Eide PK. Blood-brain barrier dysfunction in idiopathic intracranial hypertension. J Neuropathol Exp Neurol. 2019;78:808–18.
Eide PK, Hansson HA. Blood-brain barrier leakage of blood proteins in idiopathic normal pressure hydrocephalus. Brain Res. 2020;1727:1–13.
Staks T, Schuhmann-Giampieri G, Frenzel T, Weinmann HJ, Lange L, Platzek J. Pharmacokinetics, dose proportionality, and tolerability of gadobutrol after single intravenous injection in healthy volunteers. Invest Radiol. 1994;29:709–15.
Farb RI, Nicholson PJ, Peng PW, Massicotte EM, Lay C, Krings T, terBrugge KG. Spontaneous intracranial hypotension: a systematic imaging approach for CSF leak localization and management based on mri and digital subtraction myelography. AJNR Am J Neuroradiol. 2019;40:745–53.
Rahman M, Bidari SS, Quisling RG, Friedman WA. Spontaneous intracranial hypotension: dilemmas in diagnosis. Neurosurgery. 2011;69:4–14 (discussion 14).
Novotny C, Pötzi C, Asenbaum S, Peloschek P, Suess E, Hoffmann M. SPECT/CT fusion imaging in radionuclide cisternography for localization of liquor leakage sites. J Neuroimaging. 2009;19:227–34.
Roberts KF, Elbert DL, Kasten TP, Patterson BW, Sigurdson WC, Connors RE, Ovod V, Munsell LY, Mawuenyega KG, Miller-Thomas MM, et al. Amyloid-β efflux from the central nervous system into the plasma. Ann Neurol. 2014;76:837–44.
Feng W, Zhang Y, Wang Z, Xu H, Wu T, Marshall C, Gao J, Xiao M. Microglia prevent beta-amyloid plaque formation in the early stage of an Alzheimer’s disease mouse model with suppression of glymphatic clearance. Alzheimer’s Res Ther. 2020;12:125.
Pardridge WM. Drug transport in brain via the cerebrospinal fluid. Fluids Barriers CNS. 2011;8:7.
Mazur C, Powers B, Zasadny K, Sullivan JM, Dimant H, Kamme F, Hesterman J, Matson J, Oestergaard M, Seaman M, et al. Brain pharmacology of intrathecal antisense oligonucleotides revealed through multimodal imaging. JCI Insight. 2019. https://doi.org/10.1172/jci.insight.129240.
Finkel RS, Chiriboga CA, Vajsar J, Day JW, Montes J, De Vivo DC, Yamashita M, Rigo F, Hung G, Schneider E, et al. Treatment of infantile-onset spinal muscular atrophy with nusinersen: a phase 2, open-label, dose-escalation study. Lancet. 2016;388:3017–26.
Byrnes DM, Vargas F, Dermarkarian C, Kahn R, Kwon D, Hurley J, Schatz JH. Complications of intrathecal chemotherapy in adults: single-institution experience in 109 consecutive patients. J Oncol. 2019;2019:4047617.
Livshits Z, Rao RB, Smith SW. An approach to chemotherapy-associated toxicity. Emerg Med Clin North Am. 2014;32:167–203.
Nguyen AD, Nguyen TA, Zhang J, Devireddy S, Zhou P, Karydas AM, Xu X, Miller BL, Rigo F, Ferguson SM, et al. Murine knockin model for progranulin-deficient frontotemporal dementia with nonsense-mediated mRNA decay. Proc Natl Acad Sci U S A. 2018;115:E2849-e2858.
McCampbell A, Cole T, Wegener AJ, Tomassy GS, Setnicka A, Farley BJ, Schoch KM, Hoye ML, Shabsovich M, Sun L, et al. Antisense oligonucleotides extend survival and reverse decrement in muscle response in ALS models. J Clin Invest. 2018;128:3558–67.
Niu C, Prakash TP, Kim A, Quach JL, Huryn LA, Yang Y, Lopez E, Jazayeri A, Hung G, Sopher BL, et al. Antisense oligonucleotides targeting mutant Ataxin-7 restore visual function in a mouse model of spinocerebellar ataxia type 7. Sci Transl Med. 2018. https://doi.org/10.1126/scitranslmed.aap8677.
Hardcastle N, Boulis NM, Federici T. AAV gene delivery to the spinal cord: serotypes, methods, candidate diseases, and clinical trials. Expert Opin Biol Ther. 2018;18:293–307.
Edeklev CS, Halvorsen M, Lovland G, Vatnehol SAS, Gjertsen O, Nedregaard B, Sletteberg R, Ringstad G, Eide PK. Intrathecal use of gadobutrol for glymphatic MR imaging: prospective safety study of 100 patients. AJNR Am J Neuroradiol. 2019;40:1257–64.
Patel M, Atyani A, Salameh JP, McInnes M, Chakraborty S. Safety of intrathecal administration of gadolinium-based contrast agents: a systematic review and meta-analysis. Radiology. 2020. https://doi.org/10.1148/radiol.2020191373.
The authors thank dr. Øivind Gjertsen, dr. Bård Nedregaard and dr. Ruth Sletteberg from the Department of Radiology, Oslo University Hospital – Rikshospitalet, who performed some of the intrathecal gadobutrol injection procedures. In addition, the authors thank Marit Vadset from NILU-Norwegian Institute for Air Research, who analyzed a large part of the blood samples.
The work was supported by Department of Neurosurgery, Oslo university hospital, and Norwegian Institute for Air Research, Kjeller, Norway, and by the University of Oslo.
Ethics approval and consent to participate
This present study was approved by The Regional Committee for Medical and Health Research Ethics (REK) of Health Region South-East, Norway (2015/96), the Institutional Review Board of Oslo university hospital (2015/1868), and the National Medicines Agency of Norway (15/04932–7). Participants were included after providing written and oral informed consent.
Consent for publication
Geir Ringstad received a fee for speaking at the Bayer symposium at the European Congress of Radiology 2020 (Vienna, Austria). Geir Ringstad and Per Kristian Eide also have a patent pending. The other authors disclose no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hovd, M.H., Mariussen, E., Uggerud, H. et al. Population pharmacokinetic modeling of CSF to blood clearance: prospective tracer study of 161 patients under work-up for CSF disorders. Fluids Barriers CNS 19, 55 (2022). https://doi.org/10.1186/s12987-022-00352-w
- Cerebrospinal fluid
- Brain metabolites
- Intrathecal administration
- Intrathecal drugs