Synchronization in renal microcirculation unveiled with high-resolution blood flow imaging
Abstract
Internephron interaction is fundamental for kidney function. Earlier studies have shown that nephrons signal to each other, synchronize over short distances, and potentially form large synchronized 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 the frequency and phase differences in rat kidney microcirculation under different conditions. The analysis unveiled the spatial and temporal evolution of synchronized blood flow clusters of various sizes, including the formation of large (>90 vessels) and long-lived clusters (>10 periods) locked at the frequency of the tubular glomerular feedback mechanism. Administration of vasoactive agents caused significant changes in the synchronization patterns and, thus, in nephrons’ co-operative dynamics. Specifically, infusion of vasoconstrictor angiotensin II promoted stronger synchronization, while acetylcholine caused complete desynchronization. The results confirm the presence of the local synchronization 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 synchronization in pathologies development.
Editor's evaluation
The new technical advance reported here, high-resolution laser speckle contrast imaging of many microvessels simultaneously and application to the entire rat kidney surface should be of great interest to readers interested in general vascular physiology and especially renal hemodynamics research. The observed large vessel cluster synchronization and phase wave formation phenomena will help further the understanding of renal blood flow regulation in health and disease.
https://doi.org/10.7554/eLife.75284.sa0Introduction
The kidney represents a unique demand-driven, interconnected resource distribution network that is responsible for body homeostasis maintenance over a broad range of conditions, including variations in blood pressure and fluid intake and loss. With blood serving as the resource, single-nephron autoregulation mechanisms provide and regulate the demand. Based on measurements of tubule pressure responses to step changes in arterial pressure (Young and Marsh, 1981), vascular transfer functions (Chon et al., 1993; Chon et al., 2005; Cupples and Loutzenhiser, 1998; Marmarelis et al., 1999; Sakai et al., 1986; Shi et al., 2006), and renal blood flow response to arterial pressure forcing (Holstein-Rathlou et al., 1991), two critical mechanisms in renal pressure autoregulation, the myogenic mechanism and tubuloglomerular feedback (TGF), have emerged as the critical components. Their actions combine to regulate blood flow, serving to maintain the delivery of water and solutes to various regions of the nephron at levels appropriate to their dynamic ranges. The two mechanisms operate at different time scales, generating spontaneous blood flow and pressure oscillations at different frequencies: 5–10 s (0.1–0.2 Hz) for the myogenic response and 30–50 s (0.02–0.033 Hz) for the TGF (Holstein-Rathlou and Marsh, 1989; Yip et al., 1993; Cupples and Loutzenhiser, 1998; Sosnovtseva et al., 2002; Just, 2007).
Nephrons, however, do not operate as stand-alone units. Within a single kidney, all nephrons are linked via the renal vascular tree, which provides connections of different proximities - from few hundred microns for neighboring nephrons separated only by their respective afferent arterioles to nephrons only connected at the level of the renal artery, which plays a role of a single-supply source for all the nephrons. Nephrons nested in such a network are bound to communicate and affect each other to some degree. In addition to interaction through the blood flow and pressure, nephrons were found to communicate via electrical signaling (Marsh et al., 2019; Marsh et al., 2009). Such interactions are proven to play a critical role in kidney function and can lead to complex co-operative dynamics and synchronization between nephrons. Adaptive synchronization across the renal microvascular network would increase the efficiency and dynamic range of autoregulation by engaging more preglomerular resistance (Marsh et al., 2019; Zehra et al., 2021), preventing transmission of high systemic pressure to the glomeruli, which could lead to progressive glomerular and vascular injury. In addition, long-distance synchronization would argue strongly that renal autoregulation is a distributed process that can ensure an optimized oxygenation-perfusion matching and adjust to various internal or environmental conditions (Sosnovtseva et al., 2007; Mitrou et al., 2015). On the other hand, large-scale in-phase synchronization might become a mechanism leading to glomerular injury as it would substantially increase local pressure variation (Postnov et al., 2012).
Micropuncture experiments (Holstein-Rathlou, 1987; Källskog and Marsh, 1990; Yip et al., 1992; Sosnovtseva et al., 2007) showed that neighboring nephrons (originating from a common artery) adjust their TGF-mediated tubular pressure oscillations to attain a synchronized regime. Although these experiments confirmed the existence of synchronization, only pairs or triplets of nephrons could be sampled at any one time. Assessment of co-operative efforts of a larger number of nephrons required a different approach and was made possible with laser speckle contrast imaging (LSCI) technology (Holstein-Rathlou et al., 2011; Mitrou et al., 2015; Scully et al., 2013). In LSCI, media with moving light scattering particles, e.g., red blood cells, are illuminated with a near-infrared laser. The backscattered light, recorded by a camera, forms an interference pattern, which appears more or less blurred depending on the speed of the particles. This pattern is then analysed to obtain qualitative maps of particle velocity and thus a blood flow estimate (Boas and Dunn, 2010; Briers et al., 2013). First applied by Holstein-Rathlou et al. to map TGF oscillations over a large field of view, this method was later used to explore periodic activity in the myogenic frequency band (Scully et al., 2013; Mitrou et al., 2015), analyse spatial correlations in the renal blood flow (Brazhe et al., 2014), and study intrarenal drug distribution (Postnov et al., 2015a, Postnov et al., 2017). Although these results encouraged the large-scale synchronization hypothesis, the method lacked resolution, both spatial and temporal, as well as signal-to-noise ratio, to confirm it convincingly.
In this paper, we further advance renal blood flow imaging methodology and confirm the presence of synchronized clusters spanning multiple nephrons for the first time at the level of individual arterioles and venules. Our LSCI setup and data processing approach allow imaging renal microcirculation at 0.8 μm per pixel spatial resolution and imaging frequency to 160 Hz for 1024 × 1024 pixels. We apply it to study the formation of synchronous blood flow clusters in rat kidney microcirculation under normal conditions and during the administration of vasoactive agents.
Results
Synchronization patterns
An example of a high-resolution blood flow map with segmented microcirculatory vessels is shown in Figure 1A. Figure 1B shows relative blood flow (RBF) dynamics of four vessels outlined in (A) with corresponding colors (red, green, orange, and black). From the corresponding dominant frequency and phase (Figure 1C and F), it can be seen that these vessels form two frequency-locked clusters (red-green and black-orange) that are synchronous most of the time. Moreover, there are long periods when these clusters synchronize, forming a larger cluster, interrupted with an asynchronous interval (from ≈850 to ≈1050 s).

