Measuring blood oxygen saturation along a capillary vessel in human

Relating Macroscopic PET Radiomics Features to Microscopic Tumor Phenotypes Using a Stochastic Mathematical Model of Cellular Metabolism and Proliferation

Cancers can manifest large variations in tumor phenotypes due to genetic and microenvironmental factors, which has motivated the development of quantitative radiomics-based image analysis with the aim to robustly classify tumor phenotypes in vivo. Positron emission tomography (PET) imaging can be particularly helpful in elucidating the metabolic profiles of tumors. However, the relatively low resolution, high noise, and limited PET data availability make it difficult to study the relationship between the microenvironment properties and metabolic tumor phenotype as seen on the images. Most of previously proposed digital PET phantoms of tumors are static, have an over-simplified morphology, and lack the link to cellular biology that ultimately governs the tumor evolution. In this work, we propose a novel method to investigate the relationship between microscopic tumor parameters and PET image characteristics based on the computational simulation of tumor growth. We use a hybrid, multiscale, stochastic mathematical model of cellular metabolism and proliferation to generate simulated cross-sections of tumors in vascularized normal tissue on a microscopic level. The generated longitudinal tumor growth sequences are converted to PET images with realistic resolution and noise. By changing the biological parameters of the model, such as the blood vessel density and conditions for necrosis, distinct tumor phenotypes can be obtained. The simulated cellular maps were compared to real histology slides of SiHa and WiDr xenografts imaged with Hoechst 33342 and pimonidazole. As an example application of the proposed method, we simulated six tumor phenotypes that contain various amounts of hypoxic and necrotic regions induced by a lack of oxygen and glucose, including phenotypes that are distinct on the microscopic level but visually similar in PET images. We computed 22 standardized Haralick texture features for each phenotype, and identified the features that could best discriminate the phenotypes with varying image noise levels. We demonstrated that "cluster shade" and "difference entropy" are the most effective and noise-resilient features for microscopic phenotype discrimination. Longitudinal analysis of the simulated tumor growth showed that radiomics analysis can be beneficial even in small lesions with a diameter of 3.5-4 resolution units, corresponding to 8.7-10.0 mm in modern PET scanners. Certain radiomics features were shown to change non-monotonically with tumor growth, which has implications for feature selection for tracking disease progression and therapy response.

Vascular Aging in the Choroid Plexus: A 7T Ultrasmall Superparamagnetic Iron Oxide (USPIO)-MRI Study

BACKGROUND

The choroid plexus (ChP), a densely vascularized structure, has drawn increasing attention for its involvement in brain homeostasis and waste clearance. While the volumetric changes have been explored in many imaging studies, few studies have investigated the vascular degeneration associated with aging in the ChP.

PURPOSE

To investigate the sub-structural characteristics of the ChP, particularly the vascular compartment using high-resolution 7T imaging enhanced with Ferumoxytol, an ultrasmall super-paramagnetic iron oxide, which greatly increase the susceptibility contrast for vessels.

STUDY TYPE

Prospective.

SUBJECTS

Forty-nine subjects without neurological disorders (age: 21-80 years; 42 ± 17 years; 20 females).

FIELD STRENGTH/SEQUENCE

7-T with 2D and 3D T2* GRE, 3D MPRAGE T1, 2D TSE T2, and 2D FLAIR.

ASSESSMENT

The vascular and stromal compartments of the ChP were segmented using K-means clustering on post-contrast 2D GRE images. Visual and qualitative assessment of ChP vascular characteristics were conducted independently by three observers. Vascular density (Volvessel/VolChP ratio) and susceptibility change (Δχ) induced by Ferumoxytol were analyzed on 3D GRE-derived susceptibility-weighted imaging and quantitative susceptibility mapping, respectively.

STATISTICAL TESTS

Independent t-test, Mann-Whitney U test, and Chi-square test were utilized for group comparisons. The relationship between age and ChP's vascular alterations was examined using Pearson's correlation. Intra-class coefficient was calculated for inter-observer agreement. A P value <0.05 was considered statistically significant.

RESULTS

