Incorporation of an experimentally determined MTF for spatial frequency filtering and deconvolution during optical projection tomography reconstruction

Challenges and advances in optical 3D mesoscale imaging

Optical mesoscale imaging is a rapidly developing field that allows the visualisation of larger samples than is possible with standard light microscopy, and fills a gap between cell and organism resolution. It spans from advanced fluorescence imaging of micrometric cell clusters to centimetre-size complete organisms. However, with larger volume specimens, new problems arise. Imaging deeper into tissues at high resolution poses challenges ranging from optical distortions to shadowing from opaque structures. This manuscript discusses the latest developments in mesoscale imaging and highlights limitations, namely labelling, clearing, absorption, scattering, and also sample handling. We then focus on approaches that seek to turn mesoscale imaging into a more quantitative technique, analogous to quantitative tomography in medical imaging, highlighting a future role for digital and physical phantoms as well as artificial intelligence.

High resolution optical projection tomography platform for multispectral imaging of the mouse gut

Optical projection tomography (OPT) is a powerful tool for three-dimensional imaging of mesoscopic biological samples with great use for biomedical phenotyping studies. We present a fluorescent OPT platform that enables direct visualization of biological specimens and processes at a centimeter scale with high spatial resolution, as well as fast data throughput and reconstruction. We demonstrate nearly isotropic sub-28 µm resolution over more than 60 mm3 after reconstruction of a single acquisition. Our setup is optimized for imaging the mouse gut at multiple wavelengths. Thanks to a new sample preparation protocol specifically developed for gut specimens, we can observe the spatial arrangement of the intestinal villi and the vasculature network of a 3-cm long healthy mouse gut. Besides the blood vessel network surrounding the gastrointestinal tract, we observe traces of vasculature at the villi ends close to the lumen. The combination of rapid acquisition and a large field of view with high spatial resolution in 3D mesoscopic imaging holds an invaluable potential for gastrointestinal pathology research.

Hyperspectral scanning laser optical tomography

In order to study physical relationships within tissue volumes or even organism-level systems, the spatial distribution of multiple fluorescent markers needs to be resolved efficiently in three dimensions. Here, rather than acquiring discrete spectral images sequentially using multiple emission filters, a hyperspectral scanning laser optical tomography system is developed to obtain hyperspectral volumetric data sets with 2-nm spectral resolution of optically transparent mesoscopic (millimeter-centimeter) specimens. This is achieved by acquiring a series of point-scanning hyperspectral extended depth of field images at different angles and subsequently tomographically reconstructing the 3D intensity distribution for each wavelength. This technique is demonstrated to provide robust measurements via the comparison of spectral and intensity profiles of fluorescent bead phantoms. Due to its enhanced spectral resolving ability, this technique is also demonstrated to resolve largely overlapping fluorophores, as demonstrated by the 3D fluorescence hyperspectral reconstruction of a dual-labeled mouse thymus gland sample and the ability to distinguish tumorous and normal tissues of an unlabeled mouse intestine sample.

Direct inversion algorithm for focal plane scanning optical projection tomography

To achieve approximately parallel projection geometry, traditional optical projection tomography (OPT) requires the use of low numerical aperture (NA) objectives, which have a long depth-of-field at the expense of poor lateral resolution. Particularly promising methods to improve spatial resolution include ad-hoc post-processing filters that limit the effect of the system's MTF and focal-plane-scanning OPT (FPS-OPT), an alternative acquisition procedure that allows the use of higher NA objectives by limiting the effect of their shallow depth of field yet still assumes parallel projection rays during reconstruction. Here, we provide a detailed derivation that establishes the existence of a direct inversion formula for FPS-OPT. Based on this formula, we propose a point spread function-aware algorithm that is similar in form and complexity to traditional filtered backprojection (FBP). With simulations, we demonstrate that our point-spread-function aware FBP for FPS-OPT leads to more accurate images than both traditional OPT with deconvolution and FPS-OPT with naive FBP reconstruction. We further illustrate the technique on experimental zebrafish data, which shows that our approach reduces out-of-focus blur compared to a direct FBP reconstruction with FPS-OPT.

Point spread function based image reconstruction in optical projection tomography

