Lack of polymorphism on the Drosophila fourth chromosome resulting from selection

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Landmarks in the history of selective sweeps

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modelled by a realistic mutation rate and as part of a realistic distribution of fitness effects (DFE), as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modelled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false positive rates are in excess of true positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Teaser Text

Outlier-based genomic scans have proven a popular approach for identifying loci that have potentially experienced recent positive selection. However, it has previously been shown that an evolutionarily appropriate baseline model that incorporates non-equilibrium population histories, purifying and background selection, and variation in mutation and recombination rates is necessary to reduce often extreme false positive rates when performing genomic scans. Here we evaluate the power to detect recurrent selective sweeps using common SFS-based and haplotype-based methods under these increasingly realistic models. We find that while these appropriate evolutionary baselines are essential to reduce false positive rates, the power to accurately detect recurrent selective sweeps is generally low across much of the biologically relevant parameter space.

Effect of beneficial sweeps and background selection on genetic diversity in changing environments

Neutral theory predicts that the genetic diversity within a population is proportional to the effective population size. In contrast, the observed genetic diversity for various species vary in a narrow range for several orders of magnitude change in the population sizes (Lewontin's paradox). The selective sweeps and background selection, reduce the genetic variation at the linked neutral sites and have been studied considering the environment to be selectively constant. However, in a natural population, the selective environment varies with time. Here, we investigate the impact of selective sweeps and background selection on neutral genetic diversity when the selection coefficient changes periodically over time. The reduction in genetic variation due to selective sweeps is known to depend on the conditional fixation time. Here, we find that the effect of changing environment on conditional mean fixation time is most substantial for the randomly mating population than the inbreeding population with arbitrary inbreeding coefficient. We also study the effect of background selection on neutral sites when the selection coefficient of linked deleterious mutation change periodically in time. In the slowly changing environment, we find that neutral heterozygosity is significantly different, and the site frequency spectrum has a different shape than that in the static environment.

Genomic signatures of recent adaptation in a wild bumblebee

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.

Effects of Selection at Linked Sites on Patterns of Genetic Variability

Patterns of variation and evolution at a given site in a genome can be strongly influenced by the effects of selection at genetically linked sites. In particular, the recombination rates of genomic regions correlate with their amount of within-population genetic variability, the degree to which the frequency distributions of DNA sequence variants differ from their neutral expectations, and the levels of adaptation of their functional components. We review the major population genetic processes that are thought to lead to these patterns, focusing on their effects on patterns of variability: selective sweeps, background selection, associative overdominance, and Hill–Robertson interference among deleterious mutations. We emphasize the difficulties in distinguishing among the footprints of these processes and disentangling them from the effects of purely demographic factors such as population size changes. We also discuss how interactions between selective and demographic processes can significantly affect patterns of variability within genomes.

The classical hitchhiking model with continuous mutational pressure and purifying selection

Detecting selective sweeps driven by strong positive selection and localizing the targets of selection in the genome play a major role in modern population genetics and genomics. Most of these analyses are based on the classical model of genetic hitchhiking proposed by Maynard Smith and Haigh (1974, Genetical Research, 23, 23). Here, we consider extensions of the classical two-locus model. Introducing mutation at the strongly selected site, we analyze the conditions under which soft sweeps may arise. We identify a new parameter (the ratio of the beneficial mutation rate to the selection coefficient) that characterizes the occurrence of multiple-origin soft sweeps. Furthermore, we quantify the hitchhiking effect when the polymorphism at the linked locus is not neutral but maintained in a mutation-selection balance. In this case, we find a smaller relative reduction of heterozygosity at the linked site than for a neutral polymorphism. In our analysis, we use a semi-deterministic approach; i.e., we analyze the frequency process of the beneficial allele in an infinitely large population when its frequency is above a certain threshold; however, for very small frequencies in the initial phase after the onset of selection we rely on diffusion theory.

Natural variation in fecundity is correlated with species-wide levels of divergence in Caenorhabditis elegans

Life history traits underlie the fitness of organisms and are under strong natural selection. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g., a selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of a species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fecundity of 121 wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions have significantly higher fecundity than strains without evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fecundity variation. In addition, we mapped previous fecundity data from wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines that were grown in various conditions and detected eight QTL using GWA and linkage mappings. These QTL show the genetic complexity of fecundity across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fecundity than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that caused increased fecundity could underlie the selective sweeps during the worldwide expansion of C. elegans.

Population and comparative genetics of thermotolerance divergence between yeast species

Many familiar traits in the natural world-from lions' manes to the longevity of bristlecone pine trees-arose in the distant past, and have long since fixed in their respective species. A key challenge in evolutionary genetics is to figure out how and why species-defining traits have come to be. We used the thermotolerance growth advantage of the yeast Saccharomyces cerevisiae over its sister species Saccharomyces paradoxus as a model for addressing these questions. Analyzing loci at which the S. cerevisiae allele promotes thermotolerance, we detected robust evidence for positive selection, including amino acid divergence between the species and conservation within S. cerevisiae populations. Because such signatures were particularly strong at the chromosome segregation gene ESP1, we used this locus as a case study for focused mechanistic follow-up. Experiments revealed that, in culture at high temperature, the S. paradoxus ESP1 allele conferred a qualitative defect in biomass accumulation and cell division relative to the S. cerevisiae allele. Only genetic divergence in the ESP1 coding region mattered phenotypically, with no functional impact detectable from the promoter. Our data support a model in which an ancient ancestor of S. cerevisiae, under selection to boost viability at high temperature, acquired amino acid variants at ESP1 and many other loci, which have been constrained since then. Complex adaptations of this type hold promise as a paradigm for interspecies genetics, especially in deeply diverged traits that may have taken millions of years to evolve.

How Good Are Predictions of the Effects of Selective Sweeps on Levels of Neutral Diversity?

Selective sweeps are thought to play a significant role in shaping patterns of variability across genomes; accurate predictions of their effects are, therefore, important for understanding these patterns. A commonly used model of selective sweeps assumes that alleles sampled at the end of a sweep, and that fail to recombine with wild-type haplotypes during the sweep, coalesce instantaneously, leading to a simple expression for sweep effects on diversity. It is shown here that there can be a significant probability that a pair of alleles sampled at the end of a sweep coalesce during the sweep before a recombination event can occur, reducing their expected coalescent time below that given by the simple approximation. Expressions are derived for the expected reductions in pairwise neutral diversities caused by both single and recurrent sweeps in the presence of such within-sweep coalescence, although the effects of multiple recombination events during a sweep are only treated heuristically. The accuracies of the resulting expressions were checked against the results of simulations. For even moderate ratios of the recombination rate to the selection coefficient, the simple approximation can be substantially inaccurate. The selection model used here can be applied to favorable mutations with arbitrary dominance coefficients, to sex-linked loci with sex-specific selection coefficients, and to inbreeding populations. Using the results from this model, the expected differences between the levels of variability on X chromosomes and autosomes with selection at linked sites are discussed, and compared with data on a population of Drosophila melanogaster.

Sex Chromosome Evolution in Muscid Flies

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. Population genetics theory has been developed and tested to explain the rapid evolution of sex chromosomes and sex determination. However, we do not know why the sex chromosomes are divergent in some taxa and conserved in others. Addressing this question requires comparing closely related taxa with conserved and divergent sex chromosomes to identify biological features that could explain these differences. Cytological karyotypes suggest that muscid flies (e.g., house fly) and blow flies are such a taxonomic pair. The sex chromosomes appear to differ across muscid species, whereas they are conserved across blow flies. Despite the cytological evidence, we do not know the extent to which muscid sex chromosomes are independently derived along different evolutionary lineages. To address that question, we used genomic and transcriptomic sequence data to identify young sex chromosomes in two closely related muscid species, horn fly (Haematobia irritans) and stable fly (Stomoxys calcitrans). We provide evidence that the nascent sex chromosomes of horn fly and stable fly were derived independently from each other and from the young sex chromosomes of the closely related house fly (Musca domestica). We present three different scenarios that could have given rise to the sex chromosomes of horn fly and stable fly, and we describe how the scenarios could be distinguished. Distinguishing between these scenarios in future work could identify features of muscid genomes that promote sex chromosome divergence.

Drug Resistance Evolution in HIV in the Late 1990s: Hard Sweeps, Soft Sweeps, Clonal Interference and the Accumulation of Drug Resistance Mutations

The evolution of drug resistance in pathogens such as HIV is an important and widely known example in the field of evolutionary medicine. Here, we focus on a unique data set from the late 1990s with multiple viral sequences from multiple time points in 118 patients. We study patterns of evolutionary dynamics in the viral populations in these patients who were treated with Reverse Transcriptase Inhibitors and Protease Inhibitors in the late 1990s. Specifically, we aim to visualize and analyze examples of population genetic processes such as selective sweeps and clonal interference. The figures and descriptions in this paper can be used in evolution and population genetics classes. We show and analyze a wide variety of patterns, specifically: soft sweeps, hard sweeps, softening sweeps and hardening sweeps, simultaneous sweeps, accumulation of mutations and clonal interference.

The Effects on Neutral Variability of Recurrent Selective Sweeps and Background Selection

Levels of variability and rates of adaptive evolution may be affected by hitchhiking, the effect of selection on evolution at linked sites. Hitchhiking can be caused either by "selective sweeps" or by background selection, involving the spread of new favorable alleles or the elimination of deleterious mutations, respectively. Recent analyses of population genomic data have fitted models where both these processes act simultaneously, to infer the parameters of selection. Here, we investigate the consequences of relaxing a key assumption of some of these studies, that the time occupied by a selective sweep is negligible compared with the neutral coalescent time. We derive a new expression for the expected level of neutral variability in the presence of recurrent selective sweeps and background selection. We also derive approximate integral expressions for the effects of recurrent selective sweeps. The accuracy of the theoretical predictions was tested against multilocus simulations, with selection, recombination, and mutation parameters that are realistic for Drosophila melanogaster In the presence of crossing over, there is approximate agreement between the theoretical and simulation results. We show that the observed relationships between the rate of crossing over, and the level of synonymous site diversity and rate of adaptive evolution in Drosophila are probably mainly caused by background selection, whereas selective sweeps and population size changes are needed to produce the observed distortions of the site frequency spectrum.

Selective Sweeps

For almost 20 years, many inference methods have been developed to detect selective sweeps and localize the targets of directional selection in the genome. These methods are based on population genetic models that describe the effect of a beneficial allele (e.g., a new mutation) on linked neutral variation (driven by directional selection from a single copy to fixation). Here, I discuss these models, ranging from selective sweeps in a panmictic population of constant size to evolutionary traffic when simultaneous sweeps at multiple loci interfere, and emphasize the important role of demography and population structure in data analysis. In the past 10 years, soft sweeps that may arise after an environmental change from directional selection on standing variation have become a focus of population genetic research. In contrast to selective sweeps, they are caused by beneficial alleles that were neutrally segregating in a population before the environmental change or were present at a mutation-selection balance in appreciable frequency.

The Drosophila Dot Chromosome: Where Genes Flourish Amidst Repeats

The F element of the Drosophila karyotype (the fourth chromosome in Drosophila melanogaster) is often referred to as the "dot chromosome" because of its appearance in a metaphase chromosome spread. This chromosome is distinct from other Drosophila autosomes in possessing both a high level of repetitious sequences (in particular, remnants of transposable elements) and a gene density similar to that found in the other chromosome arms, ∼80 genes distributed throughout its 1.3-Mb "long arm." The dot chromosome is notorious for its lack of recombination and is often neglected as a consequence. This and other features suggest that the F element is packaged as heterochromatin throughout. F element genes have distinct characteristics (e.g, low codon bias, and larger size due both to larger introns and an increased number of exons), but exhibit expression levels comparable to genes found in euchromatin. Mapping experiments show the presence of appropriate chromatin modifications for the formation of DNaseI hypersensitive sites and transcript initiation at the 5' ends of active genes, but, in most cases, high levels of heterochromatin proteins are observed over the body of these genes. These various features raise many interesting questions about the relationships of chromatin structures with gene and chromosome function. The apparent evolution of the F element as an autosome from an ancestral sex chromosome also raises intriguing questions. The findings argue that the F element is a unique chromosome that occupies its own space in the nucleus. Further study of the F element should provide new insights into chromosome structure and function.

Support for the Dominance Theory in Drosophila Transcriptomes

Interactions among divergent elements of transcriptional networks from different species can lead to misexpression in hybrids through regulatory incompatibilities, some with the potential to generate sterility. While the possible contribution of faster-male evolution to this misexpression has been explored, the role of the hemizygous X chromosome (i.e., the dominance theory for transcriptomes) remains yet to be determined. Here, we study genome-wide patterns of gene expression in females and males of Drosophila yakuba, Drosophila santomea and their hybrids. We used attached-X stocks to specifically test the dominance theory, and we uncovered a significant contribution of recessive alleles on the X chromosome to hybrid misexpression. Our analyses also suggest a contribution of weakly deleterious regulatory mutations to gene expression divergence in genes with sex-biased expression, but only in the sex toward which the expression is biased (e.g., genes with female-biased expression when analyzed in females). In the opposite sex, we found stronger selective constraints on gene expression divergence. Although genes with a high degree of male-biased expression show a clear signal of faster-X evolution of gene expression, we also detected slower-X evolution in other gene classes (e.g., female-biased genes). This slower-X effect is mediated by significant decreases in cis- and trans-regulatory divergence. The distinct behavior of X-linked genes with a high degree of male-biased expression is consistent with these genes experiencing a higher incidence of positively selected regulatory mutations than their autosomal counterparts.

USING HITCHHIKING GENES TO STUDY ADAPTATION AND DIVERGENCE DURING SPECIATION WITHIN THE DROSOPHILA MELANOGASTER SPECIES COMPLEX

Several studies of intraspecific and interspecific DNA sequence variation from Drosophila loci have revealed a pattern of low intraspecific variation from genomic regions of low recombination. The mechanisms consistently invoked to explain these patterns are the selective sweep of advantageous mutations together with genetic hitchhiking of linked loci. To examine the effect of selective sweeps on genetic divergence during speciation, we studied two loci in different genomic regions thought to be subject to selective sweeps. We obtained DNA sequences from 1.1kb pair portions of the fourth chromosome locus cubitus interruptus Dominant (ci^D ) and from the asense locus near the telomere of the X chromosome. At ci^D , we found very low variation among multiple lines of Drosophila mauritiana and D. sechellia. This finding is consistent with an earlier report of very low variation in D. melanogaster and D. simulans at ci^D and supports the conclusion of selective sweeps and genetic hitchhiking on the nonrecombining fourth chromosome. The pattern of variation found at asense suggests that a selective sweep has occurred recently at the tip of the X chromosome in D. simulans, but not in D. melanogaster or D. mauritiana. The data from ci^D and asense are compared with data from three X chromosome loci (period, zeste, and yolk protein 2) that experience normal levels of recombination. By examining estimated genealogies and the rates at which different classes of mutations have accumulated, we conclude that selective sweeps are common occurrences on the fourth chromosome but less common near the tip of the X chromosome. An interesting pattern of low variation at ci^D among D. simulans, D. mauritiana, and D. sechellia suggests that a selective sweep may have occurred among these forms even after divergence into separate species had begun.

Genetic pattern and gene localization of polydactyly in Beijing fatty chicken

Polydactyly, a common heritable limb malformation in vertebrates, is characterized by supernumerary digits. In chickens, basic characteristics and rough dominant genes have been explored in past decades; however, the elaborate pattern of inheritance and the determinant gene remain obscure. In this study, different types of polydactylism were classified by the numbers and the shapes of toes, including the newly defined subtypes of B’ and G, for the Beijing fatty chicken, a native breed of chicken from China. Through experiments on hybridization, we demonstrated a complete dominant inheritance of polydactyly instead of an incomplete penetrance or genetic modification of the previous conjecture. In particular, by using the F2 population of the five-digit purebred line of Beijing fatty chicken backcrossed to Shiqiza chicken and by using restriction-site associated DNA based markers, we performed a genome-wide association study on the trait of polydactyly. Furthermore, whole genome resequencing strategy was applied to sweep SNPs across the whole genome. An outlier-based F_st approach was employed to search for signatures of selection, and results indicated that the determinant mutation was found in the region ranging from 8.3 Mb to 8.7 Mb, where the polydactyly candidate gene LMBR1 was located. The G/T mutation of rs80659072 was identified to be highly associated with polydactyly in our resequencing and was validated in random samples from an expanded population. Thus, we confirmed that LMBR1 was the causative gene of polydactyly in the Beijing fatty chicken by using GWAS with restriction-site associated DNA based markers and resequencing.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

What can genome-wide association studies tell us about the evolutionary forces maintaining genetic variation for quantitative traits?

Contents 21 I. 21 II. 22 III. 24 IV. 25 V. 29 30 References 30 SUMMARY: Understanding the evolutionary forces that shape genetic variation within species has long been a goal of evolutionary biology. Integrating data for the genetic architecture of traits from genome-wide association mapping studies (GWAS) along with the development of new population genetic methods for identifying selection in sequence data may allow us to evaluate the roles of mutation-selection balance and balancing selection in shaping genetic variation at various scales. Here, we review the theoretical predictions for genetic architecture and additional signals of selection on genomic sequence for the loci that affect traits. Next, we review how plant GWAS have tested for the signatures of various selective scenarios. Limited evidence to date suggests that within-population variation is maintained primarily by mutation-selection balance while variation across the landscape is the result of local adaptation. However, there are a number of inherent biases in these interpretations. We highlight these challenges and suggest ways forward to further understanding of the maintenance of variation.

Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster

Background

Evolutionary theory predicts that antagonistically selected alleles, such as those with divergent pleiotropic effects in early and late life, may often reach intermediate population frequencies due to balancing selection, an elusive process when sought out empirically. Alternatively, genetic diversity may increase as a result of positive frequency-dependent selection and genetic purging in bottlenecked populations.

Results

While experimental evolution systems with directional phenotypic selection typically result in at least local heterozygosity loss, we report that selection for increased lifespan in Drosophila melanogaster leads to an extensive genome-wide increase of nucleotide diversity in the selected lines compared to replicate control lines, pronounced in regions with no or low recombination, such as chromosome 4 and centromere neighborhoods. These changes, particularly in coding sequences, are most consistent with the operation of balancing selection and the antagonistic pleiotropy theory of aging and life history traits that tend to be intercorrelated. Genes involved in antioxidant defenses, along with multiple lncRNAs, were among those most affected by balancing selection. Despite the overwhelming genetic diversification and the paucity of selective sweep regions, two genes with functions important for central nervous system and memory, Ptp10D and Ank2, evolved under positive selection in the longevity lines.

Conclusions

Overall, the ‘evolve-and-resequence’ experimental approach proves successful in providing unique insights into the complex evolutionary dynamics of genomic regions responsible for longevity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3485-0) contains supplementary material, which is available to authorized users.

