Reproducibility of Laser Doppler Flowmetry in gingival microcirculation. A study on six different protocols

OBJECTIVES

Laser Doppler Flowmetry (LDF) is a non-invasive technique for the assessment of tissue blood flow, but increased reproducibility would facilitate longitudinal studies. The aim of the study was to assess the interday reproducibility of Laser Doppler Flowmetry (LDF) at rest, at elevated local temperatures, and with the use of the vasodilator Methyl Nicotinate (MN) in six interconnected protocols for the measurement of the blood supply to the microvascular bed of the gingiva.

METHODS

Ten healthy volunteers were included. Interweek LDF measurements with custom-made acrylic splints were performed. Six protocols were applied in separate regions of interest (ROI): 1; basal LDF, 2; LDF with thermoprobe 42 °C, 3; LDF with thermoprobe 45 °C, 4; LDF with thermoprobe 42 °C and MN, 5; LDF with thermoprobe 45 °C and MN and 6; LDF with MN.

RESULTS

Intra-individual reproducibility was assessed by the within-subject coefficient of variation (wCV) and the intraclass correlation coefficient (ICC). Basal LDF measurements demonstrated high reproducibility with wCV 11.1 in 2 min and 10.3 in 5 min. ICC was 0.9 and 0.92. wCV after heat and MN was 4.9-10.3 and ICC 0.82-0.93. The topically applied MN yielded increased blood flow.

CONCLUSION

This is the first study evaluating the reproducibility of basal LDF compared to single or multiple vasodilatory stimuli in gingiva. Multiple collector fibers probe and stabilizing acrylic splints are recommended. Vasodilatory stimulation showed a tendency toward higher reproducibility. Furthermore, MN yields vasodilation in gingiva.

Worse histopathology and prognosis in women with breast cancer diagnosed during the second trimester of pregnancy

BACKGROUND

Evidence suggests that women with breast cancer diagnosed during pregnancy (PrBC) and within 2 years of delivery (PPBC) have similar survival compared to women diagnosed not near pregnancy if adjusted for stage and subtype. To investigate whether this is true for all subtypes and for both pregnancy and post-delivery periods, we examined clinicopathologic features and survival in women with breast cancer by trimesters and 6-month post-delivery intervals.

MATERIALS AND METHODS

Women diagnosed with invasive breast cancer during 1992-2018 at ages 18-44 years were identified in the Swedish Cancer Register, with information on childbirths from the Swedish Multi-Generation Register and clinical data from Breast Cancer Quality Registers. Each woman with PrBC or PPBC was matched 1 : 2 by age and year to comparators diagnosed with breast cancer not near pregnancy. Distributions of stage, grade, and surrogate subtypes were compared. Adjusted hazard ratios (HRs) with 95% confidence intervals (CIs) for breast cancer mortality were estimated using Cox regression.

RESULTS

We identified 1430 women with PrBC and PPBC (181 during pregnancy, 499 during the first and 750 during the second year after delivery). Compared to 2860 comparators, women with PrBC and PPBC in the first year after delivery had a significantly higher proportion of luminal human epidermal growth factor receptor 2 (HER2)-positive, HER2-positive and triple-negative tumours, and more advanced stage at diagnosis. After adjustment for age, year, parity, country of birth, hospital region, subtype, and stage, women diagnosed during the second trimester had a worse prognosis than matched comparators (HR 1.8, 95% CI: 1.0-3.2).

CONCLUSIONS

Women diagnosed during pregnancy or within the first year after delivery have a worse prognosis than women diagnosed not near pregnancy due to adverse tumour biology and advanced stage at diagnosis. A worse prognosis for women diagnosed during the second trimester remained after multivariable adjustment, possibly reflecting difficulties to provide optimal treatment during ongoing pregnancy.

Skin perfusion and oxygen saturation after mastectomy and radiation therapy in breast cancer patients

The pathophysiological mechanism behind complications associated with postmastectomy radiotherapy (PMRT) and subsequent implant-based breast reconstruction are not completely understood. The aim of this study was to examine if there is a relationship between PMRT and microvascular perfusion and saturation in the skin after mastectomy and assess if there is impaired responsiveness to a topically applied vasodilator (Methyl nicotinate - MN). Skin microvascular perfusion and oxygenation >2 years after PMRT were measured using white light diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF) in the irradiated chest wall of 31 women with the contralateral breast as a control. In the non-irradiated breast, the perfusion after application of MN (median 0.84, 25th-75th centile 0.59-1.02 % RBC × mm/s) was higher compared to the irradiated chest wall (median 0.51, 25th-75th centile 0.21-0.68 % RBC × mm/s, p < 0.001). The same phenomenon was noted for saturation (median 91 %, 25th-75th centile 89-94 % compared to 89 % 25th-75th centile 77-93 %, p = 0.001). Eight of the women (26%) had a ≥10 % difference in skin oxygenation between the non-irradiated breast and the irradiated chest wall. These results indicate that late microvascular changes caused by radiotherapy of the chest wall significantly affect skin perfusion and oxygenation.

Association between cardiovascular risk profile and impaired microvascular function in a Swedish middle-aged cohort (the SCAPIS study)

AIM

The aim was to investigate the relationship between microvascular function, cardiovascular risk profile, and subclinical atherosclerotic burden.

METHODS

The study enrolled 3809 individuals, 50-65 years old, participating in the population-based observational cross-sectional Swedish CArdioPulmonary bioImage Study (SCAPIS). Microvascular function was assessed in forearm skin using an arterial occlusion and release protocol determining peak blood oxygen saturation (OxyP). Cardiovascular risk was calculated using the updated Systematic Coronary Risk Evaluation (SCORE2; 10-year risk of fatal and non-fatal CVD events). OxyP was compared with Coronary Artery Calcification Score (CACS) and to plaques in the carotid arteries.

RESULTS

Individuals with OxyP values in the lowest quartile (Q1; impaired microvascular function), had a mean SCORE2 of 5.8% compared to 3.8% in those with the highest values of OxyP (Q4), a relative risk increase of 53%. The risk of having a SCORE2 > 10% was 5 times higher for those in Q1 (OR: 4.96 95% CI: 2.76-8.93) vs. Q4 when adjusting for BMI and hsCRP. OxyP was lower in individuals with CACS > 0, and in those with both carotid plaques and CACS >0, =compared with individuals without subclinical atherosclerotic burdens (87.5 ± 5.6% and 86.9 ± 6.0%, vs 88.6 ± 5.8%, p < 0.01).

CONCLUSION

In a population without cardiovascular disease or diabetes mellitus, impaired microvascular function is associated with cardiovascular risk profiles such as higher SCORE2 risk and CACS. We suggest that OxyP may serve as a microcirculatory functional marker of subclinical atherosclerosis and CVD risk, that is not detected by structural assessments.

Absorption and reduced scattering coefficients in epidermis and dermis from a Swedish cohort study

Significance

Knowledge of optical properties is important to accurately model light propagation in tissue, but in vivo reference data are sparse.

Aim

