Genetics Population and evolutionary genetics

Estimating Divergence Parameters With Small Samples From a Large Number of Loci [Population and evolutionary genetics]

Wang, Y., Hey, J. — Tue, 26 Jan 2010 14:34:08 PST

Most methods for studying divergence with gene flow rely upon data from many individuals at few loci. Such data can be useful for inferring recent population history but they are unlikely to contain sufficient information about older events. However, the growing availability of genome sequences suggests a different kind of sampling scheme, one that may be more suited to studying relatively ancient divergence. Data sets extracted from whole-genome alignments may represent very few individuals but contain a very large number of loci. To take advantage of such data we developed a new maximum-likelihood method for genomic data under the isolation-with-migration model. Unlike many coalescent-based likelihood methods, our method does not rely on Monte Carlo sampling of genealogies, but rather provides a precise calculation of the likelihood by numerical integration over all genealogies. We demonstrate that the method works well on simulated data sets. We also consider two models for accommodating mutation rate variation among loci and find that the model that treats mutation rates as random variables leads to better estimates. We applied the method to the divergence of Drosophila melanogaster and D. simulans and detected a low, but statistically significant, signal of gene flow from D. simulans to D. melanogaster.

Mating-System Variation, Demographic History and Patterns of Nucleotide Diversity in the Tristylous Plant Eichhornia paniculata [Population and evolutionary genetics]

Ness, R. W., Wright, S. I., Barrett, S. C. H. — Tue, 26 Jan 2010 14:34:08 PST

Inbreeding in highly selfing populations reduces effective size and, combined with demographic conditions associated with selfing, this can erode genetic diversity and increase population differentiation. Here we investigate the role that variation in mating patterns and demographic history play in shaping the distribution of nucleotide variation within and among populations of the annual neotropical colonizing plant Eichhornia paniculata, a species with wide variation in selfing rates. We sequenced 10 EST-derived nuclear loci in 225 individuals from 25 populations sampled from much of the geographic range and used coalescent simulations to investigate demographic history. Highly selfing populations exhibited moderate reductions in diversity but there was no significant difference in variation between outcrossing and mixed mating populations. Population size interacted strongly with mating system and explained more of the variation in diversity within populations. Bayesian structure analysis revealed strong regional clustering and selfing populations were highly differentiated on the basis of an analysis of F_st. There was no evidence for a significant loss of within-locus linkage disequilibrium within populations, but regional samples revealed greater breakdown in Brazil than in selfing populations from the Caribbean. Coalescent simulations indicate a moderate bottleneck associated with colonization of the Caribbean from Brazil ~125,000 years before the present. Our results suggest that the recent multiple origins of selfing in E. paniculata from diverse outcrossing populations result in higher diversity than expected under long-term equilibrium.

The Spontaneous Appearance Rate of the Yeast Prion [PSI+] and Its Implications for the Evolution of the Evolvability Properties of the [PSI+] System [Population and evolutionary genetics]

Lancaster, A. K., Bardill, J. P., True, H. L., Masel, J. — Tue, 26 Jan 2010 14:34:08 PST

Epigenetically inherited aggregates of the yeast prion [PSI+] cause genomewide readthrough translation that sometimes increases evolvability in certain harsh environments. The effects of natural selection on modifiers of [PSI+] appearance have been the subject of much debate. It seems likely that [PSI+] would be at least mildly deleterious in most environments, but this may be counteracted by its evolvability properties on rare occasions. Indirect selection on modifiers of [PSI+] is predicted to depend primarily on the spontaneous [PSI+] appearance rate, but this critical parameter has not previously been adequately measured. Here we measure this epimutation rate accurately and precisely as 5.8 x 10^–7 per generation, using a fluctuation test. We also determine that genetic "mimics" of [PSI+] account for up to 80% of all phenotypes involving general nonsense suppression. Using previously developed mathematical models, we can now infer that even in the absence of opportunities for adaptation, modifiers of [PSI+] are only weakly deleterious relative to genetic drift. If we assume that the spontaneous [PSI+] appearance rate is at its evolutionary optimum, then opportunities for adaptation are inferred to be rare, such that the [PSI+] system is favored only very weakly overall. But when we account for the observed increase in the [PSI+] appearance rate in response to stress, we infer much higher overall selection in favor of [PSI+] modifiers, suggesting that [PSI+]-forming ability may be a consequence of selection for evolvability.

The Genetics of Postmating, Prezygotic Reproductive Isolation Between Drosophila virilis and D. americana [Population and evolutionary genetics]

Sweigart, A. L. — Tue, 26 Jan 2010 14:34:08 PST

Many studies have demonstrated the rapid diversification of reproductive genes that function after mating but before fertilization. This process might lead to the evolution of postmating, prezygotic barriers between species. Here, I investigate the phenotypic and genetic basis of postmating, prezygotic isolation between two closely related species of Drosophila, Drosophila virilis and D. americana. I show that a strong barrier to interspecific fertilization results in a 99% reduction in progeny production. A genetic interaction among maternal and paternal alleles at only a few loci prevents the fertilization of D. virilis females by D. americana males. These loci are autosomal and isolation acts recessively; the fertilization incompatibility is caused by at least two loci in the maternal D. virilis parent in combination with at least three loci in the paternal D. americana parent. These findings, together with results from classical experiments, suggest that male–female coevolution within D. americana may have driven postmating, prezygotic isolation between species.

Molecular and Genetic Analyses of Four Nonfunctional S Haplotype Variants Derived from a Common Ancestral S Haplotype Identified in Sour Cherry (Prunus cerasus L.) [Population and evolutionary genetics]

