Optic axis mapping with catheter-based polarization-sensitive optical coherence tomography

Birefringence offers an intrinsic contrast mechanism related to the microstructure and arrangement of fibrillary tissue components. Here we present a reconstruction strategy to recover not only the scalar amount of birefringence but also its optic axis orientation as a function of depth in tissue from measurements with catheter-based polarization sensitive optical coherence tomography. A polarization symmetry constraint, intrinsic to imaging in the backscatter direction, facilitates the required compensation for wavelength-dependent transmission through system elements, the rotating catheter, and overlying tissue layers. Applied to intravascular imaging of coronary atherosclerosis in human patients, the optic axis affords refined interpretation of plaque architecture.

Reflection artifact identification in photoacoustic imaging using multi-wavelength excitation

Photoacoustic imaging has been a focus of research for clinical applications owing to its ability for deep visualization with optical absorption contrast. However, there are various technical challenges remaining for this technique to find its place in clinics. One of the challenges is the occurrence of reflection artifacts. The reflection artifacts may lead to image misinterpretation. Here we propose a new method using multiple wavelengths for identifying and removing the reflection artifacts. By imaging the sample with multiple wavelengths, the spectral response of the features in the photoacoustic image is obtained. We assume that the spectral response of the reflection artifact is better correlated with the proper image feature of its corresponding absorber than with other features in the image. Based on this, the reflection artifacts can be identified and removed. Here, we experimentally demonstrated the potential of this method for real-time identification and correction of reflection artifacts in photoacoustic images in phantoms as well as in vivo using a handheld photoacoustic imaging probe.

Combining high wavenumber and fingerprint Raman spectroscopy for the detection of prostate cancer during radical prostatectomy

For prostate cancer (PCa) patients, radical prostatectomy (complete removal of the prostate) is the only curative surgical option. To date, there is no clinical technique allowing for real-time assessment of surgical margins to minimize the extent of residual cancer. Here, we present a tissue interrogation technique using a dual excitation wavelength Raman spectroscopy system capable of sequentially acquiring fingerprint (FP) and high wavenumber (HWN) Raman spectra. Results demonstrate the ability of the system to detect PCa in post-prostatectomy specimens. In total, 477 Raman spectra were collected from 18 human prostate slices. Each area measured with Raman spectroscopy was characterized as either normal or cancer based on histopathological analyses, and each spectrum was classified based on supervised learning using support vector machines (SVMs). Based on receiver operating characteristic (ROC) analysis, FP (area under the curve [AUC] = 0.89) had slightly superior cancer detection capabilities compared with HWN (AUC = 0.86). Optimal performance resulted from combining the spectral information from FP and HWN (AUC = 0.91), suggesting that the use of these two spectral regions may provide complementary molecular information for PCa detection. The use of leave-one-(spectrum)-out (LOO) or leave-one-patient-out (LOPO) cross-validation produced similar classification results when combining FP with HWN. Our findings suggest that the application of machine learning using multiple data points from the same patient does not result in biases necessarily impacting the reliability of the classification models.

Multi-beam scan analysis with a clinical LINAC for high resolution Cherenkov-excited molecular luminescence imaging in tissue

Cherenkov-excited luminescence scanned imaging (CELSI) is achieved with external beam radiotherapy to map out molecular luminescence intensity or lifetime in tissue. Just as in fluorescence microscopy, the choice of excitation geometry can affect the imaging time, spatial resolution and contrast recovered. In this study, the use of spatially patterned illumination was systematically studied comparing scan shapes, starting with line scan and block patterns and increasing from single beams to multiple parallel beams and then to clinically used treatment plans for radiation therapy. The image recovery was improved by a spatial-temporal modulation-demodulation method, which used the ability to capture simultaneous images of the excitation Cherenkov beam shape to deconvolve the CELSI images. Experimental studies used the multi-leaf collimator on a clinical linear accelerator (LINAC) to create the scanning patterns, and image resolution and contrast recovery were tested at different depths of tissue phantom material. As hypothesized, the smallest illumination squares achieved optimal resolution, but at the cost of lower signal and slower imaging time. Having larger excitation blocks provided superior signal but at the cost of increased radiation dose and lower resolution. Increasing the scan beams to multiple block patterns improved the performance in terms of image fidelity, lower radiation dose and faster acquisition. The spatial resolution was mostly dependent upon pixel area with an optimized side length near 38mm and a beam scan pitch of P = 0.33, and the achievable imaging depth was increased from 14mm to 18mm with sufficient resolving power for 1mm sized test objects. As a proof-of-concept, in-vivo tumor mouse imaging was performed to show 3D rendering and quantification of tissue pO2 with values of 5.6mmHg in a tumor and 77mmHg in normal tissue.

Trans-retinal cellular imaging with multimodal adaptive optics

Adaptive optics (AO), when coupled to different imaging modalities, has enabled resolution of various cell types across the entire retinal depth in the living human eye. Extraction of information from retinal cells is optimal when their optical properties, structure, and physiology are matched to the unique capabilities of each imaging modality. Despite the earlier success of multimodal AO (mAO) approaches, the full capabilities of the individual imaging modalities were often diminished rather than enhanced when integrated into multimodal platforms. Furthermore, many mAO designs added unnecessary complexity, making clinical translation difficult. In this study, we present a novel mAO system that combines two complementary approaches, scanning laser ophthalmoscopy (SLO) and optical coherence tomography (OCT), in one instrument using a simplified optical design, flexible alternation of scanning modes, and independent focus control. The mAO system imaging performance was demonstrated by visualization of cells in their mosaic arrangement across the full depth of the retina in three human subjects, including microglia, nerve fiber bundles, retinal ganglion cells and axons, and capillaries in the inner retina and foveal cones, peripheral rods, and retinal pigment epithelial cells in the outer retina. Multimodal AO is a powerful tool to capture the most complete picture of retinal health.

Erratum: Visualization of laser tattoo removal treatment effects in a mouse model by two-photon microscopy: publisher's note

[This corrects the article on p. 3735 in vol. 8, PMID: 28856046.].

Comparative study of the imaging contrasts of Mueller matrix derived parameters between transmission and backscattering polarimetry

