Spectroscopy detects skeletal muscle microvascular dysfunction during onset of sepsis in a rat fecal peritonitis model

Sepsis is a dysregulated host inflammatory response to infection potentially leading to life-threatening organ dysfunction. The objectives of this study were to determine whether early microvascular dysfunction (MVD) in skeletal muscle can be detected as dynamic changes in microvascular hemoglobin (MVHb) levels using spectroscopy and whether MVD precedes organ histopathology in septic peritonitis. Skeletal muscle of male Sprague–Dawley rats was prepared for intravital microscopy. After intraperitoneal injection of fecal slurry or saline, microscopy and spectroscopy recordings were taken for 6 h. Capillary red blood cell (RBC) dynamics and SO2 were quantified from digitized microscopy frames and MVHb levels were derived from spectroscopy data. Capillary RBC dynamics were significantly decreased by 4 h after peritoneal infection and preceded macrohemodynamic changes. At the same time, low-frequency oscillations in MVHb levels exhibited a significant increase in Power in parts of the muscle and resembled oscillations in RBC dynamics and SO2. After completion of microscopy, tissues were collected. Histopathological alterations were not observed in livers, kidneys, brains, or muscles 6 h after induction of peritonitis. The findings of this study show that, in our rat model of sepsis, MVD occurs before detectable organ histopathology and includes ~ 30-s oscillations in MVHb. Our work highlights MVHb oscillations as one of the indicators of MVD onset and provides a foundation for the use of non-invasive spectroscopy to continuously monitor MVD in septic patients.