Tsukamoto, T., Hauck, N. R., Tao, R., Jiang, N., Iezzoni, A. F. — Tue, 26 Jan 2010 14:34:09 PST

Tetraploid sour cherry (Prunus cerasus) has an S-RNase-based gametophytic self-incompatibility (GSI) system; however, individuals can be either self-incompatible (SI) or self-compatible (SC). Unlike the situation in the Solanaceae, where self-compatibility accompanying polyploidization is often due to the compatibility of heteroallelic pollen, the genotype-dependent loss of SI in sour cherry is due to the compatibility of pollen containing two nonfunctional S haplotypes. Sour cherry individuals with the S₄S₆S_36aS_36b genotype are predicted to be SC, as only pollen containing both nonfunctional S_36a and S_36b haplotypes would be SC. However, we previously found that individuals of this genotype were SI. Here we describe four nonfunctional S₃₆ variants. Our molecular analyses identified a mutation that would confer loss of stylar S function for one of the variants, and two alterations that might cause loss of pollen S function for all four variants. Genetic crosses showed that individuals possessing two nonfunctional S₃₆ haplotypes and two functional S haplotypes have reduced self-fertilization due to a very low frequency of transmission of the one pollen type that would be SC. Our finding that the underlying mechanism limiting successful transmission of genetically compatible gametes does not involve GSI is consistent with our previous genetic model for Prunus in which heteroallelic pollen is incompatible. This provides a unique case in which breakdown of SI does not occur despite the potential to generate SC pollen genotypes.

Coalescent Simulation of Intracodon Recombination [Population and evolutionary genetics]

Arenas, M., Posada, D. — Tue, 26 Jan 2010 14:34:09 PST

The coalescent with recombination is a very useful tool in molecular population genetics. Under this framework, genealogies often represent the evolution of the substitution unit, and because of this, the few coalescent algorithms implemented for the simulation of coding sequences force recombination to occur only between codons. However, it is clear that recombination is expected to occur most often within codons. Here we have developed an algorithm that can evolve coding sequences under an ancestral recombination graph that represents the genealogies at each nucleotide site, thereby allowing for intracodon recombination. The algorithm is a modification of Hudson's coalescent in which, in addition to keeping track of events occurring in the ancestral material that reaches the sample, we need to keep track of events occurring in ancestral material that does not reach the sample but that is produced by intracodon recombination. We are able to show that at typical substitution rates the number of nonsynonymous changes induced by intracodon recombination is small and that intracodon recombination does not generally result in inflated estimates of the overall nonsynonymous/synonymous substitution ratio (). On the other hand, recombination can bias the estimation of at particular codons, resulting in apparent rate variation among sites and in the spurious identification of positively selected sites. Importantly, in this case, allowing for variable synonymous rates across sites greatly reduces the false-positive rate and recovers statistical power. Finally, coalescent simulations with intracodon recombination could be used to better represent the evolution of nuclear coding genes or fast-evolving pathogens such as HIV-1.We have implemented this algorithm in a computer program called NetRecodon, freely available at http://darwin.uvigo.es.

Ploidy and the Evolution of Endosperm of Flowering Plants [Population and evolutionary genetics]

Cailleau, A., Cheptou, P.-O., Lenormand, T. — Tue, 26 Jan 2010 14:34:09 PST

In angiosperms, spermatozoa go by pair in each pollen grain and fertilize, in addition to the egg cell, one of its sister cells, called the central cell. This "double fertilization" leads to the embryo on the one hand and to its nutritive tissue, the endosperm, on the other hand. In addition, in most flowering plants, the endosperm is triploid because of a doubled maternal genetic contribution in the central cell. Most of the hypotheses trying to explain these eccentricities rest on the assumption of a male/female conflict over seed resource allocation. We investigate an alternative hypothesis on the basis of the masking of deleterious alleles. Using analytical methods, we show that a doubled maternal contribution and double fertilization tend to be favored in a wide range of conditions when deleterious mutations alter the function of the endosperm. Furthermore, we show that these conditions vary depending on whether these traits are under male or female control, which allows us to describe a new type of male/female conflict.

The Nuclear Component of a Cytonuclear Hybrid Incompatibility in Mimulus Maps to a Cluster of Pentatricopeptide Repeat Genes [Population and evolutionary genetics]

Barr, C. M., Fishman, L. — Tue, 26 Jan 2010 14:34:09 PST

Characterizing the genetic and molecular basis of hybrid incompatibilities is a first step toward understanding their evolutionary origins. We fine mapped the nuclear restorer (Rf) of cytoplasm-dependent anther sterility in Mimulus hybrids by identifying and targeting regions of the Mimulus guttatus genome containing large numbers of candidate pentatricopeptide repeat genes (PPRs). The single Mendelian locus Rf was first isolated to a 1.3-cM region on linkage group 7 that spans the genome's largest cluster of PPRs, then split into two tightly linked loci (Rf1 and Rf2) by <10 recombination events in a large (N = 6153) fine-mapping population. Progeny testing of fertile recombinants demonstrated that a dominant M. guttatus allele at each Rf locus was sufficient to restore fertility. Each Rf locus spans a physical region containing numerous PPRs with high homology to each other, suggesting recent tandem duplication or transposition. Furthermore, these PPRs have higher homology to restorers in distantly related taxa (petunia and rice) than to PPRs elsewhere in the Mimulus genome. These results suggest that the cytoplasmic male sterility (CMS)–PPR interaction is highly conserved across flowering plants. In addition, given our theoretical understanding of cytonuclear coevolution, the finding that hybrid CMS results from interactions between a chimeric mitochondrial transcript that is modified by Rf loci identified as PPRs is consistent with a history of selfish mitochondrial evolution and compensatory nuclear coevolution within M. guttatus.

