Architecture of the rat nephron-arterial network: analysis with micro-computed tomography

Mapping the blood vasculature in an intact human kidney using hierarchical phase-contrast tomography

The architecture of the kidney vasculature is essential for its function. Although structural profiling of the intact rodent kidney vasculature has been performed, it is challenging to map vascular architecture of larger human organs. We hypothesised that hierarchical phase-contrast tomography (HiP-CT) would enable quantitative analysis of the entire human kidney vasculature. Combining label-free HiP-CT imaging of an intact kidney from a 63-year-old male with topology network analysis, we quantitated vasculature architecture in the human kidney down to the scale of arterioles. Although human and rat kidney vascular topologies are comparable, vascular radius decreases at a significantly faster rate in humans as vessels branch from artery towards the cortex. At branching points of large vessels, radii are theoretically optimised to minimise flow resistance, an observation not found for smaller arterioles. Structural differences in the vasculature were found in different spatial zones of the kidney reflecting their unique functional roles. Overall, this represents the first time the entire arterial vasculature of a human kidney has been mapped providing essential inputs for computational models of kidney vascular flow and synthetic vascular architectures, with implications for understanding how the structure of individual blood vessels collectively scales to facilitate organ function.

Resting-state MRI reveals spontaneous physiological fluctuations in the kidney and tracks diabetic nephropathy in rats

The kidneys maintain fluid-electrolyte balance and excrete waste in the presence of constant fluctuations in plasma volume and systemic blood pressure. The kidneys perform these functions to control capillary perfusion and glomerular filtration by modulating the mechanisms of autoregulation. An effect of these modulations are spontaneous, natural fluctuations in glomerular perfusion. Numerous other mechanisms can lead to fluctuations in perfusion and flow. The ability to monitor these spontaneous physiological fluctuations in vivo could facilitate the early detection of kidney disease. The goal of this work was to investigate the use of resting-state magnetic resonance imaging (rsMRI) to detect spontaneous physiological fluctuations in the kidney. We performed rsMRI of rat kidneys in vivo over 10 min, applying motion correction to resolve time series in each voxel. We observed spatially variable, spontaneous fluctuations in rsMRI signal between 0 and 0.3 Hz, in frequency bands associated with autoregulatory mechanisms. We further applied rsMRI to investigate changes in these fluctuations in a rat model of diabetic nephropathy. Spectral analysis was performed on time series of rsMRI signals in the kidney cortex and medulla. The power from spectra in specific frequency bands from the cortex correlated with severity of glomerular pathology caused by diabetic nephropathy. Finally, we investigated the feasibility of using rsMRI of the human kidney in two participants, observing the presence of similar, spatially variable fluctuations. This approach may enable a range of preclinical and clinical investigations of kidney function and facilitate the development of new therapies to improve outcomes in patients with kidney disease.NEW & NOTEWORTHY This work demonstrates the development and use of resting-state MRI to detect low-frequency, spontaneous physiological fluctuations in the kidney consistent with previously observed fluctuations in perfusion and potentially due to autoregulatory function. These fluctuations are detectable in rat and human kidneys, and the power of these fluctuations is affected by diabetic nephropathy in rats.

Interacting information streams on the nephron arterial network

Blood flow and glomerular filtration in the kidney are regulated by two mechanisms acting on the afferent arteriole of each nephron. The two mechanisms operate as limit cycle oscillators, each responding to a different signal. The myogenic mechanism is sensitive to a transmural pressure difference across the wall of the arteriole, and tubuloglomerular feedback (TGF) responds to the NaCl concentration in tubular fluid flowing into the nephron's distal tubule,. The two mechanisms interact with each other, synchronize, cause oscillations in tubular flow and pressure, and form a bimodal electrical signal that propagates into the arterial network. The electrical signal enables nephrons adjacent to each other in the arterial network to synchronize, but non-adjacent nephrons do not synchronize. The arteries supplying the nephrons have the morphologic characteristics of a rooted tree network, with 3 motifs characterizing nephron distribution. We developed a model of 10 nephrons and their afferent arterioles in an arterial network that reproduced these structural characteristics, with half of its components on the renal surface, where experimental data suitable for model validation is available, and the other half below the surface, from which no experimental data has been reported. The model simulated several interactions: TGF-myogenic in each nephron with TGF modulating amplitude and frequency of the myogenic oscillation; adjacent nephron-nephron with strong coupling; non-adjacent nephron-nephron, with weak coupling because of electrical signal transmission through electrically conductive arterial walls; and coupling involving arterial nodal pressure at the ends of each arterial segment, and between arterial nodes and the afferent arterioles originating at the nodes. The model predicted full synchronization between adjacent nephrons pairs and partial synchronization among weakly coupled nephrons, reproducing experimental findings. The model also predicted aperiodic fluctuations of tubular and arterial pressures lasting longer than TGF oscillations in nephrons, again confirming experimental observations. The model did not predict complete synchronization of all nephrons.

