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.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Interaccion among systems of finite size in predictive relativistic mechanics I. General framework

We derive a Hamiltonian formulation for the three-dimensional formalism of predictive relativistic mechanics. This Hamiltonian structure is used to derive a set of dynamical equations describing the interaction among systems in perturbation theory.

Classical predictive electrodynamics of two charges with radiation: General framework. I.

Outgoing radiation is introduced in the framework of the classical predictive electrodynamics using LorentzDiracs equation as a subsidiary condition. In a perturbative scheme in the charges the first radiative self-terms of the accelerations, momentum and angular momentum of a two charge system without external field are calculated.

A Three-Dimensional Framework to Analyse the Governance of Population Registers

In June 2006, the Swiss Parliament made two important decisions with regards to public registers' governance and individuals' identification. It adopted a new law on the harmonisation of population registers in order to simplify statistical data collection and data exchange from around 4'000 decentralized registers, and it also approved the introduction of a Unique Person Identifier (UPI). The law is rather vague about the implementation of this harmonisation and even though many projects are currently being undertaken in this domain, most of them are quite technical. We believe there is a need for analysis tools and therefore we propose a conceptual framework based on three pillars (Privacy, Identity and Governance) to analyse the requirements in terms of data management for population registers.

State-to-state reaction probabilities within the quantum transition state framework

Rigorous quantum dynamics calculations of reaction rates and initial state-selected reaction probabilities of polyatomic reactions can be efficiently performed within the quantum transition state concept employing flux correlation functions and wave packet propagation utilizing the multi-configurational time-dependent Hartree approach. Here, analytical formulas and a numerical scheme extending this approach to the calculation of state-to-state reaction probabilities are presented. The formulas derived facilitate the use of three different dividing surfaces: two dividing surfaces located in the product and reactant asymptotic region facilitate full state resolution while a third dividing surface placed in the transition state region can be used to define an additional flux operator. The eigenstates of the corresponding thermal flux operator then correspond to vibrational states of the activated complex. Transforming these states to reactant and product coordinates and propagating them into the respective asymptotic region, the full scattering matrix can be obtained. To illustrate the new approach, test calculations study the D + H2(ν, j) → HD(ν′, j′) + H reaction for J = 0.

Schizotypy as an organizing framework for social and affective sciences

Schizotypy, defined in terms of commonly occurring personality traits related to the schizophrenia spectrum, has been an important construct for understanding the neurodevelopment and stress-diathesis of schizophrenia. However, as schizotypy nears its sixth decade of application, it is important to acknowledge its impressively rich literature accumulating outside of schizophrenia research. In this article, we make the case that schizotypy has considerable potential as a conceptual framework for understanding individual differences in affective and social functions beyond those directly involved in schizophrenia spectrum pathology. This case is predicated on (a) a burgeoning literature noting anomalies in a wide range of social functioning, affiliative, positive and negative emotional, expressive, and social cognitive systems, (b) practical and methodological features associated with schizotypy research that help facilitate empirical investigation, and (c) close ties to theoretical constructs of central importance to affective and social science (eg, stress diathesis, neural compensation). We highlight recent schizotypy research, ie providing insight into the nature of affective and social systems more generally. This includes current efforts to clarify the neurodevelopmental, neurobiological, and psychological underpinnings of affiliative drives, hedonic capacity, social cognition, and stress responsivity systems. Additionally, we discuss neural compensatory and resilience factors that may mitigate the expression of stress-diathesis and functional outcome, and highlight schizotypy's potential role for understanding cultural determinants of social and affective functions.

Building a Social Accounting Matrix within the ESA95 Framework: Obtaining a Dataset for Applied General Equilibrium Modelling

This research provides a description of the process followed in order to assemble a "Social Accounting Matrix" for Spain corresponding to the year 2000 (SAMSP00). As argued in the paper, this process attempts to reconcile ESA95 conventions with requirements of applied general equilibrium modelling. Particularly, problems related to the level of aggregation of net taxation data, and to the valuation system used for expressing the monetary value of input-output transactions have deserved special attention. Since the adoption of ESA95 conventions, input-output transactions have been preferably valued at basic prices, which impose additional difficulties on modellers interested in computing applied general equilibrium models. This paper addresses these difficulties by developing a procedure that allows SAM-builders to change the valuation system of input-output transactions conveniently. In addition, this procedure produces new data related to net taxation information.

