KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

Contracts, Reference Points, and Competition - Behavioral Effects of the Fundamental Transformation

In this paper we study the role of incomplete ex ante contracts for ex post trade. Previous experimental evidence indicates that a contract provides a reference point for entitlements when the terms are negotiated in a competitive market. We show that this finding no longer holds when the terms are determined in a non-competitive way. Our results imply that the presence of a "fundamental transformation" (i.e., the transition from a competitive market to a bilateral relationship) is important for a contract to become a reference point. To the best of our knowledge this behavioral aspect of the fundamental transformation has not been shown before.

Using high-resolution quantitative mapping of R1 as an index of cortical myelination

A fundamental tenet of neuroscience is that cortical functional differentiation is related to the cross-areal differences in cyto-, receptor-, and myeloarchitectonics that are observed in ex-vivo preparations. An ongoing challenge is to create noninvasive magnetic resonance (MR) imaging techniques that offer sufficient resolution, tissue contrast, accuracy and precision to allow for characterization of cortical architecture over an entire living human brain. One exciting development is the advent of fast, high-resolution quantitative mapping of basic MR parameters that reflect cortical myeloarchitecture. Here, we outline some of the theoretical and technical advances underlying this technique, particularly in terms of measuring and correcting for transmit and receive radio frequency field inhomogeneities. We also discuss new directions in analytic techniques, including higher resolution reconstructions of the cortical surface. We then discuss two recent applications of this technique. The first compares individual and group myelin maps to functional retinotopic maps in the same individuals, demonstrating a close relationship between functionally and myeloarchitectonically defined areal boundaries (as well as revealing an interesting disparity in a highly studied visual area). The second combines tonotopic and myeloarchitectonic mapping to localize primary auditory areas in individual healthy adults, using a similar strategy as combined electrophysiological and post-mortem myeloarchitectonic studies in non-human primates.

Evaluation of the caloric test by combining 3 response parameters.

The slow-phase velocity of nystagmus is one of the most sensitive parameters of vestibular function and is currently the standard for evaluating the caloric test. However, the assessment of this parameter requires recording the response by using nystagmography. The aim of this study was to evaluate whether frequency and duration of the caloric nystagmus, as measured by using a clinical test with Frenzel glasses, could predict the result of the recorded test. The retrospective analysis of 222 caloric test results recorded by means of electronystagmography has shown a good association between the 3 parameters for unilateral weakness. The asymmetry observed in the velocity can be predicted by a combination of frequency and duration. On the other hand, no relationship was observed between the parameters for directional preponderance. These results indicate that a clinical caloric test with frequency and duration as parameters can be used to predict the unilateral weakness, which would be obtained by use of nystagmography. We propose an evaluation of the caloric test on the basis of diagrams combining the 3 response parameters.

Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

Gait and foot clearance parameters obtained using shoe-worn inertial sensors in a large-population sample of older adults

In order to distinguish dysfunctional gait; clinicians require a measure of reference gait parameters for each population. This study provided normative values for widely used parameters in more than 1400 able-bodied adults over the age of 65. We also measured the foot clearance parameters (i.e., height of the foot above ground during swing phase) that are crucial to understand the complex relationship between gait and falls as well as obstacle negotiation strategies. We used a shoe-worn inertial sensor on each foot and previously validated algorithms to extract the gait parameters during 20 m walking trials in a corridor at a self-selected pace. We investigated the difference of the gait parameters between male and female participants by considering the effect of age and height factors. Besides; we examined the inter-relation of the clearance parameters with the gait speed. The sample size and breadth of gait parameters provided in this study offer a unique reference resource for the researchers.

Testing the competitive exclusion principle using various niche parameters in a native (Natrix maura) and an introduced (N. tessellata) colubrid

