THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Bucket-handle tear of the triangular fibrocartilage complex : case report of a complex peripheral injury with separation of the distal radioulnar ligaments from the articular disc.

Palmer previously proposed a classification system of triangular fibrocartilage complex (TFCC) injuries that proved to be useful in directing clinical management. However, dorsal peripheral tears (variants of class 1C) were not described and have rarely been reported in the literature since. We herewith present a rare case of bucket-handle tear of the TFCC. To our knowledge, this is the first case demonstrating partial separation of both the palmar and dorsal distal radioulnar ligaments (DRULs) from the articular disc. The particular wrist magnetic resonance (MR) arthrographic findings of this unusual complex peripheral TFCC tear (a variant of both class 1B and 1C) were nicely appreciated upon sagittal reformatted images.

Conservation of complex knotting and slipknotting patterns in proteins.

While analyzing all available protein structures for the presence of knots and slipknots, we detected a strict conservation of complex knotting patterns within and between several protein families despite their large sequence divergence. Because protein folding pathways leading to knotted native protein structures are slower and less efficient than those leading to unknotted proteins with similar size and sequence, the strict conservation of the knotting patterns indicates an important physiological role of knots and slipknots in these proteins. Although little is known about the functional role of knots, recent studies have demonstrated a protein-stabilizing ability of knots and slipknots. Some of the conserved knotting patterns occur in proteins forming transmembrane channels where the slipknot loop seems to strap together the transmembrane helices forming the channel.

Ultrastructure of the membrane attack complex of complement: detection of the tetramolecular C9-polymerizing complex C5b-8.

The ultrastructure of the membrane attack complex (MAC) of complement had been described as representing a hollow cylinder of defined dimensions that is composed of the proteins C5b, C6, C7, C8, and C9. After the characteristic cylindrical structure was identified as polymerized C9 [poly(C9)], the question arose as to the ultrastructural identity and topology of the C9-polymerizing complex C5b-8. An electron microscopic analysis of isolated MAC revealed an asymmetry of individual complexes with respect to their length. Whereas the length of one boundary (+/- SEM) was always 16 +/- 1 nm, the length of the other varied between 16 and 32 nm. In contrast, poly(C9), formed spontaneously from isolated C9, had a uniform tubule length (+/- SEM) of 16 +/- 1 nm. On examination of MAC-phospholipid vesicle complexes, an elongated structure was detected that was closely associated with the poly(C9) tubule and that extended 16-18 nm beyond the torus of the tubule and 28-30 nm above the membrane surface. The width of this structure varied depending on its two-dimensional projection in the electron microscope. By using biotinyl C5b-6 in the formation of the MAC and avidin-coated colloidal gold particles for the ultrastructural analysis, this heretofore unrecognized subunit of the MAC could be identified as the tetramolecular C5b-8 complex. Identification also was achieved by using anti-C5 Fab-coated colloidal gold particles. A similar elongated structure of 25 nm length (above the surface of the membrane) was observed on single C5b-8-vesicle complexes. It is concluded that the C5b-8 complex, which catalyzes poly(C9) formation, constitutes a structure of discrete morphology that remains as such identifiable in the fully assembled MAC, in which it is closely associated with the poly(C9) tubule.

Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

[R74W;R1070W;D1270N]: a new complex allele responsible for cystic fibrosis.

Since the beginning of population screening for CF carriers, it has become apparent that complex CFTR alleles are not uncommon. Deciphering their impact in disease pathogenesis remains a challenge for both clinicians and researchers. We report the observation of a new complex allele p.[R74W+R1070W+D1270N] found in trans with a type 1 mutation and associated with clinical diagnosis of cystic fibrosis in a one year-old Moroccan patient. This case underlines the difficulties in counseling patients with uncommon mutations and the necessity of functional studies to evaluate the structure-function relationships, since the association of several variations in cis can dramatically alter CFTR function.

Secretory IgA's complex roles in immunity and mucosal homeostasis in the gut.