Mueller matrix polarimetry is a potentially powerful tool for biomedical diagnosis. Recently, the transmission Mueller matrix microscope and backscattering Mueller matrix endoscope were developed and applied to various pathological samples. However, a comparative study of imaging contrasts of Mueller matrix derived parameters between transmission and backscattering measurements is still needed to help decide which information obtained from transmission Mueller matrix microscope can be directly applied to in vivo Mueller matrix imaging. Here, to compare the imaging contrasts of Mueller matrix derived parameters between transmission and backscattering polarimetry, we measure porcine liver tissue samples and human breast carcinoma tissue specimens. The experiments and corresponding Monte Carlo stimulation results demonstrate that the backscattering and transmission retardance-related Mueller matrix parameters have very similar contrasts to characterize the anisotropic and isotropic structures of pathological tissues, meaning that the conclusions made from Mueller matrix microscopic imaging based on retardance can also be helpful to guide the in situ backscattering Mueller matrix polarimetric diagnosis. However, the values and contrasts of depolarization-related Mueller matrix parameters have some differences between transmission and backscattering polarimetry.

Enhancing sparse-view photoacoustic tomography with combined virtually parallel projecting and spatially adaptive filtering

To fully realize the potential of photoacoustic tomography (PAT) in preclinical and clinical applications, rapid measurements and robust reconstructions are needed. Sparse-view measurements have been adopted effectively to accelerate the data acquisition. However, since the reconstruction from the sparse-view sampling data is challenging, both the effective measurement and the appropriate reconstruction should be taken into account. In this study, we present an iterative sparse-view PAT reconstruction scheme, where a concept of virtual parallel-projection matching the measurement condition is introduced to aid the "compressive sensing" in the reconstruction procedure, and meanwhile, the non-local spatially adaptive filtering exploring the a priori information of the mutual similarities in natural images is adopted to recover the unknowns in the transformed sparse domain. Consequently, the reconstructed images with the proposed sparse-view scheme can be evidently improved in comparison to those with the universal back-projection method, for the cases of same sparse views. The proposed approach has been validated by the simulations and ex vivo experiments, which exhibits desirable performances in image fidelity even from a small number of measuring positions.

Phase retrieval for arbitrary Fresnel-like linear shift-invariant imaging systems suitable for tomography

We present a generalization of the non-iterative phase retrieval in X-ray phase contrast imaging applicable for an arbitrary linear shift-invariant (LSI) imaging system with a non-negligible amount of free space propagation (termed as Fresnel-like). Our novel approach poses no restrictions on the propagation distance between optical elements of the system. In turn, the requirements are only demanded for the transfer function of the optical elements, which should be approximable by second-order Taylor polynomials. Furthermore, we show that the method can be conveniently used as an initial guess for iterative phase retrieval, resulting in faster convergence. The proposed approach is tested on synthetic and experimentally measured holograms obtained using a Bragg magnifier microscope - a representative of Fresnel-like LSI imaging systems. Finally, the algorithm is applied to a whole micro-tomographic scan of a biological specimen of a tardigrade, revealing morphological details at the spatial resolution of 300 nm - limiting resolution of the actual imaging system.

In vivo detection of UV-induced acute skin effects using optical coherence tomography

Ultraviolet (UV) rays have been identified as a carcinogen with long-term irradiation and are an important risk factor for skin cancer. Here, we report the use of optical coherence tomography/optical coherence tomography angiography (OCT/OCTA) to study acute UV-induced effects on skin in vivo. To understand the relationship between the acute effects and irradiated UV power density, three groups were irradiated with different power densities in our experiments. Furthermore, the same skin area was repeatedly scanned with OCT during UV irradiation to investigate the progress of the induced acute effects and after irradiation for observation of skin recovery. Subsequently, the OCT/OCTA results were quantitatively analyzed to acquire skin thickness and blood-vessel density for comparison. UV-induced acute effects on morphology and microcirculation can be identified from OCT/OCTA results, which showed the increases in the skin thickness and blood-vessel density and even severe damage types such as blisters. The results of quantitative analyses also illustrated that the severity of damage induced by UV irradiation can be distinguished and the skin recovery can be monitored with OCT. Our results indicate that OCT can be a promising tool for early detection of UV-induced acute skin damage.

Label-free assessment of carotid artery biochemical composition using fiber-based fluorescence lifetime imaging

Novel diagnostic tools with the ability to monitor variations in biochemical composition and provide benchmark indicators of vascular tissue maturation are needed to create functional tissue replacements. We investigated the ability of fiber-based, label-free multispectral fluorescent lifetime imaging (FLIm) to quantify the anatomical variations in biochemical composition of native carotid arteries and validated these results against biochemical assays. FLIm-derived parameters in spectral band 415-455 nm correlated with tissue collagen content (R² = 0.64) and cell number (R² = 0.61) and in spectral band 465-553 nm strongly correlated with elastin content (R² = 0.89). These results suggest that FLIm holds great potential for assessing vascular tissue maturation and functional properties based on tissue autofluorescence.

Comparison of speckleplethysmographic (SPG) and photoplethysmographic (PPG) imaging by Monte Carlo simulations and in vivo measurements

Noncontact photoplethysmography (PPG) is limited by a poor signal-to-noise ratio (SNR). A solution to this limitation is the use of alternate sources of optical contrast to generate a complementary pulsatile waveform. One such source is laser speckle contrast, which is modulated in biological tissues by the flow rate of red blood cells. Averaging a region of interest from a speckle contrast image over time allows for the calculation of a speckleplethysmogram (SPG). Similar to PPG, SPG enables monitoring of heart rate and respiratory rate. A gap in the knowledge base exists as to the precise spatiotemporal relationship between PPG and SPG signals. We have developed an eight-layer tissue model to simulate both PPG and SPG signals in a reflectance geometry via Monte Carlo methods. We modeled PPG by compression of the upper and lower blood nets due to expansion of the larger arterial layer below. The in silico PPG peak-to-peak amplitude percent was greater at 532 nm than at 860 nm (5.6% vs. 3.0%, respectively), which matches trends from the literature. We modeled SPG by changing flow speeds of red blood cells in both the capillaries and arterioles over the cardiac cycle. The in silico SPG peak-to-peak amplitude percent was 24% at 532 nm and 40% at 860 nm. In silico results are similar to in vivo results measured with a two-camera set up for simultaneous imaging of PPG and SPG. Both in silico and in vivo data suggest SPG has a much larger SNR than PPG, which may prove beneficial for noncontact, wide-field optical monitoring of cardiovascular health.