Despite the increase of animal and plant introductions worldwide and the strong augmentation of the reptile trade, few invasive snake populations have been studied. Dice snakes (Natrix tessellata) were introduced to the shores of Lake Geneva (Switzerland) in the early 1920s, and are now well established. This region of introduction was previously inhabited by Viperine snakes (N. maura). Ever since these two species have been under monitoring (which began in 1996) the Viperine snake population has shown drastic decline. We examine here the possibility of trophic competition by analysing diet composition, prey size and trophic niche overlap. Spatial distribution is also assessed in order to address the question of spatial competitive exclusion. We found very similar diets, and thus a high trophic niche overlap, indicating no partitioning of the trophic resource. No arguments in favour of spatial competitive exclusion were found. Our study suggests that trophic competition may occur between the two natricines and that it may give an explanation for the drastic decline of the Viperine snake in this area. Other pathways potentially playing a role in the exclusion of the Viperine snake are discussed.

Quantitative integration of high-resolution hydrogeophysical data through Monte-Carlo-type conditional simulations

Geophysical techniques can help to bridge the inherent gap with regard to spatial resolution and the range of coverage that plagues classical hydrological methods. This has lead to the emergence of the new and rapidly growing field of hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters and their inherent trade-off between resolution and range the fundamental usefulness of multi-method hydrogeophysical surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative integration of the resulting vast and diverse database in order to obtain a unified model of the probed subsurface region that is internally consistent with all available data. To address this problem, we have developed a strategy towards hydrogeophysical data integration based on Monte-Carlo-type conditional stochastic simulation that we consider to be particularly suitable for local-scale studies characterized by high-resolution and high-quality datasets. Monte-Carlo-based optimization techniques are flexible and versatile, allow for accounting for a wide variety of data and constraints of differing resolution and hardness and thus have the potential of providing, in a geostatistical sense, highly detailed and realistic models of the pertinent target parameter distributions. Compared to more conventional approaches of this kind, our approach provides significant advancements in the way that the larger-scale deterministic information resolved by the hydrogeophysical data can be accounted for, which represents an inherently problematic, and as of yet unresolved, aspect of Monte-Carlo-type conditional simulation techniques. We present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to corresponding field data collected at the Boise Hydrogeophysical Research Site near Boise, Idaho, USA.

Systèmes de délivrance des médicaments pour le segment antérieur: bases fondamentales et applications cliniques [Drug delivery systems to target the anterior segment of the eye: fundamental bases and clinical applications].

The development of new drug delivery systems to target the anterior segment of the eye may offer many advantages: to increase the biodisponibility of the drug, to allow the penetration of drug that cannot be formulated as solutions, to obtain constant and sustained drug release, to achieve higher local concentrations without systemic effects, to target more specifically one tissue or cell type, to reduce the frequency of instillation and therefore increase the observance and comfort of the patient while reducing side effects of frequent instillation. Several approaches are developed, aiming to increase the corneal contact time by modified formulation or reservoir systems, or by increasing the tissue permeability using iontophoresis. To date, no ocular drug delivery system is ideal for all purposes. To maximize treatment efficacy, careful evaluation of the specific pathological condition, the targeted Intraocular tissue and the location of the most severe pathology must be made before selecting the method of delivery most suitable for each individual patient.

Comprehensive analysis of proteins secreted by dermatophytes in a protein medium which to some extent mimic "in vivo" growth parameters

