Efficient automatic noise reduction of electron tomographic reconstructions based on iterative median filtering

Self-supervised noise modeling and sparsity guided electron tomography volumetric image denoising

Cryo-Electron Tomography (cryo-ET) is a revolutionary technique for visualizing macromolecular structures in near-native states. However, the physical limitations of imaging instruments lead to cryo-ET volumetric images with very low Signal-to-Noise Ratio (SNR) with complex noise, which has a side effect on the downstream analysis of the characteristics of observed macromolecules. Additionally, existing methods for image denoising are difficult to be well generalized to the complex noise in cryo-ET volumes. In this work, we propose a self-supervised deep learning model for cryo-ET volumetric image denoising based on noise modeling and sparsity guidance (NMSG), achieved by learning the noise distribution in noisy cryo-ET volumes and introducing sparsity guidance to ensure smoothness. Firstly, a Generative Adversarial Network (GAN) is utilized to learn noise distribution in cryo-ET volumes and generate noisy volumes pair from single volume. Then, a new loss function is devised to both ensure the recovery of ultrastructure and local smoothness. Experiments are done on five real cryo-ET datasets and three simulated cryo-ET datasets. The comprehensive experimental results demonstrate that our method can perform reliable denoising by training on single noisy volume, achieving better results than state-of-the-art single volume-based methods and competitive with methods trained on large-scale datasets.

[Progress in filters for denoising cryo-electron microscopy images]

Cryo-electron microscopy (cryo-EM) imaging has the unique potential to bridge the gap between cellular and molecular biology. Therefore, cryo-EM three-dimensional (3D) reconstruction has been rapidly developed in recent several years and applied widely in life science research to reveal the structures of large macromolecular assemblies and cellular complexes, which is critical to understanding their functions at all scales. Although the technical breakthrough in recent years, for example, the introduction of the direct detection device (DDD) camera and the development of cryo-EM software tools, made the three cryo-EM pioneers share the 2017 Nobel Prize, several bottleneck problems still exist that hamper the further increase of the resolution of single-particle reconstruction and hold back the application of in situ subnanometer structure determination by cryo-tomography. Radiation damage is still the key limiting factor in cryo-EM. In order to minimize the radiation damage and preserve as much resolution as possible, the imaging conditions of a low dose and weak contrast make cryo-EM images extremely noisy with very low signal-to-noise ratios (SNR), generally about 0.1. The high noise will obscure the fine details in cryo-EM images or reconstructed maps. Thus, a method to reduce the level of noise and improve the resolution has become an important issue. In this paper, we systematically reviewed and compared some robust filters in the cryo-EM field of two aspects, single-particle analysis (SPA) and cryo-electron tomography (cryo-ET), and especially studied their applications, such as, 3D reconstruction, visualization, structural analysis, and interpretation. Conventional approaches to noise reduction in cryo-EM imaging include the use of Gaussian, median, and bilateral filters, among other means. A Gaussian filter selects an appropriate filter kernel to conduct spatial convolution with a noisy image. Although noise with larger standard deviations in cryo-EM images can be suppressed and satisfactory performance is achieved in certain cases, this filter also blurs the images and over-smooths small-scale image features. This is especially detrimental when precise quantitative information needs to be extracted. Unlike a Gaussian filter, a median filter is based on the order statistics of the image and selects the median intensity in a window of the adjacent pixels to denoise the image. Although this filter is robust to outliers, it suffers from aliasing problems that possibly result in incorrect information for cryo-EM structure interpretation. A bilateral filter is a nonlinear filter that performs spatial weighted averaging and is more selective in the pixels allowing to contribute to the weighted sum, excluding the high frequency noise from the smoothing process. Thus, this filter can be used to smooth out noise while maintaining the edge details, which is similar to an anisotropic diffusion filter, and distinct from a Gaussian filter but its utility will be limited when the SNR of a cryo-EM image is very low. Generally, spatial filtering methods have the disadvantage of losing image resolution when reducing noise. A wavelet transform can exploit the wavelet's natural ability to separate a signal from noise at multiple image scales to allow for joint resolution in both the spatial and frequency domains, and thus has the potential to outperform existing methods. The modified wavelet shrinkage filter we developed can offer a remarkable improvement in image quality with a good compromise between detail preservation and noise smoothing. We expect that our review study on different filters can provide benefits to cryo-EM applications and the interpretation of biological structures.

Local Tension on Talin in Focal Adhesions Correlates with F-Actin Alignment at the Nanometer Scale

Cellular force transmission and mechanotransduction are critical in embryogenesis, normal physiology, and many diseases. Talin plays a key role in these processes by linking integrins to force-generating actomyosin. Using the previously characterized FRET-based talin tension sensor, we observed variations of tension both between and within individual focal adhesions in the same cell. Assembling and sliding adhesions showed gradients with higher talin tension toward the cell center, whereas mature, stable adhesions had uniform talin tension. Total talin accumulation was maximal in high-tension regions; by contrast, vinculin intensity was flat or maximal at the adhesion center, and actin intensity was maximal toward the cell center. To investigate mechanism, we combined talin tension imaging with cellular cryotomography to visualize the correlated actin organization at nanometer resolution. Regions of high talin tension had highly aligned linear actin filaments, whereas regions of low tension had less-well-aligned F-actin. These results reveal an orchestrated spatiotemporal relationship between talin tension, actin/vinculin localization, local actin organization, and focal adhesion dynamics.

Marker-free method for accurate alignment between correlated light, cryo-light, and electron cryo-microscopy data using sample support features

Combining fluorescence microscopy with electron cryo-tomography allows, in principle, spatial localization of tagged macromolecular assemblies and structural features within the cellular environment. To allow precise localization and scale integration between the two disparate imaging modalities, accurate alignment procedures are needed. Here, we describe a marker-free method for aligning images from light or cryo-light fluorescence microscopy and from electron cryo-microscopy that takes advantage of sample support features, namely the holes in the carbon film. We find that the accuracy of this method, as judged by prediction errors of the hole center coordinates, is better than 100 nm.

RAZA: A Rapid 3D z-crossings algorithm to segment electron tomograms and extract organelles and macromolecules

Resolving the 3D architecture of cells to atomic resolution is one of the most ambitious challenges of cellular and structural biology. Central to this process is the ability to automate tomogram segmentation to identify sub-cellular components, facilitate molecular docking and annotate detected objects with associated metadata. Here we demonstrate that RAZA (Rapid 3D z-crossings algorithm) provides a robust, accurate, intuitive, fast, and generally applicable segmentation algorithm capable of detecting organelles, membranes, macromolecular assemblies and extrinsic membrane protein domains. RAZA defines each continuous contour within a tomogram as a discrete object and extracts a set of 3D structural fingerprints (major, middle and minor axes, surface area and volume), enabling selective, semi-automated segmentation and object extraction. RAZA takes advantage of the fact that the underlying algorithm is a true 3D edge detector, allowing the axes of a detected object to be defined, independent of its random orientation within a cellular tomogram. The selectivity of object segmentation and extraction can be controlled by specifying a user-defined detection tolerance threshold for each fingerprint parameter, within which segmented objects must fall and/or by altering the number of search parameters, to define morphologically similar structures. We demonstrate the capability of RAZA to selectively extract subgroups of organelles (mitochondria) and macromolecular assemblies (ribosomes) from cellular tomograms. Furthermore, the ability of RAZA to define objects and their contours, provides a basis for molecular docking and rapid tomogram annotation.

Quantifying Variability of Manual Annotation in Cryo-Electron Tomograms

Although acknowledged to be variable and subjective, manual annotation of cryo-electron tomography data is commonly used to answer structural questions and to create a "ground truth" for evaluation of automated segmentation algorithms. Validation of such annotation is lacking, but is critical for understanding the reproducibility of manual annotations. Here, we used voxel-based similarity scores for a variety of specimens, ranging in complexity and segmented by several annotators, to quantify the variation among their annotations. In addition, we have identified procedures for merging annotations to reduce variability, thereby increasing the reliability of manual annotation. Based on our analyses, we find that it is necessary to combine multiple manual annotations to increase the confidence level for answering structural questions. We also make recommendations to guide algorithm development for automated annotation of features of interest.

Denoising of high-resolution single-particle electron-microscopy density maps by their approximation using three-dimensional Gaussian functions

Cryo-electron microscopy (cryo-EM) of frozen-hydrated preparations of isolated macromolecular complexes is the method of choice to obtain the structure of complexes that cannot be easily studied by other experimental methods due to their flexibility or large size. An increasing number of macromolecular structures are currently being obtained at subnanometer resolution but the interpretation of structural details in such EM-derived maps is often difficult because of noise at these high-frequency signal components that reduces their contrast. In this paper, we show that the method for EM density-map approximation using Gaussian functions can be used for denoising of single-particle EM maps of high (typically subnanometer) resolution. We show its denoising performance using simulated and experimental EM density maps of several complexes.

