Temporal challenges in detecting balancing selection from population genomic data

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Reproductive strategies and their consequences for divergence, gene flow, and genetic diversity in three taxa of Clarkia

Differences in reproductive strategies can have important implications for macro- and micro-evolutionary processes. We used a comparative approach through a population genetics lens to evaluate how three distinct reproductive strategies shape patterns of divergence among as well as gene flow and genetic diversity within three closely related taxa in the genus Clarkia. One taxon is a predominantly autonomous self-fertilizer and the other two taxa are predominantly outcrossing but vary in the primary pollinator they attract. In genotyping populations using genotyping-by-sequencing and comparing loci shared across taxa, our results suggest that differences in reproductive strategies in part promote evolutionary divergence among these closely related taxa. Contrary to expectations, we found that the selfing taxon had the highest levels of heterozygosity but a low rate of polymorphism. The high levels of fixed heterozygosity for a subset of loci suggests this pattern is driven by the presence of structural rearrangements in chromosomes common in other Clarkia taxa. In evaluating patterns within taxa, we found a complex interplay between reproductive strategy and geographic distribution. Differences in the mobility of primary pollinators did not translate to a difference in rates of genetic diversity and gene flow within taxa - a pattern likely due to one taxon having a patchier distribution and a less temporally and spatially reliable pollinator. Taken together, this work advances our understanding of the factors that shape gene flow and the distribution of genetic diversity within and among closely related taxa.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

The effects of inversion polymorphisms on patterns of neutral genetic diversity

The strong reduction in the frequency of recombination in heterozygotes for an inversion and a standard gene arrangement causes the arrangements to become partially isolated genetically, resulting in sequence divergence between them and changes in the levels of neutral variability at nucleotide sites within each arrangement class. Previous theoretical studies on the effects of inversions on neutral variability have assumed either that the population is panmictic or that it is divided into 2 populations subject to divergent selection. Here, the theory is extended to a model of an arbitrary number of demes connected by migration, using a finite island model with the inversion present at the same frequency in all demes. Recursion relations for mean pairwise coalescent times are used to obtain simple approximate expressions for diversity and divergence statistics for an inversion polymorphism at equilibrium under recombination and drift, and for the approach to equilibrium following the sweep of an inversion to a stable intermediate frequency. The effects of an inversion polymorphism on patterns of linkage disequilibrium are also examined. The reduction in effective recombination rate caused by population subdivision can have significant effects on these statistics. The theoretical results are discussed in relation to population genomic data on inversion polymorphisms, with an emphasis on Drosophila melanogaster. Methods are proposed for testing whether or not inversions are close to recombination-drift equilibrium, and for estimating the rate of recombinational exchange in heterozygotes for inversions; difficulties involved in estimating the ages of inversions are also discussed.

Host body size, not host population size, predicts genome-wide effective population size of parasites

The effective population size (N_e) of an organism is expected to be generally proportional to the total number of individuals in a population. In parasites, we might expect the effective population size to be proportional to host population size and host body size, because both are expected to increase the number of parasite individuals. However, among other factors, parasite populations are sometimes so extremely subdivided that high levels of inbreeding may distort these predicted relationships. Here, we used whole-genome sequence data from dove parasites (71 feather louse species of the genus Columbicola) and phylogenetic comparative methods to study the relationship between parasite effective population size and host population size and body size. We found that parasite effective population size is largely explained by host body size but not host population size. These results suggest the potential local population size (infrapopulation or deme size) is more predictive of the long-term effective population size of parasites than is the total number of potential parasite infrapopulations (i.e., host individuals).

A spatial approach to jointly estimate Wright's neighborhood size and long-term effective population size

