Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Evolution: plastic sociality in a sweat bee.

How and why do bees become social? A transplant experiment shows that sweat bees can adopt a solitary or social lifestyle in response to their environment.

Accès des patients agés aux urgences: evolution démographique et perspectives médico-ethiques [Access of elderly patients in the emergency department: demographic evolution and ethical perspectives].

Due to actual demographic evolution, emergency departments have to face a dramatic increase in admissions of elderly people. The peculiar medical and socio-demographic characteristics of these old patients emphasize the need of specific decision processes and resources allocation. An individual-based approach, related to significant ethical values, should allow better diagnostic and therapeutic attitudes. Such a way to admit, evaluate and treat older patients implies an active collaboration with patients and their relatives, but also with all medical interveners, including in particular primary care physicians.

Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

MMPBSA decomposition of the binding energy throughout a molecular dynamics simulation of Amyloid-Beta (Aß10-35) aggregation

Recent experiments with amyloid-beta (Aß) peptides indicate that the formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Aß oligomers depend on their structure, which is governed by assembly dynamics. However, a detailed knowledge of the structure of at the atomic level has not been achieved yet due to limitations of current experimental techniques. In this study, replica exchange molecular dynamics simulations are used to identify the expected diversity of dimer conformations of Aß10-35 monomers. The most representative dimer conformation has been used to track the dimer formation process between both monomers. The process has been characterized by means of the evolution of the decomposition of the binding free energy, which provides an energetic profile of the interaction. Dimers undergo a process of reorganization driven basically by inter-chain hydrophobic and hydrophilic interactions and also solvation/desolvation processes.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Host Genes Important to HIV Replication and Evolution.

Recent years have seen a significant increase in understanding of the host genetic and genomic determinants of susceptibility to HIV-1 infection and disease progression, driven in large part by candidate gene studies, genome-wide association studies, genome-wide transcriptome analyses, and large-scale in vitro genome screens. These studies have identified common variants in some host loci that clearly influence disease progression, characterized the scale and dynamics of gene and protein expression changes in response to infection, and provided the first comprehensive catalogs of genes and pathways involved in viral replication. Experimental models of AIDS and studies in natural hosts of primate lentiviruses have complemented and in some cases extended these findings. As the relevant technology continues to progress, the expectation is that such studies will increase in depth (e.g., to include host whole exome and whole genome sequencing) and in breadth (in particular, by integrating multiple data types).

The evolution of unicoloniality in ant societies : recognition, dispersal and population genetic structure

