Predicting evolution and visualizing high-dimensional fitness landscapes

The tempo and mode of an adaptive process is strongly determined by the structure of the fitness landscape that underlies it. In order to be able to predict evolutionary outcomes (even on the short term), we must know more about the nature of realistic fitness landscapes than we do today. For example, in order to know whether evolution is predominantly taking paths that move upwards in fitness and along neutral ridges, or else entails a significant number of valley crossings, we need to be able to visualize these landscapes: we must determine whether there are peaks in the landscape, where these peaks are located with respect to one another, and whether evolutionary paths can connect them. This is a difficult task because genetic fitness landscapes (as opposed to those based on traits) are high-dimensional, and tools for visualizing such landscapes are lacking. In this contribution, we focus on the predictability of evolution on rugged genetic fitness landscapes, and determine that peaks in such landscapes are highly clustered: high peaks are predominantly close to other high peaks. As a consequence, the valleys separating such peaks are shallow and narrow, such that evolutionary trajectories towards the highest peak in the landscape can be achieved via a series of valley crossingsComment: 12 pages, 7 figures. To appear in "Recent Advances in the Theory and Application of Fitness Landscapes" (A. Engelbrecht and H. Richter, eds.). Springer Series in Emergence, Complexity, and Computation, 201

Prokaryote genome fluidity is dependent on effective population size

Many prokaryote species are known to have fluid genomes, with different strains varying markedly in accessory gene content through the combined action of gene loss, gene gain via lateral transfer, as well as gene duplication. However, the evolutionary forces determining genome fluidity are not yet well understood. We here for the first time systematically analyse the degree to which this distinctive genomic feature differs between bacterial species. We find that genome fluidity is positively correlated with synonymous nucleotide diversity of the core genome, a measure of effective population size Ne. No effects of genome size, phylogeny or homologous recombination rate on genome fluidity were found. Our findings are consistent with a scenario where accessory gene content turnover is for a large part dictated by neutral evolution

Bayesian paternity analysis and mating patterns in a parasitic nematode, Trichostrongylus tenuis

Mating behaviour is a fundamental aspect of the evolutionary ecology of sexually reproducing species, but one that has been under-researched in parasitic nematodes. We analysed mating behaviour in the parasitic nematode Trichostrongylus tenuis by performing a paternity analysis in a population from a single red grouse host. Paternity of the 150 larval offspring of 25 mothers (sampled from one of the two host caeca) was assigned among 294 candidate fathers (sampled from both caeca). Each candidate father's probability of paternity of each offspring was estimated from 10-locus microsatellite genotypes. Seventy-six (51%) offspring were assigned a father with a probability of >0.8, and the estimated number of unsampled males was 136 (95% credible interval (CI) 77-219). The probability of a male from one caecum fathering an offspring in the other caecum was estimated as 0.024 (95% CI 0.003-0.077), indicating that the junction of the caeca is a strong barrier to dispersal. Levels of promiscuity (defined as the probability of two of an adult's offspring sharing only one parent) were high for both sexes. Variance in male reproductive success was moderately high, possibly because of a combination of random mating and high variance in post-copulatory reproductive success. These results provide the first data on individual mating behaviour among parasitic nematodes

A barrier to homologous recombination between sympatric strains of the cooperative soil bacterium Myxococcus xanthus

The bacterium Myxococcus xanthus glides through soil in search of prey microbes, but when food sources run out, cells cooperatively construct and sporulate within multicellular fruiting bodies. M. xanthus strains isolated from a 16 × 16-cm-scale patch of soil were previously shown to have diversified into many distinct compatibility types that are distinguished by the failure of swarming colonies to merge upon encounter. We sequenced the genomes of 22 isolates from this population belonging to the two most frequently occurring multilocus sequence type (MLST) clades to trace patterns of incipient genomic divergence, specifically related to social divergence. Although homologous recombination occurs frequently within the two MLST clades, we find an almost complete absence of recombination events between them. As the two clades are very closely related and live in sympatry, either ecological or genetic barriers must reduce genetic exchange between them. We find that the rate of change in the accessory genome is greater than the rate of amino-acid substitution in the core genome. We identify a large genomic tract that consistently differs between isolates that do not freely merge and therefore is a candidate region for harbouring gene(s) responsible for self/non-self discrimination

Contrasting patterns of population structure and gene flow facilitate exploration of connectivity in two widely distributed temperate octocorals

