Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection.

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Effect of immune pressure on hepatitis C virus evolution: insights from a single-source outbreak.

The host's immune response to hepatitis C virus (HCV) can result in the selection of characteristic mutations (adaptations) that enable the virus to escape this response. The ability of the virus to mutate at these sites is dependent on the incoming virus, the fitness cost incurred by the mutation, and the benefit to the virus in escaping the response. Studies examining viral adaptation in chronic HCV infection have shown that these characteristic immune escape mutations can be observed at the population level as human leukocyte antigen (HLA)-specific viral polymorphisms. We examined 63 individuals with chronic HCV infection who were infected from a single HCV genotype 1b source. Our aim was to determine the extent to which the host's immune pressure affects HCV diversity and the ways in which the sequence of the incoming virus, including preexisting escape mutations, can influence subsequent mutations in recipients and infection outcomes. Conclusion: HCV sequences from these individuals revealed 29 significant associations between specific HLA types within the new hosts and variations within their viruses, which likely represent new viral adaptations. These associations did not overlap with previously reported adaptations for genotypes 1a and 3a and possibly reflected a combination of constraint due to the incoming virus and genetic distance between the strains. However, these sites accounted for only a portion of the sites in which viral diversity was observed in the new hosts. Furthermore, preexisting viral adaptations in the incoming (source) virus likely influenced the outcomes in the new hosts.

Coupled reproductive ecologies in a plant-pollinator-seed predator system

