Real-time monitoring of the cellular oxygen consumption rate [OCR, attomoles of per cell per s (amol/cell/s)] can characterize the metabolic activity of a tissue, providing physiological information with clinical relevance. For instance, the metabolic rate may provide a basis for comparing normal and pathological states. Prior studies have used other methods such as respirometry of extracted mitochondria ex vivo and in vivo multiphoton microscopy to assess the OCR of tissues.12.–3 This paper focuses on assessing the OCR of the skin’s epidermis. Monitoring the epidermal OCR may characterize normal, diseased, and aged skin, detect malignancy, and monitor the effects of drugs and treatments.
Different values of OCR for various cells and tissues have been reported in the literature. Based on data from in vitro studies, Wagner et al.4 cataloged a wide range of OCR from 1 to , with a median value of 23. Balu et al.1 used multiphoton microscopy in vivo to image NADH fluorescence in the epidermis as a proxy for OCR. They reported an OCR of about in keratinocytes in the basal layer, which equals . They also reported a nonuniformity of OCR across the epidermal thickness, with more superficial keratinocytes having lower OCR.
This paper presents an alternate method using a topical optical fiber probe measurement at the skin surface to assess the blood volume and hemoglobin oxygen saturation before, during, and after manually applied pressure to blanch the superficial vascular plexus that provides oxygen to the epidermis. The refill time after release of pressure specifies the blood flow in the blanched superficial plexus. The magnitude and dynamics of changes in blood content and hemoglobin oxygen saturation yield an estimate of the epidermal OCR in units of amol . The goal of this paper is to report a pilot study that illustrates the method and cites initial values for the in vivo epidermal OCR.
Previous work has used applied pressure to remove blood from a tissue for the purpose of improving the observation of other chromophores and tissue structures. For example, Ermakov and Gellermann detected carotenoids in the skin after using pressure to remove superficial blood.5 Such reports use pressure to remove the blood, which is regarded as an interfering factor. In this paper, we use pressure to identify the superficial vascular plexus, which is easily blanched by pressure. The spectroscopy of this blanchable plexus informs about the OCR of the epidermis.
An optical fiber probe was used to spectrally measure the transport of light through in vivo forearm skin sites (Fig. 1). The probe consisted of a source fiber and a collection fiber with a center-to-center spacing of 1.88 mm (0.400-mm-core-dia fibers). The fibers were embedded in a flat plastic plate (three-dimensional-printed using polylactate) that allowed perpendicular fiber contact with the skin while distributing the pressure applied to the skin. The probe handpiece was a cylinder within which the optical probe assembly moved. As one manually pressed the probe against the skin, the outer cylinder pushed a spring against the inner probe assembly, hence applying pressure by the spring. The displacement of the spring was observable on a ruler attached to the inner assembly, which extended outside of the top of the cylinder to allow easy viewing. A reproducible pressure was applied by pressing the probe until the desired displacement was achieved. Repeated tests of the blanching force registered on a scale. The flat surface of the probe tip distributed the force, hence, the probe could also be placed lightly on the skin to monitor refill while avoiding blanching the superficial plexus of blood.
The source fiber delivered white light. The collection fiber returned light to a spectrometer (QE Pro, Ocean Optics Inc., Dunedin, Florida) which measured the light in the 300 to 1100 nm range. Light in the 500 to 1000 nm range yielded strong signals and was used for the analysis. The spectra were acquired continuously by a digital chart recorder (implemented using a MATLABTM program with a GUI interface). Spectra with light pressure were acquired for about 30 s to establish a baseline. Then, strong pressure was applied for about 10 s, which blanched the superficial blood volume, i.e., the superficial vascular plexus. Then light pressure was resumed, which allowed the superficial plexus to refill with blood.
Spectral analysis used the 500 to 1000 nm range of wavelengths. Light below 500 nm (blue light) is scattered by the epidermis and superficial papillary dermis before encountering significant blood content; therefore, spectra at such short wavelengths require a two-layer model for analysis. Light above 500 nm (green to near-infrared) penetrates into the dermis and samples the superficial plexus of microvasculature at about 200 to depth. Hence, a simple diffusion theory that assumes uniform distribution of blood could be used (see Sec. 2.3).
The analysis of the spectrum at each time point specified the apparent average blood volume fraction () and the mixed arteriovenous oxygen saturation () (see Sec. 2.3). Hence, the time courses of and were specified. The values of and after applying pressure for about 10 s () and after the refilling of the plexus () were noted: = at maximum blanching, = after refill of the plexus upon release of pressure. The data at these time points were used to calculate the OCR.
The experiment was conducted on ventral forearm skin sites: nine measurements on five subjects. For each subject, the probe was placed on the forearm, and pressure was applied then released. Then, the probe was moved to a different site on the forearm and the experiment repeated, eventually measuring three to six sites. Measurements over a visible large blood vessel were avoided. Measurements were approved by the Institutional Internal Review Board at OHSU.
Calculation of OCR
Calculation of the OCR (amol/cell/s) was accomplished by the following:
1. Blood flow
2. Oxygen extraction
3. Number of cells in epidermis per surface area
The factor is the ratio of blood in the plexus versus the apparent average blood deduced by diffusion theory from the observed spectrum. The spectral analysis uses a simple diffusion theory [see Eq. (7)] that assumes a homogeneous skin tissue with an overlying pigmented epidermis. Any observed blood absorption is assumed to be due to blood distributed evenly throughout a semi-infinite homogeneous skin. But the superficial vascular plexus has a higher local blood content () than the average () and is located at about 200 to depth (a thickness assumed here). The spectrally observed average does not equal the true local .
A Monte Carlo simulation (MCML6) was used to simulate the reflectance spectra (15 wavelengths from 500 to 1000 nm) when such a plexus was given different blood volume fractions ( to 0.05). The spectra were then analyzed to specify an apparent , using the same diffusion theory analysis that was applied to the experimental data. There was a linear relationship between and with a slope of 2.5. Hence, when spectral analysis of a measured spectrum specified to be 0.01, the true was 0.025.
The relative sizes of the arterial and venous vascular volumes are an important parameter in the calculation of OCR. The ratio of volumes, , was assumed to be 1.0 in this paper. Equation (2), however, allows for to deviate from 1.0.
The measurements yielded a reflectance spectrum, , a spectrum from a white standard (99% reflectance standard, SpectralonTM, Labsphere Inc., New Hampshire) with the probe held 1 cm above the standard, , and a spectrum with the probe held in air within a dark room, . The dark spectrum was subtracted from the other spectra, then a ratio was taken
Figure 2(a) shows the raw reflectance spectra in units of (counts): , , and applied pressure. Figure 2(b) shows the reflectance applied pressure. The factor equals 2.5 based on the average value fit for all acquired spectra on the five subjects.
The spectrum was matched by least squares fitting using a simple light transport expression based on diffusion theory [Eq. (7)]. A predicted spectrum was determined by specifying the fitting parameters: blood content (), oxygen saturation (), epidermal melanin (), scattering (), and a calibration factor (). These parameters are defined below. The parameter scales the magnitude of the reduced scattering coefficient, , whose spectral shape is described below. The expression for (counts) was7 epidermal.thickness ; is the absorption coefficient () 8 spectrum of deoxygenated whole blood;8 scattering coefficient () of skin, based on literature;9 at , which scales the shape of such that; = Rayleigh plus Mie scattering;9 = local reflectance () (diffusion theory, Farrell et al.10); center-to-center distance between source and collection fibers. , ratio of the approximate refractive indices, .
Hence, the spectral fitting yielded values for , and . After an initial analysis of all spectra on the five subjects, a common value of was determined to be 2.5. The same was used to analyze all the spectra. Subsequent analysis fit for , and only. The water content was assumed to be for all subjects. Figure 3 shows the spectral fitting of the nine spectra from five subjects.
First, the calculation of OCR is illustrated by applying the calculation to one of the pressure-release data from one subject in Fig. 4. The calculation steps are listed in Table 1. The mean values of OCR for the six pressure releases on that subject yielded an , which is close to the value of reported by Balu et al.1
Example calculation from the fourth pressure release in Fig. 3.
|measured while blanched||0.00121|
|measured after refill||0.00313|
|measured after refill||0.655|
|from Monte Carlo sims.||2.5|
|= conc. of HGb in whole blood|
Second, the analysis of the nine skin sites on five subjects is summarized in Table 2. The means and standard deviations of the OCR on each of the nine skin sites are listed. The average coefficient of variation for measurements on a site was 34%, which is site-to-site variation on the forearm during one sitting of a subject. The range of mean OCR values extended from 0.9 to . The average mean OCR value was ( skin sites). There was significant sitting-to-sitting variation on subjects 1, 2, and 5, who were measured more than once.
Summary of OCR for nine sites on five subjects.
|Site||Subject||Mean±std (n pressure releases)|
Figure 5 places the OCR data of this report on a histogram drawn using the tabulated data in the review of OCR values for various tissues and cell cultures by Wagner et al.4 A couple of the skin site data are on the edge of the histogram but most are clustered in the middle. Hence, the skin OCR values are in the more common range found by various investigators.
Topical placement of an optical fiber probe on the skin followed the dynamics of blanching by applied pressure then refilling upon release of the pressure. A spectrogram of reflectance spectra versus time was acquired and analyzed by a simple diffusion theory model to yield the time course of blood volume fraction () and mixed arteriovenous oxygen saturation of hemoglobin (). The and recordings specified the blanched value at and the refilled values at , which yielded values for the blanchable blood volume and the blood flow to this volume, which is presumed to be the superficial plexus just under the epidermis. This plexus exclusively provides oxygen for the epidermis and does not supply significant oxygen to the dermis. So, the oxygen delivery of the plexus is an indicator of the metabolic consumption of oxygen by the epidermis.
Balu et al.1 used two-photon NADH fluorescence to quantify the gradient of OCR across the epidermis. More work is needed to understand the difference between the average OCR measured by the optical fiber measurement and the OCR gradient that they observed. One of the assumptions in the fiber probe analysis is that the rate of cellular oxygen extraction from the blood is equal to the cellular OCR, which assumes that all of the oxygen used by the epidermis comes from the blood. On the other hand, some studies have reported a contribution of atmospheric oxygen to the epidermal oxygen supply.1112.–13
Gaining a better understanding of epidermal oxygen consumption may prove useful in differentiating normal from pathological states. The optical fiber probe is a low-cost portable system that offers a convenient, rapid, noninvasive assessment of the metabolic activity of the epidermis, or any pathologic growth near the skin surface. Such a probe may prove useful in clinical dermatology, which was the motivation for this study.
The authors would like to thank Dr. Sancy Leachman for arranging Dr. Liasi’s participation in the research project at OHSU.
Faezeh Talebi Liasi is a medical doctor and an alumni of the University of Washington School of Medicine. She worked on this project at the Oregon Health and Science University (OHSU) while she was still in medical school. Her research interests include optical characterization of tissue and its use in diagnostics.
Ravikant Samatham received his MS degree in mechanical engineering from the University of Nevada, Reno, in 2004 and a PhD in biomedical engineering from OHSU in 2012. He is now a senior research scientist in the Department of Dermatology at OHSU. His research focuses are second harmonic generation microscopy, optical spectroscopy, and optical imaging.
Steven Jacques received his BS degree in biology from MIT in 1972, his MSEE degree in electrical engineering and computer science, and his PhD in biophysics and medical physics from the University of California-Berkeley in 1984. He has served at Massachusetts General Hospital/Harvard Medical School, the University of Texas M.D. Anderson Cancer Center, the OHSU, and now at Tufts University. His research focuses on light propagation and interactions with tissues.