A hybrid approach to full-scale reconstruction of renal arterial network

The renal vasculature, acting as a resource distribution network, plays an important role in both the physiology and pathophysiology of the kidney. However, no imaging techniques allow an assessment of the structure and function of the renal vasculature due to limited spatial and temporal resolution. To develop realistic computer simulations of renal function, and to develop new image-based diagnostic methods based on artificial intelligence, it is necessary to have a realistic full-scale model of the renal vasculature. We propose a hybrid framework to build subject-specific models of the renal vascular network by using semi-automated segmentation of large arteries and estimation of cortex area from a micro-CT scan as a starting point, and by adopting the Global Constructive Optimization algorithm for generating smaller vessels. Our results show a close agreement between the reconstructed vasculature and existing anatomical data obtained from a rat kidney with respect to morphometric and hemodynamic parameters.

Numerical simulation of the vascular structure dependence of blood flow in the kidney

A numerical simulation was performed to clarify renal blood flow determination by the vascular structures. Large and small vessels were modeled as symmetric and asymmetric branching vessels, respectively, with simple geometries to parameterize the vascular structures. Modeling individual vessels as straight pipes, Murray's law was used to determine the vessel diameters. Blood flow in the vascular structure was calculated by network analysis based on Hagen-Poiseuille's law. Blood flow simulations for a vascular network segment demonstrated that blood flow rate and pressure vary within the same-generation vessels because of an asymmetric vessel branch while they generally tend to decrease with vessel diameter; thus, the standard deviation of flow rate relative to the mean (relative standard deviation [RSD]) increased from 0.4 to 1.0 when the number of the daughter vessels increased from 3 to 10. Blood flow simulations for an entire vascular network of a kidney showed that the vessel number and branching style, rather than Strahler order, are major parameters in successfully reproducing renal blood flow measured in published experiments. The entire vascular network could generate variation in the physiological flow rate in afferent arterioles at 0.2-0.38 in RSD, which is at least compatible with 0.16 by diameter variation within the same-generation vessels.

Numerical Modeling and Simulation of Blood Flow in a Rat Kidney: Coupling of the Myogenic Response and the Vascular Structure

A numerical simulation was carried out to investigate the blood flow behavior (i.e., flow rate and pressure) and coupling of a renal vascular network and the myogenic response to various conditions. A vascular segment and an entire kidney vascular network were modeled by assuming one single vessel as a straight pipe whose diameter was determined by Murray’s law. The myogenic response was tested on individual AA (afferent artery)–GC (glomerular capillaries)–EA (efferent artery) systems, thereby regulating blood flow throughout the vascular network. Blood flow in the vascular structure was calculated by network analysis based on Hagen–Poiseuille’s law to various boundary conditions. Simulation results demonstrated that, in the vascular segment, the inlet pressure Pinlet and the vascular structure act together on the myogenic response of each individual AA–GC–EA subsystem, such that the early-branching subsystems in the vascular network reached the well-regulated state first, with an interval of the inlet as Pinlet = 10.5–21.0 kPa, whereas the one that branched last exhibited a later interval with Pinlet = 13.0–24.0 kPa. In the entire vascular network, in contrast to the Pinlet interval (13.0–20.0 kPa) of the unified well-regulated state for all AA–GC–EA subsystems of the symmetric model, the asymmetric model exhibited the differences among subsystems with Pinlet ranging from 12.0–17.0 to 16.0–20.0 kPa, eventually achieving a well-regulated state of 13.0–18.5 kPa for the entire kidney. Furthermore, when Pinlet continued to rise (e.g., 21.0 kPa) beyond the vasoconstriction range of the myogenic response, high glomerular pressure was also related to vascular structure, where PGC of early-branching subsystems was 9.0 kPa and of late-branching one was 7.5 kPa. These findings demonstrate how the myogenic response regulates renal blood flow in vascular network system that comprises a large number of vessel elements.

Synchronization in renal microcirculation unveiled with high-resolution blood flow imaging

