Marker-free image registration of electron tomography tilt-series

Background  

Tilt series are commonly used in electron tomography as a means of collecting three-dimensional information from two-dimensional projections. A common problem encountered is the projection alignment prior to 3D reconstruction. Current alignment techniques usually employ gold particles or image derived markers to correctly align the images. When these markers are not present, correlation between adjacent views is used to align them. However, sequential pairwise correlation is prone to bias and the resulting alignment is not always optimal.     

2D fast rotational matching for image processing of biophysical data

In 3D single particle reconstruction, which involves the translational and rotational matching of a large number of electron microscopy (EM) images, the algorithmic performance is largely dependent on the efficiency and accuracy of the underlying 2D image alignment kernel. We present a novel fast rotational matching kernel for 2D images (FRM2D) that significantly reduces the cost of this alignment. The alignment problem is formulated using one translational and two rotational degrees of freedom. This allows us to take advantage of fast Fourier transforms (FFTs) in rotational space to accelerate the search of the two angular parameters, while the remaining translational parameter is explored, within a limited range, by exhaustive search. Since there are no boundary effects in FFTs of cyclic angular variables, we avoid the expensive zero padding associated with Fourier transforms in linear space. To verify the robustness of our method, efficiency and accuracy tests were carried out over a range of noise levels in realistic simulations of EM images. Performance tests against two standard alignment methods, resampling to polar coordinates and self-correlation, demonstrate that FRM2D compares very favorably to the traditional methods. FRM2D exhibits a comparable or higher robustness against noise and a significant gain in efficiency that depends on the fineness of the angular sampling and linear search range.     

Multiphase method for automatic alignment of transmission electron microscope images using markers

In order to successfully perform the 3D reconstruction in electron tomography, transmission electron microscope images must be accurately aligned or registered. So far, the problem is solved by either manually showing the corresponding fiducial markers from the set of images or automatically using simple correlation between the images on several rotations and scales. The present solutions, however, share the problem of being inefficient and/or inaccurate. We therefore propose a method in which the registration is automated using conventional colloidal gold particles as reference markers between images. We approach the problem from the computer vision viewpoint; hence, the alignment problem is divided into several subproblems: (1) finding initial matches from successive images, (2) estimating the epipolar geometry between consecutive images, (3) finding and localizing the gold particles with subpixel accuracy in each image, (4) predicting the probable matching gold particles using the epipolar constraint and its uncertainty, (5) matching and tracking the gold beads through the tilt series, and (6) optimizing the transformation parameters for the whole image set. The results show not only the reliability of the suggested method but also a high level of accuracy in alignment, since practically all the visible gold markers can be used.     

Central line bundle implementation in US intensive care units and impact on bloodstream infections

Background

Central line-associated bloodstream infections (CLABSI) represent a serious patient safety issue. To prevent these infections, bundled interventions are increasingly recommended. We examine the extent of adoption of Central Line (CL) Bundle elements throughout US intensive care units (ICU) and determine their effectiveness in preventing CLABSIs.

Methodology/Principal Findings

In this cross-sectional study, National Healthcare Safety Network (NHSN) hospitals provided the following: ICU-specific NHSN-reported rates of CLABSI/1,000 central line days; policies and compliance rates regarding bundle components; and other setting characteristics. In 250 hospitals the mean CLABSI rate was 2.1 per 1000 central line days and 49% reported having a written CL Bundle policy. However, of those that monitored compliance, only 38% reported very high compliance with the CL Bundle. Only when an ICU had a policy, monitored compliance, and had ≥95% compliance did CLABSI rates decrease. Complying with any one of three CL Bundle elements resulted in decreased CLABSI rates (β = -1.029, p = 0.015). If an ICU without good bundle compliance achieved high compliance with any one bundle element, we estimated that its CLABSI rate would decrease by 38%.

Conclusions/Significance

In NHSN hospitals across the US, the CL Bundle is associated with lower infection rates only when compliance is high. Hospitals must target improving bundle implementation and compliance as opposed to simply instituting policies.     

