 Research
 Open Access
 Published:
Are standing osmotic gradients the main driver of cerebrospinal fluid production? A computational analysis
Fluids and Barriers of the CNS volume 20, Article number: 18 (2023)
Abstract
Background
The mechanisms of cerebrospinal fluid (CSF) production by the ventricular choroid plexus (ChP) have not been fully deciphered. One prominent hypothesized mechanism is transepithelial water transport mediated by accumulation of solutes at the luminal ChP membrane that produces local osmotic gradients. However, this standing osmotic gradient hypothesis has not been systematically tested.
Methods
To assess the plausibility of the standing gradient mechanism serving as the main driver of CSF production by the ChP, we developed a threedimensional (3D) and a onedimensional (1D) computational model to quantitatively describe the associated processes in the rat ChP intermicrovillar spaces and in CSF pools between macroscopic ChP folds (1D only). The computationally expensive 3D model was used to examine the applicability of the 1D model for hypothesis testing. The 1D model was employed to predict the rate of CSF produced by the standing gradient mechanism for 200,000 parameter permutations. Model parameter values for each permutation were chosen by random sampling from distributions derived from published experimental data.
Results
Both models predict that the CSF production rate by the standing osmotic gradient mechanism is below 10% of experimentally measured values that reflect the contribution of all actual production mechanisms. The 1D model indicates that increasing the size of CSF pools between ChP folds, where diffusion dominates solute transport, would increase the contribution of the standing gradient mechanism to CSF production.
Conclusions
The models suggest that the effect of standing osmotic gradients is too small to contribute substantially to CSF production. ChP motion and movement of CSF in the ventricles, which are not accounted for in the models, would further reduce this effect, making it unlikely that standing osmotic gradients are the main drivers of CSF production.
Introduction
Cerebrospinal fluid (CSF) is produced primarily by the choroid plexus (ChP) of the four cerebral ventricles [1]. There is no consensus on how CSF is generated, but several mechanisms have been proposed, including hydrostatic pressure difference between blood and CSF [2], water transport through bloodCSF barrier tight junctions (by mechanodiffusion [3], electrodiffusion [4], or claudin mediation [5]), or by molecular ion transportermediated water translocation, e.g., through the Na^{+}K^{+}2Cl^{−} cotransporter 1 (NKCC1), Na^{+}HCO_{3}^{−} cotransporter (NBCe2), and Na^{+}K^{+}ATPase [6]. However, a prevalent hypothesis is that the ChP secretes CSF by osmosis, with water following an osmotic gradient produced by the transport of ions across the ChP epithelium [7]. This mechanism requires sufficient water permeability, along with either a global osmotic gradient across the ChP or local gradients over the epithelium. Global gradients are either small or nonexistent, with CSF reported to be either slightly hyperosmolar with respect to blood (approximately 5 mOsm) [8,9,10,11] or isosmolar [6, 12].
The origin of the idea of local osmotic gradients as the driving force behind CSF production can be traced back to experiments on water transport across the intestinal epithelium [13]. Water and solute transport from the luminal (gut) side to the serosal (blood) side was found to be a function of, primarily, cellular activity rather than a bulk osmotic pressure difference between the two sides [13]. To explain this, Diamond and Bossert proposed what came to be known as the standing gradient (SG) model [14], which was based on the work of Curran and MacIntosh who had proposed a threecompartment model of epithelial water transport [15].
Diamond and Bossert considered the intracellular and the lateral intercellular spaces as the first and the second compartment, respectively [16]. They assumed these to be separated by a selectively permeable membrane, through which solutes are actively transported from the first to the second compartment, where they produce a locally elevated solute concentration level. The ensuing local osmotic gradient across the membrane then drives water from the first to the second compartment, increasing the hydrostatic pressure therein. Consequently, the solution is transferred through an imaginary, fully permeable membrane to the third compartment, which is the basal tissue between the epithelial cells and the capillary blood vessel. The tissue was considered to have the same bulk osmolarity as the epithelial intracellular space.
Unlike the intestinal epithelium, which is “absorptive” in nature, the ChP epithelium is “secretory,” meaning that for ChP, net water flux is from the basolateral to the luminal side, i.e., from blood to the CSF space. While in absorptive epithelia a local increase in osmolyte concentration in the lateral intercellular space may serve as a driving force for water transport, in secretory epithelia, such a concentration gradient would oppose it. With the epithelial tight junctions localized on the luminal side in the ChP, the lateral intercellular spaces are on the ‘wrong’ side to enable local osmotic gradientdriven fluid transfer [1]. However, having observed the movement of water from ChP against a bulk concentration gradient, Pollay suggested that a favorable gradient could still be present, produced by increased osmolyte concentration at the luminal membrane of the ChP between microvilli, which are fingerlike protrusions of the plasma membrane [17]. Thus, local osmotic gradients in the ChP intermicrovillar spaces could be drivers of CSF production.
Since direct experimental assessment of the SG hypothesis is difficult, mathematical models have been used to aid in its evaluation. Most of this work has been done for absorptive epithelia: for example, Sackin and Boulpaep [18] suggested that slight osmolarity differences (below the limit of detection) between interstitial fluid and blood in Necturus renal proximal tubule could be responsible for the observed quasiisotonic water reabsorption. Consequently, there would be no need for an osmotic gradient in the lateral intercellular space. Schafer et al. [19] came to the same conclusion with their model for the rabbit proximal tubule. Hill argued that an abnormally high permeability would be required for SGbased isosmotic water transport in leaky epithelia [20, 21]. However, Diamond argued that these studies do not provide direct indications against the SG mechanism and called for further experimental evidence [22]. Pedley and Fischbarg [23] used a model to interpret the experimental results of Wright et al. [24] on water flux across the rabbit gallbladder epithelium, concluding that the experimental data were inconsistent with the SG model. The structural difference between the secretory ChP epithelium and absorptive epithelia requires a separate assessment of the plausibility of the SG hypothesis for the brain. We recently showed using a 1D computational model that—with the employed set of parameter values—the SG mechanism in the intermicrovillar space could not account for the CSF produced by the ChP [6]. However, the employed model did not take into account the possible effect of stagnant pools of CSF within ChP infoldings.
One aspect common to the mentioned models is their high sensitivity to input parameters and boundary conditions such as solute–solvent transport coupling at cellular membranes. This has blunted the impact of the corresponding studies on the scientific discourse. Therefore, models covering large anatomic and physiologic parameter ranges are required to allow for a more robust assessment of the SG hypothesis. In the current study, we tested whether standing osmotic gradients at the surface of the ChP epithelium could be the predominant factor for CSF production by the ChP in rats. To this end, we developed a onedimensional (1D) computational model to describe the conjugate fluid and solute transport within a subunit of intermicrovillar space (defined as a functional unit, FU) and between macroscopic ChP folds, under the assumption that standing osmotic gradients are the only mechanism contributing to CSF production. We then compared the calculated CSF production rate to measured rates of CSF production, which include the contributions of all actual mechanisms. We emphasize that, unlike the original SG model, we did not limit solute transport to the bottom 10% of the FU [14], but considered solute influx along the entire length of the intermicrovillar space. Given the wide range of values reported for the relevant anatomic and biophysical parameters and the unknown sensitivity of the CSF production rate to variations in these, we performed computations with approximately 200,000 parameter value permutations. To avoid bias against the SG hypothesis, we selected a general model setup that would over rather than underestimate the contribution of standing gradients to CSF production. Furthermore, to confirm the applicability of the 1D model for hypothesis testing, we also developed and deployed a separate threedimensional (3D) model of intermicrovillar fluid and solute transport.
Methods
Model domain
We simplified the ChP luminal membrane to a surface covered with homogeneously distributed cylindrical pins representing microvilli as shown in Fig. 1. To account for the convoluted surface of the ChP, which produces pools shielded from the CSF movement at the center of the ventricles, we defined a protected area extending from the tip of the microvilli to a distance of \({l}_{\text{prot}}\) away from them (illustrated in Figs. 1, 3). Within the protected area, only flow caused by local CSF production contributes to solute convection.
Since the microvilli are arranged in a recurring pattern, we defined a representative functional unit (FU) as the space between four adjacent microvilli (Fig. 1b, c). Microvillar radius and length are \({r}_{\text{mv}}\) and \({l}_{\text{mv}}\), respectively. The microvillar separation distance, \(p\), and the number of FUs, \(N\), are given by
and
respectively, where \(\sigma\) is the surface density of microvilli and \({A}_{\text{app}}\) is the apparent luminal surface area, i.e., the surface area of the ChP without extension by microvilli.
To derive the 1D model, we simplified the FU to a cylindrical channel with equivalent hydraulic diameter, \(d\), as shown in Fig. 1d, where
Here, \({A}_{\text{CS}}\) and \({P}_{\text{CS}}\) are the crosssectional area and perimeter, respectively, of the original FU.
Onedimensional standing gradient model
The transport of fluid and solutes in the FU are coupled. Solute transport is governed by diffusion and convection, whereas fluid transport is driven by transmembrane solute concentration differences. The fluid velocity \(u\left(z\right)\) and solute concentration \(C\left(z\right)\) along the longitudinal axis (zaxis) of the FU, as shown in Fig. 1d, are determined by Eqs. (4)–(8). This set of equations describes the conjugate transport of solute and water through the intermicrovillar space in the longitudinal direction.
Here, \(\phi \left(z\right)\) is the local solute flux, \(\rho\) is CSF density, and \(D\) is the solute diffusion coefficient in CSF. \({L}_{\text{p}}\) is the water permeability of the luminal membrane and \({C}_{0}\) is the bulk ventricular CSF osmolality. Equations (4) and (5) describe the solute mass balance in the FU and the osmotic water transport through the luminal FU surface, respectively. They correspond to the equations of the original standing gradient model [14]. The boundary conditions represented by Eqs. (6) and (7) indicate, respectively, that the base of the FU is impermeable to solutes and fluid. This is necessary because of the singledimensional characteristic of the model. However, since inflow through the FU base may occur in reality, solute flux through the lateral boundary is adapted to account for it (see Eq. (9)). Equation (8) represents the boundary condition at the FU tip, where the relative contributions of convective and diffusive transport of solutes are balanced as a function of the protected length (see Additional file 1). In the case of zero protected length, the concentration at the tip of the FU corresponds to bulk solute concentration, a condition considered in a previous model [6].
The solute flux used in Eq. (4) is calculated from the measured CSF production rate \({Q}_{\text{meas}}\) according to Eq. (9). While production may vary in time, experimentally determined \({Q}_{\text{meas}}\) and \({C}_{0}\) are timeaveraged quantities. Therefore, steadystate forms of the governing equations and boundary conditions were employed (Eqs. (4)–(8)).
Here, \({f}_{\text{b}}\) is a factor applied to the first \(d/2{l}_{\text{mv}}\) of the FU length to correct for the lack of solute flux from the FU base. We chose the length of the fluxcorrected section based on the FU aspect ratio to preempt potential numerical instabilities arising from a big jump in \(\phi (z)\). We note that the choice of a specific length has little effect on the calculated CSF production rate. \({f}_{\text{b}}\) has a value of 0.5 within the modified section (corresponding to the area ratio of the base to the side of this section) and zero everywhere else.
Threedimensional standing gradient model
The 3D model considers coupled fluid and solute transport, taking into account viscous forces that were not included in the 1D model. The fluid velocity and pressure and the solute concentration distributions are governed by Eqs. (10)–(15). This set of equations determines the same quantities as that underlying the 1D model, but does so taking into account CSF flow in all directions (rather than only in the longitudinal direction) in a more realistic FU geometry.
Here, \(\mathbf{V}=\left({u}_{x},{u}_{y},{u}_{z}\right)\) is the fluid velocity vector, \(P\) is fluid pressure, and \(\mu\) is the dynamic viscosity of CSF.
Equations (10) and (11) describe, respectively, mass conservation of the fluid in the FU and the balance between pressure and viscous forces on the fluid. Solute transport via advection and diffusion while conserving solute mass is described by Eq. (12). The transport of fluid and solute is assumed to take place under steadystate conditions. Fluid and solute enter the FU perpendicularly from the luminal surface (Fig. 1c), with the corresponding velocity and flux given by Eqs. (13) and (14), respectively. The solute flux in the 3D model does not require the correction factor \({f}_{\text{b}}\), since flux through the FU base is taken into the account through the boundary condition in Eq. (14). We note that while the overall amount of solute entering the FU is the same in the 1D and 3D models, the flux is different (Eq. (15)). This is because the surfaces through which the solutes enter are not identical. At the tip of the FU, CSF has zero gaugepressure and the same solute concentration as bulk ventricular CSF. Symmetry conditions are applied on boundaries to adjacent FUs, as the repetitive pattern of FU distribution implies zero net exchange of fluid and solute between them. The CSF production rate by the standing gradient mechanism is computed by integrating the velocity field on the entire ChP luminal surface.
Choice of model parameters
For verification of our 1D model against the original implementation of the standing gradient one [14], we employed the same model parameters as used in the original study by Diamond and Bossert (Table 1). We emphasize that this model parameter set was employed for comparison purposes only, since Diamond and Bossert did not provide results for secretory epithelia. Our subsequent calculations relied on parameters reflective of the ChP.
To test whether standing osmotic gradients in the ChP intermicrovillar spaces can account for most of the ventricular CSF production, we considered distributions of parameter values as illustrated in Fig. 2 and described in Table 2. This was done to account for uncertainties in the experimental determination of the respective parameters, as well as to reflect differences in experimental settings in cases where multiple sources of data were available. The distributions were established by assuming normal distributions and then calculating the mean and standard deviation of data reported in the literature. The distribution chosen for the diffusion coefficient of solutes in CSF is representative of Na^{+}, the most relevant ion for the SG mechanism. We computed the distribution of the ChP apparent area and luminal membrane permeability based on the distribution of other related parameters (see Additional file 1). A precise evaluation of the protected length distribution over the ChP luminal surface would necessitate a 3D in vivo scan of the entire rat ChP at singlecell resolution, followed by a cellbycell calculation of the distance from the epithelium to the bulk CSF region (see Fig. 3d for an illustration of different protected lengths). Given that such data are currently unavailable, we relied on the Waxholm Space Rat Brain Atlas (v4, RRID: SCR_017124) [25], on measurements of ChP mass and volume [26,27,28], and on a geometric analysis to estimate \({l}_{\text{prot}}\) as detailed in Additional file 1. This yielded an upper limit of 77.5 µm for the average protected length. To be conservative—longer \({l}_{\text{prot}}\) favors the SG mechanism—we chose this upper limit as the center value for a normal distribution that extends to 155 µm. For comparison of the 1D to the 3D model, we used three sets of parameter values (Table 3). One of these corresponds to the parameter set that yields the mean value of CSF production rate as per the 1D model with zero protected length. The other two sets correspond to those yielding, respectively, mean ± standard deviation production rate.
Numerical procedure
Equations (4) and (5) were solved in MATLAB R2020b using bvp4c, a fourthorder numerically accurate finite difference solver for systems of ordinary differential equations. Equations (10)–(12) were solved using the finitevolume computational fluid dynamics software ANSYS Fluent on unstructured grids consisting of approximately 131,000 to 172,000 cells (Table 3). A secondorder scheme was used for momentum and transport equation discretization. A grid independence study was performed to ensure that the calculated CSF production rate was not unduly influenced by the choice of the computational mesh.
To calculate the CSF production rate with the 1D model based on the parameter distributions described in Table 2, approximately 200,000 parameter value sets were randomly chosen while considering the likelihood of each choice. To this end, we used the Latin hypercube sampling technique [29], which generates random points within equal probability intervals of the considered distribution function, and thereby ensures proper coverage of the stochastic space. An independence study for the number of sampling points was performed to ensure that the calculated CSF production rate distribution did not change with more samples.
Choroid plexus imaging
For immunohistochemistry, anesthetized Sprague–Dawley rats (P21) were transcardially perfused with 4% paraformaldehyde, and the excised brain was immersed in the fixative at 4 °C overnight before embedding in paraffin blocks and sectioned in a microtome. Sections were deparaffinized and rehydrated in xylene and ethanol following standard protocols prior to labeling (primary antibody: antiNKCC1, 1:400, Abcam AB59791, secondary antibody: Alexa Fluor® Goat antirabbit IgG, 1:500, Life Tech A11034). Sections were mounted with ProLong Gold DAPI mounting medium (Dako) and imaged using a Zeiss LSM700 point laser (Argon Lasos RMC781272) scanning confocal microscope with a Zeiss PlanApochromat 63 × /1.4 numerical aperture oil immersion objective (Carl Zeiss, Oberkochen, Germany).
Results
To verify the implementation of the 1D SG model, we selected the same parameter values as used by Diamond and Bossert for their original standing gradient one [14]. Since the original model does not consider a protected region, we set the protected length to zero for this comparison. The solute flux was prescribed as a function of zlocation in the FU to match the conditions in [14]. Similarly, the hydraulic diameter of the FU was directly set. Figure 4 shows that the velocity and concentration distributions along the FU zaxis of the two models match very well. This indicates that Eqs. (4) and (5) with the boundary conditions (6)–(8) were solved consistently.
We assessed the applicability of the 1D SG model for testing the hypothesis that standing osmotic gradients are the main drivers of CSF production in rats by comparing its production rate predictions with those of the more intricate 3D SG model. The set of parameter values given in Table 3 was used to this end. Figure 5 shows the concentration and zvelocity distribution in the FU on two section planes as computed with the 3D model for the mean case. The concentration field appears largely onedimensional (varying only along the zaxis), while the velocity is distributed threedimensionally with the maximum at the center of the FU. Table 4 lists the production rates predicted by the 1D and 3D models, showing that the 1D model reports higher CSF production rates by the SG mechanism than the 3D model predicts.
After having verified the implementation of the 1D SG model and having shown that it likely over rather than underestimates CSF production by the SG mechanism, we employed it to calculate the CSF production rate based on the parameter distributions given in Table 2. Figure 6 shows the probability density functions of the experimentally measured CSF production rate, and of the production rate attributed to the SG mechanism as per the 1D model. The latter amounts to less than 10% of the measured total CSF production rate. We note that the experimental measurements reflect actual CSF production resulting from all contributing driving forces, while the model provides an estimate of production by the SG mechanism alone. We will argue in the next section that the calculated CSF production by SG can be considered an upper limit, and why its actual contribution to overall CSF production is likely lower.
Discussion
We aimed to test the hypothesis that standing osmotic gradients at the surface of the ChP epithelium are the main drivers of CSF production. To this end, we extended and reformulated the original standing gradient model, which was designed for absorptive epithelia, to reflect the conditions on the ChP epithelial surface. Upon verification of this new 1D SG model with the results of Diamond and Bossert [14], we compared its output with that of a new 3D SG model, showing that the 1D model attributes a higher CSF production rate to the SG mechanism for the same conditions. In other words, the 1D model likely overestimates the effect of standing osmotic gradients, which is of advantage for the purpose of plausibility testing. Finally, by performing approximately 200,000 calculations using the 1D model with different parameter value sets randomly selected based on the likelihood of each value, we found that the predicted CSF production rate by the SG mechanism, if it were the sole contributor to CSF production, is one order of magnitude below that of the actual CSF production rate. This suggests that potential intermicrovillar standing osmotic gradients, even when enhanced by macroscopic stagnant pools, are not sufficient to drive CSF production by the ChP.
Both models require the prescription of solute flux into the FU. Since solute flux into the intermicrovillar space has not been experimentally quantified, we assumed that all solutes in the ventricular space originate in intermicrovillar spaces, and that all solutes exit the ventricles driven by bulk CSF flow. Under steadystate conditions, ventricular CSF production is equal to the amount of CSF leaving the ventricles. Therefore, the cumulative solute flux from all FUs corresponds to the bulk solute concentration multiplied by the CSF production rate. This means that the actual CSF production rate is an input to the models, and the calculated production rate corresponds to the fraction that can be attributed to the SG mechanism. Since there is a wide range of measured CSF production rates in rats reported in the literature [6, 30,31,32,33], we repeated the probabilistic 1D model calculations for another value of the measured CSF production rate that marks the lower end of this range, namely 0.74 ± 0.05 µL/min [31]. We note that the nominal value used marks the upper end with 6.8 ± 0.3 µL/min (Table 2). The production rate calculated with the lower end was 0.07 ± 0.05 µL/min, thus still one order of magnitude smaller than the measured rate. This behavior is expected, since a reduction in the prescribed CSF production rate lowers solute flux into the FUs, which reduces transepithelial osmotic gradients. The prescription of solute flux into the FU in this manner likely leads to an overestimation of CSF production via standing gradients by both models, since the intermicrovillar spaces do not have to be the sole source of solutes.
Neither the 1D nor the 3D model accounts for ChP motion, which could be caused by cardiovascular and respiratory action or by head and spine movement. Any such motion would work against CSF production by the SG mechanism both by disturbing standing gradients in the intermicrovillar spaces and by flushing pools of protected CSF from ChP folds. In the 1D model, the effect of protected pools is considered by the protected length, which is defined purely based on anatomic considerations. To obtain more accurate predictions of the CSF production rate, the effective protected length should be considered instead, noting that its value is locationdependent and difficult to determine. For the purposes of this study, it is sufficient to note that the values of protected length considered here are higher than the corresponding effective values. To understand the impact of changes in protected length on CSF production, we simplified Eqs. (6)–(8) using a small Péclet number approximation (see Additional file 1), which yielded a linear relation between the production rate and protected length: an increase in the protected length while maintaining other parameters constant yields a corresponding linear increase in CSF production. Thus, while not actively contributing to CSF production, protected pools may amplify the effect of potential standing gradients in the intermicrovillar space. Since the protected length values used here are expected to be larger than the effective ones, the model likely overestimates CSF production by the SG mechanism.
The 1D model predicts a higher CSF production rate than the 3D model. This can be attributed to a difference in the area ratio of the FU solute input surface (active surface) to its total surface. In the 1D model, solutes enter the FU through the whole side surface (Fig. 1d), whereas in the 3D model, parts of the side surfaces are interfaces to neighboring FUs, i.e., no solute flux occurs through those. Consequently, solutes that diffuse in the circumferential direction and produce standing gradients away from the active surfaces cannot contribute to CSF production. This effect is not accounted for in the 1D model, where the ratio of active to total surface is equal to one. Table 4 shows that the overestimation of CSF production by the 1D model reduces as the active to total surface ratio in the 3D model increases. Consequently, it is reasonable to assume that the 1D model generally predicts a higher CSF production rate by the SG mechanism than the 3D model.
The transport processes in the ventricular space, and in particular around the ChP, are complex. They include both diffusive and advective modes, moving interfaces, and transient CSF dynamics. Neither the 1D nor the 3D model can capture the effect of all processes potentially relevant to the SG mechanism. However, the models are designed to overestimate CSF production by standing osmotic gradients. This is achieved by making simplifying assumptions on some of the model parameters: we supposed that all solutes in the ventricular space originate in the ChP intermicrovillar spaces, which leads to an overestimation of the standing gradients. Neglecting tissue and bulk fluid motion (that would disturb standing gradients) results in a larger protected length. Assuming an active to total surface ratio of one for the FU results in a high transmembrane water transport rate. Finally, considering the SG mechanism as the sole driving force for CSF production discounts the effects of other mechanisms, such as hydrostatic pressure gradients, that would reduce standing gradients. Consequently, the herereported contribution of the SG mechanism to ChP CSF production likely constitutes an upper limit.
Conclusion
We have implemented and verified a onedimensional standing osmotic gradient model of CSF production on the ChP luminal surface that also accounts for the effect of protected pools in ChP folds. Calculations with this model, based on probabilistic parameter value distributions derived from experimental measurements, suggest that potential local osmotic gradients in the intermicrovillar spaces are too small to contribute substantially to CSF production; this even though the underlying model assumptions, which include protected pools of CSF, favor the SG mechanism. ChP motion and movement of CSF in the ventricles, both not accounted for in the model, would reduce local osmotic gradients, making it unlikely that they are the main drivers of CSF production.
Availability of data and materials
Code and data supporting the conclusions of this study are available on request from the corresponding author.
Abbreviations
 1D:

