Computational Modelling for Managing Pathways to Cartilage Failure.

Over several decades the perception and therefore description of articular cartilage changed substantially. It has transitioned from being described as a relatively inert tissue with limited repair capacity, to a tissue undergoing continuous maintenance and even adaption, through a range of complex regulatory processes. Even from the narrower lens of biomechanics, the engagement with articular cartilage has changed from it being an interesting, slippery material found in the hostile mechanical environment between opposing long bones, to an intriguing example of mechanobiology in action. The progress revealing this complexity, where physics, chemistry, material science and biology are merging, has been described with increasingly sophisticated computational models. Here we describe how these computational models of cartilage as an integrated system can be combined with the approach of structural reliability analysis. That is, causal, deterministic models placed in the framework of the probabilistic approach of structural reliability analysis could be used to understand, predict, and mitigate the risk of cartilage failure or pathology. At the heart of this approach is seeing cartilage overuse and disease processes as a 'material failure', resulting in failure to perform its function, which is largely mechanical. One can then describe pathways to failure, for example, how homeostatic repair processes can be overwhelmed leading to a compromised tissue. To illustrate this 'pathways to failure' approach, we use the interplay between cartilage consolidation and lubrication to analyse the increase in expected wear rates associated with cartilage defects or meniscectomy.

Osteochondral junction leakage and cartilage joint lubrication.

BACKGROUND AND OBJECTIVES: Previous studies have shown that there is potentially interstitial fluid exchange between cartilage tissue and the subarticular spongiosa region in the case of injury or disease (e.g., osteoarthritis and osteoporosis). Interstitial flow is also required for cartilage lubrication under joint load. A key question then is how cartilage lubrication is modified by increased interstitial fluid leakage across the osteochondral junction. Thus, the purpose of this study is to develop a numerical model to investigate changes in cartilage lubrication with changes in osteochondral junction leakage. METHODS: The multi-phase coupled model includes domains corresponding to the contact gap, cartilage tissue and subchondral bone plate region (ScBP). Each of these domains are treated as poroelastic systems, with their coupling implemented through mass and pressure continuity. The effects of osteochondral junction leakage on lubrication were investigated with a parametric study on the relative permeability between the ScBP and cartilage tissue. RESULTS: Significant effects of ScBP permeability were predicted, especially during the early stage of the junction leakage development (early stage of the disease). There is a significant reduction in mixed-mode lubrication duration under the effect of increased junction leakage (the cartilage tissue mixed-mode lubrication duration is about 33% decrease for a relative permeability ratio of 0.1 between ScBP and cartilage tissue, and about 52% decrease under the osteoarthritis condition). In addition, the time for cartilage to reach steady-state consolidation is significantly reduced when ScBP permeability increases (the consolidation time reduces from roughly 2 h to 1.2 h when the relative permeability ratio increases from 0.001 to 0.1, and it reduces to 0.8 h for an advanced osteoarthritis condition). It is predicted that the initial friction coefficient could increase by over 60% when the ScBP permeability is consistent with an advanced osteoarthritis (OA) condition. CONCLUSION: Increased osteochondral junction leakage induced by joint injury and disease could result in increased cartilage surface wear rates due to more rapid interstitial fluid depressurization within articular cartilage.

Spatial composition and turnover of the main molecules in the adult glomerular basement membrane.

The glomerular basement membrane (GBM) is an important tissue structure in kidney function. It is the membrane through which filtrate and solutes must pass to reach the nephron tubules. This review focuses on the spatial location of the main extracellular matrix components of the GBM. It also attempts to explain this organization in terms of their synthesis, transport, and loss. The picture that emerges is that the collagen IV and laminin content of GBM are in a very slow dynamic disequilibrium, leading to GBM thickening with age, and in contrast, some heparan sulfate proteoglycans are in a dynamic equilibrium with a very rapid turnover (i.e. half-life measured in ~hours) and flow direction against the flow of filtrate. The highly rapid heparan sulfate turnover may serve several roles, including an unclogging mechanism for the GBM, compressive stiffness of the GBM fiber network, and/or enabling podocycte-endothelial crosstalk against the flow of filtrate.

