Long-Term Follow-Up after Laser-Assisted Pulmonary Metastasectomy Shows Complete Lung Function Recovery

Preserving maximum lung function is a fundamental goal of parenchymal-sparing pulmonary laser surgery. Long-term studies for follow-up of lung function after pulmonary laser metastasectomy are lacking. However, a sufficient postoperative lung function is essential for quality of life and reduces potential postoperative complications. In this study, we investigate the extent of loss in lung function following pulmonary laser resection after three, six, and twelve months. We conducted a retrospective analysis using a prospective database of 4595 patients, focusing on 126 patients who underwent unilateral pulmonary laser resection for lung metastases from 1996 to 2022 using a 1318 nm Nd:YAG laser or a high-power pure diode laser. Results show that from these patients, a median of three pulmonary nodules were removed, with 75% presenting central lung lesions and 25% peripheral lesions. The median preoperative FEV1 was 98% of the predicted value, decreasing to 71% postoperatively but improving to 90% after three months, 93% after six months, and 96% after twelve months. Statistical analysis using the Friedman test indicated no significant difference in FEV1 between preoperative levels and those at six and twelve months post-surgery. The findings confirm that pulmonary laser surgery effectively preserves lung function over time, with patients generally regaining their preoperative lung function within a year, regardless of the metastases' location.

A Local Iterative Approach for the Extraction of 2D Manifolds from Strongly Curved and Folded Thin-Layer Structures

Ridge surfaces represent important features for the analysis of 3-dimensional (3D) datasets in diverse applications and are often derived from varying underlying data including flow fields, geological fault data, and point data, but they can also be present in the original scalar images acquired using a plethora of imaging techniques. Our work is motivated by the analysis of image data acquired using micro-computed tomography ([Formula: see text]) of ancient, rolled and folded thin-layer structures such as papyrus, parchment, and paper as well as silver and lead sheets. From these documents we know that they are 2-dimensional (2D) in nature. Hence, we are particularly interested in reconstructing 2D manifolds that approximate the document's structure. The image data from which we want to reconstruct the 2D manifolds are often very noisy and represent folded, densely-layered structures with many artifacts, such as ruptures or layer splitting and merging. Previous ridge-surface extraction methods fail to extract the desired 2D manifold for such challenging data. We have therefore developed a novel method to extract 2D manifolds. The proposed method uses a local fast marching scheme in combination with a separation of the region covered by fast marching into two sub-regions. The 2D manifold of interest is then extracted as the surface separating the two sub-regions. The local scheme can be applied for both automatic propagation as well as interactive analysis. We demonstrate the applicability and robustness of our method on both artificial data as well as real-world data including folded silver and papyrus sheets.

Dense reconstruction of elephant trunk musculature

The elephant trunk operates as a muscular hydrostat1,2 and is actuated by the most complex musculature known in animals.3,4 Because the number of trunk muscles is unclear,5 we performed dense reconstructions of trunk muscle fascicles, elementary muscle units, from microCT scans of an Asian baby elephant trunk. Muscle architecture changes markedly across the trunk. Trunk tip and finger consist of about 8,000 extraordinarily filigree fascicles. The dexterous finger consists exclusively of microscopic radial fascicles pointing to a role of muscle miniaturization in elephant dexterity. Radial fascicles also predominate (at 82% volume) the remainder of the trunk tip, and we wonder if radial muscle fascicles are of particular significance for fine motor control of the dexterous trunk tip. By volume, trunk-shaft muscles6 comprise one-third of the numerous, small radial muscle fascicles; two-thirds of the three subtypes of large longitudinal fascicles (dorsal longitudinals, ventral outer obliques, and ventral inner obliques);7,8,9 and a small fraction of transversal fascicles. Shaft musculature is laterally, but not radially, symmetric. A predominance of dorsal over ventral radial muscles and of ventral over dorsal longitudinal muscles may result in a larger ability of the shaft to extend dorsally than ventrally10 and to bend inward rather than outward. There are around 90,000 trunk muscle fascicles. While primate hand control is based on fine control of contraction by the convergence of many motor neurons on a small set of relatively large muscles, evolution of elephant grasping has led to thousands of microscopic fascicles, which probably outnumber facial motor neurons.

Serial-section Electron Tomography and Quantitative Analysis of Microtubule Organization in 3D-reconstructed Mitotic Spindles

For the analysis of cellular architecture during mitosis, nanometer resolution is needed to visualize the organization of microtubules in spindles. Here, we present a detailed protocol that can be used to produce 3D reconstructions of whole mitotic spindles in cells grown in culture. For this, we attach mammalian cells enriched in mitotic stages to sapphire discs. Our protocol further involves cryo-immobilization by high-pressure freezing, freeze-substitution, and resin embedding. We then use fluorescence light microscopy to stage select mitotic cells in the resin-embedded samples. This is followed by large-scale electron tomography to reconstruct the selected and staged mitotic spindles in 3D. The generated and stitched electron tomograms are then used to semi-automatically segment the microtubules for subsequent quantitative analysis of spindle organization. Thus, by providing a detailed correlative light and electron microscopy (CLEM) approach, we give cell biologists a toolset to streamline the 3D visualization and analysis of spindle microtubules (http://kiewisz.shinyapps.io/asga). In addition, we refer to a recently launched platform that allows for an interactive display of the 3D-reconstructed mitotic spindles (https://cfci.shinyapps.io/ASGA_3DViewer/). Key features • High-throughput screening of mitotic cells by correlative light and electron microscopy (CLEM). • Serial-section electron tomography of selected cells. • Visualization of mitotic spindles in 3D and quantitative analysis of microtubule organization.

Simulation-based inference for efficient identification of generative models in computational connectomics

Recent advances in connectomics research enable the acquisition of increasing amounts of data about the connectivity patterns of neurons. How can we use this wealth of data to efficiently derive and test hypotheses about the principles underlying these patterns? A common approach is to simulate neuronal networks using a hypothesized wiring rule in a generative model and to compare the resulting synthetic data with empirical data. However, most wiring rules have at least some free parameters, and identifying parameters that reproduce empirical data can be challenging as it often requires manual parameter tuning. Here, we propose to use simulation-based Bayesian inference (SBI) to address this challenge. Rather than optimizing a fixed wiring rule to fit the empirical data, SBI considers many parametrizations of a rule and performs Bayesian inference to identify the parameters that are compatible with the data. It uses simulated data from multiple candidate wiring rule parameters and relies on machine learning methods to estimate a probability distribution (the 'posterior distribution over parameters conditioned on the data') that characterizes all data-compatible parameters. We demonstrate how to apply SBI in computational connectomics by inferring the parameters of wiring rules in an in silico model of the rat barrel cortex, given in vivo connectivity measurements. SBI identifies a wide range of wiring rule parameters that reproduce the measurements. We show how access to the posterior distribution over all data-compatible parameters allows us to analyze their relationship, revealing biologically plausible parameter interactions and enabling experimentally testable predictions. We further show how SBI can be applied to wiring rules at different spatial scales to quantitatively rule out invalid wiring hypotheses. Our approach is applicable to a wide range of generative models used in connectomics, providing a quantitative and efficient way to constrain model parameters with empirical connectivity data.

Evolutionary traces of miniaturization in a giant-Comparative anatomy of brain and brain nerves in Bathochordaeus stygius (Tunicata, Appendicularia)

Appendicularia comprises 70 marine, invertebrate, chordate species. Appendicularians play important ecological and evolutionary roles, yet their morphological disparity remains understudied. Most appendicularians are small, develop rapidly, and with a stereotyped cell lineage, leading to the hypothesis that Appendicularia derived progenetically from an ascidian-like ancestor. Here, we describe the detailed anatomy of the central nervous system of Bathochordaeus stygius, a giant appendicularian from the mesopelagic. We show that the brain consists of a forebrain with on average smaller and more uniform cells and a hindbrain, in which cell shapes and sizes vary to a greater extent. Cell count for the brain was 102. We demonstrate the presence of three paired brain nerves. Brain nerve 1 traces into the epidermis of the upper lip region and consists of several fibers with some supportive bulb cells in its course. Brain nerve 2 innervates oral sensory organs and brain nerve 3 innervates the ciliary ring of the gill slits and lateral epidermis. Brain nerve 3 is asymmetric, with the right nerve consisting of two neurites originating posterior to the left one that contains three neurites. Similarities and differences to the anatomy of the brain of the model species Oikopleura dioica are discussed. We interpret the small number of cells in the brain of B. stygius as an evolutionary trace of miniaturization and conclude that giant appendicularians evolved from a small, progenetic ancestor that secondarily increased in size within Appendicularia.

High-throughput segmentation, data visualization, and analysis of sea star skeletal networks

The remarkably complex skeletal systems of the sea stars (Echinodermata, Asteroidea), consisting of hundreds to thousands of individual elements (ossicles), have intrigued investigators for more than 150 years. While the general features and structural diversity of isolated asteroid ossicles have been well documented in the literature, the task of mapping the spatial organization of these constituent skeletal elements in a whole-animal context represents an incredibly laborious process, and as such, has remained largely unexplored. To address this unmet need, particularly in the context of understanding structure-function relationships in these complex skeletal systems, we present an integrated approach that combines micro-computed tomography, semi-automated ossicle segmentation, data visualization tools, and the production of additively manufactured tangible models to reveal biologically relevant structural data that can be rapidly analyzed in an intuitive manner. In the present study, we demonstrate this high-throughput workflow by segmenting and analyzing entire skeletal systems of the giant knobby star, Pisaster giganteus, at four different stages of growth. The in-depth analysis, presented herein, provides a fundamental understanding of the three-dimensional skeletal architecture of the sea star body wall, the process of skeletal maturation during growth, and the relationship between skeletal organization and morphological characteristics of individual ossicles. The widespread implementation of this approach for investigating other species, subspecies, and growth series has the potential to fundamentally improve our understanding of asteroid skeletal architecture and biodiversity in relation to mobility, feeding habits, and environmental specialization in this fascinating group of echinoderms.

A Visual Interface for Exploring Hypotheses about Neural Circuits

One of the fundamental problems in neurobiological research is to understand how neural circuits generate behaviors in response to sensory stimuli. Elucidating such neural circuits requires anatomical and functional information about the neurons that are active during the processing of the sensory information and generation of the respective response, as well as an identification of the connections between these neurons. With modern imaging techniques, both morphological properties of individual neurons as well as functional information related to sensory processing, information integration and behavior can be obtained. Given the resulting information, neurobiologists are faced with the task of identifying the anatomical structures down to individual neurons that are linked to the studied behavior and the processing of the respective sensory stimuli. Here, we present a novel interactive tool that assists neurobiologists in the aforementioned task by allowing them to extract hypothetical neural circuits constrained by anatomical and functional data. Our approach is based on two types of structural data: brain regions that are anatomically or functionally defined, and morphologies of individual neurons. Both types of structural data are interlinked and augmented with additional information. The presented tool allows the expert user to identify neurons using Boolean queries. The interactive formulation of these queries is supported by linked views, using, among other things, two novel 2D abstractions of neural circuits. The approach was validated in two case studies investigating the neural basis of vision-based behavioral responses in zebrafish larvae. Despite this particular application, we believe that the presented tool will be of general interest for exploring hypotheses about neural circuits in other species, genera and taxa.

MCRS1 modulates the heterogeneity of microtubule minus-end morphologies in mitotic spindles

Faithful chromosome segregation requires the assembly of a bipolar spindle, consisting of two antiparallel microtubule (MT) arrays having most of their minus ends focused at the spindle poles and their plus ends overlapping in the spindle midzone. Spindle assembly, chromosome alignment, and segregation require highly dynamic MTs. The plus ends of MTs have been extensively investigated but their minus-end structure remains poorly characterized. Here, we used large-scale electron tomography to study the morphology of the MT minus ends in three dimensionally reconstructed metaphase spindles in HeLa cells. In contrast to the homogeneous open morphology of the MT plus ends at the kinetochores, we found that MT minus ends are heterogeneous, showing either open or closed morphologies. Silencing the minus end-specific stabilizer, MCRS1 increased the proportion of open MT minus ends. Altogether, these data suggest a correlation between the morphology and the dynamic state of the MT ends. Taking this heterogeneity of the MT minus-end morphologies into account, our work indicates an unsynchronized behavior of MTs at the spindle poles, thus laying the groundwork for further studies on the complexity of MT dynamics regulation.

Volumetric macromolecule identification in cryo-electron tomograms using capsule networks

Background

Despite recent advances in cellular cryo-electron tomography (CET), developing automated tools for macromolecule identification in submolecular resolution remains challenging due to the lack of annotated data and high structural complexities. To date, the extent of the deep learning methods constructed for this problem is limited to conventional Convolutional Neural Networks (CNNs). Identifying macromolecules of different types and sizes is a tedious and time-consuming task. In this paper, we employ a capsule-based architecture to automate the task of macromolecule identification, that we refer to as 3D-UCaps. In particular, the architecture is composed of three components: feature extractor, capsule encoder, and CNN decoder. The feature extractor converts voxel intensities of input sub-tomograms to activities of local features. The encoder is a 3D Capsule Network (CapsNet) that takes local features to generate a low-dimensional representation of the input. Then, a 3D CNN decoder reconstructs the sub-tomograms from the given representation by upsampling.

Results

We performed binary and multi-class localization and identification tasks on synthetic and experimental data. We observed that the 3D-UNet and the 3D-UCaps had an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1-$$\end{document}F1-score mostly above 60% and 70%, respectively, on the test data. In both network architectures, we observed degradation of at least 40% in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1$$\end{document}F1-score when identifying very small particles (PDB entry 3GL1) compared to a large particle (PDB entry 4D8Q). In the multi-class identification task of experimental data, 3D-UCaps had an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1$$\end{document}F1-score of 91% on the test data in contrast to 64% of the 3D-UNet. The better \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1$$\end{document}F1-score of 3D-UCaps compared to 3D-UNet is obtained by a higher precision score. We speculate this to be due to the capsule network employed in the encoder. To study the effect of the CapsNet-based encoder architecture further, we performed an ablation study and perceived that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_1$$\end{document}F1-score is boosted as network depth is increased which is in contrast to the previously reported results for the 3D-UNet. To present a reproducible work, source code, trained models, data as well as visualization results are made publicly available.

Conclusion

Quantitative and qualitative results show that 3D-UCaps successfully perform various downstream tasks including identification and localization of macromolecules and can at least compete with CNN architectures for this task. Given that the capsule layers extract both the existence probability and the orientation of the molecules, this architecture has the potential to lead to representations of the data that are better interpretable than those of 3D-UNet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04901-w.

Three-dimensional structure of kinetochore-fibers in human mitotic spindles

During cell division, kinetochore microtubules (KMTs) provide a physical linkage between the chromosomes and the rest of the spindle. KMTs in mammalian cells are organized into bundles, so-called kinetochore-fibers (k-fibers), but the ultrastructure of these fibers is currently not well characterized. Here, we show by large-scale electron tomography that each k-fiber in HeLa cells in metaphase is composed of approximately nine KMTs, only half of which reach the spindle pole. Our comprehensive reconstructions allowed us to analyze the three-dimensional (3D) morphology of k-fibers and their surrounding MTs in detail. We found that k-fibers exhibit remarkable variation in circumference and KMT density along their length, with the pole-proximal side showing a broadening. Extending our structural analysis then to other MTs in the spindle, we further observed that the association of KMTs with non-KMTs predominantly occurs in the spindle pole regions. Our 3D reconstructions have implications for KMT growth and k-fiber self-organization models as covered in a parallel publication applying complementary live-cell imaging in combination with biophysical modeling (Conway et al., 2022). Finally, we also introduce a new visualization tool allowing an interactive display of our 3D spindle data that will serve as a resource for further structural studies on mitosis in human cells.

Hide and seek shark teeth in Random Forests: machine learning applied to Scyliorhinus canicula populations

Shark populations that are distributed alongside a latitudinal gradient often display body size differences at sexual maturity and vicariance patterns related to their number of tooth files. Previous works have demonstrated that Scyliorhinus canicula populations differ between the northeastern Atlantic Ocean and the Mediterranean Sea based on biological features and genetic analysis. In this study, we sample more than 3,000 teeth from 56 S. canicula specimens caught incidentally off Roscoff and Banyuls-sur-Mer. We investigate population differences based on tooth shape and form by using two approaches. Classification results show that the classical geometric morphometric framework is outperformed by an original Random Forests-based framework. Visually, both S. canicula populations share similar ontogenetic trends and timing of gynandric heterodonty emergence but the Atlantic population has bigger, blunter teeth, and less numerous accessory cusps than the Mediterranean population. According to the models, the populations are best differentiated based on their lateral tooth edges, which bear accessory cusps, and the tooth centroid sizes significantly improve classification performances. The differences observed are discussed in light of dietary and behavioural habits of the populations considered. The method proposed in this study could be further adapted to complement DNA analyses to identify shark species or populations based on tooth morphologies. This process would be of particular interest for fisheries management and identification of shark fossils.

Ontogeny of a tessellated surface: Carapace growth of the longhorn cowfish Lactoria cornuta

Biological armors derive their mechanical integrity in part from their geometric architectures, often involving tessellations: individual structural elements tiled together to form surface shells. The carapace of boxfish, for example, is composed of mineralized polygonal plates, called scutes, arranged in a complex geometric pattern and nearly completely encasing the body. In contrast to artificial armors, the boxfish exoskeleton grows with the fish; the relationship between the tessellation and the gross structure of the armor is therefore critical to sustained protection throughout growth. To clarify whether or how the boxfish tessellation is maintained or altered with age, we quantify architectural aspects of the tessellated carapace of the longhorn cowfish Lactoria cornuta through ontogeny (across nearly an order of magnitude in standard length) and in a high‐throughput fashion, using high‐resolution microCT data and segmentation algorithms to characterize the hundreds of scutes that cover each individual. We show that carapace growth is canalized with little variability across individuals: rather than continually adding scutes to enlarge the carapace surface, the number of scutes is surprisingly constant, with scutes increasing in volume, thickness, and especially width with age. As cowfish and their scutes grow, scutes become comparatively thinner, with the scutes at the edges (weak points in a boxy architecture) being some of the thickest and most reinforced in younger animals and thinning most slowly across ontogeny. In contrast, smaller scutes with more variable curvature were found in the limited areas of more complex topology (e.g., around fin insertions, mouth, and anus). Measurements of Gaussian and mean curvature illustrate that cowfish are essentially tessellated boxes throughout life: predominantly zero curvature surfaces comprised of mostly flat scutes, and with scutes with sharp bends used sparingly to form box edges. Since growth of a curved, tiled surface with a fixed number of tiles would require tile restructuring to accommodate the surface's changing radius of curvature, our results therefore illustrate a previously unappreciated advantage of the odd boxfish morphology: by having predominantly flat surfaces, it is the box‐like body form that in fact permits a relatively straightforward growth system of this tessellated architecture (i.e., where material is added to scute edges). Our characterization of the ontogeny and maintenance of the carapace tessellation provides insights into the potentially conflicting mechanical, geometric, and developmental constraints of this species but also perspectives into natural strategies for constructing mutable tiled architectures.

The carapace of boxfish is composed of mineralized polygonal plates, called scutes, arranged in a complex geometric pattern and nearly completely encasing the body. To clarify whether or how this armor is maintained or altered with age, we quantify architectural aspects of the carapace of the longhorn cowfish Lactoria cornuta through ontogeny, using high‐resolution microCT data and segmentation algorithms to characterize the hundreds of scutes that cover each individual.

Large abdominal mechanoreceptive sense organs in small plant-dwelling insects

The Hemiptera, with approximately 98 000 species, is one of the largest insect orders. Most species feed by sucking sap from plant tissues and are thus often vectors for economically important phytopathogens. Well known within this group are the large cicadas (Cicadomorpha: Cicadoidea: Cicadidae) because they produce extremely loud airborne sounds. Less well known are their mostly tiny relatives, the leafhoppers, spittlebugs, treehoppers and planthoppers that communicate by silent vibrational signals. While the generation of these signals has been extensively investigated, the mechanisms of their perception are poorly understood. This study provides a complete description and three-dimensional reconstruction of a large and complex array of mechanoreceptors in the first abdominal segments of the Rhododendron leafhopper Graphocephala fennahi (Cicadomorpha: Membracoidea: Cicadellidae). Further, we identify homologous organs in the spittlebug Philaenus spumarius (Cicadomorpha: Cercopoidea: Aphrophoridae) and the planthopper Issus coleoptratus (Fulgoromorpha: Fulgoroidea: Issidae). Such large abdominal sensory arrays have not been found in any other insect orders studied so far. This indicates that these sense organs, together with the signal-producing tymbal organ, constitute a synapomorphy of the Tymbalia (Hemiptera excl. Sternorrhyncha). Our results contribute to the understanding of the evolution from substrate-borne to airborne communication in insects.

Semi-automatic stitching of filamentous structures in image stacks from serial-section electron tomography

We present a software-assisted workflow for the alignment and matching of filamentous structures across a three-dimensional (3D) stack of serial images. This is achieved by combining automatic methods, visual validation, and interactive correction. After the computation of an initial automatic matching, the user can continuously improve the result by interactively correcting landmarks or matches of filaments. Supported by a visual quality assessment of regions that have been already inspected, this allows a trade-off between quality and manual labour. The software tool was developed in an interdisciplinary collaboration between computer scientists and cell biologists to investigate cell division by quantitative 3D analysis of microtubules (MTs) in both mitotic and meiotic spindles. For this, each spindle is cut into a series of semi-thick physical sections, of which electron tomograms are acquired. The serial tomograms are then stitched and non-rigidly aligned to allow tracing and connecting of MTs across tomogram boundaries. In practice, automatic stitching alone provides only an incomplete solution, because large physical distortions and a low signal-to-noise ratio often cause experimental difficulties. To derive 3D models of spindles despite dealing with imperfect data related to sample preparation and subsequent data collection, semi-automatic validation and correction is required to remove stitching mistakes. However, due to the large number of MTs in spindles (up to 30k) and their resulting dense spatial arrangement, a naive inspection of each MT is too time-consuming. Furthermore, an interactive visualisation of the full image stack is hampered by the size of the data (up to 100 GB). Here, we present a specialised, interactive, semi-automatic solution that considers all requirements for large-scale stitching of filamentous structures in serial-section image stacks. To the best of our knowledge, it is the only currently available tool which is able to process data of the type and size presented here. The key to our solution is a careful design of the visualisation and interaction tools for each processing step to guarantee real-time response, and an optimised workflow that efficiently guides the user through datasets. The final solution presented here is the result of an iterative process with tight feedback loops between the involved computer scientists and cell biologists. LAY DESCRIPTION: Electron tomography of biological samples is used for a three-dimensional (3D) reconstruction of filamentous structures, such as microtubules (MTs) in mitotic and meiotic spindles. Large-scale electron tomography can be applied to increase the reconstructed volume for the visualisation of full spindles. For this, each spindle is cut into a series of semi-thick physical sections, from which electron tomograms are acquired. The serial tomograms are then stitched and non-rigidly aligned to allow tracing and connecting of MTs across tomogram boundaries. Previously, we presented fully automatic approaches for this 3D reconstruction pipeline. However, large volumes often suffer from imperfections (ie physical distortions) caused by the image acquisition process, making it difficult to apply fully automatic approaches for matching and stitching of numerous tomograms. Therefore, we developed an interactive, semi-automatic solution that considers all requirements for large-scale stitching of microtubules in image stacks of consecutive sections. We achieved this by combining automatic methods, visual validation and interactive error correction, thus allowing the user to continuously improve the result by interactively correcting landmarks or matches of filaments. We present large-scale reconstructions of spindles in which the automatic workflow failed and where different steps of manual corrections were needed. Our approach is also applicable to other biological samples showing 3D distributions of MTs in a number of different cellular contexts.

Automatic detection of prosodic boundaries in spontaneous speech

Automatic speech recognition (ASR) and natural language processing (NLP) are expected to benefit from an effective, simple, and reliable method to automatically parse conversational speech. The ability to parse conversational speech depends crucially on the ability to identify boundaries between prosodic phrases. This is done naturally by the human ear, yet has proved surprisingly difficult to achieve reliably and simply in an automatic manner. Efforts to date have focused on detecting phrase boundaries using a variety of linguistic and acoustic cues. We propose a method which does not require model training and utilizes two prosodic cues that are based on ASR output. Boundaries are identified using discontinuities in speech rate (pre-boundary lengthening and phrase-initial acceleration) and silent pauses. The resulting phrases preserve syntactic validity, exhibit pitch reset, and compare well with manual tagging of prosodic boundaries. Collectively, our findings support the notion of prosodic phrases that represent coherent patterns across textual and acoustic parameters.

D’Arcy W. Thompson’s Cartesian transformations: a critical evaluation

The images of D’Arcy Wentworth Thompson’s book “On Growth and Form” got an iconic status and became influential for biometrics and other mathematical approaches to organismic form. In particular, this is true for those of the chapter on the theory of transformation, which even has an impact on art and humanities. Based on his approach, Thompson formulated far-reaching conclusions with a partly anti-Darwinian stance. Here, we use the example of Thompson’s transformation of crab carapaces to test to what degree the transformation of grids, landmarks, and shapes result in congruent images. For comparison, we applied the same series of tests to digitized carapaces of real crabs. Both approaches show similar results. Only the simple transformations show a reasonable form of congruence. In particular, the transformations to majoid spider crabs reveal a complicated transformation of grids with partly crossing lines. By contrast, the carapace of the lithodid species is relatively easily created despite the fact that it is no brachyuran, but evolved a spider crab-like shape convergently from a hermit crab ancestor.

Image analysis pipeline for segmentation of a biological porosity network, the lacuno-canalicular system in stingray tesserae

A prerequisite for many analysis tasks in modern comparative biology is the segmentation of 3-dimensional (3D) images of the specimens being investigated (e.g. from microCT data). Depending on the specific imaging technique that was used to acquire the images and on the image resolution, different segmentation tools are required. While some standard tools exist that can often be applied for specific subtasks, building whole processing pipelines solely from standard tools is often difficult. Some tasks may even necessitate the implementation of manual interaction tools to achieve a quality that is sufficient for subsequent analysis. In this work, we present a pipeline of segmentation tools that can be used for the semiautomatic segmentation and quantitative analysis of voids in tissue (i.e. internal structural porosity). We use this pipeline to analyze lacuno-canalicular networks in stingray tesserae from 3D images acquired with synchrotron microCT.•

The first step of this pipeline, the segmentation of the tesserae, was performed using standard marker-based watershed segmentation.

•

The efficient processing of the next two steps, that is, the segmentation of all lacunae spaces belonging to a specific tessera and the separation of these spaces into individual lacunae required recently developed, novel tools.

•

For error correction, we developed an interactive method that allowed us to quickly split lacunae that were accidentally merged, and to merge lacunae that were wrongly split.

•

Finally, the tesserae and their corresponding lacunae were subdivided into structural wedges (i.e. specific anatomical regions) using a semi-manual approach.

With this processing pipeline, analysis of a variety of interconnected structural networks (e.g. vascular or lacuno-canalicular networks) can be achieved in a comparatively high-throughput fashion. In our study system, we were able to efficiently segment more than 12,000 lacunae in high-resolution scans of nine tesserae, providing a robust data set for statistical analysis.

Co-aligned chondrocytes: Zonal morphological variation and structured arrangement of cell lacunae in tessellated cartilage

In most vertebrates the embryonic cartilaginous skeleton is replaced by bone during development. During this process, cartilage cells (chondrocytes) mineralize the extracellular matrix and undergo apoptosis, giving way to bone cells (osteocytes). In contrast, sharks and rays (elasmobranchs) have cartilaginous skeletons throughout life, where only the surface mineralizes, forming a layer of tiles (tesserae). Elasmobranch chondrocytes, unlike those of other vertebrates, survive cartilage mineralization and are maintained alive in spaces (lacunae) within tesserae. However, the functions of the chondrocytes in the mineralized tissue remain unknown. Applying a custom analysis workflow to high-resolution synchrotron microCT scans of tesserae, we characterize the morphologies and arrangements of stingray chondrocyte lacunae, using lacunar morphology as a proxy for chondrocyte morphology. We show that the cell density is comparable in unmineralized and mineralized tissue and that cells maintain similar volume even when they have been incorporated into tesserae. Our findings support previous hypotheses that elasmobranch chondrocytes, unlike those of other taxa, do not proliferate, hypertrophy or undergo apoptosis during mineralization. Tessera lacunae show zonal variation in their shapes, being flatter further from and more spherical closer to the unmineralized cartilage matrix, and larger in the center of tesserae. The lacunae show pronounced organization into parallel layers and strong orientation toward neighboring tesserae. Tesserae also exhibit local variation in lacunar density, with the density considerably higher near pores passing through the tesseral layer, suggesting pores and cells interact, and that pores may contain a nutrient source. We propose that the different lacunar types reflect the stages of the tesserae formation process, while also representing local variation in tissue architecture and cell function. Lacunae are linked by small passages (canaliculi) in the matrix to form elongated series at the tesseral periphery and tight clusters in the center of tesserae, creating a rich connectivity among cells. The network arrangement and the shape variation of chondrocytes in tesserae indicate that cells may interact within and between tesserae and manage mineralization differently from chondrocytes in other vertebrates, perhaps performing analogous roles to osteocytes in bone.

Exposure to Odors Increases Pain Threshold in Chronic Low Back Pain Patients

Objectives

Structured exposure to odors is an acknowledged therapy in patients with smell loss but has also been shown to be effective in depression. The latter might rely on connections between olfactory and emotional structures, suggesting possible effects of a similar approach in pain patients. Based on neuroanatomy, there are several interfaces between the “pain matrix” and olfactory system, such as the limbic system, hypothalamus, and mediodorsal thalamus. We aimed to investigate whether structured exposure to odors may impact perceived pain in patients with chronic low back pain.

Design

Randomized controlled parallel-group design. Subjects were tested on two occasions, at baseline and after four weeks.

Setting

Ambulatory.

Subjects

Forty-two patients with chronic low back pain

Methods

For all patients, olfactory function (using the “Sniffin’Sticks” test kit), detection, and pain thresholds for cutaneous electrical stimuli (applied to the forearm) were tested at baseline and after four weeks. Twenty-eight patients exposed themselves to four odors (rose, vanilla, chocolate, peach) every two hours over a period of four weeks (training group). Control patients (N = 14) underwent no such “olfactory training” (nontraining group).

Results

Pain thresholds were significantly increased in patients who performed olfactory training compared with patients who did not train with odors. Detection thresholds and olfactory function remained unchanged.

Conclusions

The present results indicate that regular exposure to odors increases pain thresholds in patients with chronic back pain and could be useful for general pain control in these patients. Furthermore, olfactory training in chronic pain patients might help to reduce chronification of pain by desensitization.

High-Throughput Segmentation of Tiled Biological Structures using Random-Walk Distance Transforms

Various 3D imaging techniques are routinely used to examine biological materials, the results of which are usually a stack of grayscale images. In order to quantify structural aspects of the biological materials, however, they must first be extracted from the dataset in a process called segmentation. If the individual structures to be extracted are in contact or very close to each other, distance-based segmentation methods utilizing the Euclidean distance transform are commonly employed. Major disadvantages of the Euclidean distance transform, however, are its susceptibility to noise (very common in biological data), which often leads to incorrect segmentations (i.e., poor separation of objects of interest), and its limitation of being only effective for roundish objects. In the present work, we propose an alternative distance transform method, the random-walk distance transform, and demonstrate its effectiveness in high-throughput segmentation of three microCT datasets of biological tilings (i.e., structures composed of a large number of similar repeating units). In contrast to the Euclidean distance transform, the random-walk approach represents the global, rather than the local, geometric character of the objects to be segmented and, thus, is less susceptible to noise. In addition, it is directly applicable to structures with anisotropic shape characteristics. Using three case studies—tessellated cartilage from a stingray, the dermal endoskeleton of a starfish, and the prismatic layer of a bivalve mollusc shell—we provide a typical workflow for the segmentation of tiled structures, describe core image processing concepts that are underused in biological research, and show that for each study system, large amounts of biologically-relevant data can be rapidly segmented, visualized, and analyzed.

Targeted killing of TNFR2-expressing tumor cells and Tregs by TNFR2 antagonistic antibodies in advanced Sézary syndrome

Sézary syndrome (SS) is a rare form of cutaneous T-cell lymphoma often refractory to treatment. SS is defined as adenopathy, erythroderma with high numbers of atypical T cells. This offers an opportunity for new interventions and perhaps antibody-based therapeutic by virtue of its high expression of the TNFR2 oncogene on the tumor cells and on T-regulatory cells (T_regs). Potent human-directed TNFR2 antagonistic antibodies have been created that preferentially target the TNFR2 oncogene and tumor-infiltrating TNFR2⁺ T_regs. Here we test the therapeutic potential of TNFR2 antagonists on freshly isolated lymphocytes from patients with Stage IVA SS and from healthy controls. SS patients were on a variety of end-stage multi-drug therapies. Baseline burden T_reg/T effector (T_eff) ratios and the responsiveness of tumor and infiltrating T_regs to TNFR2 antibody killing was studied. We show dose-escalating concentrations of a dominant TNFR2 antagonistic antibody killed TNFR2⁺ SS tumor cells and thus restored CD26^- subpopulations of lymphocyte cell numbers to normal. The abundant TNFR2⁺ T_regs of SS subjects are also killed with TNFR2 antagonism. Beneficial and rapid expansion of T_eff was observed. The combination of T_reg inhibition and T_eff expansion brought the high T_reg/T_eff ratio to normal. Our findings suggest a marked responsiveness of SS tumor cells and T_regs, to targeting with TNFR2 antagonistic antibodies. These results show TNFR2 antibodies are potent and efficacious in vitro.

Interactive Visualization of RNA and DNA Structures

The analysis and visualization of nucleic acids (RNA and DNA) is playing an increasingly important role due to their fundamental importance for all forms of life and the growing number of known 3D structures of such molecules. The great complexity of these structures, in particular, those of RNA, demands interactive visualization to get deeper insights into the relationship between the 2D secondary structure motifs and their 3D tertiary structures. Over the last decades, a lot of research in molecular visualization has focused on the visual exploration of protein structures while nucleic acids have only been marginally addressed. In contrast to proteins, which are composed of amino acids, the ingredients of nucleic acids are nucleotides. They form structuring patterns that differ from those of proteins and, hence, also require different visualization and exploration techniques. In order to support interactive exploration of nucleic acids, the computation of secondary structure motifs as well as their visualization in 2D and 3D must be fast. Therefore, in this paper, we focus on the performance of both the computation and visualization of nucleic acid structure. We present a ray casting-based visualization of RNA and DNA secondary and tertiary structures, which enables for the first time real-time visualization of even large molecular dynamics trajectories. Furthermore, we provide a detailed description of all important aspects to visualize nucleic acid secondary and tertiary structures. With this, we close an important gap in molecular visualization.

The Presentation of Olfactory-Trigeminal Mixed Stimuli Increases the Response to Subsequent Olfactory Stimuli

The aim of this study was to evaluate the effect of (1) the addition of trigeminal stimuli to an olfactory stimulus and (2) the congruence in the odorous mixture after repeated odor presentation. Twenty-five normosmic volunteers were enrolled and presented stimulation blocks, consisting of three habituation stimuli (H) (orange odor), one dishabituation (DH) (control condition, orange odor; congruent condition, orange odor + CO₂; incongruent condition, orange odor + l-isopulegol), and one dishabituated stimulus (D) (orange odor). Olfactory event-related potentials were analyzed. Response amplitudes differed significantly in the incongruent condition (N1P2 between H3 and D; peak to peak N1P2 at electrode positions Cz, Fz, and Pz; response amplitudes between H3 and DH). The addition of CO₂ modified the perception of orange odor, pronouncing a fruity note, whereas the addition of l-isopulegol as a DH pronounced the l-isopulegol note. This study provides evidence that incongruent trigeminal-olfactory stimulants increase the response to subsequent olfactory stimulus.

Automated segmentation of complex patterns in biological tissues: Lessons from stingray tessellated cartilage

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.

The Rewarding Effect of Pictures with Positive Emotional Connotation upon Perception and Processing of Pleasant Odors-An FMRI Study

This fMRI study was designed to investigate the effect of cross-modal conditioning in 28 female volunteers. Subjects underwent initial fMRI block design scanning during which three pleasant olfactory stimuli were presented and had to be rated with respect to intensity and pleasantness. This was followed by an odor identification task spread out over 3 days: the experimental group was rewarded for successful trials (correct odor identification) with emotionally salient photos, whilst the control group only received randomly displayed, emotionally neutral, pictures. In the final scanning session, the odors were again presented, and subjects rated pleasantness and intensity. Both pleasantness ratings and fMRI data showed effects of the rewarding procedure. Activation in nucleus accumbens and the orbitofrontal cortex confirmed the hypothesis that learnt association of odors with visual stimuli of emotionally positive valence not only increases pleasantness of the olfactory stimuli but is also reflected in the activation of brain structures relevant for hedonic and reward processing. To our knowledge, this is the first paper to report successful cross-modal conditioning of olfactory stimuli with visual clues.

Nature of the coupling between neural drive and force-generating capacity in the human quadriceps muscle

The force produced by a muscle depends on both the neural drive it receives and several biomechanical factors. When multiple muscles act on a single joint, the nature of the relationship between the neural drive and force-generating capacity of the synergistic muscles is largely unknown. This study aimed to determine the relationship between the ratio of neural drive and the ratio of muscle force-generating capacity between two synergist muscles (vastus lateralis (VL) and vastus medialis (VM)) in humans. Twenty-one participants performed isometric knee extensions at 20 and 50% of maximal voluntary contractions (MVC). Myoelectric activity (surface electromyography (EMG)) provided an index of neural drive. Physiological cross-sectional area (PCSA) was estimated from measurements of muscle volume (magnetic resonance imaging) and muscle fascicle length (three-dimensional ultrasound imaging) to represent the muscles' force-generating capacities. Neither PCSA nor neural drive was balanced between VL and VM. There was a large (r = 0.68) and moderate (r = 0.43) correlation between the ratio of VL/VM EMG amplitude and the ratio of VL/VM PCSA at 20 and 50% of MVC, respectively. This study provides evidence that neural drive is biased by muscle force-generating capacity, the greater the force-generating capacity of VL compared with VM, the stronger bias of drive to the VL.

Two new species of erect Bryozoa (Gymnolaemata: Cheilostomata) and the application of non-destructive imaging methods for quantitative taxonomy

Two new species of cheilostome Bryozoa are described from continental-slope habitats off Mauritania, including canyon and cold-water coral (mound) habitats. Internal structures of both species were visualised and quantified using microcomputed tomographic (micro-CT) methods. Cellaria bafouri n. sp. is characterised by the arrangement of zooids in alternating longitudinal rows, a smooth cryptocyst, and the presence of an ooecial plate with denticles. Smittina imragueni n. sp. exhibits many similarities with Smittina cervicornis (Pallas, 1766), but differs especially in the shape and orientation of the suboral avicularium. Observations on Smittina imragueni and material labelled as Smittina cervicornis suggest that the latter represents a species group, members of which have not yet been discriminated, possibly because of high intracolony variation and marked astogenetic changes in surface morphology. Both new species are known only from the habitats where they were collected, probably reflecting the paucity of bryozoan sampling from this geographic area and depth range. Both species are able to tolerate low oxygen concentration, which is assumed to be compensated by the high nutrient supply off Mauritania. The application of micro-CT for the semiautomatic quantification of zooidal skeletal characters was successfully tested. We were able to automatically distinguish individual zooidal cavities and acquire corresponding morphological datasets. Comparing the obtained results with conventional SEM measurements allowed ascertaining the reliability of this new method. The employment of micro-CT allows the observation and quantification of previously unseen characters that can be used in describing and differentiating species that were previously indistinguishable. Furthermore, this method might help elucidate processes of colony growth and the function of individual zooids during this process.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Registering 2D and 3D imaging data of bone during healing

PURPOSE/AIMS OF THE STUDY: Bone's hierarchical structure can be visualized using a variety of methods. Many techniques, such as light and electron microscopy generate two-dimensional (2D) images, while micro-computed tomography (µCT) allows a direct representation of the three-dimensional (3D) structure. In addition, different methods provide complementary structural information, such as the arrangement of organic or inorganic compounds. The overall aim of the present study is to answer bone research questions by linking information of different 2D and 3D imaging techniques. A great challenge in combining different methods arises from the fact that they usually reflect different characteristics of the real structure.

MATERIALS AND METHODS

We investigated bone during healing by means of µCT and a couple of 2D methods. Backscattered electron images were used to qualitatively evaluate the tissue's calcium content and served as a position map for other experimental data. Nanoindentation and X-ray scattering experiments were performed to visualize mechanical and structural properties.

RESULTS

We present an approach for the registration of 2D data in a 3D µCT reference frame, where scanning electron microscopies serve as a methodic link. Backscattered electron images are perfectly suited for registration into µCT reference frames, since both show structures based on the same physical principles. We introduce specific registration tools that have been developed to perform the registration process in a semi-automatic way.

CONCLUSIONS

By applying this routine, we were able to exactly locate structural information (e.g. mineral particle properties) in the 3D bone volume. In bone healing studies this will help to better understand basic formation, remodeling and mineralization processes.

Ligand Excluded Surface: A New Type of Molecular Surface

The most popular molecular surface in molecular visualization is the solvent excluded surface (SES). It provides information about the accessibility of a biomolecule for a solvent molecule that is geometrically approximated by a sphere. During a period of almost four decades, the SES has served for many purposes - including visualization, analysis of molecular interactions and the study of cavities in molecular structures. However, if one is interested in the surface that is accessible to a molecule whose shape differs significantly from a sphere, a different concept is necessary. To address this problem, we generalize the definition of the SES by replacing the probe sphere with the full geometry of the ligand defined by the arrangement of its van der Waals spheres. We call the new surface ligand excluded surface (LES) and present an efficient, grid-based algorithm for its computation. Furthermore, we show that this algorithm can also be used to compute molecular cavities that could host the ligand molecule. We provide a detailed description of its implementation on CPU and GPU. Furthermore, we present a performance and convergence analysis and compare the LES for several molecules, using as ligands either water or small organic molecules.

Automated stitching of microtubule centerlines across serial electron tomograms

Tracing microtubule centerlines in serial section electron tomography requires microtubules to be stitched across sections, that is lines from different sections need to be aligned, endpoints need to be matched at section boundaries to establish a correspondence between neighboring sections, and corresponding lines need to be connected across multiple sections. We present computational methods for these tasks: 1) An initial alignment is computed using a distance compatibility graph. 2) A fine alignment is then computed with a probabilistic variant of the iterative closest points algorithm, which we extended to handle the orientation of lines by introducing a periodic random variable to the probabilistic formulation. 3) Endpoint correspondence is established by formulating a matching problem in terms of a Markov random field and computing the best matching with belief propagation. Belief propagation is not generally guaranteed to converge to a minimum. We show how convergence can be achieved, nonetheless, with minimal manual input. In addition to stitching microtubule centerlines, the correspondence is also applied to transform and merge the electron tomograms. We applied the proposed methods to samples from the mitotic spindle in C. elegans, the meiotic spindle in X. laevis, and sub-pellicular microtubule arrays in T. brucei. The methods were able to stitch microtubules across section boundaries in good agreement with experts' opinions for the spindle samples. Results, however, were not satisfactory for the microtubule arrays. For certain experiments, such as an analysis of the spindle, the proposed methods can replace manual expert tracing and thus enable the analysis of microtubules over long distances with reasonable manual effort.

The segmentation of microtubules in electron tomograms using Amira

The development of automatic tools for the three-dimensional reconstruction of the microtubule cytoskeleton is crucial for large-scale analysis of mitotic spindles. Recently, we have published a method for the semiautomatic tracing of microtubules based on 3D template matching (Weber et al., J Struct Biol 178:129-138, 2012). Here, we give step-by-step instructions for the automatic tracing of microtubules emanating from centrosomes in the early mitotic Caenorhabditis elegans embryo. This approach, integrated in the visualization and data analysis software Amira, is applicable to tomographic data sets from other model systems.

Exploring cavity dynamics in biomolecular systems

Background

The internal cavities of proteins are dynamic structures and their dynamics may be associated with conformational changes which are required for the functioning of the protein. In order to study the dynamics of these internal protein cavities, appropriate tools are required that allow rapid identification of the cavities as well as assessment of their time-dependent structures.

Results

In this paper, we present such a tool and give results that illustrate the applicability for the analysis of molecular dynamics trajectories. Our algorithm consists of a pre-processing step where the structure of the cavity is computed from the Voronoi diagram of the van der Waals spheres based on coordinate sets from the molecular dynamics trajectory. The pre-processing step is followed by an interactive stage, where the user can compute, select and visualize the dynamic cavities. Importantly, the tool we discuss here allows the user to analyze the time-dependent changes of the components of the cavity structure. An overview of the cavity dynamics is derived by rendering the dynamic cavities in a single image that gives the cavity surface colored according to its time-dependent dynamics.

Conclusion

The Voronoi-based approach used here enables the user to perform accurate computations of the geometry of the internal cavities in biomolecules. For the first time, it is possible to compute dynamic molecular paths that have a user-defined minimum constriction size. To illustrate the usefulness of the tool for understanding protein dynamics, we probe the dynamic structure of internal cavities in the bacteriorhodopsin proton pump.

Anisotropic sampling of planar and two-manifold domains for texture generation and glyph distribution

We present a new method for the generation of anisotropic sample distributions on planar and two-manifold domains. Most previous work that is concerned with aperiodic point distributions is designed for isotropically shaped samples. Methods focusing on anisotropic sample distributions are rare, and either they are restricted to planar domains, are highly sensitive to the choice of parameters, or they are computationally expensive. In this paper, we present a time-efficient approach for the generation of anisotropic sample distributions that only depends on intuitive design parameters for planar and two-manifold domains. We employ an anisotropic triangulation that serves as basis for the creation of an initial sample distribution as well as for a gravitational-centered relaxation. Furthermore, we present an approach for interactive rendering of anisotropic Voronoi cells as base element for texture generation. It represents a novel and flexible visualization approach to depict metric tensor fields that can be derived from general tensor fields as well as scalar or vector fields.

Automated tracing of microtubules in electron tomograms of plastic embedded samples of Caenorhabditis elegans embryos

The ability to rapidly assess microtubule number in 3D image stacks from electron tomograms is essential for collecting statistically meaningful data sets. Here we implement microtubule tracing using 3D template matching. We evaluate our results by comparing the automatically traced centerlines to manual tracings in a large number of electron tomograms of the centrosome of the early Caenorhabditis elegans embryo. Furthermore, we give a qualitative description of the tracing results for three other types of samples. For dual-axis tomograms, the automatic tracing yields 4% false negatives and 8% false positives on average. For single-axis tomograms, the accuracy of tracing is lower (16% false negatives and 14% false positives) due to the missing wedge in electron tomography. We also implemented an editor specifically designed for correcting the automatic tracing. Besides, this editor can be used for annotating microtubules. The automatic tracing together with a manual correction significantly reduces the amount of manual labor for tracing microtubule centerlines so that large-scale analysis of microtubule network properties becomes feasible.

Voronoi-based extraction and visualization of molecular paths

Visual analysis is widely used to study the behavior of molecules. Of particular interest are the analysis of molecular interactions and the investigation of binding sites. For large molecules, however, it is difficult to detect possible binding sites and paths leading to these sites by pure visual inspection. In this paper, we present new methods for the computation and visualization of potential molecular paths. Using a novel filtering method, we extract the significant paths from the Voronoi diagram of spheres. For the interactive visualization of molecules and their paths, we present several methods using deferred shading and other state-of-the-art techniques. To allow for a fast overview of reachable regions of the molecule, we illuminate the molecular surface using a large number of light sources placed on the extracted paths. We also provide a method to compute the extension surface of selected paths and visualize it using the skin surface. Furthermore, we use the extension surface to clip the molecule to allow easy visual tracking of even deeply buried paths. The methods are applied to several proteins to demonstrate their usefulness.

Calibration of a microchannel plate based extreme ultraviolet grazing incident spectrometer at the Advanced Light Source

We present the design and calibration of a microchannel plate based extreme ultraviolet spectrometer. Calibration was performed at the Advance Light Source (ALS) at the Lawrence Berkeley National Laboratory (LBNL). This spectrometer will be used to record the single shot spectrum of radiation emitted by the tapered hybrid undulator (THUNDER) undulator installed at the LOASIS GeV-class laser-plasma-accelerator. The spectrometer uses an aberration-corrected concave grating with 1200 lines/mm covering 11-62 nm and a microchannel plate detector with a CsI coated photocathode for increased quantum efficiency in the extreme ultraviolet. A touch screen interface controls the grating angle, aperture size, and placement of the detector in vacuum, allowing for high-resolution measurements over the entire spectral range.

The injectable cyclooxygenase-2-specific inhibitor parecoxib sodium has analgesic efficacy when administered preoperatively

Preoperative administration of analgesics may prevent or reduce hyperalgesia and inhibit inflammation and pain by reducing the synthesis of prostaglandins in response to surgical injury. We evaluated in this placebo-controlled study the analgesic efficacy and safety of single doses of parecoxib sodium (20, 40, and 80 mg IV) when administered before oral surgery. Efficacy assessments were recorded during the 24-h period after completion of surgery. All doses of parecoxib sodium were consistently and significantly superior to placebo as measured by time to rescue medication, proportion of patients requiring rescue medication, patient's global assessment, and pain intensity. There were no significant differences between the Parecoxib Sodium 40- and 80-mg groups, suggesting that the analgesic effect of preoperatively administered parecoxib sodium reaches a plateau at 40 mg in this model. Forty-eight percent of the Parecoxib Sodium 40-mg group required rescue medication in the 24-h study period, compared with 93% of patients in the Placebo group. Overall, there were fewer adverse events in parecoxib sodium-treated patients compared with placebo. These findings suggest that preoperative administration of parecoxib sodium, the injectable prodrug of the cyclooxygenase-2 specific inhibitor valdecoxib, is effective, safe, and well tolerated for treating postoperative pain.