Removing Contamination-Induced Reconstruction Artifacts from Cryo-electron Tomograms

Imaging of fully hydrated, vitrified biological samples by electron tomography yields structural information about cellular protein complexes in situ. Here we present a computational procedure that removes artifacts of three-dimensional reconstruction caused by contamination present in samples during imaging by electron microscopy. Applying the procedure to phantom data and electron tomograms of cellular samples significantly improved the resolution and the interpretability of tomograms. Artifacts caused by surface contamination associated with thinning by focused ion beam, as well as those arising from gold fiducial markers and from common, lower contrast contamination, could be removed. Our procedure is widely applicable and is especially suited for applications that strive to reach a higher resolution and involve the use of recently developed, state-of-the-art instrumentation.

Efficient Extraction of Macromolecular Complexes from Electron Tomograms Based on Reduced Representation Templates

Electron tomography is the most widely applicable method for obtaining 3D information by electron microscopy. In the field of biology it has been realized that electron tomography is capable of providing a complete, molecular resolution three-dimensional mapping of entire proteoms. However, to realize this goal, information needs to be extracted efficiently from these tomograms. Owing to extremely low signal-to-noise ratios, this task is mostly carried out manually. Standard template matching approaches tend to generate large amounts of false positives. We developed an alternative method for feature extraction in biological electron tomography based on reduced representation templates, approximating the search model by a small number of anchor points used to calculate the scoring function. Using this approach we see a reduction of about 50% false positives with matched-filter approaches to below 5%. At the same time, false negatives stay below 5%, thus essentially matching the performance one would expect from human operators.

Persistent topology for cryo-EM data analysis

In this work, we introduce persistent homology for the analysis of cryo-electron microscopy (cryo-EM) density maps. We identify the topological fingerprint or topological signature of noise, which is widespread in cryo-EM data. For low signal-to-noise ratio (SNR) volumetric data, intrinsic topological features of biomolecular structures are indistinguishable from noise. To remove noise, we employ geometric flows that are found to preserve the intrinsic topological fingerprints of cryo-EM structures and diminish the topological signature of noise. In particular, persistent homology enables us to visualize the gradual separation of the topological fingerprints of cryo-EM structures from those of noise during the denoising process, which gives rise to a practical procedure for prescribing a noise threshold to extract cryo-EM structure information from noise contaminated data after certain iterations of the geometric flow equation. To further demonstrate the utility of persistent homology for cryo-EM data analysis, we consider a microtubule intermediate structure Electron Microscopy Data (EMD 1129). Three helix models, an alpha-tubulin monomer model, an alpha-tubulin and beta-tubulin model, and an alpha-tubulin and beta-tubulin dimer model, are constructed to fit the cryo-EM data. The least square fitting leads to similarly high correlation coefficients, which indicates that structure determination via optimization is an ill-posed inverse problem. However, these models have dramatically different topological fingerprints. Especially, linkages or connectivities that discriminate one model from another, play little role in the traditional density fitting or optimization but are very sensitive and crucial to topological fingerprints. The intrinsic topological features of the microtubule data are identified after topological denoising. By a comparison of the topological fingerprints of the original data and those of three models, we found that the third model is topologically favored. The present work offers persistent homology based new strategies for topological denoising and for resolving ill-posed inverse problems.

Development and integration of block operations for data invariant automation of digital preprocessing and analysis of biological and biomedical Raman spectra

High-throughput information extraction from large numbers of Raman spectra is becoming an increasingly taxing problem due to the proliferation of new applications enabled using advances in instrumentation. Fortunately, in many of these applications, the entire process can be automated, yielding reproducibly good results with significant time and cost savings. Information extraction consists of two stages, preprocessing and analysis. We focus here on the preprocessing stage, which typically involves several steps, such as calibration, background subtraction, baseline flattening, artifact removal, smoothing, and so on, before the resulting spectra can be further analyzed. Because the results of some of these steps can affect the performance of subsequent ones, attention must be given to the sequencing of steps, the compatibility of these sequences, and the propensity of each step to generate spectral distortions. We outline here important considerations to effect full automation of Raman spectral preprocessing: what is considered full automation; putative general principles to effect full automation; the proper sequencing of processing and analysis steps; conflicts and circularities arising from sequencing; and the need for, and approaches to, preprocessing quality control. These considerations are discussed and illustrated with biological and biomedical examples reflecting both successful and faulty preprocessing.

From voxels to knowledge: a practical guide to the segmentation of complex electron microscopy 3D-data