Predicting oxygen tension along the ureter.

Continuous measurement of bladder urine oxygen tension (Po2) is a method to potentially detect renal medullary hypoxia in patients at risk of acute kidney injury (AKI). To assess its practicality, we developed a computational model of the peristaltic movement of a urine bolus along the ureter and the oxygen exchange between the bolus and ureter wall. This model quantifies the changes in urine Po2 as urine transits from the renal pelvis to the bladder. The model parameters were calibrated using experimental data in rabbits, such that most of the model predictions are within ±1 SE of the reported mean in the experiment, with the average percent difference being 7.0%. Based on parametric experiments performed using a model scaled to the geometric dimensions of a human ureter, we found that bladder urine Po2 is strongly dependent on the bolus volume (i.e., bolus volume-to-surface area ratio), especially at a volume less than its physiological (baseline) volume (<0.2 mL). For the model assumptions, changes in peristaltic frequency resulted in a minimal change in bladder urine Po2 (<1 mmHg). The model also predicted that there exists a family of linear relationships between the bladder-urine Po2 and pelvic urine Po2 for different input conditions. We conclude that it may technically be possible to predict renal medullary Po2 based on the measurement of bladder urine Po2, provided that there are accurate real-time measurements of model input parameters.NEW & NOTEWORTHY Measurement of bladder urine oxygen tension has been proposed as a new method to potentially detect the risk of acute kidney injury in patients. A computational model of oxygen exchange between urine bolus and ureteral tissue shows that it may be technically possible to determine the risk of acute kidney injury based on the measurement of bladder urine oxygen tension, provided that the measurement data are properly interpreted via a computational model.

A Probabilistic Failure Risk Approach to The Problem of Articular Cartilage Lubrication.

BACKGROUND AND OBJECTIVE: The geometrical and mechanical properties that characterise the cartilage contact gap are uncertain and spatially varied. To date the effects of such uncertainties on cartilage lubrication have not been explored. Using a probabilistic approach, the purpose of this study is to numerically investigate the influence of surficial cartilage glycoaminoglycan (GAG) content on joint lubrication behaviour. Gap asperity stiffness and polymer brush border (PBB) thickness are affected by the uncertainty of surficial GAG concentration, and so their correlated effects in maintaining hydrodynamic joint lubrication are investigated. METHODS: Correlated sampling data are first generated by Monte Carlo simulation. These data are used as inputs for the cartilage contact model, which includes three distinctive features of cartilage tissue (tension-compression nonlinearity, aggrecan dependent permeability and compressive modulus) and fluid flow resistance effects of PBB on cartilage surface. The degree of hydrodynamic lubrication after thirty minutes of constant loading is used as an indicator for assessing the lubrication performance at the contact interface. RESULTS: The increase of PBB thickness with GAG concentration enhances the hydrodynamic lubrication component in the cartilage contact gap, whereas increasing the asperity stiffness with GAG concentration impairs hydrodynamic lubrication. GAG loss rate increases with the rise of GAG concentration. More aggrecan shedding through the surface could result in a thicker and denser PBB, and therefore enhance the lubrication performance in mixed-mode regime. On the other hand, higher GAG content makes the asperities stiffer, which may impede contact gap closure, and thus encourage gap fluid loss and impair the lubrication performance of cartilage. CONCLUSION: The lubrication performance of cartilage varies with the physiological conditions of the joint. Since a range of variables are internally related, the outcomes on joint lubrication are difficult to predict. A probabilistic approach accounting for the uncertainties can potentially result in more accurate evaluations of joint lubrication performance.

Estimating outflow facility parameters for the human eye using hypotensive pressure-time data.

