Properties of selected mutations and genotypic landscapes under Fisher's Geometric Model

The fitness landscape - the mapping between genotypes and fitness - determines properties of the process of adaptation. Several small genetic fitness landscapes have recently been built by selecting a handful of beneficial mutations and measuring fitness of all combinations of these mutations. Here we generate several testable predictions for the properties of these landscapes under Fisher's geometric model of adaptation (FGMA). When far from the fitness optimum, we analytically compute the fitness effect of beneficial mutations and their epistatic interactions. We show that epistasis may be negative or positive on average depending on the distance of the ancestral genotype to the optimum and whether mutations were independently selected or co-selected in an adaptive walk. Using simulations, we show that genetic landscapes built from FGMA are very close to an additive landscape when the ancestral strain is far from the optimum. However, when close to the optimum, a large diversity of landscape with substantial ruggedness and sign epistasis emerged. Strikingly, landscapes built from different realizations of stochastic adaptive walks in the same exact conditions were highly variable, suggesting that several realizations of small genetic landscapes are needed to gain information about the underlying architecture of the global adaptive landscape.Comment: 51 pages, 8 figure

Selective sweeps under dominance and inbreeding

A major research goal in evolutionary genetics is to uncover loci experiencing positive selection. One approach involves finding ‘selective sweeps’ patterns, which can either be ‘hard sweeps’ formed by de novo mutation, or ‘soft sweeps’ arising from recurrent mutation or existing standing variation. Existing theory generally assumes outcrossing populations, and it is unclear how dominance affects soft sweeps. We consider how arbitrary dominance and inbreeding via self-fertilization affect hard and soft sweep signatures. With increased self-fertilization, they are maintained over longer map distances due to reduced effective recombination and faster beneficial allele fixation times. Dominance can affect sweep patterns in outcrossers if the derived variant originates from either a single novel allele, or from recurrent mutation. These models highlight the challenges in distinguishing hard and soft sweeps, and propose methods to differentiate between scenarios

"Menage à trois": the presence/absence of thyme shapes the mutualistic interaction between the host plant [i]Medicago truncatula[/i] (Fabaceae) and its symbiotic bacterium [i]Sinorhizobium meliloti[/i]

UMR 1334 AGAP : Equipe DAVEM Diversité et adaptation de la vigne et des espèces méditerranéennesThe long-term maintenance of specialized mutualisms remains an evolutionary puzzle. Recent focus has been on factors governing the stability of these mutualisms, including sanctions by the host, partner choice, and coevolutionary constraint, that is, the genetic correlation (r(G)) between fitness of both partners. So far these studies have been typically carried out in a single environment. Here, we ask if the genetic correlation between fitness of the host plant Medicago truncatula (Fabaceae) and its bacterial symbiont Sinorhizobium meliloti is affected by the presence/absence of a monoterpene (carvacrol) leached into the soil by Thymus vulgaris-a common plant of the Mediterranean vegetation, often co-occuring with Medicago. We show that the presence of carvacrol in the soil dramatically affects fitness of the rhizobial partner and increases the magnitude of rG between plant and rhizobia fitness (r(G) = 0.02 +/- 0.05 vs. r(G) = 0.57 +/- 0.02). This finding emphasizes the importance of heterogeneity in the biotic environment for understanding the evolution of species interactions

Insertion and Deletion Processes in Recent Human History

Background: Although insertions and deletions (indels) account for a sizable portion of genetic changes within and among species, they have received little attention because they are difficult to type, are alignment dependent and their underlying mutational process is poorly understood. A fundamental question in this respect is whether insertions and deletions are governed by similar or different processes and, if so, what these differences are. Methodology/Principal Findings: We use published resequencing data from Seattle SNPs and NIEHS human polymorphism databases to construct a genomewide data set of short polymorphic insertions and deletions in the human genome (n = 6228). We contrast these patterns of polymorphism with insertions and deletions fixed in the same regions since the divergence of human and chimpanzee (n = 10546). The macaque genome is used to resolve all indels into insertions and deletions. We find that the ratio of deletions to insertions is greater within humans than between human and chimpanzee. Deletions segregate at lower frequency in humans, providing evidence for deletions being under stronger purifying selection than insertions. The insertion and deletion rates correlate with several genomic features and we find evidence that both insertions and deletions are associated with point mutations. Finally, we find no evidence for a direct effect of the local recombination rate on the insertion and deletion rate. Conclusions/Significance: Our data strongly suggest that deletions are more deleterious than insertions but that insertion