Methods
Animals. All animal procedures were approved by the University of Western Ontario Animal Care Committee of the University Council on Animal Care (Protocol #2015-105) and all methods were carried out in accordance with the guidelines and regulations set by the committee. Twenty male Sprague-Dawley rats were obtained from Charles River (Charles River Laboratories, Wilmington, MA). The rats were communally housed (2-3 rats per cage) under constant temperature and humidity with a 12-h light/dark cycle and allowed ad libitum access to food and water. At the time of use, the animals were approximately 7 weeks old. The ARRIVE guidelines were followed for reporting in vivo research 27 . Surgical preparation and physiological monitoring. Animals were randomized to a FIP or saline injection group. The rats were anesthetized with an intraperitoneal (IP) injection of sodium pentobarbital (65.64 mg/kg; Ceva Santé Animale, Libourne, Nouvelle-Aquitaine, France). The body temperature was maintained at 36.5-37.5 °C with a heat lamp and monitored using a rectal thermometer (Acorn Temp JKT; Cole-Parmer Instrument Co., Montreal, QC, Canada).
The left common carotid artery was cannulated with polyethylene 50 tubing (I.D. 0.58 mm, O.D. 0.965 mm; Becton Dickinson and Company, Sparks, MD, USA) that was pressurized with 0.9% NaCl (normal saline; pH 5.5, Na 154 mmol/L, Cl 154 mmol/L, Osmolarity 308 mOsmol/L; Baxter corporation, Mississauga, ON, Canada) with 1% heparin (Sandoz, Holzkirchen, Upper Bavaria, Germany). This allowed for monitoring of blood pressure and heart rate through a pressure transducer (Edwards Lifesciences) and Blood Pressure Analyzer (Digi-Med BPA-400; Micro-Med Inc., Louisville, KY, USA). The trachea was cannulated to allow mechanical ventilation (Inspira Advanced Safety Ventilator; Harvard Apparatus, Holliston, MA, USA) with a flow-controlled O 2 /N 2 mixture using standard settings for rats that were minimally modified based on arterial blood gas values (CG4 + cartridge, Abbott Point of Care, IL, USA; VetScan iSTAT 1 Handheld Analyzer, Abaxis, Union City, CA, USA). The right external jugular vein was cannulated using silastic tubing (I.D. 0.64 mm, O.D. 1.19 mm; Carpinteria, CA, USA). A continuous infusion of normal saline through the venous line was used for minimal fluid resuscitation (0.5 mL/h) and maintenance anesthesia. Depth of anesthesia was assessed by the palpebral reflex.
The right EDL was isolated, and its tendon was secured with a silk suture and severed distal to the ligature as previously described 11,28 . Animals randomized to the FIP group were injected with approximately 0.5 mL of a fecal slurry (approximately 0.25 g of autologous feces in normal saline) intraperitoneally through a small midline incision. The control animals received 0.5 mL of normal saline. The animal was placed on the microscope stage in the right lateral position. Suture on the EDL tendon was secured to the stage, keeping the muscle at its in situ length and orientation. The EDL was bathed with Plasma-Lyte A (37 ºC, pH 7.4, Na 140 mmol/L, Cl 98 mmol/L, Scientific Reports | (2022) 12:6339 | https://doi.org/10.1038/s41598-022-10208-w www.nature.com/scientificreports/ K 5 mmol/L, Mg 1.5 mmol/L, Acetate 27 mmol/L, Gluconate 23 mmol/L, Osmolarity 294 mOsmol/L; Baxter corporation, Mississauga, ON, Canada), isolated with an O 2 -impermeable polyvinylidene chloride film (Saran-Wrap; Dow Chemical Company, Midland, MI, USA), and gently compressed with a glass coverslip. The EDL was given 30 min to stabilize on the microscope stage and its microcirculation was recorded over a period of 6 h ( Supplementary Fig. S1 online). Gases and lactate levels in arterial blood taken from the carotid artery were measured at baseline and between each imaging period (VetScan iSTAT 1 Handheld Analyzer). Microvascular perfusion and responses are stable in our surgically exposed EDL preparation 29 .
Dual wavelength intravital video microscopy. The  Tissue collection and processing. After completion of intravital imaging in another cohort of rats, arterial blood was collected and mixed with 3% acetic acid and 1% crystal violet in a 5∶44∶1 ratio. Hemocytometer counts were averaged from 6 separate samples. Rats were overdosed with pentobarbital and perfused transaortically with phosphate-buffered saline, followed by 4% paraformaldehyde (PFA). Transaortic perfusion-fixation allowed for collection of the brain, left extensor digitorum longus, left lateral lobe of the liver and kidneys but not lungs or hearts which were bypassed by the solutions. The muscles, livers and kidneys were immersion-fixed in 4% PFA for 3 days, followed by dehydration and embedding in paraffin. The thin-sectioned tissues (5 μm) were stained with hematoxylin and eosin (H&E) and assessed for tissue damage using an upright microscope (Olympus BX-50, Hitachi SXGA ECCO HV-F22 camera, Northern Eclipse software, Olympus Corp.).

Brain immunohistochemistry.
Brains were post-fixed in 4% PFA for 24 h at 4 ºC and transferred into 30% sucrose solution at 4 ºC. Brains were frozen, sectioned at 30 µm and every 6 th coronal slice (spanning the whole brain) was immunostained for rat MHC class II with anti-rat RT1B (clone OX-6; BD Biosciences, Mississauga, ON, Canada) to detect microglial activation 30 . Specifically, free-floating sections were treated with 1% H 2 O 2 and blocked with horse serum (Vector Laboratories, Inc., Burlingame, CA, USA). The sections were incubated with anti-rat RT1B in blocking serum overnight at 4 °C, followed by 1-h incubation with biotinylated 2° antibody (BA-2000; Vector Laboratories, Inc.) in serum at room temperature. Sections were then incubated with avidin-biotin complex (Vector Laboratories, Inc.) for 1 h followed by 0.05% diaminobenzidine. After washing, sections were mounted on slides with 0.3% gelatin and air-dried, followed by dehydration, clearing in xylenes, and mounting with DePex mounting medium (Electron Microscopy Science, Hatfield, PA, USA). After examination under high magnification, stitched images of brain sections were acquired using an upright microscope (Nikon Eclipse Ni-E, Nikon DS Fi2 color camera, NIS Elements Imaging, Mississauga, ON, Canada).
Offline analysis. Functional capillary density (capillaries/mm) was determined by counting capillaries that crossed three horizontal lines on a 450 μm by 340 μm field of view and capillary RBC flow was categorized as continuous, intermittent (stopped for > 3 s at least once or reversed flow) or completely stopped during a 30-s segment of the recording. Capillary location and geometry were determined from functional images generated from the video sequence 31,32 . Analysis of RBC SO 2 and flow through capillaries in each field of view was performed as previously described 33 . Briefly, software extracted frame-by-frame light intensity data from the centre line of in-focus capillaries over the length of the sample to generate a space-time image (STI) of the passage of RBCs through each capillary. To calculate the OD of each RBC, the incident light intensity was determined from the light intensity of the plasma gaps between RBCs in a capillary. OD of each RBC at both λ was calculated as the ratio of the intensity of transmitted light through the RBC to the incident intensity. RBC hemoglobin SO 2 was determined by calculating the 438 nm/450 nm OD ratio with constants derived from in vivo calibration and displayed as SO 2 colour-coded STIs. Capillary hemodynamics were calculated from the analysis of the 450 nm STI yielding frame-by-frame data for capillary RBC velocity (μm/s), lineal density (RBC/mm)/hematocrit (%) and RBC supply rate (RBC/s) 29 . Of note, capillaries with stopped flow were assigned a velocity and supply rate of 0. www.nature.com/scientificreports/ The variation in hemoglobin content in the volume of tissue (1 mm thickness of muscle encompassing multiple capillary modules as well as arteriolar and venular trees) associated with each field of view was determined from the acquired video sequence and from the CWS spectroscopy data. The mean intensity of all pixels in each video frame of the 1-min sequence was calculated for the 450-nm λ. The resulting mean frame intensity (MFI) time series was processed to generate the change in OD relative to the first time point, i.e. ΔOD 450 (t). The CWS spectroscopy data was processed in the same manner for all λ values yielding an array of ΔOD λ (t). For this study, ΔOD values for the 6 isosbestic λ (452 nm, 500 nm, 530 nm, 545 nm, 570 nm and 584 nm; bandpass 2 nm) were calculated; note, these are the isosbestic points in the visible spectrum ( Fig. 1B). At these isosbestic λ, the extinction coefficients are independent of the hemoglobin SO 2 and hence only reflect dynamic changes in MVHb content in the sample volume. We averaged signals from the 6 isosbestic points to increase the accuracy of our MVHb measurements. Fourier analysis was applied to the MFI and spectroscopy data to quantify the Power of low-frequency oscillations observed in MVHb content. Specifically, the Fourier transforms of the MFI and spectroscopy data over time were approximated using the built-in fast Fourier transform algorithm in MATLAB (MathWorks, Natick, MA, USA). The Power was quantified as the square of the signal amplitude and the peak Power was calculated by taking the maximum of the Power spectrum in the 0.03-0.06 Hz range (corresponding periods of 16 s -32 s). The peak Power of these low-frequency oscillations from MFI and CWS recordings was compared between FIP and control groups.

Statistics.
Analysis was performed using GraphPad Prism version 7.00, (GraphPad Software; La Jolla, CA, USA). A two-way ANOVA with a Bonferroni post-hoc test was conducted to evaluate differences between the control and FIP groups over the three imaging periods. All values are reported as mean (standard error of the mean; SEM). Statistical significance was defined as P < 0.05.

Results
Animal characteristics and physiological parameters during intravital microscopy. The weights of control and FIP animals used in the experiment did not significantly differ (0.174 kg (0.007) vs 0.182 kg (0.003), P = 0.322). After surgical preparation for IVVM and during the 6-h imaging period, none of the animals died. Macrohemodynamics and physiological parameters were monitored during the study to determine the timing of their deterioration (Table 1). MAP remained stable in the control group through the 6-h imaging period. MAP was stable in the FIP group from 0.5-4 h after infection but was significantly decreased during the 4.5-6 h imaging period compared with control animals (P < 0.0001; Table 1). Thus, macrovascular changes occurred closer to the experimental endpoint in the FIP group. Lactate levels, a marker of tissue hypoperfusion, were significantly increased in rats with FIP compared with controls at 4.5 h (P = 0.0015; Table 1) and 6 h after infection (P < 0.0001). At endpoint, the feces-injected animals had systemic leukocyte counts that were significantly lower than the saline-injected animals (2.8 (1.1) vs 6.6 (0.9) × 10 6 cells/mL, respectively, P < 0.0001). Based on their decreased MAP, increased lactate levels (> 2 mmol/L) 34 , and leukopenia 35,36 , rats with FIP in our study were progressing towards a septic state at the experimental endpoint. We were, however, able to maintain stable PaO 2 values in both groups through the imaging period, which did not significantly differ between the control and FIP groups (Table 1). Also, the arterial partial pressure of O 2 to fraction of inspired O 2 ratio (PaO 2 /FiO 2 ) did not significantly change over time ( Table 1), indicating that ventilator settings used in this study did not cause acute respiratory distress syndrome, which can affect RBC SO 2 .
Functional capillary density in skeletal muscle declines by 4 h after induction of fecal peritonitis. The number of capillaries with continuously flowing RBCs was significantly decreased in the FIP group compared with control animals during the 2.5-4 h imaging period (P = 0.0002; Fig. 2) as capillaries progressed to intermittent (Fig. 3) and stopped-flow states. The decline in continuously flowing capillaries in animals with FIP continued 4.5-6 h after infection (P < 0.0001) while the number of capillaries that had stopped flow significantly increased (P < 0.0001). This reduced capillary perfusion was in concert with increased lactate levels ( Table 1) and preceded significant decrease in MAP, highlighting a mismatch between microhemodynamic and macrohemodynamic deterioration in the fecal peritonitis model of sepsis. Table 1. Mean arterial pressure, lactate, PaO 2 and PaO 2 /FiO 2 ratios. Data are reported as mean (SEM) and assessed using two-way ANOVA with Bonferroni correction. *P < 0.05 compared with control, n = 10 rats.   At baseline, RBC velocity was not significantly different between control rats and animals with FIP (Fig. 4A). However, capillary RBC velocities were significantly decreased in the FIP group 2.5-4 h (P < 0.0001) and 4.5-6 h (P < 0.0001) after infection. The distribution of capillary RBC velocities was similar in control and FIP animals in the first imaging period (Fig. 4A) but shifted towards lower velocities in the FIP group by 4.5-6 h, while the velocities in the control group became more widely distributed-an effect of surgical exposure of the muscle. Of note, we did not see increased numbers of hyperemic capillaries by 6 h after fecal injection as was observed in rats 24 h post-CLP 11 , indicating that MVD manifests differently as sepsis progresses. Capillary hematocrit was similar between the two groups during baseline recordings and 2.5-4 h after infection (Fig. 4B). At 4.5-6 h after infection, hematocrit in the capillaries that were still flowing in the FIP animals was slightly higher compared to controls (P = 0.0033). However, the hematocrit values collected during this 4.5-6 h time point did not change from baseline hematocrit in the FIP group (repeated two-way ANOVA with Tuckeys multiple comparisons test, P = 0.5808). Thus, overall, capillary hematocrit in the FIP group appeared stable. While RBC supply rates were similar between the two groups at baseline (Fig. 4C), there was a significant decline with FIP 2.5-4 h (P < 0.0001) and 4.5-6 h (P < 0.0001) after infection. RBC SO 2 was also similar at baseline between the two groups (Fig. 4D) and significantly decreased in FIP animals in the second (P < 0.0001) and third (P < 0.0001) imaging period. By 2.5-4 h after fecal injection, oscillatory flow patterns begun to emerge in capillary RBC velocity, hematocrit and supply rate at roughly 2 cycles/min (Fig. 5B). These oscillations were not readily observed in control animals, which had overall steady levels of RBCs flowing with some variability on a second-to-second basis (Fig. 5A). The ~ 30-s cycles also appeared in RBC SO 2 in capillaries of animals with fecal peritonitis (Fig. 5B), suggesting the microvasculature lost the ability to maintain constant SO 2 . In areas of the muscle characterized by oscillatory RBC flow, some capillaries had completely stationary RBCs, indicating that dysfunction within a capillary module may be heterogeneous.
Oscillations in microvascular hemoglobin levels, derived from spectroscopy data, indicate early microvascular dysfunction in skeletal muscle of rats with fecal peritonitis. The EDL is a relatively thin muscle, allowing for transillumination with negligible light scatter. ODs were calculated from light intensity data of the captured IVVM frames and spectroscopy profiles of the EDL. OD at isosbestic λ = extinction  . Through the 6-h imaging period, the mean ΔOD of 6 isosbestic λ captured with CWS from the muscle remained relatively stable at high frequencies in control animals, similar to the MFI from IVVM recordings. Graphs from a representative field of view are shown in top rows of Fig. 6A and B.
In the animals with fecal peritonitis, however, vasomotion was detected as low-frequency oscillations by 2.5-4 h after infection in both MFI and CWS recordings (bottom rows in Fig. 6A,B).
To summarize the type of data exemplified in Fig. 6 and indirectly Fig. 5, Fourier analysis was performed on the light intensity data collected with microscopy (MFI) and spectroscopy (CWS). Fourier analysis of the OD measurements from the MFI and CWS data exhibited a significant increase in Power at low frequencies in hemoglobin levels associated with vasomotion by 4 h after fecal injection (Fig. 7). Outlier analysis was performed using the ROUT method and 33 outliers out of 219 values were removed from MFI data over the three timepoints www.nature.com/scientificreports/ and 38 outliers out of 219 values were removed from the spectroscopy data. Fields of view that contained an outlier during the first, second or third imaging period were not used to preserve data matching across time. Power (amplitude 2 ) of the low-frequency oscillations from in vivo MFI signals was significantly increased in the FIP group compared to the control group in the second (P = 0.0252) and third imaging period (P = 0.0312; Fig. 7A). Likewise, the Power of low-frequency oscillations from CWS data was significantly increased in the FIP group compared to controls in the second (P = 0.0498) and third imaging period (P = 0.0066; Fig. 7B). Of note, the average peak period of these oscillations in the FIP groups was 28.8 (2.7) seconds. The Power of these oscillations was displayed as a heat map across all the fields of view for control (Fig. 7C) and FIP (Fig. 7D) rats. The oscillations appeared in several fields of view from muscles of the FIP group during the 2.5-4.5 h imaging period and were particularly prominent in rat 2 and rat 4, as labelled in the figure (Fig. 7D). However, some fields of view from FIP rats did not exhibit these low-frequency oscillations (Fig. 7C,D), which may have resulted from substantially reduced flow in these tissue volumes. These low-frequency oscillations in the spectroscopy data closely reflected oscillations in RBC velocities, hematocrit, supply rates and SO 2 in individual muscle capillaries of rats with fecal peritonitis (Fig. 5). Thus, oscillations in hemoglobin levels from the spectroscopy data appear to be one of the early indicators of microvascular derangement in sepsis. with H&E ( Fig. 8A-C). The EDL muscles from animals with FIP did not exhibit increased tissue inflammatory cell infiltration. Erythrocyte trapping within the glomerular microvasculature did not rise above control levels.
As well, livers from FIP and saline control animals did not exhibit excessive leukocyte accumulation and parenchymal cell architecture was not altered. No vacuolization of the hepatic parenchymal cells was noted in the two groups. These findings suggest that, in our rat FIP model of sepsis, microvascular dysfunction occurs despite lack of detectable organ histopathology in the early stages of the infection. Immunohistochemical labelling of brain sections for rat MHC class II did not show signs of pro-inflammatory microglial activation in brains of FIP rats with careful inspection under high magnification. While coronal slices spanning the entire brain were examined, sections encompassing the striatum and hippocampus are shown (Fig. 8D). Endothelin-1-induced brain infarct in a rat model served as positive control for this technique 30 .

Discussion
Septic microvascular dysfunction manifests as decreased arteriolar responsivity 37 , increased stopped-flow capillaries, reduced capillary RBC supply rates, reduced RBC SO 2 , and increased capillary response time to hypoxia 13 , all tied to impaired microvascular autoregulation. We believe that high-amplitude, low-frequency oscillations in microvascular hemoglobin levels may also be one of the early indicators of microvascular failure in sepsis and reflect the onset of impaired microvascular autoregulation in the early stages of severe infection. Furthermore, we demonstrate that these oscillations can be detected non-invasively in skeletal muscle remote to the infectious locus. As well, we show that these oscillations precede macrohemodynamic deterioration and histopathologic alterations in vital organs. This work highlights the potential use of non-invasive spectroscopy techniques to www.nature.com/scientificreports/ continuously monitor MVD in the early stages of sepsis using easily accessible regions of the body such as skeletal muscle. Several hypotheses have been put forth to explain the pathophysiological processes that lead to organ injury and dysfunction in sepsis. Our work centers on the idea that microvascular derangement happens early in sepsis, affecting O 2 delivery and consequently instigating tissue injury. Maldistribution of RBC flow through capillary networks in septic rat skeletal muscle, manifested as increased flow heterogeneity with stopped-and fast-flow capillaries, leads to increased diffusion distance for O 2 in the tissue 11 . Thus, the consequence of septic microvascular dysfunction is a mismatch between O 2 demand and supply. Conversely, it is also believed that early mitochondrial dysfunction in sepsis diminishes O 2 utilization and drives progression to multi-organ failure 38 . However, increased O 2 extraction in septic rat skeletal muscle suggests that mitochondrial dysfunction does not occur in the early stage of sepsis in our animal model 11 . This is in line with our results, which show reduced RBC SO 2 in skeletal muscle capillaries despite normal systemic PaO 2 in the FIP rats.
Another proposed mechanism of tissue hypoperfusion in sepsis is blockage via accumulation of adherent leukocytes in capillaries and post-capillary venules, yet leukocyte-endothelial interactions in the EDL of septic rats were observed to be reduced 6-48 h after CLP 39 . Thus, microvascular obstruction by leukocytes is an unlikely cause of the reduced muscle perfusion we observed in our feces-injected rats. Microthrombi may also contribute to capillary plugging and disseminated intravascular coagulation is a common feature of sepsis 40 . In fact, stopped-flow capillaries in the EDL from septic mice 6-7 h after induction of fecal peritonitis are marked by P-selectin-dependent platelet adhesion and fibrin formation 41 . However, in our study, RBC supply rate and SO 2 were significantly decreased by 4 h after fecal injection while the number of stopped-flow capillaries did not significantly increase till closer to the experimental endpoint. Thus, it appears that regulatory failure precedes any possible plugging of capillaries.
Tissue hypoperfusion in sepsis is believed to be tied to impaired signaling from capillaries to their supplying arterioles. It is clear that electrical communication between coupled vascular cells is diminished in early sepsis, evidenced by blunted conducted responses in mouse skeletal muscle 37,[42][43][44][45] . Diminished communication between capillaries and arterioles in O 2 demand-supply coupling may be based on ATP signalling as regulatory failure in sepsis was found to involve loss of O 2 sensing in hypoxic capillaries with an impaired RBC O 2 -dependent efflux of ATP 13 . ATP release as RBCs desaturate 46 is proposed to be a purinergic stimulus 47 for conducted hyperpolarization and the resulting functional hyperemia 15 . ATP efflux is also induced by RBC deformation 48,49 and, indeed, septic RBCs were found to be less deformable 50 . Taken together, these studies suggest that sepsis impairs conducted hyperpolarization along the capillary endothelium to the upstream arteriole, and along the arteriolar tree. We believe such impaired vascular autoregulation may underlie the observed loss of fine second-to-second control of hemoglobin delivery and the dominance of a fundamental low frequency of arteriolar vasomotion in our current study.
The goals of our study were to determine whether MVD can be detected with spectroscopy and whether spectroscopy signals qualitatively and temporally correspond to microhemodynamic changes and mean frame intensity oscillations observed with IVVM. Both methods detected temporal variability in MVHb levels in the early stages of fecal peritonitis, manifested as low-frequency OD oscillations by 4 h post-infection. The peak Power of these oscillations in spectroscopic measurements was significantly increased in FIP rats. However, in some areas of the muscle, no oscillations were observed, which may be due to capillary networks with RBC flow too low to generate a detectable oscillatory MVHb signal within the volume of tissue. This indicates that MVD is not entirely uniform across the muscle. Since the oscillations in capillary RBC flow synchronously spanned more than one 450 μm by 340 μm field of view and MVHb oscillations from spectroscopy data were derived from a volume of tissue spanning the entire EDL thickness, we believe rhythmic diameter changes occurred in larger arterioles.
The ~ 30-s microhemodynamic oscillations in the FIP group include periodically stopped flow, slower flow and reversed flow of RBCs within the vascular network. Alternating forward and reversed RBC flow through capillaries is due to pressure gradient changes in the microvascular network. In skeletal muscle, capillaries running between an arteriole and a venule are referred to as a capillary unit or module 51,52 . Of the capillaries supplied by a given terminal arteriole, roughly half drain into one venule and the other half (running in the opposite direction along muscle fibres) drain into another venule. Likewise, a given post capillary venule drains capillary modules from two or more terminal arterioles. This arrangement repeats along muscle fibres, such that a given fibre is supplied by multiple capillary modules 51 . These capillary modules are connected in series forming long columnar capillary structures, i.e. capillary fascicles 53 . The structure of the capillary fascicle appears to correspond to the muscle fascicle. Given the length of the capillary fascicle, many of the terminal arterioles supplying these modules originate from different branches of the arteriolar tree. Thus, for the microvasculature to maintain proper RBC flow direction in each capillary unit, blood pressures must be tightly controlled by coordinated regulation within multiple branches of the arteriolar tree. When regulation within the arteriolar tree is not coordinated, high pressure in one terminal arteriole may cause flow reversals in capillary modules normally supplied by adjacent terminal arterioles or low arteriolar pressures may result in flow reversals supplied by the venular tree. We observe such uncoordinated regulation as reversed flow of desaturated RBCs in muscle capillaries after induction of peritonitis (Fig. 3). This pathologic flow pattern contributed to the lower RBC SO 2 observed in our FIP rats. Since this lack of coordination of arteriolar constriction or dilation would affect modular flow over long distances, it possibly contributed to the observed ~ 30-s oscillations in our IVVM and spectroscopy data-this will be an area for future study.
It has been actively debated whether MVD precedes tissue injury in sepsis, and our work aligns with the notion that it does. As the development of acute respiratory distress syndrome 54 , hepatic dysfunction 3 , and/or acute kidney injury 4,5 is common in septic patients and skeletal muscle wasting in sepsis survivors can result in long-term mobility problems 6  www.nature.com/scientificreports/ analysis. At 24 h after induction of fecal peritonitis in mice, considerable variability was observed in tissue injury between animals and even different organs in a given animal 55 . In some reports, tissue injury was detected as early as 6 h after cecal ligation and puncture in a mouse model 16,17 . Our model differs from these studies in that there is no necrotic intestinal tissue induced by puncture as an additional inflammatory source. Furthermore, our rats were not conscious through the 6-h study period, were supported on a ventilator and continuously provided intravenous fluids over the duration of IVVM. These factors may have contributed to the lack of observed histopathology in several organs of our feces-injected rats in the early stage of the infection. Notably, the severity of the infectious insult in different animal models is hard to calibrate and directly compare, and thus the strength of our study is that MVD analysis and histology were done in the same animals. Septic patients may develop acute brain dysfunction, referred to as sepsis-associated encephalopathy, which is linked to greater mortality 56 , and sepsis survivors are at risk of long-term cognitive impairment 57 . The mechanisms that drive septic brain injury are not completely understood, and inflammation 58,59 , oxidative stress 60 , MVD 61 , and dynamic variability in MAP and arterial CO 2 levels 62 are suspected. Animal studies attempted to dissect the key inflammatory players in septic brains, with microglial activation as an important driver 63 , which has also been suggested in septic patients 59 . However, conclusions on the timing of microglial activation after induction of sepsis in animal models vary widely [64][65][66][67][68][69][70][71] . In our study, we examined inflammation in brains collected 6 h after induction of fecal peritonitis and found no evidence of microglial activation at this early time-point. It is important to understand whether peripheral MVD (such as in skeletal muscle) occurs before brain inflammation and if MVD monitoring could guide early intervention measures to reduce the risk of neurological dysfunction in sepsis survivors. Furthermore, it is unknown if MVD in the cerebrum itself causes brain damage, leading to acute loss of cognitive function or chronic long-term cognitive dysfunction.
Analogous to our spectroscopy approach to detect MVD preclinically, NIRS has been used to assess microvascular function in septic patients 25,72 . Our study supports the use of NIRS as a diagnostic tool for non-invasive continuous monitoring of MVD in sepsis. Furthermore, the analysis and interpretation of the oscillatory dynamics of MVHb signal in our study via direct comparison to in vivo microscopy observations could be applied to NIRS to garner more information from microvascular monitoring in critically ill patients. In fact, we have already begun monitoring skeletal muscle microvascular hemoglobin levels in the critically ill population using NIRS and have shown the feasibility of this approach with the presence of low frequency oscillations in some of the septic patients 26 .
Using a fecal peritonitis model of sepsis, the current study showed that dynamic changes in microvascular hemoglobin levels detected with spectroscopy correspond to microcirculatory deterioration observed with in vivo microscopy, and that MVD precedes macrohemodynamic changes and organ histopathology. This work provides a foundation for using non-invasive spectroscopy with deeper tissue penetrance, such as NIRS, as a diagnostic tool to detect and monitor MVD in septic patients. Spectroscopic detection of MVD may guide early sepsis screening and diagnostic pathways, allowing for prompt therapeutic interventions. By targeting microcirculatory preservation as well as restoration early, sepsis mortality and risk of long-term morbidity may be reduced, warranting further study in clinical populations.

Data availability
The datasets used and/or analysed in this study are available from the corresponding author on reasonable request.