Full video pulse extraction

This paper introduces a new method to automate heart-rate detection using remote photoplethysmography (rPPG). The method replaces the commonly used region of interest (RoI) detection and tracking, and does not require initialization. Instead, it combines a number of candidate pulse-signals computed in the parallel, each biased towards differently colored objects in the scene. The method is based on the observation that the temporally averaged colors of video objects (skin and background) are usually quite stable over time in typical application-driven scenarios, such as the monitoring of a subject sleeping in bed, or an infant in an incubator. The resulting system, called full video pulse extraction (FVP), allows the direct use of raw video streams for pulse extraction. Our benchmark set of diverse videos shows that FVP enables long-term sleep monitoring in visible light and in infrared, and works for adults and neonates. Although we only demonstrate the concept for heart-rate monitoring, we foresee the adaptation to a range of vital signs, thus benefiting the larger video health monitoring field.

Wavelet-based method for removing global physiological noise in functional near-infrared spectroscopy

Functional near-infrared spectroscopy (fNIRS) is a fast-developing non-invasive functional brain imaging technology widely used in cognitive neuroscience, clinical research and neural engineering. However, it is a challenge to effectively remove the global physiological noise in the fNIRS signal. The global physiological noise in fNIRS arises from multiple physiological origins in both superficial tissues and the brain. It has complex temporal, spatial and frequency characteristics, casting significant influence on the results. In the present study, we developed a novel wavelet-based method for fNIRS global physiological noise removal. The method is data-driven and does not rely on any additional hardware or subjective noise component selection procedure. It consists of two steps. Firstly, we use wavelet transform coherence to automatically detect the time-frequency points contaminated by the global physiological noise. Secondly, we decompose the fNIRS signal by using the wavelet transform, and then suppress the wavelet energy of the contaminated time-frequency points. Finally, we transform the signal back to a time series. We validated the method by using simulation and real data at both task- and resting-state. The results showed that our method can effectively remove the global physiological noise from the fNIRS signal and improve the spatial specificity of the task activation and the resting-state functional connectivity pattern.

Comparing the effective attenuation lengths for long wavelength in vivo imaging of the mouse brain

Light attenuation in thick biological tissues, caused by a combination of absorption and scattering, limits the penetration depth in multiphoton microscopy (MPM). Both tissue scattering and absorption are dependent on wavelengths, which makes it essential to choose the excitation wavelength with minimum attenuation for deep imaging. Although theoretical models have been established to predict the wavelength dependence of light attenuation in brain tissues, the accuracy of these models in experimental settings needs to be verified. Furthermore, the water absorption contribution to the tissue attenuation, especially at 1450 nm where strong water absorption is predicted to be the dominant contributor in light attenuation, has not been confirmed. Here we performed a systematic study of in vivo three-photon imaging at different excitation wavelengths, 1300 nm, 1450 nm, 1500 nm, 1550 nm, and 1700 nm, and quantified the tissue attenuation by calculating the effective attenuation length at each wavelength. The experimental data show that the effective attenuation length at 1450 nm is significantly shorter than that at 1300 nm or 1700 nm. Our results provide unequivocal validation of the theoretical estimations based on water absorption and tissue scattering in predicting the effective attenuation lengths for long wavelength in vivo imaging.

Real-time cancer detection with an integrated lensless fluorescence contact imager

Microscopic tumor cell foci left in a patient after surgery significantly increase the chance of cancer recurrence. However, fluorescence microscopes used for intraoperative navigation lack the necessary sensitivity for imaging microscopic disease and are too bulky to maneuver within the resection cavity. We have developed a scalable chip-scale fluorescence contact imager for detecting microscopic cancer in vivo and in real-time. The imager has been characterized under simulated in vivo conditions using ex vivo samples, providing strong evidence that our device can be used in vivo. Angle-selective gratings enhance the resolution of the imager without impacting its physical size. We demonstrate detection of cancer cell clusters containing as few as 25 HCC1569 breast cancer cells and 400 LNCaP prostate cancer cells with integration times of only 50 ms and 70 ms, respectively. A cell cluster recognition algorithm is used to achieve both a sensitivity and specificity of 92 % for HCC1569 cell samples, indicating the reliability of the imager. The signal-to-noise ratio (SNR) degradation with increased separation is only 1.5 dB at 250 μm. Blood scattering and absorption reduce the SNR by less than 2 dB for typical concentrations. Moreover, HER2+ breast cancer tissue taken from a patient is distinguished from normal breast tissue with an integration time of only 75 ms.

Novel analysis of fNIRS acquired dynamic hemoglobin concentrations: application in young children with autism spectrum disorder

A novel analysis of the spatial complexity of functional connectivity (SCFC) was proposed to investigate the spatial complexity of multiple dynamic functional connectivity series in an fNIRS study, using an approach combining principal component analysis and normalized entropy. The analysis was designed to describe the complex spatial features of phase synchrony based dynamic functional connectivity (dFC), which are unexplained in traditional approaches. The feasibility and validity of this method were verified in a sample of young patients with autism spectrum disorders (ASD). Our results showed that there were information exchange deficits in the right prefrontal cortex (PFC) of children with ASD, with markedly higher interregion SCFCs between the right PFC and other brain regions than those of normal controls. Furthermore, the global SCFC was significantly higher in young patients with ASD, along with considerably higher intraregion SCFCs in the prefrontal and temporal lobes which represents more diverse information exchange in these areas. The study suggests a novel method to analyze the fNIRS required dynamic hemoglobin concentrations by using concepts of SCFC. Moreover, the clinical results extend our understanding of ASD pathology, suggesting the crucial role of the right PFC during the information exchange process.