Using singleton densities to detect recent selection in Bos taurus

Many quantitative traits are subject to polygenic selection, where several genomic regions undergo small, simultaneous changes in allele frequency that collectively alter a phenotype. The widespread availability of genome data, along with novel statistical techniques, has made it easier to detect these changes. We apply one such method, the “Singleton Density Score” (SDS), to the Holstein breed of Bos taurus to detect recent selection (arising up to around 740 years ago). We identify several genes as candidates for targets of recent selection, including some relating to cell regulation, catabolic processes, neural‐cell adhesion and immunity. We do not find strong evidence that three traits that are important to humans—milk protein content, milk fat content, and stature—have been subject to directional selection. Simulations demonstrate that because B. taurus recently experienced a population bottleneck, singletons are depleted so the power of SDS methods is reduced. These results inform on which genes underlie recent genetic change in B. taurus, while providing information on how polygenic selection can be best investigated in future studies

Adaptive evolution of the symbiotic gene NORK is not correlated with shifts of rhizobial specificity in the genus Medicago

<p>Abstract</p> <p>Background</p> <p>The <it>NODULATION RECEPTOR KINASE </it>(<it>NORK</it>) gene encodes a Leucine-Rich Repeat (LRR)-containing receptor-like protein and controls the infection by symbiotic rhizobia and endomycorrhizal fungi in Legumes. The occurrence of numerous amino acid changes driven by directional selection has been reported in this gene, using a limited number of messenger RNA sequences, but the functional reason of these changes remains obscure. The <it>Medicago </it>genus, where changes in rhizobial associations have been previously examined, is a good model to test whether the evolution of <it>NORK </it>is influenced by rhizobial interactions.</p> <p>Results</p> <p>We sequenced a region of 3610 nucleotides (encoding a 392 amino acid-long region of the NORK protein) in 32 <it>Medicago </it>species. We confirm that positive selection in <it>NORK </it>has occurred within the <it>Medicago </it>genus and find that the amino acid positions targeted by selection occur in sites outside of solvent-exposed regions in LRRs, and other sites in the N-terminal region of the protein. We tested if branches of the <it>Medicago </it>phylogeny where changes of rhizobial symbionts occurred displayed accelerated rates of amino acid substitutions. Only one branch out of five tested, leading to <it>M. noeana</it>, displays such a pattern. Among other branches, the most likely for having undergone positive selection is not associated with documented shift of rhizobial specificity.</p> <p>Conclusion</p> <p>Adaptive changes in the sequence of the NORK receptor have involved the LRRs, but targeted different sites than in most previous studies of LRR proteins evolution. The fact that positive selection in <it>NORK </it>tends not to be associated to changes in rhizobial specificity indicates that this gene was probably not involved in evolving rhizobial preferences. Other explanations (<it>e.g</it>. coevolutionary arms race) must be tested to explain the adaptive evolution of <it>NORK</it>.</p

GC-Content evolution in bacterial genomes: The biased gene conversion hypothesis expands

