Natural selection shapes nucleotide polymorphism across the genome of the nematode Caenorhabditis briggsae

Accurate inference of population history in the presence of background selection

All published methods for learning about demographic history make the simplifying assumption that the genome evolves neutrally, and do not seek to account for the effects of natural selection on patterns of variation. This is a major concern, as ample work has demonstrated the pervasive effects of natural selection and in particular background selection (BGS) on patterns of genetic variation in diverse species. Simulations and theoretical work have shown that methods to infer changes in effective population size over time (Ne(t)) become increasingly inaccurate as the strength of linked selection increases. Here, we introduce an extension to the Pairwise Sequentially Markovian Coalescent (PSMC) algorithm, PSMC+, which explicitly co-models demographic history and natural selection. We benchmark our method using forward-in-time simulations with BGS and find that our approach improves the accuracy of effective population size inference. Leveraging a high resolution map of BGS in humans, we infer considerable changes in the magnitude of inferred effective population size relative to previous reports. Finally, we separately infer Ne(t) on the X chromosome and on the autosomes in diverse great apes without making a correction for selection, and find that the inferred ratio fluctuates substantially through time in a way that differs across species, showing that uncorrected selection may be an important driver of signals of genetic difference on the X chromosome and autosomes.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Genetic diversity of a short-ranged endemic terrestrial snail

The factors that influence population structure and connectivity are unknown for most terrestrial invertebrates but are of particular interest both for understanding the impacts of disturbance and for determining accurate levels of biodiversity and local endemism. The main objective of this study was to determine the historical patterns of genetic differentiation and contemporary gene flow in the terrestrial snail, Austrochloritis kosciuszkoensis (Shea & O. L. Griffiths, 2010). Snails were collected in the Mt Buffalo and Alpine National Parks in Victoria, in a bid to understand how populations of this species are connected both within continuous habitat and between adjacent, yet separate environments. Utilising both mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, the degree of population structure was determined within and between sites. Very high levels of genetic divergence were found between the Mt Buffalo and Alpine snails, with no evidence for genetic exchange detected between the two regions, indicating speciation has possibly occurred between the two regions. Our analyses of the combined mtDNA and nDNA (generated from SNPs) data have revealed patterns of genetic diversity that are consistent with a history of long-term isolation and limited connectivity. This history may be related to past cycles of changes to the climate over hundreds of thousands of years, which have, in part, caused the fragmentation of Australian forests. Within both regions, extremely limited gene flow between separate populations suggests that these land snails have very limited dispersal capabilities across existing landscape barriers, especially at Mt Buffalo: here, populations only 5 km apart from each other are genetically differentiated. The distinct genetic divergences and clearly reduced dispersal ability detected in this data explain the likely existence of at least two previously unnamed cryptic Austrochloritis species within a 30-50 km radius, and highlight the need for more concentrated efforts to understand population structure and gene flow in terrestrial invertebrates.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Cold Survival and Its Molecular Mechanisms in a Locally Adapted Nematode Population

Since Darwin, evolutionary biologists have sought to understand the drivers and mechanisms of natural trait diversity. The field advances toward this goal with the discovery of phenotypes that vary in the wild, their relationship to ecology, and their underlying genes. Here, we established resistance to extreme low temperature in the free-living nematode Caenorhabditis briggsae as an ecological and evolutionary model system. We found that C. briggsae strains of temperate origin were strikingly more cold-resistant than those isolated from tropical localities. Transcriptional profiling revealed expression patterns unique to the resistant temperate ecotype, including dozens of genes expressed at high levels even after multiple days of cold-induced physiological slowdown. Mutational analysis validated a role in cold resistance for seven such genes. These findings highlight a candidate case of robust, genetically complex adaptation in an emerging model nematode, and shed light on the mechanisms at play.

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a "faster male" model for Haldane's rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a "large-X effect" for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

A model of autowave self-organization as a hierarchy of active media in the biological evolution

Within the framework of the active media concept, we develop a biophysical model of autowave self-organization which is treated as a hierarchy of active media in the evolution of the biosphere. We also propose a mathematical model of the autowave process of speciation in a flow of mutations for the three main taxonometric groups (prokaryotes, unicellular and multicellular eukaryotes) with a naturally determined lower boundary of living matter (the appearance of prokaryotes) and an open upper boundary for the formation of new species. It is shown that the fluctuation-bifurcation description of the evolution for the formation of new taxonometric groups as a trajectory of transformation of small fluctuations into giant ones adequately reflects the process of self-organization during the formation of taxa. The major concepts of biological evolution, conditions of hierarchy formation as a fundamental manifestation of self-organization and complexity in the evolution of biological systems are considered.