Secretory IgA (SIgA) serves as the first line of defense in protecting the intestinal epithelium from enteric toxins and pathogenic microorganisms. Through a process known as immune exclusion, SIgA promotes the clearance of antigens and pathogenic microorganisms from the intestinal lumen by blocking their access to epithelial receptors, entrapping them in mucus, and facilitating their removal by peristaltic and mucociliary activities. In addition, SIgA functions in mucosal immunity and intestinal homeostasis through mechanisms that have only recently been revealed. In just the past several years, SIgA has been identified as having the capacity to directly quench bacterial virulence factors, influence composition of the intestinal microbiota by Fab-dependent and Fab-independent mechanisms, promote retro-transport of antigens across the intestinal epithelium to dendritic cell subsets in gut-associated lymphoid tissue, and, finally, to downregulate proinflammatory responses normally associated with the uptake of highly pathogenic bacteria and potentially allergenic antigens. This review summarizes the intrinsic biological activities now associated with SIgA and their relationships with immunity and intestinal homeostasis.

Analyzing complex rock slope deformation at Stampa, western Norway, by integrating geomorphology, kinematics and numerical modeling

The unstable rock slope, Stampa, above the village of Flåm, Norway, shows signs of both active and postglacial gravitational deformation over an area of 11 km2. Detailed structural field mapping, annual differential Global Navigation Satellite System (GNSS) surveys, as well as geomorphic analysis of high-resolution digital elevation models based on airborne and terrestrial laser scanning indicate that slope deformation is complex and spatially variable. Numerical modeling was used to investigate the influence of former rockslide activity and to better understand the failure mechanism. Field observations, kinematic analysis and numerical modeling indicate a strong structural control of the unstable area. Based on the integration of the above analyses, we propose that the failure mechanism is dominated by (1) a toppling component, (2) subsiding bilinear wedge failure and (3) planar sliding along the foliation at the toe of the unstable slope. Using differential GNSS, 18 points were measured annually over a period of up to 6 years. Two of these points have an average yearly movement of around 10 mm/year. They are located at the frontal cliff on almost completely detached blocks with volumes smaller than 300,000 m3. Large fractures indicate deep-seated gravitational deformation of volumes reaching several 100 million m3, but the movement rates in these areas are below 2 mm/year. Two different lobes of prehistoric rock slope failures were dated with terrestrial cosmogenic nuclides. While the northern lobe gave an average age of 4,300 years BP, the southern one resulted in two different ages (2,400 and 12,000 years BP), which represent most likely multiple rockfall events. This reflects the currently observable deformation style with unstable blocks in the northern part in between Joasete and Furekamben and no distinct blocks but a high rockfall activity around Ramnanosi in the south. With a relative susceptibility analysis it is concluded that small collapses of blocks along the frontal cliff will be more frequent. Larger collapses of free-standing blocks along the cliff with volumes > 100,000 m3, thus large enough to reach the fjord, cannot be ruled out. A larger collapse involving several million m3 is presently considered of very low likelihood.

Characterization of high osmolarity-induced vacuole fragmentation in "Saccharomyces c."

