Energetic considerations of ciliary beating and the advantage of metachronal coordination

The internal mechanism of cilia is among the most ancient biological motors on an evolutionary scale. It produces beat patterns that consist of two phases: during the effective stroke, the cilium moves approximately as a straight rod, and during the recovery stroke, it rolls close to the surface in a tangential motion. It is commonly agreed that these two phases are designed for efficient functioning: the effective stroke encounters strong viscous resistance and generates thrust, whereas the recovery stroke returns the cilium to starting position while avoiding viscous resistance. Metachronal coordination between cilia, which occurs when many of them beat close to each other, is believed to be an autonomous result of the hydrodynamical interactions in the system. Qualitatively, metachronism is perceived as a way for reducing the energy expenditure required for beating. This paper presents a quantitative study of the energy expenditure of beating cilia, and of the energetic significance of metachronism. We develop a method for computing the work done by model cilia that beat in a viscous fluid. We demonstrate that for a single cilium, beating in water, the mechanical work done during the effective stroke is approximately five times the amount of work done during the recovery stroke. Investigation of multicilia configurations shows that having neighboring cilia beat metachronally is energetically advantageous and perhaps even crucial for multiciliary functioning. Finally, the model is used to approximate the number of dynein arm attachments that are likely to occur during the effective and recovery strokes of a beat cycle, predicting that almost all of the available dynein arms should participate in generating the motion.

ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

Complex mtDNA constitutes an approximate 620-kb insertion on Arabidopsis thaliana chromosome 2: Implication of potential sequencing errors caused by large-unit repeats

Previously conducted sequence analysis of Arabidopsis thaliana (ecotype Columbia-0) reported an insertion of 270-kb mtDNA into the pericentric region on the short arm of chromosome 2. DNA fiber-based fluorescence in situ hybridization analyses reveal that the mtDNA insert is 618 ± 42 kb, ≈2.3 times greater than that determined by contig assembly and sequencing analysis. Portions of the mitochondrial genome previously believed to be absent were identified within the insert. Sections of the mtDNA are repeated throughout the insert. The cytological data illustrate that DNA contig assembly by using bacterial artificial chromosomes tends to produce a minimal clone path by skipping over duplicated regions, thereby resulting in sequencing errors. We demonstrate that fiber-fluorescence in situ hybridization is a powerful technique to analyze large repetitive regions in the higher eukaryotic genomes and is a valuable complement to ongoing large genome sequencing projects.