Estimating the genome-wide contribution of selection to temporal allele frequency change

Rapid phenotypic adaptation is often observed in natural populations and selection experiments. However, detecting the genome-wide impact of this selection is difficult since adaptation often proceeds from standing variation and selection on polygenic traits, both of which may leave faint genomic signals indistinguishable from a noisy background of genetic drift. One promising signal comes from the genome-wide covariance between allele frequency changes observable from temporal genomic data (e.g., evolve-and-resequence studies). These temporal covariances reflect how heritable fitness variation in the population leads changes in allele frequencies at one time point to be predictive of the changes at later time points, as alleles are indirectly selected due to remaining associations with selected alleles. Since genetic drift does not lead to temporal covariance, we can use these covariances to estimate what fraction of the variation in allele frequency change through time is driven by linked selection. Here, we reanalyze three selection experiments to quantify the effects of linked selection over short timescales using covariance among time points and across replicates. We estimate that at least 17 to 37% of allele frequency change is driven by selection in these experiments. Against this background of positive genome-wide temporal covariances, we also identify signals of negative temporal covariance corresponding to reversals in the direction of selection for a reasonable proportion of loci over the time course of a selection experiment. Overall, we find that in the three studies we analyzed, linked selection has a large impact on short-term allele frequency dynamics that is readily distinguishable from genetic drift.

Multi-Scale Landscape Influences on Genetic Diversity and Adaptive Traits in a Neotropical Savanna Tree

Changes in landscape structure can affect essential population ecological features, such as dispersal and recruitment, and thus genetic processes. Here, we analyze the effects of landscape metrics on adaptive quantitative traits variation, evolutionary potential, and on neutral genetic diversity in populations of the Neotropical savanna tree Tabebuia aurea. Using a multi-scale approach, we sampled five landscapes with two sites of savanna in each. To obtain neutral genetic variation, we genotyped 60 adult individuals from each site using 10 microsatellite loci. We measured seed size and mass. Seeds were grown in nursery in completely randomized experimental design and 17 traits were measured in seedlings to obtain the average, additive genetic variance (V a ) and coefficient of variation (CV a %), which measures evolvability, for each trait. We found that habitat loss increased genetic diversity (He) and allelic richness (AR), and decreased genetic differentiation among populations (F ST ), most likely due to longer dispersal distance of pollen in landscapes with lower density of flowering individuals. Habitat amount positively influenced seed size. Seeds of T. aurea are wind-dispersed and larger seeds may be dispersed to short distance, increasing genetic differentiation and decreasing genetic diversity and allelic richness. Evolvability (CV a %) in root length decreased with habitat amount. Savanna trees have higher root than shoot growth rate in the initial stages, allowing seedlings to obtain water from water tables. Landscapes with lower habitat amount may be more stressful for plant species, due to the lower plant density, edge effects and the negative impacts of agroecosystems. In these landscapes, larger roots may provide higher ability to obtain water, increasing survival and avoiding dying back because of fire. Despite the very recent agriculture expansion in Central Brazil, landscape changes are affecting neutral and adaptive variation in T. aurea. Several populations have low additive genetic variation for some traits and thus, may have limited evolvability, which may jeopardize species long-term persistence. The effect of habitat loss on highly variable neutral loci may only be detected after a certain threshold of population size is attained, that could become dangerously small masking important losses of heterozygosity endangering species conservation.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Genome structure predicts modular transcriptome responses to genetic and environmental conditions

Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species

Population genetic theory predicts that selection should be more effective when the effective population size (N_e) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger N_e. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.

Background selection as null hypothesis in population genomics: insights and challenges from Drosophila studies

