Discrete cilia modelling with singularity distributions

We discuss in detail techniques for modelling flows due to finite and infinite arrays of beating cilia. An efficient technique, based on concepts from previous ‘singularity models’ is described, that is accurate in both near and far-fields. Cilia are modelled as curved slender ellipsoidal bodies by distributing Stokeslet and potential source dipole singularities along their centrelines, leading to an integral equation that can be solved using a simple and efficient discretisation. The computed velocity on the cilium surface is found to compare favourably with the boundary condition. We then present results for two topics of current interest in biology. 1) We present the first theoretical results showing the mechanism by which rotating embryonic nodal cilia produce a leftward flow by a ‘posterior tilt,’ and track particle motion in an array of three simulated nodal cilia. We find that, contrary to recent suggestions, there is no continuous layer of negative fluid transport close to the ciliated boundary. The mean leftward particle transport is found to be just over 1 μm/s, within experimentally measured ranges. We also discuss the accuracy of models that represent the action of cilia by steady rotlet arrays, in particular, confirming the importance of image systems in the boundary in establishing the far-field fluid transport. Future modelling may lead to understanding of the mechanisms by which morphogen gradients or mechanosensing cilia convert a directional flow to asymmetric gene expression. 2) We develop a more complex and detailed model of flow patterns in the periciliary layer of the airway surface liquid. Our results confirm that shear flow of the mucous layer drives a significant volume of periciliary liquid in the direction of mucus transport even during the recovery stroke of the cilia. Finally, we discuss the advantages and disadvantages of the singularity technique and outline future theoretical and experimental developments required to apply this technique to various other biological problems, particularly in the reproductive system

On Birthing Dancing Stars: The Need for Bounded Chaos in Information Interaction

While computers causing chaos is acommon social trope, nearly the entirety of the history of computing is dedicated to generating order. Typical interactive information retrieval tasks ask computers to support the traversal and exploration of large, complex information spaces. The implicit assumption is that they are to support users in simplifying the complexity (i.e. in creating order from chaos). But for some types of task, particularly those that involve the creative application or synthesis of knowledge or the creation of new knowledge, this assumption may be incorrect. It is increasingly evident that perfect order—and the systems we create with it—support highly-structured information tasks well, but provide poor support for less-structured tasks.We need digital information environments that help create a little more chaos from order to spark creative thinking and knowledge creation. This paper argues for the need for information systems that offerwhat we term ‘bounded chaos’, and offers research directions that may support the creation of such interface

The long-term fate of permafrost peatlands under rapid climate warming

Permafrost peatlands contain globally important amounts of soil organic carbon, owing to cold conditions which suppress anaerobic decomposition. However, climate warming and permafrost thaw threaten the stability of this carbon store. The ultimate fate of permafrost peatlands and their carbon stores is unclear because of complex feedbacks between peat accumulation, hydrology and vegetation. Field monitoring campaigns only span the last few decades and therefore provide an incomplete picture of permafrost peatland response to recent rapid warming. Here we use a high-resolution palaeoecological approach to understand the longer-term response of peatlands in contrasting states of permafrost degradation to recent rapid warming. At all sites we identify a drying trend until the late-twentieth century; however, two sites subsequently experienced a rapid shift to wetter conditions as permafrost thawed in response to climatic warming, culminating in collapse of the peat domes. Commonalities between study sites lead us to propose a five-phase model for permafrost peatland response to climatic warming. This model suggests a shared ecohydrological trajectory towards a common end point: inundated Arctic fen. Although carbon accumulation is rapid in such sites, saturated soil conditions are likely to cause elevated methane emissions that have implications for climate-feedback mechanisms

Genetic relationships within and among Iberian fescues (Festuca L.) based on PCR-amplified markers

The genus Festuca comprises approximately 450 species and is widely distributed around the world. The Iberian Penninsula, with more than 100 taxa colonizing very diverse habitats, is one of its main centers of diversification. This study was conducted to assess molecular genetic variation and genetic relatedness among 91 populations of 31 taxa of Iberian fescues, based on several molecular markers (random amplified polymorphic DNA, amplified fragment length polymorphisms, and trnL sequences). The analyses showed the paraphyletic origin of the broad-leaved (subgenus Festuca, sections Scariosae and Subbulbosae, and subgenus Schedonorus) and the fine-leaved fescues (subgenus Festuca, sections Aulaxyper, Eskia, and Festuca). Schedonorus showed a weak relationship with Lolium rigidum and appeared to be the most recent of the broad-leaved clade. Section Eskia was the most ancient and Festuca the most recent of the fine-leaved clade. Festuca and Aulaxyper were the most related sections, in concordance with their taxonomic affinities. All taxa grouped into their sections, except F. ampla and F. capillifolia (section Festuca), which appeared to be more closely related to Aulaxyper and to a new independent section, respectively. Most populations clustered at the species level, but some subspecies and varieties mixed their populations. This study demonstrated the value in combining different molecular markers to uncover hidden genetic relationships between populations of Festuca

The evolution of galaxy groups and of galaxies therein