Rate of Adaptation in Large Sexual Populations [Population and evolutionary genetics]

Neher, R. A., Shraiman, B. I., Fisher, D. S. — Tue, 26 Jan 2010 14:34:09 PST

Adaptation often involves the acquisition of a large number of genomic changes that arise as mutations in single individuals. In asexual populations, combinations of mutations can fix only when they arise in the same lineage, but for populations in which genetic information is exchanged, beneficial mutations can arise in different individuals and be combined later. In large populations, when the product of the population size N and the total beneficial mutation rate U_b is large, many new beneficial alleles can be segregating in the population simultaneously. We calculate the rate of adaptation, v, in several models of such sexual populations and show that v is linear in NU_b only in sufficiently small populations. In large populations, v increases much more slowly as log NU_b. The prefactor of this logarithm, however, increases as the square of the recombination rate. This acceleration of adaptation by recombination implies a strong evolutionary advantage of sex.

The Evolution of Control and Distribution of Adaptive Mutations in a Metabolic Pathway [Population and evolutionary genetics]

Wright, K. M., Rausher, M. D. — Tue, 26 Jan 2010 14:34:09 PST

In an attempt to understand whether it should be expected that some genes tend to be used disproportionately often by natural selection, we investigated two related phenomena: the evolution of flux control among enzymes in a metabolic pathway and properties of adaptive substitutions in pathway enzymes. These two phenomena are related by the principle that adaptive substitutions should occur more frequently in enzymes with greater flux control. Predicting which enzymes will be preferentially involved in adaptive evolution thus requires an evolutionary theory of flux control. We investigated the evolution of enzyme control in metabolic pathways with two models of enzyme kinetics: metabolic control theory (MCT) and Michaelis–Menten saturation kinetics (SK). Our models generate two main predictions for pathways in which reactions are moderately to highly irreversible: (1) flux control will evolve to be highly unequal among enzymes in a pathway and (2) upstream enzymes evolve a greater control coefficient then those downstream. This results in upstream enzymes fixing the majority of beneficial mutations during adaptive evolution. Once the population has reached high fitness, the trend is reversed, with the majority of neutral/slightly deleterious mutations occurring in downstream enzymes. These patterns are the result of three factors (the first of these is unique to the MCT simulations while the other two seem to be general properties of the metabolic pathways): (1) the majority of randomly selected, starting combinations of enzyme kinetic rates generate pathways that possess greater control for the upstream enzymes compared to downstream enzymes; (2) selection against large pools of intermediate substrates tends to prevent majority control by downstream enzymes; and (3) equivalent mutations in enzyme kinetic rates have the greatest effect on flux for enzymes with high levels of flux control, and these enzymes will accumulate adaptive substitutions, strengthening their control. Prediction 1 is well supported by available data on control coefficients. Data for evaluating prediction 2 are sparse but not inconsistent with this prediction.

Regulation of Epithelial Stem Cell Replacement and Follicle Formation in the Drosophila Ovary [Population and evolutionary genetics]

Nystul, T., Spradling, A. — Tue, 26 Jan 2010 14:34:09 PST

Though much has been learned about the process of ovarian follicle maturation through studies of oogenesis in both vertebrate and invertebrate systems, less is known about how follicles form initially. In Drosophila, two somatic follicle stem cells (FSCs) in each ovariole give rise to all polar cells, stalk cells, and main body cells needed to form each follicle. We show that one daughter from each FSC founds most follicles but that cell type specification is independent of cell lineage, in contrast to previous claims of an early polar/stalk lineage restriction. Instead, key intercellular signals begin early and guide cell behavior. An initial Notch signal from germ cells is required for FSC daughters to migrate across the ovariole and on occasion to replace the opposite stem cell. Both anterior and posterior polar cells arise in region 2b at a time when ~16 cells surround the cyst. Later, during budding, stalk cells and additional polar cells are specified in a process that frequently transfers posterior follicle cells onto the anterior surface of the next older follicle. These studies provide new insight into the mechanisms that underlie stem cell replacement and follicle formation during Drosophila oogenesis.

The Power of the Methods for Detecting Interlocus Gene Conversion [Population and evolutionary genetics]

Mansai, S. P., Innan, H. — Tue, 26 Jan 2010 14:34:09 PST

Interlocus gene conversion can homogenize DNA sequences of duplicated regions with high homology. Such nonvertical events sometimes cause a misleading evolutionary interpretation of data when the effect of gene conversion is ignored. To avoid this problem, it is crucial to test the data for the presence of gene conversion. Here, we performed extensive simulations to compare four major methods to detect gene conversion. One might expect that the power increases with increase of the gene conversion rate. However, we found this is true for only two methods. For the other two, limited power is expected when gene conversion is too frequent. We suggest using multiple methods to minimize the chance of missing the footprint of gene conversion.

Gene Genealogies Strongly Distorted by Weakly Interfering Mutations in Constant Environments [Population and evolutionary genetics]

Tue, 26 Jan 2010 14:34:09 PST

