Effects of the Shape of Distribution of Mutant Effect in Nearly Neutral Mutation Models

Total Page:16

File Type:pdf, Size:1020Kb

Effects of the Shape of Distribution of Mutant Effect in Nearly Neutral Mutation Models J. Genet., Vol. 75, Number 1, April I996, pp. 33-48. (() Indian Academy of Sciences Effects of the shape of distribution of mutant effect in nearly neutral mutation models HIDENORI TACHIDA Department of Biology, Faculty of Science, Kyushu University 33, Fukuoka 812, Japan Abstract. Models of the theory of nearly neutral mutation incorporate a continuous distribu- tion of mutation effects in contrast to the theory of purely neutral mutation which allows no mutations with intermediate effects. Previous studies of one such model, namely the house-of- cards mutation model, assumed normal distribution of mutation effect. Here I study the house-of-cards mutation model in random-mating finite populations using the weak-mutation approximation, paying attention to the effects of the distribution of mutant effects. The average selection coefficient, substitution rate and average heterozygosity in the equilibrium and transient states were studied mainly by computer simulation. The main findings are: (i) Very rapid decrease of the substitution rate and very slow approach to equilibrium as selection becomes stronger are characteristics of assuming normal distribution of mutant effect. If the right tail of the mutation distribution decays more rapidly than that of the normal distribution, the decrease of substitution rate becomes slower and equilibrium is achieved more quickly. (ii) The dispersion index becomes smaller or larger than 1 depending on the time and the intensity of selection. (iii) Let N be the population size. When selection is strong the ratio of 4N times the substitution rate to the average heterozygosity, which is expected to be 1 under neutrality, is larger than 1 in earlier generations but becomes less than 1 in later generations. These findings show the importance of the distribution of mutant effect and time in determina- tion of the behaviour of various statistics frequently used in the study of molecular evolution. Keywords. Protein evolution; substitution rate; average heterozygosity; population genetics; dispersion index. 1. Introduction Understanding the mechanism of molecular evolution is one of the most important problems in biology. Kimura (1968) proposed the neutral theory, which states that the main cause of evolutionary change at the molecular level is random fixation of selectively neutral or very nearly neutral mutations rather than Darwinian selection. In its pure form (Kimura 1987 and later), the neutral theory assumes that there are two discrete classes of mutations depending on their selection coefficient, s. One class is neutral mutations whose selection coefficients satisfy N~ rsl << 1, where N~ is the effective population size. The other class is deleterious mutations whose selection coefficients satisfy -Nes >> 1, and these do not contribute to molecular evolution. The theory allows no mutations with intermediate effects. The proposal aroused a lot of contro- versy (see Lewontin 1974 and Kimura 1983 for review) but now the theory is widely accepted as an explanation for most of DNA changes in noncoding regions and synonymous sites. However, the mechanism for amino acid substitutions (replacement substitutions) is still controversial (Gillespie 1991). Ohta (1972) proposed the nearly neutral mutation theory to explain replacement substitutions. In this theory mutations with intermediate effects, i.e. Nes ~- O(1), are permitted and the distribution of mutation effect is considered somewhat continuous 33 34 Hidenori Tachida (see Ohta 1992a for review). The first model that explicitly incorporated a continuous distribution of mutation effects was that by Ohta (1977). In this model an exponential distribution was employed as a distribution of mutation effects. More precisely, the difference between the fitnesses of the mutated and original alMes was assumed to be distributed exponentially. Such models with fixed distributions for differences between the mutated and the original alleles are in a class called the shift model (Ohta and Tachida 1990). In the shift model the existence of advantageous mutations was corlsidered unlikely because, if they exist, the substitution rate becomes extremely large in large populations, which is not observed in molecular data (Kimura 1979). Ohta and Tachida (1990) introduced another type of nearly neutral mutation model, called the fixed model, in which the fitness of the mutant allele, not the difference between the fitnesses of the mutated and the original alleles as in the shift model, has a fixed distribution. This mutation model is the same as the house-of-cards model of Kingman (1978), and we call the model by this name from now on. The most distinguishing feature of this model from the shift model is that the distribution of the fitness difference between the mutated and original alMes changes as the population evolves. Thus, when the average fitness of the population is low, most mutations are advantageous; when the average fitness of the population is high most mutations become deleterious. Even though advantageous mutations are allowed in this model, the substitution rate does not become large in large populations because such populations have high average fitnesses and most mutations are deleterious. This is in accord with the traditional view of mutation, namely that mutations are mostly deleterious after evolution. Ohta and Tachida (1990) and Ohta (1992b) considered the behaviour of this mutation model in geographically structured populations while Tachida (1991) considered a simpler case of random-mating populations. Tachida's (1991) conclusions are, in summary: (i) The behaviour is mainly determined by the product c~ = 4Neo', where o" is the standard deviation of the mutant distribution. (ii) If 0~ is less than 0.2, the behaviour is almost the same as that of the neutral model. (iii) If c~ is between 0.2 and 3-5, both selection and random genetic drift contribute to the gene evolution. Although the long-term substitution rate is reduced as c~ increases, it is difficult to distinguish its behaviour from that of the neutral model using conventional neutrality tests such as that using the average and the variance of heterozygosity. (iv) If ~ is larger than 5, selection dominates the process. The average fitness quickly goes up and after that very few substitutions occur. However, it takes a long time to reach equilibrium. In all these studies, the distribution of mutant fitness is assumed to be normal. However, there is no biological reason for assuming normality of mutant fitness distribution. Thus it is necessary to investigate the effects of shape of distribution on the evolution. In the present paper I investigate the house-of-cards model in finite panmictic populations employing various mutant distributions. Specifically, I investigate the effects of the distribution shape on average fitness, substitution rate, heterozygosity and dispersion index. To save computer time, i use a jump-process approximation that can be used when the mutation rate is small. I first describe the model and explain this approximation, and then summarize results obtained by this approximation. House-of-cards mutation model 35 Appropriateness of the approximation was checked by comparison with the result obtained by the simulation study of Tachida (1991) on the Wright-Fisher model. 2. Analysis We assume the diploid Wright-Fisher model with fixed population size N (see Crow and Kimura 1970). Consider a locus with an infinite number of alleles (the infinite-allele model of Kimura and Crow 1964). Mutation always results in a new allele, labelled by a random variable S. We call S the selection coefficient of the allele. The fitness of an individual that has two alleles labelled by S x and S 2 is 1 + S 1 + S 2. For different alleles, S is independently distributed with the same density, fo(s). We designate the variance of it by o-z. This mutation model is the house-of-cards model of Kingman (1978). In a finite population the frequency changes become a multidimensional Markov process and it is difficult to deal directly with this process (but see Ethier and Kurtz 1994). So we use an approximation that is valid for small mutation rates. 2.1 Weak-mutation approximation Let u be the mutation rate per generation. IfNu is much smaller than 1 the population is expected to be mostly monomorphic, and occasionally fixations of new alleles occur very quickly. This can be approximated by a one-dimensional Markovian jump process (Zeng et al. 1989; Tachida 1991). We briefly explain this jump process. In this process the population is characterized by the selection coeff• of the allele currently fixed. Let Po (s; t) be the density of the population currently fixed by an allele with the selection coefficient s at time t. A transition from a state s to another state with a selection coefficient between r and r + dr in a very short time At occurs with a probability 2(r - s) 2NuAt ~e (r) dr, (1) 1 - exp [- 4N(r - s)] a~ because 2Nufo (r)drAt mutations within this range occur and the probability of fixation for such mutations is 2(r-s)/(1-exp[-4N(r--s)~) (Kimura 1962). Thus, if we measure time in units of 1/u generations, the transition equation for Po(S; t) is written as (Tachida 1991) dPo(S,t) = 2N { I ~ 2(s-r) d----7-- ~ 1 - exp [- 4N(s -r)j f~176 t)dr 2(r - s) t) dr }. -fT~l_exp[_4N(r_s)]fo(r)p(s , _ (2) If we scale s and r by the standard deviation o" offo (s), the equation can be written as dp(s, t) { f ~ (s-r) d~ - ~ ~o 1 - exp [- ~(s -r)] f(s)p(r' t)dr - f7 oo 1 - exp[--~(r(r-s) _s)]f(r)p(s,t)dr } ' (3) 36 Hidenori Tachida where a = 4No- and p(s, t) andf(s) are densities of the population and mutation effect respectively for scaled variables. Thus the process is characterized by a and the shape of the standardized mutation density f(s).
Recommended publications
  • 701.Full.Pdf
    THE AVERAGE NUMBER OF GENERATIONS UNTIL EXTINCTION OF AN INDIVIDUAL MUTANT GENE IN A FINITE POPULATION MOT00 KIMURA' AND TOMOKO OHTA Department of Biology, Princeton University and National Institute of Genetics, Mishima, Japan Received June 13, 1969 AS pointed out by FISHER(1930), a majority of mutant genes which appear in natural populations are lost by chance within a small number of genera- tions. For example, if the mutant gene is selectively neutral, the probability is about 0.79 that it is lost from the population during the first 7 generations. With one percent selective advantage, this probability becomes about 0.78, namely, it changes very little. In general, the probability of loss in early generations due to random sampling of gametes is very high. The question which naturally follows is how long does it take, on the average, for a single mutant gene to become lost from the population, if we exclude the cases in which it is eventually fixed (established) in the population. In the present paper, we will derive some approximation formulas which are useful to answer this question, based on the theory of KIMURAand OHTA(1969). Also, we will report the results of Monte Carlo experiments performed to check the validity of the approximation formulas. APPROXIMATION FORMULAS BASED ON DIFFUSION MODELS Let us consider a diploid population, and denote by N and Ne, respectively, its actual and effective sizes. The following treatment is based on the diffusion models of population genetics (cf. KIMURA1964). As shown by KIMURAand OHTA (1969), if p is the initial frequency of the mutant gene, the average number of generations until loss of the mutant gene (excluding the cases of its eventual fixation) is - In this formula, 1 On leave from the National Institute of Genetics, Mishima, Japan.
    [Show full text]
  • National Institute of Genetics Japan
    NATIONAL INSTITUTE OF GENETICS JAPAN ANNUAL REPORT No. 20 1969 Publi8l&ed bll THE NATIONAL INS'ITI'UTE OF GENETICS Mitlima. Simoka-ken. Japtm 1970 Annual Report of the National Institute of Genetics No. 20, 1969 Published by The National Institute of Genetics, Japan 1970 CONTENTS General statement 1 Staff..................................... .. 2 Council............................................................ 5 Association for propagation of the knowledge of genetics... 5 Projects of research for 1969 ................ 6 Researches carried out in 1969 .......................... 9 1. Cytogenetics Segregation of three types of the largest No.1 (A-I) chromosome in Rattus rattus bred in a population room. YOSIDA, T. H. and TSUCHIYA, K. 9 Polymorphism of No.9 (A-9) and No. 13 (A-13) autosomes in Rattus rattus collected in Japan. YOSIDA, T. H. and TSUCHIYA, K. and MASUJI, H.......................................... 10 Frequency of chromosome polymorphism in Rattus rattus col- lected in East and Southeast Asia. YOSIDA, T. H. and TSUCHIYA, K. 11 Fl and F2 hybrids between Asian and Oceanian black rats with different karyotypes. YOSIDA, T. H., KATO, H. and SAKAKI- BARA, K.................................................... 13 Karyotypes of black rats, Rattus rattus, obtained from Honolulu, Hawaii. YOSIDA, T. H. and TSUCHIYA, K. 14 Local differences in the frequencies of chromosome polymorphism in the house rat, Rattus rattus, in Japan. YOSIDA, T. H., TSUCHIYA, K. and SONODA, J. 15 Karyotypes of Rattus bowersii collected in Malaysia. YOSIDA, T. H. and TSUCHIYA, K. 16 Protective effects of protamine sulfate on disintegration of isolated metaphase chromosomes. KATO, H., SEKIYA, K. and YOSIDA, T. H. 17 A preliminary study on the uptake of isolated nuclei by mam­ malian cells in vitro.
    [Show full text]
  • Manipulating Underdetermination in Scientific Controversy: the Case of the Molecular Clock
    Dartmouth College Dartmouth Digital Commons Open Dartmouth: Published works by Dartmouth faculty Faculty Work 9-1-2007 Manipulating Underdetermination in Scientific Controversy: The Case of the Molecular Clock Michael Dietrich Dartmouth College Robert A. Skipper University of Cincinnati Follow this and additional works at: https://digitalcommons.dartmouth.edu/facoa Part of the Biology Commons Dartmouth Digital Commons Citation Dietrich, Michael and Skipper, Robert A., "Manipulating Underdetermination in Scientific Controversy: The Case of the Molecular Clock" (2007). Open Dartmouth: Published works by Dartmouth faculty. 16. https://digitalcommons.dartmouth.edu/facoa/16 This Article is brought to you for free and open access by the Faculty Work at Dartmouth Digital Commons. It has been accepted for inclusion in Open Dartmouth: Published works by Dartmouth faculty by an authorized administrator of Dartmouth Digital Commons. For more information, please contact [email protected]. Manipulating Underdetermination in Scientiªc Controversy: The Case of the Molecular Clock Michael R. Dietrich Dartmouth College Robert A. Skipper, Jr. University of Cincinnati Where there are cases of underdetermination in scientiªc controversies, such as the case of the molecular clock, scientists may direct the course and terms of dispute by playing off the multidimensional framework of theory evaluation. This is because assessment strategies themselves are underdetermined. Within the framework of assessment, there are a variety of trade-offs between differ- ent strategies as well as shifting emphases as speciªc strategies are given more or less weight in assessment situations. When a strategy is underdetermined, scientists can change the dynamics of a controversy by making assessments using different combinations of evaluation strategies and/or weighting what- ever strategies are in play in different ways.
    [Show full text]
  • The Average Number of Generations Until Fixation of a Mutant Gene in a Finite Population'
    THE AVERAGE NUMBER OF GENERATIONS UNTIL FIXATION OF A MUTANT GENE IN A FINITE POPULATION' MOT00 KIMURA AND TOMOKO OHTA, National Institute of Genetics, Mishima, Japan Received July 26, 1968 A mutant gene which appeared in a finite population will eventually either be lost from the population or fixed (established) in it. The mean time until either of these alternative events takes place was studied by WATTERSON(1962) and EWENS(1963). They made use of a method previously announced by DARL- ING and SIEGERT(1953), and, independently by FELLER(1954). Actually, DARL- ING and SIEGERTrefer to its application to genetics. From the standpoint of population genetics, however, it is much more desirable to determine separately the mean time until fixation and that until loss. Since the gene substitution in a population plays a key role in the evolution of the species, it may be of particular interest to know the mean time for a rare mutant gene to become fixed in a finite population, excluding the cases in which such a gene is lost from the population. In the present paper, a solution to this problem will be presented together with Monte Carlo experiments to test some of the theoretical results. Throughout this paper, the senior author (M. K.) is responsible for the mathematical treatments, while the junior author (T. 0.) is responsible for the numerical treatments based on computers. BASIC THEORY Let us consider a diploid population consisting of N individuals and having the variance effective number Ne, which may be different from the actual number (for the definition of Ne, see KIMURAand CROW1963).
    [Show full text]
  • National Institute of Genetics Japan
    NATIONAL INSTITUTE OF GENETICS JAPAN ANNUAL REPORT No. 19 1968 Published bll THE NATIONAL INSTITUTE OF GENETICS Misima. Sizuoka·ken, Japan 1969 Annual Report of the National Institute of Genetics No. 19, 1968 Published by The National Institute of Genetics, Japan 1969 CONTENTS General statement. ................................................. 1 Research member.................................................. 2 Council............................................................ 5 Projects of research for 1968....................................... 6 Researches carried out in 1968 ..................................... 9 I. Cytogenetics Scientific expedition for the study of rodents to South East Asia and Oceania. I. Members, aims and schedule of the scienti- fic expedition. YOSIDA, T. H. 9 Scientific expedition for the study of rodents to South East Asia and Oceania. II. Species and chromosome numbers of rodents collected from South East Asia and Oceania. YOSIDA, T. H., TSUCHlYA, K. and IMAI, H. T......................... 10 Scientific expedition for the study of rodents to South East Asia and Oceania. III. Chromosomal polymorphism and new karyotypes of black rat, Rattus rattus, collected in South East Asia and Oceania. YOSIDA, T. H. and TSUCHIYA, K......... 11 Scientific expedition for the study of rodents to South East Asia and Oceania. IV. Comparative idiogram analysis in several species belonging to genus Rattus collected from South East Asia and Oceania. YOSIDA, T. H. and TSUCHlYA, K......... 12 Scientific expedition for the study of rodents to South East Asia and Oceania. V. Breeding of wild rodents collected from South East Asia and Oceania. TSUCHIYA, K., YOSIDA, T. H... 13 Rat metaphase chromosomes incorporated into mouse metaphase cells. YOSIDA, T. H. and SEKIGUCHl, T...................... 15 Autoradiographic analysis of the rate of protein synthesis in diploid and tetraploid MSPC-l myeloma cells.
    [Show full text]
  • The Neoselectionist Theory of Genome Evolution
    The neoselectionist theory of genome evolution Giorgio Bernardi* Molecular Evolution Laboratory, Stazione Zoologica Anton Dohrn, Villa Comunale, 80121 Naples, Italy Edited by Tomoko Ohta, National Institute of Genetics, Mishima, Japan, and approved April 4, 2007 (received for review February 22, 2007) The vertebrate genome is a mosaic of GC-poor and GC-rich isoch- of chance in evolution. A very significant modification of the ores, megabase-sized DNA regions of fairly homogeneous base neutral theory was the nearly neutral theory of Ohta (11, 12), composition that differ in relative amount, gene density, gene who proposed that a substantial fraction of changes are caused expression, replication timing, and recombination frequency. At by random fixation of nearly neutral changes, a class that the emergence of warm-blooded vertebrates, the gene-rich, mod- ‘‘includes intermediates between neutral and advantageous, as erately GC-rich isochores of the cold-blooded ancestors underwent well as between neutral and deleterious classes’’ (13). The first a GC increase. This increase was similar in mammals and birds and four schemes of Fig. 1 display a qualitative picture of the classical was maintained during the evolution of mammalian and avian theories just mentioned. orders. Neither the GC increase nor its conservation can be ac- All classical theories proposed that natural selection acted on counted for by the random fixation of neutral or nearly neutral the ‘‘phenotype’’ [called the ‘‘classical phenotype’’ by Bernardi single-nucleotide changes (i.e., the vast majority of nucleotide and Bernardi (14), to distinguish it from the ‘‘genome pheno- substitutions) or by a biased gene conversion process occurring at type’’ presented below].
    [Show full text]
  • National Institute of Genetics Japan
    ISSN 0077-4995 NATIONAL INSTITUTE OF GENETICS JAPAN ANNUAL REPORT No. 40 1989 Published by THE NATIONAL INSTITUTE OF GENETICS Mi.ima, Sizuoka-ken, Japan 1990 Annual Report of the National Institute of Genetics No. 40, 1989 Published by The National Institute of Genetics, Japan 1990 CONTENTS General statement .........................................I Staff 3 Council and Advisory Committee. ................................ 6 Association for propagation of the knowledge of genetics 8 Projects of research for 1989 9 Research activities in 1989 14 I. Molecular Genetics Promoter selectivity of Escherichia coli RNA polymerase, I. Effect of base substitutions in promoter - 35 region on the promoter strength. KOBAYASHI, M., NAGATA, K. and ISHIHAMA, A. ... 14 The promoter selectivity of Escherichia coli RNA polymerase, II. Random screening of promoters and classification based on the promoter strength. KUBOTA, M., YAMAZAKI, Y. and ISHIHAMA, A. 15 Promoter selectivity of Escherichia coli RNA polymerase, III. Structural and functional modulations of RNA polymerase during growth phase transition. OZAKI, M., FUJITA, N., WADA, A. and ISHIHAMA, A. ..... 16 Structural modulation of ribosomes during growth phase transition of Escherichia coli. WADA, A., YAMAZAKI, Y., FUJITA, N. and ISHIHAMA, A. 17 The molecular anatomy of Escherichia coli RNA polymerase: Func­ tional mapping of the a subunit. IGARASHI, K., HAYWARD, R. S., FUJITA, N., YAMAGISHI, M. and ISHIHAMA, A. 18 Promoter selectivity of Micrococcus luteus RNA polymerase: Identification and characterization of the major sigma factor. NAKAYAMA, M., FUJITA, N., OSAWA, S. and ISHIHAMA, A..... 19 Molecular mechanisms of transcription and replication of phage Qfi: Cloning, sequencing and mapping of host factor. KAJI- TANI, M. and ISHIHAMA, A.
    [Show full text]
  • The Crafoord Prize in Biosciences 2015
    THE CRAFOORD PRIZE IN BIOSCIENCES 2015 POPULAR SCIENCE BACKGROUND The Crafoord Prize in Biosciences 2015 The Crafoord Prize in Biosciences 2015 is awarded to the geneticists Richard Lewontin, USA, and Tomoko Ohta, Japan, for their contributions to knowledge of genetic variation. Thanks to Lewontin and Ohta, we know that the genetic variation within populations of plants and animals is many times greater than science once thought – and we understand why. The late twentieth century and early twenty-first century have seen a revolution in genetic research. New techniques for DNA sequencing have resulted in the rapid growth of knowledge in many areas of biology, including ecology, evolution and systematics. The influence of genetics can also be seen in popular culture, where TV series’ crime scene investigators search for traces of DNA on a hunt for the killer. These traces are able to reveal who were at the scene of the crime because each person’s genotype is as unique as his or her fingerprints. This doesn’t only apply to humans, but to all animals and plants that reproduce sexually; there is considerable genetic variation between individuals of the same species and in the same population. This variation is so significant that we can be certain that two individuals do not have identical DNA even if they are closely related. The only exception to this is identical (monozygotic) siblings. Knowledge of the great genetic variation found within populations and each individual’s unique genotype is now a well-integrated element of our understanding of the world, as well as fundamental to science.
    [Show full text]
  • Random Genetic Drift & Gene Fixation
    Random Genetic Drift & Gene Fixation Arie Zackay July 31, 2007 Abstract Random Genetic Drift is one of the evolutionary forces that effects the distribution of alleles and changes their frequencies in populations. The main idea behind this theory is that genetic changes can happen randomly. This can be explained through the Bottleneck effect. There are two main factors that have an influence on genetic drift, the first one is the Effective Population Size Ne , which creates different genetic drift outcomes within different population sizes. The second factor that influences the rate of genetic drift is the selective advantage of a drifted allele. With these fac- tors it is possible to calculate the time and probability of an allele fixation in any chosen population. 1 Introduction In order to understand the phenomenon known to us as Random Genetic Drift, we must first of all define the meaning of Population Genetics. Population Ge- netics is the study of the changes in distribution of alleles in a certain population, considering the influence of the following evolutionary forces: natural selection, genetic drift, mutation, gene flows and selective mating. The main goals of this field of study is to determine how the frequency of a mutant gene will change over time under the influence of evolutionary forces, to examine how genetic variability is maintained, and to establish the probability that a new mutant completely replaces existing variants of genes, a phenomenon called Fixation. There are a few instigators that can cause changes in the evolution of a popu- lation. Such instigators might include: • It can be caused by a mutation with a genetic advance.
    [Show full text]
  • Dobzhansky's Evolution of Tropical Populations, and the Science and Politics Of
    CARVALHO, Tito. “A most bountiful source of inspiration:” Dobzhansky’s evolution of tropical populations, and the science and politics of genetic variation. História, Ciências, Saúde – Manguinhos, Rio de Janeiro, v.26, n.1, “A most bountiful jan.-mar. 2019, p.281-297. Abstract source of inspiration:” Theodosius Dobzhansky has been studied for how he integrated Dobzhansky’s evolution of field naturalism and laboratory experimentation in ways that helped tropical populations, and produce the Modern Synthesis, as well as how he leveraged biological expertise to support liberal and cosmopolitan the science and politics of values amidst Second World War and the Cold War. Moreover, Dobzhansky genetic variation has been central in analyses of the institutionalization of genetics in Brazil, where he spent several years. This article “A mais abundante fonte situates Dobzhansky’s Brazilian research within the science of variation and the de inspiração”: Dobzhansky politics of diversity. I conclude by raising questions about how the ways in which e sua evolução sobre as science figured in politics depended on ideas about the role of scientists in society whichwere advanced in populações dos trópicos, parallel, suggesting research on the “co- a ciência e a política da production” of natural and social orders. Keywords: evolutionary genetics; variabilidade genética transnational science; eugenics; race; tropics; Theodosius Dobzhansky (1900- 1975). Resumo Theodosius Dobzhansky tem sido estudado pelo modo como ele integrou o naturalismo de campo e a experimentação científica, que deram origem à síntese moderna, assim como a alavanca que ele deu ao conhecimento biológico para apoiar valores liberais e cosmopolitas em meio à Segunda Guerra Mundial e à Guerra Fria.
    [Show full text]
  • Evolutionary Population Genetics of Promoters: Predicting Binding Sites and Functional Phylogenies
    Evolutionary population genetics of promoters: Predicting binding sites and functional phylogenies Ville Mustonen and Michael La¨ ssig* Institut fu¨r Theoretische Physik, Universita¨t zu Ko¨ln, Zu¨lpicherstrasse 77, 50937 Cologne, Germany Edited by Tomoko Ohta, National Institute of Genetics, Mishima, Japan, and approved August 8, 2005 (received for review June 30, 2005) We study the evolution of transcription factor-binding sites in neutral background DNA. As an important step in this direction, prokaryotes, using an empirically grounded model with point the notion of a fitness landscape for binding-site sequences has been mutations and genetic drift. Selection acts on the site sequence via introduced, where the fitness of a site depends on the binding its binding affinity to the corresponding transcription factor. Cal- energy of the corresponding factor (8). The evolutionary impor- ibrating the model with populations of functional binding sites, we tance of the binding energy has also been highlighted in ref. 9, where verify this form of selection and show that typical sites are under it was shown that nucleotide substitution rates within functional substantial selection pressure for functionality: for cAMP response sites in Escherichia coli depend on the energy difference induced by protein sites in Escherichia coli, the product of fitness difference the substitution as predicted from the position weight matrix. The and effective population size takes values 2N⌬F of order 10. We biophysics of factor-DNA binding imposes stringent constraints on apply this model to cross-species comparisons of binding sites in the form of the fitness landscape (10) and has important conse- bacteria and obtain a prediction method for binding sites that uses quences for bioinformatic binding site searches (11).
    [Show full text]
  • Population Genetics, a Concise Guide
    Population Genetics vdxiaovd Population Genetics A Concise Guide John H.vvvspie THEJOHNS HOPKINS UNIVERSITYPRESS Baltimore and London ._./,.,.I.,.,_,_,,.,,... ...l,.,ll.”_.~....,.,.,., ,.*..I .....I..., ..... ,..I.......,.. ,’. , , ’ 0 1998 The JohnsHopkine University Press All rights reserved. Published 1998 Printed in the United Statesof America on acid-free paper 9876543 The JohnsHopkins University Press 2715 North Charles Street Baltimore, Maryland 21218-4363 www.press.jhu.edu Library of Congress Cataloging-in-Publication Data will be'found at the endof this book. A catalog record for this book is available from the British Library. ISBN 0-8018-5764-6 ISBN 0-8018-5755-4(pbk.) To Robin Gordon Contents List of Figures ix Preface xi 1 The Hardy-Weinberg Law 1 1.1 DNA variationin Drosophila .................... 2 1.2 Loci and alleles ............................ 5 1.3 Genotypeand allele frequencies ................... 9 1.4 Randomly matingpopulations .................... 11 1.5 Answers to problems ......................... 17 2 Genetic Drift 19 2.1 A computersimulation ........................ 20 2.2 The decay of heterozygosity ..................... 22 2.3 Mutationand drift .......................... 27 2.4 Theneutral theory .......................... 32 2.5 Effective population size ....................... 35 2.6 The coalescent ............................ 38 2.7 Binomial sampling .......................... 42 2.8 Answers to problems ......................... 47 Natural Selection 3 Natural 49 3.1 The fundamental model ......................
    [Show full text]