Voltage-sensitive dye delivery through the blood brain barrier using adenosine receptor agonist regadenoson

Optical imaging of brain activity has mostly employed genetically manipulated mice, which cannot be translated to clinical human usage. Observation of brain activity directly is challenging due to the difficulty in delivering dyes and other agents through the blood brain barrier (BBB). Using fluorescence imaging, we have demonstrated the feasibility of delivering the near-infrared voltage-sensitive dye (VSD) IR-780 perchlorate to the brain tissue through pharmacological techniques, via an adenosine agonist (regadenoson). Comparison of VSD fluorescence of mouse brains without and with regadenoson showed significantly increased residence time of the fluorescence signal in the latter case, indicative of VSD diffusion into the brain tissue. Dose and timing of regadenoson were varied to optimize BBB permeability for VSD delivery.

Video-rate all-optical ultrasound imaging

All-optical ultrasound imaging, where ultrasound is generated and detected using light, has recently been demonstrated as a viable modality that is inherently insensitive to electromagnetic interference and exhibits wide bandwidths. High-quality 2D and 3D all-optical ultrasound images of tissues have previously been presented; however, to date, long acquisition times (ranging from minutes to hours) have hindered clinical application. Here, we present the first all-optical ultrasound imaging system capable of video-rate, real-time two-dimensional imaging of biological tissue. This was achieved using a spatially extended nano-composite optical ultrasound generator, a highly sensitive fibre-optic acoustic receiver, and eccentric illumination resulting in an acoustic source exhibiting optimal directivity. This source was scanned across a one-dimensional source aperture using a fast galvo mirror, thus enabling the dynamic synthesis of source arrays comprising spatially overlapping sources at non-uniform source separation distances. The resulting system achieved a sustained frame rate of 15 Hz, a dynamic range of 30 dB, a penetration depth of at least 6 mm, a resolution of 75 µm (axial) by 100 µm (lateral), and enabled the dynamics of a pulsating ex vivo carotid artery to be captured.

Stability of gel wax based optical scattering phantoms