The consequences of selection at linked sites are multiple and widespread across the genomes of most species. Here, I first review the main concepts behind models of selection and linkage in recombining genomes, present the difficulty in parametrizing these models simply as a reduction in effective population size (N_e) and discuss the predicted impact of recombination rates on levels of diversity across genomes. Arguments are then put forward in favour of using a model of selection and linkage with neutral and deleterious mutations (i.e. the background selection model, BGS) as a sensible null hypothesis for investigating the presence of other forms of selection, such as balancing or positive. I also describe and compare two studies that have generated high-resolution landscapes of the predicted consequences of selection at linked sites in Drosophila melanogaster. Both studies show that BGS can explain a very large fraction of the observed variation in diversity across the whole genome, thus supporting its use as null model. Finally, I identify and discuss a number of caveats and challenges in studies of genetic hitchhiking that have been often overlooked, with several of them sharing a potential bias towards overestimating the evidence supporting recent selective sweeps to the detriment of a BGS explanation. One potential source of bias is the analysis of non-equilibrium populations: it is precisely because models of selection and linkage predict variation in N_e across chromosomes that demographic dynamics are not expected to be equivalent chromosome- or genome-wide. Other challenges include the use of incomplete genome annotations, the assumption of temporally stable recombination landscapes, the presence of genes under balancing selection and the consequences of ignoring non-crossover (gene conversion) recombination events.

This article is part of the themed issue ‘Evolutionary causes and consequences of recombination rate variation in sexual organisms’.

Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.

Effect of mutation mechanisms on variant composition and distribution in Caenorhabditis elegans

Genetic diversity is maintained by continuing generation and removal of variants. While examining over 800,000 DNA variants in wild isolates of Caenorhabditis elegans, we made a discovery that the proportions of variant types are not constant across the C. elegans genome. The variant proportion is defined as the fraction of a specific variant type (e.g. single nucleotide polymorphism (SNP) or indel) within a broader set of variants (e.g. all variants or all non-SNPs). The proportions of most variant types show a correlation with the recombination rate. These correlations can be explained as a result of a concerted action of two mutation mechanisms, which we named Morgan and Sanger mechanisms. The two proposed mechanisms act according to the distinct components of the recombination rate, specifically the genetic and physical distance. Regression analysis was used to explore the characteristics and contributions of the two mutation mechanisms. According to our model, ~20-40% of all mutations in C. elegans wild populations are derived from programmed meiotic double strand breaks, which precede chromosomal crossovers and thus may be the point of origin for the Morgan mechanism. A substantial part of the known correlation between the recombination rate and variant distribution appears to be caused by the mutations generated by the Morgan mechanism. Mathematically integrating the mutation model with background selection model gives a more complete depiction of how the variant landscape is shaped in C. elegans. Similar analysis should be possible in other species by examining the correlation between the recombination rate and variant landscape within the context of our mutation model.

Determining the Effect of Natural Selection on Linked Neutral Divergence across Species

A major goal in evolutionary biology is to understand how natural selection has shaped patterns of genetic variation across genomes. Studies in a variety of species have shown that neutral genetic diversity (intra-species differences) has been reduced at sites linked to those under direct selection. However, the effect of linked selection on neutral sequence divergence (inter-species differences) remains ambiguous. While empirical studies have reported correlations between divergence and recombination, which is interpreted as evidence for natural selection reducing linked neutral divergence, theory argues otherwise, especially for species that have diverged long ago. Here we address these outstanding issues by examining whether natural selection can affect divergence between both closely and distantly related species. We show that neutral divergence between closely related species (e.g. human-primate) is negatively correlated with functional content and positively correlated with human recombination rate. We also find that neutral divergence between distantly related species (e.g. human-rodent) is negatively correlated with functional content and positively correlated with estimates of background selection from primates. These patterns persist after accounting for the confounding factors of hypermutable CpG sites, GC content, and biased gene conversion. Coalescent models indicate that even when the contribution of ancestral polymorphism to divergence is small, background selection in the ancestral population can still explain a large proportion of the variance in divergence across the genome, generating the observed correlations. Our findings reveal that, contrary to previous intuition, natural selection can indirectly affect linked neutral divergence between both closely and distantly related species. Though we cannot formally exclude the possibility that the direct effects of purifying selection drive some of these patterns, such a scenario would be possible only if more of the genome is under purifying selection than currently believed. Our work has implications for understanding the evolution of genomes and interpreting patterns of genetic variation.