The aim of our study was to present in vivo skin optical properties from a large Swedish cohort including 3809 subjects using a three-layered skin model and spatially resolved diffuse reflectance spectroscopy (Periflux PF6000 EPOS).

Approach

Diffuse reflectance spectra (475 to 850 nm) at 0.4 and 1.2 mm source-detector separations were analyzed using an inverse Monte Carlo method. The model had one epidermis layer with variable thicknesses and melanin-related absorptions and two dermis layers with varying hemoglobin concentrations and equal oxygen saturations. The reduced scattering coefficient was equal across all layers.

Results

Median absorption coefficients (mm - 1 ) in the upper dermis ranged from 0.094 at 475 nm to 0.0048 at 850 nm and similarly in the lower dermis from 0.059 to 0.0035. The reduced scattering coefficient (mm - 1 ) ranged from 3.22 to 1.20, and the sampling depth (mm) ranged from 0.23 to 0.38 (0.4 mm separation) and from 0.49 to 0.68 (1.2 mm separation). There were differences in optical properties across sex, age groups, and BMI categories.

Conclusions

Reference material for skin optical properties is presented.

Workload and sex effects in comprehensive assessment of cutaneous microcirculation

INTRODUCTION

Workload and sex-related differences have been proposed as factors of importance when evaluating the microcirculation. Simultaneous assessments with diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF) enable a comprehensive evaluation of the microcirculation. The aim of the study was to compare the response between sexes in the microcirculatory parameters red blood cell (RBC) tissue fraction, RBC oxygen saturation, average vessel diameter, and speed-resolved perfusion during baseline, cycling, and recovery, respectively.

METHODS

In 24 healthy participants (aged 20 to 30 years, 12 females), cutaneous microcirculation was assessed by LDF and DRS at baseline, during a workload generated by cycling at 75 to 80 % of maximal age-predicted heart rate, and recovery, respectively.

RESULTS

Females had significantly lower RBC tissue fraction and total perfusion in forearm skin microcirculation at all phases (baseline, workload, and recovery). All microvascular parameters increased significantly during cycling, most evident in RBC oxygen saturation (34 % increase on average) and perfusion (9-fold increase in total perfusion). For perfusion, the highest speeds (>10 mm/s) increased by a factor of 31, whereas the lowest speeds (<1 mm/s) increased by a factor of 2.

CONCLUSION

Compared to a resting state, all studied microcirculation measures increased during cycling. For perfusion, this was mainly due to increased speed, and only to a minor extent due to increased RBC tissue fraction. Skin microcirculatory differences between sexes were seen in RBC concentration and total perfusion.

Speed-resolved perfusion imaging using multi-exposure laser speckle contrast imaging and machine learning

SIGNIFICANCE

Laser speckle contrast imaging (LSCI) gives a relative measure of microcirculatory perfusion. However, due to the limited information in single-exposure LSCI, models are inaccurate for skin tissue due to complex effects from e.g. static and dynamic scatterers, multiple Doppler shifts, and the speed-distribution of blood. It has been demonstrated how to account for these effects in laser Doppler flowmetry (LDF) using inverse Monte Carlo (MC) algorithms. This allows for a speed-resolved perfusion measure in absolute units %RBC × mm/s, improving the physiological interpretation of the data. Until now, this has been limited to a single-point LDF technique but recent advances in multi-exposure LSCI (MELSCI) enable the analysis in an imaging modality.

AIM

To present a method for speed-resolved perfusion imaging in absolute units %RBC × mm/s, computed from multi-exposure speckle contrast images.

APPROACH

An artificial neural network (ANN) was trained on a large simulated dataset of multi-exposure contrast values and corresponding speed-resolved perfusion. The dataset was generated using MC simulations of photon transport in randomized skin models covering a wide range of physiologically relevant geometrical and optical tissue properties. The ANN was evaluated on in vivo data sets captured during an occlusion provocation.

RESULTS

Speed-resolved perfusion was estimated in the three speed intervals 0 to 1    mm / s , 1 to 10    mm / s , and > 10    mm / s , with relative errors 9.8%, 12%, and 19%, respectively. The perfusion had a linear response to changes in both blood tissue fraction and blood flow speed and was less affected by tissue properties compared with single-exposure LSCI. The image quality was subjectively higher compared with LSCI, revealing previously unseen macro- and microvascular structures.

CONCLUSIONS

The ANN, trained on modeled data, calculates speed-resolved perfusion in absolute units from multi-exposure speckle contrast. This method facilitates the physiological interpretation of measurements using MELSCI and may increase the clinical impact of the technique.

Risk of coronary stenosis after adjuvant radiotherapy for breast cancer

Purpose

Adjuvant radiotherapy (RT) for breast cancer is associated with an increased risk of ischemic heart disease. We examined the risk of coronary artery stenosis in a large cohort of women with breast cancer receiving adjuvant RT.

Methods

A cohort of women diagnosed with breast cancer between 1992 and 2012 in three Swedish health care regions (n = 57,066) were linked to the Swedish Coronary Angiography and Angioplasty Registry (SCAAR) to identify women receiving RT who subsequently underwent a percutaneous coronary intervention (PCI) due to coronary stenosis. Cox regression analyses were performed to examine risk of a coronary intervention and competing risk analyses were performed to calculate cumulative incidence.

Results

A total of 649 women with left-sided breast cancer and 494 women with right-sided breast cancer underwent a PCI. Women who received left-sided RT had a significantly higher risk of a PCI in the left anterior descending artery (LAD) compared to women who received right-sided RT, hazard ratio (HR) 1.44 (95% confidence interval [CI] 1.21–1.77, p < 0.001). For the proximal, mid, and distal LAD, the HRs were 1.60 (95% CI 1.22–2.10), 1.38 (95% CI 1.07–1.78), and 2.43 (95% CI 1.33–4.41), respectively. The cumulative incidence of coronary events at 25 years from breast cancer diagnosis were 7.0% in women receiving left-sided RT and 4.4% in women receiving right-sided RT.

Conclusion

Implementing and further developing techniques that lower cardiac doses is important in order to reduce the risk of long-term side effects of adjuvant RT for breast cancer.

Supplementary Information

The online version of this article (10.1007/s00066-022-01927-0) contains supplementary material, which is available to authorized users.

Comprehensive imaging of microcirculatory changes in the foot during endovascular intervention - A technical feasibility study

Chronic limb-threatening ischemia (CLTI) has a major impact on patient's lives and is associated with a heavy health care burden with high morbidity and mortality. Treatment by endovascular intervention is mostly based on macrocirculatory information from angiography and does not consider the microcirculation. Despite successful endovascular intervention according to angiographic criteria, a proportion of patients fail to heal ischemic lesions. This might be due to impaired microvascular perfusion and variations in the supply to different angiosomes. Non-invasive optical techniques for microcirculatory perfusion and oxygen saturation imaging have the potential to provide the interventionist with additional information in real-time, supporting clinical decisions during the intervention. This study presents a novel multimodal imaging system, based on multi-exposure laser speckle contrast imaging and multi-spectral imaging, for continuous use during endovascular intervention. The results during intervention display spatiotemporal changes in the microcirculation compatible with expected physiological reactions during balloon dilation, with initially induced ischemia followed by a restored perfusion, and local administration of a vasodilator inducing hyperemia. We also present perioperative and postoperative follow-up measurements with a pulsatile microcirculation perfusion. Finally, cases of spatial heterogeneity in the observed oxygen saturation and perfusion are discussed. In conclusion, this technical feasibility study shows the potential of the methodology to characterize changes in microcirculation before, during, and after endovascular intervention.