Single particle analysis of filamentous and highly elongated macromolecular assemblies

The application of single particle techniques to the three-dimensional analysis of electron microscope images of elongated or filamentous macromolecular assemblies is evaluated, taking as an example the muscle thin filament. Although the thin filament contains local helical symmetry, because of the inherent variable twist along it, the helical coherence does not extend for large enough distances to allow the symmetry to be used for full reconstruction of the tropomyosin/troponin repeat along the filament. The muscle thin filament therefore represents a general case of a filamentous object in that it is not possible to exploit symmetry in a full analysis. Due to the nature of the imaging process in the electron microscope, only projections of the thin filament around its long axis are available without tilting the grid. Crucially, projection images around a single axis do not provide enough information to assign Euler angles ab initio using current methods. Tests with a model thin filament structure indicated that an out-of-plane tilt of approximately 20 degrees was needed for ab initio angular assignment of sufficient accuracy to calculate a 3D structure to a resolution of approximately 25 A. If no out-of-plane views are available, an alternative approach is to use a prior 3D model as a reference for the initial angle assignment. Tests with the thin filament model indicated that reasonably accurate angular assignment can be made using a reference containing actin, but lacking the regulatory proteins tropomyosin and troponin. We also found that an adaptation of the exact filtered back projection method is required to allow the correct weighting of projection images in which the particle has a very large axial ratio. This adaptation resulted in significant improvements in the reconstruction.     

Centripetal Disposition of Bundle Sheath Chloroplasts During the Leaf Development of Eleusine coracana

Distribution of chloroplasts in bundle sheath cells was examinedby light and electron microscopy during the leaf developmentof finger millet (Eleusine coracana Gaertn.), an NAD malic enzymetype C₄ plant with centripetal arrangement of bundle sheathchloroplasts. Young chloroplasts are almost evenly distributedalong the cell walls in bundle sheath cells of folded immatureleaves. In elongating leaves and above the elongation zone thebundle sheath chloroplasts tend to lie along the radial wallsand the walls adjacent to the vascular bundle. They furthermigrate near to the vascular bundle and finally establish acentripetal arrangement. Mitochondria, microbodies and nucleusmigrate along with the chloroplasts. Etioplasts and other organellesare centripetally located in the bundle sheath cells of etiolatedseedlings grown in the dark. Bundle sheath chloroplast, C₄ plant, chloroplast, chloroplast orientation, Eleusine coracana, finger millet     

Hill Reaction, Hydrogen Peroxide Scavenging, and Ascorbate Peroxidase Activity of Mesophyll and Bundle Sheath Chloroplasts of NADP-Malic Enzyme Type C(4) Species

Intact mesophyll and bundle sheath chloroplasts wee isolated from the NADP-malic enzyme type C₄ plants maize, sorghum (monocots), and Flaveria trinervia (dicot) using enzymic digestion and mechanical isolation techniques. Bundle sheath chloroplasts of this C₄ subgroup tend to be agranal and were previously reported to be deficient in photosystem II activity. However, following injection of intact bundle sheath chloroplasts into hypotonic medium, thylakoids had high Hill reaction activity, similar to that of mesophyll chloroplasts with the Hill oxidants dichlorophenolindophenol, p-benzoquinone, and ferricyanide (approximately 200 to 300 micromoles O₂ evolved per mg chlorophyll per hour). In comparison to that of mesophyll chloroplasts, the Hill reaction activity of bundle sheath chloroplasts of maize and sorghum was labile and lost activity during assay. Bundle sheath chloroplasts of maize also exhibited some capacity for 3-phosphoglycerate dependent O₂ evolution (29 to 58 micromoles O₂ evolved per milligram chlorophyll per hour). Both the mesophyll and bundle sheath chloroplasts were equally effective in light dependent scavenging of hydrogen peroxide. The results suggest that both chloroplast types have noncyclic electron transport and the enzymology to reduce hydrogen peroxide to water. The activities of ascorbate peroxidase from these chloroplast types was consistent with their capacity to scavenge hydrogen peroxide.     

