Bright-field and fluorescence chip-scale microscopy for biological imaging

Optical microscopy is an essential tool in biological science and one of the gold standards for medical examinations. Miniaturization of microscopes can be a crucial stepping stone towards realizing compact, cost-effective and portable platforms for biomedical research and healthcare. This thesis reports on implementations of bright-field and fluorescence chip-scale microscopes for a variety of biological imaging applications. The term “chip-scale microscopy” refers to lensless imaging techniques realized in the form of mass-producible semiconductor devices, which transforms the fundamental design of optical microscopes.

Our strategy for chip-scale microscopy involves utilization of low-cost Complementary metal Oxide Semiconductor (CMOS) image sensors, computational image processing and micro-fabricated structural components. First, the sub-pixel resolving optofluidic microscope (SROFM), will be presented, which combines microfluidics and pixel super-resolution image reconstruction to perform high-throughput imaging of fluidic samples, such as blood cells. We discuss design parameters and construction of the device, as well as the resulting images and the resolution of the device, which was 0.66 µm at the highest acuity. The potential applications of SROFM for clinical diagnosis of malaria in the resource-limited settings is discussed.

Next, the implementations of ePetri, a self-imaging Petri dish platform with microscopy resolution, are presented. Here, we simply place the sample of interest on the surface of the image sensor and capture the direct shadow images under the illumination. By taking advantage of the inherent motion of the microorganisms, we achieve high resolution (~1 µm) imaging and long term culture of motile microorganisms over ultra large field-of-view (5.7 mm × 4.4 mm) in a specialized ePetri platform. We apply the pixel super-resolution reconstruction to a set of low-resolution shadow images of the microorganisms as they move across the sensing area of an image sensor chip and render an improved resolution image. We perform longitudinal study of Euglena gracilis cultured in an ePetri platform and image based analysis on the motion and morphology of the cells. The ePetri device for imaging non-motile cells are also demonstrated, by using the sweeping illumination of a light emitting diode (LED) matrix for pixel super-resolution reconstruction of sub-pixel shifted shadow images. Using this prototype device, we demonstrate the detection of waterborne parasites for the effective diagnosis of enteric parasite infection in resource-limited settings.

Then, we demonstrate the adaptation of a smartphone’s camera to function as a compact lensless microscope, which uses ambient illumination as its light source and does not require the incorporation of a dedicated light source. The method is also based on the image reconstruction with sweeping illumination technique, where the sequence of images are captured while the user is manually tilting the device around any ambient light source, such as the sun or a lamp. Image acquisition and reconstruction is performed on the device using a custom-built android application, constructing a stand-alone imaging device for field applications. We discuss the construction of the device using a commercial smartphone and demonstrate the imaging capabilities of our system.

Finally, we report on the implementation of fluorescence chip-scale microscope, based on a silo-filter structure fabricated on the pixel array of a CMOS image sensor. The extruded pixel design with metal walls between neighboring pixels successfully guides fluorescence emission through the thick absorptive filter to the photodiode layer of a pixel. Our silo-filter CMOS image sensor prototype achieves 13-µm resolution for fluorescence imaging over a wide field-of-view (4.8 mm × 4.4 mm). Here, we demonstrate bright-field and fluorescence longitudinal imaging of living cells in a compact, low-cost configuration.

Wide Field-of-view Microscopes and Endoscopes for Time-lapse Imaging and High-throughput Screening

Wide field-of-view (FOV) microscopy is of high importance to biological research and clinical diagnosis where a high-throughput screening of samples is needed. This thesis presents the development of several novel wide FOV imaging technologies and demonstrates their capabilities in longitudinal imaging of living organisms, on the scale of viral plaques to live cells and tissues.

The ePetri Dish is a wide FOV on-chip bright-field microscope. Here we applied an ePetri platform for plaque analysis of murine norovirus 1 (MNV-1). The ePetri offers the ability to dynamically track plaques at the individual cell death event level over a wide FOV of 6 mm × 4 mm at 30 min intervals. A density-based clustering algorithm is used to analyze the spatial-temporal distribution of cell death events to identify plaques at their earliest stages. We also demonstrate the capabilities of the ePetri in viral titer count and dynamically monitoring plaque formation, growth, and the influence of antiviral drugs.