RESUME : Les dermatophytes sont les agents infectieux les plus fréquents responsables de la plupart des mycoses superficielles chez les humains et chez les animaux. Ces infections, dermatophytoses, également appelées tineas ou teignes, sont fréquentes et causent des problèmes de santé publique au niveau mondial. La capacité d'envahir et de progresser au sein des structures kératinisées est probablement liée à la sécrétion de différentes enzymes kératinolytiques, qui sont considérées comme la principale caractéristique liée à la pathogénicité de ces champignons. L'objectif de ma thèse a été premièrement de progresser dans l'identification et la caractérisation des nouvelles protéines sécrétées, afin de mieux comprendre a) la capacité globale des dermatophytes à envahir les structures kératinisées, et b) les différences dans la virulence et la spécificité d'hôte que présentent les espèces étudiées .Pour progresser dans l'identification et la caractérisation de ces nouvelles protéines, les secretomes de six espèces de dermatophytes (Trichophyton rubrum, Trichophyton violaceum, Trichophyton soudanense, Trichophyton equinum, Arthroderma vanbreuseghemii et Trichophyton tonsurans) ont été étudiés. Bien qu'il y ait un niveau globalement élevé de similitude entre les protéases sécrétées, les différentes espèces de dermatophytes sécrètent des profiles protéiques distincts lorsqu'elles sont cultivées dans les mêmes conditions de culture, et donc une signature spécifique a pu être associé à chaque espèce. Ces profiles ont été un outil avantageux pour identifier et cartographier les protéines orthologues aux six espèces et ont aussi permit la discrimination d'espèces très proches comme T. tonsurans et T. equinum qui ne peuvent pas être différenciées par l'ADN ribosomal. Ce travail également présente ce que l'on croit être la première identification global des protéines sécrétées par les dermatophytes dans des conditions de culture que incitent l'activité protéolytique extracellulaire. Ce catalogue de protéines, comprenant des endo- and exo- proteases, autres hydrolases, oxydoreductases et des protéines avec fonction inconnue, représente probablement le spectre d'enzymes qui permettent la dégradation des tissus kératinisés en composés qui peuvent être assimilés par le champignon. Les résultats suggèrent qu'un changement écologique pourrait être associé à une expression différentielle des gènes codant les protéines sécrétées, en particulier, les protéases, plutôt qu'à des divergences génétiques au niveau des gènes codant les protéines orthologues. Une sécrétion différentielle des protéines par les dermatophytes pourrait également être responsable de la variabilité inflammatoire qui causent ces agents infectieux chez les différents hôtes. Par conséquent, les protéines identifiées ici sont également importantes pour faire la lumière sur la réponse immunitaire de l'hôte au cours du processus infectieux. SUMMARY : Dermatophytes are the most common infectious agents responsible for superficial mycosis in humans and animals. Dermatophytoses, also called tineas or ringworm, are frequent and cause public health problems worldwide. The secretion of different keratinolytic enzymes is believed to be a key pathogenicity-related characteristic of these fungi. The aim of this work was first to progress in the identification and characterization of novel secreted proteins, in order to better understand a) the overall capability of dermatophytes to invade keratinised structures, and b) differences in virulence and host-specificity of the investigated species. To progress in the identification and characterization of novel proteins, the secretomes from Trichophyton rubrum, Trichophyton violaceum, Trichophyton soudanense, Trichophyton equinum, Arthroderma vanbreuseghemii and Trichophyton tonsurans were studied. Although there is a high global level of similarity among the secreted proteases, different dermatophyte species produce distinct patterns of proteins when grown in the same culture medium, and so a specific signature could be associated to each species. These patterns were useful to identify and map orthologous proteins among the six species, as well as to discriminate the closely related species T. tonsurans and T. equinum, which cannot be differentiated by ribosomal DNA. This work also presents the first in-depth identification of the major proteins secreted by dermatophytes growing under conditions promoting extracellular proteolytic activity. This catalogue of proteins, which include several endo- and exo- proteases, other hydrolases, oxydoreductases, and proteins of unknown function, probably represents the spectrum of enzymes that allow the degradation of keratinized tissues into compounds which can be assimilated by the fungus. The results suggest that ecological switching could be related to a differential expression of genes encoding secreted proteins, particularly, proteases, rather than genetic divergences of the genes encoding orthologous proteins. Differential secretion of proteins by Dermatophyte species could also be responsible for the variable inflammation caused by the infectious agent within the host. Therefore, the proteins here identified are also important to shed light into the immune response of the host during the infection process.

Outcome evaluation of ankle osteoarthritis treatments using spatio-temporal gait parameters and plantar pressure during unconstrained long distance walking : P25