Résumé de la thèseBien que le mutualisme puisse être considéré comme une relation harmonieuse entre différentes espèces, son étude révèle plutôt une exploitation réciproque où chaque partenaire tente de maximiser ses bénéfices tout en réduisant ses coûts. Dans ce contexte, l'identification des facteurs qui favorisent ou contrarient, au cours de l'évolution, une issue mutualiste est une étape majeure pour pouvoir reconstruire les étapes clés menant à l'apparition et au maintien des interactions mutualistes. Le but de ce doctorat était l'identification des traits phénotypiques qui permettent à la plante Silene latofolia (Caryophyllacée)et à son pollinisateur - prédateur de graines, la phalène Hadena bicruris (Noctuidé), d'augmenter les bénéfices nets que chacun retire de l'interaction. Ce système d'étude est particulièrement bien approprié à l'étude de ces traits, car on peut assez facilement estimer la qualité et la quantité des descendants (fitness) des deux partenaires. En effet, la femelle papillon pond un oeuf dans la fleur qu'elle pollinise et sa larve se développe dans le fruit, consommant les graines de la plante. Ainsi, sur une même plante, il est possible d'estimer les succès respectifs de la plante et du papillon à obtenir une descendance. De plus, le conflit d'intérêt autour des graines qui sont indispensables, à la fois à la plante et au papillon, peut stimuler l'évolution de traits qui limitent la surexploitation réciproque des partenaires. Dans une première étude, j'ai montré que le papillon mâle était un pollinisateur efficace de S. latifolia et qu'ainsi, il permettait à la plante d'augmenter le nombre de graines produites (i.e.bénéfice) sans pour autant augmenter la quantité de larves sur la plante. Dans ce système, les papillons pondent un seul oeuf par fleur, déposé soit à l'intérieur de la fleur, dans le tube de corolle, soit sur le pétale. Ma seconde étude montre que les plantes répondent différemment à la présence des oeufs suivant leur position. Aussi, quand l'oeuf est placé dans la fleur, la plante a davantage tendance à ne pas développer le fruit de la fleur infesté ou bien à produire des fruits plus petits que lorsque l'oeuf est placé sur le pétale. Enfin, j'ai montré que la femelle du papillon pond plus souvent sur le pétale lorsque elle visite des fleurs dotées d'un long tube de corolle, et que les larves issues de ces oeufs ont moins de chances de réussir à pénétrer dans le fruit que les larves issues des oeufs placés à l'intérieur de la fleur. Aussi, la variation observée du site de ponte pourrait être causé par la morphologie de la fleur qui contraint le papillon à pondre sur le pétale. Vu dans leur ensemble, les résultats obtenus pendant ce doctorat suggèrent que la participation des mâles à la pollination, l'absence de développement des fruits et la profondeur du tube de corolle pourraient réduire les coûts que S. latifolia subit dans son interaction avec H. bicruris. Par ailleurs, je n'ai pas détecté de mécanismes qui permettraient au papillon de réduire les coûts que la plante pourrait lui imposer. La prochaine étape serait de déterminer l'effet des traits identifiés dans ce doctorat sur la fitness globale de la plante et du papillon pour estimer pleinement leur efficacité à réduire les coûts et à favoriser une issue mutualiste. De même, il faudrait évaluer l'effet de ces traits en populations naturelles pour identifier le rôle des facteurs environnementaux sur leur efficacité.AbstractAlthough mutualisms can be regarded as harmonious relationships between the interacting partners, they are best conceptualized as reciprocal exploitations in which each partner attempts to increase its own benefits and decrease its costs. To date, identifying the factors which promote or discourage mutualistic outcomes remains a major goal to reconstruct the ecological conditions leading to mutualisms. The aim of this PhD thesis was to identify phenotypic traits that may increase the net benefits of each partner in the interaction between the plant Silene latifolia (Caryophyllaceae) and its pollinator / seed predator, the moth Hadena bicruris (Noctuidae). This study system is particularly well suited because the fitness of both interacting species can be assessed. The female moth lays its egg in the flower it pollinated, and its offspring grows in the fruit, feeding on the seeds of the plant, which allows for the follow-up of both larva and fruit fates. Furthermore, the inherent conflict of interest over the seeds as plant progeny vs. larval resource may stimulate the evolution of traits that reduce overexploitation in both the moth and plant. In a first study, I show that male moths are efficient pollinators, hence increasing seed production without increasing oviposition. The contribution of male moths to pollination might thus improve the net benefits of the interaction for the host plant. Females of the H. bicruris moth lay a single egg per flower, and place it either inside the corolla tube or on the petal. My second study shows that plants are more likely to abort the infested flower or to produce a smaller fruit when the egg was experimentally placed inside the flower compared to plants that received an egg on the petal. Finally, female moths were found to lay their eggs more frequently on the petal when visiting a flower with a deep corolla tube, and larvae hatching from these eggs less likely to successfully attack the fruit. Variation in egg position on the flower may thus be the result of a constraint imposed by floral morphology. Overall, this PhD work suggests that the pollination by male moths, flower abortion, and deep corolla tube may efficiently reduce the costs experienced by S. latifolia in its interaction with H. bicruris. Interestingly, no apparent mechanism of costs reduction was detected for the moth. Further studies should focus on the effects of these traits (i) in the long term fitness of both the plant and the insect and (ii) their interactions with environmental factors (biotic and abiotic) that may affect their efficiency in natural populations.

Kernel-based methods for change detection in remote sensing images

Nowadays, the joint exploitation of images acquired daily by remote sensing instruments and of images available from archives allows a detailed monitoring of the transitions occurring at the surface of the Earth. These modifications of the land cover generate spectral discrepancies that can be detected via the analysis of remote sensing images. Independently from the origin of the images and of type of surface change, a correct processing of such data implies the adoption of flexible, robust and possibly nonlinear method, to correctly account for the complex statistical relationships characterizing the pixels of the images. This Thesis deals with the development and the application of advanced statistical methods for multi-temporal optical remote sensing image processing tasks. Three different families of machine learning models have been explored and fundamental solutions for change detection problems are provided. In the first part, change detection with user supervision has been considered. In a first application, a nonlinear classifier has been applied with the intent of precisely delineating flooded regions from a pair of images. In a second case study, the spatial context of each pixel has been injected into another nonlinear classifier to obtain a precise mapping of new urban structures. In both cases, the user provides the classifier with examples of what he believes has changed or not. In the second part, a completely automatic and unsupervised method for precise binary detection of changes has been proposed. The technique allows a very accurate mapping without any user intervention, resulting particularly useful when readiness and reaction times of the system are a crucial constraint. In the third, the problem of statistical distributions shifting between acquisitions is studied. Two approaches to transform the couple of bi-temporal images and reduce their differences unrelated to changes in land cover are studied. The methods align the distributions of the images, so that the pixel-wise comparison could be carried out with higher accuracy. Furthermore, the second method can deal with images from different sensors, no matter the dimensionality of the data nor the spectral information content. This opens the doors to possible solutions for a crucial problem in the field: detecting changes when the images have been acquired by two different sensors.