Post-ischemic skin peak oxygen saturation is associated with cardiovascular risk factors: a Swedish cohort study

The objective of this study was to explore the associations between skin microcirculatory function and established cardiovascular risk factors in a large Swedish cohort. As part of the Swedish CArdioPulmonary bioImage Study (SCAPIS), microcirculatory data were acquired at Linköping University hospital, Linköping, Sweden during 2016-2017. The subjects, aged 50-64 years, were randomly selected from the national population register. Microcirculatory reactivity was assessed using a 5-min arterial occlusion-release protocol. Comprehensive skin microcirculatory data were continuously acquired by using a fiberoptic probe placed on the lower right arm. After exclusion of missing data (208), 1557 subjects were remaining. Among the parameters, skin microcirculatory peak oxygen saturation after occlusion release, had the strongest relationship to the cardiovascular risk factors. The linear associations between peak oxygen saturation and cardiovascular risk factors were analyzed adjusted for age and sex. We found a negative association with peak oxygen saturation (standardized regression coefficient) for blood pressure (systolic -0.05 (95% CI: -0.10;-0.003) and diastolic -0.05 (-0.10; -0.003)), BMI -0.18 (-0.23; -0.13), waist circumference (males -0.20 (-0.32; -0.16), females -0.18 (-0.25; -0.11)), prevalent diabetes -0.31 (-0.49; -0.12), hypertension -0.30 (-0.42; -0.18), dyslipidemia -0.24 (-0.40; -0.09), fasting glucose level -0.06 (-0.12; -0.01), HbA1c -0.07 (-0.12; -0.02), triglyceride level -0.09 (-0.14; -0.04), hsCRP -0.12 (-0.17; -0.07), and current smoker versus never smoked -0.50 (-0.67; -0.34). A positive association with peak oxygen saturation was found for cholesterol level 0.05 (0.005; 0.11) and HDL 0.11 (0.06; 0.17). This is the first study showing that post-ischemic skin microvascular peak oxygen saturation is associated with virtually all established cardiovascular risk factors in a population-based middle-aged cohort.

Vasomotion analysis of speed resolved perfusion, oxygen saturation, red blood cell tissue fraction, and vessel diameter: Novel microvascular perspectives

BACKGROUND

Vasomotion is the spontaneous oscillation in vascular tone in the microcirculation and is believed to be a physiological mechanism facilitating the transport of blood gases and nutrients to and from tissues. So far, Laser Doppler flowmetry has constituted the gold standard for in vivo vasomotion analysis.

MATERIALS AND METHODS

We applied vasomotion analysis to speed-resolved perfusion, oxygen saturation, red blood cell tissue (RBC) tissue fraction, and average vessel diameter from five healthy individuals at rest measured by the newly developed Periflux 6000 EPOS system over 10 minutes. Magnitude scalogram and the time-averaged wavelet spectra were divided into frequency intervals reflecting endothelial, neurogenic, myogenic, respiratory, and cardiac function.

RESULTS

Recurrent high-intensity periods of the myogenic, neurogenic, and endothelial frequency intervals were found. The neurogenic activity was considerably more pronounced for the oxygen saturation, RBC tissue fraction, and vessel diameter signals, than for the perfusion signals. In a correlation analysis we found that changes in perfusion in the myogenic, neurogenic, and endothelial frequency intervals precede changes in the other signals. Furthermore, changes in average vessel diameter were in general negatively correlated to the other signals in the same frequency intervals, indicating the importance of capillary recruitment.

CONCLUSION

We conclude that vasomotion can be observed in signals reflecting speed resolved perfusion, oxygen saturation, RBC tissue fraction, and vessel diameter. The new parameters enable new aspects of the microcirculation to be observed.

Machine learning for direct oxygen saturation and hemoglobin concentration assessment using diffuse reflectance spectroscopy

SIGNIFICANCE

Diffuse reflectance spectroscopy (DRS) is frequently used to assess oxygen saturation and hemoglobin concentration in living tissue. Methods solving the inverse problem may include time-consuming nonlinear optimization or artificial neural networks (ANN) determining the absorption coefficient one wavelength at a time.

AIM

To present an ANN-based method that directly outputs the oxygen saturation and the hemoglobin concentration using the shape of the measured spectra as input.

APPROACH

A probe-based DRS setup with dual source-detector separations in the visible wavelength range was used. ANNs were trained on spectra generated from a three-layer tissue model with oxygen saturation and hemoglobin concentration as target.

RESULTS

Modeled evaluation data with realistic measurement noise showed an absolute root-mean-square (RMS) deviation of 5.1% units for oxygen saturation estimation. The relative RMS deviation for hemoglobin concentration was 13%. This accuracy is at least twice as good as our previous nonlinear optimization method. On blood-intralipid phantoms, the RMS deviation from the oxygen saturation derived from partial oxygen pressure measurements was 5.3% and 1.6% in two separate measurement series. Results during brachial occlusion showed expected patterns.

CONCLUSIONS

The presented method, directly assessing oxygen saturation and hemoglobin concentration, is fast, accurate, and robust to noise.

Real-time video-rate perfusion imaging using multi-exposure laser speckle contrast imaging and machine learning

SIGNIFICANCE

Multi-exposure laser speckle contrast imaging (MELSCI) estimates microcirculatory blood perfusion more accurately than single-exposure LSCI. However, the technique has been hampered by technical limitations due to massive data throughput requirements and nonlinear inverse search algorithms, limiting it to an offline technique where data must be postprocessed.

AIM

To present an MELSCI system capable of continuous acquisition and processing of MELSCI data, enabling real-time video-rate perfusion imaging with high accuracy.

APPROACH

The MELSCI algorithm was implemented in programmable hardware (field programmable gate array) closely interfaced to a high-speed CMOS sensor for real-time calculation. Perfusion images were estimated in real-time from the MELSCI data using an artificial neural network trained on simulated data. The MELSCI perfusion was compared to two existing single-exposure metrics both quantitatively in a controlled phantom experiment and qualitatively in vivo.

RESULTS

The MELSCI perfusion shows higher signal dynamics compared to both single-exposure metrics, both spatially and temporally where heartbeat-related variations are resolved in much greater detail. The MELSCI perfusion is less susceptible to measurement noise and is more linear with respect to laser Doppler perfusion in the phantom experiment (R2  =  0.992).

CONCLUSIONS

The presented MELSCI system allows for real-time acquisition and calculation of high-quality perfusion at 15.6 frames per second.

Normative data and the influence of age and sex on microcirculatory function in a middle-aged cohort: results from the SCAPIS study

