Detecting cell death with optical coherence tomography and envelope statistics

A novel algorithm for multiplicative speckle noise reduction in ex vivo human brain OCT images

Optical coherence tomography (OCT) images of ex vivo human brain tissue are corrupted by multiplicative speckle noise that degrades the contrast to noise ratio (CNR) of microstructural compartments. This work proposes a novel algorithm to reduce noise corruption in OCT images that minimizes the penalized negative log likelihood of gamma distributed speckle noise. The proposed method is formulated as a majorize-minimize problem that reduces to solving an iterative regularized least squares optimization. We demonstrate the usefulness of the proposed method by removing speckle in simulated data, phantom data and real OCT images of human brain tissue. We compare the proposed method with state of the art filtering and non-local means based denoising methods. We demonstrate that our approach removes speckle accurately, improves CNR between different tissue types and better preserves small features and edges in human brain tissue.

Surgical scene generation and adversarial networks for physics-based iOCT synthesis

The development and integration of intraoperative optical coherence tomography (iOCT) into modern operating rooms has motivated novel procedures directed at improving the outcome of ophthalmic surgeries. Although computer-assisted algorithms could further advance such interventions, the limited availability and accessibility of iOCT systems constrains the generation of dedicated data sets. This paper introduces a novel framework combining a virtual setup and deep learning algorithms to generate synthetic iOCT data in a simulated environment. The virtual setup reproduces the geometry of retinal layers extracted from real data and allows the integration of virtual microsurgical instrument models. Our scene rendering approach extracts information from the environment and considers iOCT typical imaging artifacts to generate cross-sectional label maps, which in turn are used to synthesize iOCT B-scans via a generative adversarial network. In our experiments we investigate the similarity between real and synthetic images, show the relevance of using the generated data for image-guided interventions and demonstrate the potential of 3D iOCT data synthesis.

Deuterium MRS of early treatment-induced changes in tumour lactate in vitro

