No effect of recombination on the efficacy of natural selection in primates

Causes and consequences of a complex recombinational landscape in the ant Cardiocondyla obscurior

Eusocial Hymenoptera have the highest recombination rates among all multicellular animals studied so far, but it is unclear why this is and how this affects the biology of individual species. A high-resolution linkage map for the ant Cardiocondyla obscurior corroborates genome-wide high recombination rates reported for ants (8.1 cM/Mb). However, recombination is locally suppressed in regions that are enriched with TEs, that have strong haplotype divergence, or that show signatures of epistatic selection in C. obscurior The results do not support the hypotheses that high recombination rates are linked to phenotypic plasticity or to modulating selection efficiency. Instead, genetic diversity and the frequency of structural variants correlate positively with local recombination rates, potentially compensating for the low levels of genetic variation expected in haplodiploid social Hymenoptera with low effective population size. Ultimately, the data show that recombination contributes to within-population polymorphism and to the divergence of the lineages within C. obscurior.

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

The role of recombination dynamics in shaping signatures of direct and indirect selection across the Ficedula flycatcher genome†

Recombination is a central evolutionary process that reshuffles combinations of alleles along chromosomes, and consequently is expected to influence the efficacy of direct selection via Hill-Robertson interference. Additionally, the indirect effects of selection on neutral genetic diversity are expected to show a negative relationship with recombination rate, as background selection and genetic hitchhiking are stronger when recombination rate is low. However, owing to the limited availability of recombination rate estimates across divergent species, the impact of evolutionary changes in recombination rate on genomic signatures of selection remains largely unexplored. To address this question, we estimate recombination rate in two Ficedula flycatcher species, the taiga flycatcher (Ficedula albicilla) and collared flycatcher (Ficedula albicollis). We show that recombination rate is strongly correlated with signatures of indirect selection, and that evolutionary changes in recombination rate between species have observable impacts on this relationship. Conversely, signatures of direct selection on coding sequences show little to no relationship with recombination rate, even when restricted to genes where recombination rate is conserved between species. Thus, using measures of indirect and direct selection that bridge micro- and macro-evolutionary timescales, we demonstrate that the role of recombination rate and its dynamics varies for different signatures of selection.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree-based linkage mapping. This species has extended Z and W chromosomes ("neo-sex chromosomes") formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex-linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A-part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein-coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A-derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z-linked, recombination patterns.

Negative linkage disequilibrium between amino acid changing variants reveals interference among deleterious mutations in the human genome

Evolutionary forces like Hill-Robertson interference and negative epistasis can lead to deleterious mutations being found on distinct haplotypes. However, the extent to which these forces depend on the selection and dominance coefficients of deleterious mutations and shape genome-wide patterns of linkage disequilibrium (LD) in natural populations with complex demographic histories has not been tested. In this study, we first used forward-in-time simulations to predict how negative selection impacts LD. Under models where deleterious mutations have additive effects on fitness, deleterious variants less than 10 kb apart tend to be carried on different haplotypes relative to pairs of synonymous SNPs. In contrast, for recessive mutations, there is no consistent ordering of how selection coefficients affect LD decay, due to the complex interplay of different evolutionary effects. We then examined empirical data of modern humans from the 1000 Genomes Project. LD between derived alleles at nonsynonymous SNPs is lower compared to pairs of derived synonymous variants, suggesting that nonsynonymous derived alleles tend to occur on different haplotypes more than synonymous variants. This result holds when controlling for potential confounding factors by matching SNPs for frequency in the sample (allele count), physical distance, magnitude of background selection, and genetic distance between pairs of variants. Lastly, we introduce a new statistic HR(j) which allows us to detect interference using unphased genotypes. Application of this approach to high-coverage human genome sequences confirms our finding that nonsynonymous derived alleles tend to be located on different haplotypes more often than are synonymous derived alleles. Our findings suggest that interference may play a pervasive role in shaping patterns of LD between deleterious variants in the human genome, and consequently influences genome-wide patterns of LD.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Impact of Mutation Rate and Selection at Linked Sites on DNA Variation across the Genomes of Humans and Other Homininae

DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both noncoding and nonsynonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50 kb windows across the genomes of humans and nonhuman homininae. Interestingly, we find that although noncoding diversity is equally affected by these three genomic variables, nonsynonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggest that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of nonsynonymous to noncoding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Exploiting selection at linked sites to infer the rate and strength of adaptation

Genomic data encodes past evolutionary events and has the potential to reveal the strength, rate, and biological drivers of adaptation. However, jointly estimating adaptation rate (a) and adaptation strength remains challenging because evolutionary processes such as demography, linkage, and non-neutral polymorphism can confound inference. Here, we exploit the influence of background selection to reduce the fixation rate of weakly-beneficial alleles to jointly infer the strength and rate of adaptation. We develop an MK-based method (ABC-MK) to infer adaptation rate and strength, and estimate α = 0.135 in human protein-coding sequences, 72% of which is contributed by weakly-adaptive variants. We show that in this adaptation regime α is reduced ≈ 25% by linkage genome-wide. Moreover, we show that virus-interacting proteins (VIPs) undergo adaptation that is both stronger and nearly twice as frequent as the genome average (α = 0.224, 56% due to strongly-beneficial alleles). Our results suggest that while most adaptation in human proteins is weakly-beneficial, adaptation to viruses is often strongly-beneficial. Our method provides a robust framework for estimating adaptation rate and strength across species.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

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.

Genetic constraints on adaptation: a theoretical primer for the genomics era

Genetic constraints are features of inheritance systems that slow or prohibit adaptation. Several population genetic mechanisms of constraint have received sustained attention within the field since they were first articulated in the early 20th century. This attention is now reflected in a rich, and still growing, theoretical literature on the genetic limits to adaptive change. In turn, empirical research on constraints has seen a rapid expansion over the last two decades in response to changing interests of evolutionary biologists, along with new technologies, expanding data sets, and creative analytical approaches that blend mathematical modeling with genomics. Indeed, one of the most notable and exciting features of recent progress in genetic constraints is the close connection between theoretical and empirical research. In this review, we discuss five major population genetic contexts of genetic constraint: genetic dominance, pleiotropy, fitness trade-offs between types of individuals of a population, sign epistasis, and genetic linkage between loci. For each, we outline historical antecedents of the theory, specific contexts where constraints manifest, and their quantitative consequences for adaptation. From each of these theoretical foundations, we discuss recent empirical approaches for identifying and characterizing genetic constraints, each grounded and motivated by this theory, and outline promising areas for future work.

The Limits of Natural Selection in a Nonequilibrium World

Evolutionary theory predicts that factors such as a small population size or low recombination rate can limit the action of natural selection. The emerging field of comparative population genomics offers an opportunity to evaluate these hypotheses. However, classical theoretical predictions assume that populations are at demographic equilibrium. This assumption is likely to be violated in the very populations researchers use to evaluate selection's limits: populations that have experienced a recent shift in population size and/or effective recombination rates. Here we highlight theory and data analyses concerning limitations on the action of natural selection in nonequilibrial populations and argue that substantial care is needed to appropriately test whether species and populations show meaningful differences in selection efficacy. A move toward model-based inferences that explicitly incorporate nonequilibrium dynamics provides a promising approach to more accurately contrast selection efficacy across populations and interpret its significance.

Genetic Effect of Chemotherapy Exposure in Children of Testicular Cancer Survivors

PURPOSE

Cancer survivors express anxiety that chemotherapy exposure may lead to transmissible genetic damage in posttreatment children. Preclinical models suggest that chemotherapy exposure may result in considerable genomic alterations in postexposure progeny. Epidemiologic studies have not demonstrated a significant increase in congenital abnormalities in posttreatment children of cancer survivors, but the inherited genome-wide effect of chemotherapy exposure in humans is unknown.

EXPERIMENTAL DESIGN