Introduction: Ankle arthrodesis (AD) and total ankle replacement (TAR) are typical treatments for ankle osteoarthritis (AO). Despite clinical interest, there is a lack of their outcome evaluation using objective criteria. Gait analysis and plantar pressure assessment are appropriate to detect pathologies in orthopaedics but they are mostly used in lab with few gait cycles. In this study, we propose an ambulatory device based on inertial and plantar pressure sensors to compare the gait during long-distance trials between healthy subjects (H) and patients with AO or treated by AD and TAR. Methods: Our study included four groups: 11 patients with AO, 9 treated by TAR, 7 treated by AD and 6 control subjects. An ambulatory system (Physilog®, CH) was used for gait analysis; plantar pressure measurements were done using a portable insole (Pedar®-X, DE). The subjects were asked to walk 50 meters in two trials. Mean value and coefficient of variation of spatio-temporal gait parameters were calculated for each trial. Pressure distribution was analyzed in ten subregions of foot. All parameters were compared among the four groups using multi-level model-based statistical analysis. Results: Significant difference (p <0.05) with control was noticed for AO patients in maximum force in medial hindfoot and forefoot and in central forefoot. These differences were no longer significant in TAR and AD groups. Cadence and speed of all pathologic groups showed significant difference with control. Both treatments showed a significant improvement in double support and stance. TAR decreased variability in speed, stride length and knee ROM. Conclusions: In spite of a small sample size, this study showed that ankle function after AO treatments can be evaluated objectively based on plantar pressure and spatio-temporal gait parameters measured during unconstrained walking outside the lab. The combination of these two ambulatory techniques provides a promising way to evaluate foot function in clinics.

Emergence of Swiss metropole and scaling properties of urban clusters

Defining the limits of an urban agglomeration is essential both for fundamental and applied studies in quantitative and theoretical geography. A simple and consistent way for defining such urban clusters is important for performing different statistical analysis and comparisons. Traditionally, agglomerations are defined using a rather qualitative approach based on various statistical measures. This definition varies generally from one country to another, and the data taken into account are different. In this paper, we explore the use of the City Clustering Algorithm (CCA) for the agglomeration definition in Switzerland. This algorithm provides a systemic and easy way to define an urban area based only on population data. The CCA allows the specification of the spatial resolution for defining the urban clusters. The results from different resolutions are compared and analysed, and the effect of filtering the data investigated. Different scales and parameters allow highlighting different phenomena. The study of Zipf's law using the visual rank-size rule shows that it is valid only for some specific urban clusters, inside a narrow range of the spatial resolution of the CCA. The scale where emergence of one main cluster occurs can also be found in the analysis using Zipf's law. The study of the urban clusters at different scales using the lacunarity measure - a complementary measure to the fractal dimension - allows to highlight the change of scale at a given range.

Identification of Wax Esters in Latent Print Residues by Gas Chromatography-Mass Spectromertry and Their Potential Use as Aging Parameters

Recent studies show that the composition of fingerprint residue varies significantly from the same donor as well as between donors. This variability is a major drawback in latent print dating issues. This study aimed, therefore, at the definition of a parameter that is less variable from print to print, using a ratio of peak area of a target compound degrading over time divided by the summed area of peaks of more stable compounds also found in latent print residues.Gas chromatography-mass spectrometry (GC/MS) analysis of the initial lipid composition of latent prints identifies four main classes of compounds that can be used in the definition of an aging parameter: fatty acids, sterols, sterol precursors, and wax esters (WEs). Although the entities composing the first three groups are quite well known, those composing WEs are poorly reported. Therefore, the first step of the present work was to identify WE compounds present in latent print residues deposited by different donors. Of 29 WEs recorded in the chromatograms, seven were observed in the majority of samples.The identified WE compounds were subsequently used in the definition of ratios in combination with squalene and cholesterol to reduce the variability of the initial composition between latent print residues from different persons and more particularly from the same person. Finally, the influence of a latent print enhancement process on the initial composition was studied by analyzing traces after treatment with magnetic powder, 1,2-indanedione, and cyanoacrylate.