Phantoms with tuneable optical scattering properties are essential in the development and refinement of optical based imaging techniques. Mineral oil based 'gel wax' phantoms are the subject of increasing interest due to their ease and speed of manufacture, non-toxic nature, ability to cast into anatomically realistic shapes, as well as their cost-effective nature of production. The addition of scatterers such as titanium dioxide powder and monodisperse silica microspheres to the gel wax allows for the creation of phantoms with a controllable optical scattering coefficient. To enable repeated use of such phantoms, the stability of the scattering properties must be determined-a property which has yet to be investigated. We present an analysis of the stability of the reduced scattering coefficient (μ s ' ) of such phantoms over time. We conclude that due to the measurable reduction in scattering coefficient over time, gel wax phantoms embedded with silica spheres may not be suitable for repeated use over time, however gel wax-TiO2 phantoms are much more temporally stable.

Quantitative diagnosis of tissue microstructure with wide-field high spatial frequency domain imaging

Non-contact and minimally invasive endoscopic optical imaging is an invaluable diagnostic tool for tissue examination and cancer screening. The point sampling techniques with high sensitivity to the tissue microenvironment are time consuming and often not affordable in clinics. There is a major clinical need for a large field-of-view (FOV) rapid screening method to highlight subtle tissue microstructural alterations. To address this unmet need, we have developed High Spatial Frequency Domain Imaging (HSFDI)-a non-contact imaging modality that spatially maps the tissue microscopic scattering structures over a large field of view (>1cm2). Based on an analytical reflectance model of sub-diffusive light from forward-peaked highly scattering media, HSFDI quantifies the spatially-resolved parameters of the light scattering phase function (i.e., the backscattering probability and the light spreading length) from the reflectance of structured light modulated at high spatial frequencies. Enhanced signal to noise ratio (SNR) is achieved at even ultra-high modulation frequencies with single snapshot multiple frequency demodulation (SSMD). The variations in tissue microstructures, including the strength of the background (pudding) refractive index fluctuation and the prominent scattering structure (plum) morphology, can then be inferred. After validation with optical phantoms, measurements of fresh ex vivo tissue samples revealed significant contrast and differentiation of the phase function parameters between different types and disease states (normal, inflammatory, and cancerous) of tissue whereas tissue absorption and reduced scattering coefficients only show modest changes. HSFDI may provide wide-field images of microscopic structural biomarkers unobtainable with either diffuse light imaging or point-based optical sampling. Potential clinical applications include the rapid screening of excised tissue and the noninvasive examination of suspicious lesions during operation.

Real-time localization of the parathyroid gland in surgical field using Raspberry Pi during thyroidectomy: a preliminary report

We created an auto-para viewer, an autofluorescence imaging device, to localize the parathyroid glands during thyroidectomy using an inexpensive Raspberry Pi. A special emission filter in the auto-para viewer was designed to pass 1/100 of visible light and nearly all infrared light longer than 808 nm. With this emission filter, we simultaneously acquired an autofluorescence image of the parathyroid and a visible light image of the surrounding surgical field. The auto-para viewer displayed four times brighter autofluorescence of the parathyroid glands compared to the background tissues without operating room light. Additionally, it showed two times brighter autofluorescence than the background tissues simultaneously showing the surgical field illuminated by the visible light from the operating room light. The NOIR camera, using the auto-para viewer, could reduce the camera's exposure time so the parathyroid glands to be viewed in real-time, which is expected to prevent unintentional damage to the parathyroid gland during thyroidectomy.

Decoding cortical brain states from widefield calcium imaging data using visibility graph

Widefield optical imaging of neuronal populations over large portions of the cerebral cortex in awake behaving animals provides a unique opportunity for investigating the relationship between brain function and behavior. In this paper, we demonstrate that the temporal characteristics of calcium dynamics obtained through widefield imaging can be utilized to infer the corresponding behavior. Cortical activity in transgenic calcium reporter mice (n=6) expressing GCaMP6f in neocortical pyramidal neurons is recorded during active whisking (AW) and no whisking (NW). To extract features related to the temporal characteristics of calcium recordings, a method based on visibility graph (VG) is introduced. An extensive study considering different choices of features and classifiers is conducted to find the best model capable of predicting AW and NW from calcium recordings. Our experimental results show that temporal characteristics of calcium recordings identified by the proposed method carry discriminatory information that are powerful enough for decoding behavior.

Weighting function effects in a direct regularization method for image-guided near-infrared spectral tomography of breast cancer

Structural image-guided near-infrared spectral tomography (NIRST) has been developed as a way to use diffuse NIR spectroscopy within the context of image-guided quantification of tissue spectral features. A direct regularization imaging (DRI) method for NIRST has the value of not requiring any image segmentation. Here, we present a comprehensive investigational study to analyze the impact of the weighting function implied when weighting the recovery of optical coefficients in DRI based NIRST. This was done using simulations, phantom and clinical patient exam data. Simulations where the true object is known indicate that changes to this weighting function can vary the contrast by 10%, the contrast to noise ratio by 20% and the full width half maximum (FWHM) by 30%. The results from phantoms and human images show that a linear inverse distance weighting function appears optimal, and that incorporation of this function can generally improve the recovered total hemoglobin contrast of the tumor to the normal surrounding tissue by more than 15% in human cases.

Prediction of the anti-glioma therapeutic effects of temozolomide through in vivo molecular imaging of MMP expression

Currently, there is no effective way to assess the therapeutic response of temozolomide (TMZ) for the glioma. In this study, the human U87MG-fLuc glioma animal models were set up and the antitumor efficacy of TMZ was evaluated using bioluminescence imaging (BLI) and MRI. Then, bioluminescence tomography (BLT) was reconstructed using an adaptive sparsity matching pursuit (ASMP) algorithm. Second, the expression level of the MMP-750 probe was examined with or without TMZ treatment using FMI. Third, the expression of MMP2 and MMP3 was specifically examined after treatment. The results showed that TMZ effectively inhibited glioma growth. The targeted imaging of MMP-750 was decreased during the treatment of glioma with TMZ. Moreover, the MMP2 and MMP3 expression was found to correlate with the inhibition effect of TMZ. Our study indicated that the therapeutic effects of TMZ can be effectively evaluated at an early stage using molecular imaging, and MMP targeting the fluorescence probe could be utilized for the prediction and assessment of the therapeutic effects of TMZ.

Erratum: Visualization of laser tattoo removal treatment effects in a mouse model by two-photon microscopy: erratum

[This corrects the article on p. 3735 in vol. 8, PMID: 28856046.].

Volumetric non-local-means based speckle reduction for optical coherence tomography

We present a novel tomographic non-local-means based despeckling technique, TNode, for optical coherence tomography. TNode is built upon a weighting similarity criterion derived for speckle in a three-dimensional similarity window. We present an implementation using a two-dimensional search window, enabling the despeckling of volumes in the presence of motion artifacts, and an implementation using a three-dimensional window with improved performance in motion-free volumes. We show that our technique provides effective speckle reduction, comparable with B-scan compounding or out-of-plane averaging, while preserving isotropic resolution, even to the level of speckle-sized structures. We demonstrate its superior despeckling performance in a phantom data set, and in an ophthalmic data set we show that small, speckle-sized retinal vessels are clearly preserved in intensity images en-face and in two orthogonal, cross-sectional views. TNode does not rely on dictionaries or segmentation and therefore can readily be applied to arbitrary optical coherence tomography volumes. We show that despeckled esophageal volumes exhibit improved image quality and detail, even in the presence of significant motion artifacts.

Pilot study of freshly excised breast tissue response in the 300-600 GHz range

The failure to accurately define tumor margins during breast conserving surgery (BCS) results in a 20% re-excision rate. The present paper reports the investigation to evaluate the potential of terahertz imaging for breast tissue recognition within the under-explored 300 - 600 GHz range. Such a frequency window matches new BiCMOS technology capabilities and thus opens up the opportunity for near-field terahertz imaging using these devices. To assess the efficacy of this frequency band, data from 16 freshly excised breast tissue samples were collected and analyzed directly after excision. Complex refractive indices have been extracted over the as-mentioned frequency band, and amplitude frequency images show some contrast between tissue types. Principal component analysis (PCA) has also been applied to the data in an attempt to automate tissue classification. Our observations suggest that the dielectric response could potentially provide contrast for breast tissue recognition within the 300 - 600 GHz range. These results open the way for silicon-based terahertz subwavelength near field imager design, efficient up to 600 GHz to address ex vivo life-science applications.

Fundus autofluorescence beyond lipofuscin: lesson learned from ex vivo fluorescence lifetime imaging in porcine eyes

Fundus autofluorescence (FAF) imaging is a well-established method in ophthalmology; however, the fluorophores involved need more clarification. The FAF lifetimes of 20 post mortem porcine eyes were measured in two spectral channels using fluorescence lifetime imaging ophthalmoscopy (FLIO) and compared with clinical data from 44 healthy young subjects. The FAF intensity ratio of the short and the long wavelength emission (spectral ratio) was determined. Ex vivo porcine fundus fluorescence emission is generally less intense than that seen in human eyes. The porcine retina showed significantly (p<0.05) longer lifetimes than the retinal pigment epithelium (RPE): 584 ± 128 ps vs. 121 ± 55 ps 498-560 nm, 240 ± 42 ps vs. 125 ± 20 ps at 560-720 nm. Furthermore, the lifetimes of the porcine RPE were significantly shorter (121 ± 55 ps and 125 ± 20 ps) than those measured from human fundus in vivo (162 ± 14 ps and 179 ± 13 ps, respectively). The fluorescence emission of porcine retina was shifted towards a shorter wavelength compared to that of RPE and human FAF. This data shows the considerable contribution of fluorophores in the neural retina to total FAF intensity in porcine eyes.

Cellular in vivo 3D imaging of the cornea by confocal laser scanning microscopy

We present an in vivo confocal laser scanning microscopy based method for large 3D reconstruction of the cornea on a cellular level with cropped volume sizes up to 266 x 286 x 396 µm3. The microscope objective used is equipped with a piezo actuator for automated, fast and precise closed-loop focal plane control. Furthermore, we present a novel concave surface contact cap, which significantly reduces eye movements by up to 87%, hence increasing the overlapping image area of the whole stack. This increases the cuboid volume of the generated 3D reconstruction significantly. The possibility to generate oblique sections using isotropic volume stacks opens the window to slit lamp microscopy on a cellular level.

Improving mesoscopic fluorescence molecular tomography via preconditioning and regularization

Mesoscopic fluorescence molecular tomography (MFMT) is a novel imaging technique capable of obtaining 3-D distribution of molecular probes inside biological tissues at depths of a few millimeters with a resolution up to ~100 μm. However, the ill-conditioned nature of the MFMT inverse problem severely deteriorates its reconstruction performances. Furthermore, dense spatial sampling and fine discretization of the imaging volume required for high resolution reconstructions make the sensitivity matrix (Jacobian) highly correlated, which prevents even advanced algorithms from achieving optimal solutions. In this work, we propose two computational methods to respectively increase the incoherence of the sensitivity matrix and improve the convergence rate of the inverse solver. We first apply a compressed sensing (CS) based preconditioner on either the whole sensitivity matrix or sub sensitivity matrices to reduce the coherence between columns of the sensitivity matrix. Then we employed a regularization method based on the weight iterative improvement method (WIIM) to mitigate the ill-condition of the sensitivity matrix and to drive the iterative optimization process towards convergence at a faster rate. We performed numerical simulations and phantom experiments to validate the effectiveness of the proposed strategies. In both in silico and in vitro cases, we were able to improve the quality of MFMT reconstructions significantly.

A spread spectrum approach to time-domain near-infrared diffuse optical imaging using inexpensive optical transceiver modules

We introduce a compact time-domain system for near-infrared spectroscopy using a spread spectrum technique. The proof-of-concept single channel instrument utilises a low-cost commercially available optical transceiver module as a light source, controlled by a Kintex 7 field programmable gate array (FPGA). The FPGA modulates the optical transceiver with maximum-length sequences at line rates up to 10Gb/s, allowing us to achieve an instrument response function with full width at half maximum under 600ps. The instrument is characterised through a set of detailed phantom measurements as well as proof-of-concept in vivo measurements, demonstrating performance comparable with conventional pulsed time-domain near-infrared spectroscopy systems.

Three-dimensional tomographic microscopy technique with multi-frequency combination with partially coherent illuminations

We demonstrate a three-dimensional (3D) optical diffraction tomographic technique with multi-frequency combination (MFC-ODT) for the 3D quantitative phase imaging of unlabeled specimens. Three sets of through-focus intensity images are captured under an annular aperture and two circular apertures with different coherence parameters. The 3D phase optical transfer functions (POTF) corresponding to different illumination apertures are combined to obtain a synthesized frequency response, achieving high-quality, low-noise 3D reconstructions with imaging resolution up to the incoherent diffraction limit. Besides, the expression of 3D POTF for arbitrary illumination pupils is derived and analyzed, and the 3D imaging performance of annular illumination is explored. It is shown that the phase-contrast washout effect in high-NA circular apertures can be effectively addressed by introducing a complementary annular aperture, which strongly boosts the phase contrast and improves the imaging resolution. By incorporating high-NA illumination as well as high-NA detection, MFC-ODT can achieve a theoretical transverse resolution up to 200 nm and an axial resolution of 645 nm. To test the feasibility of the proposed MFC-ODT technique, the 3D refractive index reconstruction results are based on a simulated 3D resolution target and experimental investigations of micro polystyrene bead and unstained biological samples are presented. Due to its capability for high-resolution 3D phase imaging as well as the compatibility with a widely available commercial microscope, the MFC-ODT is expected to find versatile applications in biological and biomedical research.

Distinction of directional coupling in sensorimotor networks between active and passive finger movements using fNIRS

The purpose of this study is to investigate cerebral cortex activation during active movement and passive movement by using a functional near-infrared spectroscopy (fNIRS). Tasks were the flexion/extension of the right hand finger by active movement and passive movement. Oxy-hemoglobin concentration changes calculated from fNIRS and analyzed the activation and connectivity so as to understand dynamical brain relationship. The results demonstrated that the brain activation in passive movements is similar to motor execution. During active movement, the estimated causality patterns showed significant causality value from the supplementary motor area (SMA) to the primary motor cortex (M1). During the passive movement, the causality from the primary somatosensory cortex (S1) to the primary motor cortex (M1) was stronger than active movement. These results demonstrated that active and passive movements had a direct effect on the cerebral cortex but the stimulus pathway of active and passive movement is different. This study may contribute to better understanding how active and passive movements can be expressed into cortical activation by means of fNIRS.

Erratum: Methods for detecting host genetic modifiers of tumor vascular function using dynamic near-infrared fluorescence imaging: errata

[This corrects the article on p. 543 in vol. 9, PMID: 29552392.].

Simultaneous monitoring of cerebral perfusion and cytochrome c oxidase by combining broadband near-infrared spectroscopy and diffuse correlation spectroscopy

Preterm infants born with very low birth weights are at a high risk of brain injury, in part because the premature brain is believed to be prone to periods of low cerebral blood flow (CBF). Tissue damage is likely to occur if reduction in CBF is sufficient to impair cerebral energy metabolism for extended periods. Therefore, a neuromonitoring method that can detect reductions in CBF, large enough to affect metabolism, could alert the neonatal intensive care team before injury occurs. In this report, we present the development of an optical system that combines diffuse correlation spectroscopy (DCS) for monitoring CBF and broadband near-infrared spectroscopy (B-NIRS) for monitoring the oxidation state of cytochrome c oxidase (oxCCO) - a key biomarker of oxidative metabolism. The hybrid instrument includes a multiplexing system to enable concomitant DCS and B-NIRS measurements while avoiding crosstalk between the two subsystems. The ability of the instrument to monitor dynamic changes in CBF and oxCCO was demonstrated in a piglet model of neonatal hypoxia-ischemia (HI). Experiments conducted in eight animals, including two controls, showed that oxCCO exhibited a delayed response to ischemia while CBF and tissue oxygenation (StO2) responses were instantaneous. These findings suggest that simultaneous neuromonitoring of perfusion and metabolism could provide critical information regarding clinically significant hemodynamic events prior to the onset of brain injury.

Spectral-resolved cone-beam X-ray luminescence computed tomography with principle component analysis

Cone-beam X-ray luminescence computed tomography (CB-XLCT) has become a promising technique for its higher utilization of X-ray and shorter scanning time compared to the narrow-beam XLCT, but it suffers from the low-spatial resolution that results in the insufficiency to resolve the adjacent multiple probes. In multispectral CB-XLCT, multiple probes show different emission behaviors in the dimension of the spectrum. In this work, a spectral-resolved CB-XLCT method combining multispectral CB-XLCT with principle component analysis (PCA) was proposed to improve the imaging resolution. Results of digital simulation and the phantom experiment illustrated that the proposed method was capable of resolving adjacent multiple probes accurately and had better performance than the common multispectral CB-XLCT with spectrum information priori.

Optimization of wavelength selection for multispectral image acquisition: a case study of atrial ablation lesions

In vivo autofluorescence hyperspectral imaging of moving objects can be challenging due to motion artifacts and to the limited amount of acquired photons. To address both limitations, we selectively reduced the number of spectral bands while maintaining accurate target identification. Several downsampling approaches were applied to data obtained from the atrial tissue of adult pigs with sites of radiofrequency ablation lesions. Standard image qualifiers such as the mean square error, the peak signal-to-noise ratio, the structural similarity index map, and an accuracy index of lesion component images were used to quantify the effects of spectral binning, an increased spectral distance between individual bands, as well as random combinations of spectral bands. Results point to several quantitative strategies for deriving combinations of a small number of spectral bands that can successfully detect target tissue. Insights from our studies can be applied to a wide range of applications.

Terahertz time-domain spectroscopy for non-invasive assessment of water content in biological samples

We apply terahertz time-domain spectroscopy for the quantitative non-invasive assessment of the water content in biological samples, such as Carpinus caroliniana tree leaves and pork muscles. The developed experimental terahertz time-domain spectroscopy system operates both in transmission and reflection modes. The Landau-Looyenga-Lifshitz-based model is used for the calculation of the water concentration within the samples. The results of the water concentration measurements are compared with the results of the gravimetric measurements. The obtained results show that the water content in biological samples can be measured non-invasively, with a high accuracy, utilizing terahertz waves in transmission and reflection modes.

Tomographic fluorescence lifetime multiplexing in the spatial frequency domain

The ability to simultaneously recover multiple fluorophores within biological tissue (multiplexing) can have important applications for tracking parallel disease processes in vivo. Here we present a novel method for rapid and quantitative multiplexing within a scattering medium, such as biological tissue, based on fluorescence lifetime contrast. This method employs a tomographic inversion of the asymptotic (late) portion of time-resolved spatial frequency (SF) domain measurements. Using Monte Carlo simulations and phantom experiments, we show that the SF-asymptotic time domain (SF-ATD) approach provides a several-fold improvement in relative quantitation and localization accuracy over conventional SF-time domain inversion. We also show that the SF-ATD approach can exploit selective filtering of high spatial frequencies to dramatically improve reconstruction accuracy for fluorophores with subnanosecond lifetimes, which is typical of most near-infrared fluorophores. These results suggest that the SF-ATD approach will serve as a powerful new tool for whole-body lifetime multiplexing.

Fast volume-scanning light sheet microscopy reveals transient neuronal events

Light sheet fluorescence microscopy offers considerable potential to the cellular neuroscience community as it makes it possible to image extensive areas of neuronal structures, such as axons or dendrites, with a low light budget, thereby minimizing phototoxicity. However, the shallow depth of a light sheet, which is critical for achieving high contrast, well resolved images, adds a significant challenge if fast functional imaging is also required, as multiple images need to be collected across several image planes. Consequently, fast functional imaging of neurons is typically restricted to a small tissue volume where part of the neuronal structure lies within the plane of a single image. Here we describe a method by which fast functional imaging can be achieved across a much larger tissue volume; a custom-built light sheet microscope is presented that includes a synchronized galvo mirror and electrically tunable lens, enabling high speed acquisition of images across a configurable depth. We assess the utility of this technique by acquiring fast functional Ca²⁺ imaging data across a neuron's dendritic arbour in mammalian brain tissue.

Nth-order linear algorithm for diffuse correlation tomography

The current approaches to imaging the tissue blood flow index (BFI) from diffuse correlation tomography (DCT) data are either an analytical solution or a finite element method, both of which are unable to simultaneously account for the tissue heterogeneity and fully utilize the DCT data. In this study, a new imaging concept for DCT, namely NL-DCT, was created by us in which the medical images are combined with light Monte Carlo simulation to provide geometrical and heterogeneous information in tissue. Moreover, the DCT data at multiple delay time are fully utilized via iterative linear regression. The unique merit of NL-DCT in utilizing the medical images as prior information, when combined with a split Bregman algorithm for total variation minimization (Bregman-TV), was validated on a realistic human head model. Computer simulation outcomes demonstrate the accuracy and robustness of NL-DCT in localizing and separating the flow anomalies as well as the capability to preserve edges of anomalies.

Compact fiber-based multi-photon endoscope working at 1700 nm

We present the design, implementation and performance analysis of a compact multi-photon endoscope based on a piezo electric scanning tube. A miniature objective lens with a long working distance and a high numerical aperture (≈ 0.5) is designed to provide a diffraction limited spot size. Furthermore, a 1700 nm wavelength femtosecond fiber laser is used as an excitation source to overcome the scattering of biological tissues and reduce water absorption. Therefore, the novel optical system along with the unique wavelength allows us to increase the imaging depth. We demonstrate that the endoscope is capable of performing third and second harmonic generation (THG/SHG) and three-photon excitation fluorescence (3PEF) imaging over a large field of view (> 400 μm) with high lateral resolution (2.2 μm). The compact and lightweight probe design makes it suitable for minimally-invasive in-vivo imaging as a potential alternative to surgical biopsies.

3D + time blood flow mapping using SPIM-microPIV in the developing zebrafish heart

We present SPIM-μPIV as a flow imaging system, capable of measuring in vivo flow information with 3D micron-scale resolution. Our system was validated using a phantom experiment consisting of a flow of beads in a 50 μm diameter FEP tube. Then, with the help of optical gating techniques, we obtained 3D + time flow fields throughout the full heartbeat in a ∼3 day old zebrafish larva using fluorescent red blood cells as tracer particles. From this we were able to recover 3D flow fields at 31 separate phases in the heartbeat. From our measurements of this specimen, we found the net pumped blood volume through the atrium to be 0.239 nL per beat. SPIM-μPIV enables high quality in vivo measurements of flow fields that will be valuable for studies of heart function and fluid-structure interaction in a range of small-animal models.

Intravital imaging of a pulmonary endothelial surface layer in a murine sepsis model

Direct intravital imaging of an endothelial surface layer (ESL) in pulmonary microcirculation could be a valuable approach to investigate the role of a vascular endothelial barrier in various pathological conditions. Despite its importance as a marker of endothelial cell damage and impairment of the vascular system, in vivo visualization of ESL has remained a challenging technical issue. In this work, we implemented a pulmonary microcirculation imaging system integrated to a custom-design video-rate laser scanning confocal microscopy platform. Using the system, a real-time cellular-level microscopic imaging of the lung was successfully performed, which facilitated a clear identification of individual flowing erythrocytes in pulmonary capillaries. Subcellular level pulmonary ESL was identified in vivo by fluorescence angiography using a dextran conjugated fluorophore to label blood plasma and the red blood cell (RBC) exclusion imaging analysis. Degradation of ESL width was directly evaluated in a murine sepsis model in vivo, suggesting an impairment of pulmonary vascular endothelium and endothelial barrier dysfunction.

Assessing blood vessel perfusion and vital signs through retinal imaging photoplethysmography

One solution to the global challenge of increasing ocular disease is a cost-effective technique for rapid screening and assessment. Current ophthalmic imaging techniques, e.g. scanning and ocular blood flow systems, are expensive, complex to operate and utilize invasive contrast agents during assessment. The work presented here demonstrates a simple retinal imaging photoplethysmography (iPPG) system with the potential to provide screening, diagnosis, monitoring and assessment that is non-invasive, painless and radiationless. Time series of individual retinal blood vessel images, captured with an eye fundus camera, are processed using standard filtering, amplitude demodulation and principle component analysis (PCA) methods to determine the values of the heart rate (HR) and respiration rate (RR), which are in compliance with simultaneously obtained measurements using commercial pulse oximetry. It also seems possible that some information on the dynamic changes in oxygen saturation levels (SpO₂) in a retinal blood vessel may also be obtained. As a consequence, the retinal iPPG modality system demonstrates a potential avenue for rapid ophthalmic screening, and even early diagnosis, against ocular disease without the need for fluorescent or contrast agents.

Dichroism-sensitive photoacoustic computed tomography

Photoacoustic computed tomography (PACT), a fast-developing modality for deep tissue imaging, images the spatial distribution of optical absorption. PACT usually treats the absorption coefficient as a scalar. However, the absorption coefficients of many biological tissues exhibit an anisotropic property, known as dichroism or diattenuation, which depends on molecular conformation and structural alignment. Here we present a novel imaging method called dichroism-sensitive PACT (DS-PACT), which measures both the amplitude of tissue's dichroism and the orientation of the optic axis of uniaxial dichroic tissue. By modulating the polarization of linearly polarized light and measuring the alternating signals through lock-in detection, DS-PACT can boost dichroic signals from biological tissues. To validate the proposed approach, we experimentally demonstrated the performance of DS-PACT by imaging plastic polarizers and ex vivo bovine tendons deep inside scattering media. We successfully detected the orientation of the optic axis of uniaxial dichroic materials, even at a depth of 4.5 transport mean free paths. We anticipate that the proposed method will extend the capability of PACT to imaging tissue absorption anisotropy.

Burn-injured tissue detection for debridement surgery through the combination of non-invasive optical imaging techniques

The process of burn debridement is a challenging technique requiring significant skills to identify the regions that need excision and their appropriate excision depths. In order to assist surgeons, a machine learning tool is being developed to provide a quantitative assessment of burn-injured tissue. This paper presents three non-invasive optical imaging techniques capable of distinguishing four kinds of tissue-healthy skin, viable wound bed, shallow burn, and deep burn-during serial burn debridement in a porcine model. All combinations of these three techniques have been studied through a k-fold cross-validation method. In terms of global performance, the combination of all three techniques significantly improves the classification accuracy with respect to just one technique, from 0.42 up to more than 0.76. Furthermore, a non-linear spatial filtering based on the mode of a small neighborhood has been applied as a post-processing technique, in order to improve the performance of the classification. Using this technique, the global accuracy reaches a value close to 0.78 and, for some particular tissues and combination of techniques, the accuracy improves by 13%.

Incorporating reflection boundary conditions in the Neumann series radiative transport equation: application to photon propagation and reconstruction in diffuse optical imaging

We propose a formalism to incorporate boundary conditions in a Neumann-series-based radiative transport equation. The formalism accurately models the reflection of photons at the tissue-external medium interface using Fresnel's equations. The formalism was used to develop a gradient descent-based image reconstruction technique. The proposed methods were implemented for 3D diffuse optical imaging. In computational studies, it was observed that the average root-mean-square error (RMSE) for the output images and the estimated absorption coefficients reduced by 38% and 84%, respectively, when the reflection boundary conditions were incorporated. These results demonstrate the importance of incorporating boundary conditions that model the reflection of photons at the tissue-external medium interface.

Measuring water contents in animal organ tissues using terahertz spectroscopic imaging

We investigated the water contents in several organ tissues such as the liver, spleen, kidney, and brain tissue of rats using the terahertz spectroscopic imaging technique. The water contents of the tissues were determined by using a simple equation containing the absorption coefficients of fresh and lyophilized tissues and water. We compared the measured water contents with the difference in mass of tissues before and after lyophilization. All results showed a good match except for the kidney, which has several Bowman's capsules.