We have previously developed a new theory for pressure dependent outflow from the human eye, and tested the model using experimental data at intraocular pressures above normal eye pressures. In this paper, we use our model to analyze a hypotensive pressure-time dataset obtained following application of a Honan balloon. Here we show that the hypotensive pressure-time data can be successfully analyzed using our proposed pressure dependent outflow model. When the most uncertain initial data point is removed from the dataset, then parameter estimates are close to our previous parameter estimates, but clearly parameter estimates are very sensitive to assumptions. We further show that (i) for a measured intraocular pressure-time curve, the estimated model parameter for whole eye surface hydraulic conductivity is primarily a function of the ocular rigidity, and (ii) the estimated model parameter that controls the rate of decrease of outflow with increasing pressure is primarily a function of the convexity of the monotonic pressure-time curve. Reducing parameter uncertainty could be accomplished using new technologies to obtain higher quality datasets, and by gathering additional data to better define model parameter ranges for the normal eye. With additional research, we expect the pressure dependent outflow analysis described herein may find applications in the differential diagnosis, prognosis and monitoring of the glaucomatous eye.

A cardiovascular model for renal perfusion during cardiopulmonary bypass surgery.

Acute kidney injury (AKI) is a major complication following cardiac surgery requiring cardiopulmonary bypass (CPB). It is likely that poor renal perfusion contributes to the occurrence of AKI, via renal hypoxia, so it is imperative to maintain optimal renal perfusion during CPB. We have developed a straightforward cardiovascular perfusion model with parameter values calibrated against experimental and/or clinical data from several independent studies of CPB in humans and animals. Following model development and calibration, we performed a one-at-a-time parametric study to investigate the response of renal perfusion to several variables during CPB, namely pump flow (denoted CO for 'cardiac output'), renal vascular resistance, and non-renal vascular resistance. From the parametric study, we have found that all three parameters had a similarly strong influence on renal perfusion. We simulated three potential strategies for maintaining optimum renal perfusion during CPB and tested their effectiveness. The strategies were: (1) increasing the pump flow; (2) administrating noradrenaline (vasopressor); and (3) administrating fenoldopam (renal vasodilator). Simulations have revealed that administration of fenoldopam is likely to be the most effective of the three strategies. Other findings from our simulations are that increasing pump flow is less effective when central venous pressure is elevated. Further, renal autoregulation is likely inoperative during CPB, as evidenced by an unchanging renal vascular resistance with increasing CO and blood pressure. The cardiac-renal perfusion model developed in this study can be linked with other kidney models to simulate the changes in renal oxygenation during CPB.

No flow through the vitreous humor: How strong is the evidence?

When analyzing vitreal drug delivery, or the pharmacological effects of drugs on intraocular pressure, or when interpreting outflow facility measurements, it is generally accepted that the fluid in the vitreous humor is stagnant. It is accepted that for all practical purposes, the aqueous fluid exits the eye via anterior pathways only, and so there is negligible if any posteriorly directed flow of aqueous through the vitreous humor. This assumption is largely based on the interpretation of experimental data from key sources including Maurice (1957), Moseley (1984), Gaul and Brubaker (1986), Maurice (1987) and Araie et al. (1991). However, there is strong independent evidence suggesting there is a substantial fluid flow across the retinal pigment epithelium from key sources including Cantrill and Pederson (1984), Chihara and Nao-i, Tsuboi (1985), Dahrouj et al. (2014), Smith and Gardiner (2017) and Smith et al. (2019). The conflicting evidence creates a conundrum-how can both interpretations be true? This leads us to re-evaluate the evidence. We demonstrate that the data believed to be supporting no aqueous flow through the vitreous are in fact compatible with a significant normal aqueous flow. We identify strong and independent lines of evidence supporting fluid flow across the RPE, including our new outflow model for the eye. On balance it appears the current evidence favors the view that there is normally a significant aqueous flow across the RPE in vivo. This finding suggests that past and future analyses of outflow facility, interpretations of some drug distributions and the interpretation of some drug effects on eye tissues, may need to be revised.