La majorité des organelles d'une cellule adaptent leur nombre et leur taille pendant les processus de division cellulaire, de trafic vésiculaire ou suite à des changements environnementaux par des processus de fusion et de fragmentation membranaires. Ceci est valable notamment pour le golgi, les mitochondries, les péroxisomes et les lysosomes. La vacuole est le compartiment terminal de la voie endocytaire dans la levure Saccharomyces cerevisiae\ elle correspond aux lysosomes des cellules mammifères. Suite à un choc hyperosmotique, la vacuole se fragmente en plusieurs petites vésicules. Durant ce projet, cette fragmentation a été étudiée en utilisant la technique de microscopie confocale in vivo. J'ai observé que la division de la vacuole se produit d'une façon asymétrique. La première minute après le choc osmotique, les vacuoles rétrécissent et forment des longues invaginations tubulaires. Cette phase est dépendante de la protéine Vps1, un membre de la famille des protéines apparentées à la dynamine, ainsi que d'un gradient transmembranaire de protons. Pendant les 10-15 minutes qui suivent, des vésicules se détachent dans les régions où l'on observe les invaginations pendant la phase initiale. Cette deuxième phase qui mène à la fission des nouveaux compartiments vacuolaires dépend de la production du lipide PI(3,5)P2 par la protéine Fab1. J'ai établi la suite des événements du processus de fragmentation des vacuoles et propose la possibilité d'un rôle régulateur de la protéine kinase cycline-dépendante Pho85.¦En outre, j'ai tenté d'éclaircir plus spécifiquement le rôle de Vps1 pendant la fusion et fission des vacuoles. J'ai trouvé que tous les deux processus sont dépendants de l'activité GTPase de cette protéine. De plus l'association avec la membrane vacuolaire paraît régulée par le cycle d'hydrolyse du GTP. Vps1 peut lier la membrane sans la présence d'un autre facteur protéinique, ce qui permet de conclure à une interaction directe avec des lipides de la membrane. Cette interaction est au moins partiellement effectuée par le domaine GTPase, ce qui est une nouveauté pour un membre de cette famille de protéines. Une deuxième partie de Vps1, nommée insert B, est impliquée dans la liaison à la vacuole, soit par interaction directe avec la membrane, soit par régulation du domaine GTPase. En assumant que Vps1 détienne deux régions capables de liaison aux membranes, je conclus qu'elle pourrait fonctionner comme facteur de « tethering » lors de la fusion des vacuoles.¦-¦La cellule contient plusieurs sous-unités, appelées organelles, possédant chacune une fonction spécifique. Dépendant des processus qui s'y déroulent à l'intérieur, un environnement chimique spécifique est requis. Pour maintenir ces différentes conditions, les organelles sont séparées par des membranes. Lors de la division cellulaire ou en adaptation à des changements de milieu, les organelles doivent être capables de modifier leur morphologie. Cette adaptation a souvent lieu par fusion ou division des organelles. Le même principe est valable pour la vacuole dans la levure. La vacuole est une organelle qui sert principalement au stockage des aliments et à la dégradation des différents composants cellulaires. Alors que la fusion des vacuoles est un processus déjà bien décrit, la fragmentation des vacuoles a jusqu'ici été peu étudiée. Elle peut être induit par un choc osmotique: à cause de la concentration de sel élevé dans le milieu, le cytosol de la levure perd de l'eau. Par un flux d'eau de la vacuole au cytosol, la cellule est capable d'équilibrer celui-ci. Quand la vacuole perd du volume, elle doit réadapter le rapport entre surface membranaire et volume, ce qui se fait efficacement par une fragmentation d'une grande vacuole en plusieurs petites vésicules. Comment ce processus se déroule d'un point de vue morphologique n'a pas été décrit jusqu'à présent. En analysant la fragmentation vacuolaire par microscopie, j'ai trouvé que celle-ci se déroule en deux phases. Pendant la première minute suivant le choc osmotique, les vacuoles rétrécissent et forment des longues invaginations tubulaires. Cette phase dépend de la protéine Vps1, un membre de la famille des protéines apparentées à la dynamine, ainsi que du gradient transmembranaire de protons. Ce gradient s'établit par une pompe membranaire, la V-ATPase, qui transporte des protons dans la vacuole en utilisant l'énergie libérée par hydrolyse d'ATP. Après cette phase initiale, la formation de nouvelles vésicules vacuolaires dépend de la synthèse du lipide PI(3,5)P2.¦Dans la deuxième partie de l'étude, j'ai tenté de décrire comment Vps1 lie la membrane pour effectuer un remodelage de la vacuole. Vps1 est nécessaire pour la fusion et la fragmentation des vacuoles. J'ai découvert que tous les deux processus dépendent de sa capacité d'hydrolyser du GTP. Ainsi l'association avec la membrane est couplée au cycle d'hydrolyse du GTP. Vps1 peut lier la membrane sans la présence d'une autre protéine, et interagit donc très probablement avec les lipides de la membrane. Deux parties différentes de la protéine sont impliquées dans la liaison, dont une, inattendue, le domaine GTPase.¦-¦Numerous organelles undergo membrane fission and fusion events during cell division, vesicular traffic, or in response to changes in environmental conditions. Examples include Golgi (Acharya et al., 1998) mitochondria (Bleazard et al., 1999) peroxisomes (Kuravi et al., 2006) and lysosomes (Ward et al., 1997). In the yeast Saccharomyces cerevisiae the vacuole is the terminal component of the endocytic pathway and corresponds to lysosomes in mammalian cells. Yeast vacuoles fragment into multiple small vesicles in response to a hypertonic shock. This rapid and homogeneous reaction can serve as a model to study the requirements of the fragmentation process. Here, I investigated osmotically induced fragmentation by time-lapse microscopy. I observe that the small fragmentation products originate directly from the large central vacuole by asymmetric scission rather than by consecutive equal divisions and that fragmentation occurs in two distinct phases. During the first minute, vacuoles shrink and generate deep invaginations, leaving behind tubular structures. This phase requires the dynamin-like GTPase Vps1 and the vacuolar proton gradient. In the subsequent 10-15 minutes, vesicles pinch off from the tubular structures in a polarized fashion, directly generating fragmentation products of the final size. This phase depends on the production of phosphatidylinositol- 3,5-bisphosphate by the Fab1 complex. I suggest a possible regulation of vacuole fragmentation by the CDK Pho85. Based on my microscopy study I established a sequential involvement of the different fission factors.¦In addition to the morphological description of vacuole fragmentation I more specifically aimed to shed some light on the role of Vps1 in vacuole fragmentation and fusion. I find that both functions are dependent on the GTPase activity of the protein and that also the membrane association of the dynamin-like protein is coupled to the GTPase cycle. I found that Vps1 has the capacity for direct lipid binding on the vacuole and that this lipid binding is at least partially mediated through residues in the GTPase domain, a complete novelty for a dynamin family member. A second stretch located in the region of insert Β has also membrane-binding activity or regulates the association with the vacuole through the GTPase domain. Under the assumption of two membrane-binding regions I speculate on Vps1 as a possible tethering factor for vacuole fusion.