The objective of this study was to assess normative values for comprehensive forearm skin microcirculatory function: oxygen saturation, tissue fraction of red blood cells (RBCs), and speed-resolved perfusion. Furthermore, to examine the influence of age and sex on microcirculatory function. Measurements were performed using a noninvasive probe-based system, including diffuse reflectance spectroscopy and laser-Doppler flowmetry, yielding output data in absolute units. The study was conducted within the Swedish CArdioPulmonary BioImage Study (SCAPIS) and included 1,765 men and women aged 50-65 yr from the Linköping general population. Normative values are given at baseline, at the end of a 5-min occlusion of the brachial artery and during hyperemia after occlusion release. We found a consistent age distribution, in which the oldest individuals had the lowest peak oxygen saturation (P < 0.001) and the highest baseline low-speed perfusion (P < 0.001). Women had higher peak oxygen saturation (P < 0.001), lower RBC tissue fraction, in general (P < 0.001), lower baseline perfusion in all speed regions (P = 0.01), and lower peak high-speed perfusion at hyperemia (P < 0.001). The normative data can be used as reference values in future studies of disease-specific populations. The results show that age and sex are important aspects to consider in studies of microvascular function. Women and younger age were factors associated with higher peak oxygen saturation after ischemia. This is a novel parameter that reflects overall microcirculatory function associated with vascular dilation capacity.NEW & NOTEWORTHY This study expands experimental microcirculatory research to clinical use by providing normative values on microcirculatory function in a large population-based cohort. Women and younger age were factors associated with higher peak oxygen saturation after ischemia, which implies that age and sex are important aspects to consider in studies of microvascular function. This study is the first step toward using microcirculatory assessment as a tool to improve diagnosis, prognosis, and treatment in disease-specific populations.

Heart Rate Variability for Driver Sleepiness Classification in Real Road Driving Conditions

Approximately 20-30% of all road fatalities are related to driver sleepiness. A long-lasting goal in driver state research has therefore been to develop a robust sleepiness detection system. Since the alertness level is reflected in autonomous nervous system activity, it has been suggested that various heart rate variability (HRV) metrics can be used as features for driver sleepiness classification. Since the heart rate is modulated by many different factors, and not just by sleepiness, it is relevant to question the high driver sleepiness classification accuracies that have occasionally been presented in the literature. The main objective of this paper is thus to test how well a sleepiness classification system based on HRV features really is. A unique data set with 86 drivers, obtained while driving on real roads in real traffic, both in alert and sleep deprived conditions, was used to train and test a support vector machine (SVM) classifier. Subjective ratings based on the Karolinska sleepiness scale (KSS) was used as ground truth to divide the data into three classes (alert, somewhat sleepy and severely sleepy). Even though nearly all the 24 investigated HRV metrics showed significant differences between sleepiness levels, the SVM results only reached a mean accuracy of 61 %, with the worst results originating from the severely sleepy cases. In summary, the high classification performance that may arise in studies with high experimental control could not be replicated under realistic driving conditions. Future works should focus on how various confounding factors should be accounted for when using HRV based metrics as input to a driver sleepiness detection system.

Validation of speed-resolved laser Doppler perfusion in a multimodal optical system using a blood-flow phantom

The PeriFlux 6000 EPOS system combines diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF) for the assessment of oxygen saturation (expressed in percentage), red blood cell (RBC) tissue fraction (expressed as volume fraction, %RBC), and perfusion (%RBC  ×  mm  /  s) in the microcirculation. It also allows the possibility of separating the perfusion into three speed regions (0 to 1, 1 to 10, and >10  mm  /  s). We evaluate the speed-resolved perfusion components, i.e., the relative amount of perfusion within each speed region, using a blood-flow phantom. Human blood was pumped through microtubes with an inner diameter of 0.15 mm. Measured DRS and LDF spectra were compared to Monte Carlo-simulated spectra in an optimization routine, giving the best-fit parameters describing the measured spectra. The root-mean-square error for each of the three speed components (0 to 1, 1 to 10, and >10  mm  /  s, respectively) when describing the blood-flow speed in the microtubes was 2.9%, 8.1%, and 7.7%. The presented results show that the system can accurately discriminate blood perfusion originating from different blood-flow speeds, which may enable improved measurement of healthy and dysfunctional microcirculatory flow.

Abstract P2-08-05: Withdrawn

This abstract was withdrawn by the authors.

