Inferring Genome-Wide Correlations of Mutation Fitness Effects between Populations

The effect of a mutation on fitness may differ between populations depending on environmental and genetic context, but little is known about the factors that underlie such differences. To quantify genome-wide correlations in mutation fitness effects, we developed a novel concept called a joint distribution of fitness effects (DFE) between populations. We then proposed a new statistic w to measure the DFE correlation between populations. Using simulation, we showed that inferring the DFE correlation from the joint allele frequency spectrum is statistically precise and robust. Using population genomic data, we inferred DFE correlations of populations in humans, Drosophila melanogaster, and wild tomatoes. In these species, we found that the overall correlation of the joint DFE was inversely related to genetic differentiation. In humans and D. melanogaster, deleterious mutations had a lower DFE correlation than tolerated mutations, indicating a complex joint DFE. Altogether, the DFE correlation can be reliably inferred, and it offers extensive insight into the genetics of population divergence.     

Accumulation of deleterious mutations: additional Drosophila melanogaster estimates and a simulation of the effects of selection

We report an assay of egg-to-adult viability in full-sibling mutation accumulation (MA) lines derived from a completely homozygous population of Drosophila melanogaster and maintained for 210 generations. A simultaneous evaluation was also made of a large population derived from the same origin and maintained as a control for the same period. We also present computer simulations to explore the possible decline in viability of the control population due to mutation accumulation and the possible effect of selection within and between MA lines. For this purpose, we used two mutational models independent from the data analyzed and based on radically different assumptions. The first model implies a large number of mutations of small effect, whereas the second implies a much smaller number of mutations with much larger effects. The observed rate of decline in mean viability was very small but significant (0.077%). The rate of increase in among line variance (0.189 x 10(-3)) was similar to those obtained previously in the same lines. The simulation results indicated that a model of many mutations of small effect is incompatible with the evolution of the mean viability of the control and MA lines over generations, the distribution of line means after 210 generations of mutation accumulation, and the pattern of line extinction over generations. Basically, this model predicted a large drop in viability, both in the control and particularly the MA lines, that is not observed empirically. It also predicted a rate of line extinction too low in the early generations and too high in the later ones. In contrast, the model based on few mutations of large effect was generally consistent with all the observations.     

Sexual selection is ineffectual or inhibits the purging of deleterious mutations in Drosophila melanogaster

The effects of sexual selection on population mean fitness are unclear and a subject of debate. Recent models propose that, because reproductive success may be condition dependent, much of the genome may be a target of sexual selection. Under this scenario, mutations that reduce health, and thus nonsexual fitness, may also be deleterious with respect to reproductive success, meaning that sexual selection may contribute to the purging of deleterious alleles. We tested this hypothesis directly by subjecting replicate Drosophila melanogaster populations to two treatments that altered the opportunity for sexual selection and then tracked changes in the frequency of six separate deleterious alleles with recessive and visible phenotypic effects. While natural selection acted to decrease the frequency of all six mutations, the addition of sexual selection did not aid in the purging of any of them, and for three of them appears to have hampered it. Courtship and mating have harmful effects in this species and mate choice assays showed that males directed more courtship and mating behavior toward wild-type over mutant females, providing a likely explanation for sexual selection's cost. Whether this cost extends to other mutations (e.g., those lacking visible phenotypic effects) is an important topic for future research.     

Mutation accumulation in populations of varying size: large effect mutations cause most mutational decline in the rotifer Brachionus calyciflorus under UV‐C radiation

Theory predicts that fitness decline via mutation accumulation will depend on population size, but there are only a few direct tests of this key idea. To gain a qualitative understanding of the fitness effect of new mutations, we performed a mutation accumulation experiment with the facultative sexual rotifer Brachionus calyciflorus at six different population sizes under UV‐C radiation. Lifetime reproduction assays conducted after ten and sixteen UV‐C radiations showed that while small populations lost fitness, fitness losses diminished rapidly with increasing population size. Populations kept as low as 10 individuals were able to maintain fitness close to the nonmutagenized populations throughout the experiment indicating that selection was able to remove the majority of large effect mutations in small populations. Although our results also seem to imply that small populations are effectively immune to mutational decay, we caution against this interpretation. Given sufficient time, populations of moderate to large size can experience declines in fitness from accumulating weakly deleterious mutations as demonstrated by fitness estimates from simulations and, tentatively, from a long‐term experiment with populations of moderate size. There is mounting evidence to suggest that mutational distributions contain a heavier tail of large effects. Our results suggest that this is also true when the mutational spectrum is altered by UV radiation.     

Color patterns and species recognition in four closely related species of Lake Malawi cichlid

We examined interspecific female mating preferences in fourclosely related species of cichlid belonging to the Pseudotropheuszebra species complex of Lake Malawi. These species differin coloration but are similar in other respects, suggestingthat male color patterns may be important to female mate choicein species recognition. To test this hypothesis, we presented females from each species with a choice of four males, one ofher own species and three others that were each of a differentspecies. We also gave each female a choice between the threeheterospecific males only. In all four species, females showeda significant preference for conspecific males in the four-waychoice and chose the male with the most similar color patternto the conspecific male in the three-way choice. These resultsare discussed with reference to the theory of sexual selectionon color patterns as a means of sympatric speciation in cichlids.     

Phenotypic differentiation is associated with divergent sexual selection among closely related barn swallow populations

Sexual selection plays a key role in the diversification of numerous animal clades and may accelerate trait divergence during speciation. However, much of our understanding of this process comes from phylogenetic comparative studies, which rely on surrogate measures such as dimorphism that may not represent selection in wild populations. In this study, we assess sexual selection pressures for multiple male visual signals across four barn swallow (Hirundo rustica) populations. Our sample encompassed 2400 linear km and two described subspecies: European H. r. rustica (in the Czech Republic and Romania) and eastern Mediterranean H. r. transitiva (in Israel), as well as a potential area of contact (in Turkey). We demonstrate significant phenotypic differentiation in four sexual signalling axes, despite very low‐level genomic divergence and no comparable divergence in an ecological trait. Moreover, the direction of phenotypic divergence is consistent with differences in sexual selection pressures among subspecies. Thus, H. r. transitiva, which have the darkest ventral plumage of any population, experience directional selection for darker plumage. Similarly, H. r. rustica, which have the longest tail feathers of any population, experience directional selection for elongated tail feathers and disruptive selection for ventral plumage saturation. These results suggest that sexual selection is the primary driver of phenotypic differentiation in this species. Our findings add to growing evidence of phenotypic divergence with gene flow. However, to our knowledge, this is the first study to relate direct measures of the strength and targets of sexual selection to phenotypic divergence among closely related wild populations.     

Complex genetic patterns in closely related colonizing invasive species

Anthropogenic activities frequently result in both rapidly changing environments and translocation of species from their native ranges (i.e., biological invasions). Empirical studies suggest that many factors associated with these changes can lead to complex genetic patterns, particularly among invasive populations. However, genetic complexities and factors responsible for them remain uncharacterized in many cases. Here, we explore these issues in the vase tunicate Ciona intestinalis (Ascidiacea: Enterogona: Cionidae), a model species complex, of which spA and spB are rapidly spreading worldwide. We intensively sampled 26 sites (N = 873) from both coasts of North America, and performed phylogenetic and population genetics analyses based on one mitochondrial fragment (cytochrome c oxidase subunit 3-NADH dehydrogenase subunit I, COX3-ND1) and eight nuclear microsatellites. Our analyses revealed extremely complex genetic patterns in both species on both coasts. We detected a contrasting pattern based on the mitochondrial marker: two major genetic groups in C. intestinalis spA on the west coast versus no significant geographic structure in C. intestinalis spB on the east coast. For both species, geo-graphically distant populations often showed high microsatellite-based genetic affinities whereas neighboring ones often did not. In addition, mitochondrial and nuclear markers provided largely inconsistent genetic patterns. Multiple factors, including random genetic drift associated with demographic changes, rapid selection due to strong local adaptation, and varying propensity for human-mediated propagule dispersal could be responsible for the observed genetic complexities.     

