Content-aware approach for improving biomedical image analysis: an interdisciplinary study series

Biomedicine is a highly interdisciplinary research area at the interface of sciences, anatomy, physiology, and medicine. In the last decade, biomedical studies have been greatly enhanced by the introduction of new technologies and techniques for automated quantitative imaging, thus considerably advancing the possibility to investigate biological phenomena through image analysis. However, the effectiveness of this interdisciplinary approach is bounded by the limited knowledge that a biologist and a computer scientist, by professional training, have of each other’s fields. The possible solution to make up for both these lacks lies in training biologists to make them interdisciplinary researchers able to develop dedicated image processing and analysis tools by exploiting a content-aware approach. The aim of this Thesis is to show the effectiveness of a content-aware approach to automated quantitative imaging, by its application to different biomedical studies, with the secondary desirable purpose of motivating researchers to invest in interdisciplinarity. Such content-aware approach has been applied firstly to the phenomization of tumour cell response to stress by confocal fluorescent imaging, and secondly, to the texture analysis of trabecular bone microarchitecture in micro-CT scans. Third, this approach served the characterization of new 3-D multicellular spheroids of human stem cells, and the investigation of the role of the Nogo-A protein in tooth innervation. Finally, the content-aware approach also prompted to the development of two novel methods for local image analysis and colocalization quantification. In conclusion, the content-aware approach has proved its benefit through building new approaches that have improved the quality of image analysis, strengthening the statistical significance to allow unveiling biological phenomena. Hopefully, this Thesis will contribute to inspire researchers to striving hard for pursuing interdisciplinarity.

Quantitative Proteomic Analysis Reveals Metabolic Alterations, Calcium Dysregulation, And Increased Expression Of Extracellular Matrix Proteins In Laminin α2 Chain-deficient Muscle.

Congenital muscular dystrophy with laminin α2 chain deficiency (MDC1A) is one of the most severe forms of muscular disease and is characterized by severe muscle weakness and delayed motor milestones. The genetic basis of MDC1A is well known, yet the secondary mechanisms ultimately leading to muscle degeneration and subsequent connective tissue infiltration are not fully understood. In order to obtain new insights into the molecular mechanisms underlying MDC1A, we performed a comparative proteomic analysis of affected muscles (diaphragm and gastrocnemius) from laminin α2 chain-deficient dy(3K)/dy(3K) mice, using multidimensional protein identification technology combined with tandem mass tags. Out of the approximately 700 identified proteins, 113 and 101 proteins, respectively, were differentially expressed in the diseased gastrocnemius and diaphragm muscles compared with normal muscles. A large portion of these proteins are involved in different metabolic processes, bind calcium, or are expressed in the extracellular matrix. Our findings suggest that metabolic alterations and calcium dysregulation could be novel mechanisms that underlie MDC1A and might be targets that should be explored for therapy. Also, detailed knowledge of the composition of fibrotic tissue, rich in extracellular matrix proteins, in laminin α2 chain-deficient muscle might help in the design of future anti-fibrotic treatments. All MS data have been deposited in the ProteomeXchange with identifier PXD000978 (http://proteomecentral.proteomexchange.org/dataset/PXD000978).

Evaluation of follicular lymphoid depletion in the Bursa of Fabricius: an alternative methodology using digital image analysis and artificial neural networks

Fifty Bursa of Fabricius (BF) were examined by conventional optical microscopy and digital images were acquired and processed using Matlab® 6.5 software. The Artificial Neuronal Network (ANN) was generated using Neuroshell® Classifier software and the optical and digital data were compared. The ANN was able to make a comparable classification of digital and optical scores. The use of ANN was able to classify correctly the majority of the follicles, reaching sensibility and specificity of 89% and 96%, respectively. When the follicles were scored and grouped in a binary fashion the sensibility increased to 90% and obtained the maximum value for the specificity of 92%. These results demonstrate that the use of digital image analysis and ANN is a useful tool for the pathological classification of the BF lymphoid depletion. In addition it provides objective results that allow measuring the dimension of the error in the diagnosis and classification therefore making comparison between databases feasible.

Chrysotile effects on human lung cell carcinoma in culture: 3-D reconstruction and DNA quantification by image analysis

Background: Chrysotile is considered less harmful to human health than other types of asbestos fibers. Its clearance from the lung is faster and, in comparison to amphibole forms of asbestos, chrysotile asbestos fail to accumulate in the lung tissue due to a mechanism involving fibers fragmentation in short pieces. Short exposure to chrysotile has not been associated with any histopathological alteration of lung tissue. Methods: The present work focuses on the association of small chrysotile fibers with interphasic and mitotic human lung cancer cells in culture, using for analyses confocal laser scanning microscopy and 3D reconstructions. The main goal was to perform the analysis of abnormalities in mitosis of fibers-containing cells as well as to quantify nuclear DNA content of treated cells during their recovery in fiber-free culture medium. Results: HK2 cells treated with chrysotile for 48 h and recovered in additional periods of 24, 48 and 72 h in normal medium showed increased frequency of multinucleated and apoptotic cells. DNA ploidy of the cells submitted to the same chrysotile treatment schedules showed enhanced aneuploidy values. The results were consistent with the high frequency of multipolar spindles observed and with the presence of fibers in the intercellular bridge during cytokinesis. Conclusion: The present data show that 48 h chrysotile exposure can cause centrosome amplification, apoptosis and aneuploid cell formation even when long periods of recovery were provided. Internalized fibers seem to interact with the chromatin during mitosis, and they could also interfere in cytokinesis, leading to cytokinesis failure which forms aneuploid or multinucleated cells with centrosome amplification.

Characterization of micro-bubble size distribution and flow configuration in DAF contact zone by a non-intrusive image analysis system and tracer tests

This paper aims to investigate the influence of some dissolved air flotation (DAF) process variables (specifically: the hydraulic detention time in the contact zone and the supplied dissolved air concentration) and the pH values, as pretreatment chemical variables, on the micro-bubble size distribution (BSD) in a DAF contact zone. This work was carried out in a pilot plant where bubbles were measured by an appropriate non-intrusive image acquisition system. The results show that the obtained diameter ranges were in agreement with values reported in the literature (10-100mm), quite independently of the investigated conditions. The linear average diameter varied from 20 to 30mm, or equivalently, the Sauter (d(3,2)) diameter varied from 40 to 50mm. In all investigated conditions, D(50) was between 75% and 95%. The BSD might present different profile (with a bimodal curve trend), however, when analyzing the volumetric frequency distribution (in some cases with the appearance of peaks in diameters ranging from 90-100mm). Regarding volumetric frequency analysis, all the investigated parameters can modify the BSD in DAF contact zone after the release point, thus potentially causing changes in DAF kinetics. This finding prompts further research in order to verify the effect of these BSD changes on solid particle removal efficiency by DAF.

Parallel processing and image analysis in the eyes of mantis shrimps

The compound eyes of mantis shrimps, a group of tropical marine crustaceans, incorporate principles of serial and parallel processing of visual information that may be applicable to artificial imaging systems. Their eyes include numerous specializations for analysis of the spectral and polarizational properties of light, and include more photoreceptor classes for analysis of ultraviolet light, color, and polarization than occur in any other known visual system. This is possible because receptors in different regions of the eye are anatomically diverse and incorporate unusual structural features, such as spectral filters, not seen in other compound eyes. Unlike eyes of most other animals, eyes of mantis shrimps must move to acquire some types of visual information and to integrate color and polarization with spatial vision. Information leaving the retina appears to be processed into numerous parallel data streams leading into the central nervous system, greatly reducing the analytical requirements at higher levels. Many of these unusual features of mantis shrimp vision may inspire new sensor designs for machine vision

3D reconstruction and quantification of macropores using X-ray computed tomography and image analysis

Axial X-ray Computed tomography (CT) scanning provides a convenient means of recording the three-dimensional form of soil structure. The technique has been used for nearly two decades, but initial development has concentrated on qualitative description of images. More recently, increasing effort has been put into quantifying the geometry and topology of macropores likely to contribute to preferential now in soils. Here we describe a novel technique for tracing connected macropores in the CT scans. After object extraction, three-dimensional mathematical morphological filters are applied to quantify the reconstructed structure. These filters consist of sequences of so-called erosions and/or dilations of a 32-face structuring element to describe object distances and volumes of influence. The tracing and quantification methodologies were tested on a set of undisturbed soil cores collected in a Swiss pre-alpine meadow, where a new earthworm species (Aporrectodea nocturna) was accidentally introduced. Given the reduced number of samples analysed in this study, the results presented only illustrate the potential of the method to reconstruct and quantify macropores. Our results suggest that the introduction of the new species induced very limited chance to the soil structured for example, no difference in total macropore length or mean diameter was observed. However. in the zone colonised by, the new species. individual macropores tended to have a longer average length. be more vertical and be further apart at some depth. Overall, the approach proved well suited to the analysis of the three-dimensional architecture of macropores. It provides a framework for the analysis of complex structures, which are less satisfactorily observed and described using 2D imaging. (C) 2002 Elsevier Science B.V. All rights reserved.

Characterisation of the three-dimensional structure of earthworm burrow systems using image analysis and mathematical morphology

The aim of this work was to exemplify the specific contribution of both two- and three-dimensional (31)) X-ray computed tomography to characterise earthworm burrow systems. To achieve this purpose we used 3D mathematical morphology operators to characterise burrow systems resulting from the activity of an anecic (Aporrectodea noctunia), and an endogeic species (Allolobophora chlorotica), when both species were introduced either separately or together into artificial soil cores. Images of these soil cores were obtained using a medical X-ray tomography scanner. Three-dimensional reconstructions of burrow systems were obtained using a specifically developed segmentation algorithm. To study the differences between burrow systems, a set of classical tools of mathematical morphology (granulometries) were used. So-called granulometries based on different structuring elements clearly separated the different burrow systems. They enabled us to show that burrows made by the anecic species were fatter, longer, more vertical, more continuous but less sinuous than burrows of the endogeic species. The granulometry transform of the soil matrix showed that burrows made by A. nocturna were more evenly distributed than those of A. chlorotica. Although a good discrimination was possible when only one species was introduced into the soil cores, it was not possible to separate burrows of the two species from each other in cases where species were introduced into the same soil core. This limitation, partly due to the insufficient spatial resolution of the medical scanner, precluded the use of the morphological operators to study putative interactions between the two species.

Digital image analysis to standardize a photometric method in colorimetric quantification

Techniques applying digital images increasingly have been used in biology, medicine, physics, and other research areas. The image coordinates can represent light intensities values to be detected by a CCD. Based on this concept, a photometric system composed of a LED source and a digital camera as a detector was used for optical density measurements. Standards for permanganate, glucose, and protein solutions were detemined by colorimetric methods using our device. Samples of protein of Pasteurella mutocida bacteria membrane and, also, fractions of rabbit kidney membrane, rich in Na, K-ATPase, with unknown concentrations were dosed through the Hartree method using our photometric system.

Texture Image Analysis in Differentiating Malignant from Benign Adrenal Cortical Tumors in Children and Adults

Objective: To investigate the possible role of chromatin texture parameters, nuclear morphology, DNA ploidy and clinical functional status in discriminating benign from malignant adrenocortical tumors (ACT). Patients and Methods: Forty-eight cases of clinically benign (n=40) and clinically malignant (n=8) ACT with a minimum of 5-years` follow-up were evaluated for chromatin texture parameters (run length, standard deviation, configurable run length, valley, slope, peak and other 21 Markovian features that describe the distribution of the chromatin in the nucleus), nuclear morphology (nuclear area, nuclear perimeter, nuclear maximum and minumum diameter, nuclear shape), and DNA ploidy. Nuclear parameters were evaluated in Feulgen-stained 5 mu m paraffin-sections analyzed using a CAS 200 image analyzer. Results: Since ACTs present different biological features in children and adults, patients were divided into two groups: children (<= 15 years) and adults (>15 years). In the group of children DNA ploidy presented a marginal significance (p=0.05) in discriminating ACTs. None of the parameters discriminated between malignant and benign ACT in the adult group. Conclusion: ACTs are uncommon and definitive predictive criteria for malignancy remain uncertain, particularly in children. Our data point to DNA content evaluated by image analysis as a new candidate tool for this challenging task. Texture image analysis did not help to differentiate malignant from benign adrenal cortical tumors in children and adults.

Quantitative and qualitative comparison of 3.0 T and 1.5 T MR imaging of the liver in patients with diffuse parenchymal liver disease

Purpose: The purpose of our study was to compare signal characteristics and image qualities of MR imaging at 3.0 T and 1.5 T in patients with diffuse parenchymal liver disease. Materials and methods: 25 consecutive patients with diffuse parenchymal liver disease underwent abdominal MR imaging at both 3.0 T and 1.5 T within a 6-month interval. A retrospective study was conducted to obtain quantitative and qualitative data from both 3.0 T and 1.5 T MRI. Quantitative image analysis was performed by measuring the signal-to-noise ratios (SNRs) and the contrast-to-noise ratios (CNRs) by the Students t-test. Qualitative image analysis was assessed by grading each sequence on a 3- and 4-point scale, regarding the presence of artifacts and image quality, respectively. Statistical analysis consisted of the Wilcoxon signed-rank test. Results: the mean SNRs and CNRs of the liver parenchyma and the portal vein were significantly higher at 3.0 T than at 1.5 T on portal and equilibrium phases of volumetric interpolated breath-hold examination (VIBE) images (P < 0.05). The mean SNRs were significantly higher at 3.0 T than at 1.5 T on T1-weighted spoiled gradient echo (SGE) images (P < 0.05). However, there were no significantly differences on T2-weighted short-inversion-time inversion recovery (STIR) images. Overall image qualities of the 1.5 T noncontrast T1- and T2-weighted sequences were significantly better than 3.0 T (P < 0.01). In contrast, overall image quality of the 3.0 T post-gadolinium VIBE sequence was significantly better than 1.5 T (P< 0.01). Conclusions: MR imaging of post-gadolinium VIBE sequence at 3.0 T has quantitative and qualitative advantages of evaluating for diffuse parenchymal liver disease. (C) 2008 Elsevier Ireland Ltd. All rights reserved.

Effect of image analysis software on neurofunctional activation during processing of emotional human faces

Functional brain imaging techniques such as functional MRI (fMRI) that allow the in vivo investigation of the human brain have been exponentially employed to address the neurophysiological substrates of emotional processing. Despite the growing number of fMRI studies in the field, when taken separately these individual imaging studies demonstrate contrasting findings and variable pictures, and are unable to definitively characterize the neural networks underlying each specific emotional condition. Different imaging packages, as well as the statistical approaches for image processing and analysis, probably have a detrimental role by increasing the heterogeneity of findings. In particular, it is unclear to what extent the observed neurofunctional response of the brain cortex during emotional processing depends on the fMRI package used in the analysis. In this pilot study, we performed a double analysis of an fMRI dataset using emotional faces. The Statistical Parametric Mapping (SPM) version 2.6 (Wellcome Department of Cognitive Neurology, London, UK) and the XBAM 3.4 (Brain Imaging Analysis Unit, Institute of Psychiatry, Kings College London, UK) programs, which use parametric and non-parametric analysis, respectively, were used to assess our results. Both packages revealed that processing of emotional faces was associated with an increased activation in the brain`s visual areas (occipital, fusiform and lingual gyri), in the cerebellum, in the parietal cortex, in the cingulate cortex (anterior and posterior cingulate), and in the dorsolateral and ventrolateral prefrontal cortex. However, blood oxygenation level-dependent (BOLD) response in the temporal regions, insula and putamen was evident in the XBAM analysis but not in the SPM analysis. Overall, SPM and XBAM analyses revealed comparable whole-group brain responses. Further Studies are needed to explore the between-group compatibility of the different imaging packages in other cognitive and emotional processing domains. (C) 2009 Elsevier Ltd. All rights reserved.

Lower Extremities Magnetic Resonance Angiography With Blood Pressure Cuff Compression: Quantitative Dynamic Analysis

Purpose: To quantitatively evaluate changes induced by the application of a femoral blood-pressure cuff (BPC) on run-off magnetic resonance angiography (MRA). which is a method generally previously proposed to reduce venous contamination in the leg. Materials and Methods: This study was Health Insurance Portability and Accountability Act (HIPAA)- and Institutional Review Board (IRB)-compliant, We used time-resolved gradient-echo gadolinium (Gd)-enhanced MRA to measure BPC effects on arterial, venous, and soft-tissue enhancement. Seven healthy volunteers (six men) were studied with the BPC applied at the mid-femoral level unilaterally using a 1.5T MR system after intravenous injection of Gd-BOPTA. Different statistical tools were used such as the Wilcoxon signed rank test and a cubic smoothing spline fit. Results: We found that BPC application induces delayed venous filling (as previously described), but also induces significant decreases in arterial inflow, arterial enhancement, vascular-soft tissue contrast, and delayed peak enhancement (which have not been previously measured). Conclusion: The potential benefits from using a BPC for run-off MRA must be balanced against the potential pitfalls, elucidated by our findings.

Quantitative CT analysis of the glabellar and anterior nasal spine regions for the placement of implants for nasal prosthesis retention

Objectives: The aim of this study was to determine the precision of the measurements of 2 craniometric anatomic points-glabella and anterior nasal spine-in order to verify their possibility as potential locations for placing implants aimed at nasal prostheses retention. Methods: Twenty-six dry human skulls were scanned in a high-resolution spiral tomography with 1-mm axial slice thickness and 1-mm interval reconstruction using a bone tissue filter. Images obtained were stored and transferred to an independent workstation containing e-film imaging software. The measurements (in the glabella and anterior nasal fossa) were made independently by 2 observers twice for each measurement. Data were submitted to statistical analysis (parametric t test). Results: The results demonstrated no statistically significant difference between interobserver and intraobserver measurements (P > .05). The standard error was found to be between 0.49 mm and 0.84 mrn for measurements in bone protocol, indicating a high /eve/ of precision. Conclusions: The measurements obtained in anterior nasal spine and glabella were considered precise and reproducible. Mean values of such measurements pointed to the possibility of implant placement in these regions, particularly in the anterior nasal spine.