Properties of groups of galaxies depend sensitively on the algorithm for group selection, and even the most recent catalogs of groups built from redshift-space selection should suffer from projections and infalling galaxies. The cosmo-dynamical evolution of groups from initial Hubble expansion to collapse and virialization leads to a fundamental track (FT) in virial-theorem-M/L vs crossing time. The increased rates of mergers, both direct and after dynamical friction, in groups relative to clusters, explain the higher fraction of elliptical galaxies at given local number density in X-ray selected groups, relative to clusters, even when the hierarchical evolution of groups is considered. Galaxies falling into groups and clusters should later travel outwards to typically 2 virial radii, which is somewhat less than the outermost radius where observed galaxy star formation efficiencies are enhanced relative to field galaxies of same morphological type. An ongoing analysis of the internal kinematics of X-ray selected groups suggests that the radial profiles of line of sight velocity dispersion are consistent with isotropic NFW distributions for the total mass density, with higher (lower) concentrations than LambdaCDM predictions in groups of high (low) mass. The critical mass, at M200 ~ 10^13 M_sun is consistent with possible breaks in the X-ray luminosity-temperature and Fundamental Plane relations. The internal kinematics of groups indicate that the M-T relation of groups should agree with that extrapolated from clusters with no break at the group scale. The analyses of observed velocity dispersion profiles and of the FT both suggest that low velocity dispersion groups (compact and loose, X-ray emitting or undetected) are quite contaminated by chance projections.Comment: Invited review, ESO workshop "Groups of Galaxies in the Nearby Universe", held in Santiago, Chile, 5-9 December 2005, ed. I. Saviane, V. Ivanov & J. Borissova, 16 page

Nucleation of a sodium droplet on C60

We investigate theoretically the progressive coating of C60 by several sodium atoms. Density functional calculations using a nonlocal functional are performed for NaC60 and Na2C60 in various configurations. These data are used to construct an empirical atomistic model in order to treat larger sizes in a statistical and dynamical context. Fluctuating charges are incorporated to account for charge transfer between sodium and carbon atoms. By performing systematic global optimization in the size range 1<=n<=30, we find that Na_nC60 is homogeneously coated at small sizes, and that a growing droplet is formed above n=>8. The separate effects of single ionization and thermalization are also considered, as well as the changes due to a strong external electric field. The present results are discussed in the light of various experimental data.Comment: 17 pages, 10 figure

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa.

Rapidly growing biological data-including molecular sequences and fossils-hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a "Dated Tree of Life" where all node ages are directly comparable. [Bayesian phylogenetics; data mining; divide-and-conquer methods; GenBank; multilocus multispecies coalescent; next-generation sequencing; palms; primates; tree calibration.]

Validating the RedMIT/GFP-LC3 Mouse Model by Studying Mitophagy in Autosomal Dominant Optic Atrophy Due to the OPA1Q285STOP Mutation

Background: Autosomal dominant optic atrophy (ADOA) is usually caused by mutations in the essential gene, OPA1. This encodes a ubiquitous protein involved in mitochondrial dynamics, hence tissue specificity is not understood. Dysregulated mitophagy (mitochondria recycling) is implicated in ADOA, being increased in OPA1 patient fibroblasts. Furthermore, autophagy may be increased in retinal ganglion cells (RGCs) of the OPA1Q285STOPmouse model. Aims: We developed a mouse model for studying mitochondrial dynamics in order to investigate mitophagy in ADOA. Methods: We crossed the OPA1Q285STOPmouse with our RedMIT/GFP-LC3 mouse, harboring red fluorescent mitochondria and green fluorescent autophagosomes. Colocalization between mitochondria and autophagosomes, the hallmark of mitophagy, was quantified in fluorescently labeled organelles in primary cell cultures, using two high throughput imaging methods Imagestream (Amnis) and IN Cell Analyzer 1000 (GE Healthcare Life Sciences). We studied colocalization between mitochondria and autophagosomes in fixed sections using confocal microscopy. Results: We validated our imaging methods for RedMIT/GFP-LC3 mouse cells, showing that colocalization of red fluorescent mitochondria and green fluorescent autophagosomes is a useful indicator of mitophagy. We showed that colocalization increases when lysosomal processing is impaired. Further, colocalization of mitochondrial fragments and autophagosomes is increased in cultures from the OPA1Q285STOP/RedMIT/GFP-LC3 mice compared to RedMIT/GFP-LC3 control mouse cells that were wild type for OPA1. This was apparent in both mouse embryonic fibroblasts (MEFs) using IN Cell 1000 and in splenocytes using ImageStream imaging flow cytometer (Amnis). We confirmed that this represents increased mitophagic flux using lysosomal inhibitors. We also used microscopy to investigate the level of mitophagy in the retina from the OPA1Q285STOP/RedMIT/GFP-LC3 mice and the RedMIT/GFP-LC3 control mice. However, the expression levels of fluorescent proteins and the image signal-to-background ratios precluded the detection of colocalization so we were unable to show any difference in colocalization between these mice. Conclusions: We show that colocalization of fluorescent mitochondria and autophagosomes in cell cultures, but not fixed tissues from the RedMIT/GFP-LC3, can be used to detect mitophagy. We used this model to confirm that mitophagy is increased in a mouse model of ADOA. It will be useful for cell based studies of diseases caused by impaired mitochondrial dynamics

Efficiently Calculating Evolutionary Tree Measures Using SAT

We develop techniques to calculate important measures in evolutionary biology by encoding to CNF formulas and using powerful SAT solvers. Comparing evolutionary trees is a necessary step in tree reconstruction algorithms, locating recombination and lateral gene transfer, and in analyzing and visualizing sets of trees. We focus on two popular comparison measures for trees: the hybridization number and the rooted subtree-prune-and-regraft (rSPR) distance. Both have recently been shown to be NP-hard, and effcient algorithms are needed to compute and approximate these measures. We encode these as a Boolean formula such that two trees have hybridization number k (or rSPR distance k) if and only if the corresponding formula is satisfiable. We use state-of-the-art SAT solvers to determine if the formula encoding the measure has a satisfying assignment. Our encoding also provides a rich source of real-world SAT instances, and we include a comparison of several recent solvers (minisat, adaptg2wsat, novelty+p, Walksat, March KS and SATzilla).Postprint (author’s final draft