Genetic variation at neutral sites can be reduced through linkage to nearby selected sites. This pattern has been used to show the widespread effects of natural selection at shaping patterns of genetic diversity across genomes from a variety of species. However, it is not entirely clear whether natural selection has an effect on neutral divergence between species. Here we show that putatively neutral divergence between closely related species (human and chimp) and between distantly related pairs of species (humans and mice) show signatures consistent with having been affected by linkage to selected sites. Further, our theoretical models and simulations show that natural selection indirectly affecting linked neutral sites can generate these patterns. Unless substantially more of the genome is under the direct effects of purifying selection than currently believed, our results argue that natural selection has played an important role in shaping variation in levels of putatively neutral sequence divergence across the genome. Our findings further suggest that divergence-based estimates of neutral mutation rate variation across the genome as well as certain estimators of population history may be confounded by linkage to selected sites.

Determinants of genetic diversity

Genetic polymorphism varies among species and within genomes, and has important implications for the evolution and conservation of species. The determinants of this variation have been poorly understood, but population genomic data from a wide range of organisms now make it possible to delineate the underlying evolutionary processes, notably how variation in the effective population size (Ne) governs genetic diversity. Comparative population genomics is on its way to providing a solution to 'Lewontin's paradox' - the discrepancy between the many orders of magnitude of variation in population size and the much narrower distribution of diversity levels. It seems that linked selection plays an important part both in the overall genetic diversity of a species and in the variation in diversity within the genome. Genetic diversity also seems to be predictable from the life history of a species.

RhesusBase PopGateway: Genome-Wide Population Genetics Atlas in Rhesus Macaque

Although population genetics studies have significantly accelerated the evolutionary and functional interrogations of genes and regulations, limited polymorphism data are available for rhesus macaque, the model animal closely related to human. Here, we report the first genome-wide effort to identify and visualize the population genetics profile in rhesus macaque. On the basis of the whole-genome sequencing of 31 independent macaque animals, we profiled a comprehensive polymorphism map with 46,146,548 sites. The allele frequency for each polymorphism site, the haplotype structure, as well as multiple population genetics parameters were then calculated on a genome-wide scale. We further developed a specific interface, the RhesusBase PopGateway, to facilitate the visualization of these annotations, and highlighted the applications of this highly integrative platform in clarifying the selection signatures of genes and regulations in the context of the primate evolution. Overall, the updated RhesusBase provides a comprehensive monkey population genetics framework for in-depth evolutionary studies of human biology.

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and FST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.

Natural selection constrains neutral diversity across a wide range of species

The neutral theory of molecular evolution predicts that the amount of neutral polymorphisms within a species will increase proportionally with the census population size (Nc). However, this prediction has not been borne out in practice: while the range of Nc spans many orders of magnitude, levels of genetic diversity within species fall in a comparatively narrow range. Although theoretical arguments have invoked the increased efficacy of natural selection in larger populations to explain this discrepancy, few direct empirical tests of this hypothesis have been conducted. In this work, we provide a direct test of this hypothesis using population genomic data from a wide range of taxonomically diverse species. To do this, we relied on the fact that the impact of natural selection on linked neutral diversity depends on the local recombinational environment. In regions of relatively low recombination, selected variants affect more neutral sites through linkage, and the resulting correlation between recombination and polymorphism allows a quantitative assessment of the magnitude of the impact of selection on linked neutral diversity. By comparing whole genome polymorphism data and genetic maps using a coalescent modeling framework, we estimate the degree to which natural selection reduces linked neutral diversity for 40 species of obligately sexual eukaryotes. We then show that the magnitude of the impact of natural selection is positively correlated with Nc, based on body size and species range as proxies for census population size. These results demonstrate that natural selection removes more variation at linked neutral sites in species with large Nc than those with small Nc and provides direct empirical evidence that natural selection constrains levels of neutral genetic diversity across many species. This implies that natural selection may provide an explanation for this longstanding paradox of population genetics.

Full-genome evolutionary histories of selfing, splitting, and selection in Caenorhabditis

