About the Authors:
Philip Pearce
* E-mail: [email protected]
Affiliation: School of Mathematics, University of Manchester, Manchester, M13 9PL, United Kingdom
ORCID http://orcid.org/0000-0001-5788-3826
Paul Brownbill
Affiliations Maternal and Fetal Health Research Centre, Institute of Human Development, University of Manchester, Manchester, United Kingdom, Maternal and Fetal Health Research Centre, St. Mary’s Hospital, Central Manchester University Hospitals NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, M13 9WL, United Kingdom
Jiří Janáček
Affiliation: Department of Biomathematics, Institute of Physiology, v.v.i., Academy of Sciences of the Czech Republic, Prague, Czech Republic
Marie Jirkovská
Affiliation: Institute of Histology and Embryology, First Faculty of Medicine, Charles University, Albertov 4, CZ-12801 Prague 2, Czech Republic
Lucie Kubínová
Affiliation: Department of Biomathematics, Institute of Physiology, v.v.i., Academy of Sciences of the Czech Republic, Prague, Czech Republic
Igor L. Chernyavsky
Affiliation: School of Mathematics, University of Manchester, Manchester, M13 9PL, United Kingdom
Oliver E. Jensen
Affiliation: School of Mathematics, University of Manchester, Manchester, M13 9PL, United Kingdom
Introduction
Oxygen is important for both fetal and placental development throughout pregnancy [1]. Much of the research relating the structure of the placenta to fetal oxygen uptake has taken a stereological approach, using two-dimensional histological slides to calculate structural quantities such as average volumes, surface areas and lengths [2]. Recently, three-dimensional (3D) representations of the feto-placental vasculature, obtained by reconstructing confocal laser scanning microscope images, have become available [3, 4], enabling more detailed analysis of structures in specific regions of the placenta. In this paper, we use mathematical modeling in combination with 3D images to investigate blood flow and oxygen transfer in feto-placental capillaries.
The feto-placental vasculature is contained inside villous trees, which bathe in maternal blood, allowing the transfer of nutrients without mixing of the two blood supplies. As fetal blood flows from the umbilical artery to the umbilical vein, it passes through a series of vessels of varying spatial scale. At the smallest scale are the capillaries in the terminal villi, bulbous structures that are thought to be the main sites for the transfer of oxygen from mother to fetus in the final trimester of pregnancy [5]. As the primary functional unit of the mature placenta, terminal villi and the fetal capillaries therein merit careful investigation.
A specific characteristic of fetal capillaries in the terminal villi that we focus on in this study is the presence of localized dilations. Such dilations, often called sinusoids, have been observed in numerous studies [5]. The dilations are often seen in conjunction with a vasculo-syncytial membrane, where the capillary bulges against the villous surface. Due to organelle displacement, in these regions the maternal and fetal blood are separated only by a simple bilayer consisting of a capillary endothelium and a thin syncytiotrophoblast. Diverse physical explanations for the localized dilations have been hypothesized, including a reduction in feto-maternal diffusion distance, a reduction in the total resistance of the feto-placental vasculature and a deceleration in local blood flow to facilitate exchange [5].
There are a few sources linking vasculo-syncytial membranes, and the corresponding localized fetal capillary dilations, to oxygen transfer. First, there is thought to be an increased amount of fetal capillary dilation in pregnancies at high altitude, where conditions are hypoxic [6], evidenced by a reduced mean trophoblast thickness [7] and an increased percentage of villous volume occupied by fetal capillaries [8]. Second, a lack of terminal villi and vasculo–synyctial membranes, amongst other markers, characterizes a pathology known as delayed villous maturation, also known as defective villous maturation or distal villous immaturity [9–14]. There is an increased prevalence of delayed villous maturation in stillbirths where the autopsy findings were consistent with fetal hypoxia [10, 15]; it is also associated with poor neurologic outcomes due to hypoxia [9, 13]. There is thought to be an approximately 5% risk of delayed villous maturation recurring in subsequent pregnancies [15].
Past efforts to model the placenta theoretically can generally be split into those that model the maternal side with assumptions about the fetal side, and vice versa; extensive reviews of the historical literature are available in [16, 17]. In recent years, maternal flow has been modeled using Darcy’s law for flow in a porous medium [16], as a uniform flow field [18, 19], or using the Navier–Stokes equations in combination with images of fetal villous trees [20]. Where applicable, in these studies oxygen transfer to the fetal vasculature has been simplified, either by including a reaction term in the advection–diffusion equation [16], or by assuming fetal vessels to be oxygen sinks [18, 19, 21].
Oxygen diffusion from maternal to fetal blood within villous branches has been simulated in realistic 2D geometries from histology slides [22] and 3D geometries reconstructed from confocal laser scanning microscopy [23], with fetal blood vessels assumed to be perfect oxygen sinks in both cases. This assumption corresponds to the situation where all of the provided oxygen is carried away by the fetal blood, and the oxygen transfer rate in this case is expected to provide an upper bound on oxygen transfer to the fetal capillary when blood flow is taken into account.
Attempts to model blood flow in the feto-placental vasculature have simplified the microstructure in the terminal villi [24]; a problem is how to capture the microstructure in the smallest feto-placental capillaries accurately. One solution is to simulate the shape and structure of the placenta and the fetal vessels it contains [24–26]. However, to improve on existing models it is vital to understand characteristics such as the vascular resistance of real capillaries in the feto-placental vasculature. This is particularly important when modeling the transfer of a substance such as oxygen, which can be dominated by advection (due to oxygen’s affinity for hemoglobin in red blood cells) and is therefore significantly influenced by patterns of blood flow through capillaries. Spatial characteristics of oxygen transfer are also expected to be affected by the extent to which the oxygen content of the fetal blood has equilibrated with that of the maternal blood.
A multi-scale mathematical model of the entire placenta, spanning spatial scales from the umbilical cord to the terminal villi, would bring significant benefits. First, it would help identify the physical reasons for observed structures in healthy and diseased placentas. Second, it could aid clinicians in early detection and prevention of pathologies. Finally, it could enable patient-specific diagnostics using future imaging technologies, improving personalized medicine during pregnancy. This is an ambitious goal that will demand integration of contributions at many levels of spatial organization.
Three-dimensional images of fetal capillaries can be a useful tool for building theoretical models that reflect the microstructure of the placenta accurately. Because it is currently not feasible to solve the equations for blood flow and oxygen transfer in realistic geometries throughout an entire placenta, the equations can be solved in a selection of geometries to find statistical distributions of important physical properties. Such simulations are performed in the present study, using 3D images of fetal villous and capillary surfaces, albeit for a small number of examples. Several key properties that will be useful in future multi-scale models are calculated, including the vascular resistance and the total oxygen transfer rate.
It is important to be able to quantify the effects of statistical variability, experimental uncertainty and pathological placental structure on the generated statistical distribution of physical properties. This paper demonstrates two separate methods that can be used to quantify such effects. In the first method, scaling arguments are used to make predictions and assess the sensitivity of the model to changes in the geometry and the parameters. Three separate physical regimes are identified, which capture the oxygen transfer across a range of Péclet numbers: a diffusion-limited regime where oxygen transfer is determined by diffusion through the surrounding villous volume, and two flow-limited regimes, where oxygen transfer is determined by the pressure drop through the capillary and the vascular resistance. These physical regimes, and the associated scalings for the oxygen transfer, are summarized in Table 1. In the second method, an idealized axisymmetric model is used to quantify the influence of localized feto-placental capillary dilations, or a lack thereof, on oxygen transfer. As well as furthering our understanding of the effect of pathologies on blood flow and oxygen transfer, the results, which are validated by comparison with the 3D image-based modeling, provide physiological insight to help to prevent, identify and treat delayed villous maturation.
[Figure omitted. See PDF.]
Table 1. Summary of oxygen transfer scaling results for Peeff = Bu0R0/D ≫ 1.
https://doi.org/10.1371/journal.pone.0165369.t001
Methods
Image acquisition, mesh conversion and skeletonization
Three-dimensional (3D) images of fetal capillaries and villous surfaces were used from previous studies, in which the methods for specimen preparation, image acquisition and 3D image reconstruction are described in detail [3, 4]. All tissue was acquired with written consent from fully-informed women attending the General University Hospital in Prague for their antenatal care. Tissues were collected under ethical approval by the Ethics Committee of the General University Hospital in Prague. Samples from healthy placentas were taken following uniform random sampling no later than 10 min after delivery. The samples were fixed in 4% formaldehyde containing 0.5% eosin for at least 24 h, and embedded in paraffin wax; thick sections (120 μm) were cut. Stacks of optical sections 1 μm apart were imaged using a Biorad MRC 600 laser scanning confocal microscope. Surfaces of villi and capillaries were created and rendered using the image analysis program Ellipse and converted into a 3D binary image. The image was resampled as a grey image with lower resolution and treated by 3D Gaussian filtration. The triangulated isosurface of the grey value was detected.
Triangulated isosurfaces were reconstructed using a Poisson surface reconstruction algorithm in the open-source mesh processing package Meshlab (version 1.3.3, http://meshlab.sourceforge.net) to form a smoother surface. The surfaces were imported into Comsol Multiphysics (version 5.2, http://comsol.com) and converted into 3D tetrahedral meshes using the inbuilt meshing algorithm. In the process, the surface mesh was further simplified to remove small defects. Geometries were partitioned with a plane to provide a flat surface for inflow and outflow boundary conditions to be applied (Fig 1A).
[Figure omitted. See PDF.]
Fig 1. Methods.
(A) An example of how 3D images are partitioned with a plane (shown in blue) to provide a surface for boundary conditions to be applied. The red surface represents the capillary surface (endothelium, Γcap) and the grey surface represents the villous surface Γvil. In this example the capillary bifurcates within the villous branch. (B) Example of the skeletonization of a capillary from a 3D image. The skeletonization line, representing the centreline of the lumen, is shown in black. (C) The 3D images of fetal capillaries and villous surfaces used in the simulations (see S1 Images). The images shown have been partitioned by a plane on which boundary conditions are applied, for inflow (Γin), outflow (Γout) and villous tissue (Γ0). (D) Polynomial fetal oxygen–hemoglobin dissociation law from [27] (solid line) and the linear approximation between 0 and 60 mmHg, passing through the origin, found using a least-squares fit (dashed line). The gradient of the linear approximation is K = 0.019 mmHg−1. (E) Schematic diagram of the idealized axisymmetric model of a fetal capillary dilation, showing an axisymmetric tube with a localized dilation. Blood flows into the capillary through γin and leaves through γout. Oxygen is provided along the dilated section of the capillary γd, denoted by the thicker black line. The undilated sections of the capillary are labeled γu.
https://doi.org/10.1371/journal.pone.0165369.g001
For skeletonization, triangulated isosurfaces were converted into tetrahedral meshes in the open-source mesh generating package Gmsh (version 2.12, [28]) and imported into the image analysis software Avizo (Lite version 9.0, https://www.fei.com/software/avizo3d). The capillaries were skeletonized using the Avizo Skeletonization Pack (Fig 1B). This uses a distance-ordered thinning algorithm to remove voxels successively until a skeleton with the thickness of a single voxel remains; the algorithm is described in detail in [29]. The total length of each capillary was defined as the total length of the skeletonization line, including all branches. At each point along the skeletonization line, the radius of the capillary was defined as the distance between the voxel on the line and the nearest capillary surface voxel. At each point on the capillary surface, the villous distance was defined as the distance between the capillary surface voxel and the nearest villous surface voxel. The unpartitioned geometry was used for the calculation of the mean and standard deviation of the villous distance. Further details of the methods used for image analysis are provided in S1 Appendix.
Simulation of blood flow and oxygen transfer in 3D geometries
The parameters used in the simulations are summarized in Table 2. Blood flow in feto-placental capillaries has low Reynolds number, and was therefore modeled using the Stokes equations(1)where u is the flow velocity field, p is the pressure and μ is the dynamic viscosity. For simplicity, flow was assumed to be steady and non-Newtonian effects of red blood cells (RBCs) on the flow were ignored; the suspension was assumed to have constant effective viscosity μ. At small scales, the viscosity of blood is known to depend on the radius of the capillary it is flowing through, via the Fahraeus–Lindqvist effect [30]. The dynamics of blood flow in a capillary of varying radius, however, are not clear. Incorporating the effect of RBCs in 3D geometries like the ones studied in this paper would necessitate computations involving discrete RBCs, which is beyond the scope of the study.
[Figure omitted. See PDF.]
Table 2. Literature-based and calculated parameters used in the formulation of the problem and the estimation of oxygen advection enhancement by fetal hemoglobin.
https://doi.org/10.1371/journal.pone.0165369.t002
Oxygen transfer in villous branches was modeled by the advection–diffusion equation, with oxygen transfer outside capillaries assumed to be by diffusion alone. The effect of oxygen binding to hemoglobin in blood was taken into account using the method described in [18], which is summarized briefly here.
Blood is assumed to be a continuum, with hemoglobin and plasma occupying the whole volume, a reasonable assumption at normal hematocrits [33]. The total oxygen concentration is split into the concentration of oxygen dissolved in plasma c (per unit volume of blood) and the concentration of oxygen bound to hemoglobin ch (per unit volume of blood). Oxygen is assumed to diffuse only when dissolved in the plasma due to the low diffusion coefficient of RBCs [34], to give(2)where D is the diffusion coefficient of oxygen in the plasma (assumed to be constant throughout the villous branch). Oxygen is assumed to bind to hemoglobin instantaneously, which can be justified by comparing the typical time for oxygen–hemoglobin dissociation (τd ≈ 20 ms for the slowest process [35]) to the typical transit time of blood through a capillary (τt ≈ 0.3 s for a capillary of length 100μm, with centreline blood flow velocity of 300 μm/s [36]). The concentration of oxygen bound to hemoglobin is related to the partial pressure of oxygen in the plasma PO2 by the oxygen–hemoglobin dissociation law(3)Here cmax is the oxygen content of fetal blood at full saturation and the partial pressure of oxygen in the plasma has been written using Henry’s law, where khn is the Henry’s law coefficient and ρbl is the density of blood. Inserting Eq (3) into Eq (2) gives(4)Finally, the fetal oxygen–hemoglobin dissociation law given in [27] is linearized in the region between 0 mmHg and 60 mmHg (where 60 mmHg is the highest partial pressure of oxygen the fetal blood is expected to be exposed to [37]) using a least-squares fit, with the constraint that the curve passes through the origin (Fig 1D), so that(5)for K = 0.019 mmHg−1. Inserting Eq (5) into Eq (4) gives(6)where is a parameter quantifying the enhancement to advection due to oxygen’s affinity for hemoglobin.
Using typical values of the parameters, it was estimated that B ≈ 141 (see S2 Appendix). The uptake of oxygen by the placenta itself, via the outer layer of the villous surface (the syncytiotrophoblast), was neglected. The placenta is thought to consume approximately 40% of the oxygen leaving the maternal circulation [38]. However, it is not known what proportion of this consumption takes place in the terminal villi.
Surfaces where boundary conditions were applied are labeled in Fig 1C. The pressure difference between the inflow boundary Γin and the outflow boundary Γout was fixed (for flow in a rigid domain, pressure differences, rather than absolute pressures, are relevant). As is common in flow simulations [39], the normal stress was specified at these boundaries, with the tangential velocity vector ut set to zero. Thus the inlet pressure was fixed at p = ΔP, say, by specifying the normal stress n ⋅ σ ⋅ n = −ΔP, where σ = −pI + μ(∇u + (∇uT)) is the stress tensor, I is the identity tensor and n is the outer unit normal to the boundary. A no-slip condition was applied on the capillary wall Γcap, neglecting the effect of the endothelial glycocalyx. In summary, the boundary conditions applied to the velocity field were(7a) (7b) (7c)At the inflow boundary Γin, the blood was considered to be deoxygenated. At the villous surface Γvil, the concentration of oxygen was assumed to be equal to the concentration cmat of oxygen in the maternal blood in the intervillous space. There is evidence that this is reasonable from previous studies showing the maternal flow to have spatial gradients over the lengthscale of an entire placentone [16, 20]. To ensure the simulations were performed in the range of partial pressures of oxygen where the linearization of the oxygen–hemoglobin dissociation law is valid (Fig 1D), a value of cmat = 0.07 mol/m3 was taken, corresponding to a partial pressure of oxygen in the maternal blood of 52.5 mmHg. For larger partial pressures, the effect of saturation of hemoglobin may need to be taken into account. Zero diffusive flux conditions were applied at the outflow boundary Γout and the surface Γ0. In summary, the boundary conditions applied to the concentration were(7d)where n is the coordinate normal to the boundary.
Branched capillaries (Image 1 and Image 3, Fig 1C) were assumed to have one inflow and two outflow boundaries, with blood assumed to flow in the direction where the branching angle was smallest. In the unbranched capillary (Image 2, Fig 1C), Γin and Γout were designated arbitrarily; it was found in computations that reversing Γin and Γout led to less than 0.01% difference in the flow rate (to be expected from the reversibility of Stokes flow) and less than 0.1% difference in the total oxygen transfer rate.
The equations were solved numerically in Comsol Multiphysics using the direct solver MUMPS [40, 41]. Solutions were computed on dual E5-2643v3 CPU nodes with hyperthreading disabled and 768GB of RAM. Results were checked to be independent of the mesh, which consisted typically of ∼106–107 tetrahedral elements. Simulations were performed on three different 3D geometries (Fig 1C). Geometric characteristics of the corresponding capillary and villous surfaces are summarized in Table 3. The total oxygen transfer rate N from the maternal blood was calculated by integrating the advective flux leaving the domain through Γout, which was checked to be equal to the diffusive flux entering the domain across Γvil minus the diffusive flux leaving through Γin (an artifact due to the imposed boundary condition, which does not affect the results significantly for flows in which advection dominates diffusion). The integral of the diffusive flux through each boundary was computed accurately by a summation of the corresponding components of the boundary residual vector over every node on the surface [42].
[Figure omitted. See PDF.]
Table 3. Geometric and calculated properties of the three 3D images used in the paper.
https://doi.org/10.1371/journal.pone.0165369.t003
Due to the linearity of the Stokes equations, the flow field needs only one computation for each capillary to calculate the vascular resistance , where ΔP is the pressure drop between inflow and outflow and Q is the flow rate. Similarly, a single solution of the advection–diffusion Eq (6) can be used to calculate the non-dimensional concentration field, , where cin is the inflow oxygen concentration. Rearranging this equation gives the dimensional concentration field for any values of cin and cmat. Further details of the methods use in the 3D simulations are given in S1 Appendix.
Idealized axisymmetric model and optimization of capillary dilation
The parameters used in the axisymmetric model of an individual capillary dilation are summarized in Table 2. Blood flow and oxygen transfer were modeled by Eqs (1) and (6) in a simplified geometry consisting of an axisymmetric capillary of length L, with a dilated section of length λ in the centre of the capillary wall (Fig 1E), where the undilated radius was R0. To model a vasculo-syncytial membrane, oxygen was assumed to be provided only along the dilated section of the capillary, due to its proximity to the maternal blood; deoxygenated blood was assumed to flow into the capillary, driven by a pressure difference across the capillary ΔP, which was fixed to give a centreline velocity of 300 μm/s in an undilated capillary [36]. A no-slip condition was applied on the capillary wall. The diffusive flux of oxygen was set to zero on the undilated walls and the outflow boundary. To summarize, the applied boundary conditions were(8a) (8b) (8c) (8d)where the boundaries γi are illustrated in Fig 1E. Arteriolar and venular resistances and oxygen transfer upstream and downstream of the capillary were ignored for simplicity (although a modification of condition Eq (8b) on γu is discussed in S3 Appendix).
To reduce the degrees of freedom (DOF) in the optimization, the shape of the dilation R(z) was expanded as a Fourier series truncated to n terms,(9)where z is the axial coordinate. To ensure a continuous and smooth capillary wall, the conditions(10)were imposed. The conditions Eq (10) provided constraints on the parameters ak and bk. Further, to avoid errors in the simulations, the extrema of R(z) were constrained between prescribed values; it was checked that the chosen values did not affect the optimal solutions. In cases where the shape was assumed to be symmetric about the axial midpoint of the dilation, ak and bk were set to zero for odd values of k < n.
Shape optimization was performed using the fmincon function in Matlab, which uses a Sequential Quadratic Programming algorithm [43]. No further regularization was imposed beyond the implicit regularization in the choice of n (and in some cases the assumption of symmetry), which restricted the number of DOF in the shape of the dilation. At each step in the optimization process, the finite-element problem was solved numerically in Comsol Multiphysics, which output the total oxygen transfer rate N. The number of DOF in the optimization was increased successively.
The function fmincon found a local maximum in N in the parameter-space of a0, ak and bk, k < n. An important aspect of the optimization was how to check that the local maximum in N was in fact a global maximum. In cases with one DOF, the solution was checked to be optimal by comparing against a curve of N versus the maximum radius. In cases with two DOF, results were checked against contour plots of N in parameter space. In cases with successively increasing DOF, at least ten optimizations were performed, with initial conditions taken as a random sample in parameter space, constrained to ensure the dilation did not occlude the channel. If more than one local maximum was found, the shape with the largest oxygen transfer was taken to be the optimal solution. The non-invasive optimization method described here, using Comsol Multiphysics as a “black box” solver, was chosen for speed and simplicity.
When analyzing a dilation from a skeletonized capillary branch of length L, quantifying the geometry of a dilation required measuring the maximum radius Rmax and the undilated radius R0, and defining the start and end points of the dilation to calculate its length λ. Two different definitions were used: the first defined the beginning and end of the dilation as local minima in the radius, with the undilated radius being the minimum radius of the capillary branch; the second used the average radius of the capillary branch as the undilated radius and defined the edges of the dilation where the radius passed its average radius. These two separate definitions provided a range of geometrical values to be used as input to the axisymmetric model.
Results
Blood flow and oxygen transfer in real geometries: scaling and simulations
The calculated values of the vascular resistance R, the vascular resistance per unit capillary volume Rcap and the vascular resistance per unit villous volume Rvil for each capillary in Fig 1C are given in Table 3. The capillaries with the least vascular resistance are those with branches, as expected; the total vascular resistance of two capillaries in parallel, each with the same resistance, is half the resistance of each capillary. As well as branching, the most important geometrical determinants of the vascular resistance are the average radius R0 and the total length L of the capillaries; recall from Poiseuille’s law that in a straight capillary with fully developed flow, R is given by(11)
The vascular resistance of the capillaries can be used to predict which capillary is expected to receive the most oxygen from the maternal side. To this end, consider a straight capillary of radius R0 and length L, with fixed concentration cmat at the wall and fully developed Poiseuille flow entering deoxygenated. The scalings determining the oxygen transfer rate depend on the magnitude of the effective Péclet number , which gives the ratio of transport by advection to transport by diffusion, including the effect of oxygen transport by hemoglobin in red blood cells. Here D is a diffusion coefficient, u0 is a characteristic flow speed, which is proportional to the pressure drop ΔP across the capillary, and B is a numerical factor that quantifies the enhancement to advection as a result of oxygen binding to hemoglobin (see Eq (6)). When Peeff ≫ 1, advection dominates and, close to the capillary entrance, oxygen transfer is confined to a thin diffusive boundary layer near the capillary wall, where the axial component of the fluid velocity increases linearly with the distance from the wall and the curvature of the wall can be neglected [44]. Axial diffusion is expected to be of relative magnitude , and can therefore be neglected everywhere except in the region very close to the entrance [45], which we ignore for simplicity. The resulting boundary-layer equation can be solved with a similarity variable to find that the local oxygen flux per unit area along the capillary wall is proportional to , where z is the axial coordinate [44]. Integrating across the capillary surface and using Eq (11) gives the total oxygen transfer rate (labeled Nflow as it corresponds to the flow-limited regime) to be(12)where [44]; here Γ is the gamma function. The approximation Eq (12) is expected to hold for sufficiently short capillaries with Peeff ≫ 1, as long as Peeff is not large enough for oxygen transfer to be limited by diffusion through the villous volume. However, for the images in this study, if the ratio of capillary length L to radius R0 plus villous thickness d is sufficiently large (), oxygen is expected to diffuse to the centreline of the downstream region of the capillary, suppressing concentration gradients across the capillary cross-section [46] (if the simplifying assumption is made that capillary branches do not affect each other, the longest path between the inflow boundary and outflow boundary through a branched capillary could be used as L in this inequality). In this case the concentration is said to be fully developed, or equilibrated, in the downstream region of the capillary. The amount of oxygen carried away by the flow is then proportional to the amount of blood flowing through the capillary,(13)
Using the scalings Eqs (12) and (13), ignoring the effects of capillary branching and the space between capillary and villous surfaces for simplicity, it would be expected that for larger Péclet numbers, the greatest oxygen transfer would occur in Image 2 and the least in Image 3 (comparing L2/R; see Table 3). For Péclet numbers small enough for the concentration to become equilibrated downstream, the greatest oxygen transfer would be expected to occur in Image 3 and the least in Image 1 (comparing 1/R).
The predictions made above are confirmed by the numerical results (Fig 2A). At larger values of ΔP, corresponding to Peeff ∼ 103, oxygen transfer is found to be highest in Image 2 and lowest in Image 3, as expected. At smaller values of ΔP, corresponding to Peeff ∼ 102, oxygen transfer is found to be highest in Image 3 and lowest in Image 1, in line with the predicted behaviour. In this case, the oxygen distribution is insensitive to the flow of blood leaving the capillary (Fig 2C), and most of the oxygen transfer takes place at the capillary entrance and across the nearby villous surface (Fig 2B). For a higher pressure drop, oxygen transfer takes place across the entire capillary surface, although the largest amount of oxygen still enters around the capillary entrance and the nearby villous surface (Fig 2B and 2C). The majority of oxygen enters the capillary through the sections of the capillary wall closest to the villous surface (Fig 2B), illustrating the “diffusional screening” effect identified in [22, 47]. In fact, around the entrance to the villous branch, some oxygen leaves through the walls of the outflow capillaries furthest from the villous surface.
[Figure omitted. See PDF.]
Fig 2. Results of numerical simulations on 3D geometries.
(A) Oxygen flux entering each capillary versus pressure drop ΔP. The range of values of ΔP is broadly physiological, leading to flow velocities of around 300 μm/s (Table 3), which is the approximate flow velocity in normal capillaries [36]. Peeff is in the range 102–103 as indicated. (B) Distribution of oxygen flux entering through the capillary (top) and villous (bottom) surfaces, from simulations on Image 1. Results are shown for ΔP = 2 Pa (left) and ΔP = 20 Pa (right). (C) Plots of the oxygen concentration field on three planes, 15 μm apart, from simulations on Image 1. The planes are shown on the villous surface on the left. Results are given for ΔP = 2 Pa (top) and ΔP = 20 Pa (bottom). Capillary walls are denoted by black lines. The inflow and outflow capillaries are labeled on the plane at the villous entrance (left). (D) Wall shear stress at the capillary surface, from a simulation on Image 1 (ΔP = 2 Pa).
https://doi.org/10.1371/journal.pone.0165369.g002
Another way to estimate which capillary is expected to receive the most oxygen for higher values of ΔP is by solving the diffusion equation in the villous volume surrounding the capillary, with the capillary surface assumed to be deoxygenated by setting c = 0 on Γcap, which is less computationally intensive than solving for flow in the capillary and avoids the need to resolve the boundary layers in the concentration that arise at larger values of ΔP. The calculation, which is similar to those performed in [23], corresponds to the situation where all of the oxygen is carried away by the flow and oxygen transfer is limited by diffusion through the villous volume, providing an upper bound, Nmax, on the total oxygen transfer rate N. The calculated numerical values of Nmax for each capillary are given in Table 3. The values of the total oxygen transfer per unit area A, given by Nmax/A, are found to be of the same order of magnitude as the results of [23]. The value of Nmax can be related to the geometrical properties of each villous branch by considering a cylindrical capillary with a deoxygenated surface, with radius R0 and surface area A, inside a cylindrical villous branch with an oxygenated surface, with radius R0 + d. Solving the diffusion equation in the annulus between the concentric cylinders [44], the predicted maximum oxygen transfer rate is given by(14)Taking the radius of each capillary to be (as for a cylindrical capillary of volume V and length L), and using values of V, L and d from Table 3, Eq (14) gives the theoretical predictions of Nmax given in Table 3, which are of the same order of magnitude as the numerical results, and predict that Image 2 has the highest maximum oxygen transfer rate and Image 3 has the smallest, in agreement with the numerical results.
An oxygen transfer ratio Nratio = N/Nflow for each villous branch can be calculated by comparing oxygen transfer N in the villous branch to oxygen transfer Nflow in the situation where oxygen is provided directly at the surface of the capillary by imposing c = cmat on Γcap (Fig 1C). There are expected to be two main factors influencing Nratio. First, a capillary that is closer to the villous surface, on average, is expected to have a higher Nratio than a capillary that is further away, because the villous volume will have less effect on oxygen transfer in this case. Second, a capillary in which the concentration becomes more equilibrated is expected to have a higher Nratio, because in this case the oxygen concentration at the capillary surface in the downstream region of the capillary is closer to that of the villous surface. Therefore Nratio should be higher at lower values of ΔP. Image 2 is expected to have the highest Nratio because it has the highest value for the longest path length through its capillary, which is also the closest to the villous surface on average (Table 3). Images 1 and 3 are expected to have similar values of Nratio because the capillary in Image 1 is longer than the capillary in Image 3, but further away from the villous surface on average. As expected, Image 2 is found to have the highest Nratio at both ΔP = 0.5 Pa and ΔP = 20 Pa, while Images 1 and 3 are found to have roughly the same Nratio at both ΔP = 0.5 Pa and ΔP = 20 Pa (Table 3).
The average magnitude of the calculated capillary wall shear stress is found to vary between approximately 0.02 Pa and 0.25 Pa (or 0.2 dyn/cm2 and 2.5 dyn/cm2), which is close to what has recently been estimated to be exerted on the walls of the fetal villous trees by the maternal blood [20]. However, large spikes in wall shear stress are found around the inside of capillary bends, with dips around the outside of bends (Fig 2D). This could be considered counter-intuitive, as localized capillary dilations are thought to be important for the transfer of oxygen from the maternal blood, which would be enhanced by an increase in wall shear stress. A more thorough discussion of localized dilations can be found below.
It is important to understand the sensitivity of the model to changes in the geometry and the parameters. The scalings Eqs (12) and (13) suggest that, for a capillary of length L and radius R0, a uniform expansion (or contraction) of all lengthscales by magnitude k would lead to an expected increase (or decrease) in the total oxygen transfer rate N of k3 for a smaller Péclet number (where the concentration becomes equilibrated in the capillary), or k5/3 for a larger Péclet number (where the boundary-layer equation is valid, using the fact that the vascular resistance R is proportional to from Eq (11)). A uniform expansion of 10% in all lengthscales (k = 1.1), which would be expected to cause an increase in the flow rate Q of 33%, would therefore be expected to cause an increase in N of around 33% for a smaller Péclet number or 17% for a larger Péclet number. Indeed, our computations show that applying a 10% uniform expansion to Image 2 leads to the expected 33% increase in Q, with a 33% increase in oxygen transfer at ΔP = 0.5 Pa and a 15% increase at ΔP = 20 Pa. A uniform expansion (or contraction) by magnitude k in the capillary radius R0 alone would be expected to have a slightly different effect, causing an increase (or decrease) in N of k4 for a smaller Péclet number or k4/3 for a larger Péclet number. The sensitivity of the model to the parameters can also be assessed using the scalings Eqs (12) and (13) with the Formula (11) for the vascular resistance R.
Regression equation for the total oxygen transfer rate
A general regression equation relating N to ΔP across parameter space where Peeff ≫ 1 should take into account the scaling Eq (13) at lower values of ΔP, the scaling Eq (12) at higher values of ΔP and the upper bound on the oxygen transfer rate Nmax, which are summarized in Table 1. A regression equation satisfying the scalings Eqs (12) and (13) in the case where oxygen is provided directly at the capillary surface is(15)which is similar to an equation regularly used to represent the results of experiments on heat transfer in tubes [48]. Here K1 = cmat B/R and K2 = αcmat(D2 BL2/R)1/3. A correction to Eq (15) to account for the effect of the surrounding villous volume can be calculated by considering the diffusion of oxygen between a cylindrical villous surface of radius R0 + d surrounding a cylindrical capillary of radius R0, with both cylinders of equal length. The oxygen concentration is fixed on the surface of the outer cylinder, so that(16)The oxygen transfer rate per unit area is fixed on the inner cylinder, so that(17)where N is the total oxygen transfer rate and A is the surface area of the inner cylinder. Solving the axisymmetric diffusion equation in the space between the cylinders, subject to the boundary conditions Eqs (16) and (17), assuming axial diffusion to be negligible, gives the concentration cinner at the surface of the inner cylinder to be(18)where Nmax is the theoretical maximum for the oxygen transfer rate given by Eq (14). To increase the accuracy of the regression equation, the numerical values of Nmax given in Table 3 are used here. Now considering oxygen transfer into the capillary, taking the oxygen concentration at the capillary surface to be given by cinner from Eq (18), the oxygen transfer rate Nflow from Eq (15) is multiplied by due to the linearity of the advection-diffusion equation. Using this correction with Eq (15) and rearranging gives a regression equation for the oxygen transfer rate N accounting for the surrounding villous volume,(19)where Nflow is given by Eq (15) and Nmax is the upper bound on the oxygen transfer rate, calculated numerically. It is simple to adapt the regression Eq (19) to plot N versus another quantity such as the flow rate Q or the effective Péclet number Peeff.
The accuracy of the scalings Eqs (12) and (13) and the regression Eq (19) are tested across a range of Péclet numbers in Fig 3. For smaller values of ΔP, where the scaling Eq (13) holds, the data of Fig 2A for each capillary collapse close to a single straight line when NR is plotted against ΔP on log-log axes, as predicted (Fig 3A). For larger values of ΔP, when is plotted against ΔP on log-log axes, the data collapse onto straight lines of approximate gradient 1/3, confirming that the scaling Eq (12) holds for the values of ΔP considered here. The regression curves given by Eq (19) qualitatively capture the numerical results across a range of values of ΔP (Fig 3C).
[Figure omitted. See PDF.]
Fig 3. Scaling and regression.
Throughout the figure, results for Image 1, Image 2 and Image 3 are denoted by blue dots, green triangles and red squares, respectively. (A,B) Scaled oxygen transfer rate (A) NR and (B) N/(L2/R)1/3 versus pressure drop ΔP, on log-log axes. Solid lines denote numerical results and black triangles denote the gradients predicted by (A) scaling Eq (13) and (B) scaling Eq (12). (C) Oxygen transfer rate N versus ΔP with the linearized oxygen–hemoglobin dissociation law. Solid lines denote numerical results and dashed lines denote results generated by the regression Eq (19). The upper bounds Nmax on N are indicated, which have been calculated by solving the diffusion equation in the surrounding villous volume with the capillary assumed to be deoxygenated. (D) Oxygen transfer rate N with the nonlinear oxygen–hemoglobin dissociation law (solid lines) and average advection enhancement in each capillary (dashed lines) versus ΔP. Also shown on the figure is the oxygen advection enhancement parameter B (black dashed line), derived by linearizing the oxygen–hemoglobin dissociation law.
https://doi.org/10.1371/journal.pone.0165369.g003
The numerical results obtained when the nonlinear oxygen–hemoglobin dissociation law is taken into account via replacing Eq (6) with Eq (4) are given in Fig 3D. The close agreement between the results in the nonlinear case and the linearized case can be seen by comparing Fig 3C and 3D. Also plotted in Fig 3D is the average enhancement to advection in the nonlinear case, which depends on the pressure difference ΔP, compared to the linearized value B = 141. The figure shows that the linearization overestimates the average advection enhancement at smaller values of ΔP and slightly underestimates the advection enhancement at larger values of ΔP.
Predictions of the idealized axisymmetric model and optimization of capillary dilation
There is a clear dilated bend in one of the branches of the capillary in Image 1 (Fig 1C). The size of the dilation can be seen in a plot of the radius of the capillary along its centreline, before it branches (Fig 4A). In this section an idealized model in an axisymmetric geometry (Figs 1E and 4D) is used to provide a possible physical explanation for this structure and to quantify its effect on oxygen transfer in the capillary. From a family of capillary dilations, we seek the shape that optimizes oxygen transfer over the dilated portion of the capillary.
[Figure omitted. See PDF.]
Fig 4. Results from the idealized axisymmetric model of a fetal capillary dilation.
(A) Radius of the capillary (of length L = 76 μm) in Image 1 before it branches, showing the dilation and the different ways used to define its length λ. Local minima are denoted with circles (this is subjective due to noise in the data). Points where the radius of the dilation equals the average radius are denoted with squares. The dashed lines show the minimum radius Rmin = 2.8 μm, the average radius Rav = 4.7 μm and the maximum radius Rmax = 6.8 μm. (B) Increase in oxygen transfer rate for a capillary dilation compared to a straight capillary. Included in the figure are results for shapes with 1 degree of freedom (DOF), with the undilated radius being Rmin (λ = 23 μm, Peeff = 285; blue solid line) and Rav (λ = 31 μm, Peeff = 475; red dashed line). Indicated on the figure are the maximum radius and the increase in the oxygen transfer rate N for the optimal shapes in each case with 1 DOF (squares) and 2 DOF (stars). The black dash-dotted line shows the predictions of the regression Eq (15), for a capillary dilation with 1 DOF, with the undilated radius defined as Rav. (C) Optimal shapes for successively increasing DOF. Where the shapes are asymmetric about the axial midpoint of the dilation, this is labeled on the figure. The shapes correspond to the case where the undilated radius of the capillary is defined as Rav, and are plotted as a percentage of the undilated radius. (D) Top: Axial velocity magnitude, with streamlines in white. Bottom: Concentration. Results are shown for the optimal dilation with 1 DOF, in the case where the undilated radius is taken to be Rav. It has been assumed that the centreline flow velocity in a straight capillary is equal to 300 μm/s.
https://doi.org/10.1371/journal.pone.0165369.g004
For a capillary shape parameterized with one degree of freedom (DOF), an optimal maximum radius of the dilation is found, leading to the largest amount of oxygen entering the capillary (Fig 4B and 4C). Allowing an extra DOF in the shape of the dilation, while still imposing symmetry about the tube axis and the mid-plane of the dilation, gives a slightly smaller maximum radius of the optimal shape and an enhanced oxygen transfer rate (Fig 4B and 4C). Allowing further DOF in the shape of the dilation, including asymmetric shapes, does not reduce the maximum radius of the dilation (Fig 4C). The optimal shape increases the total oxygen transfer rate N into the capillary by up to 15% above the transfer in a straight capillary (Fig 4B). The model predicts an optimal maximum radius of the dilation within approximately 33% of the radius of the dilation measured from Image 1 (Fig 4B). The regression Eq (15), with ΔP fixed, L taken to be the length of the dilated section of the capillary, and the calculated values of the vascular resistance R inserted, predicts the enhancement to N in a dilated capillary with reasonable accuracy, although it does not capture the presence of an optimal maximum radius (Fig 4B). Due to the linearity of the governing equations with respect to concentration in the axisymmetric model, a change in the oxygen concentration difference between maternal and fetal blood does not affect the model predictions. The model predictions are not affected significantly by accounting for additional oxygen transfer through the villous volume upstream and downstream of the localized capillary dilation (see S3 Appendix).
The optimal maximum radius of the localized dilation arises due to two competing effects. When the size of the dilation increases, the overall resistance of the capillary is reduced, leading to an increased flow rate (pressure drop being fixed) and therefore an increase in oxygen transfer; the increased surface area also enhances transfer. However, by conservation of mass, the local flow speed reduces where the capillary radius increases, leading to a decrease in oxygen transfer (Fig 4D). If the maximum radius of the dilation is too large, the local slowdown in the flow speed overwhelms the increased flow rate and oxygen transfer decreases overall.
Discussion
To capture the inherent stochastic variability in placental structure will require sampling of numerous villous branches to obtain distributions of placental geometries. This study demonstrates how scaling arguments, 3D simulations and idealized modeling can be combined to encapsulate key measures of flow resistance and nutrient transfer. The relevant flow characteristics that have been computed are the vascular resistance R and the vascular resistance per unit villous volume Rvil (Table 3). Using R, each capillary can be coupled to a model of the feto-placental vasculature to capture the global flow characteristics through the network. Similarly, Rvil can be used to calculate total resistances through sections of a fetal villous tree of known volume.
We have applied scaling arguments to predict (at least for the three examples considered here) which capillaries are expected to receive the most oxygen from the maternal side across a range of Péclet numbers, using only the values for total capillary length and vascular resistance. The scaling laws, along with an upper bound on the oxygen transfer rate (which can be calculated by solving for diffusion in the villous volume with the capillary assumed to be deoxygenated), capture the relationship between N and the pressure drop ΔP across parameter space where the effective Péclet number is large and linearized oxygen saturation kinetics can be assumed. These results are summarized in Table 1. The regression Eq (19) has been suggested as a way to calculate N efficiently and has been shown to capture the numerical results qualitatively. The regression equation will be useful for estimating oxygen transfer in the terminal villi in future models of the entire placenta, where it will not be feasible to solve for N numerically in every feto-placental capillary. Accurate values for the fetal oxygen transfer rate can be used to improve studies assuming the fetal villi to be oxygen sinks, such as [18, 19, 21, 23].
The identified scalings, given by Eqs (12) and (13), can also help quantify the effects of statistical variability and experimental uncertainty on blood flow and oxygen transfer. For example, tissue shrinkage or expansion can occur in villous branches for several reasons: first, due to a change in the difference between maternal and fetal blood pressure [49], for example as a result of a pathology such as pre-eclampsia; second, due to active constriction of venules or arterioles upstream and downstream of fetal capillaries, respectively [50]; finally, due to shrinkage during fixation [51]. The effect of tissue shrinkage or expansion on the data in Fig 2A has been determined by examining the dependence of the scalings Eqs (12) and (13) on the geometrical parameters. The numerical value of N from Fig 2A can be multiplied by a factor depending on the magnitude of the shrinkage or expansion factor to recover the oxygen transfer rate; a similar technique has been used to correct structural quantities in stereological studies [51]. The sensitivity of the model to changes in the parameters can also be tested using Eqs (12) and (13), with the Formula (11) for the vascular resistance. Aside from the effects of changing the geometry discussed above, the scaling Eq (12) suggests that at lower values of ΔP, the results are only affected by changes to cmat, B or μ. At higher values of ΔP, the results are expected to depend on all parameters, with the results being most sensitive to changes in the diffusion coefficient D.
The predictions of the model are sensitive to changes in specific aspects of placental structure, as a result of pathologies or differences between individuals. An idealized axisymmetric model has been used to investigate the effect of localized dilations in the feto-placental vasculature on blood flow and oxygen transfer. In the model, a localized dilation of optimal shape was found to increase oxygen transfer by up to 15%. Dilated capillaries are thought to occupy more than 35% of the total villous volume at term [5], and could therefore provide a significant enhancement to fetal oxygen uptake. The model supports the hypothesis that localized dilations develop towards term to enhance oxygen transfer in the placenta without wholesale placental growth or remodeling —in the second half of gestation, placental growth increases at a slower rate than fetal growth [1], while the mean trophoblast thickness decreases until term [52]. There may be some local remodeling of the capillary wall during the formation of the localized dilations, but we do not investigate the mechanism for their formation here. The enhancement to oxygen transfer could explain the increased dilation of fetal capillaries in pregnancies at high altitude [6, 8]; at 4300 m, maternal arterial partial pressure of oxygen is reduced by almost 50% compared to sea level [53]. A possibility is that healthy placentas at sea level are suboptimal, allowing expansion of dilations in hypoxic conditions. Increased vascularization of terminal villi has also been observed in placentas at high altitude [54], which could lead to an increase in the total number of capillary dilations. The results of the model could also provide a reason for the hypoxia-related outcomes of delayed villous maturation [14]. The lack of vasculo–scyncytial membranes observed in cases of delayed villous maturation could represent a lack of localized fetal dilations and therefore an inability of the fetus to extract an appropriate amount of oxygen from the maternal blood.
As well as characterizing the sensitivity of the model to changes in the parameters, the geometry and placental structure, it is important to understand the limitations of the modeling itself when interpreting the results. In both the 3D and axisymmetric models, simplifications have been made to the equations governing blood flow and oxygen transfer. The equations have a well-known asymptotic solution in certain limits of the effective Péclet number, allowing the scalings used in the paper to be identified. However, it will be important in future studies to relax some of the assumptions that were made in employing them. In particular, the non-Newtonian effects of red blood cells on blood flow in feto-placental capillaries should be investigated. Such effects are expected to be particularly important in small vessels and it will be vital to understand how the statistical distributions of quantities such as vascular resistance are affected by their inclusion in blood flow simulations. Non-Newtonian effects can often be incorporated into models using effective parameters or simplified equations (see [55], for example), but understanding whether such methods can be used in feto-placental capillaries of varying radius will necessitate a dedicated investigation involving simulations with discrete red blood cells.
Errors have also been introduced in converting the 3D geometries to finite-element meshes, via the neglect of the glycocalyx and by the smoothing operations performed prior to conversion, which will lead to a loss of detail in the local topography. A limitation of the current study is that the 3D images were obtained from placentas that were not perfusion-fixed. The vascular resistances calculated in this study could therefore be overestimated. It is thought that the expansion of fetal capillaries with increasing arterial blood pressure could be especially prevalent in localized dilations [56], which could mean the radii of localized dilations are particularly underestimated in Fig 4A.
To summarize, in this paper we have investigated blood flow and oxygen transfer in feto-placental capillaries theoretically, using a combination of 3D simulations, physical scaling arguments and an idealized axisymmetric model. Simulations have been performed on real geometries by converting 3D images of villous and capillary surfaces to finite-element meshes. The relationship between the total oxygen transfer rate and the pressure drop through the capillary has been captured across a wide range of pressure drops by physical scaling laws and an upper bound on the oxygen transfer rate, which has been linked directly to the geometrical parameters. A regression equation has been introduced to estimate the total oxygen transfer using the vascular resistance. Scaling arguments have been shown to be an efficient tool for quantifying the effect of statistical variability and experimental uncertainty; in particular, the identified scalings have been used to estimate the effects of tissue shrinkage or expansion in villous branches on blood flow and oxygen transfer. An idealized model in an axisymmetric geometry has been used to quantify the effect of localized capillary dilations in the fetal vasculature on oxygen transfer. The model predictions support the hypothesis that localized dilations enhance oxygen transfer towards term without the need for extensive placental growth or remodeling. The combination of 3D simulations, scaling arguments and idealized modeling used in the paper provides a blueprint for quantifying variation in measures of flow resistance and oxygen transfer as further 3D images of feto-placental villous trees become available.
Supporting Information
[Figure omitted. See PDF.]
S1 Images. Three-dimensional images.
A zipped folder containing the 3D surface meshes imported into Comsol Multiphysics for the simulations (.stl files, which can be opened in Meshlab) and the 3D tetrahedral meshes used for capillary skeletonization and calculation of the average villous distances (.bdf files, which can be opened in Gmsh).
https://doi.org/10.1371/journal.pone.0165369.s001
(ZIP)
S1 Appendix. Methods used for image analysis and three-dimensional simulations.
A description of the methods used in the image analysis and 3D simulations to increase the reproducibility of the results.
https://doi.org/10.1371/journal.pone.0165369.s002
(PDF)
S2 Appendix. Estimation of the oxygen advection enhancement parameter B.
The derivation of the parameter B quantifying the enhancement to oxygen advection due to its affinity for hemoglobin.
https://doi.org/10.1371/journal.pone.0165369.s003
(PDF)
S3 Appendix. Investigation of the effect of oxygen transfer upstream and downstream of a localized dilation.
Further numerical results generated using the idealized axisymmetric model, to quantify the effect of oxygen transfer through the surrounding villous volume upstream and downstream of the localized dilation.
https://doi.org/10.1371/journal.pone.0165369.s004
(PDF)
Acknowledgments
We acknowledge helpful discussions with Simon Cotter, Ian Crocker, Alexander Heazell and Edward Johnstone.
Author Contributions
1. Conceptualization: PP PB OEJ.
2. Data curation: PP.
3. Formal analysis: PP ILC OEJ.
4. Funding acquisition: PP ILC OEJ.
5. Investigation: PP.
6. Methodology: PP ILC OEJ.
7. Project administration: PP OEJ.
8. Resources: JJ MJ LK.
9. Software: PP.
10. Supervision: PB OEJ.
11. Validation: PP.
12. Visualization: PP.
13. Writing – original draft: PP PB OEJ.
14. Writing – review & editing: PP PB JJ MJ LK ILC OEJ.
Citation: Pearce P, Brownbill P, Janáček J, Jirkovská M, Kubínová L, Chernyavsky IL, et al. (2016) Image-Based Modeling of Blood Flow and Oxygen Transfer in Feto-Placental Capillaries. PLoS ONE 11(10): e0165369. https://doi.org/10.1371/journal.pone.0165369
1. Schneider H. Ontogenic changes in the nutritive function of the placenta. Placenta. 1996;17(1):15–26. pmid:8710808
2. Mayhew TM. Stereology and the placenta: where’s the point?–a review. Placenta. 2006;27:17–25.
3. Jirkovská M, Kubínová L, Janáček J, Moravcová M, Krejčí V, Karen P. Topological properties and spatial organization of villous capillaries in normal and diabetic placentas. J Vasc Res. 2002;39(3):268–278. pmid:12097825
4. Jirkovská M, Janáček J, Kaláb J, Kubínová L. Three-dimensional arrangement of the capillary bed and its relationship to microrheology in the terminal villi of normal term placenta. Placenta. 2008;29(10):892–897. pmid:18768220
5. Benirschke K, Kaufmann P, Baergen RN. The Pathology of the Human Placenta. 5th ed. New York, NY: Springer; 2006. https://doi.org/10.1007/b137920
6. Burton GJ, Reshetnikova OS, Milovanov AP, Teleshova OV. Stereological evaluation of vascular adaptations in human placental villi to differing forms of hypoxic stress. Placenta. 1996;17(1):49–55. pmid:8710813
7. Jackson MR, Mayhew TM, Haas JD. On the factors which contribute to thinning of the villous membrane in human placentae at high altitude. I. Thinning and regional variation in thickness of trophoblast. Placenta. 1988;9(1):1–8. pmid:3362790
8. Espinoza J, Sebire NJ, McAuliffe F, Krampl E, Nicolaides KH. Placental villus morphology in relation to maternal hypoxia at high altitude. Placenta. 2001;22(6):606–608. pmid:11440551
9. Redline RW. Cerebral palsy in term infants: a clinicopathologic analysis of 158 medicolegal case reviews. Pediatr Devel Pathol. 2008;11(6):456–464. pmid:18544009
10. Higgins M, McAuliffe FM, Mooney EE. Clinical associations with a placental diagnosis of delayed villous maturation: a retrospective study. Pediatr Devel Pathol. 2011;14(4):273–279. pmid:21329421
11. Treacy A, Higgins M, Kearney JM, McAuliffe F, Mooney EE. Delayed villous maturation of the placenta: quantitative assessment in different cohorts. Pediatr Devel Pathol. 2012;16(2):63–66. pmid:23137099
12. Redline RW. Distal villous immaturity. Diagn Histopathol. 2012;18(5):189–194.
13. Harteman JC, Nikkels PGJ, Benders MJNL, Kwee A, Groenendaal F, de Vries LS. Placental pathology in full-term infants with hypoxic-ischemic neonatal encephalopathy and association with magnetic resonance imaging pattern of brain injury. J Pediatr. 2013;163(4):968–975. pmid:23891350
14. Al-Adnani M, Marnerides A, George S, Nasir A, Weber MA. “Delayed Villous Maturation” in Placental Reporting: Concordance among Consultant Pediatric Pathologists at a Single Specialist Center. Pediatr Devel Pathol. 2015;18(5):375–379. pmid:26599068
15. Stallmach T, Hebisch G, Meier K, Dudenhausen JW, Vogel M. Rescue by birth: defective placental maturation and late fetal mortality. Obst Gynecol. 2001;97(4):505–509. pmid:11275018
16. Chernyavsky IL, Jensen OE, Leach L. A mathematical model of intervillous blood flow in the human placentone. Placenta. 2010;31(1):44–52. pmid:19945160
17. Serov AS, Salafia C, Grebenkov DS, Filoche M. The role of morphology in mathematical models of placental gas exchange. J Appl Physiol. 2016;120(10):17–28. pmid:26494446
18. Serov AS, Salafia CM, Brownbill P, Grebenkov DS, Filoche M. Optimal villi density for maximal oxygen uptake in the human placenta. J Theor Biol. 2015;364:383–396. pmid:25261730
19. Serov AS, Salafia CM, Filoche M, Grebenkov DS. An analytical theory of oxygen transport in the human placenta. J Theor Biol. 2015;368:133–144. pmid:25580015
20. Lecarpentier E, Bhatt M, Bertin GI, Deloison B, Salomon L, Deloron P, et al. Computational Fluid Dynamic Simulations of Maternal Circulation: Wall Shear Stress in the Human Placenta and Its Biological Implications. PLOS ONE. 2016;11(1). pmid:26815115
21. Chernyavsky IL, Leach L, Dryden IL, Jensen OE. Transport in the placenta: homogenizing haemodynamics in a disordered medium. Phil Trans R Soc A. 2011;369(1954):4162–4182. pmid:21969671
22. Gill JS, Salafia CM, Grebenkov D, Vvedensky DD. Modeling oxygen transport in human placental terminal villi. J Theor Biol. 2011;291:33–41. pmid:21959313
23. Mayo RP, Charnock-Jones DS, Burton GJ, Oyen ML. Three-dimensional modeling of human placental terminal villi. Placenta. 2016;
24. Clark AR, Lin M, Tawhai M, Saghian R, James JL. Multiscale modelling of the feto–placental vasculature. Interface Focus. 2015;5(2):20140078. pmid:25844150
25. Yampolsky M, Salafia CM, Shlakhter O, Haas D, Eucker B, Thorp J. Modeling the variability of shapes of a human placenta. Placenta. 2008;29(9):790–797. pmid:18674815
26. Kato Y, Oyen ML, Burton GJ. Placental villous tree models for evaluating the mechanical environment in the human placenta. In: 2014 IEEE International Conference on Imaging Systems and Techniques. IEEE; 2014.
27. Käsinger W, Huch R, Huch A. In vivo oxygen dissociation curve for whole fetal blood: fitting the Adair equation and blood gas nomogram. Scand J Clin Lab Inv. 1981;41(7):701–707.
28. Geuzaine C, Remacle JF. Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities. Int J Numer Meth Eng. 2009;79(11):1309–1331.
29. Pudney C. Distance-ordered homotopic thinning: a skeletonization algorithm for 3D digital images. Comput Vis Image Und. 1998;72(3):404–413.
30. Pries AR, Neuhaus D, Gaehtgens P. Blood viscosity in tube flow: dependence on diameter and hematocrit. Am J Physiol–Heart C. 1992;263(6):H1770–H1778. pmid:1481902
31. Banerjee RK, Kwon O, Vaidya VS, Back LH. Coupled oxygen transport analysis in the avascular wall of a coronary artery stenosis during angioplasty. J Biomech. 2008;41(2):475–479. pmid:18022178
32. Késmárky G, Kenyeres P, Rábai M, Tóth K. Plasma viscosity: a forgotten variable. Clin Hemorheol Micro. 2008;39(1):243–246. pmid:18503132
33. Beard DA, Bassingthwaighte JB. Modeling advection and diffusion of oxygen in complex vascular networks. Ann Biomed Eng. 2001;29(4):298–310. pmid:11339327
34. Higgins JM, Eddington DT, Bhatia SN, Mahadevan L. Statistical dynamics of flowing red blood cells by morphological image processing. PLoS Comput Biol. 2009;5(2):e1000288. pmid:19214200
35. Yamaguchi K, Nguyen-Phu D, Scheid P, Piiper J. Kinetics of O2 uptake and release by human erythrocytes studied by a stopped-flow technique. J Applied Physiol. 1985;58(4):1215–1224. pmid:3988677
36. Marieb EN, Hoehn K. Human Anatomy & Physiology. San Francisco, CA: Pearson Benjamin Cummings; 2007.
37. Rodesch F, Simon P, Donner C, Jauniaux E. Oxygen measurements in endometrial and trophoblastic tissues during early pregnancy. Obstet Gyneco. 1992;80(2):283–285. pmid:1635745
38. Carter AM. Placental oxygen consumption. Part I: in vivo studies–a review. Placenta. 2000;21:S31–S37. pmid:10831119
39. Rannacher R. Finite Element Methods for the Incompressible Navier-Stokes Equations. In: Fundamental Directions in Mathematical Fluid Mechanics. eds. Galdi GP, Heywood JG and Rannacher RB. Basel: Birkhäuser Basel; 2000. https://doi.org/10.1007/978-3-0348-8424-2_6
40. Amestoy PR, Duff IS, L’Excellent JY, Koster J. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM J Matrix Anal A. 2001;23(1):15–41.
41. Amestoy PR, Guermouche A, L’Excellent JY, Pralet S. Hybrid scheduling for the parallel solution of linear systems. Parallel Comput. 2006;32(1):136–156.
42. Heinrich JC, Pepper DW. Intermediate Finite Element Method: Fluid Flow and Heat Transfer Applications. Philadelphia, PA: Taylor & Francis; 1999.
43. Nocedal J, Wright S. Numerical Optimization. Berlin, Germany: Springer; 1999.
44. Cussler EL. Diffusion: Mass Transfer in Fluid Systems. Cambridge, UK: Cambridge University Press; 1997.
45. Jones AS. Heat transfer in the thermal entrance region of a flat duct. J Austral Math Soc Ser B. 1975;19(02):146–160.
46. Haase AS, Chapman SJ, Tsai PA, Lohse D, Lammertink RGH. The Graetz–Nusselt problem extended to continuum flows with finite slip. J Fluid Mech. 2015;764:R3.
47. Sapoval B, Filoche M, Weibel ER. Smaller is better–but not too small: a physical scale for the design of the mammalian pulmonary acinus. Proc Natl Acad Sci USA. 2002;99(16):10411–10416. pmid:12136124
48. Mills AF. Heat Transfer. 2nd ed. New Jersey, NJ: Prentice-Hall; 1999.
49. Karimu AL, Burton GJ. The effects of maternal vascular pressure on the dimensions of the placental capillaries. BJOG. 1994;101(1):57–63. pmid:8297870
50. Brownbill P, Sibley CP. Regulation of transplacental water transfer: the role of fetoplacental venous tone. Placenta. 2006;27(6):560–567. pmid:16256193
51. Mayhew TM, Burton GJ. Methodological problems in placental morphometry: apologia for the use of stereology based on sound sampling practice. Placenta. 1988;9(6):565–581. pmid:3070535
52. Mayhew TM, Simpson RA. Quantitative evidence for the spatial dispersal of trophoblast nuclei in human placental villi during gestation. Placenta. 1994;15(8):837–844. pmid:7886024
53. Krampl E, Lees C, Bland JM, Espinoza Dorado J, Moscoso G, Campbell S. Fetal biometry at 4300 m compared to sea level in Peru. Ultrasound Obst Gyn. 2000;16(1):9–18. pmid:11084959
54. Kingdom J, Huppertz B, Seaward G, Kaufmann P. Development of the placental villous tree and its consequences for fetal growth. Eur J Obstet Gyn R B. 2000;92(1):35–43. pmid:10986432
55. Rivera RGH, Sinha K, Graham MD. Margination regimes and drainage transition in confined multicomponent suspensions. Phys Rev Lett. 2015;114(18):188101.
56. Karimu AL, Burton GJ. Compliance of the human placental villus membrane at term: The concept of the fetoplacental unit as an autoregulating gas exchange system. Placenta. 1994;15:541–558.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016 Pearce et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
During pregnancy, oxygen diffuses from maternal to fetal blood through villous trees in the placenta. In this paper, we simulate blood flow and oxygen transfer in feto-placental capillaries by converting three-dimensional representations of villous and capillary surfaces, reconstructed from confocal laser scanning microscopy, to finite-element meshes, and calculating values of vascular flow resistance and total oxygen transfer. The relationship between the total oxygen transfer rate and the pressure drop through the capillary is shown to be captured across a wide range of pressure drops by physical scaling laws and an upper bound on the oxygen transfer rate. A regression equation is introduced that can be used to estimate the oxygen transfer in a capillary using the vascular resistance. Two techniques for quantifying the effects of statistical variability, experimental uncertainty and pathological placental structure on the calculated properties are then introduced. First, scaling arguments are used to quantify the sensitivity of the model to uncertainties in the geometry and the parameters. Second, the effects of localized dilations in fetal capillaries are investigated using an idealized axisymmetric model, to quantify the possible effect of pathological placental structure on oxygen transfer. The model predicts how, for a fixed pressure drop through a capillary, oxygen transfer is maximized by an optimal width of the dilation. The results could explain the prevalence of fetal hypoxia in cases of delayed villous maturation, a pathology characterized by a lack of the vasculo-syncytial membranes often seen in conjunction with localized capillary dilations.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer