Refine
Year of publication
Document Type
- ZIB-Report (19)
- Doctoral Thesis (1)
Has Fulltext
- yes (20) (remove)
Is part of the Bibliography
- no (20)
Keywords
- model-based inversion (3)
- Monte Carlo (2)
- blood oxygen saturation (2)
- inverse problem (2)
- quantitative photoacoustic imaging (2)
- spectral unmixing (2)
- 2-photon microscopy (1)
- 2D distance map (1)
- Amira (1)
- Beton (1)
In this paper, we present a software-assisted workflow for the alignment and matching of filamentous structures across a stack of 3D serial image sections. This is achieved by a combination of automatic methods, visual validation, and interactive correction. After an initial alignment, the user can continuously improve the result by interactively correcting landmarks or matches of filaments. This is supported by a quality assessment that visualizes regions that have been already inspected and, thus, allows a trade-off between quality and manual labor.
The software tool was developed in collaboration with biologists who investigate microtubule-based spindles during cell division. To quantitatively understand the structural organization of such spindles, a 3D reconstruction of the numerous microtubules is essential. Each spindle is cut into a series of semi-thick physical sections, of which electron tomograms are acquired. The sections then need to be stitched, i.e. non-rigidly aligned; and the microtubules need to be traced in each section and connected across section boundaries. Experiments led to the conclusion that automatic methods for stitching alone provide only an incomplete solution to practical analysis needs. Automatic methods may fail due to large physical distortions, a low signal-to-noise ratio of the images, or other unexpected experimental difficulties. In such situations, semi-automatic validation and correction is required to rescue as much information as possible to derive biologically meaningful results despite of some errors related to data collection.
Since the correct stitching is visually not obvious due to the number of microtubules (up to 30k) and their dense spatial arrangement, these are difficult tasks. Furthermore, a naive inspection of each microtubule is too time consuming. In addition, interactive visualization is hampered by the size of the image data (up to 100 GB). Based on the requirements of our collaborators, we present a practical solution for the semi-automatic stitching of serial section image stacks with filamentous structures.
Quantitative photoacoustic tomography aims to recover maps of the local concentrations of tissue chromophores from multispectral images. While model-based inversion schemes are promising approaches, major challenges to their practical implementation include the unknown fluence distribution and the scale of the inverse problem. This paper describes an inversion scheme based on a radiance Monte Carlo model and an adjoint-assisted gradient optimization that incorporates fluence-dependent step sizes and adaptive moment estimation. The inversion is shown to recover absolute chromophore concentrations, blood oxygen saturation and the Grüneisen parameter from in silico 3D phantom images for different radiance approximations. The scattering coefficient was assumed to be homogeneous and known a priori.
Quantitative PA tomography of high resolution 3-D images: experimental validation in tissue phantoms
(2019)
Quantitative photoacoustic tomography aims recover the spatial distribution of absolute chromophore concentrations and their ratios from deep tissue, high-resolution images. In this study, a model-based inversion scheme based on a Monte-Carlo light transport model is experimentally validated on 3-D multispectral images of a tissue phantom acquired using an all-optical scanner with a planar detection geometry. A calibrated absorber allowed scaling of the measured data during the inversion, while an acoustic correction method was employed to compensate the effects of limited view detection. Chromophore- and fluence-dependent step sizes and Adam optimization were implemented to achieve rapid convergence. High resolution 3-D maps of absolute concentrations and their ratios were recovered with high accuracy. Potential applications of this method include quantitative functional and molecular photoacoustic tomography of deep tissue in preclinical and clinical studies.
Following axon pathfinding, growth cones transition from stochastic filopodial exploration to the formation of a limited number of synapses. How the interplay of filopodia and synapse assembly ensures robust connectivity in the brain has remained a challenging problem. Here, we developed a new 4D analysis method for filopodial dynamics and a data-driven computational model of synapse formation for R7 photoreceptor axons in developing Drosophila brains. Our live data support a 'serial synapse formation' model, where at any time point only a single 'synaptogenic' filopodium suppresses the synaptic competence of other filopodia through competition for synaptic seeding factors. Loss of the synaptic seeding factors Syd-1 and Liprin-α leads to a loss of this suppression, filopodial destabilization and reduced synapse formation, which is sufficient to cause the destabilization of entire axon terminals. Our model provides a filopodial 'winner-takes-all' mechanism that ensures the formation of an appropriate number of synapses.
Mitotic and meiotic spindles are microtubule-based structures to faithfully segregate chromosomes. Electron tomography is currently the method of choice to analyze the three-dimensional architecture of both types of spindles. Over the years, we have developed methods and software for automatic segmentation and stitching of microtubules in serial sections for large-scale reconstructions. Three-dimensional reconstruction of microtubules, however, is only the first step towards biological insight. The second step is the analysis of the structural data to derive measurable spindle properties. Here, we present a comprehensive set of techniques to quantify spindle parameters. These techniques provide quantitative analyses of specific microtubule classes and are applicable to a variety of tomographic reconstructions of spindles from different organisms.
In oocytes of many organisms, meiotic spindles form in the absence of centrosomes [1–5]. Such female meiotic spindles have a pointed appearance in metaphase with microtubules focused at acentrosomal spindle poles. At anaphase, the microtubules of acentrosomal spindles then transition to an inter- chromosomal array, while the spindle poles disappear. This transition is currently not understood. Previous studies have focused on this inter- chromosomal microtubule array and proposed a pushing model to drive chromosome segregation [6, 7]. This model includes an end-on orientation of microtubules with chromosomes. Alternatively, chromosomes were thought to associate along bundles of microtubules [8, 9]. Starting with metaphase, this second model proposed a pure lateral chromosome-to-microtubule association up to the final meiotic stages of anaphase. Here we applied large-scale electron tomography [10] of staged C. elegans oocytes in meiosis to analyze the orientation of microtubules in respect to chromosomes. We show that microtubules at metaphase I are primarily oriented laterally to the chromosomes and that microtubules switch to an end-on orientation during progression through anaphase. We further show that this switch in microtubule orientation involves a kinesin-13 microtubule depolymerase, KLP-7, which removes laterally associated microtubules around chromosomes. From this we conclude that both lateral and end-on modes of microtubule-to-chromosome orientations are successively used in C. elegans oocytes to segregate meiotic chromosomes.
Introduction – Many biological structures show recurring tiling patterns on one structural level or the other. Current image acquisition techniques are able to resolve those tiling patterns to allow quantitative analyses. The resulting image data, however, may contain an enormous number of elements. This renders manual image analysis infeasible, in particular when statistical analysis is to be conducted, requiring a larger number of image data to be analyzed. As a consequence, the analysis process needs to be automated to a large degree. In this paper, we describe a multi-step image segmentation pipeline for the automated segmentation of the calcified cartilage into individual tesserae from computed tomography images of skeletal elements of stingrays.
Methods – Besides applying state-of-the-art algorithms like anisotropic diffusion smoothing, local thresholding for foreground segmentation, distance map calculation, and hierarchical watershed, we exploit a graph-based representation for fast correction of the segmentation. In addition, we propose a new distance map that is computed only in the plane that locally best approximates the calcified cartilage. This distance map drastically improves the separation of individual tesserae. We apply our segmentation pipeline to hyomandibulae from three individuals of the round stingray (Urobatis halleri), varying both in age and size.
Results – Each of the hyomandibula datasets contains approximately 3000 tesserae. To evaluate the quality of the automated segmentation, four expert users manually generated ground truth segmentations of small parts of one hyomandibula. These ground truth segmentations allowed us to compare the segmentation quality w.r.t. individual tesserae. Additionally, to investigate the segmentation quality of whole skeletal elements, landmarks were manually placed on all tesserae and their positions were then compared to the segmented tesserae. With the proposed segmentation pipeline, we sped up the processing of a single skeletal element from days or weeks to a few hours.
Adapting trabecular structures for 3D printing: an image processing approach based on µCT data
(2017)
Materials with a trabecular structure notably combine advantages such as lightweight, reasonable strength, and permeability for fluids. This combination of advantages is especially interesting for tissue engineering in trauma surgery and orthopedics. Bone-substituting scaffolds for instance are designed with a trabecular structure in order to allow cell migration for bone ingrowth and vascularization. An emerging and recently very popular technology to produce such complex, porous structures is 3D printing. However, several technological aspects regarding the scaffold architecture, the printable resolution, and the feature size have to be considered when fabricating scaffolds for bone tissue replacement and regeneration.
Here, we present a strategy to assess and prepare realistic trabecular structures for 3D printing using image analysis with the aim of preserving the structural elements. We discuss critical conditions of the printing system and present a 3-stage approach to adapt a trabecular structure from $\mu$CT data while incorporating knowledge about the printing system. In the first stage, an image-based extraction of solid and void structures is performed, which results in voxel- and graph-based representations of the extracted structures. These representations not only allow us to quantify geometrical properties such as pore size or strut geometry and length. But, since the graph represents the geometry and the topology of the initial structure, it can be used in the second stage to modify and adjust feature size, volume and sample size in an easy and consistent way. In the final reconstruction stage, the graph is then converted into a voxel representation preserving the topology of the initial structure. This stage generates a model with respect to the printing conditions to ensure a stable and controlled voxel placement during the printing process.
The role of titanium surface nanotopography on preosteoblast morphology, adhesion and migration
(2017)
Surface structuring of titanium-based implants with appropriate nanotopographies can significantly modulate their impact on the biological behavior of cells populating these implants. Implant assisted bone tissue repair and regeneration require functional adhesion and expansion of bone progenitors. The surface nanotopography of implant materials used to support bone healing and its effect on cell behavior, in particular cell adhesion, spreading, expansion, and motility, is still not clearly understood. The aim of this study is to investigate preosteoblast proliferation, adhesion, morphology, and migration on different titanium materials with similar surface chemistry, but distinct nanotopographical features. Sonochemical treatment and anodic oxidation were employed to fabricate disordered – mesoporous titania (TMS), and ordered – titania nanotubular (TNT) topographies respectively. The morphological evaluation revealed a surface dependent shape, thickness, and spreading of cells owing to different adherence behavior. Cells were polygonal-shaped and well-spread on glass and TMS, but displayed an elongated fibroblast-like morphology on TNT surfaces. The cells on glass however, were much flatter than on nanostructured surfaces. Both nanostructured surfaces impaired cell adhesion, but TMS was more favorable for cell growth due to its support of cell attachment and spreading in contrast to TNT. Quantitative wound healing assay in combination with live-cell imaging revealed that cells seeded on TMS surfaces migrated in close proximity to neighboring cells and less directed when compared to the migratory behavior on other surfaces. The results indicate distinctly different cell adhesion and migration on ordered and disordered titania nanotopographies, providing important information that could be used in optimizing titanium-based scaffold design to foster bone tissue growth and repair.
The goal of quantitative photoacoustic tomography (qPAT) is to recover maps of the chromophore distributions from multiwavelength images of the initial pressure. Model-based inversions that incorporate the physical processes underlying the photoacoustic (PA) signal generation represent a promising approach. Monte-Carlo models of the light transport are computationally expensive, but provide accurate fluence distributions predictions, especially in the ballistic and quasi-ballistic regimes. Here, we focus on the inverse problem of 3D qPAT of blood oxygenation and investigate the application of the Monte-Carlo method in a model-based inversion scheme. A forward model of the light transport based on the MCX simulator and acoustic propagation modeled by the k-Wave toolbox was used to generate a PA image data set acquired in a tissue phantom over a planar detection geometry. The combination of the optical and acoustic models is shown to account for limited-view artifacts. In addition, the errors in the fluence due to, for example, partial volume artifacts and absorbers immediately adjacent to the region of interest are investigated. To accomplish large-scale inversions in 3D, the number of degrees of freedom is reduced by applying image segmentation to the initial pressure distribution to extract a limited number of regions with homogeneous optical parameters. The absorber concentration in the tissue phantom was estimated using a coordinate descent parameter search based on the comparison between measured and modeled PA spectra. The estimated relative concentrations using this approach lie within 5 % compared to the known concentrations. Finally, we discuss the feasibility of this approach to recover the blood oxygenation from experimental data.