An experimental study on the shallow-level migmatization of ferrogabbros from the Fuerteventura Basal Complex, Canary Islands

Spectacular shallow-level migmatization of ferrogabbroic rocks occurs in a metamorphic contact aureole of a gabbroic pluton of the Tierra Mala massif (TM) on Fuerteventura (Canary Islands). In order to improve our knowledge of the low pressure melting behavior of gabbroic rocks and to constrain the conditions of migmatization of the TM gabbros, we performed partial melting experiments on a natural ferrogabbro, which is assumed as protolith of the migmatites. The experiments were performed in an internally heated pressure vessel (IHPV) at 200 MPa, 930-1150 degreesC at relatively oxidizing conditions. Distinct amounts of water were added to the charge. From 930 to 1000 degreesC, the observed experimental phases are plagioclase (An(60-70)), clinopyroxene, amphibole (titanian magnesiohastingsites), two Fe-Ti oxides, and a basaltic, K-poor melt. Above 1000 degreesC, amphibole is no longer stable. The first melts are very rich in non-native plagioclase (>70 wt.%). This indicates that at the beginning of partial melting plagioclase is the major phase which is consumed to produce melt. In the experiments, plagioclase is stable up to high temperatures (1060 degreesC) showing increasing An content with temperature. This is not compatible with the natural migmatites, in which An-rich plagioclase is absent in the melanosomes, while amphibole is stable. Our results show that the partial melting of the natural rocks cannot be regarded as an ``in-situ'' process that occurred in a closed system. Considerable amounts of alkalis probably transported by water-rich fluids, derived from the mafic pluton underplating the TM gabbro, were necessary to drive the melting reaction out of the stability range of plagioclase. A partial melting experiment with a migmatite gabbro showing typical ``in-situ'' textures as starting material supports this assumption. Crystallization experiments performed at 1000 degreesC on a glass of the fitised ferrogabbro with different water contents added to the charge show that generally high water activities could be achieved (crystallization of amphibole), independently of the bulk water content, even in a system with very low initial bulk water content (0.3 wt.%). Increasing water contents produce plagioclase richer in An, reduces the modal proportion of plagioclase in the crystallizing assemblage and extends the melt fraction. High melt fractions of >30 wt.% could only be observed in systems with high bulk water contents (> - 2 wt.%). This indicates that the migmatites were generated under water-rich conditions (probably water-saturated), since those migmatites, which are characterized as ``in-situ'' formations, show generally high amounts of leucosomes (>30 wt.%). (C) 2003 Elsevier B.V. All rights reserved.