Neutral nucleotide diversity does not scale with population size as expected, and this "paradox of variation" is especially severe for animal mitochondria. Adaptive selective sweeps are often proposed as a major cause, but a plausible alternative is selection against large numbers of weakly deleterious mutations subject to Hill–Robertson interference. The mitochondrial genealogies of several species of whale lice (Amphipoda: Cyamus) are consistently too short relative to neutral-theory expectations, and they are also distorted in shape (branch-length proportions) and topology (relative sister-clade sizes). This pattern is not easily explained by adaptive sweeps or demographic history, but it can be reproduced in models of interference among forward and back mutations at large numbers of sites on a nonrecombining chromosome. A coalescent simulation algorithm was used to study this model over a wide range of parameter values. The genealogical distortions are all maximized when the selection coefficients are of critical intermediate sizes, such that Muller's ratchet begins to turn. In this regime, linked neutral nucleotide diversity becomes nearly insensitive to N. Mutations of this size dominate the dynamics even if there are also large numbers of more strongly and more weakly selected sites in the genome. A genealogical perspective on Hill–Robertson interference leads directly to a generalized background-selection model in which the effective population size is progressively reduced going back in time from the present.

Experimentally Increased Codon Bias in the Drosophila Adh Gene Leads to an Increase in Larval, But Not Adult, Alcohol Dehydrogenase Activity [Population and evolutionary genetics]

Hense, W., Anderson, N., Hutter, S., Stephan, W., Parsch, J., Carlini, D. B. — Tue, 26 Jan 2010 14:34:09 PST

Although most amino acids can be encoded by more than one codon, the synonymous codons are not used with equal frequency. This phenomenon is known as codon bias and appears to be a universal feature of genomes. The translational selection hypothesis posits that the use of optimal codons, which match the most abundant species of isoaccepting tRNAs, results in increased translational efficiency and accuracy. Previous work demonstrated that the experimental reduction of codon bias in the Drosophila alcohol dehydrogenase (Adh) gene led to a significant decrease in ADH protein expression. In this study we performed the converse experiment: we replaced seven suboptimal leucine codons that occur naturally in the Drosophila melanogaster Adh gene with the optimal codon. We then compared the in vivo ADH activities imparted by the wild-type and mutant alleles. The introduction of optimal leucine codons led to an increase in ADH activity in third-instar larvae. In adult flies, however, the introduction of optimal codons led to a decrease in ADH activity. There is no evidence that other selectively constrained features of the Adh gene, or its rate of transcription, were altered by the synonymous replacements. These results are consistent with translational selection for codon bias being stronger in the larval stage and suggest that there may be a selective conflict over optimal codon usage between different developmental stages.

The Genetic Signature of Conditional Expression [Population and evolutionary genetics]

Van Dyken, J. D., Wade, M. J. — Tue, 26 Jan 2010 14:34:09 PST

Conditionally expressed genes have the property that every individual in a population carries and transmits the gene, but only a fraction, , expresses the gene and exposes it to natural selection. We show that a consequence of this pattern of inheritance and expression is a weakening of the strength of natural selection, allowing deleterious mutations to accumulate within and between species and inhibiting the spread of beneficial mutations. We extend previous theory to show that conditional expression in space and time have approximately equivalent effects on relaxing the strength of selection and that the effect holds in a spatially heterogeneous environment even with low migration rates among patches. We support our analytical approximations with computer simulations and delineate the parameter range under which the approximations fail. We model the effects of conditional expression on sequence polymorphism at mutation–selection–drift equilibrium, allowing for neutral sites, and show that sequence variation within and between species is inflated by conditional expression, with the effect being strongest in populations with large effective size. As decreases, more sites are recruited into neutrality, leading to pseudogenization and increased drift load. Mutation accumulation diminishes the degree of adaptation of conditionally expressed genes to rare environments, and the mutational cost of phenotypic plasticity, which we quantify as the plasticity load, is greater for more rarely expressed genes. Our theory connects gene-level relative polymorphism and divergence with the spatial and temporal frequency of environments inducing gene expression. Our theory suggests that null hypotheses for levels of standing genetic variation and sequence divergence must be corrected to account for the frequency of expression of the genes under study.

Signatures of Recent Directional Selection Under Different Models of Population Expansion During Colonization of New Selective Environments [Population and evolutionary genetics]

Kim, Y., Gulisija, D. — Tue, 26 Jan 2010 14:34:09 PST

A major problem in population genetics is understanding how the genomic pattern of polymorphism is shaped by natural selection and the demographic history of populations. Complex population dynamics confounds patterns of variation and poses serious challenges for identifying genomic imprints of selection. We examine patterns of polymorphism using computer simulations and provide analytical predictions for hitchhiking effects under two models of adaptive niche expansion. The population split (PS) model assumes the separation of a founding population followed by directional selection in the new environment. Here, the new population undergoes a bottleneck and later expands in size. This model has been used in previous studies to account for demographic effects when testing for signatures of selection under colonization or domestication. The genotype-dependent colonization and introgression (GDCI) model is proposed in this study and assumes that a small number of migrants carrying adaptive genotype found a new population, which then grows logistically. The GDCI model also allows for constant migration between the parental and the new population. Both models predict reduction in variation and excess of high frequency of derived alleles relative to neutral expectations, with and without hitchhiking. Under comparable conditions, the GDCI model results in greater reduction in expected heterozygosity and more skew of the site frequency spectrum than the PS model. We also find that soft selective sweeps (fixation of multiple copies of a beneficial mutation) occurs less often in the GDCI model than in the PS model. This result demonstrates the importance of correctly modeling the ecological process in inferring adaptive evolution using DNA sequence polymorphism.

Estimating Breeding Values With Molecular Relatedness and Reconstructed Pedigrees in Natural Mating Populations of Common Sole, Solea Solea [Population and evolutionary genetics]

Blonk, R. J. W., Komen, H., Kamstra, A., van Arendonk, J. A. M. — Fri, 08 Jan 2010 14:15:34 PST