Interplay between extreme drift and selection intensities favors the fixation of beneficial mutations in selfing maize populations

Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.     

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

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

Fitness and entropy production in a cell population dynamics with epigenetic phenotype switching

Motivated by recent understandings in the stochastic natures of gene expression, biochemical signaling, and spontaneous reversible epigenetic switchings, we study a simple deterministic cell population dynamics in which subpopulations grow with different rates and individual cells can bi-directionally switch between a small number of different epigenetic phenotypes. Two theories in the past, the population dynamics and thermodynamics of master equations, separately defined two important concepts in mathematical terms: the fitness in the former and the (non-adiabatic) entropy production in the latter. Both of them play important roles in the evolution of the cell population dynamics. The switching sustains the variations among the subpopulation growth, thus sustains continuous natural selection. As a form of Price’s equation, the fitness increases with (i) natural selection through variations and (ii) a positive covariance between the per capita growth and switching, which represents a Lamarchian-like behavior. A negative covariance balances the natural selection in a fitness steady state --- “the red queen” scenario. At the same time the growth keeps the proportions of subpopulations away from the “intrinsic” switching equilibrium of individual cells, thus leads to a continuous entropy production. A covariance, between the per capita growth rate and the “chemical potential” of subpopulation, counteracts the entropy production. Analytical results are obtained for the limiting cases of growth dominating switching and vice versa.     

Fitness effects of new mutations in Chlamydomonas reinhardtii across two stress gradients

Most spontaneous mutations affecting fitness are likely to be deleterious, but the strength of selection acting on them might be impacted by environmental stress. Such stress‐dependent selection could expose hidden genetic variation, which in turn might increase the adaptive potential of stressed populations. On the other hand, this variation might represent a genetic load and thus lead to population extinction under stress. Previous studies to determine the link between stress and mutational effects on fitness, however, have produced inconsistent results. Here, we determined the net change in fitness in 29 genotypes of the green algae Chlamydomonas reinhardtii that accumulated mutations in the near absence of selection for approximately 1000 generations across two stress gradients, increasing NaCl and decreasing phosphate. We found mutational effects to be magnified under extremely stressful conditions, but such effects were specific both to the type of stress and to the genetic background. The detection of stress‐dependent fitness effects of mutations depended on accurately scaling relative fitness measures by generation times, thus offering an explanation for the inconsistencies among previous studies.     

Impact of demography on linked selection in two outcrossing Brassicaceae species

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.     

Interspecific interactions influence contrasting spatial genetic structures in two closely related damselfly species

Spatial genetic structure (SGS) is largely determined by colonization history, landscape and ecological characteristics of the species. Therefore, sympatric and ecologically similar species are expected to exhibit similar SGSs, potentially enabling prediction of the SGS of one species from that of another. On the other hand, due to interspecific interactions, ecologically similar species could have different SGSs. We explored the SGSs of the closely related Calopteryx splendens and Calopteryx virgo within Finland and related the genetic patterns to characteristics of the sampling localities. We observed different SGSs for the two species. Genetic differentiation even within short distances in C. splendens suggests genetic drift as an important driver. However, we also observed indication of previous gene flow (revealed by a negative relationship between genetic differentiation and increasing potential connectivity of the landscape). Interestingly, genetic diversity of C. splendens was negatively related to density of C. virgo, suggesting that interspecific interactions influence the SGS of C. splendens. In contrast, genetic differentiation between C. virgo subpopulations was low and only exhibited relationships with latitude, pointing to high gene flow, colonization history and range margin effects as the drivers of SGS. The different SGSs of the two ecologically similar species caution indirect inferences of SGS based on ecologically similar surrogate species.     

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.     

A spatial genetic structure and effects of relatedness on mate choice in a wild bird population