We developed another wide FOV imaging technique, the Talbot microscope, for the fluorescence imaging of live cells. The Talbot microscope takes advantage of the Talbot effect and can generate a focal spot array to scan the fluorescence samples directly on-chip. It has a resolution of 1.2 μm and a FOV of ~13 mm². We further upgraded the Talbot microscope for the long-term time-lapse fluorescence imaging of live cell cultures, and analyzed the cells’ dynamic response to an anticancer drug.

We present two wide FOV endoscopes for tissue imaging, named the AnCam and the PanCam. The AnCam is based on the contact image sensor (CIS) technology, and can scan the whole anal canal within 10 seconds with a resolution of 89 μm, a maximum FOV of 100 mm × 120 mm, and a depth-of-field (DOF) of 0.65 mm. We also demonstrate the performance of the AnCam in whole anal canal imaging in both animal models and real patients. In addition to this, the PanCam is based on a smartphone platform integrated with a panoramic annular lens (PAL), and can capture a FOV of 18 mm × 120 mm in a single shot with a resolution of 100─140 μm. In this work we demonstrate the PanCam’s performance in imaging a stained tissue sample.

Myocardial edema after ischemia/reperfusion is not stable and follows a bimodal pattern: Imaging and histological tissue characterization

It is widely accepted that edema occurs early in the ischemic zone and persists in stable form for at least 1 week after myocardial ischemia/reperfusion. However, there are no longitudinal studies covering from very early (minutes) to late (1 week) reperfusion stages confirming this phenomenon. This study sought to perform a comprehensive longitudinal imaging and histological characterization of the edematous reaction after experimental myocardial ischemia/reperfusion. The study population consisted of 25 instrumented Large-White pigs (30 kg to 40 kg). Closed-chest 40-min ischemia/reperfusion was performed in 20 pigs, which were sacrificed at 120 min (n = 5), 24 h (n = 5), 4 days (n = 5), and 7 days (n = 5) after reperfusion and processed for histological quantification of myocardial water content. Cardiac magnetic resonance (CMR) scans with T2-weighted short-tau inversion recovery and T2-mapping sequences were performed at every follow-up stage until sacrifice. Five additional pigs sacrificed after baseline CMR served as controls. In all pigs, reperfusion was associated with a significant increase in T2 relaxation times in the ischemic region. On 24-h CMR, ischemic myocardium T2 times returned to normal values (similar to those seen pre-infarction). Thereafter, ischemic myocardium-T2 times in CMR performed on days 4 and 7 after reperfusion progressively and systematically increased. On day 7 CMR, T2 relaxation times were as high as those observed at reperfusion. Myocardial water content analysis in the ischemic region showed a parallel bimodal pattern: 2 high water content peaks at reperfusion and at day 7, and a significant decrease at 24 h. Contrary to the accepted view, myocardial edema during the first week after ischemia/reperfusion follows a bimodal pattern. The initial wave appears abruptly upon reperfusion and dissipates at 24 h. Conversely, the deferred wave of edema appears progressively days after ischemia/reperfusion and is maximal around day 7 after reperfusion.

An external field prior for the hidden Potts model with application to cone-beam computed tomography

In images with low contrast-to-noise ratio (CNR), the information gain from the observed pixel values can be insufficient to distinguish foreground objects. A Bayesian approach to this problem is to incorporate prior information about the objects into a statistical model. A method for representing spatial prior information as an external field in a hidden Potts model is introduced. This prior distribution over the latent pixel labels is a mixture of Gaussian fields, centred on the positions of the objects at a previous point in time. It is particularly applicable in longitudinal imaging studies, where the manual segmentation of one image can be used as a prior for automatic segmentation of subsequent images. The method is demonstrated by application to cone-beam computed tomography (CT), an imaging modality that exhibits distortions in pixel values due to X-ray scatter. The external field prior results in a substantial improvement in segmentation accuracy, reducing the mean pixel misclassification rate for an electron density phantom from 87% to 6%. The method is also applied to radiotherapy patient data, demonstrating how to derive the external field prior in a clinical context.

Bayesian computational methods for spatial analysis of images

This thesis introduces a new way of using prior information in a spatial model and develops scalable algorithms for fitting this model to large imaging datasets. These methods are employed for image-guided radiation therapy and satellite based classification of land use and water quality. This study has utilized a pre-computation step to achieve a hundredfold improvement in the elapsed runtime for model fitting. This makes it much more feasible to apply these models to real-world problems, and enables full Bayesian inference for images with a million or more pixels.