Estimating the Effective Population Size from Temporal Allele Frequency Changes in Experimental Evolution

The effective population size (Ne) is a major factor determining allele frequency changes in natural and experimental populations. Temporal methods provide a powerful and simple approach to estimate short-term Ne. They use allele frequency shifts between temporal samples to calculate the standardized variance, which is directly related to Ne. Here we focus on experimental evolution studies that often rely on repeated sequencing of samples in pools (Pool-seq). Pool-seq is cost-effective and often outperforms individual-based sequencing in estimating allele frequencies, but it is associated with atypical sampling properties: Additional to sampling individuals, sequencing DNA in pools leads to a second round of sampling, which increases the variance of allele frequency estimates. We propose a new estimator of Ne, which relies on allele frequency changes in temporal data and corrects for the variance in both sampling steps. In simulations, we obtain accurate Ne estimates, as long as the drift variance is not too small compared to the sampling and sequencing variance. In addition to genome-wide Ne estimates, we extend our method using a recursive partitioning approach to estimate Ne locally along the chromosome. Since the type I error is controlled, our method permits the identification of genomic regions that differ significantly in their Ne estimates. We present an application to Pool-seq data from experimental evolution with Drosophila and provide recommendations for whole-genome data. The estimator is computationally efficient and available as an R package at https://github.com/ThomasTaus/Nest.

