Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

The evolution of the genetic load caused by recurrent mutation in small populations : genetic context… Poon, Arthur F. 2000

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata

Download

Media
831-ubc_2000-0530.pdf [ 3.78MB ]
Metadata
JSON: 831-1.0089617.json
JSON-LD: 831-1.0089617-ld.json
RDF/XML (Pretty): 831-1.0089617-rdf.xml
RDF/JSON: 831-1.0089617-rdf.json
Turtle: 831-1.0089617-turtle.txt
N-Triples: 831-1.0089617-rdf-ntriples.txt
Original Record: 831-1.0089617-source.json
Full Text
831-1.0089617-fulltext.txt
Citation
831-1.0089617.ris

Full Text

THE E V O L U T I O N OF T H E GENETIC L O A D C A U S E D B Y R E C U R R E N T MUTATION IN S M A L L POPULATIONS: GENETIC CONTEXT A N D DEMOGRAPHIC HISTORY by ART POON Hon. B . Sc., The University o f Toronto, 1998. A THESIS S U B M I T T E D IN P A R T I A L F U L F I L M E N T OF THE REQUIREMENTS FOR T H E D E G R E E OF M A S T E R OF SCIENCE in T H E F A C U L T Y OF G R A D U A T E STUDIES (Department o f Zoology) We accept this thesis as conforming to the required standard  THE UNIVERSITY OF BRITISH C O L U M B I A July 2000 © A r t Poon, 2000  UBC  6/7/2000 10:13 A M  Special Collections - Thesis Authorisation Form  In p r e s e n t i n g t h i s t h e s i s i n p a r t i a l f u l f i l m e n t o f t h e requirements f o r an a d v a n c e d d e g r e e a t t h e U n i v e r s i t y o f B r i t i s h C o l u m b i a , I a g r e e t h a t t h e L i b r a r y s h a l l make i t f r e e l y a v a i l a b l e f o r r e f e r e n c e and s t u d y . I f u r t h e r a g r e e t h a t p e r m i s s i o n f o r e x t e n s i v e c o p y i n g o f t h i s t h e s i s f o r s c h o l a r l y p u r p o s e s may be g r a n t e d by t h e head o f my d e p a r t m e n t o r by h i s o r h e r r e p r e s e n t a t i v e s . I t i s u n d e r s t o o d t h a t copying o r p u b l i c a t i o n o f t h i s t h e s i s f o r f i n a n c i a l gain s h a l l not be a l l o w e d w i t h o u t my w r i t t e n p e r m i s s i o n .  Department o f The U n i v e r s i t y o f B r i t i s h V a n c o u v e r , Canada  Columbia  Date  http://www.library.ubc.ca/spcoll/thesauth.html  Page 1 of 1  ABSTRACT  I have conducted two separate theoretical investigations on the consequences o f recurrent mutation for the extinction risk o f small populations. The continual generation o f deleterious mutations and their subsequent rise in frequency in the population cause genetic load: a reduction in the mean fitness o f the population. This evolutionary process is a function o f natural selection, mutation, and genetic drift. Predictions made by theoretical studies on the extinction risk depend on what assumptions are made about selection, mutation, and drift. The following work modifies standard assumptions that are made in models that study genetic loads i n finite populations such that the sensitivity o f the results can be evaluated. First o f all, many models assume that most mutations are unconditionally deleterious — i.e. they have the same effect in any genome. There may be, however, a large fraction o f mutations that conceal the expression o f another mutation when placed in the same genome. These are compensatory mutations. With compensatory mutation, qualitatively different conclusions are reached about the extinction risk o f a small population by genetic causes.  Secondly, the models that study the effect o f decreasing population size  on the genetic load simplify the demographic change to an instantaneous drop from infinite size. Declining populations often gradually decline i n size over several generations, however. I have studied the effect o f gradual demographic change by numerical methods. The results indicate that transient changes in the mean load caused by recessive alleles can be delayed by a gradual change in population size.  ii  TABLE OF CONTENTS  Abstract  ii  List o f Figures  iv  Acknowledgements CHAPTER I  v  The genetic load o f small populations  1  Selection, mutation, and extinction Evolutionary consequences o f small population size  1 4  Genetic loads  7  Statement o f authorship  11  C H A P T E R II  12  The effect o f compensatory mutation  2.1  Preface  12  2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10  Overview Evidence for compensatory mutations Fisher's geometrical model o f adaptation Moment analysis o f the master equation Distribution functions and numerical methods Results in the one-dimensional model Results in the ^-dimensional model Sensitivity o f the model to its assumptions Implications for the extinction risk o f small populations  12 14 17 21 28 34 39 43 49  2.11 Implications for pleiotropy C H A P T E R III 3.1 3.2 3.3 3.4 3.5 3.6 CHAPTER  51  The effect o f demographic history  54  Overview  54  Steady-state expectations o f the genetic load Marginal fitness and purging Numerical methods Numerical results Concluding remarks I V Conclusion and recommendations for further work  54 57 59 63 71 77  Bibliography  79  iii  LIST OF FIGURES  Chapter II: 2-1 2-2 2-3 2-4 2-5 2-6 2- 7  Fisher's geometrical model o f adaptation Approximation o f the fixation probability for new mutant alleles Mutant vectors in two- and three-dimensional space M e a n fixed drift load against effective population size for n - 1 Mean fixed drift load against A, for n = 1 M e a n fixed drift load against effective population size for n = 10 Mean fixed drift load against n  18 25 31 36 38 41 42  Chapter III: 3- 1 The glide and shift trajectories 61 3-2 M e a n genetic load ( x l O ) o f glide trajectories plotted against time in generations 65 3-3 Mean genetic load ( x l O ) o f glide trajectories at generation 10, plotted against s 66 3-4 M e a n genetic load ( x l O ) o f shift trajectories against time i n generations. . . 69 3-5 M e a n genetic load ( x l O ) o f shift trajectories at generation 10, plotted against s 70 -6  -6  -6  -6  'iv  ACKNOWLEDGEMENTS  I thank my advisor, S. Otto, for her patience through the development o f my thesis and for making possible my transition to graduate-level mathematics; and the rest o f my committee, M . Whitlock and R. Redfield, for their guidance and insight. I also thank the numerous professors that I have approached with theoretical problems including M . Doebeli, P. Greenwood, S. Marion, M . Barlow, B . Bergerson, and H . A . Orr; and doctoral students T. Duty and T. Johnson. For their help in programming and other computational issues, I thank A . Blachford and J. Hauser o f the Zoology Computing Unit, and M . Choptuik for use o f the Physics Department L i n u x Cluster. This thesis was supported by N S E R C P G S - A scholarship number 208606-1998 and an N S E R C grant to S. Otto.  C H A P T E R I: OVERVIEW  1.1 Selection, mutation, and extinction The mandate o f conservation biology is the preservation o f biological diversity, in part, by reducing the extinction risk o f natural populations. Although every population w i l l eventually go extinct, there are intrinsic (e.g. genetic) and extrinsic (e.g. environmental) factors that decrease the expected time to extinction. Evolutionary biology can contribute to this effort by characterizing the genetic factors that contribute to the extinction risk. Before discussing the relevance o f evolutionary processes to the extinction risk, it is necessary to first explain the demographic process o f population growth and decline. The expected time to extinction is determined by the population's ability to maintain its numbers. To model the growth and decline o f a population, it is useful to assume that its generations are discrete such that there is no mating between generations because the equations become much simpler (Hastings 1997). Consider a population o f size N. Let the expected (i.e. average) number o f offspring o f the z individual that survive to adulthood be represented by <W,>. It th  follows that the expected number o f individuals comprising the next generation, <N'>, is:  {N') = i(W ).  (1.1.1)  t  A useful quantity for describing the growth o f populations is derived from the ratio o f population sizes between generations:  1  AT  N  N  ,  v  —r  1=1  where I have used (1.1.1) to show that this ratio is equal to the expected number o f surviving offspring averaged over all individuals in the parental generation. A population i n which R is less than one, such that the number o f offspring that w i l l survive to adulthood is less than the number o f parents, w i l l soon go extinct. Furthermore, i f R only slightly exceeds one, then a substantial extinction risk remains because the actual (or realized) population size, TV, may by chance be less than expected over a number o f generations. Such variation in R reflects the cumulative effect o f fluctuations i n the realized number o f offspring, W , over all t  individuals in the population. The fluctuations around each average, <W,>, are either caused by chance (i.e. demographic stochasticity) or by environmental changes (Lande 1988). The mean time to extinction is maximized when R is considerably greater than one, when the carrying capacity o f the habitat is very large, and when there is little variation over time in R.  The summary variable R retains no information about individual variation i n W nor h  about the variation o f the mean values, <W,>, among individuals. A s we shall see, the latter is necessary for studying the evolutionary dynamics o f R. In fact, W is considered to be synonymous with Darwinian fitness (Crow and K i m u r a 1970) and reflects the performance of an individual i n its environment. Similarly, R is referred to as the "mean absolute fitness" of a population. The number o f descendants that an individual w i l l produce on average, <Wj>, is determined by how well-adapted it is to its environment. Variation among individuals i n the degree o f adaptation, and hence in the expected number o f descendants left, is the sine qua non o f natural selection. Yet there is no response to natural selection, and thus  2  evolution, unless fitness variation in the population is inherited by subsequent generations. The principal mechanism o f biological inheritance is genetic.  Therefore, i f variation in  <W,> among individuals has a heritable genetic basis, then there w i l l be evolution by natural selection. Indeed, genetic material that, when expressed, tends to produce a phenotype (i.e. characteristic or form) that outperforms any other can be expected to be found i n a progressively greater proportion o f the population over time. A s a result, the mean fitness, R, of the population w i l l often increase except when it is at the carrying capacity o f its habitat. It follows that selection on variation in fitness can reduce the risk o f extinction.  Genetic variance is constantly being introduced into a population by recurrent mutation. A mutation is a change in the genetic sequence that is brought about by the inherent inaccuracy o f our genomic replication machinery, as well as other genetic processes including transposable genetic elements (Nuzhdin and Mackay 1995). Every biological organism exhibits some measurable rate o f mutation (Drake et al. 1999). Although any mutation in a gene sequence in a strict sense creates a different version (i. e. allele) o f that gene, I am mainly concerned with those that modify the performance o f the subsequent gene product. Because o f the inherent complexity o f biological systems, predicting the effect that a particular mutation w i l l have on fitness is a difficult task. However, both the results o f experimental studies on mutation (Mukai 1964; Fry et al. 1999) and our understanding of how genes and proteins work imply that most mutations that affect fitness reduce individual fitness (i.e. are deleterious). It follows that an individual carrying a large number o f mutations, relative to the mean number i n the population, w i l l on average contribute fewer offspring to the next generation, reducing the transmission probability o f those mutant  3  alleles. This variance among individuals in <W,> introduced by mutation causes a net reduction in R. The magnitude o f this reduction is determined by the number o f deleterious mutant alleles in the population (Haldane 1937). This number reflects a balance between the generation o f mutations and their elimination by selection. This dynamic equilibrium is often referred to as the mutation-selection balance. Ih a small population, the mean sojourn time of a new deleterious mutant allele (i.e. the mean number o f generations until its descendants are eliminated from the population) is also affected by genetic drift. The details o f genetic drift w i l l be clarified in the next section. This phenomenon is particularly relevant to problems in conservation biology, as the populations that are studied often have been reduced to small numbers (Caughley 1994).  1.2 Evolutionary consequences of small population size  If the offspring are equivalent with respect to juvenile mortality, then the subset that survive to become the next adult generation w i l l be a random sample o f those offspring. There is some chance that any particular parent w i l l be over- or under-represented in this sample o f offspring. Therefore the expected number o f offspring, <W>, is a different quantity than the realized number, W. Random fluctuations i n the relative success o f genetically distinct groups in a population vary the expected pattern o f inheritance and thus cause genetic drift. When the size o f this sample between generations is small, these fluctuations increase in magnitude. Because the sample is by definition independent o f the biological character o f each individual, genetic drift is a process that is independent o f natural selection. Put another way, an offspring that is genetically well-adapted to its  4  environment has no guarantee o f surviving to reproductive maturity because o f accidents that occur over time.  The variance o f allele frequencies i n the sample o f offspring caused by genetic drift is a decreasing function o f the population size, N. However, this relationship makes several assumptions about the reproductive character o f the population, such as equal sex ratio and random mating. I f these assumptions are violated, then drift in the population w i l l behave as it would in an ideal population o f another size in which the assumptions hold true. The size of this ideal population is the effective population size. Whereas the original quantity, called the census population size, is the countable number o f individuals, the latter is an evolutionary quantity that describes the effect o f drift on the genetic state o f that population (Crow and K i m u r a 1970). When the sex ratio is unequal, for example, the mean number o f gametes contributed to the next generation is no longer equal among individuals. The effective population size consequently becomes smaller than the census size (Crow and K i m u r a 1970). Indeed, empirical estimates o f the effective size o f populations is roughly one-tenth the census size on average (Frankham 1995).  Genetic drift can interfere with natural selection by causing individuals with many deleterious alleles to become over-represented in subsequent generations and as a result cause a net reduction in R.  There is, in fact, some probability that sampling error can  eventually cause every individual i n the population to have descended from a common ancestor that carried the most deleterious mutant alleles i n its generation. For a haploid asexual population, this process is analogous to Muller's ratchet (Muller 1964; Felsenstein  5  1974), in which the class of individuals with the least number o f mutations goes extinct by chance. It is now necessary for our purposes, however, to concentrate on the number o f gene copies i n the population rather than the number o f individuals. When the population has more than one allele at a locus (i.e. a polymorphic locus), it is more convenient to collect gene copies together according to their allelic form (i. e. identity by state) and characterize the genetic state o f a population by the frequencies o f each allele. There is a chance that drift w i l l eventually cause a deleterious allele to become the only allele in the population. This state is called fixation. The average number o f generations required for a new mutant allele to be fixed in a population is determined by the intensity o f genetic drift, which i n turn is dependent on the effective population size over time.  Kondrashov (1995) inferred from the available empirical data on the rate and severity of deleterious mutation that the accumulation o f slightly deleterious mutant alleles in the genome should be sufficient to cause the extinction o f most populations. This, o f course, assumes that the accumulation o f mutant alleles corresponds to an accumulation o f unconditionally deleterious effects on fitness. O n the contrary, it is possible that new mutant alleles that cancel out the deleterious effects o f previous mutations are frequently generated. Mutations that are not required to be in combination with other genes to increase fitness are called beneficial. If a mutation w i l l only increase fitness in the presence o f certain deleterious alleles, it is called a compensatory mutation. In the second chapter, I w i l l investigate the effect o f compensatory mutation on the extinction risk o f small populations that accumulate many deleterious mutations. Moreover, the effects o f genetic drift are often studied in a population held at a constant size (e.g. K i m u r a 1955; Lande 1994). Endangered  6  populations, however, are often in decline or have only recently declined to their current level because o f human-mediated disturbance (Ehrlich and Ehrlich 1981; W i l s o n 1988). Very few studies have investigated how the demographic history o f a population affects its genetic state and subsequently its extinction risk (Kirkpatrick and Jarne 2000). This is the subject o f the third chapter. Hence, the following two chapters represent further contributions to an existing literature on the application o f evolutionary theory to the problem of extinction risk (e.g. L y n c h and Gabriel 1990; Lande 1988, 1994; Burger and L y n c h 1995; Frankham 1995; Hedrick and M i l l e r 1992). A common feature o f these studies is the genetic load, which acts as an index describing the genetic component o f a population's extinction risk. The next section o f this chapter w i l l be devoted to identifying the various types o f genetic load that are directly caused by recurrent mutation in finite populations.  1.3 Genetic load Genetic variation is the raw material from which all evolutionary adaptations are derived. A lack o f genetic variation may even contribute to the extinction risk o f a population, as adaptation to a changing environment is required for survival (Burger and L y n c h 1995). In addition, there is some evidence that populations with a lack o f genetic variation are less resistant to disease and parasitism (Frankham 1995), although we do not know whether the level o f resistance depends on genetic variation itself or reflects the fixation o f deleterious alleles. Nevertheless, many features o f our genetic system that maintain genetic variation are unavoidably associated with a cost, paid as a reduction in the mean fitness o f the population (Crow 1970). We refer to each o f these fitness reductions as a "genetic load," following Muller's (1950) use o f the term "load" to describe the deleterious  7  impact o f mutation on populations. Load is measured relative to an ideal population that lacks the genetic feature o f interest. Mathematically, this may be written:  where Ro is the mean absolute fitness o f the ideal population (Crow 1970). Although mutation generates genetic variation, for example, most new mutant alleles that affect fitness have deleterious effects on fitness (Mukai 1964). Thus, recurrent mutation causes, on average, a reduction i n fitness that is called the mutation load (Muller 1950; Morton et al. 1956) and is measured relative to the mean fitness o f a population that does not experience mutation. If the effects o f drift are negligible, as is the case in an infinitely-large population, then selection w i l l act against mutation to establish an equilibrium frequency for deleterious alleles (Haldane 1937), such that mutant alleles are being removed from the population as quickly as they are being generated. In this thesis, I w i l l limit the term "mutation load" to refer to the expected load in a large population where drift is absent. Mutation load, then, is assumed to be present in a population o f any given size, in addition to which other genetic loads w i l l cause further changes i n the mean fitness.  In finite populations, genetic drift w i l l move the frequencies o f mutant alleles away from the equilibria determined by mutation and selection, by causing an allele to be over- or under-represented in subsequent generations. This in turn can create an additional change in fitness called the drift load. The drift load can be apportioned into two components, each associated with a different feature o f the probability distribution o f allele frequencies. First  8  of all, there is a portion o f the drift load that is caused by asymmetries in the spread o f the frequencies o f deleterious unfixed alleles about their deterministic equilibria, due to genetic drift. W e refer to the difference between the expected mutation load and the total reduction in fitness caused by variation in the frequency o f polymorphic genes as the segregating drift load. The segregating drift load can contribute the greater fraction o f the drift load in moderately large populations, in which the rate o f fixation is low and allelic polymorphism is maintained. A s population size decreases, however, allele frequency distributions become more U-shaped (Wright 1937) because new mutant alleles are rapidly eliminated or fixed in the population so that few loci are polymorphic with alleles that affect fitness. Deleterious alleles that are fixed i n a population contribute to the second component o f the drift load, which I w i l l refer to as the fixed drift load. Because the fixation probability o f deleterious mutations increases exponentially with decreasing population size (Crow and K i m u r a 1970), fixed mutations make a larger contribution to the drift load in smaller populations. Furthermore, fixed mutations accumulate over time, which can contribute to an incremental reduction in fitness and population size known as the mutational meltdown process (Gabriel et al. 1993; L y n c h and Gabriel 1990). The accumulation o f fixed drift load may eventually exceed both the mutation and segregating drift loads in small populations (Lande 1994). Therefore, it is reasonable and useful simplification to focus on the fixed drift load to assess the genetic contribution to the extinction risk o f a small population. The model presented in Chapter II takes advantage o f this simplification to model the effect o f compensatory mutations on the genetic load. Because the mutation load does not change much as a function o f N, it acts as a baseline from which further reductions in mean fitness occur. The segregating drift load is assumed to be small relative to the fixed drift load when the  9  population is held at a small size. O n the other hand, chapter III presents numerical results that do not require such simplifications so the genetic load is treated as a single quantity including mutation load and drift load, both segregating and fixed. Both chapters demonstrate how the expected genetic load o f a small population is affected by changes in standard assumptions o f previous evolutionary models.  10  S T A T E M E N T OF AUTHORSHIP  The following chapter is adapted from a co-authored manuscript to appear in the journal Evolution. I appear as main author and my supervisor, S. Otto, is senior (second) author. I was responsible for the majority o f the analysis and for preparing the manuscript. S. Otto was responsible for the initial solution o f the one-dimensional analysis.  We are in agreement that the contributions o f the thesis author are as stated above:  Senior Author  Thesis Author  11  Chapter II: THE EFFECTS OFCOMPENSATORY  MUTATION  2.1 Preface This chapter appears in the journal Evolution under the title "Compensating for our load o f mutations: freezing the meltdown o f small populations" (Poon and Otto, in press), as a manuscript co-authored with my advisor, S. Otto. The contributions o f each author to the manuscript are described on a separate document that precedes this chapter, as per the requirements o f the Special Collections and University Archives Division. I was responsible for the majority o f the analysis and writing. Redundancies i n the manuscript i n the context o f this thesis have been removed.  2.2 Overview Do small populations, in which there is a substantial probability that slightly deleterious mutations w i l l fix, have a high risk o f going extinct as a result? The majority of theoretical studies that have addressed this question affirm that many asexual and some sexual populations with l o w effective population sizes (see Section 1.2) w i l l inevitably go extinct through mutation accumulation (Butcher 1995; Gabriel et al. 1993; Lande 1994, 1998; L y n c h et al. 1993, 1995a, 1995b; L y n c h and Gabriel 1990; Schultz and L y n c h 1997). Even populations as large as our own species, given what we have measured o f mutation, may be in genetic peril (Crow 1997; Eyre-Walker and Keightley 1999; Kondrashov 1995). Thus it appears that deleterious mutation presents a genuine hazard to the conservation o f species and the well-being o f human populations. However, this threat may have been  12  overestimated as an artifact o f how the effect o f mutation on the mean fitness o f populations has been modeled. The earliest work, such as Haldane's (1937) pioneering investigations into the effect o f variation in fitness in an infinite population, simplified the biology o f mutations by ignoring the phenotypic origin o f their effects on fitness (Hartl and Taubes 1998). Subsequent investigators have done much to improve the biological realism within this classical framework by introducing rates o f reversal, epistasis, and effect variation (Butcher 1995; Kondrashov 1994; Lande 1998; Schultz and L y n c h 1997). Nevertheless, most previous investigations have assumed that mutations can be classified as either beneficial or deleterious regardless o f the current phenotype. Such mutations are said to be unconditionally beneficial or deleterious (Wagner and Gabriel 1990).  The effect o f a mutation at biochemical, physiological, morphological and behavioral levels may, however, depend on its genomic context (Wright 1968). In other words, a mutation may be deleterious in some genomes and advantageous i n others. Interactions within and among gene products, furthermore, provide opportunities for one mutation to compensate for the effects o f another at some level o f the phenotype. Hence, new mutations may restore fitness losses incurred by previous mutations without requiring true reversals i n the gene sequence (Burch and Chao 1999; Hartl and Taubes 1996; K i m u r a 1990; Ohta 1992; Wagner and Gabriel 1990). These are known as compensatory or suppressor mutations. There is an abundance o f experimental evidence, reviewed below, for mutations that conceal maladaptive mutant phenotypes at molecular, biochemical, and organismal levels. With few exceptions (e.g. Hartl and Taubes 1998), we do not yet know how compensatory mutations might influence predictions made by evolutionary models because context-dependence is not  13  a part o f the classical framework. Therefore it would be useful to determine how changing the way we model mutations affects the outcome o f models addressing biological problems, such as the extinction risk o f small populations. To do so requires a framework that defines a mutation by its effect on the phenotype, subsequently treating the fitness effect as a result o f the phenotypic alteration. Fisher's geometrical model o f adaptation (Fisher 1930), which I w i l l describe in Section 2.4, is a good candidate that has been applied i n a handful o f studies to evolutionary models o f mutation (Haiti and Taubes 1998; Orr 1998; Wagner and Gabriel 1990). One study in particular (Hartl and Taubes 1998) has derived an analytical estimate for the long-term mean fitness reduction caused by fixed mutations, which is called the fixed drift load (see Section 1.3). For reasons discussed previously, this is an important quantity regarding the extinction risk o f small populations. I have conducted an analysis o f a model similar to Hartl and Taubes' (1998), representing the evolution o f a finite population under mutation, selection, and drift. In this analysis, a probability distribution describing the amount by which a population needs to adapt is summarized by its moments (e.g. its central tendency and spread; Kendall and Stuart 1963). Furthermore, I have supplemented the results o f the analysis with computer simulations to evaluate the precision o f approximations made in the analysis. Differences between Hartl and Taubes' (1998) analysis and the model presented here have substantial effects on fixed drift load predictions, which I w i l l investigate at the conclusion o f this chapter.  2.3 Evidence for Compensatory Mutations There are two types o f empirical evidence that support the postulate that compensatory mutations are a common biological phenomenon. First, there is an extensive  14  body o f experimental work identifying "suppressor" mutations — an older term synonymous with compensatory mutation — that conceal the expression o f a mutant gene (Hartman and Roth 1973; Jarvik and Botstein 1975). These mutations may occur within the mutant gene itself or i n other genes. Compensation by additional mutations in the mutant gene (i.e. intragenic compensation) is often accomplished by completely or partially restoring the original structural or functional conformation o f the gene product (e.g. H o u and Schimmel 1992; Mateu and Fersht 1999; K i m et al. 1994; Hanson et al. 1993). A recent study o f human hemoglobin ( K i m et al. 1994) provides a good example o f this process. The authors modeled the energetic folding o f a mutant hemoglobin protein and identified a second basepair mutation that best recovered normal quaternary structure. That these mutations were compensatory was then verified empirically by site-directed mutagenesis. Compensatory intragenic mutations are particularly well understood i n the case o f R N A structures ( r R N A , Clark et al. 1984, Hancock et al. 1988, Springer et al. 1995; m R N A , Stephan and Kirby 1993; t R N A , H o u and Schimmel 1992). For example, mutations in stem structures o f t R N A s can be compensated for by a matching mutation in the opposite strand (Cedergren et al. 1981; Steinberg and Cedergren 1994). Compensation by mutations i n genes other than that containing the first mutation — which w i l l herein be referred to as the target gene — (i.e. extragenic or intergenic compensation) is often mediated by some biochemical relationship between gene products. For example, a compensatory mutation in a gene involved in translation can conceal the expression o f the target gene by causing occasional errors in the translation o f the target gene's m R N A transcript, such that the subsequent gene product is wildtype (Hartman and Roth 1973; Waterston and Brenner 1978; Murgola 1985; Y u and Spreitzer 1992; E l Mezaine et al. 1998); frequently, these compensatory mutations turn out  15  to be point-mutations in certain t R N A genes. Alternatively, compensation can be mediated through biochemical pathways (Hartman and Roth 1973). When a mutant enzyme does not function sufficiently, the build-up o f unmetabolized substrate or the lack o f metabolized product may have deleterious effects. Subsequent mutations that shunt flux towards an alternate pathway (Maringanti and Imlay 1999; Dickinson et al. 1995; Gachotte et al. 1997), block flux into the impaired pathway (Manning et al. 1999), or deactivate negative extragenic regulators (Matsuno and Sonenshein 1999) may compensate for this blockage.  A second source o f evidence for compensatory mutation comes from studies in which direct measures o f fitness compensation are made.  Many classical studies find that  deleterious mutations kept in laboratory populations for extended periods have decreasing fitness effects over time (Haldane 1957; Muller 1938). Furthermore, there are more recent experimental studies i n which populations with reduced fitness due to the effects o f a particular allele recover through compensatory effects at other loci (Burch and Chao 1999; Elena et al. 1998; Liang et al. 1998). It is worth discussing the results from Burch and Chao (1999) because they are particularly relevant to our model. Burch and Chao (1999) fixed a highly deleterious mutation in a strain o f the bacterial virus 06 and then propagated several  descendant lineages, varying the effective population size among the lineages by allowing only a small sample to propagate (i.e. bottlenecking). The mean fitness o f these bottlenecked lineages recovered over time at a rate inversely proportional to their effective population size. Lineages that recovered the fitness loss in one step were presumed to carry the original allele that was regenerated by a reverse mutation i n the deleterious allele. Other lineages that recovered fitness in a gradual step-wise fashion were inferred to have fixed compensatory  16  mutations at other sites i n the genome. Burch and Chao's (1999) experiment provides evidence that compensatory mutations are more common than a reverse mutation, and therefore that the reverse mutations are relatively rare in small populations that have fewer genomes in which to accumulate mutations. It follows that small populations may depend more heavily on compensatory mutations for fitness recovery.  2.4 Fisher's Geometrical Model of Adaptation Fisher's (1930) geometrical model o f adaptation provides a framework in which unconditionally deleterious and compensatory mutations can be studied together. This geometrical model was originally devised by R. A . Fisher to promote a micromutational view of adaptation, in which adaptation is achieved by successive substitutions o f mutations with small effects on the phenotype (Leigh 1987). Fisher's model has made several recent appearances in the literature (Burch and Chao 1999; Hartl and Taubes 1998; Orr 1998; Peck et al. 1997). Its original premise is that an organism's phenotype can be broken down into some number o f components, n, that evolve independently o f one another. I f each component is represented by an axis, then an ^-dimensional rectangular coordinate space (i.e. «-space) can be constructed that represents the full range o f possible phenotypes. We let the scale o f each axis be equal to the fitness reduction caused by a given change i n phenotype from the optimum. Suppose that a given population is genetically monomorphic and that the fixation o f mutant alleles occurs instantaneously. Let a point i n the «-space labeled A represent the current average phenotype o f the population. If there exists a single optimallyadapted phenotype, represented by a point i n the «-space labeled O (see F i g . 2-1), then the  17  Figure 2-1. - Fisher's geometrical model o f adaptation. The model for two dimensions (n = 2) is shown. The population state is represented by point A that resides at some distance \AO\ from the optimum O, which is set at the origin for convenience. The fitness o f a population is a function o f \AO\ with a maximum dX\AO\-Q. Mutations are represented by vectors in the ndimensional character space. (/) Mutations that increase \AO\ are deleterious. (//') Conversely, those decreasing \AO\ are advantageous, (iii) Mutations between points on the circle So do not change \AO\ and are neutral, (iv) The probability that a mutation w i l l be advantageous is a decreasing function o f its length. Mutation vectors longer than the diameter o f So cannot be advantageous. The effect o f a mutation on fitness depends on its context. A mutation has an advantageous effect on one side o f the optimum (ii) and a deleterious effect on the other (v).  18  population's degree o f adaptiveness is some function o f the Euclidean distance between these two points:  where a, and o, are the coordinate projections along the z' axis such that A = (a\, ai, th  and 0 = (o\,02,  ct ) n  o ). For convenience, we place the optimum at the origin such that o, = n  0 for all /. It is useful to use a spherical coordinate system instead o f the Cartesian coordinate system, such that \AO\ becomes the radius o f a hypersphere, So, that is centered at O. Because every point on the surface o f So is equidistant from the origin, each corresponds to a phenotype with the same fitness value.  A mutation is represented by a vector in this «-space that originates from the point A and terminates at some other point A' (Fig. 2-1). The length o f a mutation vector corresponds to the amount o f phenotypic change that it has caused along some combination o f components. Mutations o f a fixed length w i l l terminate on the surface o f a hypersphere SA centered at A; this sphere w i l l be useful for explaining the derivation o f the probability distribution functions that characterize the mutation process in this system. The fitness effect associated with a mutation depends on where A' is located relative to the surface o f So (see Fig. 2-1, i-iii). Without this contextual information, one cannot determine the fitness effects of a particular mutation even i f one knows its size and orientation i n Fisher's model. Mutations between points on the surface o f So w i l l have a neutral effect on fitness because there is no change in distance to the origin, although each brings about a different change in phenotype. A mutation that brings a population inside So reduces the distance to the origin  19  and thus has an advantageous effect on fitness. Conversely, a mutation that moves the population outside So has a deleterious effect. It is useful to note that a correspondence exists between the size o f mutations and their expected effect on fitness. Small mutations are more likely to be favorable because the local curvature o f So becomes flatter, such that the probability o f being advantageous approaches 1 / 2 (Fisher 1930). Although the probability that a mutation is favorable decreases with increasing mutation size (see F i g . 2-1, iv), advantageous mutations o f larger effect are more likely to be incorporated into an adaptive walk because they confer a greater fitness advantage than their slightly advantageous counterparts and are therefore less likely to be lost by drift (Kimura 1983). When mutations become too large, however, this effect disappears because large mutations tend to overshoot the optimum and the mean fitness advantage decreases. Mutations larger than the diameter o f So are never advantageous. Unlike most classical population genetic models, the genetic context o f a mutation defines its effect on fitness in Fisher's model. For example, we can visualize a vector with some fixed direction which is adaptive on one side o f the origin and maladaptive on another (see Fig. 2-1, ii and v). Thus, in Fisher's model there are no mutations that are intrinsically and unconditionally deleterious or advantageous. It is this quality o f the model that makes it useful to our purposes.  The number o f orthogonal axes, n, in Fisher's model corresponds to the number o f phenotypic aspects by which an organism must adapt (Fisher 1930) assuming that only one optimal phenotype exists. Consequently, the model makes certain postulates about the number o f ways that an organism can have a suboptimal fitness. Specifically, there are 2"~  l  more ways to have a fitness o f l - 2 e than l-e (e > 0). This becomes evident i f one considers  20  that the surface area o f an ^-dimensional hypersphere increases at a rate proportional to the (fl-l)-th power o f its radius. Therefore n, to some extent, quantifies the rarity o f optimal versus nearly optimal states. To better understand this we make use o f Fisher's (1930) original analogy o f a microscope, except that it is now focused with n dials. The accumulation o f error i n the tuning o f each dial results in a large number o f alternate settings of the dials that could all lead to the same degree o f poor focus.  2.5 Moment analysis of the master equation A s a population fixes mutations, it w i l l wander in «-space at varying distances to the optimum. W e use z' to denote the Euclidean distance o f the population from the optimum. Throughout this analysis, we w i l l assume that the mean environmental effect is zero and that the population size is small enough that the segregating drift load is negligibly small relative to the fixed drift load. The state o f a population can be represented by a probability distribution p such that p(z', t) is the probability that the population state is within the interval [z, z'+dz') at time t. We can describe the dynamics o f this probability distribution over time with the following differential equation:  dt  (Kimura 1965), where z is a state variable similar to z' that is used to integrate over the entire range o f previous states, <U/is the genomic rate o f fixation for new alleles, and g(z', z) is the  transition probability from state z to z'. The first term represents the rate that the population  21  state leaves the interval [z', z'+dz') upon mutation, and the second term integrates the cumulative probability that the population state w i l l enter the interval [z', z'+dz') from another interval [z, z+dz). W e assume that the occurrence o f mutant alleles is rare enough and/or that recombination is frequent enough that the fixation o f a mutation is not interfered with by the simultaneous segregation o f other alleles. If p(z', t) converges on the same distribution i n the long term regardless o f its initial state, then a steady state exists for this differential equation. The steady state is very useful because it describes the probability distribution i n the long term regardless o f the initial conditions. Because we are primarily interested in finding a steady state we can ignore the details o f the approach to a steady state which allows us to redefine the time scale through a change o f variables:  T — t • fAf  where T measures time i n units based on the average waiting time between fixation events. The differential equation then becomes:  (2.5.1)  which implies that a new mutation fixes with every time step dr. To use this equation requires that we define the transition probability density function g(z', z). Because each transition represents the fixation of a new allele, g{z', z) depends on the probability density function o f mutations from z to z' — which we w i l l denote as m(z', z) — weighted by the  22  fixation probability u(s) o f the new allele, where s is the selection coefficient that depends on  z and z:  g(z',z) =  m(z',z)u(s)  (2.5.2)  The denominator o f (2.5.2) is a normalizing constant such that J  g(z',z)dz' = 1. Note that  the probability that a mutation moves the population from z ' t o z, m{z', z), w i l l depend on the number o f orthogonal axes, n. For example, a mutation w i l l have a greater chance o f bringing a population closer to the origin o f a plane (n = 2) than the origin o f a threedimensional volume (n - 3), because a larger proportion o f mutations w i l l lie within the hypersphere So- This dependence on n w i l l be elaborated upon below, where simulation methods are presented.  ,  The resulting phenotypic change, once mapped to a change in the Euclidean distance from the origin, can be translated via the fitness function into a change in fitness. Because we are only concerned with changes o f Euclidean distance, all calculations from this point are independent o f the number o f dimensions, n, in the model. Given a fitness function w, we can calculate the selection coefficient:  s=  w(z') - w(z) w(z)  (2.5.3)  23  Note that s is a relative quantity, such that the same change i n fitness causes far stronger selection when fitness is low (w(z) «  1) than when it is high (w(z) ~ 1). The fixation  probability o f a mutation i n a diploid population with an effective size N is a function o f s: e  l-e~  2s  u  ( ) = :  1  (2.5.4)  =4177  s  —  e  (Crow and K i m u r a 1970). This function is well approximated in the neighborhood o f s = 0 (see Fig. 2-2) by:  1  u(s;s<0) = — e~ 2N  2NM  u(s;s>0) = 2s + —e~ ' 2N 2N  s  To simplify the analysis, I choose a linear function to describe absolute fitness:  w(z') =  l - B  (2.5.5)  where z' = B defines the point beyond which fitness is zero and where the optimal fitness was set, without loss o f generality, to one. Because  dw(z') dz'  (2.5.6)  B  24  F i g u r e 2-2. - Approximation o f the fixation probability for new mutant alleles. The exact and approximate fixation probabilities are plotted against the selection coefficient, s. The dashed line corresponds to the exact formula (Eq. 2.5.4; C r o w and K i m u r a 1970) for effective population size JV= 10. Our approximation, represented by the solid line, provides a good fit i n the vicinity o f s = 0, and overestimates the probability o f fixation in the region \s\ > 0.05.  25  the intensity o f selection equals the reciprocal o f the upper bound.  We are interested in finding a steady state o f the differential equation (2.5.1), so we set the time-derivative on the left-hand side to zero, giving:  0 = -p (z') + [j(z,z) (z)dz eq  (2.5.7)  Peq  wherep^z') represents the probability that the population state w i l l be i n the range [z' z'+dz') at steady state. Although it is difficult to solve (2.5.7) forp q(z^) directly, we can e  characterize the distribution, Peqiz*), by its moments. The n moment o f a distribution h(x) is th  given by the integral:  E[x ] = j:jx-a) f(x)dx n  (2.5.8)  n  where a is usually set to 0 to produce the non-central moments (Kendall and Stuart 1963). I proceed by using the method o f moment generating functions (Burger 1991) which takes advantage o f the expansion o f the exponential:  26  e =l+—+— +— +. 1! 2! 3! x  (Abramowitz and Stegun 1965) together with (2.5.8) to decompose functions into their moments. A n y density function h(x) o f an absolutely continuous random variable x can be transformed by integration into a moment generating function h(k) = \ e h(x)dx, where k is a real number. The moment generating function can be expanded to yield all moments o f the distribution described by h(x), and it also uniquely determines h(x). A p p l y i n g the method of moment generating functions to p (z'): eq  0=  j e ' (z')dz'+ ki  Peq  J J  e 'g(z',z) (z)dzdz', h  Peq  we produce  0=  \  g(k\z)p {z)dz eq  If we assume that the population is not too far from the optimum at steady state, we can take a Taylor series expansion o f g(k;z) around z = 0 to produce a polynomial o f z that we integrate to give the moments of p : eq  Vo W + Yx(k)m + y (k)m +...) 0 = - & , ( * ) + (• x  2  2  27  (2.5.9)  where m\ is the i moment o f the distribution p (z) and yfjc) is the j eq  coefficient o f the Taylor  series expansion o f g(k;z) around z = 0. Taking derivatives o f (2.5.9) with respect to k gives us a system o f equations that we limit to terms o f 0(z ) and less: 2  -m + (y L + Y[ L™, + 72 L m ) = 0 x  0  0  2  (2.5.10)  Simulations indicate that the higher moments are small and may be safely neglected in most cases. A l l analytical results and approximations were compared to simulation results. It remains to determine what form the mutational distribution function m(z', z) assumes in this analysis. I w i l l discuss this while presenting details o f the simulation methods i n the next section.  2.6 Distribution functions and numerical methods The evolution o f a population by the fixation o f new mutations is a simple Markovian stochastic process. Hence, no further information is required to determine the probability that a random variable assumes a particular state when the state at some previous moment i n time is given (Gillespie 1992). It is relatively simple to simulate a M a r k o v process with a Monte Carlo simulation program. I developed Monte Carlo simulations o f the onedimensional and ^-dimensional models i n order to evaluate their analysis. Uniformlydistributed random numbers over the interval (0, 1) were generated with shuffling to avoid low-order correlations (ran2 from Press 1992). These were i n turn used to generate random deviates o f other probability distributions by transformation whenever possible, or else by the  28  rejection method (Press 1992). In each cycle, a mutation vector was generated with some length r and orientation 6, each being drawn at random from probability distributions that differed between the one- and ^-dimensional models. In the one-dimensional simulations, mutation lengths were drawn from an exponential distribution with mean 1 / A:  (2.6.1)  I chose (2.6.1) for its relative mathematical simplicity (Mukai et al. 1972; Ohta 1977). Although the general consensus is that the distribution o f mutant effects on fitness is probably L-shaped (Keightley 1994; Ohta 1998), there is no simple correspondence between the phenotypic and fitness effects o f a mutation, and we cannot be certain what distribution is most representative. With n=\,  mutations were oriented towards or away from the optimum  with equal probability. That is, the direction o f the mutation (6>i) was set to +1 half o f the time and -1 otherwise. A mutation vector with length r less than z terminates on the same side o f the optimum. If a mutation vector is oriented towards the optimum but is so large that it overshoots the optimum, then the resulting z is the absolute difference between its length r andz. Hence,  (2.6.2)  z = Z + 9r . x  29  The selection coefficient (2.5.3) corresponding to a shift from z to z'then determined the probability o f fixation for this mutation. If the mutation was fixed, z became z ' a s the population changed state, and the process began another cycle.  When evolution takes place in more than one dimension, each mutation can always be described by a triangle with one vertex at the origin, such that the two adjacent sides become radii z ' a n d z and the remaining side becomes the mutation vector (Fig. 2-3; see also Hartl and Taubes 1998). Thus, regardless o f how many dimensions in) there are, the process can be summarized i n two. B y this scheme, any mutant vector is completely characterized by its length r and orientation 6, defined as the angle between the mutant vector and the distance vector from the optimum (Fig. 2-3). To determine the probability distribution o f r in the ndimensional model, I made use o f the fact that all mutation vectors o f equal length terminated on a hypersphere SA- Recall that according to the one-dimensional model, the probability that one o f these vectors is o f length r is described by (2.6.1). A l l that is needed then is to extrapolate (2.6.1) to account for all mutation vectors o f length r (a "bottom-up" procedure). The surface area o f a hypersphere increases at a rate proportional to r ~ with respect to the n  radius r. Using this relation we weight (2.6.1) by r "  Xe~ r'  -1  X  and normalize to get:  lr  \;Xe' r -'dr Xr  n  (2.6.3)  T(n)  30  r  F i g u r e 2-3. - Mutant vectors in two- and three-dimensional space. A n y mutation can be completely described by a triangle o f vectors regardless o f the number o f dimensions n. z represents the distance from the optimum before mutation. Each mutation is a vector o f some length r and angle 9 relative to z. Determining the next distance to the optimum, z', is then a simple matter o f trigonometric relations. A number o f vectors share the same values r and 9, which can be quantified by the sphere with radius r sm9 projected onto the plane orthogonal to the distance vector z. This quantity can then be translated into a probability function o f 9 (see text).  31  This is the gamma probability density function (Press 1992). For n>\,  (2.6.3) produces a  bell-shaped curve with a mean mutation length equal to n I A. I w i l l discuss the relative merits o f bottom-up against alternative methods o f modeling mutations i n Section 2.9. For large n and small A, implementing (2.6.3) i n a simulation becomes very inefficient because there is an excess o f unrealistically large mutations with a near-zero probability o f fixation. Consequently, I chose to truncate the distribution by the same upper bound as that defining the fitness function, B , assuming that mutations o f length greater than B were very unlikely to fix. Rescaling (2.6.3) accordingly gives us an incomplete gamma probability density function (Press 1992) for mutation length:  AV V A  (2.6.4)  B  A s n increases and A decreases, a greater proportion o f the original bell-shaped curve is truncated until only a monotonically-increasing left tail remains.  Having described the distribution o f mutation lengths, we turn now to a description o f the probability distribution o f 6. Increasing n expands the number o f positions associated  with mutation vectors o f the same length r and angle 6, which i n turn modifies the  probability distribution o f 6. To illustrate, consider that when n = 2 there are only two  vectors associated with each value o f 6 (0 < 6 <ri). When n = 3, there is a cone o f vectors  32  sharing an angle 6 around radius z, tracing a circular path on a plane orthogonal to the original position vector, z (see Fig. 2-3). Similarly, with increasing n this path becomes a hypersphere o f n - 1 dimensions. We can take advantage o f this trend to calculate the probability distribution o f 0 as a function o f n, apportioning the surface area o f the mutant ndimensional hypersphere SA (see Fig. 1) into (n - l)-dimensional hyperspheres associated with values o f 6. Because the radius of the hypersphere orthogonal to radius z is equal to r sin 6, the probability distribution function for 9 is:  f(9) = ~7n  rn (2 rsin 0)  n-2  (2.6.5)  Jo ^c'sin^fl  where 0 < 6 < n (Hartl and Taubes 1998; Rice 1990) and where c' is a normalizing constant equal to  {* sin- ddd 2  Jo  V^m(n-l)]'  where T denotes the gamma function (Abramowitz and Stegun 1965). Equation (2.6.5) produces a distribution similar to a Gaussian function that becomes narrower with increasing n. I used a rejection algorithm to generate pseudo-random numbers from this distribution (Press 1992). This method becomes less efficient, however, as the distribution becomes  33  narrower, causing the simulations to take longer when n is large. For large n (n > 10), (2.6.5) is well-approximated by a transformed Gaussian function centered at n 12 (Hartl and Taubes 1996; Fisher 1930), which implies that most mutations point i n a direction perpendicular to the direction o f the optimum. We have used this approximation to generate random deviates of 6 when n is large to reduce computational time. Given a mutation vector (r, 6) and the current distance from the optimum z, we can calculate the next distance z' from the optimum by the classic L a w o f Cosines (see Fig. 2-3):  (2.6.6)  Again, the mutation is fixed with probability u(s), a function o f the selection coefficient associated with a phenotypic change from z to z' (Eq. 2.5.3). This completes the derivation of the transition probability g(z', z) in the ^-dimensional model.  2.7 Results in the one-dimensional model  Here I w i l l describe the one-dimensional analysis as a means o f introduction to the problem and as a segue into the general n-dimensional analysis discussed below. Recall that the transition probability g(z', z) is generated from the probability density function for mutational effects on the phenotype, m(z', z), weighted by the fixation probability, u{s). The first step is to determine the probability density function, m(z', z), based on an exponential distribution o f mutation lengths. Because distance to the optimum {i.e. z or z) is an absolute quantity, every transformation from z to z' is associated with two points on the real line on  34  either side o f the origin. Assuming exponentially distributed mutation lengths and treating advantageous and deleterious mutations separately, we have:  M\z\-\z'\)  -M\z\+\z'\)  M\z'\-\z\)  -A(|z'|+|z|)  mi  (14 N; M<N) = f(«"  mi  (14 14 k1>N)=2 ~' (e  ) [advantageous]  ) [deleterious]  Partitioning the denominator o f (2.5.2) into two integrals over z>  z and z' < z and  substituting the above equations produces a complicated solution for g(z', z) i n one dimension.  First moment estimates from a second-order moment analysis o f this model correspond relatively well to simulation results (see Fig. 2-4). The analytical solution for the first moment, which represents the mean distance from the optimum at steady state, can be approximated by the sum o f two rational polynomials:  m, ~  B(2N + BX)  BA  2N + (2N + BX)  2  2  (2.7.1)  27V(27V + 5 A )  where we have omitted the subscript from N for clarity. Equation (2.7.1) assumes that e  e  2N  »  N , and so gives a good approximation o f the exact solution when N> 10. When 3  compared to simulation results, we find that this second-order system makes an accurate prediction o f the first moment for low values o f'X(Fig. 2-5). However, for large values o f A,  35  0.05  T  0  20  40  60  80  100  120  N  Figure 2-4. - Mean fixed drift load against effective population size for n = 1. The mean fixed drift load is equal to the first moment of the distribution Peq(x) when B=\. Solid squares indicate the average of five simulation runs evaluated at A = 10, each consisting of 1000 fixed mutation steps. (Throughout, the first 50 steps were ignored to minimize the influence of starting conditions.) Crosses indicate the average of five simulation runs evaluated at A = 1, each consisting of 1000 fixed mutation steps. The solid curve represents equation (2.7.1) evaluated at A = 10, and the dotted curve represents the upper limit approximation at A = 0 (2.7.2).  36  (2.7.1) predicts that the first moment is nearly independent o f A, while simulation results show that the mean distance to the optimum decreases with A (Fig. 2-5). Consequently, (2.7.1) w i l l overestimate the mean drift load when A is large (i.e. when the average phenotypic effect o f a mutation is small).  Given that the first moment obtains a maximum value at A = 0 in simulation results, we can approximate an upper bound estimate o f the mean distance to the optimum by setting A = 0 such that (2.7.1) becomes:  (2.7.2)  B  1 + 27V  B y substituting (2.7.2) into the absolute fitness function (2.5.5), we find that the upper bound estimate o f the expected drift load from fixed mutations is:  L =— —  (2.7.3)  \ + 2N  Note that this upper bound estimate is independent o f the selection intensity, 1 / B (see E q . 2.5.6). Because (2.7.3) is a reciprocal function o f N, the mean drift load due to fixed mutations becomes very small for N> 10. For reasonably-sized populations, the probability of fixation (2.5.4) for deleterious mutations becomes very small and most such mutations are readily compensated.  37  0.012  0.01  0.008  •o cd o c co  •  • <  0.006  CD  E  0.004  •  simulation equation  0.002  10  20  30  40  50  lambda  Figure 2-5. - Mean fixed drift load against A for n = 1. The mean fixed drift load is equal to the first moment of the distribution Peq(x) when 5 = 1 . For this plot, the effective population size is set to N = 50. Points indicate the average of 10 simulation runs, each consisting of 1000 fixed mutation steps. Standard errors are too small to appear on the graph. Simulation results indicate that the first moment decreases monotonically with A. The dashed line represents the second-order approximation (2.7.1).  38  2.8 Results in the n-dimensional model The mutation probability function m(z', z) becomes very complicated when Fisher's space consists o f more than one dimension. Given a population residing at a distance z from the optimum, we must consider the cumulative probability o f every mutation vector o f some length r and orientation 6 that terminates at the surface o f a hypersphere with radius z' i n order to calculate m(z', z). It is simpler to index each mutation by the angle (f> internal to So  formed by the radii extending to A and to A':  m(z',z) = C / ; A e x p ( - A ^ z + (z'f 1  - 2z^'cos</>)sin' ,  2  040  which contains (2.6.6); C is a normalizing coefficient. A n exact solution o f this integral appears to be mathematically intractable. However, simulation results indicate that a maximum with respect to A exists for the first moment o f the distribution p (z') at A = 0. eq  Therefore, we set A = 0 to obtain an upper bound estimate o f the fixed drift load for the case of exponentially distributed mutation lengths along each axis. This assumption allows mutations in a population at a distance z from the optimum to jump to any point within the nspace with equal probability. The rc-space can then be thought o f as a "uniform probability volume" divided into a series o f concentric shells around the optimum like an onion. It follows that the probability density function o f mutations terminating on the surface o f a hypersphere centered at the optimum with radius z is:  39  (2.8.1)  where I have bounded the distribution to exclude lethal mutations (z' < B). Note that (2.8.1) is independent o f the original position z. Equation (2.8.1) is equal to the limit o f (2.6.4) as A approaches zero. To prove this, I w i l l use the following asymptotic expansion for the denominator o f (2.6.4):  r ( n ) - T(n, Afl) = e - " ( A B ) " f :  r 0  *  }  UB)  J  (Eq. 6.2.5 from Press 1992). Substituting this expansion into (2.6.4) and taking the limit, we have:  lim  £  e (XB) m  n  -Xz'  A^O  ,„  e \B) x  n  /  1  n  1  / +  W  (XBy  f\n—\  Afl  A  + ...  (n + \)n  /\n-\  n(z') B"  The moment analysis o f an ^-dimensional model using this approximation produces results which are very close to those from the simulations when A is near zero (see Figs.  40  0.35 -equation  0.3  0.2  A O  A  exponential,  O  gamma,  X=0.1 X=0.1  C  _  simulation,  X=0.1  0.25 TJ CO O  •  0.1 5 i 0.1  A O  A  0.05  A  0 20  40  60  80  1 oo  N  Figure 2-6. - Mean fixed drift load against effective population size for n = 10. The dashed curve represents the upper limit equation (2.8.3), where A = 0. Solid squares indicate simulation results with A = 0.1 using the original mutational model, which assumes exponential distributions of mutational effects along each axis that accumulate to produce a bell-shaped gamma distribution for total mutation length (2.6.4). We were unable to obtain simulation runs for values of N greater than 30, because mutations that fixed were extremely rare under these conditions. Open triangles represent the mean of five simulation runs consisting of 1000 fixed mutation steps each, with an exponential distribution (A = 0.1) describing the total mutation length. Open diamonds represent the mean of five simulation runs consisting of 1000 fixed mutation steps each, with an L-shaped gamma distribution (A = 0.1, a = 0.5) describing the total mutation length.  41  0.6 0.5 -I 0.4 x  T3  ro o 0.3 c  £  s  0.2 \  /  A  Jl-  A  -equation  A  0.1  •  simulation,  A  simulation,  X  approximation  ^=1  X=20  0 0  10  15  20  25  # dimensions (n)  Figure 2-7. - Mean fixed drift load against n. For this plot, the effective population size is set to N= 10. The dashed curve represents the upper limit approximation (2.8.3), where A = 0. Solid squares indicate simulation results for A = 1. We were unable to obtain results for n > 10, because fixation events became extremely rare under these conditions. Crosses indicate the means of five simulation runs consisting of 500 mutation steps for a truncated distribution at A = 0 that generated only values of r that were likely to fix. Truncating the distribution decreased computation time but in exchange slightly overestimated the load for larger values of n. Closed triangles indicate the average of five simulation runs for A = 20, each consisting of 1000 fixed mutation steps.  42  2-6, 2-7). The upper limit solution for the first moment can be approximated by the sum o f a simple rational polynomial and a vanishing term:  nB  n + 2N  J O  (2N) en  2N  ^  (2.8.2)  r(n)-T(n,2N)  The neglected terms are very nearly zero i n the domain N> n. Even for N<n,  simulations  indicate that nB I (n + 2N) is a good approximation for the mean distance to the optimum for low A and provides an upper bound for m\ as A increases. Note that the first term reduces appropriately to (2.7.2) i n one dimension (n = 1). Consequently, the mean equilibrium drift load caused by the fixation o f new mutations is:  n + 2N  when A is small, which overestimates the load when small mutations are common (higher A). Equation (2.8.3) is again a reciprocal function of N, such that the mean drift load due to fixed mutations rapidly decreases with effective population size. A l s o , (2.8.3) is a hyperbolic function o f the number o f dimensions n.  2.9 Sensitivity of the model to its assumptions The features o f Fisher's geometric model o f adaptation that make it mathematically useful also make it biologically oversimplified. One frequent criticism o f Fisher's model is that its spherical symmetry is too idealized to apply to real organisms. A l l o f the orthogonal  43  axes are standardized to fit under the same fitness function and mutational distribution. With respect to fitness, we can change the scale for each axis such that the same displacement experiences the same selection intensity. However, this would alter the mutation probability associated with a displacement along each axis. Consequently, the spherical geometry o f Fisher's model cannot perfectly capture both the fitness effects o f mutations and their frequency distribution. The model also assumes that the mutational distribution is symmetrical along each axis, with an equal probability o f going towards and away from the optimum. If the number o f potentially compensating loci varies among traits, for example, then this assumption would not hold. Hence Fisher's model is almost surely an imperfect representation o f a biological system. O n the other hand, models that ignore compensatory mutations are also imperfect representations. Given that neither Fisher's model nor a model lacking compensatory mutations perfectly captures reality, modeling both is still instructive because it allows us to determine the sensitivity o f evolutionary predictions to two extreme alternatives.  There is good reason to believe that the results obtained from this formulation are robust enough to changes i n the assumptions to be useful for understanding the fixed drift load in small populations. There is a handful o f theoretical studies that have either examined a related problem in a roughly similar fashion (Lynch and Gabriel 1990; Haiti and Taubes 1998) or addressed a different evolutionary problem i n a model analogous to Fisher's geometric model (Robertson 1970). It is interesting that the results o f all these studies bear a non-trivial resemblance to those presented here. Robertson (1970) studied the load caused by the drift o f two alleles maintained near a polymorphic equilibrium q by heterozygote  44  superiority. The departure o f allele frequencies from the equilibrium could be represented by the variance o f Wright's (1937) distribution, giving an expected load of:  (2.9.1)  l + 4Nq(l-q)  (Robertson, 1970). This model is superficially analogous to Fisher's model for n = 1, in that there is an optimum located on a real line, over which selection and drift act antagonistically to determine a steady-state distribution. Note that equation (2.9.1) and the result obtained by the one-dimensional model, (2.7.3), are similar in their inverse dependence on population size.  A simulation study o f the mutational meltdown o f small populations by L y n c h and Gabriel (1990), like Fisher's model, allows for compensatory mutations. In this model, Af individuals experience a Poisson distribution (Feller 1951) o f mutations whose effects on fitness have a constant mean and variance. Although mutations are on average deleterious, there is a constant probability that a beneficial mutation w i l l appear i n this model, with the restriction that individuals with a fitness greater than one are not allowed. Consequently, the model is similar to ones investigating the incorporation o f unconditionally beneficial mutations (Lande 1994; Schultz and L y n c h 1997), except that beneficial mutations only appear once fitness has decayed to some extent. This differs from Fisher's model, where the probability o f a beneficial mutation continues to increase as the population moves further from the optimum and approaches 1 / 2 . L y n c h and Gabriel (1990) found that allowing for beneficial mutations caused the mean extinction time o f the population to increase by orders  45  of magnitude. The mutational meltdown is not frozen, however, unless the fraction o f beneficial mutations is sufficiently high (Schultz and L y n c h 1997).  More recently, Hartl and Taubes (1998) estimated the fixed drift load in their study o f adaptation in Fisher's model as:  (2.9.2)  Although we have studied the same model here, different assumptions have been made about the nature o f mutations. We have evaluated the effect o f each difference by running simulations under various alternative assumptions. First, Hartl and Taubes (1998) used an absolute selection coefficient s = w(z') — w(z) whereas the probability o f fixation for a new mutation is really a function o f the relative selection coefficient s = (w(z') - w(z)) I w(z) (Eq. 2.5.3; i.e. selective advantage, Fisher 1930). Effectively, their assumption was equivalent to relaxing selection away from the optimum. According to our simulations, changing s from an absolute to a relative measure in Hartl and Taubes' (1998) study produces results that can be approximated by:  L=  n n + SN  (2.9.3)  Note that whereas the absolute selection coefficient causes the load estimate to have a linear dependency on the number o f dimensions, n, (2.8.3) and (2.9.3) are nonlinear functions  46  which asymptote at 1 as n gets large. The remaining difference between (2.8.3) and (2.9.3) results from what is assumed about selection and mutation. Hartl and Taubes (1998) used a quadratic fitness curve instead o f a linear one (Eq. 2.5.5). Simulations indicate that the term 8N in (2.9.3) becomes roughly AN with a linear fitness function, suggesting that the fixed drift load is lowered when the fitness function curves downward away from the optimum (i.e. with synergistic epistasis; Schultz and L y n c h 1997). Furthermore, I had assembled a probability distribution for mutation length, r, from component distributions along each axis in a bottom-up fashion, whereas Hartl and Taubes (1998) assumed that the total mutation length in n dimensions was uniformly distributed (see below). This difference explains the remaining discrepancy between (2.9.3) and (2.8.3). That both estimates o f the fixed drift load are inversely related to the effective population size, however, implies that this result from Fisher's model is fairly robust.  Comparing the results from Hartl and Taubes (1998) and this study highlights our ignorance about how mutations affect phenotype. Recall that I have assumed that the distribution for mutation length along each axis is exponential and that the total mutation length is derived from these component effects (a bottom-up derivation). However, the conventional procedure with respect to Fisher's model is to select a specific shape to the total distribution o f mutational lengths, leaving the component distributions along each axis unspecified (a top-down derivation) (Kimura 1983; Orr 1998; Hartl and Taubes 1998). Little data exist concerning the shape of the distribution o f mutational effects on phenotype. There is a general consensus, however, that the distribution o f deleterious mutant effects on fitness is L-shaped (Mackay et al. 1992; L y m a n et al. 1992; Keightley 1994; but see Garcia-  47  Dorado 1997). For mathematical reasons, it is difficult to derive a probability density function o f mutant effects on phenotype that is L-shaped using the bottom-up approach (Orr 2000). Nevertheless, the mutant effects on fitness may still have an L-shaped distribution because the phenotypic changes must first be mapped onto some fitness function. This mapping depends, however, on the current state o f the population, as well as on the number of dimensions o f the system in Fisher's model. Recall that mutations become less likely to be oriented directly towards or away from the optimum as the population moves further away from the optimum and as n increases (Figs. 2-1, 2-3; E q . 2.5.3; Hartl and Taubes 1998). In other words, mutation vectors become increasingly likely to be tangential to the surface o f the hypersphere 5*0. Consequently, mutations with vanishingly small effects on fitness may increase in frequency. In this case, a roughly L-shaped distribution o f fitness effects may be seen even when the distribution o f phenotypic effects is bell-shaped (as in the bottom-up derivation).  O n the other hand, very large populations (such as the Drosophila populations used to generate the expectation o f an L-shaped distribution o f mutational effects on fitness) are expected to reside at or close to an optimum. If this is true, then the observed distribution o f fitness effects should be similar i n shape to the distribution o f total mutation length for populations near an optimum. Consequently, it is worth investigating mutations whose total length follows an L-shaped distribution. We have run simulations using an exponential distribution and an L-shaped gamma distribution (that can be obtained from (2.6.3) by setting n = 1 and 0.5, respectively, with A = 0.1 in both cases) to describe the distribution o f total mutation lengths i n n dimensions (Fig. 2-6). In both cases, the dependency o f the fixed drift  48  load on the effective population size, N, is qualitatively similar to our original result (Fig. 2-6). Therefore, results from Fisher's model remain fairly robust under varying mutational models.  2.10 Implications for the extinction risk of small populations The most important result o f this model is that changing the way mutations are modeled results i n qualitatively different predictions regarding the extinction risk o f small populations. The inevitable extinction o f small populations is a common feature o f genetic models developed within the classical framework o f mutations whose effects on fitness are unconditionally deleterious and cannot be compensated (Lynch et al. 1993; Lande 1994). A l l o w i n g for compensatory mutations, however, we find that mean fitness does not suffer an inexorable decline but rather reaches a steady-state level. This "freezing" o f the mutational meltdown occurs because, as the population moves further from the optimum, a higher fraction o f mutations increase fitness. Although it seems unlikely that all mutations can be compensated (as in Fisher's model), it is equally unlikely that compensatory mutations can be safely ignored (as i n the classic framework). Whether reality lies nearer the level o f compensation assumed in Fisher's model or in the classic framework is an empirical matter that has yet to be addressed. Such data are needed to resolve whether small populations are at risk o f genetic meltdown and to determine the rate o f this meltdown.  Some previous models o f the accumulation o f mutations i n small populations focus on unconditionally deleterious mutations with a constant fitness effect s occurring at a constant rate /I (Lynch et al. 1993; L y n c h and Gabriel 1990, but see discussion above).  49  Similar models draw unconditionally deleterious fitness effects from a continuous distribution (Lande 1994) or both deleterious and advantageous mutations at different rates (Schultz and Lande 1997; Lande 1998). In these examples, a population is unable to halt or even slow its fitness decline once it has begun, because then the influx o f deleterious mutations overwhelms that o f beneficial mutations. I f deleterious mutations are progressively less likely to accumulate as mean fitness decreases — as occurs when there is synergistic epistasis among deleterious mutations (Kondrashov 1994; Schultz and L y n c h 1997) — the population reaches a point beyond which its mean fitness no longer decreases. Schultz and L y n c h (1997) argue, however, that an unrealistically high level o f synergistic epistasis is required. Interestingly, the meltdown is similarly frozen i n Fisher's model, in which a population at the optimum has no beneficial mutations available to it and an increasing number as it declines in fitness. B y evaluating the problem i n terms o f Fisher's model, we are explicitly considering the context-dependence o f a mutation's effect on fitness.  If there is a substantial class o f mutations that are readily compensated, we expect . that the drift load caused by the fixation o f new mutations i n small populations is unlikely to be sufficient to cause the extinction o f all but the smallest populations. I have shown that the load caused by the fixation o f deleterious mutations i n the range N > n is too small to affect a population with a modest reproductive excess. There may exist, however, a sub-class o f mutations that cannot be compensated. Such mutations cannot be modeled using Fisher's model and w i l l contribute to a mutational meltdown (Lynch and Gabriel 1990), albeit at a slower rate. O f course, other risks o f extinction including demographic stochasticity and ecological degradation are more than sufficient to place many small populations in peril  50  (Lande 1988, 1993). Furthermore, small populations may be unable to adapt quickly enough to a changing world, which would also place them at risk (Burger and L y n c h 1995). The pervasiveness o f anthropogenic disturbance o f habitats may reduce the population size o f species to a point where all these forces can act i n concert to cause extinction. Therefore, although our analysis lends insight into the effects o f mutation on small populations when compensatory mutations are allowed, it would be unwise to base decisions regarding the conservation o f endangered populations solely on these results.  2.11 Implications for pleiotropy Hartl and Taubes (1998) have discussed the relationship between the number o f dimensions, n, i n Fisher's adaptive space and concrete biological properties, such as the number o f loci underlying a quantitative trait. It is also possible that n reflects the extent o f pleiotropy in the average mutation (Orr 1998). Because this analysis substantially alters predictions o f how n affects the fixed drift load, it is worth discussing the biological significance o f n. To avoid confusion, I w i l l denote the number o f axes i n Fisher's model as « F and use n to denote the hypothetical number o f components that make up an organism's phenotype. Although there is an indefinitely large number o f measurable characters in any organism, this does not necessarily imply that n is also large. This is analogous to a principal components analysis o f multiple morphological measurements, where most o f the variation present can be accounted for by a much smaller number o f axes. W e can partition n into discrete sets and apply Fisher's model to each set, effectively making different assumptions about how mutation and selection act on these different trait combinations. A t one extreme, letting «F = n implies that the average mutation w i l l nearly always have a non-zero effect on  51  every phenotypic component, which is biologically equivalent to being completely pleiotropic. A t the other extreme, n can be completely partitioned into n individual sets, each represented by a one-dimensional version o f Fisher's model (i.e. n?= 1). This assumes that every mutation w i l l affect only one of the n phenotypic components at a time; there are no pleiotropic effects on any other components. The subsequent fixed drift load is predicted by the product o f n and the one-dimensional result, 1 / (1 + 2N) i f loads combine additively across sets. This load without pleiotropy, n I (1 + 2N), is always greater than the load with complete pleiotropy, nl(n  + 2N), for n> 1. This implies that the pleiotropy inherent i n  Fisher's model substantially reduces the predicted drift load caused by new mutations, because every mutation can compensate several components at once.  Alternatively, let n represent how an organism's phenotype w i l l respond to mutation and selection, rather than being a description o f the phenotype itself (Orr 2000). Increasing n within Fisher's model (with «  = F  n) w i l l decrease the chance that a new mutation bring a  population closer to the optimum (Fisher 1930). If the organism's phenotype evolves as i f it were a single character (n = 1), then half o f the mutations that occur point in an opposite direction to the remaining half. When n is greater than one, then mutations do not affect all characters equally. Because the chance that maladapted characters w i l l be little affected by a new mutation increases with n (see Fig. 2-3), the delay between deleterious and compensatory mutation events also increases. This explains why the fixed drift load is an increasing function o f n within Fisher's model.  52  I have reviewed the experimental literature that provides good evidence that compensatory mutations occur (Section 2.3), but there is a lack o f empirical measures for what the overall extent o f compensation might be. From Fisher's model, we can predict how n is related to measurable quantities, such as the fixed drift load and the average degree o f compensation among mutations, that can be used to extrapolate estimates o f n from empirical data. Such interpretations are only as relevant, however, as Fisher's model is an accurate representation o f biology. That predictions concerning the extinction risk o f small populations depend qualitatively on this issue emphasizes the importance o f further research on compensatory mutation.  53  CHAPTER  III:  T H E EFFECTS OF DEMOGRAPHIC HISTORY  3.1 Overview Most o f the theoretical studies that I have discussed so far have been concerned with the evolutionary consequences o f a constant small population size in relation to the genetic load (e.g. L y n c h and Gabriel 1990; Lande 1994). This "small-population paradigm" (sensu Caughley 1994) is common because it is mathematically simpler to model allele frequency evolution in a population o f constant size. The evolutionary consequences o f population decline, however, do not readily accommodate theoretical study (Caughley 1994). It remains to be understood how changes i n the size o f a population (i.e. its demographic history) affect allele frequency evolution and the corresponding genetic load. Because fluctuations in the allele frequency accumulate over time, alleles in a finite population w i l l rarely occur at the equilibrium frequency determined by mutation-selection balance (see Section 1.1). The corresponding change i n the genetic load is one consequence o f the demographic history. In this chapter, I w i l l investigate how the demographic history o f a population that has recently declined from a very large size affects the transient genetic load. In other words, consider a population that is currently at a small size, but was at a very large size some known number of generations ago. With all else being equal, how does the size o f the population in the intervening generations affect the current genetic load?  3.2 Steady-state expectations of genetic load  54  To understand the effect o f demographic history on the genetic load, we must place results o f the model in the context o f the steady-state load expected i n populations o f constant infinite or finite size. The genetic load caused by deleterious alleles at equilibrium in an infinite population for an arbitrary degree o f dominance, h, is given by Crow (1970) as being approximately:  L = n + qhs(\- ju)  (3.2.1)  where s is the selection coefficient o f the deleterious allele, ju is the forward mutation rate to the allele, and q denotes the equilibrium allele frequency under mutation and selection (Haldane 1937). When the deleterious allele is completely recessive such that h = 0, (3.2.1) simplifies to fi. If h is appreciably greater than 0, then q can be approximated by fi I (hs) and (3.2.1) can be approximated by 2^u; on the other hand, i f h is very nearly 0, then (3.2.1)  lies between jl and 2^u and the exact equilibrium solution for q is required (Crow 1970). Note that the load caused by a deleterious allele in an infinite population is independent o f the strength o f selection against the homozygote, s. This is equivalent to stating that every mutation with some heterozygous expression, whether lethal or nearly neutral, contributes equally to the genetic load i n an infinitely large population.  Numerical results obtained by K i m u r a et al. (1963) indicated that the steady-state genetic load i n a finite population caused by a moderately deleterious allele (s < 0.1) increases above the equilibrium (3.2.1) with smaller constant population size, regardless o f  55  whether the allele was additive or recessive (h = 0.1). Moreover, the steady-state load i n either case was never less than the value attained in an infinite population (Kimura et al. 1963).  Their most interesting result, however, was that the steady-state load caused by a  slightly deleterious mutation was greater than the load caused by a comparatively severe mutation until the population size became small enough that drift overcame selection against the severe mutation. For s > 0.5, however, the required population size was very small (TV < 10). With the exception o f some special cases, the steady-state distribution o f allele frequency obtained by Wright (1945) and K i m u r a  (1955)  derive the moments o f the genetic load (Kimura et al.  cannot be explicitly integrated to  1963).  In the case o f lethal alleles (e.g.  h = 0,s= 1), however, the steady-state allele frequency distribution can be approximated by:  - ^\xp(-2Nx )x ^ 2  2  (Crow and K i m u r a  1970).  4N  (3.2.2)  Integrating the product of the load caused by a completely-  recessive lethal allele, L = x , over the distribution 2  -  (3.2.2)  gives the mean load:  r q + 2 ^ , 2 ^  *\  ni+2Nfi)  (3.2.3)  j  where T(-) denotes the gamma function. The second term is vanishingly small for all u. and JV> 0,  such that  (3.2.3)  is well-approximated by u,. Hence, the mean steady-state load caused  by a completely recessive lethal allele is essentially independent o f population size. Together, these results imply that mildly deleterious alleles (s < 0.1) are more important in  56  the long-term effect o f a reduction in size on the genetic load (Kimura et al. 1963). Results from the steady-state distribution, however, conceal the transient evolution o f the genetic load that can be revealed by iterative calculations o f the distribution o f genetic load. Before I discuss the basis o f such calculations, however, I w i l l discuss a process that has an important effect on the transient load.  3.3 M a r g i n a l fitness and purging  A change i n allele frequency does not only alter the genetic load, but it can also alter the marginal fitness o f a deleterious allele. The marginal fitness o f an allele is the fitness averaged over all possible zygotic combinations between the allele and other alleles in the population. W i t h inbreeding, the marginal fitness o f an allele is more influenced by its effect in homozygotes. During a transition from panmixis to an inbred mating system, the onset o f inbreeding may therefore be accompanied by a change in the marginal fitness o f deleterious alleles. The effect that a change o f an allele's marginal fitness has on its frequency in subsequent generations depends on the dominance relationship between alleles i n the population. Because homozygotic combinations become more common with inbreeding, the marginal fitness o f a recessive deleterious allele declines with inbreeding. The reduction in the mean fitness o f a population that accompanies this decline is called inbreeding depression (Charlesworth and Charlesworth 1987). In the generations following the onset o f inbreeding, the overabundance o f homozygotes o f the recessive deleterious allele may cause more copies of the allele to be eliminated by selection from the population than expected under random mating (i.e. before inbreeding). This phenomenon is called purging (Templeton and Read 1984). One motivation for the large amount o f literature (Lacy and Ballou 1998; F u 1999;  57  F u et al. 1998; Hedrick 1994; Barrett and Charlesworth 1991) devoted to studying purging is the possibility that prolonged inbreeding may be a useful technique for improving the "genetic quality" o f small populations (Templeton and Read 1984; Hedrick 1994; Hedrick and M i l l e r 1992). However, because the elimination o f alleles causes a temporary reduction in the mean fitness o f the population, the inbreeding process may itself cause extinction (Hedrick 1994; F u etal. 1997).  Because the marginal fitness o f a deleterious allele also changes with the frequency o f the allele in the population, purging is not unique to changes i n the mating system but can also be produced by genetic drift (Wang et al. 1999; Kirkpatrick and Jarne 2000). A recessive deleterious allele that becomes over-represented i n the population after genetic drift forms more homozygotes, thus exposing the allele to increased selection (Lande and Schemske 1985; Wang et al. 1999; Kirkpatrick and Jarne 2000). This phenomenon has recently been studied in theoretical models o f finite populations initially sampled from an infinitely large ancestral population, and subsequently held at a constant size (Wang et al. 1999) or allowed to return to the original population size (Kirkpatrick and Jarne 2000). Although such bottleneck models can be associated with a founding event (i.e. the colonization o f an unoccupied habitat) or a sudden and very severe truncation in population size, it is less useful for studying the dynamics o f load in gradually declining populations. For mathematical reasons, the analytical methods available for the study o f the evolution o f genetic loads due to recurrent mutation are not easily extended to declining populations. Therefore, I have chosen to conduct a numerical investigation o f the transient genetic load in a declining population.  58  3.4 Numerical methods I have carried out exact calculations i n a Wright-Fisher model o f allele frequency evolution i n a randomly-mating diploid population as described by Ewens (1979), in which the dot product o f the probability distribution o f allele frequencies at the t-th generation, P , t  and the matrix o f transition probabilities, M , gives the distribution o f allele frequencies in the next generation. Each distribution P is a vector o f 2N + 1 elements, and the matrix M i s a t  t  rectangular (2A +i+l)x(2A +l) matrix. Modifying the size o f the matrix and the associated /  /  t  t  vectors permits changes in population size (Tyvand 1993). The /, k-th entry o f M i s equal to the transition probability from k to i alleles in the population according to the binomial formula:  . 2AT, -i +1  1-  P{i\k) =  (3.4.1)  2N.tJ  (Feller 1950) where the allele frequency k/(2N) may be modified by selection and mutation. I assumed equal forward and backward rates o f mutation at a typical value o f 1 0 . Thus, -6  successive generations were calculated with evolutionary processes occurring i n the following order: selection, mutation, genetic drift. A t the end o f each cycle, the nth moment of the probability distribution o f genetic load at that generation was calculated from P  x  according to the formula:  (L') = _  k ,2 - - k k V p(k,t). :(1-: 2N(t) 2N(t) 2N(t)  59  (3-4.2)  where p(k, t) is the k-ih. element o f P . The mean and variance o f load are easily obtained t  from the first and second moments obtained from (3.5.1) (Kendall and Stuart 1963).  I have modeled a step-wise decline in population size. The population size is assumed to be initially large, declines to a size o f N* for T generations, and then declines again to its final size of N'. I explored two different classes o f trajectories (Figure 3-1). The first ("glide" trajectories) set % = 10, N' = 10, and varied the intermediate iV* from 10 to °°. The second ("shift" trajectories) set N* = 100, N' = 10, and varied x from 1 to 10. Assuming a very large initial population, the allele frequency distribution in generation 0 was the equilibrium allele frequency determined by mutation and selection. The exact numerical solution o f the equilibrium frequency was used instead o f Haldane's (1937) approximations because the degree o f inaccuracy in the approximate solution varies with the parameters s and h and can obscure the results.  I used computer algorithms written in the C programming language to perform exact calculations o f the Wright-Fisher model. The number o f calculations required per generation was on the order o f N since each allele frequency class k I (2N ), k= {0, 1,2, ..., 2iVt}in the 2  t  vector P o f length 2JV +l could make some contribution to every class comprising the t  t  distribution vector in the next generation,  P +it  Considering the number o f calculations  required to take mutation, selection, and drift into account, the total number o f calculations becomes prohibitively large for large populations. Approximate methods that lower the number o f calculations required for large population dynamics, such as combining adjacent  60  F i g u r e 3-1. - The glide and shift trajectories. In each case, population size is plotted against elapsed time in generations. It is assumed that the first generation is sampled from an infinitely large population. The glide trajectories, represented by the left plot, are step-wise trajectories that summarize the effect of historical population size preceding the final generation to a single value, N*. Upon the final generation (T = 10), the population makes a second decline in numbers to If = 10. The shift trajectories, represented by the right plot, are step-wise trajectories that summarize the effect o f the rate o f decline to a single variable, x, that corresponds to the generation upon which the second decline from N* to i V occurs.  shift  glide x=10  N  i  i  +  +  •  i  time  time  61  i  +  i  »  allele frequency classes into "bins" whose behavior is determined by the mean allele frequency o f its member classes, failed to achieve a sufficient degree o f accuracy to be useful. Consequently, the dynamics o f the allele frequency distribution were explored solely for small populations (N < 100). Recall that the effect o f drift on the allele frequency distribution is inversely proportional to population size. Because the distribution changes slowly while the population is still large, a disproportionately greater part o f allele frequency evolution occurs while the population is o f small size. Given that the decline is recent, it is a reasonable approximation to omit the generations when the population is large so that the population immediately declines from an infinite to small size (N< 100).  The calculations were not quite exact because o f limits on the accuracy o f computer operations with 16-bit processing. The binomial distribution (3.4.1) is prohibitively difficult to implement for large N(N>  50) because several factors become vanishingly small (e.g.  high powers o f x) or excessively large (i.e. the factorial AH). One method for circumventing this problem is by decomposing a function into parts so that extremely large or small numbers do not appear. For example, the binomial coefficients in (3.4.1) can each be  a\  \  a  (a-b)\b\  b  r  expressed as an iterative product (e.g.  a  •x  a-I b-1  x...x  a-b + 1. 1  ). Such  iterative calculations, however, can be time-consuming. Alternatively, performing the calculations on a logarithmic scale and subsequently back-transforming the results can also prevent the use o f extremely large or small numbers. I used this second method for calculating the binomial coefficients, taking advantage o f the fact that the factorial o f an integer n can be written as a gamma function:  62  w! = r(w + l) = J z V ' d z 0  (Abramowitz and Stegun 1972). B y applying an approximation obtained from Press (1992, p. 213), I obtained an approximate formula for the binomial distribution that I used for all population sizes:  x' (1 - x) ~' = exp{ln T(N +1) - In F(N - i +1) - In r(i +1) + i In x + (N - i) ln(l - x)} N  V 'J  where the notation in the binomial formula (3.4.1) has been changed to avoid cluttering the right-hand side. This formula allowed fast calculations o f the binomial distribution that were sufficiently accurate in comparison to exact calculations for small population size.  3.5 N u m e r i c a l results The numerical results o f the calculations from the glide trajectories are presented i n Figures 3-2 and 3-3. The general qualitative result o f these figures is that the mean genetic load caused by a mutant allele depends not only on the dominance and strength o f selection but also on the demographic history o f the declining population, represented by N*. Only the results for two members o f the glide trajectories are shown for clarity: fast corresponds to the population with the smallest size before x = 10 (TV* = 10), and slow, corresponds to the population with the largest size before x = 10 (N* = 100) studied. In the additive model (Fig. 3 -2a), the load caused by a lethal ( 5 = 1 ) allele remains at the equilibrium value, 2jx (Eq. 3.2.1), regardless o f N*. This is supported by the fact that the initial sample from the infinite  63  ancestral population (i.e. generation 1) has the same mean load irrespective o f N* and s and that the mean additive load is a linear function o f allele frequency when mutant homozygotes are absent (Eq. 3.4.2). In contrast, the mean load o f non-lethal additive alleles increases above 2\i at a rate that is greater for smaller N* (Fig. 3-2o). Note that the load for s = 0.5 increases more rapidly but is eventually overtaken by the load for 5 = 0.1 regardless o f N*. Relative to the load o f recessive alleles (Figs. 3-2b, c), however, the transient mean load o f additive alleles is very slow to increase from 2/1.  Drift assumes a more complicated role in the transient evolution o f partially recessive alleles in declining populations. This is illustrated by Figure 3-2Z>. Unlike the additive case, the initial sample o f size iV* from the infinite ancestral population results i n a larger mean genetic load with decreasing N*. The mean load increases at a higher rate for weaklyselected alleles (s = 0.1) than lethal alleles when partially recessive, regardless of N*. In the fast population, the purging o f lethal alleles eventually overcomes this increase. The slow population also experiences a slight reduction i n mean load after increase. When the slow population declines to N' - 10, drift becomes stronger and the mean load jumps upwards as a result. This is similar to the initial sample o f  = 10 individuals from the ancestral  population (i.e. fast) except that the distribution o f load has been changed slightly over nine generations o f being at N* = 100. Because the decline to N' = 10 is recent i n the slow population, its mean load due to partially recessive lethals is greater than that o f the fast population, in which selection has had the opportunity to purge lethal alleles.  64  h=0.1  h=0.5  — ° —  s=0.1,  fast  —o—s=1.0, fast —o—  -s=0.1,  slow  — -o— -s=1.0, slow  1  2  3  4  5 6 7 8 generations  9 10  2  3  4  5 6 7 8 generations  9  10  h=0.0 120.0  100.0  •s=1.0, fast •s=0.1, fast o - - •s=1.0, slow —o - - •s=0.1, slow  & - h=0.1 * - h=0.5  4 5 6 7 generations  9  10  F i g u r e 3-2. - Mean genetic load (xlfT ) of glide trajectories plotted against time in generations. 6  Individual plots correspond to (a) additive, (b) partially recessive, and (c) completely recessive models of selection, respectively. Note that the scale of genetic load is much larger in (c) than for (a) and (b). The equilibrium load in an infinite population is marked by an arrow and the value corresponding to the degree of dominance. The fast trajectories represent populations that immediately declines from being infinitely large to N* = 10. The slow trajectories represent populations that decline from being infinitely large to N* = 100 and remain at that size until the final generation, whereupon they decline to Af = 10. In the additive lethal case, the fast and slow trajectories overlap. Note that in each case at generation 10 the population size was the same (AT = 10) across trajectories when load was calculated, but that all selection (and purging) occurred in previous generations (when N = A *). 7  65  F i g u r e 3-3. - Mean genetic load (xlO ) of glide trajectories at generation 10, plotted against s. -6  Individual plots correspond to (a) additive, (b) partially recessive, and (c) completely recessive models of selection, respectively. For clarity, only the results for N* = {10, 20, 100} are plotted. Note that the scale of the mean load varies among plots. The effect of demographic history is greatest for completely recessive alleles, with the greatest divergence caused by s = 0.1 at the 10 generation. The genetic load caused by completely recessive lethal alleles in the N* = 20 trajectory exceeds both the smallest (N* = 10) and largest (N* = 100) trajectories studied. The additive plot at s = 0.0001 goes below 2p because the corresponding equilibrium genetic load is approximately equal to 1.96. th  66  Figure 3-2c shows numerical results for completely recessive alleles. The overall transient genetic load is roughly an order o f magnitude greater than the load in the previous two cases o f additive and partial dominance (Fig. 3-2c). The initial sample o f N* individuals from the infinite ancestral population causes a larger jump in the mean load for N* = 10 (fast) than N* = \00 (slow). This jump is very large overall in comparison to the partially recessive case because the equilibrium frequency o f lethal alleles (~\i/(hs); Haldane 1937) is greater in the ancestral population when completely recessive, causing the initial spread o f the allele frequency distribution by drift to extend to higher values. Again, the purging o f lethal alleles in the fast population eventually overcomes the increase o f load due to drift so that the mean load is greater i n the slow population (Fig. 3-2c). In the short-term, the mean load caused by lethal alleles is greater than that caused by non-lethal (s = 0.1) alleles. Eventually, the mean load o f a completely recessive lethal allele converges in the long-term on the steady-state value jx that is less than the steady-state load caused by non-lethal alleles (Section 3.2). This example illustrates that the transient dynamics can be substantially different than steady state expectations.  Figure 3-3 again presents results from the glide trajectories except that now the abscissa is a log-transformed axis o f the selection coefficient, s, and the ordinate is the mean genetic load at the last generation (when A ' = 10). In the additive model (Fig. 3-3a), the 7  mean genetic load o f both lethal and weakly-selected (s < 10 ) alleles is equivalent to the corresponding equilibrium load in an infinite population, irrespective o f the historical population size, N*. Moderately-selected alleles ( 1 0  -3  < s < 0.5) cause an increase in the  mean load for lower N*. Again, the mean load at generation 10 is greater overall for partially  67  recessive alleles (Fig. 3-3/3), being greatest for intermediate selection coefficients for all N*. For partially or completely recessive alleles, purging in populations with smaller N* causes the mean load plots to cross over at approximately s = 0.55, indicating that selection must be severe for the population with a historically smaller size to have the lesser mean load by the final generation. For completely recessive alleles (Fig. 3-3c), the selection coefficient causing the greatest mean load increases with N*. The plots no longer cross over at a single point as is the case for partially recessive alleles.  Figures 3-4 and 3-5 illustrate the results from the shift trajectories, i n which populations share the same historical population size (N* = 100) and final size JV = 10 but differ in the time at which the transition from N* to N' occurs, x. Figure 3-4 is analogous to Figure 3-2c, as only results for lethal alleles are shown. Because the mean load remains at 2p in the additive lethal case for all N*, there are no differences among shift trajectories, and I have elected to omit this plot from Figure 3-4. For partially or completely recessive alleles (Fig. 3-4a, b), the evolutionary trend o f mean load i n a population that declines in later generations roughly mimics the trend experienced by a population that declines early (e.g. x = 1). The mean load immediately increases upon the decline in numbers at generation x and eventually decreases due to purging. Although the mean load changes while the population is at N* = 100 (see x = 10 trajectory), this is not reflected i n the maximum mean attained by different trajectories. For example, in the completely recessive case (Fig. 3-4b), the mean load while the population is at N* increases over time but populations having spent more time at N* reach lower maximum values o f mean load.  68  1.8 -  1.6  -i  1  1  1  2  1  3  1  4  1  1  5  6  1  7  1  8  1  9  —  10  generations  h=0.0  b. 140.0  0  1  J  ,  ,  1  2  , 3  , 4  ,  ,  5  6  , 7  , 8  , 9  10  generations  F i g u r e 3-4. - Mean genetic load (xlO ) of shift trajectories against time in generations. -6  Individual plots correspond to (a) partially recessive, and (b) completely recessive models of selection, respectively. Only the case for lethal alleles (s=l .0) is shown. The mean genetic load does not change from 2 p for additive lethal alleles. Each shift trajectory represents a population that immediately declines from infinite size and remains at N* = 100 until the x-th generation, upon which it declines to its final population size, N* = 10. The mean load jumps upwards upon the second decline to JV by a similar amount at x = 10 as when x = 1 despite having spent a number of generations at N* = 100. Note, however, that the maximum mean load attained declines gradually with increasing x.  69  1.92 -  1  1  0.0001  1  0.001  — i  r-i  0.01  0.1  selection coefficient  r  1  F i g u r e 3-5 - Mean genetic load (xlO ) of shift trajectories at generation 10, plotted against s. -6  Individual plots correspond to (a) additive, (b) partially recessive, and (c) completely recessive models of selection, respectively. The shift trajectories for x = 1 and x = 10 are equivalent to the N* = 10 and N* = 100 plots of Fig. 3-3. Trajectories for x = 9 are similar to x = 10 for all degrees of dominance. For x = 5, however, there are substantial departures of the mean load for severely deleterious and recessive alleles. This results from a delay of the increase in the load following the second decline of the population to TV = 10, as illustrated in Fig. 3-4.  70  Figure 3-5 is analogous to Figure 3-3. When the mutant allele is additive, the mean genetic load o f a population that historically declined from N* to A ' is a decreasing function 7  o f x (Fig. 3-5a). Thus given a population's current size, the load is least when its historic size was large and the decline in size was rapid and recent. This relationship between the load and x remains true for partially or completely recessive alleles that are mildly deleterious (s > 0.5), but changes into a non-linear relationship for severely deleterious alleles (s < 0.5) (Fig. 3-5b, c). In the case o f lethal alleles, the load caused by shift trajectories with intermediate values o f x can exceed the load caused by upper (x = 10) and lower (x = 1) histories, which reflects the delayed increase o f the mean load revealed by Figure 3-4.  3.6 C o n c l u d i n g remarks Genetic drift causes the mean genetic load o f a population to increase either by causing deleterious alleles to fix or by increasing the variation o f the unfixed allele frequency distribution. The long-term increase in the mean load due to the fixation o f alleles by drift is a decreasing function o f s because the fixation probability o f deleterious alleles decreases exponentially with increasing s. Thus, slightly deleterious alleles are more important in the long term with respect to the transient load o f a small population. Because lethal alleles cannot be fixed in the population, changes in the mean genetic load caused by lethals must be brought about by the second process. Increasing variation in allele frequency due to drift can be converted to an increase in the mean genetic load i f the load is a non-linear function o f allele frequency (i.e. i f the allele is recessive; E q . 3.4.2). Hence, there is no change in the mean load caused by lethal alleles in the additive case (e.g. F i g 3.2a) because the alleles  71  cannot become fixed and because increasing variation does not change the mean load. Variation in the frequency o f completely recessive lethal alleles, on the other hand, causes a greater increase i n mean load than slightly deleterious alleles i n the short term (Fig 3-2c). The amount o f variation in allele frequency caused by genetic drift between generations is greater for smaller population size. If a population remains at a small size for several generations, then the cumulative allele frequency variation eventually causes the purging o f severely deleterious alleles to overcome the increase o f mean load due to drift. Thus, purging acts in opposition to drift. This rise and fall in the mean load due to recessive lethal alleles is a transient phenomenon that occurs at small population size (Fig. 3-4). Hence, recessive lethal alleles are more important than slightly deleterious alleles i n determining the transient genetic load o f a small population in the short term. This confirms numerical results obtained by previous investigators (Wang et al. 1999; Kirkpatrick and Jarne 2000).  Wang et al. (1999) performed an extensive simulation study on the transient dynamics o f the genetic load for a population bottlenecked from an infinite initial size to a constant size (N = 50). For 2A chromosome sets containing 2900 mutable loci, new 7  mutations were each generated with a selection coefficient, s, and degree o f dominance, h, drawn from distributions chosen to reflect experimental evidence. For the case in which s was exponentially distributed and h was uniformly distributed between 0 and exp(-ks) 12, where A: is a constant (Wang et al. 1998), their results indicated that the genetic load increased at a decaying rate over time. Because the mean degree o f dominance decreased with increasing strength of selection, lethal alleles were highly recessive on average and were rapidly purged from the population. A s a result, slightly deleterious mutations became more  72  important for the long-term genetic load o f the small population. Given the negative correlation between s and h, however, highly deleterious mutations were more important in the short-term. Kirkpatrick and Jarne (2000) performed a similar numerical study in which a severe bottleneck from infinite size was immediately followed by rapid growth o f the population. The genetic load o f individual loci with specific values o f s and h was calculated for different bottleneck sizes. Purging o f recessive deleterious alleles caused the genetic load to temporarily drop below the equilibrium value before returning to that equilibrium as the population continued to grow.  There has not been, to my knowledge, a numerical study o f the evolution o f the genetic load in a population whose decline has occurred gradually. Consider a population o f current small size that has recently declined from a very large ancestral population, where the decline may have either been severe so that the population was historically small in size, or gradual so that it is only recently small in size. If completely recessive lethal alleles are rare, then the current genetic load is greater i f the population decline was severe so that deleterious alleles have had more opportunity to fix in the population. O n the other hand, i f completely recessive lethal alleles are abundant, then the current genetic load is greater i n a population that has declined gradually because the purging o f recessive lethals is delayed. Therefore, estimating the genetic load in a declining population requires knowledge o f the joint distribution o f 5 and h for mutations in the ancestral population. The general consensus is that most mutant alleles are slightly deleterious (Mukai 1964; Keightley 1995) with lethal mutations being comparatively rare. Because mostly recessive alleles make a disproportionately large contribution to the transient genetic load, however, the cumulative  73  effect o f demographic history w i l l depend more on the distribution o f selection coefficients within this dominance class. In reviewing the results o f several mutation accumulation experiments (e.g. Simmons and C r o w 1977), Caballero and Keightley (1994) observed that mutant alleles with large effects tended to recessivity, whereas alleles with small effects exhibited no tendency.  I have made several simplifying assumptions that may reduce the accuracy o f genetic load estimates obtained from the numerical model. Because the model addresses a single locus in isolation, there is no direct means to represent the effects o f linkage and epistasis. Linkage between deleterious alleles may have a substantial effect on the effectiveness o f purging (Charlesworth et al. 1992; Wang et al. 1999). Synergistic epistasis among deleterious alleles, such that the total effect exceeds the sum o f individual effects, w i l l cause purging to become more efficient (Wang et al. 1999). Furthermore, I have neglected the distinction between the census size and the effective size o f a population (Section 1.2). Finally, the variance i n genetic load w i l l cause empirical estimates to deviate from the mean values calculated by this model; this is the subject o f further study.  I know o f no experimental study on the genetic load in gradually declining populations. We know from mutation accumulation experiments such as M u k a i et al. (1972) that deleterious mutations are abundant, but less is known about how the genetic load caused by mutations depends on population size. Keightley et al. (1998) reviewed two contradictory experiments (Shabalina et al. 1997; Gilligan et al. 1997) i n which populations o f Drosophila recently caught from the w i l d were maintained as random-mating laboratory populations o f  74  constant size. Although fitness assays were performed to monitor the decline o f the populations relative to their controls, selection between families was minimized. Because this impedes the elimination o f deleterious alleles, the genetic load estimates were greater than would be obtained in equivalent natural populations. Shabalina et al. (1997) observed a decrease in fitness over time in populations held at N - 200 when fitness was assayed under harsh conditions. Gilligan et al. (1997), on the other hand, did not observe a lower fitness in their small population compared to their large populations or controls.  There is a rich empirical literature devoted to verifying the biological importance o f purging due to inbreeding. The recovery o f fitness in inbred lines was initially observed by Bowman and Falconer (1960) in mice and subsequently by Templeton and Read (1984) in Speke's gazelle, although the latter has recently been criticized for containing possible biases (Willis and Wiese 1997; Ballou 1997). Several experimental studies have since attempted to induce purging by artificially increasing the rate o f inbreeding in natural populations, with both affirmative (Carr and Dudash 1997; Barrett and Charlesworth 1991), negative (Dudash et al. 1997; M c C a l l et al. 1994), or mixed results (Lacy and Ballou 1998). Others have compared the magnitude o f inbreeding depression among lines derived from ancestral populations with a known history o f inbreeding (Latta and Ritland 1994; Weeks et al. 1999). The evidence tends to portray purging by inbreeding as a relatively weak process.  The empirical study o f purging by drift, on the other hand, is less extensive. Several studies (reviewed in Hauser et al. 1994) have contrasted the inbreeding depression i n large and small populations, but none o f them was able to detect a significant effect o f population  75  size. Brakefield and Saccheri (1994), however, observed a recovery o f fitness in laboratory populations o f satyrine butterflies (Bicyclus ariynand) founded from one or more pairs. The divergence o f morphological and molecular markers among bottlenecked lines was consistent with strong genetic drift. Lines founded from a single pair and allowed to increase to a carrying capacity o f K= 300 immediately experienced a dramatic (>50%) loss o f fecundity, followed by a gradual recovery in subsequent generations (Brakefield and Saccheri 1994). Similar results have been obtained by Bryant et al. (1990) for serially-bottlenecked populations o f house flies and by Fowler and Whitlock (1999) for populations o f D. melanogaster each founded by a single pair. Using a different experimental design, Ehiobu et al. (1989) adjusted the number o f generations in serially bottlenecked populations o f four and twenty D. melanogaster such that the same inbreeding coefficient was obtained by the final generation as caused by a single full-sib mating. A full-sib mating became a third treatment. The three treatments allowed them to distinguish between the effects o f drift and inbreeding on inbreeding depression. Ehiobu et al. (1989) observed a significantly greater larval viability i n their N=20 lines than the full-sib treatment, suggesting purging had occured. The N = 4 lines were not significantly different than the full-sib treatment, which may reflect a lesser opportunity for purging i n fewer generations or the fixation o f deleterious alleles by a greater intensity o f drift. These results, together with the results o f the numerical model described above, suggest that population that has declined gradually i n size may have a greater genetic load than a population that has declined severely.  76  CHAPTER  IV  CONCLUSIONS A N D RECOMMENDATION FOR FURTHER WORK  The evolutionary theory o f genetic load has developed over several decades, from the pioneering work o f Haldane (1937) and Muller (1950), into a rich literature where the influence of recombination (e.g. Otto and Feldman 1997), mating systems (e.g. Lande and Schemske 1985; Barrett and Charlesworth 1991), and mutation rate (Felsenstein 1974) on load have been considered. Application o f load theory to problems i n conservation biology has stimulated new interest into the effects o f drift on the genetic load (Lynch and Gabriel 1990; Lande 1994), as drift is increasingly important in smaller populations. In particular, the mutational meltdown process (Lynch and Gabriel 1990), by which the mean fitness o f a population is continually eroded by the fixation o f deleterious mutations by drift, is predicted to cause the extinction o f sexual populations as large as N= 10 (Lynch et al. 1995b). If mutations frequently compensate 3  for each other's effects, however, then the erosion o f fitness is halted as compensatory mutations become increasingly advantageous (Chap. II). The importance o f drift also emphasizes the undeveloped state o f theoretical models o f drift in relation to the genetic load (Kirkpatrick and Jarne 2000), especially for the transient dynamics following a change i n population size. Most of the studies on this matter have consisted o f numerical simulation and have calculated load for bottlenecked populations at a constant size (Wang et al. 1999) or undergoing rapid growth (Kirkpatrick and Jarne 2000). Exact numerical evaluation o f recursive formulae modeling the evolution o f allele frequency distributions in a declining population indicate that processes observed in the aforementioned studies (i.e. immediate boost in load, and subsequent purging) can be delayed, such that populations that have declined gradually may have a greater genetic load than others that have undergone a severe decline i n numbers (Chap. III).  77  Fisher's geometrical model o f adaptation (Fisher 1930), applied in Chapter II to an analysis o f compensatory mutation, is a good first step towards understanding the effects o f compensatory mutations on load, but more empirically grounded models should be explored next.  Although the compensatory and mutational meltdown models can produce qualitatively  different predictions on the extinction risk o f small populations, this divergence depends on parameters for which there is an insufficient amount o f empirical data. For example, the parameter n representing the number o f dimensions in the model is not a measurable biological quantity. Further theoretical work is required to derive models containing parameters that can be derived from empirical results. The parameter n is related to the probability that the next mutation w i l l be compensatory, which is a quantity that may be estimable from experimental data. Naturally, progress in this problem w i l l also require further experiments on the accumulation o f mutations in small populations (e.g. Gilligan et al. 1997). Further development of the analytic theory o f drift and genetic load is also required, especially for the load in declining populations.  78  BIBLIOGRAPHY  Abrarnowitz, M . , and I. A . Stegun. 1965. Handbook o f mathematical functions, with formulas, graphs, and mathematical tables. Dover Publications, N e w York. Ames, W . F. 1992. Numerical methods for partial differential equations. 3 ed. Academic Press, Boston. r d  Ballou, J. D . 1997. Ancestral Inbreeding Only M i n i m a l l y Affects Inbreeding Depression in Mammalian Populations. J. Hered. 88(3): 169-178. Barrett, S. C . H . , and D . Charlesworth. 1991. Effects o f a change i n the level o f inbreeding on the genetic load. Nature 352(6335): 522-524. Barton, N . H . , and M . Turelli. 1989. Evolutionary quantitative genetics: how little do we know? A n n . Rev. Genet. 23: 337-370. Bowman, J. C , and D . S. Falconer. 1960. Inbreeding depression and heterosis o f litter size i n mice. Genet. Res. 1: 262-274. Brakefield, P. M . , and I. J. Saccheri. 1994. Guidelines i n conservation genetics and the use o f the population cage experiments with butterflies to investigate the effects o f genetic drift and inbreeding. Pp. 165-179 i n V . Loeschcke, J. Tomiuk, and S. K . Jain, eds. Conservation Genetics. Birkhauser, Berlin. Bryant, E . H . , L . M . Meffert, and S. A . McCommas. 1990. Fitness rebound i n serially bottlenecked populations o f the house fly. A m . Nat. 136: 542-549. Burch, C . L . , and L . Chao. 1999. Evolution by small steps and rugged landscapes i n the R N A virus <P6. Genetics 151: 921-927. Burger, R . 1991. Moments, cumulants, and polygenic dynamics. J. Math. B i o l . 30: 199-213. Burger, R., and M . Lynch. 1995. Evolution and extinction in a changing environment: A quantitative-genetic analysis. Evolution 49(1): 151-163. Butcher, D . 1995. Muller's ratchet, epistasis and mutation effects. Genetics 141: 431-437. Caballero, A . , and P. D . Keightley. 1994. A pleiotropic nonadditive model o f variation i n quantitative traits. Genetics 138: 883-900. Carr, D . E . , and M . R. Dudash. 1997. The effects o f five generations o f enforced selfiing on potential male and female function i n Mimulusguttatus. Evolution 51(6): 1797-1807. Caughley, G . 1994. Directions i n conservation biology. J. A p p l . E c o l . 63: 215-244. Cedergren, R . J., D . Sankoff, B . LaRue, and H . Grosjean. 1981. The evolving t R N A molecule. Crit. Rev. Biochem. M o l . B i o l . 11(1): 35-104. 79  Charlesworth, B . 1980. Evolution in Age-Structured Populations. Cambridge University Press, Cambridge. Charlesworth, D . , and B . Charlesworth. 1987. Inbreeding depression and its evolutionary consequences. A n n . Rev. E c o l . Syst. 18: 237-268. ., M . T. Morgan, and B . Charlesworth. 1992. The effect o f linkage and population size on inbreeding depression due to mutational load. Genet. Res. 59: 49-61. Couvet, D . , and J. Ronfort. 1994. Mutation load depending on variance in reproductive success and mating system. E X S 68: 55-68. Clark, C . G . , B . W . Tague, V . C . Ware, and S. A . Gerbi. 1984. Xenopus laevis 28S ribosomal R N A : a secondary structure model and its evolutionary and functional implications. Nucleic Acids Res. 12(15): 6197-6220. Crow, J. F. 1948. Alternative hypotheses o f hybrid vigor. Genetics 33: 477-487. . 1970. Genetic Loads and the Cost o f Natural Selection. Pp. 128-177 i n K . Kojima, ed. Mathematical Topics in Population Genetics. Springer, N e w Y o r k . . 1997. The high spontaneous mutation rate: is it a health risk? Proc. Natl. Acad. Sci. U S A 94: 8380-8386. Crow, J. F., and M . Kimura. 1970. An Introduction to Population Genetics Theory. Burgess Publishing Company. Minneapolis. Dickinson, R. J., M . A . Sobanski and M . J. E . Hewlins. 1995. In Saccharomyces cerevisiae deletion o f phosphoglucose isomerase can be suppressed by increased activities o f enzymes o f the hexose monophosphate pathway. Microbiology 141(2): 385-391. Drake, J. W . , B . Charlesworth, D . Charlesworth, and J. F . Crow. 1998. Rates o f spontaneous mutation. Genetics 148(4): 1667-1686. Dobzhansky, T. 1957. Genetic loads in natural populations. Science 126: 191-194. Dudash, M . R., and D . E . Carr. 1998. Genetics underlying inbreeding depression i n Mimulus with contrasting mating systems. Nature 393: 682-684. ., ., and C . B . Fenster. 1997. Five generations o f enforced selfing and outcrossing in Mimulus guttatus: inbreeding depression variation at the population and family level. Evolution 51: 54-65. E l Mezaine, A . , S. K . Lehtinen, N . Hance, L . G . Nijtmans, D . Dunbar, I. J. Holt, and H . T. Jacobs. 1998. A t R N A suppressor mutation in human mitochondria. Nature Genet. 18(4): 350353.  80  Elena, S., M . Davila, I. S. Novella, J. J. Holland, E . Domingo, and A . M o y a . 1998. Evolutionary dynamics o f fitness recovery from the debilitating effects o f Muller's ratchet. Evolution 52: 309-314. Ehrlich, P. R., and A . H . Ehrlich. 1981. Extinction: the causes and consequences o f the disappearance o f species. Random House, N e w York. Ewens, W . J. 1965. The adequacy o f the diffusion approximation to certain distributions in genetics. Biometrics 21: 386-394. . 1979. Mathematical population genetics. Springer, N e w York. Eyre-Walker, A . , and P. D . Keightley. 1999. H i g h genomic deleterious mutation rates in hominids. Nature 397: 344-347. Falconer, D . S. 1989. Introduction to Quantitative Genetics. 3 Technical, Hong K o n g .  rd  ed. Longman Scientific &  Felsenstein, J. 1974. The evolutionary advantage o f recombination. Genetics 78: 737-756. Feller, W . 1951. Diffusion processes in genetics. Pp. 227'-246 in Proceedings oj the Second Berkeley Symposium on Mathematical Statistics and Probability. University o f California Press, Berkeley. . 1968. An Introduction to Probability Theory and Its Applications, Vol. 1. 3 Wiley & Sons, Inc., N e w York.  rd  ed. John  Fisher, R. A . 1930. Genetical Theory of Natural Selection. The Clarendon press, Oxford. Fowler, K . , and M . C . Whitlock. 1999. The distribution o f phenotypic variance with inbreeding. Evolution 53(4): 1143-1156. Frankham, R. 1995. Conservation genetics. Annu. Rev. Genet. 29: 305-327. Frankel, O. H . , and M . E . Soule. 1981. Conservation and evolution. Cambridge University Press, Cambridge. Fry, J. D . , P. D . Keightley, S. L . Heinsohn, and S. V . Nuzhdin. 1999. N e w estimates o f the rates and effects o f mildly deleterious mutation in Drosophila melanogaster. Gabriel, W . , M . Lynch, and R. Burger. 1993. Muller's ratchet and mutational meltdowns. Evolution 47: 1744-1757. Gachotte, D . , C . A . Pierson, N . D . Lees, R. Barbuch, C . Koegel, and M . Bard. 1997. A yeast sterol auxotroph (erg25) is rescued by addition o f azole antifungals and reduced levels o f heme. Proc. Natl. Acad. Sci. U S A 94: 11173-11178. Gillespie, D . T. 1992. Markov processes : an introduction for physical scientists. Academic Press, Boston. 81  Gilligan, D . M . , L . M . Woodworm, M . E . Montgomery, D . A . Briscoe, and R. Frankham. 1997. Is mutation accumulation a threat to the survival o f endangered populations? Cons. B i o l . 11(5): 1235-1241. Garcia-Dorado, A . 1997. The rate and effects distribution o f viability mutation in Drosophila: M i n i m u m distance estimation. Evolution 51 (4): 1130-1139. Gomulkiewicz, R., R. D . Holt, and M . Barfield. 1999. The effects o f density dependence and immigration on local adaptation and niche evolution in a black-hole sink environment. Theor. Popul. B i o l . 55(3): 283-296. Haigh, J. 1978. The accumulation o f deleterious genes i n a population — Muller's Ratchet. Theor. Popul. B i o l . 14: 251-267. Haldane, J. B . S. 1937. The effect o f variation on fitness. A m . Nat. 71: 337-349. . 1957. The cost o f natural selection. J. Genet. 55: 511-524. Hancock, J. M . , D . Tautz, and G . A . Dover. 1988. Evolution o f the secondary structures and compensatory mutations o f the ribosomal R N A s o f Drosophila melanogaster. M o l . B i o l . E v o l . 5: 393-414. Hanson, D . K . , D . M . Tiede, S. L . Nance, C - H . Chang, and M . Schiffer. 1993. Site-specific and compensatory mutations imply unexpected pathways for proton delivery to the Q - B binding site of the photosynthetic reaction center. Proc. Natl. Acad. Sci. U S A 90(19): 8929-8933. Hartl, D . L . , and C . H . Taubes. 1996. Compensatory nearly neutral mutations: selection without adaptation. J. Theor. B i o l . 182: 303-309. . 1998. Towards a theory o f evolutionary adaptation. Genetica 103: 525-533. Hartman, P. E . , and J. R. Roth. 1973. Mechanisms o f suppression. A d v . Genet. 17: 1-105. Hastings, A . 1997. Population biology: concepts and models. Springer, N e w York. Hauser, T. P., C . Damgaard, and V . Loeschcke. Effects o f inbreeding i n small plant populations: Expectations and implications for conservation. Pp. 115-129 in V . Loeschcke, J. Tomiuk, and S. K . Jain, eds. Conservation Genetics. Birkhauser, Berlin. Hedrick, P. W . 1994. Purging inbreeding depression and the probability o f extinction: full-sib mating. Heredity 73: 363-372. ., and P. S. Miller. 1992. Conservation genetics: techniques and fundamentals. A p p l . 2(1): 30-46.  Ecol.  Higgs, P. G . 1998. Compensatory neutral mutations and the evolution o f R N A . Genetica 103: 91-101.  82  I  Holt, R. D . 1997. O n the evolutionary stability o f sink populations. E v o l . E c o l . 11(6): 723731. Hou, Y . M . , and P. Schimmel. 1992. Functional compensation o f a recognition-defective transfer R N A by a distal base pair substitution. Biochemistry 31: 10310-10314. Jarvik, J., and D . Botstein. 1975. Conditional-lethal mutations that suppress genetic defects in morphogenesis by altering structural proteins. Proc. Natl. A c a d . Sci. U S A 72(7): 2738-2742. Jucovic, M . , and A . R. Poteete. 1998. Delineation o f an evolutionary salvage pathway by compensatory mutations o f a defective lysozyme. Prot. S c i . 7: 2200-2209. Keightley, P. D . 1994. The distribution o f mutation effects on viability in Drosophila melanogaster. Genetics 138: 1315-1322. ., A . Caballero, and A . Garcia-Dorado. 1998. Population genetics: Surviving under mutation pressure. Curr. B i o l . 8: R235-R237. Kendall, M . G . , and A . Stuart. 1963. The Advanced Theory o f Statistics. Charles Griffin & Company Limited, London. K i m , H . - W . , T.-J. Shen, D . P. Sun, N . T. H o , M . Madrid, M . F. Tam, M . Z o u , P. F. Cottam, and C. H o . 1994. Restoring allosterism with compensatory mutations in hemoglobin. Proc. Natl. Acad. Sci. U S A 91: 11547-11551. Kimura, M . 1955. Solution o f a Process o f Random Genetic Drift with a Continuous Model. Proc. Natl. A c a d . Sci. U S A 41: 144-150. . 1964. Diffusion Models in Population Genetics. J. A p p l . Prob. 1: 177-232. . 1965. A stochastic model concerning the maintenance o f genetic variability in quantitative characters. Proc. Natl. Acad. Sci. U S A 54:731-736. . 1983. The neutral theory o f molecular evolution. Cambridge University Press, Cambridge [Cambridgeshire] ; N e w York. . 1990. Some models o f neutral evolution, compensatory evolution, and the shifting balance process. Theor. Popl. B i o l . 37: 150-158. ., T. Maruyama, and J. F. Crow. 1963. The mutation load in small populations. Genetics 48: 1303-1312. Kirkpatrick, M . , and P. Jarne. 2000. The Effects o f a Bottleneck on the Inbreeding Depression and the Genetic Load. A m . Nat. 155(2): 154-167. Kondrashov, A . S. 1994. Muller's ratchet under epistatic selection. Genetics 136: 1469-1473. . 1995. Contamination o f the genome by very slightly deleterious mutations: why have we not died 100 times over? J. Theor. B i o l . 175: 583-594. 83  Lacy, R. C . and J. D . Ballou. 1998. Effectiveness o f selection i n reducing the genetic load in populations o f Peromyscus polionotus during generations o f inbreeding. Lande, R. 1988. Genetics and demography in biological conservation. Science 241: 14551460. . 1993. Risks o f population extinction from demographic and environmental stochasticity and random catastrophes. A m . Nat. 142: 911-927. —  . 1994. Risk o f population extinction from fixation o f new deleterious mutations.  Evolution 48: 1460-1469. . 1995. Mutation and conservation. Cons. B i o l . 9: 782-791. . 1998. Risk o f population extinction from fixation o f deleterious and reverse mutations. Genetica 103: 21-27. ., and D . W . Schemske. 1985. The evolution o f self-fertilization and inbreeding depression i n plants. I. Genetic models. Evolution 39: 24-40. Latta, R., and K . Ritland. 1994. The relationship between inbreeding depression and prior inbreeding among populations o f four Mimulus taxa. Evolution 48(3): 806-817. Latter, B . D . H . , J. C . Mulley, D . Reid, and L . Pascoe. 1995. Reduced Genetic Load Revealed by Slow Inbreeding i n Drosophila melanogaster. Genetics 139: 287-298. Leigh, E . G . 1987. Ronald Fisher and the development o f evolutionary theory. II. Influences o f new variation on evolutionary process. Oxford University Press, Oxford. Lenormand, T., T. Guillemaud, D . Bourguet, and M . Raymond. 1998. Appearance and sweep of a gene duplication: Adaptive response and potential for new functions i n the mosquito Culex pipiens. Evolution 52: 1705-1712. Lerner, I. M . 1954. Genetic homeostasis. Oliver and Boyd. Edinburgh. Liang, C , L . Rong, M . Laughrea, L . Kleiman, and M . A . Wainberg. 1998. Compensatory point mutations i n the human immunodeficiency virus type 1 G a g region that are distal from deletion mutations i n the dimerization initiation site can restore viral replication. J. V i r o l . 72: 66296636.  Lighthill, M . J. 1958. An introduction to Fourier analysis and generalised functions. Cambridge University Press, Cambridge. Lyman, R . F., F. Lawrence, S. V . Nuzhdin, and T. F . C . Mackay. 1996. Effects o f Single PElement Insertions on Bristle Number and Viability i n Drosophila melanogaster. Genetics 143: 277-292.  84  L y n c h , M . , R. Burger, D . Butcher, and W . Gabriel. 1993. The mutational meltdown in asexual populations. J. Hered. 84: 339-344. ., J. Conery, and R. Burger. 1995a. Mutation accumulation and the extinction o f small populations. A m . Nat. 146:489-518. ., J. Conery, and R. Burger. 1995b. Mutational meltdowns in sexual populations. Evolution 49: 1067-1080. ., and W . Gabriel. 1990. Mutation Load and the Survival o f Small Populations. Evolution 44: 1725-1737. Mackay, T. F . C , R. F. Lyman, and M . S. Jackson. 1992. Effects o f P Element Insertions on Quantitative Traits in Drosophila melanogaster. Genetics 130: 315-332. Manning, K . , M . A l - D h a l i m y , M . Finegold, and M . Grompe. 1999. In vivo suppressor mutations correct a murine model o f hereditary tyrosinemia type I. Proc. A c a d . Natl. Sci. U S A 96(21): 11928-11933. Maringanti, S., and J. A . Imlay. 1999. A n intracellular iron chelator pleiotropically suppresses enzymatic and growth defects o f superoxide dismutase-deficient Escherichia coli. J. Bacteriol. 181(12): 3792-3802. Maruyama, T., and P. A . Fuerst. 1985. Population bottlenecks and nonequilibrium models in population genetics. II. Number o f alleles in a small populations that was formed by a recent bottleneck. Genetics 111: 675-689. Mateu, M . G . , and A . R. Fersht. 1999. Mutually compensatory mutations during evolution o f the tetramerization domain o f tumor suppressor p53 lead to impaired hetero-oligomerization. Proc. Natl. A c a d . Sci. U S A 96: 3595-3599. Matsuno K . , and A . L . Sonenshein. 1999. Role o f S p o V G i n asymmetric septation i n Bacillus subtilis. J. Bacteriol. 181(11): 3392-3401. M c C a l l , C , D . M . Waller, and T. Mitchell-Olds. 1994. Effects o f serial inbreeding on fitness components i n Impatiens capensis. Evolution 48: 818-827. Morton, N . E . , J. F. Crow, and H . J. Muller. 1956. A n estimate o f the mutational damage in man from data on consanguineous marriages. Proc. Natl. A c a d . Sci. U S A 42: 855-863. M u k a i , T. 1964. The genetic structure o f natural populations o f Drosophila melanogaster. I. Spontaneous mutation rate o f polygenes controlling viability. Genetics 50: 1-19. , S. I. Chigusa, L . E . Mettler, and J. F. Crow. 1972. Mutation rate and dominance o f genes affecting viability i n Drosophila melanogaster. Genetics 72: 335-355. Muller, H . J. 1938. Reversibility in evolution considered from the standpoint o f genetics. B i o l . Rev. Camb. Philos. Soc. 14: 261-280. 85  . 1950. Our load o f mutations. Amer. J. H u m . Genet. 2: 111-176. . 1964. The relation o f recombination to mutational advance. Mutat. Res. 1: 2-9. Murgola, E . J. 1985. t R N A , suppression, and the code. A n n . Rev. Genet. 19: 57-80. Nei, M . , T. Maruyama, and R. Chakrabotry. 1975. The bottleneck effect and genetic variability in populations. Evolution 29: 1-10. Nuzhdin, S. V . , and T. F. Mackay. 1995. The genomic rate o f transposable element movement in Drosophila melanogaster [letter]. M o l . B i o l . E v o l . 12(1): 180-181. Ohta, T. 1977. Extensions to the neutral mutation random drift hypothesis. Pp. 148-167 in M . Kimura, ed. Molecular evolution and polymorphism. National Institute o f Genetics, Mishima, Japan. . 1992. The nearly neutral theory o f molecular evolution. A n n . Rev. E c o l . Syst. 23: 263-286. . 1998. Evolution by nearly-neutral mutations. Genetica 103: 83-90. Otto, S. P., and M . W . Feldman. 1997. Deleterious mutations, variable epistatic interactions, and the evolution o f recombination. Theor. Popul. B i o l . 51(2): 134-147. Otto, S. P., and M . C . Whitlock. 1997. The Probability o f Fixation in Populations o f Changing Size. Genetics 146: 723-733. Orr, H . A . 1998. The population genetics o f adaptation: The distribution o f factors fixed during adaptive evolution. Evolution 52: 935-949. Orr, H . A . 2000. Adaptation and the cost o f complexity. Evolution (in press). Peck, J. R., G . Barreau, and S. C . Heath. 1997. Imperfect genes, Fisherian mutation and the evolution o f sex. Genetics 145: 1171-1199. Press, W . H . 1992. Numerical recipes in C : the art o f scientific computing. Cambridge University Press, N e w York. Rice, S. 1990. A geometric model for the evolution o f development. J. Theor. B i o l . 143: 319342. Robertson, A . 1952. The effects o f inbreeding on the variation due to recessive genes. Genetics 37: 187-207. . 1970. The reduction in fitness from genetic drift at heterotic loci i n small populations. Genet. Res., Camb. 15: 257-259.  86  Schultz, S. T., and M . Lynch. 1997. Mutation and extinction: The role o f variable mutational effects, synergistic epistasis, beneficial mutations, and degree o f outcrossing. Evolution 51: 1363-1371. Shabalina, S. A . , L . Y . Yampolsky, and A . S. Kondrashov. 1997. Rapid decline o f fitness in panmictic populations o f Drosophila melanogaster maintained under relaxed natural selection. Proc. Natl. A c a d . Sci. U S A 94: 13034-13039. Simmons, M . J., and J. F. Crow. 1977. Mutations affecting fitness in Drosophila populations. Annu. Rev. Genet. 11: 49-78. Springer, M . S., L . J. Hollar, and A . Burk. 1995. Compensatory substitutions and the evolution of the mitochondrial 12S r R N A gene in mammals. M o l . B i o l . E v o l . 12(6): 1138-1150. Steinberg, S., and R. Cedergren. 1994. Structural compensation in atypical mitochondrial t R N A s . Nature Struct. B i o l . 1: 507-510. Stephan, W . 1996. The rate o f compensatory evolution. Genetics 144(1): 419-426. ., and D . A . Kirby. 1993. R N A folding i n Drosophila shows a distance effect for compensatory fitness interactions. Genetics 135(1): 97-103. Templeton, A . R., and B . Read. 1964. Factors eliminating inbreeding depression i n a captive herd o f Speke's gazelle. Zoo B i o l . 3: 177-199. Tyvand, P . A . A n Exact Algebraic Theory o f Genetic Drift in Finite D i p l o i d Populations with Random Mating. J. Theor. B i o l . 163: 315-331. van Kampen, N . G . 1976. The expansion o f the master equation. A d v . Chem. Phys. 34: 245309. . 1992. Stochastic processes in physics and chemistry. Elsevier Science Publishers B . V . , Amsterdam. Wagner, G . P., and W . Gabriel. 1990. Quantitative Variation in Finite Parthenogenetic Populations What Stops Muller's Ratchet in the Absence o f Recombination? Evolution 44: 715731. Wang., J., A . Caballero, P. D . Keightley, and W . G . H i l l . 1998. Bottleneck Effect on Genetic Variance: A Theoretical Investigation o f the Role o f Dominance. Genetics 150: 435-447. ., W . G . H i l l , D . Charlesworth, and B . Charlesworth. 1999. Dynamics o f inbreeding depression due to deleterious mutations in small populations: mutation parameters and inbreeding rate. Genet. Res., Camb. 74: 165-178. Waterston, R. H . , and S. Brenner. 1978. A suppressor mutation i n the nematode acting on specific alleles o f many genes. Nature 275: 715-719.  87  Weeks, S. C , V . Marcus, and B . R. Crosser. 1999. Inbreeding depression in a self-compatible, androdioecious crustacean, Eulimnadia texana. Evolution 53(2): 472-483. W i l l i s , K . , and R. Wiese. 1997. Elimination o f inbreeding depression from captive populations: Speke's gazelle revisited. Zoo B i o l . 16: 9-16. Wilson, E . O. 1988. Biodiversity. National Academy Press, Washington, D . C . Wright, S. 1922. Coefficients o f inbreeding and relationship. A m . Nat 56: 330-338. . 1931. Evolution in Mendelian populations. Genetics 16: 97-159. . 1937. The distribution o f gene frequencies in populations. Proc. Natl. Acad. Sci. U S A 23: 307-320. . 1945. The differential equation o f the distribution o f gene frequencies. Proc. Natl. Acad. Sci. U S A 31: 382-389.  . 1969. Evolution and the Genetics of Populations. Vol.2. The Theory of Gene Frequencies. University o f Chicago Press, Chicago. Y u , W . , and R. J. Spreitzer. 1992. Chloroplast heteroplasmicity is stabilized by an ambersuppressor tryptophan t R N A C U A . Proc. Natl. A c a d . S c i . U S A 89: 3904-3907.  88  

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0089617/manifest

Comment

Related Items