Modern 3D electron microscopy approaches have recently allowed unprecedented insight into the 3D ultrastructural organization of cells and tissues, enabling the visualization of large macromolecular machines, such as adhesion complexes, as well as higher-order structures, such as the cytoskeleton and cellular organelles in their respective cell and tissue context. Given the inherent complexity of cellular volumes, it is essential to first extract the features of interest in order to allow visualization, quantification, and therefore comprehension of their 3D organization. Each data set is defined by distinct characteristics, e.g., signal-to-noise ratio, crispness (sharpness) of the data, heterogeneity of its features, crowdedness of features, presence or absence of characteristic shapes that allow for easy identification, and the percentage of the entire volume that a specific region of interest occupies. All these characteristics need to be considered when deciding on which approach to take for segmentation.

The six different 3D ultrastructural data sets presented were obtained by three different imaging approaches: resin embedded stained electron tomography, focused ion beam- and serial block face- scanning electron microscopy (FIB-SEM, SBF-SEM) of mildly stained and heavily stained samples, respectively. For these data sets, four different segmentation approaches have been applied: (1) fully manual model building followed solely by visualization of the model, (2) manual tracing segmentation of the data followed by surface rendering, (3) semi-automated approaches followed by surface rendering, or (4) automated custom-designed segmentation algorithms followed by surface rendering and quantitative analysis. Depending on the combination of data set characteristics, it was found that typically one of these four categorical approaches outperforms the others, but depending on the exact sequence of criteria, more than one approach may be successful. Based on these data, we propose a triage scheme that categorizes both objective data set characteristics and subjective personal criteria for the analysis of the different data sets.

Local regularization of tilt projections reduces artifacts in electron tomography

Electron tomography produces very high resolution 3D image volumes useful for investigating the structure and function of cellular components. Unfortunately, unavoidable discontinuities and physical constraints in the acquisition geometry lead to a range of artifacts that can affect the reconstructed image. In particular, highly electron dense regions, such as gold nanoparticles, can hide proximal biological structures and degrade the overall quality of the reconstructed tomograms. In this work we introduce a pre-reconstruction non-conservative non-linear isotropic diffusion (NID) filter that automatically identifies and reduces local irregularities in the tilt projections. We illustrate the improvement in quality obtained using this approach for reconstructed tomograms generated from samples of malaria parasite-infected red blood cells. A quantitative and qualitative evaluation for our approach on both simulated and real data is provided.

Robust membrane detection based on tensor voting for electron tomography

Electron tomography enables three-dimensional (3D) visualization and analysis of the subcellular architecture at a resolution of a few nanometers. Segmentation of structural components present in 3D images (tomograms) is often necessary for their interpretation. However, it is severely hampered by a number of factors that are inherent to electron tomography (e.g. noise, low contrast, distortion). Thus, there is a need for new and improved computational methods to facilitate this challenging task. In this work, we present a new method for membrane segmentation that is based on anisotropic propagation of the local structural information using the tensor voting algorithm. The local structure at each voxel is then refined according to the information received from other voxels. Because voxels belonging to the same membrane have coherent structural information, the underlying global structure is strengthened. In this way, local information is easily integrated at a global scale to yield segmented structures. This method performs well under low signal-to-noise ratio typically found in tomograms of vitrified samples under cryo-tomography conditions and can bridge gaps present on membranes. The performance of the method is demonstrated by applications to tomograms of different biological samples and by quantitative comparison with standard template matching procedure.

Three-dimensional reconstructions of actin filaments capped by Arp2/3 complex

The primary function of Arp2/3 complex is the generation of free barbed ends by nucleating new filaments from the sides of pre-existing filaments. The pathway of branch formation is complex and involves nucleation promoting factors, actin monomers and nucleotides. A less prominent function of Arp2/3 complex is capping of actin filament pointed ends. Here we show, using electron microscopy, electron tomography, and image reconstruction of negatively-stained samples at ∼2-3nm resolution, that Arp2/3 complex bound to the pointed ends of actin filaments has a conformation similar to that in the branch junction with the Arps arranged in an actin-filament like configuration. This is direct evidence for the existence of two distinct activation pathways for Arp2/3 complex, one in the context of branch formation, one in the context of pointed-end capping, with essentially the same conformational end point.

Multiscale geometric modeling of macromolecules I: Cartesian representation

This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the polarized curvature, for the prediction of protein binding sites.

Microscopy basics and the study of actin-actin-binding protein interactions