S/HIC: Robust Identification of Soft and Hard Sweeps Using Machine Learning

Detecting the targets of adaptive natural selection from whole genome sequencing data is a central problem for population genetics. However, to date most methods have shown sub-optimal performance under realistic demographic scenarios. Moreover, over the past decade there has been a renewed interest in determining the importance of selection from standing variation in adaptation of natural populations, yet very few methods for inferring this model of adaptation at the genome scale have been introduced. Here we introduce a new method, S/HIC, which uses supervised machine learning to precisely infer the location of both hard and soft selective sweeps. We show that S/HIC has unrivaled accuracy for detecting sweeps under demographic histories that are relevant to human populations, and distinguishing sweeps from linked as well as neutrally evolving regions. Moreover, we show that S/HIC is uniquely robust among its competitors to model misspecification. Thus, even if the true demographic model of a population differs catastrophically from that specified by the user, S/HIC still retains impressive discriminatory power. Finally, we apply S/HIC to the case of resequencing data from human chromosome 18 in a European population sample, and demonstrate that we can reliably recover selective sweeps that have been identified earlier using less specific and sensitive methods.

Extraordinary conservation of entire chromosomes in insects over long evolutionary periods

Comparison of the genomes of different Drosophila species has shown that six different chromosomes, the so-called ''Muller elements," constitute the building blocks for all Drosophila species. Here, we confirm previous results suggesting that this conservation of the Muller elements extends far beyond Drosophila, to at least tephritid fruit flies, thought to have diverged from drosophilids 60-70 mYr ago. Less than 10 percent of genes differ in chromosome location between the two insect groups. Within chromosomes, however, the order is highly scrambled, as expected from the comparison between Drosophila species. The data also support the notion that the sex chromosomes of tephritid flies originated from an ancestor of the dot chromosome 4 of Drosophila. Overall, therefore, no new chromosome has been created for perhaps a billion generations over the two evolutionary lines. This stability at the chromosome level, which appears to extend to all Diptera including mosquitoes, is in stark contrast to other groups such as mammals, birds, fish and plants, in which chromosome numbers and organization vary enormously among species that have diverged over much fewer generations.