This is the final version of the article. Available from Springer Nature via the DOI in this record.Connectivity is an important component of metapopulation dynamics in marine systems and can influence population persistence, migration rates and conservation decisions associated with Marine Protected Areas (MPAs). In this study, we compared the genetic diversity, gene flow and population structure of two octocoral species, Eunicella verrucosa and Alcyonium digitatum, in the northeast Atlantic (ranging from the northwest of Ireland and the southern North Sea, to southern Portugal), using two panels of thirteen and eight microsatellite loci, respectively. Our results identified regional genetic structure in E. verrucosa partitioned between populations from southern Portugal, northwest Ireland, and Britain/France; subsequent hierarchical analysis of population structure also indicated reduced gene flow between southwest Britain and northwest France. However, over a similar geographical area, A. digitatum showed little evidence of population structure, suggesting high gene flow and/or a large effective population size; indeed, the only significant genetic differentiation detected in A. digitatum occurred between North Sea samples and those from the English Channel/northeast Atlantic. In both species the vast majority of gene flow originated from sample sites within regions, with populations in southwest Britain being the predominant source of contemporary exogenous genetic variants for the populations studied. Unsurprisingly, historical patterns of gene flow appeared more complex, though again southwest Britain appeared an important source of genetic variation for both species. Our findings have major conservation implications, particularly for E. verrucosa, a protected species in UK waters and listed by the IUCN as ‘Vulnerable’, and for the designation and management of European MPAs.We thank Natural England (project No. RP0286, contract No. SAE 03-02-146), the NERC (grant No. NE/L002434/1) and the University of Exeter for funding this research. Additional funding for sample collection, travel and microsatellite development was provided by the EU Framework 7 ASSEMBLE programme, agreement no. 227799, and NERC grant No. NBAF-362

Effective polyploidy causes phenotypic delay and influences bacterial evolvability

Whether mutations in bacteria exhibit a noticeable delay before expressing their corresponding mutant phenotype was discussed intensively in the 1940s to 1950s, but the discussion eventually waned for lack of supportive evidence and perceived incompatibility with observed mutant distributions in fluctuation tests. Phenotypic delay in bacteria is widely assumed to be negligible, despite the lack of direct evidence. Here, we revisited the question using recombineering to introduce antibiotic resistance mutations into E. coli at defined time points and then tracking expression of the corresponding mutant phenotype over time. Contrary to previous assumptions, we found a substantial median phenotypic delay of three to four generations. We provided evidence that the primary source of this delay is multifork replication causing cells to be effectively polyploid, whereby wild-type gene copies transiently mask the phenotype of recessive mutant gene copies in the same cell. Using modeling and simulation methods, we explored the consequences of effective polyploidy for mutation rate estimation by fluctuation tests and sequencing-based methods. For recessive mutations, despite the substantial phenotypic delay, the per-copy or per-genome mutation rate is accurately estimated. However, the per-cell rate cannot be estimated by existing methods. Finally, with a mathematical model, we showed that effective polyploidy increases the frequency of costly recessive mutations in the standing genetic variation (SGV), and thus their potential contribution to evolutionary adaptation, while drastically reducing the chance that de novo recessive mutations can rescue populations facing a harsh environmental change such as antibiotic treatment. Overall, we have identified phenotypic delay and effective polyploidy as previously overlooked but essential components in bacterial evolvability, including antibiotic resistance evolution

Rapid and widespread de novo evolution of kin discrimination

Diverse forms of kin discrimination, broadly defined as alteration of social behavior as a function of genetic relatedness among interactants, are common among social organisms from microbes to humans. However, the evolutionary origins and causes of kin-discriminatory behavior remain largely obscure. One form of kin discrimination observed in microbes is the failure of genetically distinct colonies to merge freely upon encounter. Here, we first use natural isolates of the highly social bacterium Myxococcus xanthus to show that colony-merger incompatibilities can be strong barriers to social interaction, particularly by reducing chimerism in multicellular fruiting bodies that develop near colony-territory borders. We then use experimental laboratory populations to test hypotheses regarding the evolutionary origins of kin discrimination. We show that the generic process of adaptation, irrespective of selective environment, is sufficient to repeatedly generate kin-discriminatory behaviors between evolved populations and their common ancestor. Further, we find that kin discrimination pervasively evolves indirectly between allopatric replicate populations that adapt to the same ecological habitat and that this occurs generically in many distinct habitats. Patterns of interpopulation discrimination imply that kin discrimination phenotypes evolved via many diverse genetic mechanisms and mutation-accumulation patterns support this inference. Strong incompatibility phenotypes emerged abruptly in some populations but strengthened gradually in others. The indirect evolution of kin discrimination in an asexual microbe is analogous to the indirect evolution of reproductive incompatibility in sexual eukaryotes and linguistic incompatibility among human cultures, the commonality being indirect, noncoordinated divergence of complex systems evolving in isolation