Actin is a multifunctional eukaryotic protein with a globular monomer form that polymerizes into a thin, linear microfilament in cells. Through interactions with various actin-binding proteins (ABPs), actin plays an active role in many cellular processes, such as cell motility and structure. Microscopy techniques are powerful tools for determining the role and mechanism of actin-ABP interactions in these processes. In this article, we describe the basic concepts of fluorescent speckle microscopy, total internal reflection fluorescence microscopy, atomic force microscopy, and cryoelectron microscopy and review recent studies that utilize these techniques to visualize the binding of actin with ABPs.

The architectural relationship of components controlling mast cell endocytosis

Eukaryotic cells use multiple routes for receptor internalization. Here, we examine the topographical relationships of clathrin-dependent and clathrin-independent endocytic structures on the plasma membranes of leukemia-derived mast cells. The high affinity IgE receptor (FcεRI) utilizes both pathways, whereas transferrin receptor serves as a marker for the classical clathrin-mediated endocytosis pathway. Both receptors were tracked by live-cell imaging in the presence or absence of inhibitors that established their differential dependence on specific endocytic adaptor proteins. The topology of antigen-bound FcεRI, clathrin, dynamin, Arf6 and Eps15-positive structures were analyzed by 2D and 3D immunoelectron microscopy techniques, revealing their remarkable spatial relationships and unique geometry. We conclude that the mast cell plasma membrane has multiple specialized domains for endocytosis. Their close proximity might reflect shared components, such as lipids and adaptor proteins, that facilitate inward membrane curvature. Intersections between these specialized domains might represent sorting stations that direct cargo to specific endocytic pathways.

Geometric modeling of subcellular structures, organelles, and multiprotein complexes

Recently, the structure, function, stability, and dynamics of subcellular structures, organelles, and multiprotein complexes have emerged as a leading interest in structural biology. Geometric modeling not only provides visualizations of shapes for large biomolecular complexes but also fills the gap between structural information and theoretical modeling, and enables the understanding of function, stability, and dynamics. This paper introduces a suite of computational tools for volumetric data processing, information extraction, surface mesh rendering, geometric measurement, and curvature estimation of biomolecular complexes. Particular emphasis is given to the modeling of cryo-electron microscopy data. Lagrangian-triangle meshes are employed for the surface presentation. On the basis of this representation, algorithms are developed for surface area and surface-enclosed volume calculation, and curvature estimation. Methods for volumetric meshing have also been presented. Because the technological development in computer science and mathematics has led to multiple choices at each stage of the geometric modeling, we discuss the rationales in the design and selection of various algorithms. Analytical models are designed to test the computational accuracy and convergence of proposed algorithms. Finally, we select a set of six cryo-electron microscopy data representing typical subcellular complexes to demonstrate the efficacy of the proposed algorithms in handling biomolecular surfaces and explore their capability of geometric characterization of binding targets. This paper offers a comprehensive protocol for the geometric modeling of subcellular structures, organelles, and multiprotein complexes.

Developing a denoising filter for electron microscopy and tomography data in the cloud

The low radiation conditions and the predominantly phase-object image formation of cryo-electron microscopy (cryo-EM) result in extremely high noise levels and low contrast in the recorded micrographs. The process of single particle or tomographic 3D reconstruction does not completely eliminate this noise and is even capable of introducing new sources of noise during alignment or when correcting for instrument parameters. The recently developed Digital Paths Supervised Variance (DPSV) denoising filter uses local variance information to control regional noise in a robust and adaptive manner. The performance of the DPSV filter was evaluated in this review qualitatively and quantitatively using simulated and experimental data from cryo-EM and tomography in two and three dimensions. We also assessed the benefit of filtering experimental reconstructions for visualization purposes and for enhancing the accuracy of feature detection. The DPSV filter eliminates high-frequency noise artifacts (density gaps), which would normally preclude the accurate segmentation of tomography reconstructions or the detection of alpha-helices in single-particle reconstructions. This collaborative software development project was carried out entirely by virtual interactions among the authors using publicly available development and file sharing tools.

A 3D image filter for parameter-free segmentation of macromolecular structures from electron tomograms