A microcomputer based system for three-dimensional reconstructions from tomographic or histologic sections

The reconstructions of three-dimensional (3-D) objects from serial two-dimensional (2-D) images can contribute to the understanding of many biologic structures, from organelles to organs and tissues. The 3-D reconstruction of sections can be divided into several major tasks: image acquisition, alignment of slices, internal object definition, object reconstruction and rotation of the completed image. A fast, versatile, interactive system was devised for the reconstruction of 3-D objects from serial 2-D images using a low-cost microcomputer, original programs and commercial software. The system allows reconstruction from any serial images, e.g., electron micrographs, histologic sections or computed tomograms. A photographic image or a microscopic field is acquired into the computer memory using a video digitizer. Slices are superimposed and aligned to each other using an operator-interactive program. A contour-(edge-) finding algorithm isolates an object of interest from the background image by "subtraction" of the image from an overlaid, slightly shifted identical image. Contours for each slice are input to a reconstruction procedure, which calculates the x, y and z coordinates of every point in a slice and the thickness and number of slices. It then calculates the illumination for every point using a given point source of light and an intensity-fading coefficient. Finally, the points are represented by cubes to provide dimension and reflective surfaces. A cube of appropriate shade and color represents in 2-D the equivalent of a 3-D object; this results in a very effective 3-D image. The reconstruction is rotated by recalculating the positions of every point defining the object and rebuilding the image.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure of actin filament bundles from microvilli of sea urchin eggs

The detailed substructure of actin filament bundles in microvilli of fertilized sea urchin eggs has been studied by analysing electron microscope images of negatively stained specimens. Transverse stripes which repeat about every 130 Å along the axis of a bundle, as previously observed by Burgess & Schroeder (1977), reflect the positions of cross-bridges that connect the filaments into a bundle. Analysis of optical transforms of the micrographs reveals that there are approximately 14 actin monomers between cross-overs of the two long-pitch helical strands of the actin filaments, with three cross-bridges in this interval. The structure is basically similar to that of the hexagonally packed bundles prepared in vitro from high speed supernatants of sea urchin eggs by Kane (1975) and analyzed by DeRosier et al. (1977). One clear difference, however, is that the in vivo microvillar filament bundles are supercoiled, giving rise to long axial repeats of 1500 to 2000 Å.Computationally filtered images of regions that were only slightly supercoiled reveal the relative alignment of filaments within the bundles and show that crossbridges appear to interact with four actin monomers, apparently linking two actin monomers on one strand of one filament to the nearest two monomers on a neighbouring filament. However, the cross-bridges are not spaced at equal intervals corresponding to four actin subunits, presumably because of the lack of hexagonal symmetry in the individual filaments, which have about 14 actin monomers between cross-overs. Instead, the cross-bridges are arranged quasiequivalently along the longitudinal axis of the bundles, in steps of four or five actin subunit spacings (28 Å each).     

A Strategy for Electron Tomographic Data Collection and Crystallographic Reconstruction of Biological Bundles

Structures of highly ordered biological bundles have unique features which call for special experimental and computational methods in electron cryomicroscopy. They can be considered as three-dimensional quasi-crystals and reconstructed using a crystallographic approach. However, they are neither “infinitely” large with respect to the borders of the bundle, nor are they a single unit cell in thickness along the viewing direction. Also, because of their shape, bundles do not generally have a preferred azimuthal orientation, which poses challenges for orientation estimation and refinement. We developed a strategy for recording and processing electron cryomicroscopic images that differs from classical two-dimensional crystalline reconstruction techniques. These developments allowed us to merge data from tomographic tilt series of ice-embedded acrosomal bundles. The goal is to determine accurately amplitudes and phases at the diffraction maxima in terms ofhklindices, and compute a three-dimensional map from the diffraction data.     

The effect of overabundant projection directions on 3D reconstruction algorithms