Two testicular cancer survivors cured with chemotherapy who had children pre- and postexposure without sperm banking were identified. Familial germline whole genome sequencing (WGS) was performed for these families, and analytic methods were utilized to identify de novo alterations, including mutations, recombinations, and structural rearrangements in the pre- and postexposure offspring.

RESULTS

No increase in de novo germline mutations in postexposure children compared with their preexposure siblings was found. Furthermore, there were no increased short insertion/deletions, recombination frequency, or structural rearrangements in these postexposure children.

CONCLUSIONS

In two families of male cancer survivors, there was no transmissible genomic impact of significant mutagenic exposure in postexposure children. This study may provide possible reassuring evidence for patients undergoing chemotherapy who are unable to have pretreatment sperm cryopreservation. Expanded cohorts that utilize WGS to identify environmental exposure effects on the inherited genome may inform the generalizability of these results. Clin Cancer Res; 22(9); 2183-9. ©2015 AACR.

Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits

BACKGROUND

Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each).

RESULTS

For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination.

CONCLUSIONS

Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

Tissue-Specific Evolution of Protein Coding Genes in Human and Mouse

Protein-coding genes evolve at different rates, and the influence of different parameters, from gene size to expression level, has been extensively studied. While in yeast gene expression level is the major causal factor of gene evolutionary rate, the situation is more complex in animals. Here we investigate these relations further, especially taking in account gene expression in different organs as well as indirect correlations between parameters. We used RNA-seq data from two large datasets, covering 22 mouse tissues and 27 human tissues. Over all tissues, evolutionary rate only correlates weakly with levels and breadth of expression. The strongest explanatory factors of purifying selection are GC content, expression in many developmental stages, and expression in brain tissues. While the main component of evolutionary rate is purifying selection, we also find tissue-specific patterns for sites under neutral evolution and for positive selection. We observe fast evolution of genes expressed in testis, but also in other tissues, notably liver, which are explained by weak purifying selection rather than by positive selection.

Highly variable recombinational landscape modulates efficacy of natural selection in birds

Determining the rate of protein evolution and identifying the causes of its variation across the genome are powerful ways to understand forces that are important for genome evolution. By using a multitissue transcriptome data set from great tit (Parus major), we analyzed patterns of molecular evolution between two passerine birds, great tit and zebra finch (Taeniopygia guttata), using the chicken genome (Gallus gallus) as an outgroup. We investigated whether a special feature of avian genomes, the highly variable recombinational landscape, modulates the efficacy of natural selection through the effects of Hill-Robertson interference, which predicts that selection should be more effective in removing deleterious mutations and incorporating beneficial mutations in high-recombination regions than in low-recombination regions. In agreement with these predictions, genes located in low-recombination regions tend to have a high proportion of neutrally evolving sites and relaxed selective constraint on sites subject to purifying selection, whereas genes that show strong support for past episodes of positive selection appear disproportionally in high-recombination regions. There is also evidence that genes located in high-recombination regions tend to have higher gene expression specificity than those located in low-recombination regions. Furthermore, more compact genes (i.e., those with fewer/shorter introns or shorter proteins) evolve faster than less compact ones. In sum, our results demonstrate that transcriptome sequencing is a powerful method to answer fundamental questions about genome evolution in nonmodel organisms.

Population size and the rate of evolution

Does evolution proceed faster in larger or smaller populations? The relationship between effective population size (Ne) and the rate of evolution has consequences for our ability to understand and interpret genomic variation, and is central to many aspects of evolution and ecology. Many factors affect the relationship between Ne and the rate of evolution, and recent theoretical and empirical studies have shown some surprising and sometimes counterintuitive results. Some mechanisms tend to make the relationship positive, others negative, and they can act simultaneously. The relationship also depends on whether one is interested in the rate of neutral, adaptive, or deleterious evolution. Here, we synthesize theoretical and empirical approaches to understanding the relationship and highlight areas that remain poorly understood.

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR-retrotransposons, the rates of synonymous and non-synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non-synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter-specific asymmetric evolution.

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.