2D GRE images demonstrated superior contrast and accurate delineation of ChP substructures (ICC = 0.86). Older subjects exhibited a significantly smaller vascular density (16.5 ± 4.34%) and lower Δχ (22.10 ± 12.82 ppb) compared to younger subjects (24.85 ± 6.84% and 34.64 ± 12.69 ppb). Vascular density and mean Δχ within the ChP negatively correlated with age (r = -0.48, and r = -0.45).

DATA CONCLUSION

Ferumoxytol-enhanced 7T images can demonstrate ChP alterations in elderly with decreased vascular density and expansion of nonvascular compartment.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

In vivo structural and functional imaging of human nailbed microvasculature using photoacoustic microscopy

Monitoring microvascular structure and function is of great significance for the diagnosis of many diseases. In this study, we demonstrate the feasibility of OR-PAM to nailbed microcirculation detection as a new, to the best of our knowledge, application scenario in humans. We propose a dual-wavelength optical-resolution photoacoustic microscopy (OR-PAM) with improved local-flexible coupling to image human nailbed microvasculature. Microchip lasers with 532 nm wavelength are employed as the pump sources. The 558 nm laser is generated from the 532 nm laser through the stimulated Raman scattering effect. The flowing water, circulated by a peristaltic pump, maintains the acoustic coupling between the ultrasonic transducer and the sample. These designs improve the sensitivity, practicality, and stability of the OR-PAM system for human in vivo experiments. The imaging of the mouse ear demonstrates the ability of our system to acquire structural and functional information. Then, the system is applied to image human nailbed microvasculature. The imaging results reveal that the superficial capillaries are arranged in a straight sagittal pattern, approximately parallel to the long axis of the finger. The arterial and venular limbs are distinguished according to their oxygen saturation differences. Additionally, the images successfully discover the capillary loops with single or multiple twists, the oxygen release at the end of the capillary loop, and the changes when the nailbed is abnormal.

Vascular-centric mapping of in vivo blood oxygen saturation in preclinical models

Assessing intravascular blood oxygen saturation (SO₂) is crucial for characterizing in vivo microenvironmental changes in preclinical models of injury and disease. However, most conventional optical imaging techniques for mapping in vivo SO₂ assume or compute a single value of the optical path-length in tissue. This is especially detrimental when mapping in vivo SO₂ in experimental disease or wound healing models that are characterized by vascular and tissue remodeling. Therefore, to circumvent this limitation we developed an in vivo SO₂ mapping technique that utilizes hemoglobin-based intrinsic optical signal (IOS) imaging combined with a vascular-centric estimation of optical path-lengths. In vivo arterial and venous SO₂ distributions derived with this approach closely matched those reported in the literature, while those derived using the single path-length (i.e. conventional) approach did not. Moreover, in vivo cerebrovascular SO₂ strongly correlated (R² > 0.7) with changes in systemic SO₂ measured with a pulse oximeter during hypoxia and hyperoxia paradigms. Finally, in a calvarial bone healing model, in vivo SO₂ assessed over four weeks was spatiotemporally correlated with angiogenesis and osteogenesis (R² > 0.6). During the early stages of bone healing (i.e. day 10), angiogenic vessels surrounding the calvarial defect exhibited mean SO₂ that was elevated by10 % (p < 0.05) relative to that observed at a later stage (i.e., day 26), indicative of their role in osteogenesis. These correlations were not evident with the conventional SO₂ mapping approach. The feasibility of our wide field-of-view in vivo SO₂ mapping approach illustrates its potential for characterizing the microvascular environment in applications ranging from tissue engineering to cancer.

Eye-specific 3D modeling of factors influencing oxygen concentration in the lamina cribrosa