The experimental process of collecting images from macromolecules in an electron microscope is such that it does not allow for prior specification of the angular distribution of the projection images. As a consequence, an uneven distribution of projection directions may occur. Concerns have been raised recently about the behavior of 3D reconstruction algorithms for the case of unevenly distributed projections. It has been illustrated on experimental data that in the case of a heavily uneven distribution of projection directions some algorithms tend to elongate the reconstructed volumes along the overloaded direction so much as to make a quantitative biological analysis impossible. In answer to these concerns we have developed a strategy for quantitative comparison and optimization of 3D reconstruction algorithms. We apply this strategy to quantitatively analyze algebraic reconstruction techniques (ART) with blobs, simultaneous iterative reconstruction techniques (SIRT) with voxels, and weighted backprojection (WBP). We show that the elongation artifacts that had been previously reported can be strongly reduced. With our specific choices for the free parameters of the three algorithms, WBP reconstructions tend to be inferior to those obtained with either SIRT or ART and the results obtained with ART are comparable to those with SIRT, but at a very small fraction of the computational cost of SIRT.     

The impact of multiple protein sequence alignment on phylogenetic estimation

Multiple sequence alignment is typically the first step in estimating phylogenetic trees, with the assumption being that as alignments improve, so will phylogenetic reconstructions. Over the last decade or so, new multiple sequence alignment methods have been developed to improve comparative analyses of protein structure, but these new methods have not been typically used in phylogenetic analyses. In this paper, we report on a simulation study that we performed to evaluate the consequences of using these new multiple sequence alignment methods in terms of the resultant phylogenetic reconstruction. We find that while alignment accuracy is positively correlated with phylogenetic accuracy, the amount of improvement in phylogenetic estimation that results from an improved alignment can range from quite small to substantial. We observe that phylogenetic accuracy is most highly correlated with alignment accuracy when sequences are most difficult to align, and that variation in alignment accuracy can have little impact on phylogenetic accuracy when alignment error rates are generally low. We discuss these observations and implications for future work.     

Image correlation method for DNA sequence alignment

The complexity of searches and the volume of genomic data make sequence alignment one of bioinformatics most active research areas. New alignment approaches have incorporated digital signal processing techniques. Among these, correlation methods are highly sensitive. This paper proposes a novel sequence alignment method based on 2-dimensional images, where each nucleic acid base is represented as a fixed gray intensity pixel. Query and known database sequences are coded to their pixel representation and sequence alignment is handled as object recognition in a scene problem. Query and database become object and scene, respectively. An image correlation process is carried out in order to search for the best match between them. Given that this procedure can be implemented in an optical correlator, the correlation could eventually be accomplished at light speed. This paper shows an initial research stage where results were "digitally" obtained by simulating an optical correlation of DNA sequences represented as images. A total of 303 queries (variable lengths from 50 to 4500 base pairs) and 100 scenes represented by 100 x 100 images each (in total, one million base pair database) were considered for the image correlation analysis. The results showed that correlations reached very high sensitivity (99.01%), specificity (98.99%) and outperformed BLAST when mutation numbers increased. However, digital correlation processes were hundred times slower than BLAST. We are currently starting an initiative to evaluate the correlation speed process of a real experimental optical correlator. By doing this, we expect to fully exploit optical correlation light properties. As the optical correlator works jointly with the computer, digital algorithms should also be optimized. The results presented in this paper are encouraging and support the study of image correlation methods on sequence alignment.     

Markov random field based automatic image alignment for electron tomography

We present a method for automatic full-precision alignment of the images in a tomographic tilt series. Full-precision automatic alignment of cryo electron microscopy images has remained a difficult challenge to date, due to the limited electron dose and low image contrast. These facts lead to poor signal to noise ratio (SNR) in the images, which causes automatic feature trackers to generate errors, even with high contrast gold particles as fiducial features. To enable fully automatic alignment for full-precision reconstructions, we frame the problem probabilistically as finding the most likely particle tracks given a set of noisy images, using contextual information to make the solution more robust to the noise in each image. To solve this maximum likelihood problem, we use Markov Random Fields (MRF) to establish the correspondence of features in alignment and robust optimization for projection model estimation. The resulting algorithm, called Robust Alignment and Projection Estimation for Tomographic Reconstruction, or RAPTOR, has not needed any manual intervention for the difficult datasets we have tried, and has provided sub-pixel alignment that is as good as the manual approach by an expert user. We are able to automatically map complete and partial marker trajectories and thus obtain highly accurate image alignment. Our method has been applied to challenging cryo electron tomographic datasets with low SNR from intact bacterial cells, as well as several plastic section and X-ray datasets.     