The characterization of functional elements in genomes relies on the identification of the footprints of natural selection. In this quest, taking into account neutral evolutionary processes such as mutation and genetic drift is crucial because these forces can generate patterns that may obscure or mimic signatures of selection. In mammals, and probably in many eukaryotes, another such confounding factor called GC-Biased Gene Conversion (gBGC) has been documented. This mechanism generates patterns identical to what is expected under selection for higher GC-content, specifically in highly recombining genomic regions. Recent results have suggested that a mysterious selective force favouring higher GC-content exists in Bacteria but the possibility that it could be gBGC has been excluded. Here, we show that gBGC is probably at work in most if not all bacterial species. First we find a consistent positive relationship between the GC-content of a gene and evidence of intra-genic recombination throughout a broad spectrum of bacterial clades. Second, we show that the evolutionary force responsible for this pattern is acting independently from selection on codon usage, and could potentially interfere with selection in favor of optimal AU-ending codons. A comparison with data from human populations shows that the intensity of gBGC in Bacteria is comparable to what has been reported in mammals. We propose that gBGC is not restricted to sexual Eukaryotes but also widespread among Bacteria and could therefore be an ancestral feature of cellular organisms. We argue that if gBGC occurs in bacteria, it can account for previously unexplained observations, such as the apparent non-equilibrium of base substitution patterns and the heterogeneity of gene composition within bacterial genomes. Because gBGC produces patterns similar to positive selection, it is essential to take this process into account when studying the evolutionary forces at work in bacterial genomes

Microsatellite diversity and broad scale geographic structure in a model legume: building a set of nested core collection for studying naturally occurring variation in Medicago truncatula

BACKGROUND: Exploiting genetic diversity requires previous knowledge of the extent and structure of the variation occurring in a species. Such knowledge can in turn be used to build a core-collection, i.e. a subset of accessions that aim at representing the genetic diversity of this species with a minimum of repetitiveness. We investigate the patterns of genetic diversity and population structure in a collection of 346 inbred lines representing the breadth of naturally occurring diversity in the Legume plant model Medicago truncatula using 13 microsatellite loci distributed throughout the genome. RESULTS: We confirm the uniqueness of all these genotypes and reveal a large amount of genetic diversity and allelic variation within this autogamous species. Spatial genetic correlation was found only for individuals originating from the same population and between neighbouring populations. Using a model-based clustering algorithm, we identified four main genetic clusters in the set of individuals analyzed. This stratification matches broad geographic regions. We also identified a set of "admixed" individuals that do not fit with this population structure scheme. CONCLUSION: The stratification inferred is discussed considering potential historical events like expansion, refuge history and admixture between neighbouring groups. Information on the allelic richness and the inferred population structure are used to build a nested core-collection. The set of inbred lines and the core collections are publicly available and will help coordinating efforts for the study of naturally occurring variation in the growing Medicago truncatula community

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions

Purifying Selection and Molecular Adaptation in the Genome of Verminephrobacter, the Heritable Symbiotic Bacteria of Earthworms

While genomic erosion is common among intracellular symbionts, patterns of genome evolution in heritable extracellular endosymbionts remain elusive. We study vertically transmitted extracellular endosymbionts (Verminephrobacter, Betaproteobacteria) that form a beneficial, species-specific, and evolutionarily old (60–130 Myr) association with earthworms. We assembled a draft genome of Verminephrobacter aporrectodeae and compared it with the genomes of Verminephrobacter eiseniae and two nonsymbiotic close relatives (Acidovorax). Similar to V. eiseniae, the V. aporrectodeae genome was not markedly reduced in size and showed no A–T bias. We characterized the strength of purifying selection (ω = dN/dS) and codon usage bias in 876 orthologous genes. Symbiont genomes exhibited strong purifying selection (ω = 0.09 ± 0.07), although transition to symbiosis entailed relaxation of purifying selection as evidenced by 50% higher ω values and less codon usage bias in symbiont compared with reference genomes. Relaxation was not evenly distributed among functional gene categories but was overrepresented in genes involved in signal transduction and cell envelope biogenesis. The same gene categories also harbored instances of positive selection in the Verminephrobacter clade. In total, positive selection was detected in 89 genes, including also genes involved in DNA metabolism, tRNA modification, and TonB-dependent iron uptake, potentially highlighting functions important in symbiosis. Our results suggest that the transition to symbiosis was accompanied by molecular adaptation, while purifying selection was only moderately relaxed, despite the evolutionary age and stability of the host association. We hypothesize that biparental transmission of symbionts and rare genetic mixing during transmission can prevent genome erosion in heritable symbionts