Does human capital stimulate investment in physical capital?: Evidence from a cost system framework

The direct effect of human capital on economic growth has been widely analysed in the economic literature. This paper, however, focuses on its indirect effect as a stimulus for private investment in physical capital. The methodological framework used is the duality theory, estimating a cost system aggregated with human capital. Empirical evidence is given for Spain for the period 1980-2000. We provide evidence on the indirect effect of human capital in making private capital investment more attractive. Among the main explanations forthis process, we observe that higher worker skill levels enable higher returns to be extracted from investment in physical capital.

The connectome viewer toolkit: an open source framework to manage, analyze, and visualize connectomes.

Advanced neuroinformatics tools are required for methods of connectome mapping, analysis, and visualization. The inherent multi-modality of connectome datasets poses new challenges for data organization, integration, and sharing. We have designed and implemented the Connectome Viewer Toolkit - a set of free and extensible open source neuroimaging tools written in Python. The key components of the toolkit are as follows: (1) The Connectome File Format is an XML-based container format to standardize multi-modal data integration and structured metadata annotation. (2) The Connectome File Format Library enables management and sharing of connectome files. (3) The Connectome Viewer is an integrated research and development environment for visualization and analysis of multi-modal connectome data. The Connectome Viewer's plugin architecture supports extensions with network analysis packages and an interactive scripting shell, to enable easy development and community contributions. Integration with tools from the scientific Python community allows the leveraging of numerous existing libraries for powerful connectome data mining, exploration, and comparison. We demonstrate the applicability of the Connectome Viewer Toolkit using Diffusion MRI datasets processed by the Connectome Mapper. The Connectome Viewer Toolkit is available from http://www.cmtk.org/

Molecular genetic framework for protophloem formation.

The phloem performs essential systemic functions in tracheophytes, yet little is known about its molecular genetic specification. Here we show that application of the peptide ligand CLAVATA3/EMBRYO SURROUNDING REGION 45 (CLE45) specifically inhibits specification of protophloem in Arabidopsis roots by locking the sieve element precursor cell in its preceding developmental state. CLE45 treatment, as well as viable transgenic expression of a weak CLE45(G6T) variant, interferes not only with commitment to sieve element fate but also with the formative sieve element precursor cell division that creates protophloem and metaphloem cell files. However, the absence of this division appears to be a secondary effect of discontinuous sieve element files and subsequent systemically reduced auxin signaling in the root meristem. In the absence of the formative sieve element precursor cell division, metaphloem identity is seemingly adopted by the normally procambial cell file instead, pointing to possibly independent positional cues for metaphloem formation. The protophloem formation and differentiation defects in brevis radix (brx) and octopus (ops) mutants are similar to those observed in transgenic seedlings with increased CLE45 activity and can be rescued by loss of function of a putative CLE45 receptor, BARELY ANY MERISTEM 3 (BAM3). Conversely, a dominant gain-of-function ops allele or mild OPS dosage increase suppresses brx defects and confers CLE45 resistance. Thus, our data suggest that delicate quantitative interplay between the opposing activities of BAM3-mediated CLE45 signals and OPS-dependent signals determines cellular commitment to protophloem sieve element fate, with OPS acting as a positive, quantitative master regulator of phloem fate.

A competency framework in cardiothoracic surgery for training and revalidation - an international comparison.

The conventional methods of education, certification and recertification in cardiothoracic surgery face a paradigm shift in line with recent innovations in diagnostics and therapeutics. The attributes of a competent clinician entail proficiency in knowledge, communication, teamwork, management, health advocacy, professionalism and technical skills. This article investigates the skills required for a cardiothoracic surgeon to be competent. The relevant practice of certification and recertification across various regions has also been explored. Validated and competency-based curricula should be designed to develop core competencies to successfully integrate them into practice. Challenges to the implementation of such curricula and potential solutions are explored. Patient safety remains the ultimate aim to ensure excellence of both competency and performance.