Recombination rates and genomic shuffling in human and chimpanzee--a new twist in the chromosomal speciation theory

A long-standing question in evolutionary biology concerns the effect of recombination in shaping the genomic architecture of organisms and, in particular, how this impacts the speciation process. Despite efforts employed in the last decade, the role of chromosomal reorganizations in the human-chimpanzee speciation process remains unresolved. Through whole-genome comparisons, we have analyzed the genome-wide impact of genomic shuffling in the distribution of human recombination rates during the human-chimpanzee speciation process. We have constructed a highly refined map of the reorganizations and evolutionary breakpoint regions in the human and chimpanzee genomes based on orthologous genes and genome sequence alignments. The analysis of the most recent human and chimpanzee recombination maps inferred from genome-wide single-nucleotide polymorphism data revealed that the standardized recombination rate was significantly lower in rearranged than in collinear chromosomes. In fact, rearranged chromosomes presented significantly lower recombination rates than chromosomes that have been maintained since the ancestor of great apes, and this was related with the lineage in which they become fixed. Importantly, inverted regions had lower recombination rates than collinear and noninverted regions, independently of the effect of centromeres. Our observations have implications for the chromosomal speciation theory, providing new evidences for the contribution of inversions in suppressing recombination in mammals.

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.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Genome evolution in outcrossing versus selfing versus asexual species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes, and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. However, while asexual species essentially conform to theoretical predictions, the situation is less simple in selfing species. We discuss the possible reasons to potentially explain this paradox. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Direct and indirect consequences of meiotic recombination: implications for genome evolution

There is considerable variation within eukaryotic genomes in the local rate of crossing over. Why is this and what effect does it have on genome evolution? On the genome scale, it is known that by shuffling alleles, recombination increases the efficacy of selection. By contrast, the extent to which differences in the recombination rate modulate the efficacy of selection between genomic regions is unclear. Recombination also has direct consequences on the origin and fate of mutations: biased gene conversion and other forms of meiotic drive promote the fixation of mutations in a similar way to selection, and recombination itself may be mutagenic. Consideration of both the direct and indirect effects of recombination is necessary to understand why its rate is so variable and for correct interpretation of patterns of genome evolution.

Quantifying the variation in the effective population size within a genome

The effective population size (N(e)) is one of the most fundamental parameters in population genetics. It is thought to vary across the genome as a consequence of differences in the rate of recombination and the density of selected sites due to the processes of genetic hitchhiking and background selection. Although it is known that there is intragenomic variation in the effective population size in some species, it is not known whether this is widespread or how much variation in the effective population size there is. Here, we test whether the effective population size varies across the genome, between protein-coding genes, in 10 eukaryotic species by considering whether there is significant variation in neutral diversity, taking into account differences in the mutation rate between loci by using the divergence between species. In most species we find significant evidence of variation. We investigate whether the variation in N(e) is correlated to recombination rate and the density of selected sites in four species, for which these data are available. We find that N(e) is positively correlated to recombination rate in one species, Drosophila melanogaster, and negatively correlated to a measure of the density of selected sites in two others, humans and Arabidopsis thaliana. However, much of the variation remains unexplained. We use a hierarchical Bayesian analysis to quantify the amount of variation in the effective population size and show that it is quite modest in all species-most genes have an N(e) that is within a few fold of all other genes. Nonetheless we show that this modest variation in N(e) is sufficient to cause significant differences in the efficiency of natural selection across the genome, by demonstrating that the ratio of the number of nonsynonymous to synonymous polymorphisms is significantly correlated to synonymous diversity and estimates of N(e), even taking into account the obvious nonindependence between these measures.

Genomic determinants of protein evolution and polymorphism in Arabidopsis