Captive populations where natural mating in groups is used to obtain offspring typically yield unbalanced population structures with highly skewed parental contributions and unknown pedigrees. Consequently, for genetic parameter estimation, relationships need to be reconstructed or estimated using DNA marker data. With missing parents and natural mating groups, commonly used pedigree reconstruction methods are not accurate and lead to loss of data. Relatedness estimators, however, infer relationships between all animals sampled. In this study, we compared a pedigree relatedness method and a relatedness estimator ("molecular relatedness") method using accuracy of estimated breeding values. A commercial data set of common sole, Solea solea, with 51 parents and 1953 offspring ("full data set") was used. Due to missing parents, for 1338 offspring, a pedigree could be reconstructed with 10 microsatellite markers ("reduced data set"). Cross-validation of both methods using the reduced data set showed an accuracy of estimated breeding values of 0.54 with pedigree reconstruction and 0.55 with molecular relatedness. Accuracy of estimated breeding values increased to 0.60 when applying molecular relatedness to the full data set. Our results indicate that pedigree reconstruction and molecular relatedness predict breeding values equally well in a population with skewed contributions to families. This is probably due to the presence of few large full-sib families. However, unlike methods with pedigree reconstruction, molecular relatedness methods ensure availability of all genotyped selection candidates, which results in higher accuracy of breeding value estimation.

Evolution at a High Imposed Mutation Rate: Adaptation Obscures the Load in Phage T7 [Population and evolutionary genetics]

Springman, R., Keller, T., Molineux, I. J., Bull, J. J. — Fri, 08 Jan 2010 14:15:34 PST

Evolution at high mutation rates is expected to reduce population fitness deterministically by the accumulation of deleterious mutations. A high enough rate should even cause extinction (lethal mutagenesis), a principle motivating the clinical use of mutagenic drugs to treat viral infections. The impact of a high mutation rate on long-term viral fitness was tested here. A large population of the DNA bacteriophage T7 was grown with a mutagen, producing a genomic rate of 4 nonlethal mutations per generation, two to three orders of magnitude above the baseline rate. Fitness—viral growth rate in the mutagenic environment—was predicted to decline substantially; after 200 generations, fitness had increased, rejecting the model. A high mutation load was nonetheless evident from (i) many low- to moderate-frequency mutations in the population (averaging 245 per genome) and (ii) an 80% drop in average burst size. Twenty-eight mutations reached high frequency and were thus presumably adaptive, clustered mostly in DNA metabolism genes, chiefly DNA polymerase. Yet blocking DNA polymerase evolution failed to yield a fitness decrease after 100 generations. Although mutagenic drugs have caused viral extinction in vitro under some conditions, this study is the first to match theory and fitness evolution at a high mutation rate. Failure of the theory challenges the quantitative basis of lethal mutagenesis and highlights the potential for adaptive evolution at high mutation rates.

Human Triallelic Sites: Evidence for a New Mutational Mechanism? [Population and evolutionary genetics]

Hodgkinson, A., Eyre-Walker, A. — Fri, 08 Jan 2010 14:15:34 PST

Most SNPs in the human genome are biallelic; however, there are some sites that are triallelic. We show here that there are approximately twice as many triallelic sites as we would expect by chance. This excess does not appear to be caused by natural selection or mutational hotspots. Instead we propose that a new mutation can induce another mutation either within the same individual or subsequently during recombination. We provide evidence for this model by showing that the rarer two alleles at triallelic sites tend to cluster on phylogenetic trees of human haplotypes. However, we find no association between the density of triallelic sites and the rate of recombination, which leads us to suggest that triallelic sites might be generated by the simultaneous production of two new mutations within the same individual on the same genetic background. Under this model we estimate that simultaneous mutation contributes ~3% of all distinct SNPs. We also show that there is a twofold excess of adjacent SNPs. Approximately half of these seem to be generated simultaneously since they have identical minor allele frequencies. We estimate that the mutation of adjacent nucleotides accounts for a little less than 1% of all SNPs.

Bayesian Computation and Model Selection Without Likelihoods [Population and evolutionary genetics]

Leuenberger, C., Wegmann, D. — Fri, 08 Jan 2010 14:15:35 PST

Until recently, the use of Bayesian inference was limited to a few cases because for many realistic probability models the likelihood function cannot be calculated analytically. The situation changed with the advent of likelihood-free inference algorithms, often subsumed under the term approximate Bayesian computation (ABC). A key innovation was the use of a postsampling regression adjustment, allowing larger tolerance values and as such shifting computation time to realistic orders of magnitude. Here we propose a reformulation of the regression adjustment in terms of a general linear model (GLM). This allows the integration into the sound theoretical framework of Bayesian statistics and the use of its methods, including model selection via Bayes factors. We then apply the proposed methodology to the question of population subdivision among western chimpanzees, Pan troglodytes verus.

Recurrent Selection on the Winters sex-ratio Genes in Drosophila simulans [Population and evolutionary genetics]

Kingan, S. B., Garrigan, D., Hartl, D. L. — Fri, 08 Jan 2010 14:15:35 PST