Comparison of wounds created by non-bladed trocars and pyramidal tip trocars in the pig

Wounds made by the Endopath nonbladed obturator, the Step trocar, and conventional pyramidal tip trocars were compared. The endopath nonbladed obturator and the Step trocar made wounds by separating tissue fibers, whereas the pyramidal tip trocar cut tissue fibers. The wounds of the Endopath nonbladed obturator and the Step trocar were similar in length but were narrower than wounds made by the pyramidal tip trocar. Further studies are needed to determine whether the wounds made by the Endopath nonbladed obturator and the Step trocar will have fewer complications than conventional pyramidal tip trocars.

TNF-alpha, not CD154 (CD40L), plays a major role in SEB-dependent, CD4(+) T cell-induced endothelial cell activation in vitro

CD4(+) T cell effector molecules, in particular TNF-alpha and CD154, activate endothelial cells. However, the relative contributions of TNF-alpha and CD154 in mediating endothelial cell activation during complex Ag-driven CD4(+) T cell-endothelial cell interactions are not known. We utilized an in vitro model of CD4(+) T cell-endothelial cell interactions to characterize the contributions of TNF-alpha and CD154 in mediating upregulation of adhesion molecules CD54, CD62E, and CD106 on human umbilical vein endothelial cells (HUVEC). HUVEC were first treated with IFN-gamma to upregulate MHC Class II expression. IFN-gamma minimally effects HUVEC adhesion molecule expression but renders them capable of MHC class II restricted interactions with CD4(+) T cells. Coculturing MHC class II+ HUVEC and CD4(+) T cells with the superantigen SEB induces a rapid and marked upregulation of CD54, CD62E, and CD106 expression on HUVEC, as shown by FACS analysis. To study the effector molecules mediating SEB-driven, CD4(+) T cell-dependent endothelial cell activation, similar experiments were performed in the presence of neutralizing anti-CD154, anti-TNF-alpha, or anti-IL1 antibodies, as well as combinations of these antibodies. In contrast to the anti-CD154 or anti-IL-1 antibodies, the anti-TNF-alpha mAb markedly inhibited SEB-dependent, CD4(+) T cell-induced HUVEC activation. We conclude that TNF-alpha, not CD154, plays the major role in SEB-driven, CD4(+) T cell-induced endothelial cell activation in vitro.