Quantitative underwater 3D motion analysis using submerged video cameras: accuracy analysis and trajectory reconstruction

In this study we aim at investigating the applicability of underwater 3D motion capture based on submerged video cameras in terms of 3D accuracy analysis and trajectory reconstruction. Static points with classical direct linear transform (DLT) solution, a moving wand with bundle adjustment and a moving 2D plate with Zhang's method were considered for camera calibration. As an example of the final application, we reconstructed the hand motion trajectories in different swimming styles and qualitatively compared this with Maglischo's model. Four highly trained male swimmers performed butterfly, breaststroke and freestyle tasks. The middle fingertip trajectories of both hands in the underwater phase were considered. The accuracy (mean absolute error) of the two calibration approaches (wand: 0.96 mm – 2D plate: 0.73 mm) was comparable to out of water results and highly superior to the classical DLT results (9.74 mm). Among all the swimmers, the hands' trajectories of the expert swimmer in the style were almost symmetric and in good agreement with Maglischo's model. The kinematic results highlight symmetry or asymmetry between the two hand sides, intra- and inter-subject variability in terms of the motion patterns and agreement or disagreement with the model. The two outcomes, calibration results and trajectory reconstruction, both move towards the quantitative 3D underwater motion analysis.     

Fiducial-less alignment of cryo-sections

Cryo-electron tomography of vitreous sections is currently the most promising technique for visualizing arbitrary regions of eukaryotic cells or tissue at molecular resolution. Despite significant progress in the sample preparation techniques over the past few years, the three dimensional reconstruction using electron tomography is not as simple as in plunge frozen samples for various reasons, but mainly due to the effects of irradiation on the sections and the resulting poor alignment. Here, we present a new algorithm, which can provide a useful three-dimensional marker model after investigation of hundreds to thousands of observations calculated using local cross-correlation throughout the tilt series. The observations are chosen according to their coherence to a particular model and assigned to virtual markers. Through this type of measurement a merit figure can be calculated, precisely estimating the quality of the reconstruction. The merit figures of this alignment method are comparable to those obtained with plunge frozen samples using fiducial gold markers. An additional advantage of the algorithm is the implicit detection of areas in the sections that behave as rigid bodies and can thus be properly reconstructed.     