The nematode Caenorhabditis briggsae is a model for comparative developmental evolution with C. elegans. Worldwide collections of C. briggsae have implicated an intriguing history of divergence among genetic groups separated by latitude, or by restricted geography, that is being exploited to dissect the genetic basis to adaptive evolution and reproductive incompatibility; yet, the genomic scope and timing of population divergence is unclear. We performed high-coverage whole-genome sequencing of 37 wild isolates of the nematode C. briggsae and applied a pairwise sequentially Markovian coalescent (PSMC) model to 703 combinations of genomic haplotypes to draw inferences about population history, the genomic scope of natural selection, and to compare with 40 wild isolates of C. elegans. We estimate that a diaspora of at least six distinct C. briggsae lineages separated from one another approximately 200,000 generations ago, including the “Temperate” and “Tropical” phylogeographic groups that dominate most samples worldwide. Moreover, an ancient population split in its history approximately 2 million generations ago, coupled with only rare gene flow among lineage groups, validates this system as a model for incipient speciation. Low versus high recombination regions of the genome give distinct signatures of population size change through time, indicative of widespread effects of selection on highly linked portions of the genome owing to extreme inbreeding by self-fertilization. Analysis of functional mutations indicates that genomic context, owing to selection that acts on long linkage blocks, is a more important driver of population variation than are the functional attributes of the individually encoded genes.

The effects of purifying selection on patterns of genetic differentiation between Drosophila melanogaster populations

Using the data provided by the Drosophila Population Genomics Project, we investigate factors that affect the genetic differentiation between Rwandan and French populations of D. melanogaster. By examining within-population polymorphisms, we show that sites in long introns (especially those >2000 bp) have significantly lower π (nucleotide diversity) and more low-frequency variants (as measured by Tajima's D, minor allele frequencies, and prevalence of variants that are private to one of the two populations) than short introns, suggesting a positive relationship between intron length and selective constraint. A similar analysis of protein-coding polymorphisms shows that 0-fold (degenerate) sites in more conserved genes are under stronger purifying selection than those in less conserved genes. There is limited evidence that selection on codon bias has an effect on differentiation (as measured by F_ST) at 4-fold (degenerate) sites, and 4-fold sites and sites in 8–30 bp of short introns ⩽65 bp have comparable F_ST values. Consistent with the expected effect of purifying selection, sites in long introns and 0-fold sites in conserved genes are less differentiated than those in short introns and less conserved genes, respectively. Genes in non-crossover regions (for example, the fourth chromosome) have very high F_ST values at both 0-fold and 4-fold degenerate sites, which is probably because of the large reduction in within-population diversity caused by tight linkage between many selected sites. Our analyses also reveal subtle statistical properties of F_ST, which arise when information from multiple single nucleotide polymorphisms is combined and can lead to the masking of important signals of selection.

Characterization of genetic diversity in the nematode Pristionchus pacificus from population-scale resequencing data

The hermaphroditic nematode Pristionchus pacificus is an established model system for comparative studies with Caenorhabditis elegans in developmental biology, ecology, and population genetics. In this study, we present whole-genome sequencing data of 104 P. pacificus strains and the draft assembly of the obligate outcrossing sister species P. exspectatus. We characterize genetic diversity within P. pacificus and investigate the population genetic processes shaping this diversity. P. pacificus is 10 times more diverse than C. elegans and exhibits substantial population structure that allows us to probe its evolution on multiple timescales. Consistent with reduced effective recombination in this self-fertilizing species, we find haplotype blocks that span several megabases. Using the P. exspectatus genome as an outgroup, we polarized variation in P. pacificus and found a site frequency spectrum (SFS) that decays more rapidly than expected in neutral models. The SFS at putatively neutral sites is U shaped, which is a characteristic feature of pervasive linked selection. Based on the additional findings (i) that the majority of nonsynonymous variation is eliminated over timescales on the order of the separation between clades, (ii) that diversity is reduced in gene-rich regions, and (iii) that highly differentiated clades show very similar patterns of diversity, we conclude that purifying selection on many mutations with weak effects is a major force shaping genetic diversity in P. pacificus.

An ABC method for estimating the rate and distribution of effects of beneficial mutations

Determining the distribution of adaptive mutations available to natural selection is a difficult task. These are rare events and most of them are lost by chance. Some theoretical works propose that the distribution of newly arising beneficial mutations should be close to exponential. Empirical data are scarce and do not always support an exponential distribution. Analysis of the dynamics of adaptation in asexual populations of microorganisms has revealed that these can be summarized by two effective parameters, the effective mutation rate, U_e, and the effective selection coefficient of a beneficial mutation, S_e. Here, we show that these effective parameters will not always reflect the rate and mean effect of beneficial mutations, especially when the distribution of arising mutations has high variance, and the mutation rate is high. We propose a method to estimate the distribution of arising beneficial mutations, which is motivated by a common experimental setup. The method, which we call One Biallelic Marker Approximate Bayesian Computation, makes use of experimental data consisting of periodic measures of neutral marker frequencies and mean population fitness. Using simulations, we find that this method allows the discrimination of the shape of the distribution of arising mutations and that it provides reasonable estimates of their rates and mean effects in ranges of the parameter space that may be of biological relevance.

Patterns of Selection in Plant Genomes

Plants show a wide range of variation in mating system, ploidy level, and demographic history, allowing for unique opportunities to investigate the evolutionary and genetic factors affecting genome-wide patterns of positive and negative selection. In this review, we highlight recent progress in our understanding of the extent and nature of selection on plant genomes. We discuss differences in selection as they relate to variation in demography, recombination, mating system, and ploidy. We focus on the population genetic consequences of these factors and argue that, although variation in the magnitude of purifying selection is well documented, quantifying rates of positive selection and disentangling the relative importance of recombination, demography, and ploidy are ongoing challenges. Large-scale comparative studies that examine the relative and joint importance of these processes, combined with explicit models of population history and selection, are key and feasible goals for future work.

The Drosophila early ovarian transcriptome provides insight to the molecular causes of recombination rate variation across genomes

Background

Evidence in yeast indicates that gene expression is correlated with recombination activity and double-strand break (DSB) formation in some hotspots. Studies of nucleosome occupancy in yeast and mice also suggest that open chromatin influences the formation of DSBs. In Drosophila melanogaster, high-resolution recombination maps show an excess of DSBs within annotated transcripts relative to intergenic sequences. The impact of active transcription on recombination landscapes, however, remains unexplored in a multicellular organism. We then investigated the transcription profile during early meiosis in D. melanogaster females to obtain a glimpse at the relevant transcriptional dynamics during DSB formation, and test the specific hypothesis that DSBs preferentially target transcriptionally active genomic regions.

Results

Our study of transcript profiles of early- and late-meiosis using mRNA-seq revealed, 1) significant differences in gene expression, 2) new genes and exons, 3) parent-of-origin effects on transcription in early-meiosis stages, and 4) a nonrandom genomic distribution of transcribed genes. Importantly, genomic regions that are more actively transcribed during early meiosis show higher rates of recombination, and we ruled out DSB preference for genic regions that are not transcribed.