Spatially continuous patterns of genetic differentiation, which are common in nature, are often poorly described by existing population genetic theory or methods that assume panmixia or discrete, clearly definable populations. There is therefore a need for statistical approaches in population genetics that can accommodate continuous geographic structure, and that ideally use georeferenced individuals as the unit of analysis, rather than populations or subpopulations. In addition, researchers are often interested describing the diversity of a population distributed continuously in space, and this diversity is intimately linked to the dispersal potential of the organism. A statistical model that leverages information from patterns of isolation-by-distance to jointly infer parameters that control local demography (such as Wright's neighborhood size), and the long-term effective size (Ne) of a population would be useful. Here, we introduce such a model that uses individual-level pairwise genetic and geographic distances to infer Wright's neighborhood size and long-term Ne. We demonstrate the utility of our model by applying it to complex, forward-time demographic simulations as well as an empirical dataset of the Red Sea clownfish (Amphiprion bicinctus). The model performed well on simulated data relative to alternative approaches and produced reasonable empirical results given the natural history of clownfish. The resulting inferences provide important insights into the population genetic dynamics of spatially structure populations.

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

Host association, environment, and geography underlie genomic differentiation in a major forest pest

Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high‐throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance‐based redundancy analyses (dbRDA), we show that differentiation of thousands of genome‐wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host‐associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.

Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps

Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.

A Change in Conservation Status of Pachyphytum caesium (Crassulaceae), a Threatened Species from Central Mexico Based on Genetic Studies

Simple Summary

Species decline has accelerated in recent decades, with rare species often being the first to go extinct, mainly due to low population sizes. This study worked with Pachyphytum cesium (Crassulaceae), an endemic species from central Mexico threatened by human activities and climate change. To increase our knowledge of the species, we analyze the diversity and genetic structure of all known populations of P. caesium to change their current genetic status and propose conservation strategies. The results indicate that this species presents low to moderate levels of genetic diversity and incipiently structured populations due to low genetic flow. We found that genetic parameters are essential to the conservation strategies and status vulnerability reclassification. Based on our results, we propose reclassifying the conservation status of the species in danger of extinction; hence a series of conservation strategies are provided to prevent its disappearance.

Abstract

Genetic studies in rare species are important to determine their patterns of genetic diversity among populations and implement conservation plans aimed to reduce extinction risks. Pachyphytum caesium is an endemic species with extremely small populations in central Mexico. This work analyzes the diversity and genetic structure of Pachyphytum cesium (Crassulaceae) to change the conservation status and propose conservation strategies. Twelve dominant ISSR markers were used to describe the genetic diversity of all known populations. Additionally, we carried out two simulations to explore how the loss of individuals or the local populations extinction affect the genetics parameters of the species. The average results indicated moderate levels of genetic diversity (number of alleles = 89.7 ± 3.9, expected heterozygosity = 0.212 ± 0.0, and percentage of polymorphic loci = 56.1 ± 3.9), parameters that decreased significantly when simulations were performed in P. caesium. Additionally, a genetic structure of the populations was found with low gene flow (Nm = 0.92). Genetic parameters are negatively affected by the loss of individuals and the local extinction of populations. Based on our results, we propose to reclassify the conservation status of the species in danger of extinction, hence a series of conservation strategies are provided to prevent its disappearance.

Genealogical structure changes as range expansions transition from pushed to pulled

Range expansions accelerate evolution through multiple mechanisms, including gene surfing and genetic drift. The inference and control of these evolutionary processes ultimately rely on the information contained in genealogical trees. Currently, there are two opposing views on how range expansions shape genealogies. In invasion biology, expansions are typically approximated by a series of population bottlenecks producing genealogies with only pairwise mergers between lineages-a process known as the Kingman coalescent. Conversely, traveling wave models predict a coalescent with multiple mergers, known as the Bolthausen-Sznitman coalescent. Here, we unify these two approaches and show that expansions can generate an entire spectrum of coalescent topologies. Specifically, we show that tree topology is controlled by growth dynamics at the front and exhibits large differences between pulled and pushed expansions. These differences are explained by the fluctuations in the total number of descendants left by the early founders. High growth cooperativity leads to a narrow distribution of reproductive values and the Kingman coalescent. Conversely, low growth cooperativity results in a broad distribution, whose exponent controls the merger sizes in the genealogies. These broad distribution and non-Kingman tree topologies emerge due to the fluctuations in the front shape and position and do not occur in quasi-deterministic simulations. Overall, our results show that range expansions provide a robust mechanism for generating different types of multiple mergers, which could be similar to those observed in populations with strong selection or high fecundity. Thus, caution should be exercised in making inferences about the origin of non-Kingman genealogies.

Inference with selection, varying population size, and evolving population structure: application of ABC to a forward-backward coalescent process with interactions

Genetic data are often used to infer demographic history and changes or detect genes under selection. Inferential methods are commonly based on models making various strong assumptions: demography and population structures are supposed a priori known, the evolution of the genetic composition of a population does not affect demography nor population structure, and there is no selection nor interaction between and within genetic strains. In this paper, we present a stochastic birth-death model with competitive interactions and asexual reproduction. We develop an inferential procedure for ecological, demographic, and genetic parameters. We first show how genetic diversity and genealogies are related to birth and death rates, and to how individuals compete within and between strains. This leads us to propose an original model of phylogenies, with trait structure and interactions, that allows multiple merging. Second, we develop an Approximate Bayesian Computation framework to use our model for analyzing genetic data. We apply our procedure to simulated data from a toy model, and to real data by analyzing the genetic diversity of microsatellites on Y-chromosomes sampled from Central Asia human populations in order to test whether different social organizations show significantly different fertilities.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

Temporal migration rates affect the genetic structure of populations in the biennial Erysimum mediohispanicum with reproductive asynchrony

Migration is a process with important implications for the genetic structure of populations. However, there is an aspect of migration seldom investigated in plants: migration between temporally isolated groups of individuals within the same geographic population. The genetic implications of temporal migration can be particularly relevant for semelparous organisms, which are those that reproduce only once in a lifetime after a certain period of growth. In this case, reproductive asynchrony in individuals of the same population generates demes of individuals differing in their developmental stage (non-reproductive and reproductive). These demes are connected by temporal migrants, that is, individuals that become annually asynchronous with respect to the rest of individuals of their same deme. Here, we investigated the extent of temporal migration and its effects on temporal genetic structure in the biennial plant Erysimum mediohispanicum. To this end, we conducted two independent complementary approaches. First, we empirically estimated temporal migration rates and temporal genetic structure in four populations of E. mediohispanicum during three consecutive years using nuclear microsatellites markers. Second, we developed a demographic genetic simulation model to assess genetic structure for different migration scenarios differing in temporal migration rates and their occurrence probabilities. We hypothesized that genetic structure decreased with increasing temporal migration rates due to the homogenizing effect of migration. Empirical and modelling results were consistent and indicated a U-shape relationship between genetic structure and temporal migration rates. Overall, they indicated the existence of temporal genetic structure and that such genetic structure indeed decreased with increasing temporal migration rates. However, genetic structure increased again at high temporal migration rates. The results shed light into the effects of reproductive asynchrony on important population genetic parameters. Our study contributes to unravel the complexity of some processes that may account for genetic diversity and genetic structure of natural populations.

Range-wide genetic structure in the thorn-tailed rayadito suggests limited gene flow towards peripheral populations

Understanding the population genetic consequences of habitat heterogeneity requires assessing whether patterns of gene flow correspond to landscape configuration. Studies of the genetic structure of populations are still scarce for Neotropical forest birds. We assessed range-wide genetic structure and contemporary gene flow in the thorn-tailed rayadito (Aphrastura spinicauda), a passerine bird inhabiting the temperate forests of South America. We used 12 microsatellite loci to genotype 582 individuals from eight localities across a large latitudinal range (30°S–56°S). Using population structure metrics, multivariate analyses, clustering algorithms, and Bayesian methods, we found evidence for moderately low regional genetic structure and reduced gene flow towards the range margins. Genetic differentiation increased with geographic distance, particularly in the southern part of the species’ distribution where forests are continuously distributed. Populations in the north seem to experience limited gene flow likely due to forest discontinuity, and may comprise a demographically independent unit. The southernmost population, on the other hand, is genetically depauperate and different from all other populations. Different analytical approaches support the presence of three to five genetic clusters. We hypothesize that the genetic structure of the species follows a hierarchical clustered pattern.

A Few Stickleback Suffice for the Transport of Alleles to New Lakes

Threespine stickleback populations provide a striking example of local adaptation to divergent habitats in populations that are connected by recurrent gene flow. These small fish occur in marine and freshwater habitats throughout the Northern Hemisphere, and in numerous cases the smaller freshwater populations have been established “de novo” from marine colonists. Independently evolved freshwater populations exhibit similar phenotypes that have been shown to derive largely from the same standing genetic variants. Geographic isolation prevents direct migration between the freshwater populations, strongly suggesting that these shared locally adaptive alleles are transported through the marine population. However it is still largely unknown how gene flow, recombination, and selection jointly impact the standing variation that might fuel this adaptation. Here we use individual-based, spatially explicit simulations to determine the levels of gene flow that best match observed patterns of allele sharing among habitats in stickleback. We aim to better understand how gene flow and local adaptation in large metapopulations determine the speed of adaptation and re-use of standing genetic variation. In our simulations we find that repeated adaptation uses a shared set of alleles that are maintained at low frequency by migration-selection balance in oceanic populations. This process occurs over a realistic range of intermediate levels of gene flow that match previous empirical population genomic studies in stickleback. Examining these simulations more deeply reveals how lower levels of gene flow leads to slow, independent adaptation to different habitats, whereas higher levels of gene flow leads to significant mutation load – but an increased probability of successful population genomic scans for locally adapted alleles. Surprisingly, we find that the genealogical origins of most freshwater adapted alleles can be traced back to the original generation of marine individuals that colonized the lakes, as opposed to subsequent migrants. These simulations provide deeper context for existing studies of stickleback evolutionary genomics, and guidance for future empirical studies in this model. More broadly, our results support existing theory of local adaptation but extend it by more completely documenting the genealogical history of adaptive alleles in a metapopulation.

Spatial Population Genetics: It's About Time

Many important questions about the history and dynamics of organisms have a geographical component: How many are there, and where do they live? How do they move and interbreed across the landscape? How were they moving a thousand years ago, and where were the ancestors of a particular individual alive today? Answers to these questions can have profound consequences for our understanding of history, ecology, and the evolutionary process. In this review, we discuss how geographic aspects of the distribution, movement, and reproduction of organisms are reflected in their pedigree across space and time. Because the structure of the pedigree is what determines patterns of relatedness in modern genetic variation, our aim is to thus provide intuition for how these processes leave an imprint in genetic data. We also highlight some current methods and gaps in the statistical toolbox of spatial population genetics.

Isolation by phenology synergizes isolation by distance across a continuous landscape

Pollen is generally dispersed over short distances, which promotes population genetic structure across continuous two-dimensional space. Quantitative genetic variance in flowering time structures mating pools in the temporal dimension, at least with respect to the phenology loci. We asked if these two phenomena, isolation by distance (IBD) and isolation by phenology (IBP), synergistically promote genetic structure. We constructed an individual-based model that tracked genotype frequencies at flowering time and neutral loci across a uniform landscape, over multiple generations, under four mating schemes: panmixia, IBD only, IBP only, and IBP × IBD. IBD × IBP divided the population into spatial clusters of early-, mid-, and late-flowering genotypes and strongly increased its quantitative genetic variance. Flowering time did not cluster under IBP, but its genetic variance increased moderately. IBD induced mild spatial structure in a nonassortative reference trait but did not change its variance. Importantly, the spatial correlation of genotypes at neutral loci was twice as strong under IBD × IBP compared with IBD alone. IBD × IBP also drew neutral loci into gametic disequilibrium with flowering time loci, structuring them temporally. Temporal and spatial mating pool structure promotes local differentiation. This trend would facilitate adaptation on small spatial scales.

Selection, Linkage, and Population Structure Interact To Shape Genetic Variation Among Threespine Stickleback Genomes

The outcome of selection on genetic variation depends on the geographic organization of individuals and populations as well as the organization of loci within the genome. Spatially variable selection between marine and freshwater habitats has had a significant and heterogeneous impact on patterns of genetic variation across the genome of threespine stickleback fish. When marine stickleback invade freshwater habitats, more than a quarter of the genome can respond to divergent selection, even in as little as 50 years. This process largely uses standing genetic variation that can be found ubiquitously at low frequency in marine populations, can be millions of years old, and is likely maintained by significant bidirectional gene flow. Here, we combine population genomic data of marine and freshwater stickleback from Cook Inlet, Alaska, with genetic maps of stickleback fish derived from those same populations to examine how linkage to loci under selection affects genetic variation across the stickleback genome. Divergent selection has had opposing effects on linked genetic variation on chromosomes from marine and freshwater stickleback populations: near loci under selection, marine chromosomes are depauperate of variation, while these same regions among freshwater genomes are the most genetically diverse. Forward genetic simulations recapitulate this pattern when different selective environments also differ in population structure. Lastly, dense genetic maps demonstrate that the interaction between selection and population structure may impact large stretches of the stickleback genome. These findings advance our understanding of how the structuring of populations across geography influences the outcomes of selection, and how the recombination landscape broadens the genomic reach of selection.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

Local PCA Shows How the Effect of Population Structure Differs Along the Genome

Population structure leads to systematic patterns in measures of mean relatedness between individuals in large genomic data sets, which are often discovered and visualized using dimension reduction techniques such as principal component analysis (PCA). Mean relatedness is an average of the relationships across locus-specific genealogical trees, which can be strongly affected on intermediate genomic scales by linked selection and other factors. We show how to use local PCA to describe this intermediate-scale heterogeneity in patterns of relatedness, and apply the method to genomic data from three species, finding in each that the effect of population structure can vary substantially across only a few megabases. In a global human data set, localized heterogeneity is likely explained by polymorphic chromosomal inversions. In a range-wide data set of Medicago truncatula, factors that produce heterogeneity are shared between chromosomes, correlate with local gene density, and may be caused by linked selection, such as background selection or local adaptation. In a data set of primarily African Drosophila melanogaster, large-scale heterogeneity across each chromosome arm is explained by known chromosomal inversions thought to be under recent selection and, after removing samples carrying inversions, remaining heterogeneity is correlated with recombination rate and gene density, again suggesting a role for linked selection. The visualization method provides a flexible new way to discover biological drivers of genetic variation, and its application to data highlights the strong effects that linked selection and chromosomal inversions can have on observed patterns of genetic variation.

The IICR and the non-stationary structured coalescent: towards demographic inference with arbitrary changes in population structure

In the last years, a wide range of methods allowing to reconstruct past population size changes from genome-wide data have been developed. At the same time, there has been an increasing recognition that population structure can generate genetic data similar to those produced under models of population size change. Recently, Mazet et al. (Heredity 116:362-371, 2016) showed that, for any model of population structure, it is always possible to find a panmictic model with a particular function of population size changes, having exactly the same distribution of T₂ (the coalescence time for a sample of size two) as that of the structured model. They called this function IICR (Inverse Instantaneous Coalescence Rate) and showed that it does not necessarily correspond to population size changes under non-panmictic models. Besides, most of the methods used to analyse data under models of population structure tend to arbitrarily fix that structure and to minimise or neglect population size changes. Here, we extend the seminal work of Herbots (PhD thesis, University of London, 1994) on the structured coalescent and propose a new framework, the Non-Stationary Structured Coalescent (NSSC) that incorporates demographic events (changes in gene flow and/or deme sizes) to models of nearly any complexity. We show how to compute the IICR under a wide family of stationary and non-stationary models. As an example we address the question of human and Neanderthal evolution and discuss how the NSSC framework allows to interpret genomic data under this new perspective.

Is the sky the limit? On the expansion threshold of a species' range

More than 100 years after Grigg's influential analysis of species' borders, the causes of limits to species' ranges still represent a puzzle that has never been understood with clarity. The topic has become especially important recently as many scientists have become interested in the potential for species' ranges to shift in response to climate change-and yet nearly all of those studies fail to recognise or incorporate evolutionary genetics in a way that relates to theoretical developments. I show that range margins can be understood based on just two measurable parameters: (i) the fitness cost of dispersal-a measure of environmental heterogeneity-and (ii) the strength of genetic drift, which reduces genetic diversity. Together, these two parameters define an 'expansion threshold': adaptation fails when genetic drift reduces genetic diversity below that required for adaptation to a heterogeneous environment. When the key parameters drop below this expansion threshold locally, a sharp range margin forms. When they drop below this threshold throughout the species' range, adaptation collapses everywhere, resulting in either extinction or formation of a fragmented metapopulation. Because the effects of dispersal differ fundamentally with dimension, the second parameter-the strength of genetic drift-is qualitatively different compared to a linear habitat. In two-dimensional habitats, genetic drift becomes effectively independent of selection. It decreases with 'neighbourhood size'-the number of individuals accessible by dispersal within one generation. Moreover, in contrast to earlier predictions, which neglected evolution of genetic variance and/or stochasticity in two dimensions, dispersal into small marginal populations aids adaptation. This is because the reduction of both genetic and demographic stochasticity has a stronger effect than the cost of dispersal through increased maladaptation. The expansion threshold thus provides a novel, theoretically justified, and testable prediction for formation of the range margin and collapse of the species' range.

Genomic Signatures of Sexual Conflict

Sexual conflict is a specific class of intergenomic conflict that describes the reciprocal sex-specific fitness costs generated by antagonistic reproductive interactions. The potential for sexual conflict is an inherent property of having a shared genome between the sexes and, therefore, is an extreme form of an environment-dependent fitness effect. In this way, many of the predictions from environment-dependent selection can be used to formulate expected patterns of genome evolution under sexual conflict. However, the pleiotropic and transmission constraints inherent to having alleles move across sex-specific backgrounds from generation to generation further modulate the anticipated signatures of selection. We outline methods for detecting candidate sexual conflict loci both across and within populations. Additionally, we consider the ability of genome scans to identify sexually antagonistic loci by modeling allele frequency changes within males and females due to a single generation of selection. In particular, we highlight the need to integrate genotype, phenotype, and functional information to truly distinguish sexual conflict from other forms of sexual differentiation.

Regional and local environmental conditions do not shape the response to warming of a marine habitat-forming species

The differential response of marine populations to climate change remains poorly understood. Here, we combine common garden thermotolerance experiments in aquaria and population genetics to disentangle the factors driving the population response to thermal stress in a temperate habitat-forming species: the octocoral Paramuricea clavata. Using eight populations separated from tens of meters to hundreds of kilometers, which were differentially impacted by recent mortality events, we identify 25 °C as a critical thermal threshold. After one week of exposure at this temperature, seven of the eight populations were affected by tissue necrosis and after 30 days of exposure at this temperature, the mean % of affected colonies increased gradually from 3 to 97%. We then demonstrate the weak relation between the observed differential phenotypic responses and the local temperature regimes experienced by each population. A significant correlation was observed between these responses and the extent of genetic drift impacting each population. Local adaptation may thus be hindered by genetic drift, which seems to be the main driver of the differential response. Accordingly, conservation measures should promote connectivity and control density erosion in order to limit the impact of genetic drift on marine populations facing climate change.

Population genetic and field-ecological analyses return similar estimates of dispersal over space and time in an endangered amphibian

The explosive growth of empirical population genetics has seen a proliferation of analytical methods leading to a steady increase in our ability to accurately measure key population parameters, including genetic isolation, effective population size, and gene flow, in natural systems. Assuming they yield similar results, population genetic methods offer an attractive complement to, or replacement of, traditional field‐ecological studies. However, empirical assessments of the concordance between direct field‐ecological and indirect population genetic studies of the same populations are uncommon in the literature. In this study, we investigate genetic isolation, rates of dispersal, and population sizes for the endangered California tiger salamander, Ambystoma californiense, across multiple breeding seasons in an intact vernal pool network. We then compare our molecular results to a previously published study based on multiyear, mark–recapture data from the same breeding sites. We found that field and genetic estimates of population size were only weakly correlated, but dispersal rates were remarkably congruent across studies and methods. In fact, dispersal probability functions derived from genetic data and traditional field‐ecological data were a significant match, suggesting that either method can be used effectively to assess population connectivity. These results provide one of the first explicit tests of the correspondence between landscape genetic and field‐ecological approaches to measuring functional population connectivity and suggest that even single‐year genetic samples can return biologically meaningful estimates of natural dispersal and gene flow.

Imbalanced positive selection maintains the functional divergence of duplicated DIHYDROKAEMPFEROL 4-REDUCTASE genes

Gene duplication could be beneficial by functional division but might increase the risk of genetic load. The dynamics of duplicated paralogs number could involve recombination, positive selection, and functional divergence. Duplication of DIHYDROFLAVONOL 4-REDUCTASE (DFR) has been reported in several organisms and may have been retained by escape from adaptive conflict (EAC). In this study, we screened the angiosperm DFR gene focusing on a diversified genus Scutellaria to investigate how these duplicated genes are retained. We deduced that gene duplication involved multiple independent events in angiosperms, but the duplication of DFR was before the divergence of Scutellaria. Asymmetric positive selective pressures resulted in different evolutionary rates between the duplicates. Different numbers of regulatory elements, differential codon usages, radical amino acid changes, and differential gene expressions provide evidences of functional divergence between the two DFR duplicates in Scutellaria, implying adaptive subfunctionalization between duplicates. The discovery of pseudogenes accompanying a reduced replacement rate in one DFR paralogous gene suggested possibly leading to “loss of function” due to dosage imbalance after the transient adaptive subfunctionalization in the early stage of duplication. Notwithstanding, episodic gene duplication and functional divergence may be relevant to the diversification of ecological function of DFR gene in Scutellaria.

Can multi-generational exposure to ocean warming and acidification lead to the adaptation of life history and physiology in a marine metazoan?

Ocean warming and acidification are concomitant global drivers that are currently threatening the survival of marine organisms. How species will respond to these changes depends on their capacity for plastic and adaptive responses. Little is known about the mechanisms that govern plasticity and adaptability or how global changes will influence these relationships across multiple generations. Here, we exposed the emerging model marine polychaete Ophryotrocha labronica to conditions simulating ocean warming and acidification, in isolation and in combination over five generations to identify: (i) how multiple versus single global change drivers alter both juvenile and adult life-history traits; (ii) the mechanistic link between adult physiological and fitness-related life-history traits; and (iii) whether the phenotypic changes observed over multiple generations are of plastic and/or adaptive origin. Two juvenile (developmental rate; survival to sexual maturity) and two adult (average reproductive body size; fecundity) life-history traits were measured in each generation, in addition to three physiological (cellular reactive oxygen species content, mitochondrial density, mitochondrial capacity) traits. We found that multi-generational exposure to warming alone caused an increase in juvenile developmental rate, reactive oxygen species production and mitochondrial density, decreases in average reproductive body size and fecundity, and fluctuations in mitochondrial capacity, relative to control conditions. Exposure to ocean acidification alone had only minor effects on juvenile developmental rate. Remarkably, when both drivers of global change were present, only mitochondrial capacity was significantly affected, suggesting that ocean warming and acidification act as opposing vectors of stress across multiple generations.

Epigenetic Inheritance across the Landscape

The study of epigenomic variation at the landscape-level in plants may add important insight to studies of adaptive variation. A major goal of landscape genomic studies is to identify genomic regions contributing to adaptive variation across the landscape. Heritable variation in epigenetic marks, resulting in transgenerational plasticity, can influence fitness-related traits. Epigenetic marks are influenced by the genome, the environment, and their interaction, and can be inherited independently of the genome. Thus, epigenomic variation likely influences the heritability of many adaptive traits, but the extent of this influence remains largely unknown. Here, we summarize the relevance of epigenetic inheritance to ecological and evolutionary processes, and review the literature on landscape-level patterns of epigenetic variation. Landscape-level patterns of epigenomic variation in plants generally show greater levels of isolation by distance and isolation by environment then is found for the genome, but the causes of these patterns are not yet clear. Linkage between the environment and epigenomic variation has been clearly shown within a single generation, but demonstrating transgenerational inheritance requires more complex breeding and/or experimental designs. Transgenerational epigenetic variation may alter the interpretation of landscape genomic studies that rely upon phenotypic analyses, but should have less influence on landscape genomic approaches that rely upon outlier analyses or genome–environment associations. We suggest that multi-generation common garden experiments conducted across multiple environments will allow researchers to understand which parts of the epigenome are inherited, as well as to parse out the relative contribution of heritable epigenetic variation to the phenotype.

Excess of Deleterious Mutations around HLA Genes Reveals Evolutionary Cost of Balancing Selection

Deleterious mutations are expected to evolve under negative selection and are usually purged from the population. However, deleterious alleles segregate in the human population and some disease-associated variants are maintained at considerable frequencies. Here, we test the hypothesis that balancing selection may counteract purifying selection in neighboring regions and thus maintain deleterious variants at higher frequency than expected from their detrimental fitness effect. We first show in realistic simulations that balancing selection reduces the density of polymorphic sites surrounding a locus under balancing selection, but at the same time markedly increases the population frequency of the remaining variants, including even substantially deleterious alleles. To test the predictions of our simulations empirically, we then use whole-exome sequencing data from 6,500 human individuals and focus on the most established example for balancing selection in the human genome, the major histocompatibility complex (MHC). Our analysis shows an elevated frequency of putatively deleterious coding variants in nonhuman leukocyte antigen (non-HLA) genes localized in the MHC region. The mean frequency of these variants declined with physical distance from the classical HLA genes, indicating dependency on genetic linkage. These results reveal an indirect cost of the genetic diversity maintained by balancing selection, which has hitherto been perceived as mostly advantageous, and have implications both for the evolution of recombination and also for the epidemiology of various MHC-associated diseases.

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation-interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.

Comparative spatial genetics and epigenetics of plant populations: heuristic value and a proof of concept

Despite the recent upsurge of interest on natural epigenetic variation of nonmodel organisms, factors conditioning the spatial structure of epigenetic diversity in wild plant populations remain virtually unexplored. We propose that information on processes shaping natural epigenetic variation can be gained using the spatial structure of genetic diversity as null model. Departures of epigenetic isolation-by-distance (IBD) patterns from genetic IBD patterns for the same sample, particularly differences in slope of similarity-distance regressions, will reflect the action of factors that operate specifically on epigenetic variation, including imperfect transgenerational inheritance and responsiveness to environmental factors of epigenetic marks. As a proof of concept, we provide a comparative analysis of spatial genetic and epigenetic structure of 200 mapped individuals of the perennial herb Helleborus foetidus. Plants were fingerprinted using nuclear microsatellites, amplified fragment length polymorphisms (AFLP) and methylation-sensitive AFLP markers. Expectations from individual-level IBD patterns were tested by means of kinship-distance regressions. Both genetic and epigenetic similarity between H. foetidus individuals conformed to theoretical expectations under individual-level IBD models. Irrespective of marker type, there were significant negative linear relationships between the kinship coefficient for plant pairs and their spatial separation. Regression slopes were significantly steeper for epigenetic markers. Epigenetic similarity between individuals was much greater than genetic similarity at shortest distances, such epigenetic 'kinship excess' tending to decrease as plant separation increased. Results suggest that moderate-to-high heritability and responsiveness to local environments are major drivers of epigenetic spatial structure in H. foetidus, and illustrate the heuristic value of comparing genetic and epigenetic spatial structure for formulating and testing hypotheses on forces shaping epigenetic diversity in wild plant populations.

Mating system divergence affects the distribution of sequence diversity within and among populations of recently diverged subspecies of Clarkia xantiana (Onagraceae)

PREMISE OF THE STUDY

The population biology of outcrossing and self-fertilizing taxa is thought to differ because of the advantage that selfers have in colonizing unoccupied sites where mates and pollinators may be limiting (Baker's Law). This reduced tendency for outcrossers to colonize new sites, along with their greater dependence on pollinators to disperse pollen, has the potential to differently influence the genetic diversity and structure of outcrossing and selfing populations.

METHODS

We conducted a comparative population genetic study of two sister outcrossing and selfing subspecies of Clarkia xantiana that have very recently diverged. We used DNA sequence variation (>40 kb from eight nuclear loci) from large samples of individuals from 14 populations to assess geographic patterns of genetic diversity and make inferences about the demographic and colonization histories of each subspecies.

KEY RESULTS

We show that sequence variation is strongly reduced across all selfing populations. The demographic history of selfing populations exhibits recent colonization bottlenecks, whereas such bottlenecks are rarely observed for the outcrosser. The greater effect of genetic drift in the selfer has resulted in strong population genetic structure, but with no pattern of isolation by distance. By contrast, the stronger effect of gene flow in the outcrosser has resulted in considerably less structure, but a significant pattern of isolation by distance.

CONCLUSIONS

Taken together, our results suggest that selfing populations are not at migration-drift equilibrium, are affected by strong episodes of genetic drift during colonization, and experience little or no subsequent gene flow from other populations after those founder events.

On the importance of being structured: instantaneous coalescence rates and human evolution--lessons for ancestral population size inference?

Most species are structured and influenced by processes that either increased or reduced gene flow between populations. However, most population genetic inference methods assume panmixia and reconstruct a history characterized by population size changes. This is potentially problematic as population structure can generate spurious signals of population size change through time. Moreover, when the model assumed for demographic inference is misspecified, genomic data will likely increase the precision of misleading if not meaningless parameters. For instance, if data were generated under an n-island model (characterized by the number of islands and migrants exchanged) inference based on a model of population size change would produce precise estimates of a bottleneck that would be meaningless. In addition, archaeological or climatic events around the bottleneck's timing might provide a reasonable but potentially misleading scenario. In a context of model uncertainty (panmixia versus structure) genomic data may thus not necessarily lead to improved statistical inference. We consider two haploid genomes and develop a theory that explains why any demographic model with structure will necessarily be interpreted as a series of changes in population size by inference methods ignoring structure. We formalize a parameter, the inverse instantaneous coalescence rate, and show that it is equivalent to a population size only in panmictic models, and is mostly misleading for structured models. We argue that this issue affects all population genetics methods ignoring population structure which may thus infer population size changes that never took place. We apply our approach to human genomic data.

Deeply divergent sympatric mitochondrial lineages of the earthworm Lumbricus rubellus are not reproductively isolated

Background

The accurate delimitation of species is essential to numerous areas of biological research. An unbiased assessment of the diversity, including the cryptic diversity, is of particular importance for the below ground fauna, a major component of global biodiversity. On the British Isles, the epigeic earthworm Lumbricus rubellus, which is a sentinel species in soil ecotoxicology, consists of two cryptic taxa that are differentiated in both the nuclear and the mitochondrial (mtDNA) genomes. Recently, several deeply divergent mtDNA lineages were detected in mainland Europe, but whether these earthworms also constitute cryptic species remains unclear. This information is important from an evolutionary perspective, but it is also essential for the interpretation and the design of ecotoxicological projects. In this study, we used genome-wide RADseq data to assess the reproductive isolation of the divergent mitochondrial lineages of L. rubellus that occur in sympatry in multiple localities in Central Europe.

Results

We identified five divergent (up to 16 % net p-distance) mitochondrial lineages of L. rubellus in sympatry. Because the clustering of the RADseq data was according to the population of origin and not the mtDNA lineage, reproductive isolation among the mtDNA lineages was not likely. Although each population contained multiple mtDNA lineages, subdivisions within the populations were not observed for the nuclear genome. The lack of fixed differences and sharing of the overwhelming majority of nuclear polymorphisms between localities, indicated that the populations did not constitute allopatric species. The nucleotide diversity within the populations was high, 0.7–0.8 %.

Conclusions

The deeply divergent mtDNA sympatric lineages of L. rubellus in Central Europe were not reproductively isolated groups. The earthworm L. rubellus, which is represented by several mtDNA lineages in continental Europe, apparently is a single highly polymorphic species rather than a complex of several cryptic species. This study demonstrated the critical importance of the use of multilocus nuclear data for the unbiased assessment of cryptic diversity and for the delimitation of species in soil invertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0488-9) contains supplementary material, which is available to authorized users.

Ecology has contrasting effects on genetic variation within species versus rates of molecular evolution across species in water beetles

Comparative analysis is a potentially powerful approach to study the effects of ecological traits on genetic variation and rate of evolution across species. However, the lack of suitable datasets means that comparative studies of correlates of genetic traits across an entire clade have been rare. Here, we use a large DNA-barcode dataset (5062 sequences) of water beetles to test the effects of species ecology and geographical distribution on genetic variation within species and rates of molecular evolution across species. We investigated species traits predicted to influence their genetic characteristics, such as surrogate measures of species population size, latitudinal distribution and habitat types, taking phylogeny into account. Genetic variation of cytochrome oxidase I in water beetles was positively correlated with occupancy (numbers of sites of species presence) and negatively with latitude, whereas substitution rates across species depended mainly on habitat types, and running water specialists had the highest rate. These results are consistent with theoretical predictions from nearly-neutral theories of evolution, and suggest that the comparative analysis using large databases can give insights into correlates of genetic variation and molecular evolution.

Combining genetic and demographic data for the conservation of a Mediterranean marine habitat-forming species

The integration of ecological and evolutionary data is highly valuable for conservation planning. However, it has been rarely used in the marine realm, where the adequate design of marine protected areas (MPAs) is urgently needed. Here, we examined the interacting processes underlying the patterns of genetic structure and demographic strucuture of a highly vulnerable Mediterranean habitat-forming species (i.e. Paramuricea clavata (Risso, 1826)), with particular emphasis on the processes of contemporary dispersal, genetic drift, and colonization of a new population. Isolation by distance and genetic discontinuities were found, and three genetic clusters were detected; each submitted to variations in the relative impact of drift and gene flow. No founder effect was found in the new population. The interplay of ecology and evolution revealed that drift is strongly impacting the smallest, most isolated populations, where partial mortality of individuals was highest. Moreover, the eco-evolutionary analyses entailed important conservation implications for P. clavata. Our study supports the inclusion of habitat-forming organisms in the design of MPAs and highlights the need to account for genetic drift in the development of MPAs. Moreover, it reinforces the importance of integrating genetic and demographic data in marine conservation.

Nucleotide variation and balancing selection at the Ckma gene in Atlantic cod: analysis with multiple merger coalescent models

High-fecundity organisms, such as Atlantic cod, can withstand substantial natural selection and the entailing genetic load of replacing alleles at a number of loci due to their excess reproductive capacity. High-fecundity organisms may reproduce by sweepstakes leading to highly skewed heavy-tailed offspring distribution. Under such reproduction the Kingman coalescent of binary mergers breaks down and models of multiple merger coalescent are more appropriate. Here we study nucleotide variation at the Ckma (Creatine Kinase Muscle type A) gene in Atlantic cod. The gene shows extreme differentiation between the North (Canada, Greenland, Iceland, Norway, Barents Sea) and the South (Faroe Islands, North-, Baltic-, Celtic-, and Irish Seas) with F_ST > 0.8 between regions whereas neutral loci show no differentiation. This is evidence of natural selection. The protein sequence is conserved by purifying selection whereas silent and non-coding sites show extreme differentiation. The unfolded site-frequency spectrum has three modes, a mode at singleton sites and two high frequency modes at opposite frequencies representing divergent branches of the gene genealogy that is evidence for balancing selection. Analysis with multiple-merger coalescent models can account for the high frequency of singleton sites and indicate reproductive sweepstakes. Coalescent time scales vary with population size and with the inverse of variance in offspring number. Parameter estimates using multiple-merger coalescent models show that times scales are faster than under the Kingman coalescent.

Evolutionary Algorithm with Diversity-Reference Adaptive Control in Dynamic Environments

Evolutionary algorithms (EAs) can be used to find solutions in dynamic environments. In such cases, after a change in the environment, EAs can either be restarted or they can take advantage of previous knowledge to resume the evolutionary process. The second option tends to be faster and demands less computational effort. The preservation or growth of population diversity is one of the strategies used to advance the evolutionary process after modifications to the environment. We propose a new adaptive method to control population diversity based on a model-reference. The EA evolves the population whereas a control strategy, independently, handles the population diversity. Thus, the adaptive EA evolves a population that follows a diversity-reference model. The proposed model, called the Diversity-Reference Adaptive Control Evolutionary Algorithm (DRAC), aims to maintain or increase the population diversity, thus avoiding premature convergence, and assuring exploration of the solution space during the whole evolutionary process. We also propose a diversity models based on the dynamics of heterozygosity of the population, as models to be tracked by the diversity control. The performance of DRAC showed promising results when compared with the standard genetic algorithm and six other adaptive evolutionary algorithms in 14 different experiments with three different types of environments.