Elevated production of lactate is a key characteristic of aberrant tumour cell metabolism and can be non-invasively measured as an early marker of tumour response using deuterium (² H) MRS. Following treatment, changes in the ² H-labelled lactate signal could identify tumour cell death or impaired metabolic function, which precede morphological changes conventionally used to assess tumour response. In this work, the association between apoptotic cell death, extracellular lactate concentration, and early treatment-induced changes in the ² H-labelled lactate signal was established in an in vitro tumour model. Experiments were conducted at 7 T on acute myeloid leukaemia (AML) cells, which had been treated with 10 μg/mL of the chemotherapeutic agent cisplatin. At 24 and 48 h after cisplatin treatment the cells were supplied with 20 mM of [6,6'-² H₂ ]glucose and scanned over 2 h using a two-dimensional ² H MR spectroscopic imaging sequence. The resulting signals from ² H-labelled glucose, lactate, and water were quantified using a spectral fitting algorithm implemented on the Oxford Spectroscopy Analysis MATLAB toolbox. After scanning, the cells were processed for histological stains (terminal deoxynucleotidyl transferase UTP nick end labelling and haematoxylin and eosin) to assess apoptotic area fraction and cell morphology respectively, while a colorimetric assay was used to measure extracellular lactate concentrations in the supernatant. Significantly lower levels of ² H-labelled lactate were observed in the 48 h treated cells compared with the untreated and 24 h treated cells, and these changes were significantly correlated with an increase in apoptotic fraction and a decrease in extracellular lactate. By establishing the biological processes associated with treatment-induced changes in the ² H-labelled lactate signal, these findings suggest that ² H MRS of lactate may be valuable in evaluating early tumour response.

Volumetric tumor delineation and assessment of its early response to radiotherapy with optical coherence tomography

Texture analyses of optical coherence tomography (OCT) images have shown initial promise for differentiation of normal and tumor tissues. This work develops a fully automatic volumetric tumor delineation technique employing quantitative OCT image speckle analysis based on Gamma distribution fits. We test its performance in-vivo using immunodeficient mice with dorsal skin window chambers and subcutaneously grown tumor models. Tumor boundaries detection is confirmed using epi-fluorescence microscopy, combined photoacoustic-ultrasound imaging, and histology. Pilot animal study of tumor response to radiotherapy demonstrates high accuracy, objective nature, novelty of the proposed method in the volumetric separation of tumor and normal tissues, and the sensitivity of the fitting parameters to radiation-induced tissue changes. Overall, the developed methodology enables hitherto impossible longitudinal studies for detecting subtle tissue alterations stemming from therapeutic insult.

A Miniaturized Platform for Multiplexed Drug Response Imaging in Live Tumors

Simple Summary

We have developed an implantable microdevice that is placed into a live tumor, and can directly image how effective various chemotherapy drugs are at inducing cell death, without having to remove or process the tumor tissue. Currently drug optimization is performed by assessing tumor shrinkage after treating a patient with systemic doses of a chemotherapy agent; this only evaluates a single treatment at a time and typically takes weeks-months before an optimal treatment strategy is found (if found at all) for a specific patient. In contrast, using the technology presented here, a personalized cancer treatment strategy can potentially be optimized and tailored to a specific patient’s tumor characteristics within several hours, without requiring surgical tissue removal or prolonged trials of potentially ineffective chemotherapies.

Abstract

By observing the activity of anti-cancer agents directly in tumors, there is potential to greatly expand our understanding of drug response and develop more personalized cancer treatments. Implantable microdevices (IMD) have been recently developed to deliver microdoses of chemotherapeutic agents locally into confined regions of live tumors; the tissue can be subsequently removed and analyzed to evaluate drug response. This method has the potential to rapidly screen multiple drugs, but requires surgical tissue removal and only evaluates drug response at a single timepoint when the tissue is excised. Here, we describe a “lab-in-a-tumor” implantable microdevice (LIT-IMD) platform to image cell-death drug response within a live tumor, without requiring surgical resection or tissue processing. The LIT-IMD is inserted into a live tumor and delivers multiple drug microdoses into spatially discrete locations. In parallel, it locally delivers microdose levels of a fluorescent cell-death assay, which diffuses into drug-exposed tissues and accumulates at sites of cell death. An integrated miniaturized fluorescence imaging probe images each region to evaluate drug-induced cell death. We demonstrate ability to evaluate multi-drug response over 8 h using murine tumor models and show correlation with gold-standard conventional fluorescence microscopy and histopathology. This is the first demonstration of a fully integrated platform for evaluating multiple chemotherapy responses in situ. This approach could enable a more complete understanding of drug activity in live tumors, and could expand the utility of drug-response measurements to a wide range of settings where surgery is not feasible.

Monitoring Quantitative Ultrasound Parameter Changes in a Cell Pellet Model of Cell Starvation

Although it has previously been shown that the spectral analysis of ultrasound backscatter data is sensitive to the cellular changes caused by apoptosis, the sensitivity of spectral analysis to oncosis or ischemic cell death had not previously been studied. Whereas many anticancer treatments induce apoptosis, others induce cell starvation, or oncosis. HT-29 colon adenocarcinoma cells were formed into pellets and covered in phosphate-buffered saline at room temperature for 56 h. The pellets were imaged every 8 h with high-frequency (55 MHz) ultrasound and the raw radio-frequency data processed. The lack of nutrients available to the cells induced cell death due to oncosis. The attenuation slope, speed of sound, spectral slope, and midband fit were estimated at each of the eight time points to identify changes as the cells died due to starvation. The spectral slope decreased monotonically over the 56 h, whereas the attenuation slope showed an increase between 1 and 48 h, followed by a slight decrease between 48 and 56 h. The midband fit did not vary over time. The speed of sound increased from 1514 to 1532 m/s over the first 24 h, after which time it plateaued. These in vitro results indicate different trends in ultrasound parameter changes from those of in vitro apoptotic cells, suggesting that these different methods of cell death can be identified not only by morphological markers, but also by specific ultrasound signatures.

Influence of eye biometrics and corneal micro-structure on noncontact tonometry

Purpose

Tonometry is widely used as the main screening tool supporting glaucoma diagnosis. Still, its accuracy could be improved if full knowledge about the variation of the corneal biomechanical properties was available. In this study, Optical Coherence Tomography (OCT) speckle statistics are used to infer the organisation of the corneal micro-structure and hence, to analyse its influence on intraocular pressure (IOP) measurements.

Methods

Fifty-six subjects were recruited for this prospective study. Macro and micro-structural corneal parameters as well as subject age were considered. Macro-structural analysis included the parameters that are associated with the ocular anatomy, such as central corneal thickness (CCT), corneal radius, axial length, anterior chamber depth and white-to-white corneal diameter. Micro-structural parameters which included OCT speckle statistics were related to the internal organisation of the corneal tissue and its physiological changes during lifetime. The corneal speckle obtained from OCT was modelled with the Generalised Gamma (GG) distribution that is characterised with a scale parameter and two shape parameters.

Results

In macro-structure analysis, only CCT showed a statistically significant correlation with IOP (R² = 0.25, p<0.001). The scale parameter and the ratio of the shape parameters of GG distribution showed statistically significant correlation with IOP (R² = 0.19, p<0.001 and R² = 0.17, p<0.001, respectively). For the studied group, a weak, although significant correlation was found between age and IOP (R² = 0.053, p = 0.04). Forward stepwise regression showed that CCT and the scale parameter of the Generalised Gamma distribution can be combined in a regression model (R² = 0.39, p<0.001) to study the role of the corneal structure on IOP.

Conclusions

We show, for the first time, that corneal micro-structure influences the IOP measurements obtained from noncontact tonometry. OCT speckle statistics can be employed to learn about the corneal micro-structure and hence, to further calibrate the IOP measurements.

Assessment of corneal properties based on statistical modeling of OCT speckle

A new approach to assess the properties of the corneal micro-structure in vivo based on the statistical modeling of speckle obtained from Optical Coherence Tomography (OCT) is presented. A number of statistical models were proposed to fit the corneal speckle data obtained from OCT raw image. Short-term changes in corneal properties were studied by inducing corneal swelling whereas age-related changes were observed analyzing data of sixty-five subjects aged between twenty-four and seventy-three years. Generalized Gamma distribution has shown to be the best model, in terms of the Akaike's Information Criterion, to fit the OCT corneal speckle. Its parameters have shown statistically significant differences (Kruskal-Wallis, p < 0.001) for short and age-related corneal changes. In addition, it was observed that age-related changes influence the corneal biomechanical behaviour when corneal swelling is induced. This study shows that Generalized Gamma distribution can be utilized to modeling corneal speckle in OCT in vivo providing complementary quantified information where micro-structure of corneal tissue is of essence.

Assesment of apoptosis induced changes in scattering using optical coherence tomography

The aim of this study is to identify changes in scattering with optical coherence tomography (OCT) and relate these measurements with mitochondrial changes during the initiation of apoptosis. Human retinal pigment epithelial cells were cultured and apoptosis was induced using 10% alcohol. Using the attenuation coefficient and backscattering, changes were measured during cell death in a cell-pellet and monolayer respectively. To confirm apoptosis, fluorescent activated cell sorting was used. Mitochondrial activity during apoptosis was assessed using an oxidative stress assay and fluorescent confocal microscopy. Pelleted apoptotic cells measured with OCT showed a clear rise while untreated cells showed a very small increase in attenuation coefficient. Monolayered apoptotic cells displayed a distinct increase, while untreated cells showed a small increase in the backscattering. Apoptosis was confirmed by FACS experiments. Mitochondrial changes during the onset of apoptosis were also measured. The results demonstrate that apoptotic cell death could be monitored in real-time by OCT. Changes in the scattering after induction of apoptosis are likely to be related to changes in the intracellular morphology. Oxidative stress-induced mitochondrial swelling could be responsible for the initial increase, while cell blebbing and secondary necrosis subsequently for the observed decrease in scattering.

Probability density function formalism for optical coherence tomography signal analysis: a controlled phantom study

The distribution of backscattered intensities as described by the probability density function (PDF) of tissue-scattered light contains information that may be useful for tissue assessment and diagnosis, including characterization of its pathology. In this Letter, we examine the PDF description of the light scattering statistics in a well characterized tissue-like particulate medium using optical coherence tomography (OCT). It is shown that for low scatterer density, the governing statistics depart considerably from a Gaussian description and follow the K distribution for both OCT amplitude and intensity. The PDF formalism is shown to be independent of the scatterer flow conditions; this is expected from theory, and suggests robustness and motion independence of the OCT amplitude (and OCT intensity) PDF metrics in the context of potential biomedical applications.

Age-related changes of the corneal speckle by Optical Coherence Tomography

Ageing of the population has several consequences for the human visual system. The corneal tissue which is responsible for two thirds of the total optical power is not an exception, suffering structural and anatomical changes with the age. In this study, we have explored the usefulness of Optical Coherence Tomography (OCT) speckle to study the age-related changes in corneal stroma. The Generalized Gamma distribution has been suggested as a model for OCT corneal speckle and was used here to investigate OCT envelope statistics in three different age groups (65 subjects in total). The scaling and the two shape parameters have shown significant statistical differences (p <; 0.05) among the age groups, showing the possibility to track the corneal changes using envelope statistics. This study shows that Generalized Gamma distribution has a great potential to fit the corneal speckle and could be a helpful complement to evaluate the structure and elastic properties of the cornea in vivo.

Optical coherence tomography spectral analysis for detecting apoptosis in vitro and in vivo

Apoptosis is a form of programmed cell death characterized by a series of predictable morphological changes at the subcellular level, which modify the light-scattering properties of cells. We present a spectroscopic optical coherence tomography (OCT) technique to detect changes in subcellular morphology related to apoptosis in vitro and in vivo. OCT data were acquired from acute myeloid leukemia (AML) cells treated with cisplatin over a 48-h period. The backscatter spectrum of the OCT signal acquired from the cell samples was characterized by calculating its in vitro integrated backscatter (IB) and spectral slope (SS). The IB increased with treatment duration, while the SS decreased, with the most significant changes occurring after 24 to 48 h of treatment. These changes coincided with striking morphological transformations in the cells and their nuclei. Similar trends in the spectral parameter values were observed in vivo in solid tumors grown from AML cells in mice, which were treated with chemotherapy and radiation. Our results provide a strong foundation from which future experiments may be designed to further understand the effect of cellular morphology and kinetics of apoptosis on the OCT signal and demonstrate the feasibility of using this technique in vivo.

Longitudinal label-free tracking of cell death dynamics in living engineered human skin tissue with a multimodal microscope

We demonstrate real-time, longitudinal, label-free tracking of apoptotic and necrotic cells in living tissue using a multimodal microscope. The integrated imaging platform combines multi-photon microscopy (MPM, based on two-photon excitation fluorescence), optical coherence microscopy (OCM), and fluorescence lifetime imaging microscopy (FLIM). Three-dimensional (3-D) co-registered images are captured that carry comprehensive information of the sample, including structural, molecular, and metabolic properties, based on light scattering, autofluorescence intensity, and autofluorescence lifetime, respectively. Different cell death processes, namely, apoptosis and necrosis, of keratinocytes from different epidermal layers are longitudinally monitored and investigated. Differentiation of the two cell death processes in a complex living tissue environment is enabled by quantitative image analysis and high-confidence classification processing based on the multidimensional, cross-validating imaging data. These results suggest that despite the limitations of each individual label-free modality, this multimodal imaging approach holds the promise for studies of different cell death processes in living tissue and in vivo organs.

Speckle statistics in OCT images: Monte Carlo simulations and experimental studies

The speckle pattern of an optical coherence tomography (OCT) image carries potentially useful sample information that may assist in tissue characterization. Recent biomedical results in vivo indicate that the distribution of signal intensities within an OCT tissue image is well described by a log-normal-like (Gamma) function. To fully understand and exploit this finding, an OCT Monte Carlo model that accounts for speckle effects was developed. The resultant Monte Carlo speckle statistics predictions agree well with experimental OCT results from a series of control phantoms with variable scattering properties; the Gamma distribution provides a good fit to the theoretical and experimental results. The ability to quantify subresolution tissue features via OCT speckle analysis may prove useful in diagnostic photomedicine.

In vitro detection of apoptosis using oscillating and pulsed gradient diffusion magnetic resonance imaging

Cellular apoptosis, a common pathway towards tumor regression, is induced by many radiotherapy and chemotherapy regimens. Imaging methods that can detect apoptosis may be able to assess treatment response earlier than typical tumor volume measurements. In this paper, a wide range of diffusion experiments and a simple model of diffusion in tissues were used to probe the microstructural effects of apoptosis. Experiments were conducted on acute myeloid leukemia cell pellets, where apoptosis was induced by treatment with the chemotherapeutic agent cisplatin. Seventy-two hours following cisplatin treatment, pulsed and oscillating gradient diffusion measurements were utilized to assess effects across a broad range of structural scales. The presence of apoptosis, which was histologically confirmed by TUNEL (terminal deoxynucleotidyl transferase UTP nick end labelling) staining, significantly changed diffusion properties. To describe these changes, the data were fit to the parallel plane model, which characterizes the effects of restricted diffusion in terms of three parameters: d, the restricted size, Dfree , the intrinsic, free diffusion coefficient of the solvent, and Drest , the long time or "restricted" diffusion coefficient. Apoptotic samples exhibited significant decreases in parameters d and Dfree and a significant increase in Drest . These changes appear consistent with the established morphological effects of apoptosis. In particular, the decrease in d may be a result of the combined effects of cell shrinkage, nuclear fragmentation and membrane blebbing, the decrease in Dfree may relate to cytosolic condensation, while the increase in Drest can be attributed to increased membrane permeability and extracellular volume fraction. By non-invasively detecting apoptosis, the methods reported in this study have the potential to improve upon current MRI methods for monitoring therapeutic response. Furthermore, these methods may offer sufficient specificity to differentiate between apoptosis and other modes of cell death, such as oncosis or necrosis.

Texture analysis of optical coherence tomography speckle for characterizing biological tissues in vivo

We demonstrate a method for differentiating tissue disease states using the intrinsic texture properties of speckle in optical coherence tomography (OCT) images of normal and tumor tissues obtained in vivo. This approach fits a gamma distribution function to the nonlog-compressed OCT image intensities, thus allowing differentiation of normal and tumor tissues in an ME-180 human cervical cancer mouse xenograft model. Quantitative speckle intensity distribution analysis thus shows promise for identifying tissue pathologies, with potential for early cancer detection in vivo.

Quantitative measurement of optical attenuation coefficients of cell lines CNE1, CNE2, and NP69 using optical coherence tomography

The radiotherapy-related types of nasopharyngeal carcinoma (NPC) have been established, which give the most effective treatment for NPC patients using the individual therapy. To diagnose the types of NPC, we assess the general NPC cell lines CNE1, CNE2 and normal nasopharyngeal cell line NP69 using optical coherence tomography (OCT) in two steps: firstly, the OCT images of the three different types of cell pellets are captured. Secondly, by fitting Beer's law to the averaged A-scans in these OCT datasets, the attenuation coefficients (μ t ) of the cells can be extracted. The median attenuation coefficients (interquartile range) of CNE1, CNE2, and NP69 are 5.58 mm(-1) (IQR 5.55 to 5.65 mm(-1)), 5.91 mm(-1) (IQR 5.82 to 5.88 mm(-1)), and 8.96 mm(-1) (IQR 8.80 to 9.47 mm(-1)), respectively. The distinguishable quantitative OCT analysis (by μ t ) shows that the types of NPC could potentially be differentiated in real time and noninvasive.

Multimodality molecular imaging of apoptosis in oncology

OBJECTIVE

The purposes of this review are to describe the signaling pathways of and the cellular changes that occur with apoptosis and other forms of cell death, summarize tracers and modalities used for imaging of apoptosis, delineate the relation between apoptosis and inhibition of protein translation, and describe spectroscopic technologies that entail high-frequency ultrasound and infrared and midinfrared light in characterizing the intracellular events of apoptosis.

CONCLUSION

Apoptosis is a highly orchestrated set of biochemical and morphologic cellular events. These events present many potential targets for the imaging of apoptosis in vivo. Imaging of apoptosis can facilitate early assessment of anticancer treatment before tumor shrinkage, which may increase the effectiveness of delivery of chemotherapy and radiation therapy and speed drug development.

Detecting apoptosis using dynamic light scattering with optical coherence tomography

A dynamic light scattering technique is implemented using optical coherence tomography (OCT) to measure the change in intracellular motion as cells undergo apoptosis. Acute myeloid leukemia cells were treated with cisplatin and imaged at a frame rate of 166 Hz using a 1300 nm swept-source OCT system at various times over a period of 48 h. Time correlation analysis of the speckle intensities indicated a significant increase in intracellular motion 24 h after treatment. This rise in intracellular motion correlated with histological findings of irregularly shaped and fragmented cells indicative of cell membrane blebbing and fragmentation.