Conclusions

Our results provide evidence in a multicellular organism that transcription during the initial phases of meiosis increases the likelihood of DSB and give insight into the molecular determinants of recombination rate variation across the D. melanogaster genome. We propose that a model where variation in gene expression plays a role altering the recombination landscape across the genome could provide a molecular, heritable and plastic mechanism to observed patterns of recombination variation, from the high level of intra-specific variation to the known influence of environmental factors and stress conditions.

Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

Genomic signatures of selection at linked sites: unifying the disparity among species

Population genetics theory supplies powerful predictions about how natural selection interacts with genetic linkage to sculpt the genomic landscape of nucleotide polymorphism. Both the spread of beneficial mutations and the removal of deleterious mutations act to depress polymorphism levels, especially in low-recombination regions. However, empiricists have documented extreme disparities among species. Here we characterize the dominant features that could drive differences in linked selection among species--including roles for selective sweeps being 'hard' or 'soft'--and the concealing effects of demography and confounding genomic variables. We advocate targeted studies of closely related species to unify our understanding of how selection and linkage interact to shape genome evolution.

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.

Snpdat: easy and rapid annotation of results from de novo snp discovery projects for model and non-model organisms

Background

Single nucleotide polymorphisms (SNPs) are the most abundant genetic variant found in vertebrates and invertebrates. SNP discovery has become a highly automated, robust and relatively inexpensive process allowing the identification of many thousands of mutations for model and non-model organisms. Annotating large numbers of SNPs can be a difficult and complex process. Many tools available are optimised for use with organisms densely sampled for SNPs, such as humans. There are currently few tools available that are species non-specific or support non-model organism data.

Results

Here we present SNPdat, a high throughput analysis tool that can provide a comprehensive annotation of both novel and known SNPs for any organism with a draft sequence and annotation. Using a dataset of 4,566 SNPs identified in cattle using high-throughput DNA sequencing we demonstrate the annotations performed and the statistics that can be generated by SNPdat.

Conclusions