Citation Format: Jaraj D, Ahlgren J, Arnesson L-G, Einbeigi Z, Höijer J, Klintman M, Malmström P, Vikhe Patil E, Sund M, Fredriksson I, Bergh J, Pettersson A. Withdrawn [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P2-08-05.

Machine learning in multiexposure laser speckle contrast imaging can replace conventional laser Doppler flowmetry

Laser speckle contrast imaging (LSCI) enables video rate imaging of blood flow. However, its relation to tissue blood perfusion is nonlinear and depends strongly on exposure time. By contrast, the perfusion estimate from the slower laser Doppler flowmetry (LDF) technique has a relationship to blood perfusion that is better understood. Multiexposure LSCI (MELSCI) enables a perfusion estimate closer to the actual perfusion than that using a single exposure time. We present and evaluate a method that utilizes contrasts from seven exposure times between 1 and 64 ms to calculate a perfusion estimate that resembles the perfusion estimate from LDF. The method is based on artificial neural networks (ANN) for fast and accurate processing of MELSCI contrasts to perfusion. The networks are trained using modeling of Doppler histograms and speckle contrasts from tissue models. The importance of accounting for noise is demonstrated. Results show that by using ANN, MELSCI data can be processed to LDF perfusion with high accuracy, with a correlation coefficient R  =  1.000 for noise-free data, R  =  0.993 when a moderate degree of noise is present, and R  =  0.995 for in vivo data from an occlusion-release experiment.

In vivo characterization of light scattering properties of human skin in the 475- to 850-nm wavelength range in a Swedish cohort

We have determined in vivo optical scattering properties of normal human skin in 1734 subjects, mostly with fair skin type, within the Swedish CArdioPulmonary bioImage Study. The measurements were performed with a noninvasive system, integrating spatially resolved diffuse reflectance spectroscopy and laser Doppler flowmetry. Data were analyzed with an inverse Monte Carlo algorithm, accounting for both scattering, geometrical, and absorbing properties of the tissue. The reduced scattering coefficient was found to decrease from 3.16 ± 0.72 to 1.13 ± 0.27 mm-1 (mean ± SD) in the 475- to 850-nm wavelength range. There was a negative correlation between the reduced scattering coefficient and age, and a significant difference between men and women in the reduced scattering coefficient as well as in the fraction of small scattering particles. This large study on tissue scattering with mean values and normal variation can serve as a reference when designing diagnostic techniques or when evaluating the effect of therapeutic optical systems.

Abstract P4-09-08: A targeted breast cancer radiosensitivity gene expression panel

Background: A majority of patients with early breast cancer is operated with breast conserving surgery (BCS) and adjuvant radiotherapy (RT) is administered to prevent ipsilateral breast tumor recurrence (IBTR), including a new ipsilateral cancer. The EBCTCG meta-analysis showed a majority of patients treated with surgery only to be recurrence free at 10 years, and more than 10% to suffer an IBTR despite RT, thus implying considerable over- and under treatment. A wide range of prognosticators, including multigene tests, are well established, but we lack predictive factors for RT, which is the aim in the present study.

Patients and methods: Fresh frozen tissue from 340 patients operated with BCS with or without RT and with or without IBTR was collected (without IBTR N=196, with IBTR n=144). Patients were stratified according to estrogen receptor (ER) status and RT, and divided into a training cohort (N=172) and a validation cohort (N=168). The training cohort was analyzed with whole transcriptome analysis (Illumina HT12 v4) and top discriminating genes for IBTR (N=155) were selected based on a random forest machine learning algorithm with recursive feature elimination and cross-validation. Further, genes described in the literature as associated with radioresistance were included in the panel to a total of 248 genes. A custom nCounter (Nanostring Technologies) gene expression panel was designed and both the training and validation cohorts were analyzed with the custom panel. Single-sample classifiers using a k-top scoring pairs algorithm were trained in the training cohort and validated in the validation cohort. Area under the curve (AUC) with a receiver operator characteristics (ROC) analysis were calculated and p-values were calculated with a log-rank test. All calculations were done using the R statistical environment.

Results: Our classifiers were prognostic for IBTR in the validation cohort among ER+ patients given RT (AUC 0.67, p=0.005), ER+ patients not given RT (AUC=0.89, p=0.015) and ER- patients given RT (AUC=0.78, p&lt;0.001), while the number of ER- patients not given RT was too small for subgroup analysis (N=4). We also created a sequential algorithm were a first classifier was applied to test the risk of IBTR without RT. If low, the tumor was classified as “surgery only”. If classified as high, a second classifier was applied to test the risk of recurrence when given RT. If the risk was predicted low after RT, the tumor was classified as “radiosensitive”. If high, the tumor was classified as “radioresistant”. Among ER+ patients in the validation cohort, the “radiosensitive” tumors had an excellent effect of RT (p&lt;0.001), the “radioresistant” had no effect of RT (p=0.4) and a very high risk of recurrence (55% at 10 years). The tumors predicted as “surgery only” had no effect of RT (p=0.4), and a lower risk of recurrence than the “radioresistant” patients (25% at 10 years).

Conclusions: Our targeted radiosensitivity gene expression panel could identify patients of high or low risk of LR, with or without RT. The most promising was however that it seems as the panel could be used as a predictive marker, i.e., finding patients that do, or do not, respond to RT. Further refinement and testing of the panel and models is ongoing.

Citation Format: Sjöström M, Staaf J, Edén P, Wärnberg F, Bergh J, Malmström P, Fernö M, Niméus E, Fredriksson I. A targeted breast cancer radiosensitivity gene expression panel [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P4-09-08.

A 15.6 frames per second 1-megapixel multiple exposure laser speckle contrast imaging setup

A multiple exposure laser speckle contrast imaging (MELSCI) setup for visualizing blood perfusion was developed using a field programmable gate array (FPGA), connected to a 1000 frames per second (fps) 1-megapixel camera sensor. Multiple exposure time images at 1, 2, 4, 8, 16, 32 and 64 milliseconds were calculated by cumulative summation of 64 consecutive snapshot images. The local contrast was calculated for all exposure times using regions of 4 × 4 pixels. Averaging of multiple contrast images from the 64-millisecond acquisition was done to improve the signal-to-noise ratio. The results show that with an effective implementation of the algorithm on an FPGA, contrast images at all exposure times can be calculated in only 28 milliseconds. The algorithm was applied to data recorded during a 5 minutes finger occlusion. Expected contrast changes were found during occlusion and the following hyperemia in the occluded finger, while unprovoked fingers showed constant contrast during the experiment. The developed setup is capable of massive data processing on an FPGA that enables processing of MELSCI data in 15.6 fps (1000/64 milliseconds). It also leads to improved frame rates, enhanced image quality and enables the calculation of improved microcirculatory perfusion estimates compared to single exposure time systems.

Evaluation of a pointwise microcirculation assessment method using liquid and multilayered tissue simulating phantoms

A fiber-optic probe-based instrument, designed for assessment of parameters related to microcirculation, red blood cell tissue fraction (fRBC), oxygen saturation (SO2), and speed resolved perfusion, has been evaluated using state-of-the-art tissue phantoms. The probe integrates diffuse reflectance spectroscopy (DRS) at two source-detector separations and laser Doppler flowmetry, using an inverse Monte Carlo method for identifying the parameters of a multilayered tissue model. Here, we characterize the accuracy of the DRS aspect of the instrument using (1) liquid blood phantoms containing yeast and (2) epidermis-dermis mimicking solid-layered phantoms fabricated from polydimethylsiloxane, titanium oxide, hemoglobin, and coffee. The root-mean-square (RMS) deviations for fRBC for the two liquid phantoms were 11% and 5.3%, respectively, and 11% for the solid phantoms with highest hemoglobin signatures. The RMS deviation for SO2 was 5.2% and 2.9%, respectively, for the liquid phantoms, and 2.9% for the solid phantoms. RMS deviation for the reduced scattering coefficient (μs'), for the solid phantoms was 15% (475 to 850 nm). For the liquid phantoms, the RMS deviation in average vessel diameter (D) was 1  μm. In conclusion, the skin microcirculation parameters fRBC and SO2, as well as, μs' and D are estimated with reasonable accuracy.

Vessel packaging effect in laser speckle contrast imaging and laser Doppler imaging

Laser speckle-based techniques are frequently used to assess microcirculatory blood flow. Perfusion estimates are calculated either by analyzing the speckle fluctuations over time as in laser Doppler flowmetry (LDF), or by analyzing the speckle contrast as in laser speckle contrast imaging (LSCI). The perfusion estimates depend on the amount of blood and its speed distribution. However, the perfusion estimates are commonly given in arbitrary units as they are nonlinear and depend on the magnitude and the spatial distribution of the optical properties in the tissue under investigation. We describe how the spatial confinement of blood to vessels, called the vessel packaging effect, can be modeled in LDF and LSCI, which affect the Doppler power spectra and speckle contrast, and the underlying bio-optical mechanisms for these effects. As an example, the perfusion estimate is reduced by 25% for LDF and often more than 50% for LSCI when blood is located in vessels with an average diameter of 40  μm, instead of being homogeneously distributed within the tissue. This significant effect can be compensated for only with knowledge of the average diameter of the vessels in the tissue.

Risk of death from breast cancer after treatment for ductal carcinoma in situ

Background

Studies to date have failed to demonstrate any survival benefit from preventing local recurrence after treatment for ductal breast carcinoma in situ (DCIS). Patient- and tumour-related risk factors for death from breast cancer in women with a primary DCIS were analysed here in a large case–control study.

Methods

A nested case–control study was conducted in a population-based cohort of women with primary DCIS between 1992 and 2012. Women who later died from breast cancer were identified. Four controls per case were selected randomly by incidence density sampling. Medical records and pathology reports were retrieved. Conditional logistic regression was used to calculate odds ratios (ORs) and 95 per cent confidence intervals for risk of death from breast cancer.

Results

From a cohort of 6964 women, 96 who died from breast cancer were identified and these were compared with a group of 318 controls. Tumour size over 25 mm or multifocal DCIS (OR 2·55, 95 per cent c.i. 1·53 to 4·25), a positive or uncertain margin status (OR 3·91, 1·59 to 9·61) and detection outside the screening programme (OR 2·12, 1·16 to 3·86) increased the risk of death from breast cancer. The risks were not affected by age or type of treatment. In the multivariable analysis, tumour size (OR 1·95, 1·06 to 3·67) and margin status (OR 2·69, 1·15 to 7·11) remained significant.

Conclusion

In the present study, large tumour size and positive or uncertain margin status were associated with a higher risk of death from breast cancer after treatment for primary DCIS. More extensive treatment was not associated with lower risk, which may be due to confounding by indication, or indicate that some DCIS has an inherent potential for metastatic spread.

Skin microvascular endothelial dysfunction is associated with type 2 diabetes independently of microalbuminuria and arterial stiffness

Skin and kidney microvascular functions may be affected independently in diabetes mellitus. We investigated skin microcirculatory function in 79 subjects with diabetes type 2, where 41 had microalbuminuria and 38 not, and in 41 age-matched controls. The oxygen saturation, fraction of red blood cells and speed-resolved microcirculatory perfusion (% red blood cells × mm/s) divided into three speed regions: 0-1, 1-10 and above 10 mm/s, were assessed during baseline and after local heating of the foot with a new device integrating diffuse reflectance spectroscopy and laser Doppler flowmetry. Arterial stiffness was assessed as carotid-femoral pulse wave velocity. Subjects with diabetes and microalbuminuria had significantly higher carotid-femoral pulse wave velocity compared to subjects without microalbuminuria and to controls. The perfusion for speeds 0-1 mm/s and red blood cell tissue fraction were reduced in subjects with diabetes at baseline and after heating, independent of microalbuminuria. These parameters were correlated to HbA1c. In conclusion, the reduced nutritive perfusion and red blood cell tissue fraction in type 2 diabetes were related to long-term glucose control but independent of microvascular changes in the kidneys and large-vessel stiffness. This may be due to different pathogenic pathways in the development of nephropathy, large-vessel stiffness and cutaneous microvascular impairment.

Abstract P3-17-03: DCIS and the risk of breast cancer death - A case control study

Introduction

The risk of breast cancer death after a primary ductal carcinoma in situ (DCIS) is less than 2 % after 10 years. Whereas in situ recurrences do not influence survival, a 17-fold elevated risk of breast cancer specific mortality has been shown for invasive recurrences. Adjuvant radiotherapy (RT) effectively reduces recurrences after breast conserving surgery (BCS) for DCIS, but no studies have been able to demonstrate a survival benefit from adjuvant RT treatment or from choosing mastectomy instead of BCS. Here patient and tumour related risk factors for breast cancer death in women with a pure primary DCIS were studied.

Patients and methods

Women registered with a primary DCIS, between 1992-2012 in three of Sweden´s health care regions with a population of approximately 5.2 million, were enrolled in a nested case-control study. Out of 6,964 women with DCIS, 96 patients who later died from breast cancer were identified. Four controls per case (n=318) were randomly selected by incidence density sampling. We retrieved medical records and pathology reports and calculated OR with 95% CIs for various variables using conditional logistic regression.

Results

Of the 96 cases, 10 patients developed distant metastasis without a known local recurrence. In 56 patients death was preceded by an invasive ipsilateral recurrence and in 3 patients by a recurrent ipsilateral DCIS. Seven patients had invasive breast events in both the ipsilateral and the contralateral breast. Seventeen patients had contralateral invasive breast cancer and 3 patients contralateral DCIS.

Multifocality and tumour size over 25mm (OR 2.6 (1.6 to 4.2)), positive or uncertain margin status (OR 2.8 (1.6 to 4.9)) and detection outside screening (OR 2.1 (1.2 to 3.9)) increased the risk of breast cancer death in univariate analysis, when adjusted for age and year of diagnosis. Suspicion of micro-invasion and nuclear grade 3 was associated with a nonsignificant increased risk, OR 1.8 (0.6 to 5.0) and 2.6 (0.9-6.5), respectively. The risk was not affected by age or treatment. Tumour size and margin status remained significant in the multivariable analysis, when adjusted for treatment and for contralateral breast cancer (OR 2.0 (1.2 to 3.7)).

Conclusion

In the present study, large tumours and positive or uncertain margin status were significant risk factors for later breast cancer death after a primary DCIS. More extensive treatment was not related to a lower risk. The significance of tumour biology and nuclear grade will be further examined and evaluated.

Citation Format: Wadsten C, Garmo H, Fredriksson I, Sund M, Wärnberg F. DCIS and the risk of breast cancer death - A case control study [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P3-17-03.

On the equivalence and differences between laser Doppler flowmetry and laser speckle contrast analysis

Laser Doppler flowmetry (LDF) and laser speckle contrast analysis (LASCA) both utilize the spatiotemporal properties of laser speckle patterns to assess microcirculatory blood flow in tissue. Although the techniques analyze the speckle pattern differently, there is a close relationship between them. We present a theoretical overview describing how the LDF power spectrum and the LASCA contrast can be calculated from each other, and how both these can be calculated from an optical Doppler spectrum containing various degrees of Doppler shifted light. The theoretical relationships are further demonstrated using time-resolved speckle simulations. A wide range of Monte Carlo simulated tissue models is then used to show how perfusion estimates for LDF and LASCA are affected by changes in blood concentration and speed distribution, as well as by geometrical and optical properties. We conclude that perfusion estimates from conventional single exposure time LASCA are in general more sensitive to changes in optical and geometrical properties and are less accurate in the prediction of real perfusion changes, especially speed changes. Since there is a theoretical one-to-one relationship between Doppler power spectrum and contrast, one can conclude that those drawbacks with the LASCA technique can be overcome using a multiple exposure time setup.

Abstract P6-10-08: Pregnancy associated breast cancer (PABC); no evidence of patients' or doctors' delays

Background:

A small, but not negligible, proportion of breast cancers in young women are detected in association with childbearing. While pregnancy usually is a period of intense medical observation, signs and symptoms of a malignancy may be overlooked or misinterpreted as pregnancy-related, resulting in diagnostic and treatment delays. Also, a delayed diagnosis in pregnant women has been suggested as a reason for the more advanced disease and poorer outcomes in women with pregnancy-associated breast cancer.

Material and Methods:

For the purpose of the present study, pregnancy-associated breast cancer (PABC) was defined as an invasive breast tumor diagnosed during pregnancy and up to two years post-delivery (non-PABC cases were diagnosed outside this time window, or nulliparous). Based on a systematic review of medical records for women aged 15-44 years at diagnosis with PABC and non-PABC identified in Swedish health care registers, chart information was retrieved by trained nurses for a total of 570 women (285 PABC women and 285 age and hospital matched non-PABC women) treated at 11 hospitals across Sweden between 1992 and 2009. Median waiting times from initial signs or symptoms in days to start of treatment, and time periods within, were computed using the Kaplan-Meier method and compared using the logrank test for the Kaplan-Meier curves. Dates on first symptoms were available for 122 matched PABC/non-PABC pairs, in total 244 patients. Full dates to assess and compare times between first health care contact – diagnosis – start of treatment, were available for 246 PABC/non-PABC pairs, in total 492 women.

Objective:

To examine and compare lengths of several defined waiting times within the time period from initial symptoms to start of treatment in women diagnosed with PABC and non-PABC.

Results:

Patient delay-time between first symptom and first point of contact with health care provider.

Median time between first symptoms and first contact with health care was 36 days and 45 days for women with PABC and non-PABC, respectively (logrank test p-value 0.48).

Time between first health care contact and diagnosis

Median time between first contact and diagnosis was 7 days for both PABC and non-PABC women (logrank test p-value 0.16).

Time between diagnosis and start of treatment

The median waiting time from date of diagnosis to initiation of treatment was shorter in women with PABC (22 days) compared to non-PABC women (26 days) (logrank test p-value 0.14).

Time between first contact and start of treatment

The median delay of start of treatment from first contact with a health care provider was 34 days in PABC women and 37 days in non-PABC women. (logrank test p-value 0.14).

Conclusions:

Patients' delay and the time between first contact with health care and start of treatment was shorter in women with PABC compared to non-PABC. Taken together, the present results do not support the notion that diagnostic and treatment delays are more common in women diagnosed with breast cancer during or shortly after pregnancy.

Citation Format: Lambe M, Fredriksson I, Johansson ALV. Pregnancy associated breast cancer (PABC); no evidence of patients' or doctors' delays. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P6-10-08.

Oxygen saturation, red blood cell tissue fraction and speed resolved perfusion - A new optical method for microcirculatory assessment

We have developed a new fiber-optic system that combines diffuse reflectance spectroscopy (DRS) and laser Doppler Flowmetry (LDF) for a multi-modal assessment of the microcirculation. Quantitative data is achieved with an inverse Monte Carlo algorithm based on an individually adaptive skin model. The output parameters are calculated from the model and given in absolute units: hemoglobin oxygen saturation (%), red blood cell (RBC) tissue fraction (%), and the speed resolved RBC perfusion separated into three speed regions; 0-1mm/s, 1-10mm/s and above 10mm/s (% mm/s). The aim was to explore microcirculatory parameters using the new optical method, integrating DRS and LDF in a joint skin model, during local heating of the dorsal foot and venous and arterial occlusion of the forearm in 23 healthy subjects (age 20-28years). There were differences in the three speed regions in regard to blood flow changes due to local heating, where perfusion for high speeds increased the most. There was also a high correlation between changes in oxygenation and changes in perfusion for higher speeds. Oxygen saturation at baseline was 44% on foot, increasing to 83% at plateau after heating. The larger increase in perfusion for higher speeds than for lower speeds together with the oxygenation increase during thermal provocation, shows a local thermoregulatory blood flow in presumably arteriolar dermal vessels. In conclusion, there are improved possibilities to assess microcirculation using integrated DRS and LDF in a joint skin model by enabling both oxygenation and speed resolved blood flow assessment simultaneously and in the same skin site. Output parameters in absolute units may also yield new insights about the microcirculatory system.

Microcirculation assessment using an individualized model for diffuse reflectance spectroscopy and conventional laser Doppler flowmetry

Microvascular assessment would benefit from co-registration of blood flow and hemoglobin oxygenation dynamics during stimulus response tests. We used a fiber-optic probe for simultaneous recording of white light diffuse reflectance (DRS; 475-850 nm) and laser Doppler flowmetry (LDF; 780 nm) spectra at two source-detector distances (0.4 and 1.2 mm). An inverse Monte Carlo algorithm, based on a multiparameter three-layer adaptive skin model, was used for analyzing DRS data. LDF spectra were conventionally processed for perfusion. The system was evaluated on volar forearm recordings of 33 healthy subjects during a 5-min systolic occlusion protocol. The calibration scheme and the optimal adaptive skin model fitted DRS spectra at both distances within 10%. During occlusion, perfusion decreased within 5 s while oxygenation decreased slowly (mean time constant 61 s; dissociation of oxygen from hemoglobin). After occlusion release, perfusion and oxygenation increased within 3 s (inflow of oxygenized blood). The increased perfusion was due to increased blood tissue fraction and speed. The supranormal hemoglobin oxygenation indicates a blood flow in excess of metabolic demands. In conclusion, by integrating DRS and LDF in a fiber-optic probe, a powerful tool for assessment of blood flow and oxygenation in the same microvascular bed has been presented.

Inverse Monte Carlo in a multilayered tissue model: merging diffuse reflectance spectroscopy and laser Doppler flowmetry

The tissue fraction of red blood cells (RBCs) and their oxygenation and speed-resolved perfusion are estimated in absolute units by combining diffuse reflectance spectroscopy (DRS) and laser Doppler flowmetry (LDF). The DRS spectra (450 to 850 nm) are assessed at two source-detector separations (0.4 and 1.2 mm), allowing for a relative calibration routine, whereas LDF spectra are assessed at 1.2 mm in the same fiber-optic probe. Data are analyzed using nonlinear optimization in an inverse Monte Carlo technique by applying an adaptive multilayered tissue model based on geometrical, scattering, and absorbing properties, as well as RBC flow-speed information. Simulations of 250 tissue-like models including up to 2000 individual blood vessels were used to evaluate the method. The absolute root mean square (RMS) deviation between estimated and true oxygenation was 4.1 percentage units, whereas the relative RMS deviations for the RBC tissue fraction and perfusion were 19% and 23%, respectively. Examples of in vivo measurements on forearm and foot during common provocations are presented. The method offers several advantages such as simultaneous quantification of RBC tissue fraction and oxygenation and perfusion from the same, predictable, sampling volume. The perfusion estimate is speed resolved, absolute (% RBC×mm/s), and more accurate due to the combination with DRS.

Inverse Monte Carlo for estimation of scattering and absorption in liquid optical phantoms

A spectroscopic probe with multiple detecting fibers was used for quantifying absorption and scattering in liquid optical phantoms. The phantoms were mixtures of Intralipid and red and blue food dyes. Intensity calibration for the detecting fibers was undertaken using either a microsphere suspension (absolute calibration) or a uniform detector illumination (relative calibration between detectors). Two different scattering phase functions were used in an inverse Monte Carlo algorithm. Data were evaluated for residual spectra (systematic deviations and magnitude) and accuracy in estimation of scattering and absorption. Spectral fitting was improved by allowing for a 10% intensity relaxation in the optimization algorithm. For a multi-detector setup, non-systematic residual spectrum was only found using the more complex Gegenbauer-kernel phase function. However, the choice of phase function did not influence the accuracy in the estimation of absorption and scattering. Similar estimation accuracy as in the multi-detector setup was also obtained using either two relative calibrated detectors or one absolute calibrated detector at a fiber separation of 0.46 mm.

Inverse Monte Carlo method in a multilayered tissue model for diffuse reflectance spectroscopy

Model based data analysis of diffuse reflectance spectroscopy data enables the estimation of optical and structural tissue parameters. The aim of this study was to present an inverse Monte Carlo method based on spectra from two source-detector distances (0.4 and 1.2 mm), using a multilayered tissue model. The tissue model variables include geometrical properties, light scattering properties, tissue chromophores such as melanin and hemoglobin, oxygen saturation and average vessel diameter. The method utilizes a small set of presimulated Monte Carlo data for combinations of different levels of epidermal thickness and tissue scattering. The path length distributions in the different layers are stored and the effect of the other parameters is added in the post-processing. The accuracy of the method was evaluated using Monte Carlo simulations of tissue-like models containing discrete blood vessels, evaluating blood tissue fraction and oxygenation. It was also compared to a homogeneous model. The multilayer model performed better than the homogeneous model and all tissue parameters significantly improved spectral fitting. Recorded in vivo spectra were fitted well at both distances, which we previously found was not possible with a homogeneous model. No absolute intensity calibration is needed and the algorithm is fast enough for real-time processing.

Model-based quantitative laser Doppler flowmetry in skin

Laser Doppler flowmetry (LDF) can be used for assessing the microcirculatory perfusion. However, conventional LDF (cLDF) gives only a relative perfusion estimate for an unknown measurement volume, with no information about the blood flow speed distribution. To overcome these limitations, a model-based analysis method for quantitative LDF (qLDF) is proposed. The method uses inverse Monte Carlo technique with an adaptive three-layer skin model. By analyzing the optimal model where measured and simulated LDF spectra detected at two different source-detector separations match, the absolute microcirculatory perfusion for a specified speed region in a predefined volume is determined. qLDF displayed errors<12% when evaluated using simulations of physiologically relevant variations in the layer structure, in the optical properties of static tissue, and in blood absorption. Inhomogeneous models containing small blood vessels, hair, and sweat glands displayed errors<5%. Evaluation models containing single larger blood vessels displayed significant errors but could be dismissed by residual analysis. In vivo measurements using local heat provocation displayed a higher perfusion increase with qLDF than cLDF, due to nonlinear effects in the latter. The qLDF showed that the perfusion increase occurred due to an increased amount of red blood cells with a speed>1 mm∕s.

Reduced arteriovenous shunting capacity after local heating and redistribution of baseline skin blood flow in type 2 diabetes assessed with velocity-resolved quantitative laser Doppler flowmetry

OBJECTIVE

To compare the microcirculatory velocity distribution in type 2 diabetic patients and nondiabetic control subjects at baseline and after local heating.

RESEARCH DESIGN AND METHODS

The skin blood flow response to local heating (44 degrees C for 20 min) was assessed in 28 diabetic patients and 29 control subjects using a new velocity-resolved quantitative laser Doppler flowmetry technique (qLDF). The qLDF estimates erythrocyte (RBC) perfusion (velocity x concentration), in a physiologically relevant unit (grams RBC per 100 g tissue x millimeters per second) in a fixed output volume, separated into three velocity regions: v <1 mm/s, v 1-10 mm/s, and v >10 mm/s.

RESULTS

The increased blood flow occurs in vessels with a velocity >1 mm/s. A significantly lower response in qLDF total perfusion was found in diabetic patients than in control subjects after heat provocation because of less high-velocity blood flow (v >10 mm/s). The RBC concentration in diabetic patients increased sevenfold for v between 1 and 10 mm/s, and 15-fold for v >10 mm/s, whereas no significant increase was found for v <1 mm/s. The mean velocity increased from 0.94 to 7.3 mm/s in diabetic patients and from 0.83 to 9.7 mm/s in control subjects.

CONCLUSIONS

The perfusion increase occurs in larger shunting vessels and not as an increase in capillary flow. Baseline diabetic patient data indicated a redistribution of flow to higher velocity regions, associated with longer duration of diabetes. A lower perfusion was associated with a higher BMI and a lower toe-to-brachial systolic blood pressure ratio.

Simulation of reflected light intensity changes during navigation and radio-frequency lesioning in the brain

An electrode with adjacent optical fibers for measurements during navigation and radio frequency lesioning in the brain is modeled for Monte Carlo simulations of light transport in brain tissue. Relative reflected light intensity at 780 nm, I780, from this electrode and probes with identical fiber configuration are simulated using the intensity from native white matter as reference. Models are made of homogeneous native and coagulated gray, thalamus, and white matter as well as blood. Dual layer models, including models with a layer of cerebrospinal fluid between the fibers and the brain tissue, are also made. Simulated I780 was 0.16 for gray matter, 0.67 for coagulate gray matter, 0.36 for thalamus, 0.39 for coagulated thalamus, unity for white matter, 0.70 for coagulated white matter, and 0.24 for blood. Thalamic matter is also found to reflect more light than gray matter and less than white matter in clinical studies. In conclusion, the reflected light intensity can be used to differentiate between gray and white matter during navigation. Furthermore, coagulation of light gray tissue, such as the thalamus, might be difficult to detect using I780, but coagulation in darker gray tissue should result in a rapid increase of I780.

Forced detection Monte Carlo algorithms for accelerated blood vessel image simulations

Two forced detection (FD) variance reduction Monte Carlo algorithms for image simulations of tissue-embedded objects with matched refractive index are presented. The principle of the algorithms is to force a fraction of the photon weight to the detector at each and every scattering event. The fractional weight is given by the probability for the photon to reach the detector without further interactions. Two imaging setups are applied to a tissue model including blood vessels, where the FD algorithms produce identical results as traditional brute force simulations, while being accelerated with two orders of magnitude. Extending the methods to include refraction mismatches is discussed.

Forced detection Monte Carlo algorithms for accelerated blood vessel image simulations

Two forced detection (FD) variance reduction Monte Carlo algorithms for image simulations of tissue-embedded objects with matched refractive index are presented. The principle of the algorithms is to force a fraction of the photon weight to the detector at each and every scattering event. The fractional weight is given by the probability for the photon to reach the detector without further interactions. Two imaging setups are applied to a tissue model including blood vessels, where the FD algorithms produce identical results as traditional brute force simulations, while being accelerated with two orders of magnitude. Extending the methods to include refraction mismatches is discussed.

Spectral determination of a two-parametric phase function for polydispersive scattering liquids

A method for determining a two-parametric Gegenbauer-kernel phase function that accurately describes the diffuse reflectance from a polydispersive scattering media at small source-detector separations (0.23 to 1.2 mm), is presented. The method involves spectral collimated transmission measurements, spatially resolved spectral diffuse reflectance (SRDR) measurements, and inverse Monte Carlo technique. Both absolute calibration (using a monodispersive polystyrene microsphere suspension) and relative calibration (eliminating differences between fibers) of SRDR spectra yielded comparable results. When applied to water dilutions of milk, simulated and measured spectra deviated less than 6.5% and 2.5% for the absolute and relative calibration case, respectively, even for the closest fiber separation. Corresponding values for milk including ink as an absorber, were 13.4% and 7.3%.