Onedimensional
 3D:

Threedimensional
 CSF:

Cerebrospinal fluid
 ChP:

Choroid plexus
 FU:

Functional unit
 NBCe2:

Na^{+}HCO_{3}^{−} cotransporter
 NKCC1:

Na^{+}K^{+}2Cl^{−} cotransporter 1
 SD:

Standard deviation
 SG:

Standing gradient
References
Macaulay N, Keep RF, Zeuthen T. Cerebrospinal fluid production by the choroid plexus: a century of barrier research revisited. Fluids Barriers CNS. 2022;19(1):11–22.
Johanson CE, Duncan JA, Klinge PM, Brinker T, Stopa EG, Silverberg GD. Multiplicity of cerebrospinal fluid functions: new challenges in health and disease. Cerebrospinal Fluid Res. 2008;5(1):10.
Hill AE, ShacharHill B. A new approach to epithelial isotonic fluid transport: an osmosensor feedback model. J Membr Biol. 2006;210(2):77–90.
Fischbarg J. Fluid transport across leaky epithelia: central role of the tight junction and supporting role of aquaporins. Physiol Rev. 2010;90(4):1271–90.
Furuse M, Fujita K, Hiiragi T, Fujimoto K, Tsukita S. Claudin1 and 2: novel integral membrane proteins localizing at tight junctions with no sequence similarity to occludin. J Cell Biol. 1998;141(7):1539–50.
Oernbo EK, Steffensen AB, Razzaghi Khamesi P, ToftBertelsen TL, Barbuskaite D, Vilhardt F, Gerkau NJ, Tritsaris K, Simonsen AH, Lolansen SD, Andreassen SN, Hasselbalch SG, Zeuthen T, Rose CR, Kurtcuoglu V, MacAulay N. Membrane transporters control cerebrospinal fluid formation independently of conventional osmosis to modulate intracranial pressure. Fluids Barriers CNS. 2022;19(1):65.
Damkier HH, Brown PD, Praetorius J. Cerebrospinal fluid secretion by the choroid plexus. Physiol Rev. 2013;93(4):1847–92.
Bradbury MWB, Kleeman CR. The effect of chronic osmotic disturbance on the concentrations of cations in cerebrospinal fluid. J Physiol. 1969;204(1):181–93.
Davson H, Segal MB. Physiology of the CSF and bloodbrain barriers. Boca Raton: CRC Press; 1996.
DePasquale M, Patlak CS, Cserr HF. Brain ion and volume regulation during acute hypernatremia in Brattleboro rats. Am J Physiol Renal Physiol. 1989;256(6):F1059–66.
Pollay M, Curl F. Secretion of cerebrospinal fluid by the ventricular ependyma of the rabbit. Am J Physiol Legacy Content. 1967;213(4):1031–8.
Hendry EB. The osmotic pressure and chemical composition of human body fluids. Clin Chem. 1962;8(3):246–65.
Reid EW. Intestinal absorption of solutions. J Physiol. 1902;28(3):241–56.
Diamond JM, Bossert WH. Standinggradient osmotic flow: a mechanism for coupling of water and solute transport in epithelia. J Gen Physiol. 1967;50(8):2061–83.
Curran PF, Macintosh JR. A model system for biological water transport. Nature. 1962;193(4813):347–8.
Friedman MH. Epithelial transport. Principles and models of biological transport. Berlin: Springer; 2008. p. 193–234.
Pollay M. Formation of cerebrospinal fluid: Relation of studies of isolated choroid plexus to the standing gradient hypothesis. J Neurosurg. 1975;42(6):665–73.
Sackin H, Boulpaep EL. Models for coupling of salt and water transport; proximal tubular reabsorption in Necturus kidney. J Gen Physiol. 1975;66(6):671–733.
Schafer JA, Patlak CS, Andreoli TE. Fluid absorption and active and passive ion flows in the rabbit superficial pars recta. J Am Soc Nephrol. 1977;11(4):784–800.
Hill AE. Solutesolvent coupling in epithelia: a critical examination of the standinggradient osmotic flow theory. Proc R Soc Lond B. 1975;190(1098):99–114.
Hill AE. Saltwater coupling in leaky epithelia. J Membr Biol. 1980;56(3):177–82.
Diamond JM. Osmotic water flow in leaky epithelia. J Membr Biol. 1979;51(3):195–216.
Pedley TJ, Fischbarg J. Unstirred layer effects in osmotic water flow across gallbladder epithelium. J Membr Biol. 1980;54(2):89–102.
Wright EM, Smulders AP, Tormey JD. The role of the lateral intercellular spaces and solute polarization effects in the passive flow of water across the rabbit gallbladder. J Membr Biol. 1972;7(1):198–219.
Papp EA, Leergaard TB, Calabrese E, Johnson GA, Bjaalie JG. Waxholm Space atlas of the Sprague Dawley rat brain. Neuroimage. 2014;97:374–86.
Quay WB. Regional and quantitative differences in the postweaning development of choroid plexuses in the rat brain. Brain Res. 1972;36(1):37–45.
Keep RF, Jones HC, Cawkwell RD. A morphometric analysis of the development of the fourth ventricle choroid plexus in the rat. Dev Brain Res. 1986;27(1):77–85.
Neurohr GE, Amon A. Relevance and regulation of cell density. Trends Cell Biol. 2020;30(3):213–25.
McKay MD, Beckman RJ, Conover WJ. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics. 1979;21(2):239–45.
Han ME, Kim HJ, Lee YS, Kim DH, Choi JT, Pan CS, Yoon S, Baek SY, Kim BS, Kim JB, Oh SO. Regulation of cerebrospinal fluid production by caffeine consumption. BMC Neurosci. 2009;10(1):110.
Karimy JK, Kahle KT, Kurland DB, Yu E, Gerzanich V, Simard JM. A novel method to study cerebrospinal fluid dynamics in rats. J Neurosci Methods. 2015;241:78–84.
Bothwell SW, Omileke D, Patabendige A, Spratt NJ. CSF secretion is not altered by NKCC1 nor TRPV4 antagonism in healthy rats. Brain Sci. 2021;11(9):1117.
Murtha LA, Yang Q, Parsons MW, Levi CR, Beard DJ, Spratt NJ, McLeod DD. Cerebrospinal fluid is drained primarily via the spinal canal and olfactory route in young and aged spontaneously hypertensive rats. Fluids Barriers CNS. 2014;11(1):12.
Rosenberg GA, Kyner WT, Estrada E. Bulk flow of brain interstitial fluid under normal and hyperosmolar conditions. Am J Physiol Renal Physiol. 1980;238(1):42–9.
Goodman JA, Kroenke CD, Bretthorst GL, Ackerman JJH, Neil JJ. Sodium ion apparent diffusion coefficient in living rat brain. Magn Reson Med. 2005;53(5):1040–5.
Bandaranayake NM, Nemoto EM, Stezoski SW. Rat brain osmolality during barbiturate anesthesia and global brain ischemia. Stroke. 1978;9(3):249–54.
Lui ACP, Polis TZ, Cicutti NJ. Densities of cerebrospinal fluid and spinal anaesthetic solutions in surgical patients at body temperature. Can J Anaesth. 1998;45(4):297.
Acknowledgements
We thank Kjeld Møllgård, Department of Cellular and Molecular Medicine, University of Copenhagen for providing the brain slices used for immunohistochemistry, Anca Stoica and Trine ToftBertelsen, Department of Neuroscience, University of Copenhagen for the immunostaining and protocol thereof. We would also like to thank Thomas Zeuthen, University of Copenhagen, for feedback on the model development and are grateful to both him and Victoria Makrides, University of Zurich, for suggestions on the manuscript.
Funding
This study was financed, in part, by the Swiss National Science Foundation through project 205321_182683 (to VK), by the Fidelity Bermuda Foundation, and by the Novo Nordisk Foundation, Tandem program (NNF17OC0024718, to NM).
Author information
Authors and Affiliations
Contributions
PRK: software, methodology, investigation, data curation, formal analysis, visualization, validation, writing—original draft. VC: validation, writing—original draft. EKH: data curation, writing—review and editing. NMA: funding acquisition, conceptualization, data curation, writing—review and editing. VK: funding acquisition, conceptualization, project administration, supervision, formal analysis, writing—review and editing. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All animal experiments conformed to the legislations for animal protection and care in the European Community Council Directive (2010/63/EU) and followed all ethical regulations under the animal permission no. 201615020100944.
Consent for publication
Not applicable.
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.
Supplementary Information
Additional file 1.
Derivation of the functional unit tip boundary condition, estimation of choroid plexus surface area and luminal membrane permeability, Péclet number calculation, and estimation of the protected length
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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Razzaghi Khamesi, P., Charitatos, V., Heerfordt, E.K. et al. Are standing osmotic gradients the main driver of cerebrospinal fluid production? A computational analysis. Fluids Barriers CNS 20, 18 (2023). https://doi.org/10.1186/s12987023004192
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12987023004192
Keywords
 CSF production
 Choroid plexus
 Computational modeling