SNPdat provides users with a simple tool for annotation of genomes that are either not supported by other tools or have a small number of annotated SNPs available. SNPdat can also be used to analyse datasets from organisms which are densely sampled for SNPs. As a command line tool it can easily be incorporated into existing SNP discovery pipelines and fills a niche for analyses involving non-model organisms that are not supported by many available SNP annotation tools. SNPdat will be of great interest to scientists involved in SNP discovery and analysis projects, particularly those with limited bioinformatics experience.

Adaptive evolution: evaluating empirical support for theoretical predictions

Adaptive evolution is shaped by the interaction of population genetics, natural selection and underlying network and biochemical constraints. Variation created by mutation, the raw material for evolutionary change, is translated into phenotypes by flux through metabolic pathways and by the topography and dynamics of molecular networks. Finally, the retention of genetic variation and the efficacy of selection depend on population genetics and demographic history. Emergent high-throughput experimental methods and sequencing technologies allow us to gather more evidence and to move beyond the theory in different systems and populations. Here we review the extent to which recent evidence supports long-established theoretical principles of adaptation.

A coalescent model of background selection with recombination, demography and variation in selection coefficients

There is increasing evidence that background selection, the effects of the elimination of recurring deleterious mutations by natural selection on variability at linked sites, may be a major factor shaping genome-wide patterns of genetic diversity. To accurately quantify the importance of background selection, it is vital to have computationally efficient models that include essential biological features. To this end, a structured coalescent procedure is used to construct a model of background selection that takes into account the effects of recombination, recent changes in population size and variation in selection coefficients against deleterious mutations across sites. Furthermore, this model allows a flexible organization of selected and neutral sites in the region concerned, and has the ability to generate sequence variability at both selected and neutral sites, allowing the correlation between these two types of sites to be studied. The accuracy of the model is verified by checking against the results of forward simulations. These simulations also reveal several patterns of diversity that are in qualitative agreement with observations reported in recent studies of DNA sequence polymorphisms. These results suggest that the model should be useful for data analysis.

Temperature-dependent behaviours are genetically variable in the nematode Caenorhabditis briggsae

Temperature-dependent behaviours in Caenorhabditis elegans, such as thermotaxis and isothermal tracking, are complex behavioural responses that integrate sensation, foraging and learning, and have driven investigations to discover many essential genetic and neural pathways. The ease of manipulation of the Caenorhabditis model system also has encouraged its application to comparative analyses of phenotypic evolution, particularly contrasts of the classic model C. elegans with C. briggsae. And yet few studies have investigated natural genetic variation in behaviour in any nematode. Here we measure thermotaxis and isothermal tracking behaviour in genetically distinct strains of C. briggsae, further motivated by the latitudinal differentiation in C. briggsae that is associated with temperature-dependent fitness differences in this species. We demonstrate that C. briggsae performs thermotaxis and isothermal tracking largely similar to that of C. elegans, with a tendency to prefer its rearing temperature. Comparisons of these behaviours among strains reveal substantial heritable natural variation within each species that corresponds to three general patterns of behavioural response. However, intraspecific genetic differences in thermal behaviour often exceed interspecific differences. These patterns of temperature-dependent behaviour motivate further development of C. briggsae as a model system for dissecting the genetic underpinnings of complex behavioural traits.

Recombination modulates how selection affects linked sites in Drosophila

Recombination rate in Drosophila species shapes the impact of selection in the genome and is positively correlated with nucleotide diversity.

One of the most influential observations in molecular evolution has been a strong association between local recombination rate and nucleotide polymorphisms across the genome. This is interpreted as evidence for ubiquitous natural selection. The alternative explanation, that recombination is mutagenic, has been rejected by the absence of a similar association between local recombination rate and nucleotide divergence between species. However, many recent studies show that recombination rates are often very different even in closely related species, questioning whether an association between recombination rate and divergence between species has been tested satisfactorily. To circumvent this problem, we directly surveyed recombination across approximately 43% of the D. pseudoobscura physical genome in two separate recombination maps and 31% of the D. miranda physical genome, and we identified both global and local differences in recombination rate between these two closely related species. Using only regions with conserved recombination rates between and within species and accounting for multiple covariates, our data support the conclusion that recombination is positively related to diversity because recombination modulates Hill–Robertson effects in the genome and not because recombination is predominately mutagenic. Finally, we find evidence for dips in diversity around nonsynonymous substitutions. We infer that at least some of this reduction in diversity resulted from selective sweeps and examine these dips in the context of recombination rate.