Mapping hippocampal and ventricular change in Alzheimer disease

We developed an anatomical mapping technique to detect hippocampal and ventricular changes in Alzheimer disease (AD). The resulting maps are sensitive to longitudinal changes in brain structure as the disease progresses. An anatomical surface modeling approach was combined with surface-based statistics to visualize the region and rate of atrophy in serial MRI scans and isolate where these changes link with cognitive decline. Fifty-two high-resolution MRI scans were acquired from 12 AD patients (age: 68.4 ± 1.9 years) and 14 matched controls (age: 71.4 ± 0.9 years), each scanned twice (2.1 ± 0.4 years apart). 3D parametric mesh models of the hippocampus and temporal horns were created in sequential scans and averaged across subjects to identify systematic patterns of atrophy. As an index of radial atrophy, 3D distance fields were generated relating each anatomical surface point to a medial curve threading down the medial axis of each structure. Hippocampal atrophic rates and ventricular expansion were assessed statistically using surface-based permutation testing and were faster in AD than in controls. Using color-coded maps and video sequences, these changes were visualized as they progressed anatomically over time. Additional maps localized regions where atrophic changes linked with cognitive decline. Temporal horn expansion maps were more sensitive to AD progression than maps of hippocampal atrophy, but both maps correlated with clinical deterioration. These quantitative, dynamic visualizations of hippocampal atrophy and ventricular expansion rates in aging and AD may provide a promising measure to track AD progression in drug trials.

Mapping hippocampal and ventricular change in Alzheimer disease

We developed an anatomical mapping technique to detect hippocampal and ventricular changes in Alzheimer disease (AD). The resulting maps are sensitive to longitudinal changes in brain structure as the disease progresses. An anatomical surface modeling approach was combined with surface-based statistics to visualize the region and rate of atrophy in serial MRI scans and isolate where these changes link with cognitive decline. Fifty-two high-resolution MRI scans were acquired from 12 AD patients (age: 68.4 +/- 1.9 years) and 14 matched controls (age: 71.4 +/- 0.9 years), each scanned twice (2.1 +/- 0.4 years apart). 3D parametric mesh models of the hippocampus and temporal horns were created in sequential scans and averaged across subjects to identify systematic patterns of atrophy. As an index of radial atrophy, 3D distance fields were generated relating each anatomical surface point to a medial curve threading down the medial axis of each structure. Hippocampal atrophic rates and ventricular expansion were assessed statistically using surface-based permutation testing and were faster in AD than in controls. Using color-coded maps and video sequences, these changes were visualized as they progressed anatomically over time. Additional maps localized regions where atrophic changes linked with cognitive decline. Temporal horn expansion maps were more sensitive to AD progression than maps of hippocampal atrophy, but both maps correlated with clinical deterioration. These quantitative, dynamic visualizations of hippocampal atrophy and ventricular expansion rates in aging and AD may provide a promising measure to track AD progression in drug trials. (C) 2004 Elsevier Inc. All rights reserved.

Association between biomechanical structural stresses of atherosclerotic carotid plaques and subsequent ischaemic cerebrovascular events - A longitudinal in vivo magnetic resonance imaging-based finite element study