Signatures of positive selection: from selective sweeps at individual loci to subtle allele frequency changes in polygenic adaptation

In the past 15 years, numerous methods have been developed to detect selective sweeps underlying adaptations. These methods are based on relatively simple population genetic models, including one or two loci at which positive directional selection occurs, and one or two marker loci at which the impact of selection on linked neutral variation is quantified. Information about the phenotype under selection is not included in these models (except for fitness). In contrast, in the quantitative genetic models of adaptation, selection acts on one or more phenotypic traits, such that a genotype-phenotype map is required to bridge the gap to population genetics theory. Here I describe the range of population genetic models from selective sweeps in a panmictic population of constant size to evolutionary traffic when simultaneous sweeps at multiple loci interfere, and I also consider the case of polygenic selection characterized by subtle allele frequency shifts at many loci. Furthermore, I present an overview of the statistical tests that have been proposed based on these population genetics models to detect evidence for positive selection in the genome.

Biogeography of Drosophila (Diptera: Drosophilidae) in East and Southeast Asia

The causes of high biological diversity in biodiversity hotspots have long been a major subject of study in conservation biology. To investigate this matter, we conducted a phylogeographic study of five Drosophila (Diptera: Drosophilidae) species from East and Southeast Asia: Drosophila albomicans Duda, D. formosana Duda, D. immigrans Sturtevant, D. melanogaster Meigen, and D. simulans Sturtevant. We collected 185 samples from 28 localities in eight countries. From each collected individual, we sequenced the autosomal extra sex comb gene (esc) and seven mitochondrial genes, including nicotinamide adenine dinucleotide hydrate-reductase dehydrogenase subunit 4 (ND4), ND4L, tRNA-His, tRNA-Pro, tRNA-Thr, partial ND5, and partial ND6. Phylogenetic analyses using maximum- likelihood and Bayesian methods revealed interesting population structure and identified the existence of two distinct D. formosana lineages (Southeast Asian and Taiwanese populations). Genetic differentiation among groups of D. immigrans suggests the possibility of endemic speciation in Taiwan. In contrast, D. melanogaster remained one extensively large population throughout East and Southeast Asia, including nearby islets. A molecular clock was used to estimate divergence times, which were compared with past geographical events to infer evolutionary scenarios. Our findings suggest that interglacial periods may have caused population isolation, thus enhancing population differentiation more strongly for some of the Drosophila species. The population structure of each Drosophila species in East and Southeast Asia has been influenced by past geographic events.