Individuals within a species differ in the DNA sequences of their genes. This sequence variation affects how well individuals survive or reproduce and is transmitted to their offspring. Genes near each other on individual chromosomes tend to be passed to offspring together—neighboring genes are unlikely to be separated by exchanges of genetic material derived from different parents during meiotic recombination. When genes are inherited together, however, the evolutionary forces acting on one gene can interfere with variation at its neighbors. Thus, variation at multiple genes can be lost if natural selection acts on one gene in close proximity. Recombination can prevent or reduce this loss of variation, but previous tests of this phenomenon failed to account for recombination rate differences between species. In this study, we show that some parts of the genome differ in recombination rate between two species of fruit fly, Drosophila pseudoobscura and D. miranda. Avoiding an assumption made in previous studies, we then examine sequence variation within and between fly species in those parts of the genome that have conserved recombination rates. Based on the results, we conclude that recombination indeed preserves variation within species that would otherwise have been eliminated by natural selection.

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.

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.

Birth, death and subfunctionalization in the Arabidopsis genome

Arabidopsis thaliana is now a model system, not just for plant biology but also for comparative genomics. The completion of the sequences of two closely related species, Arabidopsis lyrata and Brassica rapa, is complemented by genomic comparisons among A. thaliana accessions and mutation accumulation lines. Together these genomic data document the birth of new genes via gene duplication, transposon exaptation and de novo formation of new genes from noncoding sequence. Most novel loci exhibit low expression, and are undergoing pseudogenization or subfunctionalization. Comparatively, A. thaliana has lost large amounts of sequence through deletion, particularly of transposable elements. Intraspecific genomic variation indicates high rates of deletion mutations and deletion polymorphisms across accessions, shedding light on the history of Arabidopsis genome architecture.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

Nematode model systems in evolution and development

The free-living nematode Caenorhabditis elegans is one of the most important model organisms in all areas of modern biology. Using the knowledge about C. elegans as a baseline, nematodes are now intensively studied in evolution and development. Evolutionary developmental biology or for short, 'evo-devo' has been developed as a new research discipline during the last two decades to investigate how changes in developmental processes and mechanisms result in the modification of morphological structures and phenotypic novelty. In this article, we review the concepts that make nematode evo-devo a successful approach to evolutionary biology. We introduce selected model systems for nematode evo-devo and provide a detailed discussion of four selected case studies. The most striking finding of nematode evo-devo is the magnitude of developmental variation in the context of a conserved body plan. Detailed investigation of early embryogenesis, gonad formation, vulva development, and sex determination revealed that molecular mechanisms evolve rapidly, often in the context of a conserved body plan. These studies highlight the importance of developmental systems drift and neutrality in evolution.

The effects of deleterious mutations on evolution at linked sites

The process of evolution at a given site in the genome can be influenced by the action of selection at other sites, especially when these are closely linked to it. Such selection reduces the effective population size experienced by the site in question (the Hill-Robertson effect), reducing the level of variability and the efficacy of selection. In particular, deleterious variants are continually being produced by mutation and then eliminated by selection at sites throughout the genome. The resulting reduction in variability at linked neutral or nearly neutral sites can be predicted from the theory of background selection, which assumes that deleterious mutations have such large effects that their behavior in the population is effectively deterministic. More weakly selected mutations can accumulate by Muller's ratchet after a shutdown of recombination, as in an evolving Y chromosome. Many functionally significant sites are probably so weakly selected that Hill-Robertson interference undermines the effective strength of selection upon them, when recombination is rare or absent. This leads to large departures from deterministic equilibrium and smaller effects on linked neutral sites than under background selection or Muller's ratchet. Evidence is discussed that is consistent with the action of these processes in shaping genome-wide patterns of variation and evolution.

The genetic basis of evolution, adaptation and speciation

A primary question in biology concerns the genetic basis of the evolution of novel traits, often in response to environmental changes, and how this can subsequently cause species isolation. This topic was the focus of the symposium on the Genetics of Speciation and Evolution at the annual meeting of the Canadian Society for Ecology and Evolution, held in Banff in May 2011. The presentations revealed some of the rapid advances being made in understanding the genetic basis of adaptation and speciation, as well as the elegant interplay between an organism's genetic complement and the environment that organism experiences.