Introduction Societies of ants, bees, wasps and termites dominate many terrestrial ecosystems (Wilson 1971). Their evolutionary and ecological success is based upon the regulation of internal conflicts (e.g. Ratnieks et al. 2006), control of diseases (e.g. Schmid-Hempel 1998) and individual skills and collective intelligence in resource acquisition, nest building and defence (e.g. Camazine 2001). Individuals in social species can pass on their genes not only directly trough their own offspring, but also indirectly by favouring the reproduction of relatives. The inclusive fitness theory of Hamilton (1963; 1964) provides a powerful explanation for the evolution of reproductive altruism and cooperation in groups with related individuals. The same theory also led to the realization that insect societies are subject to internal conflicts over reproduction. Relatedness of less-than-one is not sufficient to eliminate all incentive for individual selfishness. This would indeed require a relatedness of one, as found among cells of an organism (Hardin 1968; Keller 1999). The challenge for evolutionary biology is to understand how groups can prevent or reduce the selfish exploitation of resources by group members, and how societies with low relatedness are maintained. In social insects the evolutionary shift from single- to multiple queens colonies modified the relatedness structure, the dispersal, and the mode of colony founding (e.g. (Crozier & Pamilo 1996). In ants, the most common, and presumably ancestral mode of reproduction is the emission of winged males and females, which found a new colony independently after mating and dispersal flights (Hölldobler & Wilson 1990). The alternative reproductive tactic for ant queens in multiple-queen colonies (polygyne) is to seek to be re-accepted in their natal colonies, where they may remain as additional reproductives or subsequently disperse on foot with part of the colony (budding) (Bourke & Franks 1995; Crozier & Pamilo 1996; Hölldobler & Wilson 1990). Such ant colonies can contain up to several hundred reproductive queens with an even more numerous workforce (Cherix 1980; Cherix 1983). As a consequence in polygynous ants the relatedness among nestmates is very low, and workers raise brood of queens to which they are only distantly related (Crozier & Pamilo 1996; Queller & Strassmann 1998). Therefore workers could increase their inclusive fitness by preferentially caring for their closest relatives and discriminate against less related or foreign individuals (Keller 1997; Queller & Strassmann 2002; Tarpy et al. 2004). However, the bulk of the evidence suggests that social insects do not behave nepotistically, probably because of the costs entailed by decreased colony efficiency or discrimination errors (Keller 1997). Recently, the consensus that nepotistic behaviour does not occur in insect colonies was challenged by a study in the ant Formica fusca (Hannonen & Sundström 2003b) showing that the reproductive share of queens more closely related to workers increases during brood development. However, this pattern can be explained either by nepotism with workers preferentially rearing the brood of more closely related queens or intrinsic differences in the viability of eggs laid by queens. In the first chapter, we designed an experiment to disentangle nepotism and differences in brood viability. We tested if workers prefer to rear their kin when given the choice between highly related and unrelated brood in the ant F. exsecta. We also looked for differences in egg viability among queens and simulated if such differences in egg viability may mistakenly lead to the conclusion that workers behave nepotistically. The acceptance of queens in polygnous ants raises the question whether the varying degree of relatedness affects their share in reproduction. In such colonies workers should favour nestmate queens over foreign queens. Numerous studies have investigated reproductive skew and partitioning of reproduction among queens (Bourke et al. 1997; Fournier et al. 2004; Fournier & Keller 2001; Hammond et al. 2006; Hannonen & Sundström 2003a; Heinze et al. 2001; Kümmerli & Keller 2007; Langer et al. 2004; Pamilo & Seppä 1994; Ross 1988; Ross 1993; Rüppell et al. 2002), yet almost no information is available on whether differences among queens in their relatedness to other colony members affects their share in reproduction. Such data are necessary to compare the relative reproductive success of dispersing and non-dispersing individuals. Moreover, information on whether there is a difference in reproductive success between resident and dispersing queens is also important for our understanding of the genetic structure of ant colonies and the dynamics of within group conflicts. In chapter two, we created single-queen colonies and then introduced a foreign queens originating from another colony kept under similar conditions in order to estimate the rate of queen acceptance into foreign established colonies, and to quantify the reproductive share of resident and introduced queens. An increasing number of studies have investigated the discrimination ability between ant workers (e.g. Holzer et al. 2006; Pedersen et al. 2006), but few have addressed the recognition and discrimination behaviour of workers towards reproductive individuals entering colonies (Bennett 1988; Brown et al. 2003; Evans 1996; Fortelius et al. 1993; Kikuchi et al. 2007; Rosengren & Pamilo 1986; Stuart et al. 1993; Sundström 1997; Vásquez & Silverman in press). These studies are important, because accepting new queens will generally have a large impact on colony kin structure and inclusive fitness of workers (Heinze & Keller 2000). In chapter three, we examined whether resident workers reject young foreign queens that enter into their nest. We introduced mated queens into their natal nest, a foreign-female producing nest, or a foreign male-producing nest and measured their survival. In addition, we also introduced young virgin and mated queens into their natal nest to examine whether the mating status of the queens influences their survival and acceptance by workers. On top of polgyny, some ant species have evolved an extraordinary social organization called 'unicoloniality' (Hölldobler & Wilson 1977; Pedersen et al. 2006). In unicolonial ants, intercolony borders are absent and workers and queens mix among the physically separated nests, such that nests form one large supercolony. Super-colonies can become very large, so that direct cooperative interactions are impossible between individuals of distant nests. Unicoloniality is an evolutionary paradox and a potential problem for kin selection theory because the mixing of queens and workers between nests leads to extremely low relatedness among nestmates (Bourke & Franks 1995; Crozier & Pamilo 1996; Keller 1995). A better understanding of the evolution and maintenance of unicoloniality requests detailed information on the discrimination behavior, dispersal, population structure, and the scale of competition. Cryptic genetic population structure may provide important information on the relevant scale to be considered when measuring relatedness and the role of kin selection. Theoretical studies have shown that relatedness should be measured at the level of the `economic neighborhood', which is the scale at which intraspecific competition generally takes place (Griffin & West 2002; Kelly 1994; Queller 1994; Taylor 1992). In chapter four, we conducted alarge-scale study to determine whether the unicolonial ant Formica paralugubris forms populations that are organised in discrete supercolonies or whether there is a continuous gradation in the level of aggression that may correlate with genetic isolation by distance and/or spatial distance between nests. In chapter five, we investigated the fine-scale population structure in three populations of F. paralugubris. We have developed mitochondria) markers, which together with the nuclear markers allowed us to detect cryptic genetic clusters of nests, to obtain more precise information on the genetic differentiation within populations, and to separate male and female gene flow. These new data provide important information on the scale to be considered when measuring relatedness in native unicolonial populations.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Evolutionary dynamics of interlinked public goods traits: an experimental study of siderophore production in Pseudomonas aeruginosa.

Public goods cooperation is common in microbes, and there is much interest in understanding how such traits evolve. Research in recent years has identified several important factors that shape the evolutionary dynamics of such systems, yet few studies have investigated scenarios involving interactions between multiple public goods. Here, we offer general predictions about the evolutionary trajectories of two public goods traits having positive, negative or neutral regulatory influence on one another's expression, and we report on a test of some of our predictions in the context of Pseudomonas aeruginosa's production of two interlinked iron-scavenging siderophores. First, we confirmed that both pyoverdine and pyochelin siderophores do operate as public goods under appropriate environmental conditions. We then tracked their production in lines experimentally evolved under different iron-limitation regimes known to favour different siderophore expression profiles. Under strong iron limitation, where pyoverdine represses pyochelin, we saw a decline in pyoverdine and a concomitant increase in pyochelin - consistent with expansion of pyoverdine-defective cheats derepressed for pyochelin. Under moderate iron limitation, pyochelin declined - again consistent with an expected cheat invasion scenario - but there was no concomitant shift in pyoverdine because cross-suppression between the traits is unidirectional only. Alternating exposure to strong and moderate iron limitation caused qualitatively similar though lesser shifts compared to the constant-environment regimes. Our results confirm that the regulatory interconnections between public goods traits can significantly modulate the course of evolution, yet also suggest how we can start to predict the impacts such complexities will have on phenotypic divergence and community stability.

Generation and nonlinear evolution of shore-oblique/transverse sand bars

The coupling between topography, waves and currents in the surf zone may selforganize to produce the formation of shore-transverse or shore-oblique sand bars on an otherwise alongshore uniform beach. In the absence of shore-parallel bars, this has been shown by previous studies of linear stability analysis, but is now extended to the finite-amplitude regime. To this end, a nonlinear model coupling wave transformation and breaking, a shallow-water equations solver, sediment transport and bed updating is developed. The sediment flux consists of a stirring factor multiplied by the depthaveraged current plus a downslope correction. It is found that the cross-shore profile of the ratio of stirring factor to water depth together with the wave incidence angle primarily determine the shape and the type of bars, either transverse or oblique to the shore. In the latter case, they can open an acute angle against the current (upcurrent oriented) or with the current (down-current oriented). At the initial stages of development, both the intensity of the instability which is responsible for the formation of the bars and the damping due to downslope transport grow at a similar rate with bar amplitude, the former being somewhat stronger. As bars keep on growing, their finite-amplitude shape either enhances downslope transport or weakens the instability mechanism so that an equilibrium between both opposing tendencies occurs, leading to a final saturated amplitude. The overall shape of the saturated bars in plan view is similar to that of the small-amplitude ones. However, the final spacings may be up to a factor of 2 larger and final celerities can also be about a factor of 2 smaller or larger. In the case of alongshore migrating bars, the asymmetry of the longshore sections, the lee being steeper than the stoss, is well reproduced. Complex dynamics with merging and splitting of individual bars sometimes occur. Finally, in the case of shore-normal incidence the rip currents in the troughs between the bars are jet-like while the onshore return flow is wider and weaker as is observed in nature.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Organic matter processing and soil evolution in a braided river system

Traditionally, braided river research has considered flow, sediment transport processes and, recently, vegetation dynamics in relation to river morphodynamics. However, if considering the development of woody vegetated patches over a time scale of decades, we must consider the extent to which soil forming processes, particularly related to soil organic matter, impact the alluvial geomorphic-vegetation system. Here we quantify the soil organic matter processing (humification) that occurs on young alluvial landforms. We sampled different geomorphic units, ranging from the active river channel to established river terraces in a braided river system. For each geomorphic unit, soil pits were used to sample sediment/soil layers that were analysed in terms of grain size (<2mm) and organic matter quantity and quality (RockEval method). A principal components analysis was used to identify patterns in the dataset. Results suggest that during the succession from bare river gravels to a terrace soil, there is a transition from small amounts of external organic matter supply provided by sedimentation processes (e.g. organic matter transported in suspension and deposited on bars), to large amounts of autogenic in situ organic matter production due to plant colonisation. This appears to change the time scale and pathways of alluvial succession (bio-geomorphic succession). However, this process is complicated by: the ongoing possibility of local sedimentation, which can serve to isolate surface layers via aggradation from the exogenic supply; and erosion which tends to create fresh deposits upon which organic matter processing must re-start. The result is a complex pattern of organic matter states as well as a general lack of any clear chronosequence within the active river corridor. This state reflects the continual battle between deposition events that can isolate organic matter from the surface, erosion events that can destroy accumulating organic matter and the early ecosystem processes necessary to assist the co-evolution of soil and vegetation. A key question emerges over the extent to which the fresh organic matter deposited in the active zone is capable of significantly transforming the local geochemical environment sufficiently to accelerate soil development.

Microbial communities in a porphyry copper tailings impoundment and their impact on the geochemical dynamics of the mine waste

The distribution and diversity of acidophilic bacteria of a tailings impoundment at the La Andina copper mine, Chile, was examined. The tailings have low sulfide (1.7% pyrite equivalent) and carbonate (1.4% calcite equivalent) contents and are stratified into three distinct zones: a surface (0-70-80 cm) `oxidation zone' characterized by low-pH (2.5-4), a `neutralization zone' (70-80 to 300-400 cm) and an unaltered `primary zone' below 400 cm. A combined cultivation-dependent and biomolecular approach (terminal restriction enzyme fragment length polymorphism and 16S rRNA clone library analysis) was used to characterize the indigenous prokaryotic communities in the mine tailings. Total cell counts showed that the microbial biomass was greatest in the top 125 cm of the tailings. The largest numbers of bacteria (10(9) g(-1) dry weight of tailings) were found at the oxidation front (the junction between the oxidation and neutralization zones), where sulfide minerals and oxygen were both present. The dominant iron-/sulfur-oxidizing bacteria identified at the oxidation front included bacteria of the genus Leptospirillum (detected by molecular methods), and Gram-positive iron-oxidizing acidophiles related to Sulfobacillus (identified both by molecular and cultivation methods). Acidithiobacillus ferrooxidans was also detected, albeit in relatively small numbers. Heterotrophic acidophiles related to Acidobacterium capsulatum were found by molecular methods, while another Acidobacterium-like bacterium and an Acidiphilium sp. were isolated from oxidation zone samples. A conceptual model was developed, based on microbiological and geochemical data derived from the tailings, to account for the biogeochemical evolution of the Piuquenes tailings impoundment.

Population viscosity can promote the evolution of altruistic sterile helpers and eusociality.

Because it increases relatedness between interacting individuals, population viscosity has been proposed to favour the evolution of altruistic helping. However, because it increases local competition between relatives, population viscosity may also act as a brake for the evolution of helping behaviours. In simple models, the kin selected fecundity benefits of helping are exactly cancelled out by the cost of increased competition between relatives when helping occurs after dispersal. This result has lead to the widespread view, especially among people working with social organisms, that special conditions are required for the evolution of altruism. Here, we re-examine this result by constructing a simple population genetic model where we analyse whether the evolution of a sterile worker caste (i.e. an extreme case of altruism) can be selected for by limited dispersal. We show that a sterile worker caste can be selected for even under the simplest life-cycle assumptions. This has relevant consequences for our understanding of the evolution of altruism in social organisms, as many social insects are characterized by limited dispersal and significant genetic population structure.