Recent results from Drosophila suggest that positive selection has a substantial impact on genomic patterns of polymorphism and divergence. However, species with smaller population sizes and/or stronger population structure may not be expected to exhibit Drosophila-like patterns of sequence variation. We test this prediction and identify determinants of levels of polymorphism and rates of protein evolution using genomic data from Arabidopsis thaliana and the recently sequenced Arabidopsis lyrata genome. We find that, in contrast to Drosophila, there is no negative relationship between nonsynonymous divergence and silent polymorphism at any spatial scale examined. Instead, synonymous divergence is a major predictor of silent polymorphism, which suggests variation in mutation rate as the main determinant of silent variation. Variation in rates of protein divergence is mainly correlated with gene expression level and breadth, consistent with results for a broad range of taxa, and map-based estimates of recombination rate are only weakly correlated with nonsynonymous divergence. Variation in mutation rates and the strength of purifying selection seem to be major drivers of patterns of polymorphism and divergence in Arabidopsis. Nevertheless, a model allowing for varying negative and positive selection by functional gene category explains the data better than a homogeneous model, implying the action of positive selection on a subset of genes. Genes involved in disease resistance and abiotic stress display high proportions of adaptive substitution. Our results are important for a general understanding of the determinants of rates of protein evolution and the impact of selection on patterns of polymorphism and divergence.

Effective population size and the efficacy of selection on the X chromosomes of two closely related Drosophila species

The prevalence of natural selection relative to genetic drift is of central interest in evolutionary biology. Depending on the distribution of fitness effects of new mutations, the importance of these evolutionary forces may differ in species with different effective population sizes. Here, we survey population genetic variation at 105 orthologous X-linked protein coding regions in Drosophila melanogaster and its sister species D. simulans, two closely related species with distinct demographic histories. We observe significantly higher levels of polymorphism and evidence for stronger selection on codon usage bias in D. simulans, consistent with a larger historical effective population size on average for this species. Despite these differences, we estimate that <10% of newly arising nonsynonymous mutations have deleterious fitness effects in the nearly neutral range (i.e., −10 < N_es < 0) in both species. The inferred distributions of fitness effects and demographic models translate into surprisingly high estimates of the fraction of “adaptive” protein divergence in both species (∼85–90%). Despite evidence for different demographic histories, differences in population size have apparently played little role in the dynamics of protein evolution in these two species, and estimates of the adaptive fraction (α) of protein divergence in both species remain high even if we account for recent 10-fold growth. Furthermore, although several recent studies have noted strong signatures of recurrent adaptive protein evolution at genes involved in immunity, reproduction, sexual conflict, and intragenomic conflict, our finding of high levels of adaptive protein divergence at randomly chosen proteins (with respect to function) suggests that many other factors likely contribute to the adaptive protein divergence signature in Drosophila.