Internephron interaction is fundamental for kidney function. Earlier studies have shown that nephrons signal to each other, synchronise over short distances, and potentially form large synchronised clusters. Such clusters would play an important role in renal autoregulation, but due to the technological limitations, their presence is yet to be confirmed. In the present study, we introduce an approach for high-resolution laser speckle imaging of renal blood flow and apply it to estimate frequency and phase differences in rat kidney microcirculation under different conditions. The analysis unveiled spatial and temporal evolution of synchronised blood flow clusters of various sizes, including the formation of large (>90 vessels) long-lived clusters (>10 periods) locked at the frequency of the tubular glomerular feedback mechanism. Administration of vasoactive agents caused significant changes in the synchronisation patterns and, thus, in nephrons' co-operative dynamics. Specifically, infusion of vasoconstrictor angiotensin II promoted stronger synchronisation, while acetylcholine caused complete desynchronisation. The results confirm the presence of the local synchronisation in the renal microcirculatory blood flow and that it changes depending on the condition of the vascular network and the blood pressure, which will have further implications for the role of such synchronisation in pathologies development.

Evaluation of 2D super-resolution ultrasound imaging of the rat renal vasculature using ex vivo micro-computed tomography

Super-resolution ultrasound imaging (SRUS) enables in vivo microvascular imaging of deeper-lying tissues and organs, such as the kidneys or liver. The technique allows new insights into microvascular anatomy and physiology and the development of disease-related microvascular abnormalities. However, the microvascular anatomy is intricate and challenging to depict with the currently available imaging techniques, and validation of the microvascular structures of deeper-lying organs obtained with SRUS remains difficult. Our study aimed to directly compare the vascular anatomy in two in vivo 2D SRUS images of a Sprague-Dawley rat kidney with ex vivo μCT of the same kidney. Co-registering the SRUS images to the μCT volume revealed visually very similar vascular features of vessels ranging from ~ 100 to 1300 μm in diameter and illustrated a high level of vessel branching complexity captured in the 2D SRUS images. Additionally, it was shown that it is difficult to use μCT data of a whole rat kidney specimen to validate the super-resolution capability of our ultrasound scans, i.e., validating the actual microvasculature of the rat kidney. Lastly, by comparing the two imaging modalities, fundamental challenges for 2D SRUS were demonstrated, including the complexity of projecting a 3D vessel network into 2D. These challenges should be considered when interpreting clinical or preclinical SRUS data in future studies.

Tubuloglomerular Feedback Synchronization in Nephrovascular Networks

To perform their functions, the kidneys maintain stable blood perfusion in the face of fluctuations in systemic BP. This is done through autoregulation of blood flow by the generic myogenic response and the kidney-specific tubuloglomerular feedback (TGF) mechanism. The central theme of this paper is that, to achieve autoregulation, nephrons do not work as single units to manage their individual blood flows, but rather communicate electrically over long distances to other nephrons via the vascular tree. Accordingly, we define the nephrovascular unit (NVU) to be a structure consisting of the nephron, glomerulus, afferent arteriole, and efferent arteriole. We discuss features that require and enable distributed autoregulation mediated by TGF across the kidney. These features include the highly variable topology of the renal vasculature which creates variability in circulation and the potential for mismatch between tubular oxygen demand and delivery; the self-sustained oscillations in each NVU arising from the autoregulatory mechanisms; and the presence of extensive gap junctions formed by connexins and their properties that enable long-distance transmission of TGF signals. The existence of TGF synchronization across the renal microvascular network enables an understanding of how NVUs optimize oxygenation-perfusion matching while preventing transmission of high systemic pressure to the glomeruli, which could lead to progressive glomerular and vascular injury.

Renal microvasculature in the adult pipid frog, Xenopus laevis: A scanning electron microscope study of vascular corrosion casts

We studied the opisthonephric (mesonephric) kidneys of adult male and female Xenopus laevis using scanning electron microscopy (SEM) of vascular corrosion casts and light microscopy of paraplast embedded tissue sections. Both techniques displayed glomeruli from ventral to mid-dorsal regions of the kidneys with single glomeruli located dorsally close beneath the renal capsule. Glomeruli in general were fed by a single afferent arteriole and drained via a single thinner efferent arteriole into peritubular vessels. Light microscopy and SEM of vascular corrosion casts revealed sphincters at the origins of afferent arterioles, which arose closely, spaced from their parent renal arteries. The second source of renal blood supply via renal portal veins varied interindividually in branching patterns with vessels showing up to five branching orders before they became peritubular vessels. Main trunks and their first- and second-order branches revealed clear longish endothelial cell nuclei imprint patterns oriented parallel to the vessels longitudinal axis, a pattern characteristic for arteries. Peritubular vessels had irregular contours and were never seen as clear cylindrical structures. They ran rather parallel, anastomosed with neighbors and changed into renal venules and veins, which finally emptied into the ventrally located posterior caval vein. A third source of blood supply of the peritubular vessels by straight terminal portions of renal arteries (vasa recta) was not found.