Renal oxygenation: From data to insight.

Computational models have made a major contribution to the field of physiology. As the complexity of our understanding of biological systems expands, the need for computational methods only increases. But collaboration between experimental physiologists and computational modellers (ie theoretical physiologists) is not easy. One of the major challenges is to break down the barriers created by differences in vocabulary and approach between the two disciplines. In this review, we have two major aims. Firstly, we wish to contribute to the effort to break down these barriers and so encourage more interdisciplinary collaboration. So, we begin with a "primer" on the ways in which computational models can help us understand physiology and pathophysiology. Second, we aim to provide an update of recent efforts in one specific area of physiology, renal oxygenation. This work is shedding new light on the causes and consequences of renal hypoxia. But as importantly, computational modelling is providing direction for experimental physiologists working in the field of renal oxygenation by: (a) generating new hypotheses that can be tested in experimental studies, (b) allowing experiments that are technically unfeasible to be simulated in silico, or variables that cannot be measured experimentally to be estimated, and (c) providing a means by which the quality of experimental data can be assessed. Critically, based on our experience, we strongly believe that experimental and theoretical physiology should not be seen as separate exercises. Rather, they should be integrated to permit an iterative process between modelling and experimentation.

What Makes the Kidney Susceptible to Hypoxia?

Per gram of tissue, the kidneys are among our most highly perfused organs. Yet the renal cortex and, in particular, the renal medulla are susceptible to hypoxia. In turn, hypoxia is a major pathophysiological feature of both acute kidney injury and chronic kidney disease. We identify seven factors that render the kidney susceptible to hypoxia: (1) the large metabolic demand imposed by active reabsorption of sodium; (2) limitations on oxygen delivery to cortical tissue imposed by the density of peritubular capillaries; (3) the poor capacity for angiogenesis in the adult kidney; (4) the limited ability of the renal vasculature to dilate in response to hypoxia; (5) diffusive oxygen shunting between arteries and veins in the cortex and descending and ascending vasa recta in the medulla; (6) the physiological requirement for low medullary blood flow to facilitate urinary concentration; and (7) the topography of vascular-tubular arrangements in the outer medulla that limit oxygen delivery to the thick ascending limb of Henle's loop. Recent collaborative efforts between anatomists, physiologists, and mathematicians have improved our understanding of the roles of these factors in both physiological regulation of intrarenal oxygenation and development of renal hypoxia under pathophysiological conditions. We are also better able to understand these apparent maladaptations in the context of evolution. That is, they can be explained by the combined effects of historical contingency (our ancestral life in the sea) and selection pressures imposed by the multiple functions of the kidney to regulate extracellular fluid volume, retain water, and control erythrocyte production.

Analysis of the critical determinants of renal medullary oxygenation.

We have previously developed a three-dimensional computational model of oxygen transport in the renal medulla. In the present study, we used this model to quantify the sensitivity of renal medullary oxygenation to four of its major known determinants: medullary blood flow (MBF), medullary oxygen consumption rate (V̇o2,M), hemoglobin (Hb) concentration in the blood, and renal perfusion pressure. We also examined medullary oxygenation under special conditions of hydropenia, extracellular fluid volume expansion by infusion of isotonic saline, and hemodilution during cardiopulmonary bypass. Under baseline (normal) conditions, the average medullary tissue Po2 predicted for the whole renal medulla was ~30 mmHg. The periphery of the interbundle region in the outer medulla was identified as the most hypoxic region in the renal medulla, which demonstrates that the model prediction is qualitatively accurate. Medullary oxygenation was most sensitive to changes in renal perfusion pressure followed by Hb, MBF, and V̇o2,M, in that order. The medullary oxygenation also became sensitized by prohypoxic changes in other parameters, leading to a greater fall in medullary tissue Po2 when multiple parameters changed simultaneously. Hydropenia did not induce a significant change in medullary oxygenation compared with the baseline state, while volume expansion resulted in a large increase in inner medulla tissue Po2 (by ~15 mmHg). Under conditions of cardiopulmonary bypass, the renal medulla became severely hypoxic, due to hemodilution, with one-third of the outer stripe of outer medulla tissue having a Po2 of <5 mmHg.