Selfish genes, such as meiotic drive elements, propagate themselves through a population without increasing the fitness of host organisms. X-linked (or Y-linked) meiotic drive elements reduce the transmission of the Y (X) chromosome and skew progeny and population sex ratios, leading to intense conflict among genomic compartments. Drosophila simulans is unusual in having a least three distinct systems of X chromosome meiotic drive. Here, we characterize naturally occurring genetic variation at the Winters sex-ratio driver (Distorter on the X or Dox), its progenitor gene (Mother of Dox or MDox), and its suppressor gene (Not Much Yang or Nmy), which have been previously mapped and characterized. We survey three North American populations as well as 13 globally distributed strains and present molecular polymorphism data at the three loci. We find that all three genes show signatures of selection in North America, judging from levels of polymorphism and skews in the site-frequency spectrum. These signatures likely result from the biased transmission of the driver and selection on the suppressor for the maintenance of equal sex ratios. Coalescent modeling indicates that the timing of selection is more recent than the age of the alleles, suggesting that the driver and suppressor are coevolving under an evolutionary "arms race." None of the Winters sex-ratio genes are fixed in D. simulans, and at all loci we find ancestral alleles, which lack the gene insertions and exhibit high levels of nucleotide polymorphism compared to the derived alleles. In addition, we find several "null" alleles that have mutations on the derived Dox background, which result in loss of drive function. We discuss the possible causes of the maintenance of presence–absence polymorphism in the Winters sex-ratio genes.

Genetic Testing of the Hypothesis That Hybrid Male Lethality Results From a Failure in Dosage Compensation [Population and evolutionary genetics]

Barbash, D. A. — Fri, 08 Jan 2010 14:15:35 PST

Several recent studies have suggested that F₁ hybrid male lethality in crosses between Drosophila melanogaster and D. simulans is due to a failure in dosage compensation, caused by incompatibilities between D. simulans dosage compensation proteins and the D. melanogaster X chromosome. Contrary to the predictions of this hypothesis, mutations in four essential D. melanogaster dosage compensation genes are shown here to moderately increase rather than decrease hybrid male viability.

Phylodynamics of Infectious Disease Epidemics [Population and evolutionary genetics]

Tue, 08 Dec 2009 12:05:15 PST

We present a formalism for unifying the inference of population size from genetic sequences and mathematical models of infectious disease in populations. Virus phylogenies have been used in many recent studies to infer properties of epidemics. These approaches rely on coalescent models that may not be appropriate for infectious diseases. We account for phylogenetic patterns of viruses in susceptible–infected (SI), susceptible–infected–susceptible (SIS), and susceptible–infected–recovered (SIR) models of infectious disease, and our approach may be a viable alternative to demographic models used to reconstruct epidemic dynamics. The method allows epidemiological parameters, such as the reproductive number, to be estimated directly from viral sequence data. We also describe patterns of phylogenetic clustering that are often construed as arising from a short chain of transmissions. Our model reproduces the moments of the distribution of phylogenetic cluster sizes and may therefore serve as a null hypothesis for cluster sizes under simple epidemiological models. We examine a small cross-sectional sample of human immunodeficiency (HIV)-1 sequences collected in the United States and compare our results to standard estimates of effective population size. Estimated prevalence is consistent with estimates of effective population size and the known history of the HIV epidemic. While our model accurately estimates prevalence during exponential growth, we find that periods of decline are harder to identify.

Exact Tests for Hardy-Weinberg Proportions [Population and evolutionary genetics]

Engels, W. R. — Tue, 08 Dec 2009 12:05:15 PST

Exact conditional tests are often required to evaluate statistically whether a sample of diploids comes from a population with Hardy–Weinberg proportions or to confirm the accuracy of genotype assignments. This requirement is especially common when the sample includes multiple alleles and sparse data, thus rendering asymptotic methods, such as the common ²-test, unreliable. Such an exact test can be performed using the likelihood ratio as its test statistic rather than the more commonly used probability test. Conceptual advantages in using the likelihood ratio are discussed. A substantially improved algorithm is described to permit the performance of a full-enumeration exact test on sample sizes that are too large for previous methods. An improved Monte Carlo algorithm is also proposed for samples that preclude full enumeration. These algorithms are about two orders of magnitude faster than those currently in use. Finally, methods are derived to compute the number of possible samples with a given set of allele counts, a useful quantity for evaluating the feasibility of the full enumeration procedure. Software implementing these methods, ExactoHW, is provided.

Identification of the Major Sex-Determining Region of Turbot (Scophthalmus maximus) [Population and evolutionary genetics]

Tue, 08 Dec 2009 12:05:15 PST

Sex determination in fish is a labile character in evolutionary terms. The sex-determining (SD) master gene can differ even between closely related fish species. This group is an interesting model for studying the evolution of the SD region and the gonadal differentiation pathway. The turbot (Scophthalmus maximus) is a flatfish of great commercial value, where a strong sexual dimorphism exists for growth rate. Following a QTL and marker association approach in five families and a natural population, we identified the main SD region of turbot at the proximal end of linkage group (LG) 5, close to the SmaUSC-E30 marker. The refined map of this region suggested that this marker would be 2.6 cM and 1.4 Mb from the putative SD gene. This region appeared mostly undifferentiated between males and females, and no relevant recombination frequency differences were detected between sexes. Comparative genomics of LG5 marker sequences against five model species showed no similarity of this chromosome to the sex chromosomes of medaka, stickleback, and fugu, but suggested a similarity to a sex-associated QTL from Oreochromis spp. The segregation analysis of the closest markers to the SD region demonstrated a ZW/ZZ model of sex determination in turbot. A small proportion of families did not fit perfectly with this model, which suggests that other minor genetic and/or environmental factors are involved in sex determination in this species.

The Genetic Basis of Phenotypic Adaptation II: The Distribution of Adaptive Substitutions in the Moving Optimum Model [Population and evolutionary genetics]

Kopp, M., Hermisson, J. — Tue, 08 Dec 2009 12:05:15 PST