Rho kinase inhibitors reduce voltage-dependent Ca2+ channel signaling in aortic and renal microvascular smooth muscle cells

Voltage-dependent L-type Ca²⁺ channels (L-VDCCs) and the RhoA/Rho kinase pathway are two predominant intracellular signaling pathways that regulate renal microvascular reactivity. Traditionally, these two pathways have been thought to act independently; however, recent evidence suggests that these pathways could be convergent. We hypothesized that Rho kinase inhibitors can influence L-VDCC signaling. The effects of Rho kinase inhibitors Y-27632 or RKI-1447 on KCl-induced depolarization or the L-VDCC agonist Bay K8644 were assessed in afferent arterioles using an in vitro blood-perfused rat juxtamedullary nephron preparation. Superfusion of KCl (30-90 mM) led to concentration-dependent vasoconstriction of afferent arterioles. Administration of Y-27632 (1, 5, and 10 µM) or RKI-1447 (0.1, 1, and 10 µM) significantly increased the starting diameter by 16-65%. KCl-induced vasoconstriction was markedly attenuated with 5 and 10 µM Y-27632 and with 10 µM RKI-1447 (P < 0.05 vs. KCl alone). Y-27632 (5 µM) also significantly attenuated Bay K8644-induced vasoconstriction (P < 0.05). Changes in intracellular Ca²⁺ concentration ([Ca²⁺]_i) were estimated by fura-2 fluorescence during KCl-induced depolarization in cultured A7r5 cells and in freshly isolated preglomerular microvascular smooth muscle cells. Administration of 90 mM KCl significantly increased fura-2 fluorescence in both cell types. KCl-mediated elevation of [Ca²⁺]_i in A7r5 cells was suppressed by 1-10 µM Y-27632 (P < 0.05), but 10 µM Y-27632 was required to suppress Ca²⁺ responses in preglomerular microvascular smooth muscle cells. RKI-1447, however, significantly attenuated KCl-mediated elevation of [Ca²⁺]_i. Y-27632 markedly inhibited Bay K8644-induced elevation of [Ca²⁺]_i in both cell types. The results of the present study indicate that the Rho kinase inhibitors Y-27632 and RKI-1447 can partially inhibit L-VDCC function and participate in L-VDCC signaling.

Myoendothelial communication in the renal vasculature and the impact of drugs used clinically to treat hypertension

The renal vasculature has many peculiarities including highly irregular branching. Renal blood flow must sustain adequate perfusion and maintain a high glomerular filtration. Renal autoregulation helps control renal blood flow. The local autoregulatory mechanism, tubuloglomerular feedback, elicits a vasoconstriction that can be found not only in neighboring nephrons but over large areas of the kidney indicating that the renal vasculature supports strong conduction of vascular responses. The basis for conduction is intercellular communication through gap junctions. The endothelium is strongly coupled and serves as a vascular conduction highway leading the signal to the vascular smooth muscle cells through myoendothelial coupling. Extensive intercellular coupling is also found in renin secreting cells where gap junctions seem to tie the cells together to improve control of renin secretion. Lack of coupling leads to dysregulation of renin secretion and hypertension. However, the activity of the renin-angiotensin system also controls gap junction expression in the kidney. Treatment reducing angiotensin II activity, as used in hypertension treatment, can affect expression of renal and vascular gap junction.

The nephron-arterial network and its interactions

Tubuloglomerular feedback and the myogenic mechanism form an ensemble in renal afferent arterioles that regulate single-nephron blood flow and glomerular filtration. Each mechanism generates a self-sustained oscillation, the mechanisms interact, and the oscillations synchronize. The synchronization generates a bimodal electrical signal in the arteriolar wall that propagates retrograde to a vascular node, where it meets similar electrical signals from other nephrons. Each signal carries information about the time-dependent behavior of the regulatory ensemble. The converging signals support synchronization of the nephrons participating in the information exchange, and the synchronization can lead to formation of nephron clusters. We review the experimental evidence and the theoretical implications of these interactions and consider additional interactions that can limit the size of nephron clusters. The architecture of the arterial tree figures prominently in these interactions.

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.