The investigation of fluid flow in cartilage contact gap.

Synovial fluid flow in articular joint capsule plays an important role during mixed mode lubrication. However, the actual fluid flow behaviour during cartilage contact has not been fully understood so far. This is due to the difficulties in measuring the gap permeability using conventional experimental techniques. The problem becomes further complicated with consideration of the cartilage surface roughness. Here a validated numerical study was developed to quantify the gap permeability of lateral synovial fluid flow. Both macro- and micro-scale gap flow models were created based on Darcy's law at the macro-scale and the Navier-stokes equation at the micro-scale. To generate model inputs, the cartilage topography was numerically synthesised based on the experimental measurements of bovine medial tibia cartilage surface roughness using Dektak Stylus Profilers. The experimental results show that the average roughness height Ra is 1.97 µm and root-mean-square roughness height Rq is 2.44 µm, while the correlation lengths of the secondary and tertiary undulations are round 100 µm and 20 µm, respectively. The numerical results indicate that the contact gap height and fluid pressure gradient are two critical parameters which significantly affect the gap permeability. As the contact gap closes, there is a decrease in gap permeability, and most importantly, the gap permeability is also very sensitive to the fluid pressure gradient. Furthermore, with gap closure, the permeability of the contact gap gradually approaches that of the cartilage tissue, at which point the contact gap is functional closed. This occurs at a contact gap height around 1 µm and fluid pressure gradient below 5 × 105 Pa/m in this study.

Estimating three-dimensional outflow and pressure gradients within the human eye.

In this paper we set the previously reported pressure-dependent, ordinary differential equation outflow model by Smith and Gardiner for the human eye, into a new three-dimensional (3D) porous media outflow model of the eye, and calibrate model parameters using data reported in the literature. Assuming normal outflow through anterior pathways, we test the ability of 3D flow model to predict the pressure elevation with a silicone oil tamponade. Then assuming outflow across the retinal pigment epithelium is normal, we test the ability of the 3D model to predict the pressure elevation in Schwartz-Matsuo syndrome. For the first time we find the flow model can successfully model both conditions, which helps to build confidence in the validity and accuracy of the 3D pressure-dependent outflow model proposed here. We employ this flow model to estimate the translaminar pressure gradient within the optic nerve head of a normal eye in both the upright and supine postures, and during the day and at night. Based on a ratio of estimated and measured pressure gradients, we define a factor of safety against acute interruption of axonal transport at the laminar cribrosa. Using a completely independent method, based on the behaviour of dynein molecular motors, we compute the factor of safety against stalling the dynein molecule motors, and so compromising retrograde axonal transport. We show these two independent methods for estimating factors of safety agree reasonably well and appear to be consistent. Taken together, the new 3D pressure-dependent outflow model proves itself to capable of providing a useful modeling platform for analyzing eye behaviour in a variety of physiological and clinically useful contexts, including IOP elevation in Schwartz-Matsuo syndrome and with silicone oil tamponade, and potentially for risk assessment for optic glaucomatous neuropathy.

Universal hidden order in amorphous cellular geometries.

Partitioning space into cells with certain extreme geometrical properties is a central problem in many fields of science and technology. Here we investigate the Quantizer problem, defined as the optimisation of the moment of inertia of Voronoi cells, i.e., similarly-sized 'sphere-like' polyhedra that tile space are preferred. We employ Lloyd's centroidal Voronoi diagram algorithm to solve this problem and find that it converges to disordered states associated with deep local minima. These states are universal in the sense that their structure factors are characterised by a complete independence of a wide class of initial conditions they evolved from. They moreover exhibit an anomalous suppression of long-wavelength density fluctuations and quickly become effectively hyperuniform. Our findings warrant the search for novel amorphous hyperuniform phases and cellular materials with unique physical properties.