3D image reconstruction of large cellular volumes by electron tomography (ET) at high (≤ 5 nm) resolution can now routinely resolve organellar and compartmental membrane structures, protein coats, cytoskeletal filaments, and macromolecules. However, current image analysis methods for identifying in situ macromolecular structures within the crowded 3D ultrastructural landscape of a cell remain labor-intensive, time-consuming, and prone to user-bias and/or error. This paper demonstrates the development and application of a parameter-free, 3D implementation of the bilateral edge-detection (BLE) algorithm for the rapid and accurate segmentation of cellular tomograms. The performance of the 3D BLE filter has been tested on a range of synthetic and real biological data sets and validated against current leading filters-the pseudo 3D recursive and Canny filters. The performance of the 3D BLE filter was found to be comparable to or better than that of both the 3D recursive and Canny filters while offering the significant advantage that it requires no parameter input or optimisation. Edge widths as little as 2 pixels are reproducibly detected with signal intensity and grey scale values as low as 0.72% above the mean of the background noise. The 3D BLE thus provides an efficient method for the automated segmentation of complex cellular structures across multiple scales for further downstream processing, such as cellular annotation and sub-tomogram averaging, and provides a valuable tool for the accurate and high-throughput identification and annotation of 3D structural complexity at the subcellular level, as well as for mapping the spatial and temporal rearrangement of macromolecular assemblies in situ within cellular tomograms.

Electron tomography of cells

The electron microscope has contributed deep insights into biological structure since its invention nearly 80 years ago. Advances in instrumentation and methodology in recent decades have now enabled electron tomography to become the highest resolution three-dimensional (3D) imaging technique available for unique objects such as cells. Cells can be imaged either plastic-embedded or frozen-hydrated. Then the series of projection images are aligned and back-projected to generate a 3D reconstruction or 'tomogram'. Here, we review how electron tomography has begun to reveal the molecular organization of cells and how the existing and upcoming technologies promise even greater insights into structural cell biology.

A supramolecular assembly formed by influenza A virus genomic RNA segments

The influenza A virus genome consists of eight viral RNAs (vRNAs) that form viral ribonucleoproteins (vRNPs). Even though evidence supporting segment-specific packaging of vRNAs is accumulating, the mechanism ensuring selective packaging of one copy of each vRNA into the viral particles remains largely unknown. We used electron tomography to show that the eight vRNPs emerge from a common ‘transition zone’ located underneath the matrix layer at the budding tip of the virions, where they appear to be interconnected and often form a star-like structure. This zone appears as a platform in 3D surface rendering and is thick enough to contain all known packaging signals. In vitro, all vRNA segments are involved in a single network of intermolecular interactions. The regions involved in the strongest interactions were identified and correspond to known packaging signals. A limited set of nucleotides in the 5′ region of vRNA 7 was shown to interact with vRNA 6 and to be crucial for packaging of the former vRNA. Collectively, our findings support a model in which the eight genomic RNA segments are selected and packaged as an organized supramolecular complex held together by direct base pairing of the packaging signals.

Automated spectral smoothing with spatially adaptive penalized least squares

A variety of data smoothing techniques exist to address the issue of noise in spectroscopic data. The vast majority, however, require parameter specification by a knowledgeable user, which is typically accomplished by trial and error. In most situations, optimized parameters represent a compromise between noise reduction and signal preservation. In this work, we demonstrate a nonparametric regression approach to spectral smoothing using a spatially adaptive penalized least squares (SAPLS) approach. An iterative optimization procedure is employed that permits gradual flexibility in the smooth fit when statistically significant trends based on multiscale statistics assuming white Gaussian noise are detected. With an estimate of the noise level in the spectrum the procedure is fully automatic with a specified confidence level for the statistics. Potential application to the heteroscedastic noise case is also demonstrated. Performance was assessed in simulations conducted on several synthetic spectra using traditional error measures as well as comparisons of local extrema in the resulting smoothed signals to those in the true spectra. For the simulated spectra, a best case comparison with the Savitzky-Golay smoothing via an exhaustive parameter search was performed while the SAPLS method was assessed for automated application. The application to several dissimilar experimentally obtained Raman spectra is also presented.

Three-dimensional structure of Rubella virus factories

Viral factories are complex structures in the infected cell where viruses compartmentalize their life cycle. Rubella virus (RUBV) assembles factories by recruitment of rough endoplasmic reticulum (RER), mitochondria and Golgi around modified lysosomes known as cytopathic vacuoles or CPVs. These organelles contain active replication complexes that transfer replicated RNA to assembly sites in Golgi membranes. We have studied the structure of RUBV factory in three dimensions by electron tomography and freeze-fracture. CPVs contain stacked membranes, rigid sheets, small vesicles and large vacuoles. These membranes are interconnected and in communication with the endocytic pathway since they incorporate endocytosed BSA-gold. RER and CPVs are coupled through protein bridges and closely apposed membranes. Golgi vesicles attach to the CPVs but no tight contacts with mitochondria were detected. Immunogold labelling confirmed that the mitochondrial protein p32 is an abundant component around and inside CPVs where it could play important roles in factory activities.