Reprint of "Fiducial-less alignment of cryo-sections" [J. Struct. Biol. 159 (2007) 413-423

Cryo-electron tomography of vitreous sections is currently the most promising technique for visualizing arbitrary regions of eukaryotic cells or tissue at molecular resolution. Despite significant progress in the sample preparation techniques over the past few years, the three dimensional reconstruction using electron tomography is not as simple as in plunge frozen samples for various reasons, but mainly due to the effects of irradiation on the sections and the resulting poor alignment. Here, we present a new algorithm, which can provide a useful three-dimensional marker model after investigation of hundreds to thousands of observations calculated using local cross-correlation throughout the tilt series. The observations are chosen according to their coherence to a particular model and assigned to virtual markers. Through this type of measurement a merit figure can be calculated, precisely estimating the quality of the reconstruction. The merit figures of this alignment method are comparable to those obtained with plunge frozen samples using fiducial gold markers. An additional advantage of the algorithm is the implicit detection of areas in the sections that behave as rigid bodies and can thus be properly reconstructed.     

Reduction of surgical site infections after implementation of a bundle of care

Background

Surgical Site Infections (SSI) are relatively frequent complications after colorectal surgery and are associated with substantial morbidity and mortality.

Objective

Implementing a bundle of care and measuring the effects on the SSI rate.

Design

Prospective quasi experimental cohort study.

Methods

A prospective surveillance for SSI after colorectal surgery was performed in the Amphia Hospital, Breda, from January 1, 2008 until January 1, 2012. As part of a National patient safety initiative, a bundle of care consisting of 4 elements covering the surgical process was introduced in 2009. The elements of the bundle were perioperative antibiotic prophylaxis, hair removal before surgery, perioperative normothermia and discipline in the operating room. Bundle compliance was measured every 3 months in a random sample of surgical procedures.

Results

Bundle compliance improved significantly from an average of 10% in 2009 to 60% in 2011. 1537 colorectal procedures were performed during the study period and 300 SSI (19.5%) occurred. SSI were associated with a prolonged length of stay (mean additional length of stay 18 days) and a significantly higher 6 months mortality (Adjusted OR: 2.71, 95% confidence interval 1.76–4.18). Logistic regression showed a significant decrease of the SSI rate that paralleled the introduction of the bundle. The adjusted Odds ratio of the SSI rate was 36% lower in 2011 compared to 2008.

Conclusion

The implementation of the bundle was associated with improved compliance over time and a 36% reduction of the SSI rate after adjustment for confounders. This makes the bundle an important tool to improve patient safety.     

Reconstruction of 3-Dimensional Histology Volume and its Application to Study Mouse Mammary Glands

Histology volume reconstruction facilitates the study of 3D shape and volume change of an organ at the level of macrostructures made up of cells. It can also be used to investigate and validate novel techniques and algorithms in volumetric medical imaging and therapies. Creating 3D high-resolution atlases of different organs^1,2,3 is another application of histology volume reconstruction. This provides a resource for investigating tissue structures and the spatial relationship between various cellular features. We present an image registration approach for histology volume reconstruction, which uses a set of optical blockface images. The reconstructed histology volume represents a reliable shape of the processed specimen with no propagated post-processing registration error. The Hematoxylin and Eosin (H&E) stained sections of two mouse mammary glands were registered to their corresponding blockface images using boundary points extracted from the edges of the specimen in histology and blockface images. The accuracy of the registration was visually evaluated. The alignment of the macrostructures of the mammary glands was also visually assessed at high resolution.This study delineates the different steps of this image registration pipeline, ranging from excision of the mammary gland through to 3D histology volume reconstruction. While 2D histology images reveal the structural differences between pairs of sections, 3D histology volume provides the ability to visualize the differences in shape and volume of the mammary glands.     

Shaping the dynamics of a bidirectional neural interface

Progress in decoding neural signals has enabled the development of interfaces that translate cortical brain activities into commands for operating robotic arms and other devices. The electrical stimulation of sensory areas provides a means to create artificial sensory information about the state of a device. Taken together, neural activity recording and microstimulation techniques allow us to embed a portion of the central nervous system within a closed-loop system, whose behavior emerges from the combined dynamical properties of its neural and artificial components. In this study we asked if it is possible to concurrently regulate this bidirectional brain-machine interaction so as to shape a desired dynamical behavior of the combined system. To this end, we followed a well-known biological pathway. In vertebrates, the communications between brain and limb mechanics are mediated by the spinal cord, which combines brain instructions with sensory information and organizes coordinated patterns of muscle forces driving the limbs along dynamically stable trajectories. We report the creation and testing of the first neural interface that emulates this sensory-motor interaction. The interface organizes a bidirectional communication between sensory and motor areas of the brain of anaesthetized rats and an external dynamical object with programmable properties. The system includes (a) a motor interface decoding signals from a motor cortical area, and (b) a sensory interface encoding the state of the external object into electrical stimuli to a somatosensory area. The interactions between brain activities and the state of the external object generate a family of trajectories converging upon a selected equilibrium point from arbitrary starting locations. Thus, the bidirectional interface establishes the possibility to specify not only a particular movement trajectory but an entire family of motions, which includes the prescribed reactions to unexpected perturbations.