We consider a population that adapts to a gradually changing environment. Our aim is to describe how ecological and genetic factors combine to determine the genetic basis of adaptation. Specifically, we consider the evolution of a polygenic trait that is under stabilizing selection with a moving optimum. The ecological dynamics are defined by the strength of selection, tilde;, and the speed of the optimum, v; the key genetic parameters are the mutation rate and the variance of the effects of new mutations, . We develop analytical approximations within an "adaptive-walk" framework and describe how selection acts as a sieve that transforms a given distribution of new mutations into the distribution of adaptive substitutions. Our analytical results are complemented by individual-based simulations. We find that (i) the ecological dynamics have a strong effect on the distribution of adaptive substitutions and their impact depends largely on a single composite measure =v/(tilde;³), which combines the ecological and genetic parameters; (ii) depending on , we can distinguish two distinct adaptive regimes: for large the adaptive process is mutation limited and dominated by genetic constraints, whereas for small it is environmentally limited and dominated by the external ecological dynamics; (iii) deviations from the adaptive-walk approximation occur for large mutation rates, when different mutant alleles interact via linkage or epistasis; and (iv) in contrast to predictions from previous models assuming constant selection, the distribution of adaptive substitutions is generally not exponential.

X-Linked Variation in Immune Response in Drosophila melanogaster [Population and evolutionary genetics]

Hill-Burns, E. M., Clark, A. G. — Tue, 08 Dec 2009 12:05:15 PST

This study quantifies the effects of naturally occurring X-linked variation on immune response in Drosophila melanogaster to assess associations between immunity genotypes and innate immune response. We constructed a set of 168 X-chromosomal extraction lines, incorporating X chromosomes from a natural population into co-isogenic autosomal backgrounds, and genotyped the lines at 88 SNPs in 20 X-linked immune genes. We find that genetic variation in many of the genes is associated with immune response phenotypes, including bacterial load and immune gene expression. Many of the associations act in a sex-specific or sexually antagonistic manner, supporting the theory that with the selective pressures facing genes on the X chromosome, sexually antagonistic variation may be more easily maintained.

Measuring and Detecting Molecular Adaptation in Codon Usage Against Nonsense Errors During Protein Translation [Population and evolutionary genetics]

Gilchrist, M. A., Shah, P., Zaretzki, R. — Tue, 08 Dec 2009 12:05:15 PST

Codon usage bias (CUB) has been documented across a wide range of taxa and is the subject of numerous studies. While most explanations of CUB invoke some type of natural selection, most measures of CUB adaptation are heuristically defined. In contrast, we present a novel and mechanistic method for defining and contextualizing CUB adaptation to reduce the cost of nonsense errors during protein translation. Using a model of protein translation, we develop a general approach for measuring the protein production cost in the face of nonsense errors of a given allele as well as the mean and variance of these costs across its coding synonyms. We then use these results to define the nonsense error adaptation index (NAI) of the allele or a contiguous subset thereof. Conceptually, the NAI value of an allele is a relative measure of its elevation on a specific and well-defined adaptive landscape. To illustrate its utility, we calculate NAI values for the entire coding sequence and across a set of nonoverlapping windows for each gene in the Saccharomyces cerevisiae S288c genome. Our results provide clear evidence of adaptation to reduce the cost of nonsense errors and increasing adaptation with codon position and expression. The magnitude and nature of this adaptation are also largely consistent with simulation results in which nonsense errors are the only selective force driving CUB evolution. Because NAI is derived from mechanistic models, it is both easier to interpret and more amenable to future refinement than other commonly used measures of codon bias. Further, our approach can also be used as a starting point for developing other mechanistically derived measures of adaptation such as for translational accuracy.

Adaptive Divergence in Experimental Populations of Pseudomonas fluorescens. IV. Genetic Constraints Guide Evolutionary Trajectories in a Parallel Adaptive Radiation [Population and evolutionary genetics]

McDonald, M. J., Gehrig, S. M., Meintjes, P. L., Zhang, X.-X., Rainey, P. B. — Wed, 18 Nov 2009 08:39:31 PST

The capacity for phenotypic evolution is dependent upon complex webs of functional interactions that connect genotype and phenotype. Wrinkly spreader (WS) genotypes arise repeatedly during the course of a model Pseudomonas adaptive radiation. Previous work showed that the evolution of WS variation was explained in part by spontaneous mutations in wspF, a component of the Wsp-signaling module, but also drew attention to the existence of unknown mutational causes. Here, we identify two new mutational pathways (Aws and Mws) that allow realization of the WS phenotype: in common with the Wsp module these pathways contain a di-guanylate cyclase-encoding gene subject to negative regulation. Together, mutations in the Wsp, Aws, and Mws regulatory modules account for the spectrum of WS phenotype-generating mutations found among a collection of 26 spontaneously arising WS genotypes obtained from independent adaptive radiations. Despite a large number of potential mutational pathways, the repeated discovery of mutations in a small number of loci (parallel evolution) prompted the construction of an ancestral genotype devoid of known (Wsp, Aws, and Mws) regulatory modules to see whether the types derived from this genotype could converge upon the WS phenotype via a novel route. Such types—with equivalent fitness effects—did emerge, although they took significantly longer to do so. Together our data provide an explanation for why WS evolution follows a limited number of mutational pathways and show how genetic architecture can bias the molecular variation presented to selection.

Testing for Spatially Divergent Selection: Comparing QST to FST [Population and evolutionary genetics]

Whitlock, M. C., Guillaume, F. — Wed, 18 Nov 2009 08:39:31 PST