The evolutionary history of the CD209 (DC-SIGN) family in humans and non-human primates

The CD209 gene family that encodes C-type lectins in primates includes CD209 (DC-SIGN), CD209L (L-SIGN) and CD209L2. Understanding the evolution of these genes can help understand the duplication events generating this family, the process leading to the repeated neck region and identify protein domains under selective pressure. We compiled sequences from 14 primates representing 40 million years of evolution and from three non-primate mammal species. Phylogenetic analyses used Bayesian inference, and nucleotide substitutional patterns were assessed by codon-based maximum likelihood. Analyses suggest that CD209 genes emerged from a first duplication event in the common ancestor of anthropoids, yielding CD209L2 and an ancestral CD209 gene, which, in turn, duplicated in the common Old World primate ancestor, giving rise to CD209L and CD209. K(A)/K(S) values averaged over the entire tree were 0.43 (CD209), 0.52 (CD209L) and 0.35 (CD209L2), consistent with overall signatures of purifying selection. We also assessed the Toll-like receptor (TLR) gene family, which shares with CD209 genes a common profile of evolutionary constraint. The general feature of purifying selection of CD209 genes, despite an apparent redundancy (gene absence and gene loss), may reflect the need to faithfully recognize a multiplicity of pathogen motifs, commensals and a number of self-antigens

Trait decoupling promotes evolutionary diversification of the trophic and acoustic system of damselfishes.

Trait decoupling, wherein evolutionary release of constraints permits specialization of formerly integrated structures, represents a major conceptual framework for interpreting patterns of organismal diversity. However, few empirical tests of this hypothesis exist. A central prediction, that the tempo of morphological evolution and ecological diversification should increase following decoupling events, remains inadequately tested. In damselfishes (Pomacentridae), a ceratomandibular ligament links the hyoid bar and lower jaws, coupling two main morphofunctional units directly involved in both feeding and sound production. Here, we test the decoupling hypothesis by examining the evolutionary consequences of the loss of the ceratomandibular ligament in multiple damselfish lineages. As predicted, we find that rates of morphological evolution of trophic structures increased following the loss of the ligament. However, this increase in evolutionary rate is not associated with an increase in trophic breadth, but rather with morphofunctional specialization for the capture of zooplanktonic prey. Lineages lacking the ceratomandibular ligament also shows different acoustic signals (i.e. higher variation of pulse periods) from others, resulting in an increase of the acoustic diversity across the family. Our results support the idea that trait decoupling can increase morphological and behavioural diversity through increased specialization rather than the generation of novel ecotypes.

Segmental copy number variation shapes tissue transcriptomes.

Copy number variation (CNV) is a key source of genetic diversity, but a comprehensive understanding of its phenotypic effect is only beginning to emerge. We have generated a CNV map in wild mice and classical inbred strains. Genome-wide expression data from six major organs show not only that expression of genes within CNVs tend to correlate with copy number changes, but also that CNVs influence the expression of genes in their vicinity, an effect that extends up to half a megabase. Genes within CNVs show lower expression and more specific spatial expression patterns than genes mapping elsewhere. Our analyses reveal differential constraint on copy number changes of genes expressed in different tissues. Dosage alterations of brain-expressed genes are less frequent than those of other genes and are buffered by tighter transcriptional regulation. Our study provides initial evidence that CNVs shape tissue transcriptomes on a global scale.

Signal search analysis server.

