On the mechanistic nature of epistasis in a canonical cis -regulatory element

Understanding the relation between genotype and phenotype remains a major challenge. The difficulty of predicting individual mutation effects, and particularly the interactions between them, has prevented the development of a comprehensive theory that links genotypic changes to their phenotypic effects. We show that a general thermodynamic framework for gene regulation, based on a biophysical understanding of protein-DNA binding, accurately predicts the sign of epistasis in a canonical cis-regulatory element consisting of overlapping RNA polymerase and repressor binding sites. Sign and magnitude of individual mutation effects are sufficient to predict the sign of epistasis and its environmental dependence. Thus, the thermodynamic model offers the correct null prediction for epistasis between mutations across DNA-binding sites. Our results indicate that a predictive theory for the effects of cis-regulatory mutations is possible from first principles, as long as the essential molecular mechanisms and the constraints these impose on a biological system are accounted for

The Inexorable Spread of a Newly Arisen Neo-Y Chromosome

A newly arisen Y-chromosome can become established in one part of a species range by genetic drift or through the effects of selection on sexually antagonistic alleles. However, it is difficult to explain why it should then spread throughout the species range after this initial episode. As it spreads into new populations, it will actually enter females. It would then be expected to perform poorly since it will have been shaped by the selective regime of the male-only environment from which it came. We address this problem using computer models of hybrid zone dynamics where a neo-XY chromosomal race meets the ancestral karyotype. Our models consider that the neo-Y was established by the fusion of an autosome with the ancestral X-chromosome (thereby creating the Y and the ‘fused X’). Our principal finding is that sexually antagonistic effects of the Y induce indirect selection in favour of the fused X-chromosomes, causing their spread. The Y-chromosome can then spread, protected behind the advancing shield of the fused X distribution. This mode of spread provides a robust explanation of how newly arisen Y-chromosomes can spread. A Y-chromosome would be expected to accumulate mutations that would cause it to be selected against when it is a rare newly arrived migrant. The Y can spread, nevertheless, because of the indirect selection induced by gene flow (which can only be observed in models comprising multiple populations). These results suggest a fundamental re-evaluation of sex-chromosome hybrid zones. The well-understood evolutionary events that initiate the Y-chromosome's degeneration will actually fuel its range expansion

Limits to the Rate of Adaptive Substitution in Sexual Populations

In large populations, many beneficial mutations may be simultaneously available and may compete with one another, slowing adaptation. By finding the probability of fixation of a favorable allele in a simple model of a haploid sexual population, we find limits to the rate of adaptive substitution, , that depend on simple parameter combinations. When variance in fitness is low and linkage is loose, the baseline rate of substitution is , where is the population size, is the rate of beneficial mutations per genome, and is their mean selective advantage. Heritable variance in log fitness due to unlinked loci reduces by under polygamy and under monogamy. With a linear genetic map of length Morgans, interference is yet stronger. We use a scaling argument to show that the density of adaptive substitutions depends on , , , and only through the baseline density: . Under the approximation that the interference due to different sweeps adds up, we show that , implying that interference prevents the rate of adaptive substitution from exceeding one per centimorgan per 200 generations. Simulations and numerical calculations confirm the scaling argument and confirm the additive approximation for ; for higher , the rate of adaptation grows above , but only very slowly. We also consider the effect of sweeps on neutral diversity and show that, while even occasional sweeps can greatly reduce neutral diversity, this effect saturates as sweeps become more common—diversity can be maintained even in populations experiencing very strong interference. Our results indicate that for some organisms the rate of adaptive substitution may be primarily recombination-limited, depending only weakly on the mutation supply and the strength of selection

A simple mathematical model of gradual Darwinian evolution: Emergence of a Gaussian trait distribution in adaptation along a fitness gradient

We consider a simple mathematical model of gradual Darwinian evolution in continuous time and continuous trait space, due to intraspecific competition for common resource in an asexually reproducing population in constant environment, while far from evolutionary stable equilibrium. The model admits exact analytical solution. In particular, Gaussian distribution of the trait emerges from generic initial conditions.Comment: 21 pages, 2 figures, as accepted to J Math Biol 2013/03/1

Thermal adaptation and clinal mitochondrial DNA variation of European anchovy

Natural populations of widely distributed organisms often exhibit genetic clinal variation over their geographical ranges. The European anchovy, Engraulis encrasicolus, illustrates this by displaying a two-clade mitochondrial structure clinally arranged along the eastern Atlantic. One clade has low frequencies at higher latitudes, whereas the other has an anti-tropical distribution, with frequencies decreasing towards the tropics. The distribution pattern of these clades has been explained as a consequence of secondary contact after an ancient geographical isolation. However, it is not unlikely that selection acts on mitochondria whose genes are involved in relevant oxidative phosphorylation processes. In this study, we performed selection tests on a fragment of 1044 bp of the mitochondrial cytochrome b gene using 455 individuals from 18 locations. We also tested correlations of six environmental features: temperature, salinity, apparent oxygen utilization and nutrient concentrations of phosphate, nitrate and silicate, on a compilation of mitochondrial clade frequencies from 66 sampling sites comprising 2776 specimens from previously published studies. Positive selection in a single codon was detected predominantly (99%) in the anti-tropical clade and temperature was the most relevant environmental predictor, contributing with 59% of the variance in the geographical distribution of clade frequencies. These findings strongly suggest that temperature is shaping the contemporary distribution of mitochondrial DNA clade frequencies in the European anchovy.FCT [SFRH/BD/36600/2007, IF/00043/2012

Habitat Association and Seasonality in a Mosaic and Bimodal Hybrid Zone between Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae)