NextGen sequencing reveals short double crossovers contribute disproportionately to genetic diversity in Toxoplasma gondii

Background

Toxoplasma gondii is a widespread protozoan parasite of animals that causes zoonotic disease in humans. Three clonal variants predominate in North America and Europe, while South American strains are genetically diverse, and undergo more frequent recombination. All three northern clonal variants share a monomorphic version of chromosome Ia (ChrIa), which is also found in unrelated, but successful southern lineages. Although this pattern could reflect a selective advantage, it might also arise from non-Mendelian segregation during meiosis. To understand the inheritance of ChrIa, we performed a genetic cross between the northern clonal type 2 ME49 strain and a divergent southern type 10 strain called VAND, which harbors a divergent ChrIa.

Results

NextGen sequencing of haploid F1 progeny was used to generate a genetic map revealing a low level of conventional recombination, with an unexpectedly high frequency of short, double crossovers. Notably, both the monomorphic and divergent versions of ChrIa were isolated with equal frequency. As well, ChrIa showed no evidence of being a sex chromosome, of harboring an inversion, or distorting patterns of segregation. Although VAND was unable to self fertilize in the cat, it underwent successful out-crossing with ME49 and hybrid survival was strongly associated with inheritance of ChrIII from ME49 and ChrIb from VAND.

Conclusions

Our findings suggest that the successful spread of the monomorphic ChrIa in the wild has not been driven by meiotic drive or related processes, but rather is due to a fitness advantage. As well, the high frequency of short double crossovers is expected to greatly increase genetic diversity among progeny from genetic crosses, thereby providing an unexpected and likely important source of diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1168) contains supplementary material, which is available to authorized users.

The relations between recombination rate and patterns of molecular variation and evolution in Drosophila

Genetic recombination affects levels of variability and the efficacy of selection because natural selection acting at one site affects evolutionary processes at linked sites. The variation in local recombination rates across the Drosophila genome provides excellent material for testing hypotheses concerning the evolutionary consequences of recombination. The current state of knowledge from studies of Drosophila genomics and population genetics is reviewed here. Selection at linked sites has influenced the relations between recombination rates and patterns of molecular variation and evolution, such that higher rates of recombination are associated with both higher levels of variability and a greater efficacy of selection. It seems likely that background selection against deleterious mutations is a major factor contributing to these patterns in genome regions in which crossing over is rare or absent, whereas selective sweeps of positively selected mutations probably play an important role in regions with crossing over.

Dynamics of adaptive alleles in divergently selected body weight lines of chickens

By studying genomic changes over time in populations subjected to strong artificial directional selection, we can gain insights to the dynamics of beneficial alleles originating from the founder population or emerging as novel mutations undergoing ongoing selection. The Virginia lines are a chicken resource population generated by long-term bi-directional, single-trait selection for juvenile body weight. We studied genome-wide allele frequency changes from generation 40 to 53 using genome-wide genotypes from directional and relaxed selection lines. Overall, there were small changes in allele frequencies at individual loci over the studied time period; but, on average, the changes were greater in lines with larger phenotypic changes. This is consistent with previous findings that much of the response to selection over the first 40 years of selection was attributable to utilization of standing genetic variation at many loci in the genome, indicating a mostly polygenic architecture for body weight. Over the course of the selection experiment, the largest phenotypic response to selection was observed in the high-weight selected line, and in this line we detected a single locus where the allele frequency changed rapidly during a late stage of the experiment. This locus likely contains a novel, beneficial mutation that appeared between generations 40 and 45 and was driven to fixation in 5 to 10 generations. This result illustrates the dependence of continued long-term selection response on standing genetic variation at many loci as well as strong, novel, beneficial mutations.

Genome-wide analysis in chicken reveals that local levels of genetic diversity are mainly governed by the rate of recombination

Background

Polymorphism is key to the evolutionary potential of populations. Understanding which factors shape levels of genetic diversity within genomes forms a central question in evolutionary genomics and is of importance for the possibility to infer episodes of adaptive evolution from signs of reduced diversity. There is an on-going debate on the relative role of mutation and selection in governing diversity levels. This question is also related to the role of recombination because recombination is expected to indirectly affect polymorphism via the efficacy of selection. Moreover, recombination might itself be mutagenic and thereby assert a direct effect on diversity levels.

Results

We used whole-genome re-sequencing data from domestic chicken (broiler and layer breeds) and its wild ancestor (the red jungle fowl) to study the relationship between genetic diversity and several genomic parameters. We found that recombination rate had the largest effect on local levels of nucleotide diversity. The fact that divergence (a proxy for mutation rate) and recombination rate were negatively correlated argues against a mutagenic role of recombination. Furthermore, divergence had limited influence on polymorphism.

Conclusions

Overall, our results are consistent with a selection model, in which regions within a short distance from loci under selection show reduced polymorphism levels. This conclusion lends further support from the observations of strong correlations between intergenic levels of diversity and diversity at synonymous as well as non-synonymous sites. Our results also demonstrate differences between the two domestic breeds and red jungle fowl, where the domestic breeds show a stronger relationship between intergenic diversity levels and diversity at synonymous and non-synonymous sites. This finding, together with overall lower diversity levels in domesticates compared to red jungle fowl, seem attributable to artificial selection during domestication.

Witnessing Phenotypic and Molecular Evolution in the Fruit Fly