Detecting positive selection within genomes: the problem of biased gene conversion

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio (d(N)/d(S)), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d(N)/d(S) (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.

Contrasting methods of quantifying fine structure of human recombination

There has been considerable excitement over the ability to construct linkage maps based only on genome-wide genotype data for single nucleotide polymorphic sites (SNPs) in a population sample. These maps, which are derived from estimates of linkage disequilibrium (LD), rely on population genetics theory to relate the decay of LD to the local rate of recombination, but other population processes also come into play. Here we contrast these LD maps to the classically derived, pedigree-based human recombination maps. The LD maps have a level of resolution greatly exceeding that of the pedigree maps, and at this fine scale, sperm typing allows a means of validation. While at a gross level both the pedigree maps and the sperm typing methods generally agree with LD maps, there are significant local differences between them, and the fact that these maps measure different genetic features should be remembered when using them for other genetic inferences.

Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.

Microsatellite Diversity in Natural Populations of Ascomycetous Fungus, Emericella Nidulans, from Different Climatic-Edaphic Conditions in Israel

The genetic divergence of Israeli populations of the soil ascomycetous fungusEmericella nidulanswas studied on regional and local scales using fifteen microsatellite (SSR) markers. The study was performed in the framework of the "Evolution Canyon" research program at the Institute of Evolution, University of Haifa, in three "Evolution Canyons" (ECs): EC I (Mt. Carmel), EC II (western Upper Galilee), and EC III (the southern Negev desert). The first two canyons (EC I and EC II) are located in the northern part of Israel at a distance of 38 km apart; EC III is located southward at a distance of nearly 350 km from the northern ECs. In each canyon,E. nidulansstrains were isolated from opposite slopes and, in EC III, from the valley bottom. All three EC populations ofE. nidulanswere found to be genetically distinct. The estimated genetic divergences correspond to geographical distances and ecological differences between the three studied canyons. On a regional scale, simple sequence repeat (SSR) polymorphism tends to increase with severity of ecological conditions. In general, both environmental parameters (soil moisture and temperature) and genetic factors (predicted number of repeats in SSR markers, distance from marker to centromere, codon evolutionary chronologies, and hydrophobic vs. hydrophilic character of encoded amino acid) influenced genetic diversity ofE. nidulanspopulations.

Comparative genomics based on massive parallel transcriptome sequencing reveals patterns of substitution and selection across 10 bird species

Next-generation sequencing technology provides an attractive means to obtain large-scale sequence data necessary for comparative genomic analysis. To analyse the patterns of mutation rate variation and selection intensity across the avian genome, we performed brain transcriptome sequencing using Roche 454 technology of 10 different non-model avian species. Contigs from de novo assemblies were aligned to the two available avian reference genomes, chicken and zebra finch. In total, we identified 6499 different genes across all 10 species, with approximately 1000 genes found in each full run per species. We found evidence for a higher mutation rate of the Z chromosome than of autosomes (male-biased mutation) and a negative correlation between the neutral substitution rate (d(S)) and chromosome size. Analyses of the mean d(N)/d(S) ratio (omega) of genes across chromosomes supported the Hill-Robertson effect (the effect of selection at linked loci) and point at stochastic problems with omega as an independent measure of selection. Overall, this study demonstrates the usefulness of next-generation sequencing for obtaining genomic resources for comparative genomic analysis of non-model organisms.

Mutation and the evolution of recombination

Under the classical view, selection depends more or less directly on mutation: standing genetic variance is maintained by a balance between selection and mutation, and adaptation is fuelled by new favourable mutations. Recombination is favoured if it breaks negative associations among selected alleles, which interfere with adaptation. Such associations may be generated by negative epistasis, or by random drift (leading to the Hill-Robertson effect). Both deterministic and stochastic explanations depend primarily on the genomic mutation rate, U. This may be large enough to explain high recombination rates in some organisms, but seems unlikely to be so in general. Random drift is a more general source of negative linkage disequilibria, and can cause selection for recombination even in large populations, through the chance loss of new favourable mutations. The rate of species-wide substitutions is much too low to drive this mechanism, but local fluctuations in selection, combined with gene flow, may suffice. These arguments are illustrated by comparing the interaction between good and bad mutations at unlinked loci under the infinitesimal model.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

Protein rates of evolution are predicted by double-strand break events, independent of crossing-over rates

Theory predicts that, owing to reduced Hill–Robertson interference, genomic regions with high crossing-over rates should experience more efficient selection. In Saccharomyces cerevisiae a negative correlation between the local recombination rate, assayed as meiotic double-strand breaks (DSBs), and the local rate of protein evolution has been considered consistent with such a model. Although DSBs are a prerequisite for crossing-over, they need not result in crossing-over. With recent high-resolution crossover data, we now return to this issue comparing two species of yeast. Strikingly, even allowing for crossover rates, both the rate of premeiotic DSBs and of noncrossover recombination events predict a gene's rate of evolution. This both questions the validity of prior analyses and strongly suggests that any correlation between crossover rates and rates of protein evolution could be owing to slow-evolving genes being prone to DSBs or a direct effect of DSBs on sequence evolution. To ask if classical theory of recombination has any relevance, we determine whether crossover rates predict rates of protein evolution, controlling for noncrossover DSB events, gene ontology (GO) class, gene expression, protein abundance, nucleotide content, and dispensability. We find that genes with high crossing-over rates have low rates of protein evolution after such control, although any correlation is weaker than that previously reported considering meiotic DSBs as a proxy. The data are consistent both with recombination enhancing the efficiency of purifying selection and, independently, with DSBs being associated with low rates of evolution.

Complex evolution of a Y-chromosomal double homeobox 4 (DUX4)-related gene family in hominoids

The human Y chromosome carries four human Y-chromosomal euchromatin/heterochromatin transition regions, all of which are characterized by the presence of interchromosomal segmental duplications. The Yq11.1/Yq11.21 transition region harbours a peculiar segment composed of an imperfectly organized tandem-repeat structure encoding four members of the double homeobox (DUX) gene family. By comparative fluorescence in situ hybridization (FISH) analysis we have documented the primary appearance of Y-chromosomal DUX genes (DUXY) on the gibbon Y chromosome. The major amplification and dispersal of DUXY paralogs occurred after the gibbon and hominid lineages had diverged. Orthologous DUXY loci of human and chimpanzee show a highly similar structural organization. Sequence alignment survey, phylogenetic reconstruction and recombination detection analyses of human and chimpanzee DUXY genes revealed the existence of all copies in a common ancestor. Comparative analysis of the circumjacent beta-satellites indicated that DUXY genes and beta-satellites evolved in concert. However, evolutionary forces acting on DUXY genes may have induced amino acid sequence differences in the orthologous chimpanzee and human DUXY open reading frames (ORFs). The acquisition of complete ORFs in human copies might relate to evolutionary advantageous functions indicating neo-functionalization. We propose an evolutionary scenario in which an ancestral tandem array DUX gene cassette transposed to the hominoid Y chromosome followed by lineage-specific chromosomal rearrangements paved the way for a species-specific evolution of the Y-chromosomal members of a large highly diverged homeobox gene family.

Pervasive hitchhiking at coding and regulatory sites in humans

Much effort and interest have focused on assessing the importance of natural selection, particularly positive natural selection, in shaping the human genome. Although scans for positive selection have identified candidate loci that may be associated with positive selection in humans, such scans do not indicate whether adaptation is frequent in general in humans. Studies based on the reasoning of the MacDonald-Kreitman test, which, in principle, can be used to evaluate the extent of positive selection, suggested that adaptation is detectable in the human genome but that it is less common than in Drosophila or Escherichia coli. Both positive and purifying natural selection at functional sites should affect levels and patterns of polymorphism at linked nonfunctional sites. Here, we search for these effects by analyzing patterns of neutral polymorphism in humans in relation to the rates of recombination, functional density, and functional divergence with chimpanzees. We find that the levels of neutral polymorphism are lower in the regions of lower recombination and in the regions of higher functional density or divergence. These correlations persist after controlling for the variation in GC content, density of simple repeats, selective constraint, mutation rate, and depth of sequencing coverage. We argue that these results are most plausibly explained by the effects of natural selection at functional sites -- either recurrent selective sweeps or background selection -- on the levels of linked neutral polymorphism. Natural selection at both coding and regulatory sites appears to affect linked neutral polymorphism, reducing neutral polymorphism by 6% genome-wide and by 11% in the gene-rich half of the human genome. These findings suggest that the effects of natural selection at linked sites cannot be ignored in the study of neutral human polymorphism.

An empirical test for branch-specific positive selection

The use of phylogenetic analysis to predict positive selection specific to human genes is complicated by the very close evolutionary relationship with our nearest extant primate relatives, chimpanzees. To assess the power and limitations inherent in use of maximum-likelihood (ML) analysis of codon substitution patterns in such recently diverged species, a series of simulations was performed to assess the impact of several parameters of the evolutionary model on prediction of human-specific positive selection, including branch length and d(N)/d(S) ratio. Parameters were varied across a range of values observed in alignments of 175 transcription factor (TF) genes that were sequenced in 12 primate species. The ML method largely lacks the power to detect positive selection that has occurred since the most recent common ancestor between humans and chimpanzees. An alternative null model was developed on the basis of gene-specific evaluation of the empirical distribution of ML results, using simulated neutrally evolving sequences. This empirical test provides greater sensitivity to detect lineage-specific positive selection in the context of recent evolutionary divergence.