Our goal was to identify the factors with the strongest influence on the minimum lamina cribrosa (LC) oxygen concentration as potentially indicative of conditions increasing hypoxia risk. Because direct measurement of LC hemodynamics and oxygenation is not yet possible, we developed 3D eye-specific LC vasculature models. The vasculature of a normal monkey eye was perfusion-labeled post-mortem. Serial cryosections through the optic nerve head were imaged using fluorescence and polarized light microscopy to visualize the vasculature and collagen, respectively. The vasculature within a 450 μm-thick region containing the LC - identified from the collagen, was segmented, skeletonized, and meshed for simulations. Using Monte Carlo sampling, 200 vascular network models were generated with varying vessel diameter, neural tissue oxygen consumption rate, inflow hematocrit, and blood pressures (arteriole, venule, anterior boundary, and posterior boundary). Factors were varied over ranges of baseline ±20% with uniform probability. For each model we first obtained the blood flow, and from this the neural tissue oxygen concentration. ANOVA was used to identify the factors with the strongest influence on the minimum (10th percentile) oxygen concentration in the LC. The three most influential factors were, in ranked order, vessel diameter, neural tissue oxygen consumption rate, and arteriole pressure. There was a strong interaction between vessel diameter and arteriole pressure whereby the impact of one factor was larger when the other factor was small. Our results show that, for the eye analyzed, conditions that reduce vessel diameter, such as vessel compression due to elevated intraocular pressure or gaze-induced tissue deformation, may particularly contribute to decreased LC oxygen concentration. More eyes must be analyzed before generalizing.

Skin Complications of Diabetes Mellitus Revealed by Polarized Hyperspectral Imaging and Machine Learning

Aging and diabetes lead to protein glycation and cause dysfunction of collagen-containing tissues. The accompanying structural and functional changes of collagen significantly contribute to the development of various pathological malformations affecting the skin, blood vessels, and nerves, causing a number of complications, increasing disability risks and threat to life. In fact, no methods of non-invasive assessment of glycation and associated metabolic processes in biotissues or prediction of possible skin complications, e.g., ulcers, currently exist for endocrinologists and clinical diagnosis. In this publication, utilizing emerging photonics-based technology, innovative solutions in machine learning, and definitive physiological characteristics, we introduce a diagnostic approach capable of evaluating the skin complications of diabetes mellitus at the very earlier stage. The results of the feasibility studies, as well as the actual tests on patients with diabetes and healthy volunteers, clearly show the ability of the approach to differentiate diabetic and control groups. Furthermore, the developed in-house polarization-based hyperspectral imaging technique accomplished with the implementation of the artificial neural network provides new horizons in the study and diagnosis of age-related diseases.

Quantification of ocular surface microcirculation by computer assisted video microscopy and diffuse reflectance spectroscopy

In piglets we tested the applicability of digital video microscopy and diffuse reflectance spectroscopy for non-invasive assessments of limbal and bulbar conjunctival microcirculation. A priori we postulated that the metabolic rate is higher in limbal as compared to bulbar conjunctiva, and that this difference is reflected in microvascular structure or function between the two locations. Two study sites, Oslo University Hospital (OUH), Norway and Cleveland Clinic (CC), USA, used the same video microscopy and spectroscopy techniques to record limbal and bulbar microcirculation in sleeping piglets. Recordings were analyzed with custom-made software to quantify functional capillary density, capillary flow velocity and microvascular oxygen saturation in measuring volumes of approximately 0.1 mm³. The functional capillary density was higher in limbus than in bulbar conjunctiva at both study sites (OUH: 18.1 ± 2.9 versus 12.2 ± 2.9 crossings per mm line, p < 0.01; CC: 11.3 ± 3.0 versus 7.1 ± 2.8 crossings per mm line, p < 0.01). Median categorial capillary blood flow velocity was higher in bulbar as compared with limbal recordings (CC: 3 (1-3) versus 1 (0-3), p < 0.01). Conjunctival microvascular oxygen saturation was 88 ± 5.9% in OUH versus 94 ± 7.5% in CC piglets. Non-invasive digital video microscopy and diffuse reflectance spectroscopy can be used to obtain data from conjunctival microcirculation in piglets. Limbal conjunctival microcirculation has a larger capacity for oxygen delivery as compared with bulbar conjunctiva.

Single wavelength measurements of absorption coefficients based on iso-pathlength point