This multi-day exercise is designed for a college Genetics and Evolution laboratory to demonstrate concepts of inheritance and phenotypic and molecular evolution using a live model organism, Drosophila simulans. Students set up an experimental fruit fly population consisting of ten white eyed flies and one red eyed fly. Having red eyes is advantageous compared to having white eyes, allowing students to track the spread of this advantageous trait over several generations. Ultimately, the students perform PCR and gel electrophoresis at two neutral markers, one located in close proximity to the eye-color locus, and one located at the other end of the chromosome. Students observe that most flies have red eyes, and these red-eyed flies have lost variation at the near marker, but maintained variation at the far marker, hence observing a "selective sweep" and the "hitchhiking" of a nearby neutral variant. Students literally observe phenotypic and molecular evolution in their classroom!

The many landscapes of recombination in Drosophila melanogaster

Recombination is a fundamental biological process with profound evolutionary implications. Theory predicts that recombination increases the effectiveness of selection in natural populations. Yet, direct tests of this prediction have been restricted to qualitative trends due to the lack of detailed characterization of recombination rate variation across genomes and within species. The use of imprecise recombination rates can also skew population genetic analyses designed to assess the presence and mode of selection across genomes. Here we report the first integrated high-resolution description of genomic and population variation in recombination, which also distinguishes between the two outcomes of meiotic recombination: crossing over (CO) and gene conversion (GC). We characterized the products of 5,860 female meioses in Drosophila melanogaster by genotyping a total of 139 million informative SNPs and mapped 106,964 recombination events at a resolution down to 2 kilobases. This approach allowed us to generate whole-genome CO and GC maps as well as a detailed description of variation in recombination among individuals of this species. We describe many levels of variation in recombination rates. At a large-scale (100 kb), CO rates exhibit extreme and highly punctuated variation along chromosomes, with hot and coldspots. We also show extensive intra-specific variation in CO landscapes that is associated with hotspots at low frequency in our sample. GC rates are more uniformly distributed across the genome than CO rates and detectable in regions with reduced or absent CO. At a local scale, recombination events are associated with numerous sequence motifs and tend to occur within transcript regions, thus suggesting that chromatin accessibility favors double-strand breaks. All these non-independent layers of variation in recombination across genomes and among individuals need to be taken into account in order to obtain relevant estimates of recombination rates, and should be included in a new generation of population genetic models of the interaction between selection and linkage.

Geographic variation in Drosophila: From molecules to morphology and back

The examination of spatial variation can act as a substitute for temporal variation in population studies. Of particular use in the effort to understand the selective forces that govern spatial variation are the sibling species Drosophila melanogaster and D. simulans. Recent work at the level of both molecules and morphology has uncovered a great deal of spatial variation within and between these two species. Here we summarize these data, with reference to what they tell us about the history and nature of selection operating in these species.

Patterns of neutral diversity under general models of selective sweeps

Two major sources of stochasticity in the dynamics of neutral alleles result from resampling of finite populations (genetic drift) and the random genetic background of nearby selected alleles on which the neutral alleles are found (linked selection). There is now good evidence that linked selection plays an important role in shaping polymorphism levels in a number of species. One of the best-investigated models of linked selection is the recurrent full-sweep model, in which newly arisen selected alleles fix rapidly. However, the bulk of selected alleles that sweep into the population may not be destined for rapid fixation. Here we develop a general model of recurrent selective sweeps in a coalescent framework, one that generalizes the recurrent full-sweep model to the case where selected alleles do not sweep to fixation. We show that in a large population, only the initial rapid increase of a selected allele affects the genealogy at partially linked sites, which under fairly general assumptions are unaffected by the subsequent fate of the selected allele. We also apply the theory to a simple model to investigate the impact of recurrent partial sweeps on levels of neutral diversity and find that for a given reduction in diversity, the impact of recurrent partial sweeps on the frequency spectrum at neutral sites is determined primarily by the frequencies rapidly achieved by the selected alleles. Consequently, recurrent sweeps of selected alleles to low frequencies can have a profound effect on levels of diversity but can leave the frequency spectrum relatively unperturbed. In fact, the limiting coalescent model under a high rate of sweeps to low frequency is identical to the standard neutral model. The general model of selective sweeps we describe goes some way toward providing a more flexible framework to describe genomic patterns of diversity than is currently available.

The dynamics of HIV-1 adaptation in early infection

Human immunodeficiency virus type 1 (HIV-1) undergoes a severe population bottleneck during sexual transmission and yet adapts extremely rapidly to the earliest immune responses. The bottleneck has been inferred to typically consist of a single genome, and typically eight amino acid mutations in viral proteins spread to fixation by the end of the early chronic phase of infection in response to selection by CD8(+) T cells. Stochastic simulation was used to examine the effects of the transmission bottleneck and of potential interference among spreading immune-escape mutations on the adaptive dynamics of the virus in early infection. If major viral population genetic parameters are assigned realistic values that permit rapid adaptive evolution, then a bottleneck of a single genome is not inconsistent with the observed pattern of adaptive fixations. One requirement is strong selection by CD8(+) T cells that decreases over time. Such selection may reduce effective population sizes at linked loci through genetic hitchhiking. However, this effect is predicted to be minor in early infection because the transmission bottleneck reduces the effective population size to such an extent that the resulting strong selection and weak mutation cause beneficial mutations to fix sequentially and thus avoid interference.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Signatures of selection in the genomes of commercial and non-commercial chicken breeds

Identifying genomics regions that are affected by selection is important to understand the domestication and selection history of the domesticated chicken, as well as understanding molecular pathways underlying phenotypic traits and breeding goals. While whole-genome approaches, either high-density SNP chips or massively parallel sequencing, have been successfully applied to identify evidence for selective sweeps in chicken, it has been difficult to distinguish patterns of selection and stochastic and breed specific effects. Here we present a study to identify selective sweeps in a large number of chicken breeds (67 in total) using a high-density (58 K) SNP chip. We analyzed commercial chickens representing all major breeding goals. In addition, we analyzed non-commercial chicken diversity for almost all recognized traditional Dutch breeds and a selection of representative breeds from China. Based on their shared history or breeding goal we in silico grouped the breeds into 14 breed groups. We identified 396 chromosomal regions that show suggestive evidence of selection in at least one breed group with 26 of these regions showing strong evidence of selection. Of these 26 regions, 13 were previously described and 13 yield new candidate genes for performance traits in chicken. Our approach demonstrates the strength of including many different populations with similar, and breed groups with different selection histories to reduce stochastic effects based on single populations.