Stimulation of erythropoietin release by hypoxia and hypoxemia: similar but different.

Erythropoietin is released from the kidney in response to tissue hypoxia. Montero and Lundby found that increases in plasma erythropoietin induced by reducing arterial oxygen content in healthy humans were independent of arterial oxygen tension. Their observations accord with the established physiology of kidney oxygenation and can be predicted by a computational model of renal oxygen transport. However, model simulations indicate that the interpretation implicit in the title of their paper may be an oversimplification.

A model of oxygen transport in the rat renal medulla.

The renal medulla is prone to hypoxia. Medullary hypoxia is postulated to be a leading cause of acute kidney injury, so there is considerable interest in predicting the oxygen tension in the medulla. Therefore we have developed a computational model for blood and oxygen transport within a physiologically normal rat renal medulla, using a multilevel modeling approach. For the top-level model we use the theory of porous media and advection-dispersion transport through a realistic three-dimensional representation of the medulla's gross anatomy to describe blood flow and oxygen transport throughout the renal medulla. For the lower-level models, we employ two-dimensional reaction-diffusion models describing the distribution of oxygen through tissue surrounding the vasculature. Steady-state model predictions at the two levels are satisfied simultaneously, through iteration between the levels. The computational model was validated by simulating eight sets of experimental data regarding renal oxygenation in rats (using 4 sets of control groups and 4 sets of treatment groups, described in 4 independent publications). Predicted medullary tissue oxygen tension or microvascular oxygen tension for control groups and for treatment groups that underwent moderate perturbation in hemodynamic and renal functions is within ±2 SE values observed experimentally. Diffusive shunting between descending and ascending vasa recta is predicted to be only 3% of the oxygen delivered. The validation tests confirm that the computational model is robust and capable of capturing the behavior of renal medullary oxygenation in both normal and early-stage pathological states in the rat.

Effect of collagen length distribution and timing for repair on the active TGF-ß concentration in tendon.

The composition of extracellular matrix (ECM) in tendon depends on the secretion profile of resident cells known as tenocytes. For tissues with a mechanical role like tendon, mechanical strain is known to play an important role in determining the secretion profile of resident cells. Previously we explored the idea of estimating average concentrations of ECM molecules as a function of tendon strain magnitude and number of loading cycles. Specifically, we developed a model of the mechanical fatigue damage of tendon collagen fibers and introduced elementary cell responses (ECRs) by which local cellular-level responses to the strain environment, combined with the fatigue damage model, were scaled up to predict tissue-level responses. Using this approach, we demonstrated that the proposed model is capable of estimating average concentrations of ECM molecules that qualitatively accord with experimental observations. In this study, we increase model realism by extending this approach to consider the implications of a non-uniform collagen fiber distribution, and the influence of time delay on repair of damaged collagen fibers. Using this approach, we focus the study on the average tenocyte secretion profile for active transforming growth factor beta (TGF-ß), and discover that increasing fiber length dispersion and/or increasing repair delay leads to increasing active TGF-ß concentrations, and reduced sensitivity of average concentration profile of TGF-ß to tendon strain.

Three-dimensional morphometric analysis of the renal vasculature.