Signal search analysis is a general method to discover and characterize sequence motifs that are positionally correlated with a functional site (e.g. a transcription or translation start site). The method has played an instrumental role in the analysis of eukaryotic promoter elements. The signal search analysis server provides access to four different computer programs as well as to a large number of precompiled functional site collections. The programs offered allow: (i) the identification of non-random sequence regions under evolutionary constraint; (ii) the detection of consensus sequence-based motifs that are over- or under-represented at a particular distance from a functional site; (iii) the analysis of the positional distribution of a consensus sequence- or weight matrix-based sequence motif around a functional site; and (iv) the optimization of a weight matrix description of a locally over-represented sequence motif. These programs can be accessed at: http://www.isrec.isb-sib.ch/ssa/.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Modeling complex wells with the multi-scale finite-volume method

In this paper, an extension of the multi-scale finite-volume (MSFV) method is devised, which allows to Simulate flow and transport in reservoirs with complex well configurations. The new framework fits nicely into the data Structure of the original MSFV method,and has the important property that large patches covering the whole well are not required. For each well. an additional degree of freedom is introduced. While the treatment of pressure-constraint wells is trivial (the well-bore reference pressure is explicitly specified), additional equations have to be solved to obtain the unknown well-bore pressure of rate-constraint wells. Numerical Simulations of test cases with multiple complex wells demonstrate the ability of the new algorithm to capture the interference between the various wells and the reservoir accurately. (c) 2008 Elsevier Inc. All rights reserved.

Mating triggers dynamic immune regulations in wood ant queens.

Mating can affect female immunity in multiple ways. On the one hand, the immune system may be activated by pathogens transmitted during mating, sperm and seminal proteins, or wounds inflicted by males. On the other hand, immune defences may also be down-regulated to reallocate resources to reproduction. Ants are interesting models to study post-mating immune regulation because queens mate early in life, store sperm for many years, and use it until their death many years later, while males typically die after mating. This long-term commitment between queens and their mates limits the opportunity for sexual conflict but raises the new constraint of long-term sperm survival. In this study, we examine experimentally the effect of mating on immunity in wood ant queens. Specifically, we compared the phenoloxidase and antibacterial activities of mated and virgin Formica paralugubris queens. Queens had reduced levels of active phenoloxidase after mating, but elevated antibacterial activity 7 days after mating. These results indicate that the process of mating, dealation and ovary activation triggers dynamic patterns of immune regulation in ant queens that probably reflect functional responses to mating and pathogen exposure that are independent of sexual conflict.

Topological locking restrains replication fork reversal.

Two-dimensional agarose gel electrophoresis, psoralen cross-linking, and electron microscopy were used to study the effects of positive supercoiling on fork reversal in isolated replication intermediates of bacterial DNA plasmids. The results obtained demonstrate that the formation of Holliday-like junctions at both forks of a replication bubble creates a topological constraint that prevents further regression of the forks. We propose that this topological locking of replication intermediates provides a biological safety mechanism that protects DNA molecules against extensive fork reversals.

Patient versus general practitioner perception of problems with treatment adherence in type 2 diabetes: from adherence to concordance.

OBJECTIVES: To determine the prevalence of problems with treatment adherence among type-2 diabetic patients with regards to medication, dietary advice, and physical activity; to identify the associated clinical and psychosocial factors; and to investigate the degree of agreement between patient-perceived and GP-perceived adherence. METHODS: Consecutive patients were solicited during visits to 39 GPs. In total, 521 patients self-reported on treatment adherence, anxiety and depression, and disease perception. The GPs reported clinical and laboratory data and patients' adherence. A multivariate analysis identified the factors associated with adherence problems. RESULTS: Problems of adherence to medication, dietary advice, and physical activity recommendations were reported by 17%, 62%, and 47% of the patients, respectively. Six independent factors were found associated with adherence problems: young age, body-mass index (BMI) > 30 kg/m(2), glycosylated haemoglobin (HbA(1c)) > 8%, single life, depression, and perception of medication as a constraint. Agreement between patients' and GPs' assessments of treatment problems reached 70%. CONCLUSION: In type 2 diabetes, problems with dietary advice or physical activity are far more frequent than problems with medication, and not all physicians are fully aware of patients' problems. More active listening and shared decision-making should enhance adherence and improve outcomes.