Understanding why some hybrid zones are bimodal and others unimodal can aid in identifying barriers to gene exchange following secondary contact. The hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi contains a mix of allopatric parental populations and inter-mingled bimodal and unimodal sympatric populations, and provides an ideal system to examine the roles of local selection and gene flow between populations in maintaining bimodality. However, it is first necessary to confirm, over a larger spatial scale, previously identified associations between population composition and season and habitat. Here we use cline-fitting of one morphological and one song trait along two valley transects, and intervening mountains, to confirm previously identified habitat associations (mountain versus valley) and seasonal changes in population composition. As expected from previous findings of studies on a smaller spatial scale, C. jacobsi dominated mountain habitats and mixed populations dominated valleys, and C. brunneus became more prevalent in August. Controlling for habitat and incorporating into the analysis seasonal changes in cline parameters and the standard errors of parental trait values revealed wider clines than previous studies (best estimates of 6.4 to 24.5 km in our study versus 2.8 to 4.7 km in previous studies) and increased percentage of trait variance explained (52.7% and 61.5% for transects 1 and 2 respectively, versus 17.6%). Revealing such strong and consistent patterns within a complex hybrid zone will allow more focused examination of the causes of variation in bimodality in mixed populations, in particular the roles of local selection versus habitat heterogeneity and gene flow between differentiated populations

Red Queen Dynamics with Non-Standard Fitness Interactions

Antagonistic coevolution between hosts and parasites can involve rapid fluctuations of genotype frequencies that are known as Red Queen dynamics. Under such dynamics, recombination in the hosts may be advantageous because genetic shuffling can quickly produce disproportionately fit offspring (the Red Queen hypothesis). Previous models investigating these dynamics have assumed rather simple models of genetic interactions between hosts and parasites. Here, we assess the robustness of earlier theoretical predictions about the Red Queen with respect to the underlying host-parasite interactions. To this end, we created large numbers of random interaction matrices, analysed the resulting dynamics through simulation, and ascertained whether recombination was favoured or disfavoured. We observed Red Queen dynamics in many of our simulations provided the interaction matrices exhibited sufficient ‘antagonicity’. In agreement with previous studies, strong selection on either hosts or parasites favours selection for increased recombination. However, fast changes in the sign of linkage disequilibrium or epistasis were only infrequently observed and do not appear to be a necessary condition for the Red Queen hypothesis to work. Indeed, recombination was often favoured even though the linkage disequilibrium remained of constant sign throughout the simulations. We conclude that Red Queen-type dynamics involving persistent fluctuations in host and parasite genotype frequencies appear to not be an artefact of specific assumptions about host-parasite fitness interactions, but emerge readily with the general interactions studied here. Our results also indicate that although recombination is often favoured, some of the factors previously thought to be important in this process such as linkage disequilibrium fluctuations need to be reassessed when fitness interactions between hosts and parasites are complex

Hybridization and pre-zygotic reproductive barriers in Plasmodium

Host-parasite interactio

Asymmetric reproductive isolation between terminal forms of the salamander ring species Ensatina eschscholtzii revealed by fine-scale genetic analysis of a hybrid zone

<p>Abstract</p> <p>Background</p> <p>Ring species, exemplified by salamanders of the <it>Ensatina eschscholtzii </it>complex, represent a special window into the speciation process because they allow the history of species formation to be traced back in time through the geographically differentiated forms connecting the two terminal forms of the ring. Of particular interest is the nature and extent of reproductive isolation between the geographically terminal forms, in this case <it>E. e. eschscholtzii </it>and <it>E. e. klauberi</it>. Previous studies have documented infrequent hybridization at the end of the ring. Here, we report the first fine-scale genetic analysis of a hybrid zone between the terminal forms in southern California using individual-based Bayesian analyses of multilocus genetic data to estimate levels and direction of hybridization and maximum-likelihood analysis of linkage disequilibrium and cline shape to make inferences about migration and selection in the hybrid zone.</p> <p>Results</p> <p>The center of the hybrid zone has a high proportion of hybrids, about half of which were classified as F1s. Clines are narrow with respect to dispersal, and there are significant deviations from Hardy-Weinberg equilibrium as well as nonrandom associations (linkage disequilibria) between alleles characteristic of each parental type. There is cytonuclear discordance, both in terms of introgression and the geographic position of mitochondrial versus nuclear clines. Genetic disequilibrium is concentrated on the <it>eschscholtzii </it>side of the zone. Nearly all hybrids possess <it>klauberi </it>mtDNA, indicating that most hybrids are formed from female <it>klauberi </it>mating with male <it>eschscholtzii </it>or male hybrids (but not vice versa).</p> <p>Conclusions</p> <p>Our results are consistent with a tension zone trapped at an ecotone, with gene combinations characteristic of <it>klauberi </it>showing up on the <it>eschscholtzii </it>side of the zone due to asymmetric hybridization. We suggest that the observed asymmetry is best explained by increased discriminatory power of <it>eschscholtzii </it>females, or asymmetric postzygotic isolation. The relatively high frequency of hybrids, particularly F1s, contrasts with other contacts between the terminal forms, and with other contacts between other divergent <it>Ensatina </it>lineages, highlighting the diverse outcomes of secondary contact within a single species complex.</p