Vascular topology and morphology are critical in the regulation of blood flow and the transport of small solutes, including oxygen, carbon dioxide, nitric oxide, and hydrogen sulfide. Renal vascular morphology is particularly challenging, since many arterial walls are partially wrapped by the walls of veins. In the absence of a precise characterization of three-dimensional branching vascular geometry, accurate computational modeling of the intrarenal transport of small diffusible molecules is impossible. An enormous manual effort was required to achieve a relatively precise characterization of rat renal vascular geometry, highlighting the need for an automated method for analysis of branched vasculature morphology to allow characterization of the renal vascular geometry of other species, including humans. We present a semisupervised method for three-dimensional morphometric analysis of renal vasculature images generated by computed tomography. We derive quantitative vascular attributes important to mass transport between arteries, veins, and the renal tissue and present methods for their computation for a three-dimensional vascular geometry. To validate the algorithm, we compare automated vascular estimates with subjective manual measurements for a portion of rabbit kidney. Although increased image resolution can improve outcomes, our results demonstrate that the method can quantify the morphological characteristics of artery-vein pairs, comparing favorably with manual measurements. Similar to the rat, we show that rabbit artery-vein pairs become less intimate along the course of the renal vasculature, but the total wrapped mass transfer coefficient increases and then decreases. This new method will facilitate new quantitative physiological models describing the transport of small molecules within the kidney.

Estimating outflow facility through pressure dependent pathways of the human eye.

We develop and test a new theory for pressure dependent outflow from the eye. The theory comprises three main parameters: (i) a constant hydraulic conductivity, (ii) an exponential decay constant and (iii) a no-flow intraocular pressure, from which the total pressure dependent outflow, average outflow facilities and local outflow facilities for the whole eye may be evaluated. We use a new notation to specify precisely the meaning of model parameters and so model outputs. Drawing on a range of published data, we apply the theory to animal eyes, enucleated eyes and in vivo human eyes, and demonstrate how to evaluate model parameters. It is shown that the theory can fit high quality experimental data remarkably well. The new theory predicts that outflow facilities and total pressure dependent outflow for the whole eye are more than twice as large as estimates based on the Goldman equation and fluorometric analysis of anterior aqueous outflow. It appears likely that this discrepancy can be largely explained by pseudofacility and aqueous flow through the retinal pigmented epithelium, while any residual discrepancy may be due to pathological processes in aged eyes. The model predicts that if the hydraulic conductivity is too small, or the exponential decay constant is too large, then intraocular eye pressure may become unstable when subjected to normal circadian changes in aqueous production. The model also predicts relationships between variables that may be helpful when planning future experiments, and the model generates many novel testable hypotheses. With additional research, the analysis described here may find application in the differential diagnosis, prognosis and monitoring of glaucoma.

Micro-computed tomographic analysis of the radial geometry of intrarenal artery-vein pairs in rats and rabbits: Comparison with light microscopy.

We assessed the utility of synchrotron-radiation micro-computed tomography (micro-CT) for quantification of the radial geometry of the renal cortical vasculature. The kidneys of nine rats and six rabbits were perfusion fixed and the renal circulation filled with Microfil. In order to assess shrinkage of Microfil, rat kidneys were imaged at the Australian Synchrotron immediately upon tissue preparation and then post fixed in paraformaldehyde and reimaged 24 hours later. The Microfil shrank only 2-5% over the 24 hour period. All subsequent micro-CT imaging was completed within 24 hours of sample preparation. After micro-CT imaging, the kidneys were processed for histological analysis. In both rat and rabbit kidneys, vascular structures identified in histological sections could be identified in two-dimensional (2D) micro-CT images from the original kidney. Vascular morphology was similar in the two sets of images. Radial geometry quantified by manual analysis of 2D images from micro-CT was consistent with corresponding data generated by light microscopy. However, due to limited spatial resolution when imaging a whole organ using contrast-enhanced micro-CT, only arteries ≥100 and ≥60 µm in diameter, for the rat and rabbit respectively, could be assessed. We conclude that it is feasible and valid to use micro-CT to quantify vascular geometry of the renal cortical circulation in both the rat and rabbit. However, a combination of light microscopic and micro-CT approaches are required to evaluate the spatial relationships between intrarenal arteries and veins over an extensive range of vessel size.