Blood flow data example.
(A) Normalized flow map with microcirculatory vessels. Yellow and blue colors correspond to high and low blood flow, respectively. (B) Relative blood flow (RBF) dynamics in four vessels, which are outlined in (A) with corresponding colors. (C) Dominant frequency in the tubuloglomerular feedback (TGF) band of the blood flow oscillations is shown in (B). (D) Phase value at the dominant frequency. One can see that there are two vessel pairs in which blood flow is synchronized for more than 90% of the time. Synchronization between the pairs is also well observed but breaks down when TGF activity change its frequency at 850–1050 s.
We generalized this approach to the full-field blood flow imaging of the kidney surface. To simplify the visualisation of frequencies and phases, we generate maps where colored circles represent the segmented vessels (outlined in magenta in Figure 1A), with the circle center located in the center of the mass of the vessel. The circle’s color represents the peak TGF frequency observed in the corresponding vessel, while the size of the circles is identical and does not hold any information. Figure 2 shows examples of high (A), moderate (B), and low (C) instantaneous synchronization degree (see Methods, Eq. 1), with S = 0.97, 0.65, and 0.23, respectively. In the first case, flow in almost all of the identified vessels oscillates at the same TGF frequency, forming a cluster of >100 vessels within the field of view. This cluster is likely to be even larger since it can spread outside the view field and in-depth in the kidney cortex. In the case of moderate synchronization degree, several frequency clusters with ≈ 10–60 vessels can be identified, while for the low S, there are multiple groups of two to three vessels displaying the same frequency but no distinct pattern over the field of view. All three regimes were observed in the same animal during control, vasoconstricted (angiotensin II [AngII] infusion), and vasodilated (acetylcholine [ACh] infusion) conditions, respectively. Averaged over the whole observation period and all animals (N = 5, 20 min per condition), synchronization degree has moderate to high () values during AngII infusion, low to moderate () during ACh, and low to high in control (). The average TGF frequencies were measured to be Hz in control, during the AngII infusion, and during the ACh infusion.

Examples of instantaneous frequency patterns.
Each circle location corresponds to a segmented vessel with color-coded frequency. Bottom panels show normalized blood flow in each segmented vessel over 5 min centered around the time corresponding to the frequency pattern snapshot. (A) High synchronization degree - large cluster covers the whole field of view. (B) Moderate synchronization degree - most of the vessels are split between two clusters. (C) Low synchronization degree - no clear clustering pattern, some vessels do not have pronounced tubuloglomerular feedback activity (frequency color-coded as black). All data were captured in the same animal but under different conditions: control (A), angiotensin II infusion (B), and acetylcholine infusion (C). Black-colored circles represent vessels where TGF activity was considered too weak (less than 10% prominence of the activity peak).
To visualize clustering tendencies, we calculate the probability of a randomly chosen vessel at any given moment of time to have the same dominant frequency as X% of the vessels in the field of view. From Figure 3, it can be seen that during AngII infusion, vessels are significantly more likely to be synchronized with 20–50% of the field of view than in control or during ACh infusion. In the latter condition, probability of a vessel being synchronized with less than 15% is significantly higher than during AngII infusion (p<0.05). Such behavior is observed both with no restriction on minimum synchronization duration (A) and when only frequency-locking for three TGF periods and longer is taken into account.

Probability of a randomly chosen vessel displaying the same dominant tubuloglomerular feedback (TGF) frequency as X% of the vessels in the field of view.
(A) Any duration of frequency matching is considered, and (B) dominant frequency should match for at least three TGF periods. It can be seen that infusion of angiotensin II significantly (p<0.05) increases the prevalence of clusters covering 20–50% of the vessels - chances that a randomly chosen vessel belongs to such a cluster at any given moment are ≈50%, while for control and acetylcholine infusion they are ≈20%. N = 5 animals were used to create these graphs. Paired t-test was used to produce p-values. p-values smaller than 0.05 are considered to be significant and marked with ‘*’.
Phase waves and spatial localization
Another distinct feature that we have observed is that the phase of oscillations within a cluster is space-dependent - it is mostly the same for the closely positioned vessels and gradually changes with distance. Figure 4 illustrates how phase within clusters evolves in space and time, forming the phase waves. Note that both (A) and (B) panels are from the same animal as was shown in Figure 3 in control and vasoconstricted conditions, respectively, and that only vessels belonging to the largest cluster are shown.

Phase waves.
(A) and (B) show the spatio-temporal evolution of tubuloglomerular feedback activity phase in vessels that belong to the largest cluster observed during control and angiotensin II infusion. Both observations are from the same animal as shown in Figure 2A, B. Arrows indicate the wave direction.
Change of the wave direction by ≈90° in the same animal with constant vascular structure suggests different synchronization centers, where phase waves originate. Propagation speed is and 0.30 mm/s for (A) and (B), respectively, and spatial period of the phase waves is ≈ 8.85 and 3.95 mm (respective phase difference 0.71 and 1.59 rad was observed over 1 mm distance). The possibility of such phase waves was previously hinted at in our earlier LSCI study, where spatial correlations in renal blood flow were analyzed (Brazhe et al., 2014). See Figure 4—videos 1–3 for more examples of phase waves dynamics.
Phase difference distribution (Figure 5A), calculated over all of the frequency matching vessels in all animals, shows the prevalence of in-phase synchronization (, rad). This result is in good agreement with experimental micropuncture observations, where it is explained by the presence of fast electrical coupling acting over short distances. In all conditions, the larger phase differences are less prevalent, with the anti-phase synchronization (, rad) observed over just 1% of the time in control and during ACh infusion. However, AngII infusion increases this number to 5%, showing a statistically significant difference with other conditions. Phase difference grows with distance, as can be seen from Figure 5B, reflecting the presence of phase waves. While change is relatively small in the control and vasodilated conditions (≈0.3 and 0.25 rad/mm), it is strongly enhanced in the vasoconstricted condition, reaching, on average, ≈1 rad over 1 mm distance. As we showed with mathematical modeling (Postnov et al., 2012), such difference can be explained by strengthened hemodynamic coupling, which the increased vascular tone should cause.