As a result of the shallow depth of focus of the optical imaging system, the use of standard filtered back projection in optical projection tomography causes space-variant tangential blurring that increases with the distance to the rotation axis. We present a novel optical tomographic image reconstruction technique that incorporates the point spread function of the imaging lens in an iterative reconstruction. The technique is demonstrated using numerical simulations, tested on experimental optical projection tomography data of single fluorescent beads, and applied to high-resolution emission optical projection tomography imaging of an entire zebrafish larva. Compared to filtered back projection our results show greatly reduced radial and tangential blurring over the entire [Formula: see text] mm² field of view, and a significantly improved signal to noise ratio.

Inverse filtering approach to measure directional in-plane modulation transfer function using a sphere phantom for a digital tomosynthesis system

We propose a method to measure the directional in-plane modulation transfer function (MTF) of a digital tomosynthesis system using a sphere phantom. To assess the spatial resolution of an in-plane image of the tomosynthesis system, projection data of a sphere phantom were generated within a limited data acquisition range of 40°, and reconstructed by the FDK algorithm. To measure the in-plane MTF, we divided the Fourier transform of the reconstructed sphere phantom by that of the ideal sphere phantom, and then performed plane integral along the fz-direction. When dividing, small values in the denominator can introduce estimation errors, and these errors were reduced by the proposed method. To evaluate the performance of the proposed method, the in-plane MTF estimated by simulation and experimental data was compared to the ideal in-plane MTF generated by computer simulations using a point object. For quantitative evaluation, we measured frequency values at half-maximum and full-maximum of the directional in-plane MTF along the three different directions (i.e., f_0° -, f_30° -, and f_60° -directions) and compared them with those of the ideal in-plane MTF. Although the sphere phantom has been regarded as an inappropriate object due to the anisotropic characteristics of tomosynthesis image, our results show that the proposed method has a reliable estimation performance, demonstrating the sphere phantom is still suitable for measuring the directional in-plane MTF for a digital tomosynthesis system.

Image resolution and deconvolution in optical tomography

We present a frequency domain analysis of the image resolution of optical tomography systems. The result of our analysis is a description of the spatially-variant resolution in optical tomographic image after reconstruction as a function of the properties of the imaging system geometry. We validate our model using optical projection tomography (OPT) measurements of fluorescent beads embedded in agarose gel. Our model correctly describes both the radial and tangential resolution of the measured images. In addition, we present a correction of the tomographic images for the spatially-varying resolution using a deconvolution algorithm. The resulting corrected tomographic reconstruction shows a homogeneous and isotropic pixel-limited resolution across the entire image. Our method is applied to OPT measurements of a zebrafish, showing improved resolution. Aside from allowing image correction and providing a resolution measure for OPT systems, our model provides a powerful tool for the design of optical tomographic systems.

Applicability, usability, and limitations of murine embryonic imaging with optical coherence tomography and optical projection tomography

We present an analysis of imaging murine embryos at various embryonic developmental stages (embryonic day 9.5, 11.5, and 13.5) by optical coherence tomography (OCT) and optical projection tomography (OPT). We demonstrate that while OCT was capable of rapid high-resolution live 3D imaging, its limited penetration depth prevented visualization of deeper structures, particularly in later stage embryos. In contrast, OPT was able to image the whole embryos, but could not be used in vivo because the embryos must be fixed and cleared. Moreover, the fixation process significantly altered the embryo morphology, which was quantified by the volume of the eye-globes before and after fixation. All of these factors should be weighed when determining which imaging modality one should use to achieve particular goals of a study.

Visualising apoptosis in live zebrafish using fluorescence lifetime imaging with optical projection tomography to map FRET biosensor activity in space and time

Fluorescence lifetime imaging (FLIM) combined with optical projection tomography (OPT) has the potential to map Förster resonant energy transfer (FRET) readouts in space and time in intact transparent or near transparent live organisms such as zebrafish larvae, thereby providing a means to visualise cell signalling processes in their physiological context. Here the first application of FLIM OPT to read out biological function in live transgenic zebrafish larvae using a genetically expressed FRET biosensor is reported. Apoptosis, or programmed cell death, is mapped in 3-D by imaging the activity of a FRET biosensor that is cleaved by Caspase 3, which is a key effector of apoptosis. Although apoptosis is a naturally occurring process during development, it can also be triggered in a variety of ways, including through gamma irradiation. FLIM OPT is shown here to enable apoptosis to be monitored over time, in live zebrafish larvae via changes in Caspase 3 activation following gamma irradiation at 24 hours post fertilisation. Significant apoptosis was observed at 3.5 hours post irradiation, predominantly in the head region.