The joint effects of background selection and genetic recombination on local gene genealogies

Background selection, the effects of the continual removal of deleterious mutations by natural selection on variability at linked sites, is potentially a major determinant of DNA sequence variability. However, the joint effects of background selection and genetic recombination on the shape of the neutral gene genealogy have proved hard to study analytically. The only existing formula concerns the mean coalescent time for a pair of alleles, making it difficult to assess the importance of background selection from genome-wide data on sequence polymorphism. Here we develop a structured coalescent model of background selection with recombination and implement it in a computer program that efficiently generates neutral gene genealogies for an arbitrary sample size. We check the validity of the structured coalescent model against forward-in-time simulations and show that it accurately captures the effects of background selection. The model produces more accurate predictions of the mean coalescent time than the existing formula and supports the conclusion that the effect of background selection is greater in the interior of a deleterious region than at its boundaries. The level of linkage disequilibrium between sites is elevated by background selection, to an extent that is well summarized by a change in effective population size. The structured coalescent model is readily extendable to more realistic situations and should prove useful for analyzing genome-wide polymorphism data.

Temporal variations in the hepatitis C virus intrahost population during chronic infection

The intrahost evolution of hepatitis C virus (HCV) holds keys to understanding mechanisms responsible for the establishment of chronic infections and to development of a vaccine and therapeutics. In this study, intrahost variants of two variable HCV genomic regions, HVR1 and NS5A, were sequenced from four treatment-naïve chronically infected patients who were followed up from the acute stage of infection for 9 to 18 years. Median-joining network analysis indicated that the majority of the HCV intrahost variants were observed only at certain time points, but some variants were detectable at more than one time point. In all patients, these variants were found organized into communities or subpopulations. We hypothesize that HCV intrahost evolution is defined by two processes: incremental changes within communities through random mutation and alternations between coexisting communities. The HCV population was observed to incrementally evolve within a single community during approximately the first 3 years of infection, followed by dispersion into several subpopulations. Two patients demonstrated this pattern of dispersion for the rest of the observation period, while HCV variants in the other two patients converged into another single subpopulation after ∼9 to 12 years of dispersion. The final subpopulation in these two patients was under purifying selection. Intrahost HCV evolution in all four patients was characterized by a consistent increase in negative selection over time, suggesting the increasing HCV adaptation to the host late in infection. The data suggest specific staging of HCV intrahost evolution.

The influence of recombination on SNP diversity in chickens

The association between recombination rate and diversity, but not divergence is considered to be driven mainly by natural selection: fixation of positively selected variants and associated hitchhiking effects and/or background selection eliminating deleterious alleles. In the present study, we investigated the relationship between recombination rate, SNP diversity and interspecies divergence for 29 loci in chickens. We found that recombination rate is positively correlated with nucleotide diversity but is not correlated with interspecies divergence. It appears that variation in recombination rate explains over 30% of the variation in levels of diversity among 29 loci. Our data suggested that natural selection is a main factor in shaping SNP diversity in chickens. Since SNP diversity is significantly lower at Z-linked than at autosomal loci, we argued that genetic hitchhiking might be more important than background selection in producing the observed correlation.

Molecular population genomics: a short history

Population genomics is the study of the amount and causes of genome-wide variability in natural populations, a topic that has been under discussion since Darwin. This paper first briefly reviews the early development of molecular approaches to the subject: the pioneering unbiased surveys of genetic variability at multiple loci by means of gel electrophoresis and restriction enzyme mapping. The results of surveys of levels of genome-wide variability using DNA resequencing studies are then discussed. Studies of the extent to which variability for different classes of variants (non-synonymous, synonymous and non-coding) are affected by natural selection, or other directional forces such as biased gene conversion, are also described. Finally, the effects of deleterious mutations on population fitness and the possible role of Hill–Robertson interference in shaping patterns of sequence variability are discussed.

Pervasive adaptive protein evolution apparent in diversity patterns around amino acid substitutions in Drosophila simulans

In Drosophila, multiple lines of evidence converge in suggesting that beneficial substitutions to the genome may be common. All suffer from confounding factors, however, such that the interpretation of the evidence—in particular, conclusions about the rate and strength of beneficial substitutions—remains tentative. Here, we use genome-wide polymorphism data in D. simulans and sequenced genomes of its close relatives to construct a readily interpretable characterization of the effects of positive selection: the shape of average neutral diversity around amino acid substitutions. As expected under recurrent selective sweeps, we find a trough in diversity levels around amino acid but not around synonymous substitutions, a distinctive pattern that is not expected under alternative models. This characterization is richer than previous approaches, which relied on limited summaries of the data (e.g., the slope of a scatter plot), and relates to underlying selection parameters in a straightforward way, allowing us to make more reliable inferences about the prevalence and strength of adaptation. Specifically, we develop a coalescent-based model for the shape of the entire curve and use it to infer adaptive parameters by maximum likelihood. Our inference suggests that ∼13% of amino acid substitutions cause selective sweeps. Interestingly, it reveals two classes of beneficial fixations: a minority (approximately 3%) that appears to have had large selective effects and accounts for most of the reduction in diversity, and the remaining 10%, which seem to have had very weak selective effects. These estimates therefore help to reconcile the apparent conflict among previously published estimates of the strength of selection. More generally, our findings provide unequivocal evidence for strongly beneficial substitutions in Drosophila and illustrate how the rapidly accumulating genome-wide data can be leveraged to address enduring questions about the genetic basis of adaptation.