Formation and reduction of streak artefacts in electron tomography

We have analysed the formation of streak artefacts in the reconstruction based on the filtered back projection algorithm in electron tomography (ET) and accordingly applied an adaptive interpolation technique to artefact reduction. In the adaptive interpolation to recover the missing information, the edge positions in a projection curve were tracked to reduce the interpolation error. A simulation was used to demonstrate the effectiveness of the artefact reduction. Furthermore, image reconstruction of integrated circuit specimens in the ET experiments with the ultra-high voltage electron microscope show that the strong streak artefacts can be reduced effectively by our artefact reduction technique.

Methods for segmentation and interpretation of electron tomographic reconstructions

Electron tomography has become a powerful tool for revealing the molecular architecture of biological cells and tissues. In principle, electron tomography can provide high-resolution mapping of entire proteomes. The achievable resolution (3-8 nm) is capable of bridging the gap between live-cell imaging and atomic resolution structures. However, the relevant information is not readily accessible from the data and needs to be identified, extracted, and processed before it can be used. Because electron tomography imaging and image acquisition technologies have enjoyed major advances in the last few years and continue to increase data throughput, the need for approaches that allow automatic and objective interpretation of electron tomograms becomes more and more urgent. This chapter provides an overview of the state of the art in this field and attempts to identify the major bottlenecks that prevent approaches for interpreting electron tomography data to develop their full potential.

Confidence intervals for fitting of atomic models into low-resolution densities

The fitting of high-resolution structures into low-resolution densities obtained from techniques such as electron microscopy or small-angle X-ray scattering can yield powerful new insights. While several algorithms for achieving optimal fits have recently been developed, relatively little effort has been devoted to developing objective measures for judging the quality of the resulting fits, in particular with regard to the danger of overfitting. Here, a general method is presented for obtaining confidence intervals for atomic coordinates resulting from fitting of atomic resolution domain structures into low-resolution densities using well established statistical tools. It is demonstrated that the resulting confidence intervals are sufficiently accurate to allow meaningful statistical tests and to provide tools for detecting potential overfitting.

TOMOBFLOW: feature-preserving noise filtering for electron tomography

Background

Noise filtering techniques are needed in electron tomography to allow proper interpretation of datasets. The standard linear filtering techniques are characterized by a tradeoff between the amount of reduced noise and the blurring of the features of interest. On the other hand, sophisticated anisotropic nonlinear filtering techniques allow noise reduction with good preservation of structures. However, these techniques are computationally intensive and are difficult to be tuned to the problem at hand.

Results

TOMOBFLOW is a program for noise filtering with capabilities of preservation of biologically relevant information. It is an efficient implementation of the Beltrami flow, a nonlinear filtering method that locally tunes the strength of the smoothing according to an edge indicator based on geometry properties. The fact that this method does not have free parameters hard to be tuned makes TOMOBFLOW a user-friendly filtering program equipped with the power of diffusion-based filtering methods. Furthermore, TOMOBFLOW is provided with abilities to deal with different types and formats of images in order to make it useful for electron tomography in particular and bioimaging in general.

Conclusion

TOMOBFLOW allows efficient noise filtering of bioimaging datasets with preservation of the features of interest, thereby yielding data better suited for post-processing, visualization and interpretation. It is available at the web site .

The unique architecture of Bunyamwera virus factories around the Golgi complex

Viral factories are novel structures built by viruses in infected cells. During their construction organelles are recruited and build a large scaffold for viral replication and morphogenesis. We have studied how a bunyavirus uses the Golgi to build the factory. With the help of confocal and 3D ultrastructural imaging together with molecular mapping in situ and in vitro we have characterized a tubular structure that harbours the viral replication complexes in a globular domain. Numerous ribonucleoproteins were released from purified tubes disrupted in vitro. Actin and myosin I were identified by peptide mass fingerprinting in isolated tubes while actin and the viral NSm non-structural protein were detected in the tubes' internal proteinaceous scaffold by immunogold labelling. Studies with NSm deletion mutants and drugs affecting actin showed that both NSm and actin are key factors for tube and virus assembly in Golgi. Three-dimensional reconstructions based on oriented serial sections of infected cells showed that tubes anchor cell organelles to Golgi stacks and make contacts with intracellular viruses. We propose that this new structure, unique among enveloped viruses, assembles in association with the most stable component of Golgi stacks, the actin-containing matrix scaffold, connecting viral replication and morphogenesis inside viral factories.