Assessment of optical CT as a future QA tool for synchrotron x-ray microbeam therapy

Synchrotron microbeam radiation therapy (MRT) is an advanced form of radiotherapy for which it is extremely difficult to provide adequate quality assurance. This may delay or limit its clinical uptake, particularly in the paediatric patient populations for whom it could be especially suitable. This study investigates the extent to which new developments in 3D dosimetry using optical computed tomography (CT) can visualise MRT dose distributions, and assesses what further developments are necessary before fully quantitative 3D measurements can be achieved. Two experiments are reported. In the first cylindrical samples of the radiochromic polymer PRESAGE(®) were irradiated with different complex MRT geometries including multiport treatments of collimated 'pencil' beams, interlaced microplanar arrays and a multiport treatment using an anthropomorphic head phantom. Samples were scanned using transmission optical CT. In the second experiment, optical CT measurements of the biologically important peak-to-valley dose ratio (PVDR) were compared with expected values from Monte Carlo simulations. The depth-of-field (DOF) of the optical CT system was characterised using a knife-edge method and the possibility of spatial resolution improvement through deconvolution of a measured point spread function (PSF) was investigated. 3D datasets from the first experiment revealed excellent visualisation of the 50 μm beams and various discrepancies from the planned delivery dose were found. The optical CT PVDR measurements were found to be consistently 30% of the expected Monte Carlo values and deconvolution of the microbeam profiles was found to lead to increased noise. The reason for the underestimation of the PVDR by optical CT was attributed to lack of spatial resolution, supported by the results of the DOF characterisation. Solutions are suggested for the outstanding challenges and the data are shown already to be useful in identifying potential treatment anomalies.

Measurement of the modulation transfer function of an X-ray microscope based on multiple Fourier orders analysis of a Siemens star

Soft X-ray tomography (SXT) is becoming a powerful imaging technique to analyze eukaryotic whole cells close to their native state. Central to the analysis of the quality of SXT 3D reconstruction is the estimation of the spatial resolution and Depth of Field of the X-ray microscope. In turn, the characterization of the Modulation Transfer Function (MTF) of the optical system is key to calculate both parameters. Consequently, in this work we introduce a fully automated technique to accurately estimate the transfer function of such an optical system. Our proposal is based on the preprocessing of the experimental images to obtain an estimate of the input pattern, followed by the analysis in Fourier space of multiple orders of a Siemens Star test sample, extending in this way its measured frequency range.

Helical optical projection tomography

A new technique termed Helical Optical Projection Tomography (hOPT) has been developed with the aim to overcome some of the limitations of current 3D optical imaging techniques. hOPT is based on Optical Projection Tomography (OPT) with the major difference that there is a translation of the sample in the vertical direction during the image acquisition process, requiring a new approach to image reconstruction. Contrary to OPT, hOPT makes possible to obtain 3D-optical images of intact long samples without imposing limits on the sample length. This has been tested using hOPT to image long murine tissue samples such as spinal cords and large intestines. Moreover, 3D-reconstructed images of the colon of DSS-treated mice, a model for Inflammatory Bowel Disease, allowed the identification of the structural alterations. Finally, the geometry of the hOPT device facilitates the addition of a Selective Plane Illumination Microscopy (SPIM) arm, providing the possibility of delivering high resolution images of selected areas together with complete volumetric information.

Simultaneous angular multiplexing optical projection tomography at shifted focal planes

We describe an angular multiplexing technique for optical projection tomography that improves resolution, signal-to-noise ratio, and imaging speed by ameliorating the trade-off between spatial resolution and depth of field and improving the light collection efficiency. Here we demonstrate that imaging at two orthogonal angular projections simultaneously, focused on shifted planes in the sample, improves the average spatial resolution by ~20% and the light collection efficiency by a factor of ~4, thereby enabling increased acquisition speed and reduced light dose.