Characterizing the nature of beneficial changes to the genome is essential to our understanding of adaptation. To do so, researchers identify and analyze footprints that beneficial changes leave in patterns of genetic variation within and between species. In order to teach us about adaptive evolution, these footprints need to be specific to positive selection as well as rich enough to allow for reliable inferences. Here, we identify such a footprint: a pronounced trough in the average levels of genetic diversity surrounding amino acid substitutions throughout the D. simulans genome. Based on this pattern, we infer that approximately 13% of amino acid substitutions were beneficial, a minority of which (3%) conferred a large selective advantage of nearly 0.5% and the majority of which (10%) conferred a much smaller advantage of about 0.01%. These findings offer insights into the distribution of selection effects driving beneficial changes to the D. simulans genome and suggest how the widely varying estimates obtained in previous studies of Drosophila may be reconciled. Moreover, the approach that we introduce is readily applicable to other taxa and thus should help to gain important insights into how the rate and strength of adaptive evolution vary depending on life-history, population size, and ecology.

Population genomics of transposable elements in Drosophila melanogaster

Transposable elements (TEs) are the primary contributors to the genome bulk in many organisms and are major players in genome evolution. A clear and thorough understanding of the population dynamics of TEs is therefore essential for full comprehension of the eukaryotic genome evolution and function. Although TEs in Drosophila melanogaster have received much attention, population dynamics of most TE families in this species remains entirely unexplored. It is not clear whether the same population processes can account for the population behaviors of all TEs in Drosophila or whether, as has been suggested previously, different orders behave according to very different rules. In this work, we analyzed population frequencies for a large number of individual TEs (755 TEs) in five North American and one sub-Saharan African D. melanogaster populations (75 strains in total). These TEs have been annotated in the reference D. melanogaster euchromatic genome and have been sampled from all three major orders (non-LTR, LTR, and TIR) and from all families with more than 20 TE copies (55 families in total). We find strong evidence that TEs in Drosophila across all orders and families are subject to purifying selection at the level of ectopic recombination. We showed that strength of this selection varies predictably with recombination rate, length of individual TEs, and copy number and length of other TEs in the same family. Importantly, these rules do not appear to vary across orders. Finally, we built a statistical model that considered only individual TE-level (such as the TE length) and family-level properties (such as the copy number) and were able to explain more than 40% of the variation in TE frequencies in D. melanogaster.

Molecular evolution of a Y chromosome to autosome gene duplication in Drosophila

In contrast to the rest of the genome, the Y chromosome is restricted to males and lacks recombination. As a result, Y chromosomes are unable to respond efficiently to selection, and newly formed Y chromosomes degenerate until few genes remain. The rapid loss of genes from newly formed Y chromosomes has been well studied, but gene loss from highly degenerate Y chromosomes has only recently received attention. Here, we identify and characterize a Y to autosome duplication of the male fertility gene kl-5 that occurred during the evolution of the testacea group species of Drosophila. The duplication was likely DNA based, as other Y-linked genes remain on the Y chromosome, the locations of introns are conserved, and expression analyses suggest that regulatory elements remain linked. Genetic mapping reveals that the autosomal copy of kl-5 resides on the dot chromosome, a tiny autosome with strongly suppressed recombination. Molecular evolutionary analyses show that autosomal copies of kl-5 have reduced polymorphism and little recombination. Importantly, the rate of protein evolution of kl-5 has increased significantly in lineages where it is on the dot versus Y linked. Further analyses suggest this pattern is a consequence of relaxed purifying selection, rather than adaptive evolution. Thus, although the initial fixation of the kl-5 duplication may have been advantageous, slightly deleterious mutations have accumulated in the dot-linked copies of kl-5 faster than in the Y-linked copies. Because the dot chromosome contains seven times more genes than the Y and is exposed to selection in both males and females, these results suggest that the dot suffers the deleterious effects of genetic linkage to more selective targets compared with the Y chromosome. Thus, a highly degenerate Y chromosome may not be the worst environment in the genome, as is generally thought, but may in fact be protected from the accumulation of deleterious mutations relative to other nonrecombining regions that contain more genes.

Genome-wide effects of long-term divergent selection

To understand the genetic mechanisms leading to phenotypic differentiation, it is important to identify genomic regions under selection. We scanned the genome of two chicken lines from a single trait selection experiment, where 50 generations of selection have resulted in a 9-fold difference in body weight. Analyses of nearly 60,000 SNP markers showed that the effects of selection on the genome are dramatic. The lines were fixed for alternative alleles in more than 50 regions as a result of selection. Another 10 regions displayed strong evidence for ongoing differentiation during the last 10 generations. Many more regions across the genome showed large differences in allele frequency between the lines, indicating that the phenotypic evolution in the lines in 50 generations is the result of an exploitation of standing genetic variation at 100s of loci across the genome.

Nonrecombining genes in a recombination environment: the Drosophila "dot" chromosome

Rate of recombination is a powerful variable affecting several aspects of molecular variation and evolution. A nonrecombining portion of the genome of most Drosophila species, the "dot" chromosome or F element, exhibits very low levels of variation and unusual codon usage. One lineage of Drosophila, the willistoni/saltans groups, has the F element fused to a normally recombining E element. Here, we present polymorphism data for genes on the F element in two Drosophila willistoni and one D. insularis populations, genes previously studied in D. melanogaster. The D. willistoni populations were known to be very low in inversion polymorphism, thus minimizing the recombination suppression effect of inversions. We first confirmed, by in situ hybridization, that D. insularis has the same E + F fusion as D. willistoni, implying this was a monophyletic event. A clear gradient in codon usage exists along the willistoni F element, from the centromere distally to the fusion with E; estimates of recombination rates parallel this gradient and also indicate D. insularis has greater recombination than D. willistoni. In contrast to D. melanogaster, genes on the F element exhibit moderate levels of nucleotide polymorphism not distinguishable from two genes elsewhere in the genome. Although some linkage disequilibrium (LD) was detected between polymorphic sites within genes (generally <500 bp apart), no long-range LD between F element loci exists in the two willistoni group species. In general, the distribution of allele frequencies of F element genes display the typical pattern of expectations of neutral variation at equilibrium. These results are consistent with the hypothesis that recombination allows the accumulation of nucleotide variation as well as allows selection to act on synonymous codon usage. It is estimated that the fusion occurred ∼20 Mya and while the F element in the willistoni lineage has evolved "normal" levels and patterns of nucleotide variation, equilibrium may not have been reached for codon usage.