Fully automated high-performance signal-to-noise ratio enhancement based on an iterative three-point zero-order Savitzky-Golay filter

The automated processing of data from high-throughput and real-time collection procedures is becoming a pressing problem. Currently the focus is shifting to automated smoothing techniques where, unlike background subtraction techniques, very few methods exist. We have developed a filter based on the widely used and conceptually simple moving average method or zero-order Savitzky-Golay filter and its iterative relative, the Kolmogorov-Zurbenko filter. A crucial difference, however, between these filters and our implementation is that our fully automated smoothing filter requires no parameter specification or parameter optimization. Results are comparable to, or better than, Savitzky-Golay filters with optimized parameters and superior to the automated iterative median filter. Our approach, because it is based on the highly familiar moving average concept, is intuitive, fast, and straightforward to implement and should therefore be of immediate and considerable practical use in a wide variety of spectroscopy applications.

The structural basis of actin filament branching by the Arp2/3 complex

The actin-related protein 2/3 (Arp2/3) complex mediates the formation of branched actin filaments at the leading edge of motile cells and in the comet tails moving certain intracellular pathogens. Crystal structures of the Arp2/3 complex are available, but the architecture of the junction formed by the Arp2/3 complex at the base of the branch was not known. In this study, we use electron tomography to reconstruct the branch junction with sufficient resolution to show how the Arp2/3 complex interacts with the mother filament. Our analysis reveals conformational changes in both the mother filament and Arp2/3 complex upon branch formation. The Arp2 and Arp3 subunits reorganize into a dimer, providing a short-pitch template for elongation of the daughter filament. Two subunits of the mother filament undergo conformational changes that increase stability of the branch. These data provide a rationale for why branch formation requires cooperative interactions among the Arp2/3 complex, nucleation-promoting factors, an actin monomer, and the mother filament.

Computational resources for cryo-electron tomography in Bsoft

The Bsoft package [Heymann, J.B., Belnap, D.M., 2007. Bsoft: image processing and molecular modeling for electron microscopy. J. Struct. Biol. 157, 3-18] has been enhanced by adding utilities for processing electron tomographic (ET) data; in particular, cryo-ET data characterized by low contrast and high noise. To handle the high computational load efficiently, a workflow was developed, based on the database-like parameter handling in Bsoft, aimed at minimizing user interaction and facilitating automation. To the same end, scripting elements distribute the processing among multiple processors on the same or different computers. The resolution of a tomogram depends on the precision of projection alignment, which is usually based on pinpointing fiducial markers (electron-dense gold particles). Alignment requires accurate specification of the tilt axis, and our protocol includes a procedure for determining it to adequate accuracy. Refinement of projection alignment provides information that allows assessment of its precision, as well as projection quality control. We implemented a reciprocal space algorithm that affords an alternative to back-projection or real space algorithms for calculating tomograms. Resources are also included that allow resolution assessment by cross-validation (NLOO2D); denoising and interpretation; and the extraction, mutual alignment, and averaging of tomographic sub-volumes.

Bilateral edge filter: photometrically weighted, discontinuity based edge detection

Edge-detection algorithms have the potential to play an increasingly important role both in single particle analysis (for the detection of randomly oriented particles), and in tomography (for the segmentation of 3D volumes). However, the majority of traditional linear filters are significantly affected by noise as well as artefacts, and offer limited selectivity. The Bilateral edge filter presented here is an adaptation of the Bilateral filter [Jiang, W., Baker, M.L., Wu, Q., Bajaj, C., Chiu, W., 2003. Applications of a bilateral denoising filter in biological electron microscopy. J. Struct. Biol. 144, 114-122] designed for enhanced edge detection. It uses photometric weighting to identify significant discontinuities (representing edges), minimizing artefacts and noise. Compared with common edge-detectors (LoG, Marr-Hildreth) the Bilateral edge filter yielded significantly better results. Indeed data was of a similar quality to that of the Canny edge-detector, which is considered as a leading standard in edge detection [Basu, M., 2002. Gaussian-based edge-detection methods-a survey. IEEE Trans. Syst. Man Cybern. C Appl. Rev. 32, 252-260]. Compared to the Canny edge-detector the Bilateral edge-detector has the advantages that it only requires the adjustment of a single parameter, is theoretically faster for reasonably sized images, and can be used in selective contrast enhancement of images. The simplicity and speed of the filter for single particle and tomographic analysis are discussed.