Inbreeding depression, as commonly found in natural populations, should favour the evolution of inbreeding avoidance mechanisms. If natal dispersal, the first and probably most effective mechanism, does not lead to a complete separation of males and females from a common origin, a small-scale genetic population structure may result and other mechanisms to avoid inbreeding may exist. We studied the genetic population structure and individual mating patterns in blue tits (Parus caeruleus). The population showed a local genetic structure in two out of four years: genetic relatedness between individuals (estimated from microsatellite markers) decreased with distance. This pattern was mainly caused by immigrants to the study area; these, if paired with fellow immigrants, were more related than expected by chance. Since blue tits did not avoid inbreeding with their social partner, we examined if individuals preferred less related partners at later stages of the mate choice process. We found no evidence that females or males avoided inbreeding through extra-pair copulations or through mate desertion and postbreeding dispersal. Although the small-scale genetic population structure suggests that blue tits could use a simple rule of thumb to select less related mates, females did not generally prefer more distantly breeding extra-pair partners. However, the proportion of young fathered by an extra-pair male in mixed paternity broods depended on the genetic relatedness with the female. This suggests that there is a fertilization bias towards less related copulation partners and that blue tits are able to reduce the costs of inbreeding through a postcopulatory process.     

Extending the range of additivity in using inclusive fitness

Inclusive fitness is a concept widely utilized by social biologists as the quantity organisms appear designed to maximize. However, inclusive fitness theory has long been criticized on the (uncontested) grounds that other quantities, such as offspring number, predict gene frequency changes accurately in a wider range of mathematical models. Here, we articulate a set of modeling assumptions that extend the range of scenarios in which inclusive fitness can be applied. We reanalyze recent formal analyses that searched for, but did not find, inclusive fitness maximization. We show (a) that previous models have not used Hamilton''s definition of inclusive fitness, (b) a reinterpretation of Hamilton''s definition that makes it usable in this context, and (c) that under the assumption of probabilistic mixing of phenotypes, inclusive fitness is indeed maximized in these models. We also show how to understand mathematically, and at an individual level, the definition of inclusive fitness, in an explicit population genetic model in which exact additivity is not assumed. We hope that in articulating these modeling assumptions and providing formal support for inclusive fitness maximization, we help bridge the gap between empiricists and theoreticians, which in some ways has been widening, demonstrating to mathematicians why biologists are content to use inclusive fitness, and offering one way to utilize inclusive fitness in general models of social behavior.     

Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species

The wild currant tomato Solanum pimpinellifolium inhabits a wide range of abiotic habitats across its native range of Ecuador and Peru. Although it has served as a key genetic resource for the improvement of domestic cultivars, little is known about the genetic basis of traits underlying local adaptation in this species, nor what abiotic variables are most important for driving differentiation. Here we use redundancy analysis (RDA) and other multivariate statistical methods (structural equation modelling [SEM] and generalized dissimilarity modelling [GDM]) to quantify the relationship of genomic variation (6,830 single nucleotide polymorphisms [SNPs]) with climate and geography, among 140 wild accessions. RDA, SEM and GDM each identified environment as explaining more genomic variation than geography, suggesting that local adaptation to heterogeneous abiotic habitats may be an important source of genetic diversity in this species. Environmental factors describing temporal variation in precipitation and evaporative demand explained the most SNP variation among accessions, indicating that these forces may represent key selective agents. Lastly, by studying how SNP–environment associations vary throughout the genome (44,064 SNPs), we mapped the location and investigated the functions of loci putatively contributing to climatic adaptations. Together, our findings indicate an important role for selection imposed by the abiotic environment in driving genomic differentiation between populations.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

The spatial population structure of a harpacticoid copepod community in spring

The spatial population patterns of an assemblage of meiobenthic harpacticoid copepod species were analysed for two sets of samples resulting from a systematic survey of an intertidal, estuarine, sandy beach in spring. The fauna was characterised by mesopsammic species, the 6 most abundant of which showed a great degree of interassociation, such that they could not be considered to distribute themselves independently. Populations of each of the 6 species were overdispersed, such that a negative binomial distribution adequately fitted all population samples. Spatial segregation of closely related species was indicated during the late spring sample, so leading to a patchy population distribution. It is proposed that during early spring, during periods of extreme disturbance in the seasonal r-selecting environment, the species populations do not interact, but that during periods of less extreme perturbation an interactive community evolves.