Localization of phase and synchronization in space.
(A) Phase differences prevalence for synchronized vessels. (B) Phase differences distribution over the distance between vessels. (C) Synchronization duration normalized by the total observation time for different distances between vessels. N = 5 animals were used to create these graphs. Paired t-test was used to produce p-values. p-values smaller than 0.05 are considered to be significant and marked with ‘*.’ Values between 0.05 and 0.1 are shown as p<0.1, where relevant, to highlight a trend in the data.
Since synchronization in renal blood flow is extended far beyond a pair of nephrons and is unlike synchronization in relatively homogeneous media, one would expect some space localization due to the topological features of the vascular network. Figure 5C illustrates how average synchronization duration changes with distance. It is clear that synchronization is stronger in the vasoconstricted condition, with average duration reaching ≈40±13% of the observation time for neighboring vessels and gradually reducing to 24±6% for vessels located at more than 1 mm distance. In control, synchronization duration varies greatly, but on average, it also reduces with distance, although at a slower rate than during AngII infusion—from ≈30% to ≈21%. During the ACh infusion the synchronization duration is ≈15-13%, with only 2% reduction over 1mm of distance. Higher synchronization duration in control and during AngII infusion compared with ACh infusion suggests a long-distance nephron-to-nephron communication or a common driving force.
Discussion
In this study, we have designed a methodology for high-resolution blood flow imaging in renal microcirculation and applied it to study the synchronization of TGF oscillations in control, vasoconstricted (AngII infusion), and vasodilated (ACh infusion) conditions. Our data confirm that blood flow in renal microcirculation tends to demonstrate clustered, frequency-locked activity, with the clustering size and tendencies changing depending on the animal condition. Synchronization tends to be stronger, and cluster size is larger during AngII infusion, with synchronization degree and average synchronization duration ranging from ≈40 ± 13 to of the observed time depending on the distance between vessels. During the ACh infusion, on the contrary, synchronization seems to be disrupted, with () and average synchronization duration of the observed time. Finally, in the normotensive condition, we observed mixed behavior with highly variable synchronization degree and duration ranging from ≈30 ± 27 to . A possible explanation for the stronger synchronization during the AngII infusion will be a stronger hemodynamic coupling due to increased vascular resistance. It is also supported by the increased number of vessels synchronized in anti-phase during the infusion, which, as we predicted using mathematical modeling (Postnov et al., 2012), is a natural consequence of stronger hemodynamic coupling.
We also observed phase waves that travel over the TGF-frequency clustered vessels (Figures 4 and 5B). Interestingly, the direction of the phase waves appears to be not predetermined by the vascular structure - depending on the flow dynamics and other, yet unknown factors, it can change even within the same animal (an example is shown in Figure 4). Such behavior might be related to diffusive interaction between nephrons or the formation of synchronization centers (pacemakers) at different locations. The presence of phase waves supports the deterministic nature of synchronization rather than random entrainment at the same dominant TGF frequency. Similar effects are well studied in excitable media with diffusion interaction mechanisms such as brain (Dahlem et al., 2010; Shibata and Bures, 1972) and heart (Alonso et al., 2016) tissues. Such similarity might suggest the presence of a diffusive mechanism in the inter-nephron interaction or a long-distance fast electrical signaling, e.g., via conducted vasoreactivity. However, a detailed exploration of the topic would require additional experiments combining multi-scale blood flow and structural imaging methods.
In our experiments, we imaged 1.5 × 1.5 mm2 field of view in five animals with individual segmented vessels on average, each of which is 1o–30 µm in diameter. While it provides an estimate of synchronization in the nephrons activity, direct translation from individual segmented vessels to nephrons is challenging and requires further exploration. Factors that are critical to consider are: (i) penetration depth of LSCI when applied to renal imaging, (ii) type of the vessels in the field of view, and (iii) topology of the nephro-vascular network. While in theory, LSCI can collect the blood flow signal from as deep as 300–400 μm, in practice, visually resolvable signal typically comes from top 50 to 150 μm of the vascular structure (Davis et al., 2014; Zheng et al., 2021). Considering high vascular density close to the renal surface, it would mean that LSCI is likely limited to imaging vessels originating from ≈10000 nephrons in outer 30% of rat renal cortex (Letts et al., 2017), which would result in ≈40 nephrons in the 1.5 × 1.5 µm field of view. The larger number of segmented vessels can be explained by their mixed type - afferent and efferent arterioles as well as venules are likely to be segmented. Distinguishing vessel types in LSCI images will require further exploration and registration with high-resolution structural imaging. It, however, does not mean that the observed clusters were limited to, at the most, 40 nephrons. When considering renal vascular topology, it is to be expected that within the 1.5 × 1.5 µm field of view, we observe arterioles that arise from different non-terminal arteries (Marsh et al., 2017). Depending on the branching order, each of such arteries can branch into 10 to several hundreds of nephrons, but only a small number of these nephrons will have arterioles reaching close enough to the surface to be segmented from LSCI images. Thus, when a synchronous cluster is observed with LSCI, it is likely to extend several branching orders in depth and reach the size of hundreds and even thousands of nephrons.
While the exact role of inter-nephron communication, co-operative dynamics, and synchronization in kidney-related pathology development is still unclear and requires further exploration, it is evidently altered by the blood pressure and vascular tone. Strong local coupling and in-phase synchronization, while being not evident at the renal artery level (Postnov et al., 2015a), are likely to increase pressure variation at the level of afferent arterioles (Postnov et al., 2012; Postnov et al., 2016a), thus increasing chances of local damage and aggravating pathological condition.
Materials and methods
Animal preparation
Request a detailed protocolAll experimental protocols were approved by the Danish National Animal Experiments Inspectorate (License 2015-15-0201-00463) and were conducted according to the American Physiological Society guidelines. Male Sprague Dawley rats (Taconic, Denmark) with average weight ≈290 g (n = 5) were used. Before starting surgical procedures, animals were anesthetized in a chamber with 8% sevoflurane. During the surgery, sevoflurane concentration was reduced to a final concentration of ≈2%. Two catheters were inserted in the right jugular vein to allow continuous systemic infusion of drugs and saline. Another catheter was inserted in the carotid artery to measure mean arterial pressure with a pressure transducer (Statham P23-dB, Gould, Oxnard, CA). Then tracheotomy was performed, after which the rat was placed on a servo-controlled heating table maintaining body temperature at 37°C and connected to a mechanical animal ventilator (60 breaths/min; 8 ml/kg bodyweight). To avoid secondary heartbeat and breathing artifacts, Nimbex (muscle relaxant, Sigma) was administered in a concentration of 0.85 mg/ml, first as a bolus injection of 0.5 ml, followed by a continuous intravenous infusion at a rate of 20 μl/min. The left kidney was then exposed, and the left ureter was catheterized to ensure free urine flow. To reduce motion artifacts and avoid drying the kidney surface during the experiment, we placed the kidney in a plastic fixation holder, covered it with warm agarose solution (1% Agarose, Sigma, 99% saline), and put a thin (0.1 mm) cover glass on top of the kidney. Metal thread (40 µm in diameter) was bent in a ‘U’ shape and positioned on top of the cover glass at the flattest location of the kidney surface, marking the region of interest for imaging procedures. Following the surgical procedures, the animal was left to stabilize for 20 min. Experiments were continued only if the mean arterial pressure remained within 100–120 mm Hg during the control period. At the end of the experiment, animals were euthanized by overdose of sevoflurane, followed by cervical dislocation.
LSCI data acquisition
Request a detailed protocolSchematics of the imaging setup are shown in Figure 6A, B. To assess microcirculation in the kidney cortex, we built a high resolution LSCI setup. A single-mode fiber-coupled laser diode (785 nm, LP785-SF100, Thorlabs, USA) controlled with a laser driver and temperature controller (CLD1011LP, Thorlabs, USA) was used to deliver coherent light onto kidney surface with a power density of approximately 10 mW/cm2, providing an optimal signal-to-noise ratio (Postnov et al., 2015b). Backscattered light was collected by a zoom imaging lens (VZM 1000i, Edmund Optics) at ×5 magnification and recorded with a CMOS camera (Basler acA2000-165umNIR, 2048 × 1088 pixels, 5.5 μm pixel size). Exposure time was set to 5 ms, which provides optimal sensitivity to slow flows, such as in capillary perfused tissue (Yuan et al., 2005). A subset of 1024 × 1024 pixels was used for recording, from a 1.5 × 1.5 mm region of the renal surface. In addition, a linear polarizing filter was placed in front of the objective to reduce artifacts from reflected light. While the imaging setup in this configuration allows imaging resolution 1 of ≈0.8 × 0.8 µm2/px at the frame-rate of 160 frames per second (fps), in our experiments, we found it optimal, in terms of the field-of-view and data storage, to acquire images at 50 fps and resolution of ≈1.5 × 1.5 µm2/px. The system provides an image resolution of 2.25 µm2/px compared with previously published studies on renal synchronization, which used the system with the image resolution of 45–400 µm2/px (Holstein-Rathlou et al., 2011; Scully et al., 2013; Mitrou et al., 2015). Similarly, the system has a higher temporal resolution (50 fps instead of 25) and utilizes a better light source (Postnov et al., 2015b). Combined and considering noise reduction achieved by spatial and temporal averaging, these factors improve the expected signal-to-noise ratio per area unit by at least 30 times. The examples of data collected with the previously used system and the system presented in this study can be seen in Figure 1—figure supplement 1. Further information on the role of high-resolution imaging for identifying vessel-specific oscillations frequency and phase can be found in our recent study (Lee et al., 2022).
Imaging was performed in a non-stimulated state (control) and following administration of AngII and ACh, respectively. After collecting 20 min (2400 frames) of the baseline data, we initiated continuous administration of AngII at a concentration of 4 ng/ml and an infusion rate of 20 μl/min to cause systemic vasoconstriction. The infusion lasted for 30 min, out of which the first 10 min were allocated for blood flow to stabilize, and the following 20 min were recorded. 15 min after completion of the AngII infusion, we infused ACh at a concentration of 0.0375 mg/ml and a rate of 20 μl/min, first permitting 10 min for blood flow to stabilize, and then recorded subsequent 20 min. Across all experiments, average arterial pressure was 112 ± 2, 127 ± 6, and 100 ± 4 during the control, vasoconstrited (AngII), and vasodilated (ACh) conditions, respectively. Averaged across all vessels blood flow index has decreased by 17.3±2.3% in the vasoconstricted condition and increased by 27.7±13.3% in the vasodilated condition. For the details on blood flow index measurements across the conditions see Figure 1—figure supplement 1.
Data analysis
Image registration
Request a detailed protocolTo allow high-resolution LSCI of the renal microcirculation, we needed to reduce motion artifacts, as any lateral motion larger than 5–10 µm will prevent accurate estimation of the blood flow and further segmentation of microcirculatory vessels. Unlike brain imaging, when working with the kidney, there is no bone tissue that can be fixed to reduce respiratory motion, and applying even slight pressure on the kidney might result in abnormal blood flow due to the occlusion of the small vessels on its surface. At the same time, raw laser speckle images, or contrast images without temporal averaging, are not suitable for automated registration due to the absence of clear intensity landmarks (Postnov et al., 2016b). To resolve this issue, we placed a ‘U’ shaped metal marker on the cover glass, which moves along with the kidney, as described above. As the first step of analyzing the data, the marker is segmented in all frames via thresholding and then used to estimate the translation type geometrical transformation required to register images. Estimated geometrical transformation is then applied to the raw laser speckle images prior to performing the contrast analysis.
Contrast analysis
Request a detailed protocolRegistered laser speckle images were processed to calculate temporal contrast , where and are SD and mean of pixel intensity over 25 frames (Postnov et al., 2016b). Contrast values were then converted to the blood flow index as , which are then used in the ensuing analysis.
Vessels segmentation
Request a detailed protocolTo segment individual microcirculatory vessels, we calculated averaged in time blood flow index (BFI) images and applied adaptive thresholding (MATLAB) to them. Automated segmentation was followed by manual clean-up, where we removed artifacts and occasional large surface vessels. We then calculated blood flow dynamics for each segmented individual microcirculatory vessel by averaging BFI values in pixels belonging to this vessel.
Synchronization analysis
Request a detailed protocolTo study synchronization patterns between microcirculatory vessels and, thus, obtain insight into inter-nephron communication, we apply continuous wavelet transform analysis (Morse wavelet, MATLAB) to segmented vessels’ BFI. We identified the frequency and phase of dominant periodic activity in the 0.015–0.05 Hz frequency band associated with the TGF mechanism. Vessels with less than 10% prominence of the activity peak were discarded and not used for synchronization analysis. In this study, we consider blood flow in different segmented vessels to be synchronized whenever their dominant frequencies match. To quantify blood flow synchronization over the field of view, we analyzed phase differences between synchronized vessels, average synchronization duration, and its dependency on the distance between vessels, and the probability of the vessel’s blood flow to be synchronized with N% of the vessels in the field of view. To provide a ‘single-value’ characterization of the synchronization at a given moment of time, we also introduced the synchronization degree parameter :
where is a number of frequency matching pairs of vessels, is a total number of observed segmented vessels, and is the time. represents the relation of the observed number of frequency matching pairs to their maximum possible number. Thus corresponds to all segmented blood vessels having the same dominant frequency, while to all segmented blood vessels having a different dominant frequencies. However, it is important to notice that is impossible to reach due to the discrete nature of the measured data and the analysis. In our case, the 0.015–0.05 Hz range is split into 18 fixed values, so that if there were more than 18 vessels, it became unavoidable for some of them to have an identical dominant frequencies. In practice, for 100 vessels with randomly chosen dominant frequencies, the minimum observed would be ≈0.25±0.06, which can be confirmed with a simple computational experiment.
Statistical analysis
Request a detailed protocolPaired t-test was applied to compare results between control, AngII infusion, and ACh infusion. p-values greater than 0.05 are reported as not significant. Results were expressed as mean ± SD unless indicated otherwise.
Data availability
The data underlying this article are available at public data repository (Dryad): https://doi.org/10.5061/dryad.g79cnp5r2.
-
Dryad Digital RepositoryData from: Synchronization in renal microcirculation unveiled with high-resolution blood flow imaging.https://doi.org/10.5061/dryad.g79cnp5r2
References
-
Nonlinear physics of electrical wave propagation in the heart: a review Reports on Progress in PhysicsPhysical Society (Great Britain) 79:096601.https://doi.org/10.1088/0034-4885/79/9/096601
-
Laser speckle contrast imaging in biomedical opticsJournal of Biomedical Optics 15:011109.https://doi.org/10.1117/1.3285504
-
Laser speckle contrast imaging: theoretical and practical limitationsJournal of Biomedical Optics 18:066018.https://doi.org/10.1117/1.JBO.18.6.066018
-
On the efficacy of linear system analysis of renal autoregulation in ratsIEEE Transactions on Bio-Medical Engineering 40:8–20.https://doi.org/10.1109/10.204766
-
Interactions of TGF-dependent and myogenic oscillations in tubular pressureAmerican Journal of Physiology. Renal Physiology 288:F298–F307.https://doi.org/10.1152/ajprenal.00164.2004
-
Dynamic autoregulation in the in vitro perfused hydronephrotic rat kidneyAmerican Journal of Physiology-Renal Physiology 275:F126–F130.https://doi.org/10.1152/ajprenal.1998.275.1.F126
-
Imaging depth and multiple scattering in laser speckle contrast imagingJournal of Biomedical Optics 19:086001.https://doi.org/10.1117/1.JBO.19.8.086001
-
Oscillations of tubular pressure, flow, and distal chloride concentration in ratsThe American Journal of Physiology 256:F1007–F1014.https://doi.org/10.1152/ajprenal.1989.256.6.F1007
-
Tubuloglomerular feedback dynamics and renal blood flow autoregulation in ratsAmerican Journal of Physiology-Renal Physiology 260:F53–F68.https://doi.org/10.1152/ajprenal.1991.260.1.F53
-
Nephron blood flow dynamics measured by laser speckle contrast imagingAmerican Journal of Physiology-Renal Physiology 300:F319–F329.https://doi.org/10.1152/ajprenal.00417.2010
-
Mechanisms of renal blood flow autoregulation: dynamics and contributionsAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 292:R1–R17.https://doi.org/10.1152/ajpregu.00332.2006
-
TGF-initiated vascular interactions between adjacent nephrons in the rat kidneyThe American Journal of Physiology 259:F60–F64.https://doi.org/10.1152/ajprenal.1990.259.1.F60
-
Multi-scale laser speckle contrast imaging of microcirculatory vasoreactivityBiomedical Optics Express 13:2312–2322.https://doi.org/10.1364/BOE.451014
-
Nephron morphometry in mice and rats using tomographic microscopyAmerican Journal of Physiology. Renal Physiology 312:F210–F229.https://doi.org/10.1152/ajprenal.00207.2016
-
Nonlinear analysis of renal autoregulation in rats using principal dynamic modesAnnals of Biomedical Engineering 27:23–31.https://doi.org/10.1114/1.222
-
Electrotonic vascular signal conduction and nephron synchronizationAmerican Journal of Physiology. Renal Physiology 296:F751–F761.https://doi.org/10.1152/ajprenal.90669.2008
-
Architecture of the rat nephron-arterial network: analysis with micro-computed tomographyAmerican Journal of Physiology. Renal Physiology 313:F351–F360.https://doi.org/10.1152/ajprenal.00092.2017
-
The nephron-arterial network and its interactionsAmerican Journal of Physiology. Renal Physiology 316:F769–F784.https://doi.org/10.1152/ajprenal.00484.2018
-
Laser speckle contrast imaging reveals large-scale synchronization of cortical autoregulation dynamics influenced by nitric oxideAmerican Journal of Physiology. Renal Physiology 308:661–670.https://doi.org/10.1152/ajprenal.00022.2014
-
Dynamics of nephron-vascular networkBulletin of Mathematical Biology 74:2820–2841.https://doi.org/10.1007/s11538-012-9781-6
-
Laser speckle imaging of intra organ drug distributionBiomedical Optics Express 6:5055–5062.https://doi.org/10.1364/BOE.6.005055
-
Improved detectability of microcirculatory dynamics by laser speckle flowmetryJournal of Biophotonics 8:790–794.https://doi.org/10.1002/jbio.201500152
-
Modeling of Kidney Hemodynamics: Probability-Based Topology of an Arterial NetworkPLOS Computational Biology 12:e1004922.https://doi.org/10.1371/journal.pcbi.1004922
-
Estimation of vessel diameter and blood flow dynamics from laser speckle imagesBiomedical Optics Express 7:2759–2768.https://doi.org/10.1364/BOE.7.002759
-
A simple method to ensure homogeneous drug distribution during intrarenal infusionAmerican Journal of Physiology. Renal Physiology 312:F543–F548.https://doi.org/10.1152/ajprenal.00417.2016
-
Frequency domain analysis of renal autoregulation in the ratThe American Journal of Physiology 250:364–373.https://doi.org/10.1152/ajprenal.1986.250.2.F364
-
Detecting physiological systems with laser speckle perfusion imaging of the renal cortexAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 304:929–939.https://doi.org/10.1152/ajpregu.00002.2013
-
Tubuloglomerular feedback-dependent modulation of renal myogenic autoregulation by nitric oxideAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 290:R982–R991.https://doi.org/10.1152/ajpregu.00346.2005
-
Reverberation of cortical spreading depression along closed-loop pathways in rat cerebral cortexJournal of Neurophysiology 35:381–388.https://doi.org/10.1152/jn.1972.35.3.381
-
Bimodal oscillations in nephron autoregulationPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 66:061909.https://doi.org/10.1103/PhysRevE.66.061909
-
Synchronization among mechanisms of renal autoregulation is reduced in hypertensive ratsAmerican Journal of Physiology. Renal Physiology 293:1545–1555.https://doi.org/10.1152/ajprenal.00054.2007
-
Dynamics of TGF-initiated nephron-nephron interactions in normotensive rats and SHRThe American Journal of Physiology 262:980–988.https://doi.org/10.1152/ajprenal.1992.262.6.F980
-
Mechanisms of temporal variation in single-nephron blood flow in ratsThe American Journal of Physiology 264:427–434.https://doi.org/10.1152/ajprenal.1993.264.3.F427
-
Pulse wave propagation in rat renal tubules: implications for GFR autoregulationThe American Journal of Physiology 240:446–458.https://doi.org/10.1152/ajprenal.1981.240.5.F446
-
Tubuloglomerular Feedback Synchronization in Nephrovascular NetworksJournal of the American Society of Nephrology 32:1293–1304.https://doi.org/10.1681/ASN.2020040423
-
Depth resolution in multifocus laser speckle contrast imagingOptics Letters 46:5059–5062.https://doi.org/10.1364/OL.436334
Decision letter
-
Martin R PollakSenior and Reviewing Editor; Harvard Medical School, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your work entitled "Synchronization in renal microcirculation unveiled with high-resolution blood flow imaging" for further consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and overseen by Martin Pollak as Senior Editor. The reviewers have opted to remain anonymous.
There are some issues that need to be addressed. In particular, we would like to see some evidence regarding how this method compares with previous methodologies (see Reviewer 1). We would also like to see the other specific points raised by the reviewers addressed to our satisfaction.
Reviewer #1 (Recommendations for the authors):
The key observations of the study are the synchronization of oscillations in nephron blood flow oscillation frequency at frequencies previously demonstrated to be governed by tubuloglomerular feedback (TGF). Using laser speckle microscopy, the authors observe similar frequency of blood flow oscillation across many nephrons, spanning an area over 1 mm2. The authors then characterize this synchronization of oscillation frequency in response to systemic administration of angiotensin II and acetylcholine to induce systemic vasoconstriction and vasodilation and observe that angiotensin II increases differences in phase, while acetylcholine decreases the duration of synchronization. The authors claim that this manuscript presents a novel laser speckle microscopy system to enable characterization of blood flow across a larger area of the kidney, yet the system is not characterized for resolution or field of analysis nor is it quantitatively compared to previous systems. Therefore, the merits of the manuscript rest on the characterization of nephron synchronization, which is largely correlative and lacking in data presentation and controls.
– In the abstract, the authors claim that this manuscript introduces "an approach for high-resolution laser speckle imaging." As motivation for this work, in the introduction, the authors state, "The methods lacked resolution, both spatial and temporal, as well as signal-to-noise ratio, to confirm it convincingly. In this paper we further advance renal blood flow imaging methodology…" However, the authors do not characterize their new system, in terms of resolution (spatial and temporal) and signal-to-noise ratio. Therefore, the advances claimed cannot be evaluated against previous methods.
– In Figure 2, the length scale of subpanels A-C is not indicated. It's mentioned in the discussion to be 1.5 x 1.5 mm, but it is not indicated on the figure or in the caption. How are the size of the dots determined? Also, given the thresholding and image processing required to generate this data set, it is difficult to understand how the dots presented in the figure relate back to the structure of the kidney. Raw data and images with the nephrons indicated (as in Figure 1) should be provided.
– Similarly, it is difficult to understand how representative the data provided are. While the analysis has been performed on multiple rats (as indicated in the data of Figure 3), there is no information on the variation from measurement to measurement vs. animal to animal, regarding segmentation, synchronization, or effects of administered compound.
– In all conditions, the authors administer the vasoconstriction agent followed by vasodilation. Does this ordering affect the synchronization?
– In the discussion, the authors claim that TGF-frequency clustered waves propagate in a manner where "the direction of the phase is not predetermined by the vascular structure – depending on the flow dynamics and other, yet unknown factors…" What data justify the claim regarding structure and flow dynamics?
Reviewer #2 (Recommendations for the authors):
In this study, Postnov et al., are reporting the application of a new technical advance, further improvement of their recently developed imaging technology using LSCI. In the present work, they performed high-resolution LSCI imaging measurement and analysis of renal blood flow dynamics in multiple individual small renal vessels simultaneously in the entire surface of intact rat kidneys. They found synchronized clusters of blood vessels spanning multiple (>40) nephrons, increased or decreased level of synchronization in response to AngII or Ach, respectively, and the development of space-dependent phase waves. From these results the authors concluded that strengthened hemodynamic coupling due to increased vascular tone, alternating synchronization (pacemaker) centers likely played important roles in the observed phenomena.
Strengths:
The paper is well-written, and the study appears to be well performed. The results are clearly presented and support the conclusions. The presently applied imaging technical advance will be helpful to better understand the detailed, nephron-to-nephron level functional features and role of renal hemodynamics in physiological and disease conditions. However, there are a few points that could be more clarified.
Weaknesses:
1. It would be more convincing if vasoconstricted (Ang II infusion) and vasodilated conditions (Ach infusion) would be validated experimentally by quantitative vessel diameter measurements rather than simply assumed based on the expected vasoactivity of the applied drugs. This is especially important due to the rather low doses applied.
2. Since the LSCI signal strength depends on healthy blood flow rates, vasoconstrictor states may limit the sensitivity of this technique to detect the detailed, single vessel hemodynamic oscillations and the large cluster synchronization features described here. It would be interesting to know the level of total renal blood flow alterations in the presently applied vasoactive treatment conditions.
3. As it is understandable for a technical advance type of report, the study is rather descriptive, and the study design is lacking mechanistic examination of the observed synchronization and phase waves phenomena. Although the presence of fast electric or diffusive coupling that generates long-distance nephron-to-nephron communication was discussed. The applied vasoactive compounds, AngII and Ach have numerous extra-renal effects, for example on nerve activity which may be involved in cluster synchronization. It would be interesting to examine how the large cluster synchronization and phase waves are altered in kidneys one week after denervation. Since a vasoconstrictor state was associated with strengthened synchronization, one is wondering if stronger synchronization could be a sign of increased sympathetic nerve activity.
4. The physiological significance of nephron cluster synchronizations should be clearly stated.
In addition to the many strengths detailed above, I have a few suggestions for further improvement.
1. Vasoconstricted (Ang II infusion) and vasodilated conditions (Ach infusion) need to be validated experimentally by quantitative vessel diameter measurements. Vessel dimensions could be measured using the LSCI images or other higher resolution imaging modalities that could be simultaneously applied.
2. It would be important to show if in the AngII condition the reduced blood flow and reduced sensitivity of the technique could contribute, at least in part, to the observed large clustering. Perhaps including an original, continuous tracing of RBF oscillations before and after AngII treatment would alleviate this concern.
3. Renal denervation could address mechanistic details of the observed large clustering and phase wave phenomena.
4. One or two additional sentences in the introduction should state the physiological significance of nephron synchronization.
https://doi.org/10.7554/eLife.75284.sa1Author response
Reviewer #1 (Recommendations for the authors):
The key observations of the study are the synchronization of oscillations in nephron blood flow oscillation frequency at frequencies previously demonstrated to be governed by tubuloglomerular feedback (TGF). Using laser speckle microscopy, the authors observe similar frequency of blood flow oscillation across many nephrons, spanning an area over 1 mm2. The authors then characterize this synchronization of oscillation frequency in response to systemic administration of angiotensin II and acetylcholine to induce systemic vasoconstriction and vasodilation and observe that angiotensin II increases differences in phase, while acetylcholine decreases the duration of synchronization. The authors claim that this manuscript presents a novel laser speckle microscopy system to enable characterization of blood flow across a larger area of the kidney, yet the system is not characterized for resolution or field of analysis nor is it quantitatively compared to previous systems. Therefore, the merits of the manuscript rest on the characterization of nephron synchronization, which is largely correlative and lacking in data presentation and controls.
– In the abstract, the authors claim that this manuscript introduces "an approach for high-resolution laser speckle imaging." As motivation for this work, in the introduction, the authors state, "The methods lacked resolution, both spatial and temporal, as well as signal-to-noise ratio, to confirm it convincingly. In this paper we further advance renal blood flow imaging methodology…" However, the authors do not characterize their new system, in terms of resolution (spatial and temporal) and signal-to-noise ratio. Therefore, the advances claimed cannot be evaluated against previous methods.
We have added clarification in the text and a supplementary figure (Figure 1—figure supplement 1) that compares measurements taken with our system versus measurements taken with the system used for the renal LSCI before. Further information on the importance of high-resolution imaging and vessel masking approach can be found in our recent study (see Lee at all 2022).
Changes:
Text added (Lines 253-263): “The system provides an image resolution of 2.25µ2 compared to previously published studies on renal synchronization, which used the system with the image resolution of 45 – 400µ2 {holstein2011nephron,scully2013detecting,mitrou2015laser}. Similarly, the system has a higher temporal resolution (50 frames per second instead of 25) and utilizes a better light source {postnov2015improved}. Combined and considering noise reduction achieved by spatial and temporal averaging, these factors improve the expected signal to noise ratio per area unit by at least 30 times. The examples of data collected with the previously used system and the system presented in this study can be seen in Supplementary Figure 1. Further information on the role of high-resolution imaging for identifying vessel-specific oscillations frequency and phase can be found in our recent study{lee2022multi}.”
– In Figure 2, the length scale of subpanels A-C is not indicated. It's mentioned in the discussion to be 1.5 x 1.5 mm, but it is not indicated on the figure or in the caption. How are the size of the dots determined? Also, given the thresholding and image processing required to generate this data set, it is difficult to understand how the dots presented in the figure relate back to the structure of the kidney. Raw data and images with the nephrons indicated (as in Figure 1) should be provided.
– Similarly, it is difficult to understand how representative the data provided are. While the analysis has been performed on multiple rats (as indicated in the data of Figure 3), there is no information on the variation from measurement to measurement vs. animal to animal, regarding segmentation, synchronization, or effects of administered compound.
Colour-coded circles in Figure 2 are a way to provide clear visualization of the dynamics identified in separate vessels. The size of the circles is identical and holds no information, while the location is defined by the centre of mass of corresponding vessels. All the panels presented in Figures 2 and 4 serve to demonstrate examples of observed dynamics and correspond to different recordings from the same animal shown in Figure 1 A. Thus, we already show the data indicating the positions as asked by the reviewer.
The information on variation between animals is provided in Figures 3 and 5 and is discussed in detail in the text. The average and deviation in figures 3 and 5 are calculated over all 5 animals.
We have clarified the text and scale bars to all relevant figures.
Changes:
Scale bars added to figures 2 and 4.
Text added (Lines 95-100): “To simplify the visualisation of frequencies and phases, we generate maps where coloured circles represent the segmented vessels (outlined in magenta in Figure~{Figure 1}(A)), with the circle centre located in the centre of the mass of the vessel. The circle's colour represents the peak TGF frequency observed in the corresponding vessel, while the size of the circles is identical and does not hold any information.”
– In all conditions, the authors administer the vasoconstriction agent followed by vasodilation. Does this ordering affect the synchronization?
We should highlight that the study aims not to see the effect of ACh after infusion of AngII, but to observe synchronization in three different conditions (normal, vasoconstricted, vasodilated). The infusions order protocol was designed for this purpose considering the washout time of the drugs. AngII action is expected to be entirely gone (and receptors recover) within 25 minutes of the resting time after the infusion is over. Changing the order should not change the result as long as the washout times taken into account.
– In the discussion, the authors claim that TGF-frequency clustered waves propagate in a manner where "the direction of the phase is not predetermined by the vascular structure – depending on the flow dynamics and other, yet unknown factors…" What data justify the claim regarding structure and flow dynamics?
It is justified by the direction of the phase waves changing in the same animal depending on time and condition (meaning that it changes within the same vascular structure). Figure 2 shows an example: in the top panels, the phase wave travels from bottom-left to top-right, while in the bottom panels, it travels from top left to bottom right. As we write in the discussion, however, confirming and understanding the mechanism of the phase waves requires further analysis and experiments that would combine multi-scale blood flow and structural imaging.
Changes:
Text changed (Line 177): "is not predetermined" to "appears to be not predetermined".
Text added (Lines 186-187): "However, a detailed exploration of the topic would require additional experiments combining multi-scale blood flow and structural imaging methods."
Reviewer #2 (Recommendations for the authors):
In this study, Postnov et al., are reporting the application of a new technical advance, further improvement of their recently developed imaging technology using LSCI. In the present work, they performed high-resolution LSCI imaging measurement and analysis of renal blood flow dynamics in multiple individual small renal vessels simultaneously in the entire surface of intact rat kidneys. They found synchronized clusters of blood vessels spanning multiple (>40) nephrons, increased or decreased level of synchronization in response to AngII or Ach, respectively, and the development of space-dependent phase waves. From these results the authors concluded that strengthened hemodynamic coupling due to increased vascular tone, alternating synchronization (pacemaker) centers likely played important roles in the observed phenomena.
Strengths:
The paper is well-written, and the study appears to be well performed. The results are clearly presented and support the conclusions. The presently applied imaging technical advance will be helpful to better understand the detailed, nephron-to-nephron level functional features and role of renal hemodynamics in physiological and disease conditions. However, there are a few points that could be more clarified.
Weaknesses:
1. It would be more convincing if vasoconstricted (Ang II infusion) and vasodilated conditions (Ach infusion) would be validated experimentally by quantitative vessel diameter measurements rather than simply assumed based on the expected vasoactivity of the applied drugs. This is especially important due to the rather low doses applied.
We do confirm the effect of the vasoactive drugs by measuring average arterial blood pressure via carotid artery catheter. As it is stated at the end of the Methods section 4.2 the average arterial pressure was observed to be 112±2 in control, 127±6 during vasoconstriction and 100±4 during vasodilation which agrees well with our expectations. The effect is further confirmed by changes in blood flow index measured with LSCI.
Regarding the diameter measurements – reliable measurements of star-vessels diameter with LSCI do not appear to be feasible due to specifics of renal vascular geometry or unless the field of view is further reduced which will make synchronization studies impossible. It could be done if LSCI is combined with multi-photon microscopy, which is out of the scope for the present paper but is part of studies that we plan in future.
2. Since the LSCI signal strength depends on healthy blood flow rates, vasoconstrictor states may limit the sensitivity of this technique to detect the detailed, single vessel hemodynamic oscillations and the large cluster synchronization features described here. It would be interesting to know the level of total renal blood flow alterations in the presently applied vasoactive treatment conditions.
We thank reviewer for a suggestion and have added supplementary figure S2, which shows average blood flow index in all three conditions with and without normalization. From it you can see that blood flow index is reduced by 20% during AngII infusion and increased by 30% during Ach infusion.
We should mention, however that while LSCI sensitivity depends on the relation between particles speed and the exposure time, the effect can be neglected in our study. We use exposure time = 5ms, which is optimal for slow flows (e.g. capillary perfused tissue) and matches our application well. Even with vasoactive agents induced the blood flow index change is +-30%, while noticeable loss of sensitivity is expected for flows which are much faster or slower (e.g. 10 times) (see Yuan et al., 2005 [38])
Changes:
Text added (Lines 247-248): “Exposure time was set to 5 ms, which provides optimal sensitivity to slow flows, such as in capillary perfused tissue{yuan2005determination}.”
Text added (Lines 273-276): “Averaged across all vessels blood flow index has decreased by 17.3±2.3% in the vasoconstricted condition and increased by 27.7±13.3% in the vasodilated condition. For the details on blood flow index measurements across the conditions see Supplementary Figure 2.”
Figure added: Figure 1 —figure supplement 2.
3. As it is understandable for a technical advance type of report, the study is rather descriptive, and the study design is lacking mechanistic examination of the observed synchronization and phase waves phenomena. Although the presence of fast electric or diffusive coupling that generates long-distance nephron-to-nephron communication was discussed. The applied vasoactive compounds, AngII and Ach have numerous extra-renal effects, for example on nerve activity which may be involved in cluster synchronization. It would be interesting to examine how the large cluster synchronization and phase waves are altered in kidneys one week after denervation. Since a vasoconstrictor state was associated with strengthened synchronization, one is wondering if stronger synchronization could be a sign of increased sympathetic nerve activity.
We are thankful to reviewer for the suggestion. It is indeed crucial to study synchronization under other conditions including the denervation. While it is out of the scope in the present study, we do aim to cover it in the future studies.
4. The physiological significance of nephron cluster synchronizations should be clearly stated.
We have clarified it in the text. Although we should bring it to the attention that since the topic is not well explored such discussion has a theoretical and speculative nature, requiring multiple studies to be done to come to specific conclusions.
Changes:
Text added (Lines 53-62): “Adaptive synchronization across the renal microvascular network would increase the efficiency and dynamic range of autoregulation by engaging more pre-glomerular resistance
{marsh2019nephron,zehra2021tubuloglomerular}, preventing transmission of high systemic pressure to the glomeruli, which could lead to progressive glomerular and vascular injury. In addition, long-distance synchronization would argue strongly that renal autoregulation is a distributed process that can ensure an optimized oxygenation-perfusion matching and adjust to various internal or environmental conditions {sosnovtseva2007synchronization, mitrou2015laser}. On the other hand, large scale in-phase synchronization might become a mechanism leading to glomerular injury as it would substantially increase local pressure variation {postnov2012dynamics}.”
In addition to the many strengths detailed above, I have a few suggestions for further improvement.
1. Vasoconstricted (Ang II infusion) and vasodilated conditions (Ach infusion) need to be validated experimentally by quantitative vessel diameter measurements. Vessel dimensions could be measured using the LSCI images or other higher resolution imaging modalities that could be simultaneously applied.
Please see Reply 1 above
2. It would be important to show if in the AngII condition the reduced blood flow and reduced sensitivity of the technique could contribute, at least in part, to the observed large clustering. Perhaps including an original, continuous tracing of RBF oscillations before and after AngII treatment would alleviate this concern.
Please see Reply 2 above
3. Renal denervation could address mechanistic details of the observed large clustering and phase wave phenomena.
Please see Reply 3 above
4. One or two additional sentences in the introduction should state the physiological significance of nephron synchronization.
Please see Reply 4 above
https://doi.org/10.7554/eLife.75284.sa2Article and author information
Author details
Funding
Novo Nordisk Fonden
- Dmitry Postnov
Lundbeckfonden
- Dmitry Postnov
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
All experimental protocols were approved by the Danish National Animal Experiments Inspectorate (License 2015-15-0201-00463) and were conducted according to the American Physiological Society guidelines.
Senior and Reviewing Editor
- Martin R Pollak, Harvard Medical School, United States
Version history
- Preprint posted: October 16, 2021 (view preprint)
- Received: November 4, 2021
- Accepted: May 6, 2022
- Accepted Manuscript published: May 6, 2022 (version 1)
- Version of Record published: May 17, 2022 (version 2)
- Version of Record updated: May 19, 2022 (version 3)
Copyright
© 2022, Postnov et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 603
- Page views
-
- 118
- Downloads
-
- 5
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
- Medicine
Background:
Although there are several efficacious vaccines against COVID-19, vaccination rates in many regions around the world remain insufficient to prevent continued high disease burden and emergence of viral variants. Repurposing of existing therapeutics that prevent or mitigate severe COVID-19 could help to address these challenges. The objective of this study was to determine whether prior use of bisphosphonates is associated with reduced incidence and/or severity of COVID-19.
Methods:
A retrospective cohort study utilizing payer-complete health insurance claims data from 8,239,790 patients with continuous medical and prescription insurance January 1, 2019 to June 30, 2020 was performed. The primary exposure of interest was use of any bisphosphonate from January 1, 2019 to February 29, 2020. Bisphosphonate users were identified as patients having at least one bisphosphonate claim during this period, who were then 1:1 propensity score-matched to bisphosphonate non-users by age, gender, insurance type, primary-care-provider visit in 2019, and comorbidity burden. Main outcomes of interest included: (a) any testing for SARS-CoV-2 infection; (b) COVID-19 diagnosis; and (c) hospitalization with a COVID-19 diagnosis between March 1, 2020 and June 30, 2020. Multiple sensitivity analyses were also performed to assess core study outcomes amongst more restrictive matches between BP users/non-users, as well as assessing the relationship between BP-use and other respiratory infections (pneumonia, acute bronchitis) both during the same study period as well as before the COVID outbreak.
Results:
A total of 7,906,603 patients for whom continuous medical and prescription insurance information was available were selected. A total of 450,366 bisphosphonate users were identified and 1:1 propensity score-matched to bisphosphonate non-users. Bisphosphonate users had lower odds ratios (OR) of testing for SARS-CoV-2 infection (OR = 0.22; 95%CI:0.21–0.23; p<0.001), COVID-19 diagnosis (OR = 0.23; 95%CI:0.22–0.24; p<0.001), and COVID-19-related hospitalization (OR = 0.26; 95%CI:0.24–0.29; p<0.001). Sensitivity analyses yielded results consistent with the primary analysis. Bisphosphonate-use was also associated with decreased odds of acute bronchitis (OR = 0.23; 95%CI:0.22–0.23; p<0.001) or pneumonia (OR = 0.32; 95%CI:0.31–0.34; p<0.001) in 2019, suggesting that bisphosphonates may protect against respiratory infections by a variety of pathogens, including but not limited to SARS-CoV-2.
Conclusions:
Prior bisphosphonate-use was associated with dramatically reduced odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations. Prospective clinical trials will be required to establish a causal role for bisphosphonate-use in COVID-19-related outcomes.
Funding:
This study was supported by NIH grants, AR068383 and AI155865, a grant from MassCPR (to UHvA) and a CRI Irvington postdoctoral fellowship, CRI2453 (to PH).
-
- Medicine
Background:
To systematically identify cell types in the human ligament, investigate how ligamental cell identities, functions, and interactions participated in the process of ligamental degeneration, and explore the changes of ligamental microenvironment homeostasis in the disease progression.
Methods:
Using single-cell RNA sequencing and spatial RNA sequencing of approximately 49,356 cells, we created a comprehensive cell atlas of healthy and degenerated human anterior cruciate ligaments. We explored the variations of the cell subtypes’ spatial distributions and the different processes involved in the disease progression, linked them with the ligamental degeneration process using computational analysis, and verified findings with immunohistochemical and immunofluorescent staining.
Results:
We identified new fibroblast subgroups that contributed to the disease, mapped out their spatial distribution in the tissue and revealed two dynamic trajectories in the process of the degenerative process. We compared the cellular interactions between different tissue states and identified important signaling pathways that may contribute to the disease.
Conclusions:
This cell atlas provides the molecular foundation for investigating how ligamental cell identities, biochemical functions, and interactions contributed to the ligamental degeneration process. The discoveries revealed the pathogenesis of ligamental degeneration at the single-cell and spatial level, which is characterized by extracellular matrix remodeling. Our results provide new insights into the control of ligamental degeneration and potential clues to developing novel diagnostic and therapeutic strategies.
Funding:
This study was funded by the National Natural Science Foundation of China (81972123, 82172508, 82372490) and 1.3.5 Project for Disciplines of Excellence of West China Hospital Sichuan University (ZYJC21030, ZY2017301).