MOTIVATION: Homology search for RNAs can use secondary structure information to increase power by modeling base pairs, as in covariance models, but the resulting computational costs are high. Typical acceleration strategies rely on at least one filtering stage using sequence-only search. RESULTS: Here we present the multi-segment CYK (MSCYK) filter, which implements a heuristic of ungapped structural alignment for RNA homology search. Compared to gapped alignment, this approximation has lower computation time requirements (O(N⁴) reduced to O(N³), and space requirements (O(N³) reduced to O(N²). A vector-parallel implementation of this method gives up to 100-fold speed-up; vector-parallel implementations of standard gapped alignment at two levels of precision give 3- and 6-fold speed-ups. These approaches are combined to create a filtering pipeline that scores RNA secondary structure at all stages, with results that are synergistic with existing methods.

Live fluorescence microscopy has the unique capability to probe dynamic processes, linking molecular components and their localization with function. A key goal of microscopy is to increase spatial and temporal resolution while simultaneously permitting identification of multiple specific components. We demonstrate a new microscope platform, OMX, that enables subsecond, multicolor four-dimensional data acquisition and also provides access to subdiffraction structured illumination imaging. Using this platform to image chromosome movement during a complete yeast cell cycle at one 3D image stack per second reveals an unexpected degree of photosensitivity of fluorophore-containing cells. To avoid perturbation of cell division, excitation levels had to be attenuated between 100 and 10,000× below the level normally used for imaging. We show that an image denoising algorithm that exploits redundancy in the image sequence over space and time allows recovery of biological information from the low light level noisy images while maintaining full cell viability with no fading.

Many important physiological processes operate at time and space scales far beyond those accessible to atom-realistic simulations, and yet discrete stochastic rather than continuum methods may best represent finite numbers of molecules interacting in complex cellular spaces. We describe and validate new tools and algorithms developed for a new version of the MCell simulation program (MCell3), which supports generalized Monte Carlo modeling of diffusion and chemical reaction in solution, on surfaces representing membranes, and combinations thereof. A new syntax for describing the spatial directionality of surface reactions is introduced, along with optimizations and algorithms that can substantially reduce computational costs (e.g., event scheduling, variable time and space steps). Examples for simple reactions in simple spaces are validated by comparison to analytic solutions. Thus we show how spatially realistic Monte Carlo simulations of biological systems can be far more cost-effective than often is assumed, and provide a level of accuracy and insight beyond that of continuum methods.

Conventional acquisition of three-dimensional (3D) microscopy data requires sequential z scanning and is often too slow to capture biological events. We report an aberration-corrected multifocus microscopy method capable of producing an instant focal stack of nine 2D images. Appended to an epifluorescence microscope, the multifocus system enables high-resolution 3D imaging in multiple colors with single-molecule sensitivity, at speeds limited by the camera readout time of a single image.

In order to study anatomy of organisms with high-resolution there is an increasing demand to image large specimen in three dimensions (3D). Confocal microscopy is able to produce high-resolution 3D images, but these are limited by its relatively small field of view compared to the size of large biological specimens. To overcome this drawback, motorized stages moving the sample are used to create a tiled scan of the whole specimen. The physical coordinates provided by the microscope stage are not precise enough to allow reconstruction (”Stitching”) of the whole image from individual image stacks.
We developed an algorithm, as well as an ImageJ plug-in, based on the Fourier Shift Theorem that computes all possible translations (x, y, z) between two 3D images at once, yielding the best overlap in terms of the cross correlation measure. Apart from the obvious gain in computation time it has the advantage that it cannot be trapped in local minima as it simply computes all possible solutions. Computing the overlap between two adjacent image stacks is fast (12 seconds for two 512x512x89 images on a Intel ® Core2Duo with 2.2GHz) making it suitable for real time use, i.e. computing the output image during acquisition of the individual image stacks.
To compensate the possible shading- and brightness differences we apply a smooth linear intensity transition between the overlapping stacks. Additionally we extended the to generic 3D registration using gradient based rotation detection on top of the phase correlation method. We demonstrate the performance of our 3D stitching plug-in on several tiled confocal images and show an example of its application for 3D registration.

Gap junctions (GJs) represent connexin-rich membrane domains that connect interiors of adjoining cells in mammalian tissues. How fast GJs can respond to bacterial pathogens has not been known previously. Using Bessel beam plane illumination and confocal spinning disk microscopy, we found fast ( 500 ms) formation of connexin-depleted regions (CDRs) inside GJ plaques between cells exposed to AB5 toxins. CDR formation appears as a fast redistribution of connexin channels within GJ plaques with minor changes in outline or geometry. CDR formation does not depend on membrane trafficking or submembrane cytoskeleton and has no effect on GJ conductance. However, CDR responses depend on membrane lipids, can be modified by cholesterol-clustering agents and extracellular K(+) ion concentration, and influence cAMP signaling. The CDR response of GJ plaques to bacterial toxins is a phenomenon observed for all tested connexin isoforms. Through signaling, the CDR response may enable cells to sense exposure to AB5 toxins. CDR formation may reflect lipid-phase separation events in the biological membrane of the GJ plaque, leading to increased connexin packing and lipid reorganization. Our data demonstrate very fast dynamics (in the millisecond-to-second range) within GJ plaques, which previously were considered to be relatively stable, long-lived structures.

Gap junctions (GJs) represent connexin-rich membrane domains that connect interiors of adjoining cells in mammalian tissues. How fast GJs can respond to bacterial pathogens has not been known previously. Using Bessel beam plane illumination and confocal spinning disk microscopy, we found fast (~500 ms) formation of connexin-depleted regions (CDRs) inside GJ plaques between cells exposed to AB5 toxins. CDR formation appears as a fast redistribution of connexin channels within GJ plaques with minor changes in outline or geometry. CDR formation does not depend on membrane trafficking or submembrane cytoskeleton and has no effect on GJ conductance. However, CDR responses depend on membrane lipids, can be modified by cholesterol-clustering agents and extracellular K(+) ion concentration, and influence cAMP signaling. The CDR response of GJ plaques to bacterial toxins is a phenomenon observed for all tested connexin isoforms. Through signaling, the CDR response may enable cells to sense exposure to AB5 toxins. CDR formation may reflect lipid-phase separation events in the biological membrane of the GJ plaque, leading to increased connexin packing and lipid reorganization. Our data demonstrate very fast dynamics (in the millisecond-to-second range) within GJ plaques, which previously were considered to be relatively stable, long-lived structures.

Cortical information processing is under state-dependent control of subcortical neuromodulatory systems. Although this modulatory effect is thought to be mediated mainly by slow nonsynaptic metabotropic receptors, other mechanisms, such as direct synaptic transmission, are possible. Yet, it is currently unknown if any such form of subcortical control exists. Here, we present direct evidence of a strong, spatiotemporally precise excitatory input from an ascending neuromodulatory center. Selective stimulation of serotonergic median raphe neurons produced a rapid activation of hippocampal interneurons. At the network level, this subcortical drive was manifested as a pattern of effective disynaptic GABAergic inhibition that spread throughout the circuit. This form of subcortical network regulation should be incorporated into current concepts of normal and pathological cortical function.

The comprehensive reconstruction of cell lineages in complex multicellular organisms is a central goal of developmental biology. We present an open-source computational framework for the segmentation and tracking of cell nuclei with high accuracy and speed. We demonstrate its (i) generality by reconstructing cell lineages in four-dimensional, terabyte-sized image data sets of fruit fly, zebrafish and mouse embryos acquired with three types of fluorescence microscopes, (ii) scalability by analyzing advanced stages of development with up to 20,000 cells per time point at 26,000 cells min(-1) on a single computer workstation and (iii) ease of use by adjusting only two parameters across all data sets and providing visualization and editing tools for efficient data curation. Our approach achieves on average 97.0% linkage accuracy across all species and imaging modalities. Using our system, we performed the first cell lineage reconstruction of early Drosophila melanogaster nervous system development, revealing neuroblast dynamics throughout an entire embryo.