In optical sensing, to reveal the chemical composition of tissues, the main challenge is isolating absorption from scattering. Most techniques use multiple wavelengths, which adds an error due to the optical pathlength differences. We suggest using a unique measurement angle for cylindrical tissues, the iso-pathlength (IPL) point, which depends on tissue geometry only (specifically the effective radius). We present a method for absorption assessment from a single wavelength at multiple measurement angles. The IPL point presented similar optical pathlengths for different tissues, both in simulation and experiments, hence it is optimal. Finally, in vivo measurements validated our proposed method.

High spatiotemporal vessel-specific hemodynamic mapping with multi-echo single-vessel fMRI

High-resolution fMRI enables noninvasive mapping of the hemodynamic responses from individual penetrating vessels in animal brains. Here, a 2D multi-echo single-vessel fMRI (MESV-fMRI) method has been developed to map the fMRI signal from arterioles and venules with a 100 ms sampling rate at multiple echo times (TE, 3-30 ms) and short acquisition windows (<1 ms). The T₂*-weighted signal shows the increased extravascular effect on venule voxels as a function of TE. In contrast, the arteriole voxels show an increased fMRI signal with earlier onset than venules voxels at the short TE (3 ms) with increased blood inflow and volume effects. MESV-fMRI enables vessel-specific T2* mapping and presents T2*-based fMRI time courses with higher contrast-to-noise ratios (CNRs) than the T2*-weighted fMRI signal at a given TE. The vessel-specific T2* mapping also allows semi-quantitative estimation of the oxygen saturation levels (Y) and their changes (ΔY) at a given blood volume fraction upon neuronal activation. The MESV-fMRI method enables vessel-specific T2* measurements with high spatiotemporal resolution for better modeling of the fMRI signal based on the hemodynamic parameters.

Dynamic microcirculation PIPE model for functional neuroimaging, non-neuroimaging, and coherent hemodynamics spectroscopy: blood volume and flow velocity variations, and vascular autoregulation

We present a dynamic microcirculation PIPE model for functional neuroimaging, non-neuroimaging, and coherent hemodynamics spectroscopy. The temporal evolution of the concentration and oxygen saturation of hemoglobin in tissue, comprised of the contributions from the arterioles, capillaries, and venules of microvasculature, is determined by time-resolved hemodynamic and metabolic variations in blood volume, flow velocity, and oxygen consumption with a fluid mechanics treatment. Key parameters regarding microcirculation can be assessed, including the effective blood transit times through the capillaries and the venules, and the rate constant of oxygen release from hemoglobin to tissue. The vascular autoregulation can further be quantified from the relationship between the resolved blood volume and flow velocity variations. The PIPE model shows excellent agreement with the experimental cerebral and cutaneous coherent hemodynamics spectroscopy (CHS) and fMRI-BOLD data. It further identifies the impaired cerebral autoregulation distinctively in hemodialysis patients compared to healthy subjects measured by CHS. This new dynamic microcirculation PIPE model provides a valuable tool for brain and other functional studies with hemodynamic-based techniques. It is instrumental in recovering physiological parameters from analyzing and interpreting the signals measured by hemodynamic-based neuroimaging and non-neuroimaging techniques such as functional near-infrared spectroscopy (fNIRS) and functional magnetic resonance imaging (fMRI) in response to brain activation, physiological challenges, or physical maneuvers.

Hyperspectral imaging of human skin aided by artificial neural networks

We developed a compact, hand-held hyperspectral imaging system for 2D neural network-based visualization of skin chromophores and blood oxygenation. State-of-the-art micro-optic multichannel matrix sensor combined with the tunable Fabry-Perot micro interferometer enables a portable diagnostic device sensitive to the changes of the oxygen saturation as well as the variations of blood volume fraction of human skin. Generalized object-oriented Monte Carlo model is used extensively for the training of an artificial neural network utilized for the hyperspectral image processing. In addition, the results are verified and validated via actual experiments with tissue phantoms and human skin in vivo. The proposed approach enables a tool combining both the speed of an artificial neural network processing and the accuracy and flexibility of advanced Monte Carlo modeling. Finally, the results of the feasibility studies and the experimental tests on biotissue phantoms and healthy volunteers are presented.