Circumscription of the Malvales and relationships to other Rosidae: evidence from rbcL sequence data

The order Malvales remains poorly circumscribed, despite its seemingly indisputable core constituents: Bombacaceae, Malvaceae, Sterculiaceae, and Tiliaceae. We conducted a two-step parsimony analysis on 125 rbcL sequences to clarify the composition of Malvales, to determine the relationships of some controversial families, and to identify the placement of the Malvales within Rosidae. We sampled taxa that have been previously suggested to be within, or close to, Malvales (83 sequences), plus additional rosids (26 sequences) and nonrosid eudicots (16 sequences) to provide a broader framework for the analysis. The resulting trees strongly support the monophyly of the core malvalean families, listed above. In addition, these data serve to identify a broader group of taxa that are closely associated with the core families. This expanded malvalean clade is composed of four major subclades: (1) the core families (Bombacaceae, Malvaceae, Sterculiaceae, Tiliaceae); (2) Bixaceae, Cochlospermaceae, and Sphaerosepalaceae (Rhopalocarpaceae); (3) Thymelaeaceae sensu lato (s.l.); and (4) Cistaceae, Dipterocarpaceae s.l., Sarcolaenaceae (Chlaenaceae), and Muntingia. In addition, Neurada (Neuradaceae or Rosaceae) falls in the expanded malvalean clade but not clearly within any of the four major subclades. This expanded malvalean clade is sister to either the expanded capparalean clade of Rodman et al. or the sapindalean clade of Gadek et al. Members of Elaeocarpaceae, hypothesized by most authors as a sister group to the four core malvalean families, are shown to not fall close to these taxa. Also excluded as members of, or sister groups to, the expanded malvalean clade were the families Aextoxicaceae, Barbeyaceae, Cannabinaceae, Cecropiaceae, Dichapetalaceae, Elaeagnaceae, Euphorbiaceae s.l., Huaceae, Lecythidaceae, Moraceae s.l., Pandaceae, Plagiopteraceae, Rhamnaceae, Scytopetalaceae, Ulmaceae, and Urticaceae.