Background: High-resolution magnetic resonance (MR) imaging has been used for MR imaging-based structural stress analysis of atherosclerotic plaques. The biomechanical stress profile of stable plaques has been observed to differ from that of unstable plaques; however, the role that structural stresses play in determining plaque vulnerability remains speculative. Methods: A total of 61 patients with previous history of symptomatic carotid artery disease underwent carotid plaque MR imaging. Plaque components of the index artery such as fibrous tissue, lipid content and plaque haemorrhage (PH) were delineated and used for finite element analysis-based maximum structural stress (M-C Stress) quantification. These patients were followed up for 2 years. The clinical end point was occurrence of an ischaemic cerebrovascular event. The association of the time to the clinical end point with plaque morphology and M-C Stress was analysed. Results: During a median follow-up duration of 514 days, 20% of patients (n=12) experienced an ischaemic event in the territory of the index carotid artery. Cox regression analysis indicated that M-C Stress (hazard ratio (HR): 12.98 (95% confidence interval (CI): 1.32-26.67, pZ0.02), fibrous cap (FC) disruption (HR: 7.39 (95% CI: 1.61e33.82), p Z 0.009) and PH (HR: 5.85 (95% CI: 1.27e26.77), p Z 0.02) are associated with the development of subsequent cerebrovascular events. Plaques associated with future events had higher M-C Stress than those which had remained asymptomatic (median (interquartile range, IQR): 330 kPa (229e494) vs. 254 kPa (166-290), p Z0.04). Conclusions: High biomechanical structural stresses, in addition to FC rupture and PH, are associated with subsequent cerebrovascular events.

Spatial longitudinal coherence length of a thermal source and its influence on lensless ghost imaging

The spatial longitudinal coherence length (SLCL), which is determined by the size of and the distance from the source, is introduced to investigate the longitudinal resolution of lensless ghost imaging. Its influence is discussed quantitatively by simulation. The discrepancy of position sensitivity between Scareelli et al. [Appl. Phys. Lett. 88, 061106 (2006)] and Basano and Ottonello [Appl. Phys. Lett. 88, 091109 (2006)] is clarified. (C) 2008 Optical Society of America.

Test-retest variability in structural and functional parameters of glaucoma damage in the glaucoma imaging longitudinal study.

PURPOSE:

To determine the test-retest variability in perimetric, optic disc, and macular thickness parameters in a cohort of treated patients with established glaucoma.

PATIENTS AND METHODS:

In this cohort study, the authors analyzed the imaging studies and visual field tests at the baseline and 6-month visits of 162 eyes of 162 participant in the Glaucoma Imaging Longitudinal Study (GILS). They assessed the difference, expressed as the standard error of measurement, of Humphrey field analyzer II (HFA) Swedish Interactive Threshold Algorithm fast, Heidelberg retinal tomograph (HRT) II, and retinal thickness analyzer (RTA) parameters between the two visits and assumed that this difference was due to measurement variability, not pathologic change. A statistically significant change was defined as twice the standard error of measurement.

RESULTS:

In this cohort of treated glaucoma patients, it was found that statistically significant changes were 3.2 dB for mean deviation (MD), 2.2 for pattern standard deviation (PSD), 0.12 for cup shape measure, 0.26 mm for rim area, and 32.8 microm and 31.8 microm for superior and inferior macular thickness, respectively. On the basis of these values, it was estimated that the number of potential progression events detectable in this cohort by the parameters of MD, PSD, cup shape measure, rim area, superior macular thickness, and inferior macular thickness was 7.5, 6.0, 2.3, 5.7, 3.1, and 3.4, respectively.

CONCLUSIONS:

The variability of the measurements of MD, PSD, and rim area, relative to the range of possible values, is less than the variability of cup shape measure or macular thickness measurements. Therefore, the former measurements may be more useful global measurements for assessing progressive glaucoma damage.

Living Renal Allograft Transplantation: Diffusion-weighted MR Imaging in Longitudinal Follow-up of the Donated and the Remaining Kidney.

Purpose To determine whether diffusion-weighted (DW) magnetic resonance (MR) imaging in living renal allograft donation allows monitoring of potential changes in the nontransplanted remaining kidney of the donor because of unilateral nephrectomy and changes in the transplanted kidney before and after transplantation in donor and recipient, respectively, and whether DW MR parameters are correlated in the same kidney before and after transplantation. Materials and Methods The study protocol was approved by the local ethics committee; written informed consent was obtained. Thirteen healthy kidney donors and their corresponding recipients prospectively underwent DW MR imaging (multiple b values) in donors before donation and in donors and recipients at day 8 and months 3 and 12 after donation. Total apparent diffusion coefficient (ADCT) values were determined; contribution of microcirculation was quantified in perfusion fraction (FP). Longitudinal changes of diffusion parameters were compared (repeated-measures one-way analysis of variance with post hoc pairwise comparisons). Correlations were tested (linear regression). Results ADCT values in nontransplanted kidney of donors increased from a preexplantation value of (188 ± 9 [standard deviation]) to (202 ± 11) × 10(-5) mm(2)/sec in medulla and from (199 ± 11) to (210 ± 13) × 10(-5) mm(2)/sec in cortex 1 week after donation (P < .004). Medullary, but not cortical, ADCT values stayed increased up to 1 year. ADCT values in allografts in recipients were stable. Compared with values obtained before transplantation in donors, the corticomedullary difference was reduced in allografts (P < .03). Cortical ADCT values correlated with estimated glomerular filtration rate in recipients (R = 0.56, P < .001) but not donors. Cortical ADCT values in the same kidney before transplantation in donors correlated with those in recipients on day 8 after transplantation (R = 0.77, P = .006). FP did not show significant changes. Conclusion DW MR imaging depicts early adaptations in the remaining nontransplanted kidney of donors after nephrectomy. All diffusion parameters remained constant in allograft recipients after transplantation. This method has potential monitoring utility, although assessment of clinical relevance is needed. © RSNA, 2013 Online supplemental material is available for this article.

Renal Blood Oxygenation Level-dependent Imaging in Longitudinal Follow-up of Donated and Remaining Kidneys

Purpose To determine renal oxygenation changes associated with uninephrectomy and transplantation in both native donor kidneys and transplanted kidneys by using blood oxygenation level-dependent (BOLD) MR imaging. Materials and Methods The study protocol was approved by the local ethics committee. Thirteen healthy kidney donors and their corresponding recipients underwent kidney BOLD MR imaging with a 3-T imager. Written informed consent was obtained from each subject. BOLD MR imaging was performed in donors before uninephrectomy and in donors and recipients 8 days, 3 months, and 12 months after transplantation. R2* values, which are inversely related to tissue partial pressure of oxygen, were determined in the cortex and medulla. Longitudinal R2* changes were statistically analyzed by using repeated measures one-way analysis of variance with post hoc pair-wise comparisons. Results R2* values in the remaining kidneys significantly decreased early after uninephrectomy in both the medulla and cortex (P < .003), from 28.9 sec(-1) ± 2.3 to 26.4 sec(-1) ± 2.5 in the medulla and from 18.3 sec(-1) ± 1.5 to 16.3 sec(-1) ± 1.0 in the cortex, indicating increased oxygen content. In donors, R2* remained significantly decreased in both the medulla and cortex at 3 (P < .01) and 12 (P < .01) months. In transplanted kidneys, R2* remained stable during the first year after transplantation, with no significant change. Among donors, cortical R2* was found to be negatively correlated with estimated glomerular filtration rate (R = -0.47, P < .001). Conclusion The results suggest that BOLD MR imaging may potentially be used to monitor renal functional changes in both remaining and corresponding transplanted kidneys. (©) RSNA, 2016.

Longitudinal diffusion tensor imaging in dementia with Lewy bodies and Alzheimer's disease

Copyright © 2016 Elsevier Ltd. All rights reserved. Acknowledgements The study was supported by the NIHR Biomedical Research Unit in Dementia and the Biomedical Research Centre awarded to Cambridge University Hospitals NHS Foundation Trust and the University of Cambridge, and the NIHR Biomedical Research Unit in Dementia and the Biomedical Research Centre awarded to Newcastle upon Tyne Hospitals NHS Foundation Trust and Newcastle University. Elijah Mak was in receipt of a Gates Cambridge PhD studentship.

Anatomical MR imaging of long bones : comparative performance of MRI at 1.5 T and 3 T

The current gold standard for the design of orthopaedic implants is 3D models of long bones obtained using computed tomography (CT). However, high-resolution CT imaging involves high radiation exposure, which limits its use in healthy human volunteers. Magnetic resonance imaging (MRI) is an attractive alternative for the scanning of healthy human volunteers for research purposes. Current limitations of MRI include difficulties of tissue segmentation within joints and long scanning times. In this work, we explore the possibility of overcoming these limitations through the use of MRI scanners operating at a higher field strength. We quantitatively compare the quality of anatomical MR images of long bones obtained at 1.5 T and 3 T and optimise the scanning protocol of 3 T MRI. FLASH images of the right leg of five human volunteers acquired at 1.5 T and 3 T were compared in terms of signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). The comparison showed a relatively high CNR and SNR at 3 T for most regions of the femur and tibia, with the exception of the distal diaphyseal region of the femur and the mid diaphyseal region of the tibia. This was accompanied by an ~65% increase in the longitudinal spin relaxation time (T1) of the muscle at 3 T compared to 1.5 T. The results suggest that MRI at 3 T may be able to enhance the segmentability and potentially improve the accuracy of 3D anatomical models of long bones, compared to 1.5 T. We discuss how the total imaging times at 3 T can be kept short while maximising the CNR and SNR of the images obtained.