Q_ST is a standardized measure of the genetic differentiation of a quantitative trait among populations. The distribution of Q_ST's for neutral traits can be predicted from the F_ST for neutral marker loci. To test for the neutral differentiation of a quantitative trait among populations, it is necessary to ask whether the Q_ST of that trait is in the tail of the probability distribution of neutral traits. This neutral distribution can be estimated using the Lewontin–Krakauer distribution and the F_ST from a relatively small number of marker loci. We develop a simulation method to test whether the Q_ST of a given trait is consistent with the null hypothesis of selective neutrality over space. The method is most powerful with small mean F_ST, strong selection, and a large number (>10) of measured populations. The power and type I error rate of the new method are far superior to the traditional method of comparing Q_ST and F_ST.

Population Differentiation as an Indicator of Recent Positive Selection in Humans: An Empirical Evaluation [Population and evolutionary genetics]

Wed, 18 Nov 2009 08:39:31 PST

We have evaluated the extent to which SNPs identified by genomewide surveys as showing unusually high levels of population differentiation in humans have experienced recent positive selection, starting from a set of 32 nonsynonymous SNPs in 27 genes highlighted by the HapMap1 project. These SNPs were genotyped again in the HapMap samples and in the Human Genome Diversity Project–Centre d'Etude du Polymorphisme Humain (HGDP–CEPH) panel of 52 populations representing worldwide diversity; extended haplotype homozygosity was investigated around all of them, and full resequence data were examined for 9 genes (5 from public sources and 4 from new data sets). For 7 of the genes, genotyping errors were responsible for an artifactual signal of high population differentiation and for 2, the population differentiation did not exceed our significance threshold. For the 18 genes with confirmed high population differentiation, 3 showed evidence of positive selection as measured by unusually extended haplotypes within a population, and 7 more did in between-population analyses. The 9 genes with resequence data included 7 with high population differentiation, and 5 showed evidence of positive selection on the haplotype carrying the nonsynonymous SNP from skewed allele frequency spectra; in addition, 2 showed evidence of positive selection on unrelated haplotypes. Thus, in humans, high population differentiation is (apart from technical artifacts) an effective way of enriching for recently selected genes, but is not an infallible pointer to recent positive selection supported by other lines of evidence.

The Population Genetics of Adaptation: Multiple Substitutions on a Smooth Fitness Landscape [Population and evolutionary genetics]

Unckless, R. L., Orr, H. A. — Wed, 18 Nov 2009 08:39:31 PST

Much recent work in the theoretical study of adaptation has focused on the so-called strong selection–weak mutation (SSWM) limit, wherein adaptation is due to new mutations of definite selective advantage. This work, in turn, has focused on the first step (substitution) during adaptive evolution. Here we extend this theory to allow multiple steps during adaptation. We find analytic solutions to the probability that adaptation follows a certain path during evolution as well as the probability that adaptation arrives at a given genotype regardless of the path taken. We also consider the probability of parallel adaptation and the proportion of the total increase in fitness caused by the first substitution. Our key assumption is that there is no epistasis among beneficial mutations.

Closed-Form Two-Locus Sampling Distributions: Accuracy and Universality [Population and evolutionary genetics]

Jenkins, P. A., Song, Y. S. — Wed, 18 Nov 2009 08:39:31 PST

Sampling distributions play an important role in population genetics analyses, but closed-form sampling formulas are generally intractable to obtain. In the presence of recombination, there is no known closed-form sampling formula that holds for an arbitrary recombination rate. However, we recently showed that it is possible to obtain useful closed-form sampling formulas when the population-scaled recombination rate is large. Specifically, in the case of the two-locus infinite-alleles model, we considered an asymptotic expansion of the sampling formula in inverse powers of and obtained closed-form expressions for the first few terms in the expansion. In this article, we generalize this result to an arbitrary finite-alleles mutation model and show that, up to the first few terms in the expansion that we are able to compute analytically, the functional form of the asymptotic sampling formula is common to all mutation models. We carry out an extensive study of the accuracy of the asymptotic formula for the two-locus parent-independent mutation model and discuss in detail a concrete application in the context of the composite-likelihood method. Furthermore, using our asymptotic sampling formula, we establish a simple sufficient condition for a given two-locus sample configuration to have a finite maximum-likelihood estimate (MLE) of . This condition is the first analytic result on the classification of the MLE of and is instantaneous to check in practice, provided that one-locus probabilities are known.

The Sheltered Genetic Load Linked to the S Locus in Plants: New Insights From Theoretical and Empirical Approaches in Sporophytic Self-Incompatibility [Population and evolutionary genetics]

Llaurens, V., Gonthier, L., Billiard, S. — Wed, 18 Nov 2009 08:39:31 PST

Inbreeding depression and mating systems evolution are closely linked, because the purging of deleterious mutations and the fitness of individuals may depend on outcrossing vs. selfing rates. Further, the accumulation of deleterious mutations may vary among genomic regions, especially for genes closely linked to loci under balancing selection. Sporophytic self-incompatibility (SSI) is a common genetic mechanism in angiosperm that enables hermaphrodite plants to avoid selfing and promote outcrossing. The SSI phenotype is determined by the S locus and may depend on dominance relationships among alleles. Since most individuals are heterozygous at the S locus and recombination is suppressed in the S-locus region, it has been suggested that deleterious mutations could accumulate at genes linked to the S locus, generating a "sheltered load." In this article, we first theoretically investigate the conditions generating sheltered load in SSI. We show that deleterious mutations can accumulate in linkage with specific S alleles, and particularly if those S alleles are dominant. Second, we looked for the presence of sheltered load in Arabidopsis halleri using CO₂ gas treatment to overcome self-incompatibility. By examining the segregation of S alleles and measuring the relative fitness of progeny, we found significant sheltered load associated with the most dominant S allele (S15) of three S alleles tested. This sheltered load seems to be expressed at several stages of the life cycle and to have a larger effect than genomic inbreeding depression.