Role of biased gene conversion in one-locus neutral theory and genome evolution

The implications of biased gene conversion acting on selectively neutral alleles are investigated for a single diallelic locus in a finite population. Even a very slight conversion bias can significantly alter fixation probabilities. We argue that most newly arising mutants will be at a conversion disadvantage, resulting in a potentially greatly decreased substitution rate of new alleles compared with predictions from strict neutral theory. Thus, conversion bias potential allows for conservation of particular alleles without having to invoke selection. Conversely, we also show that bias can be important in the maintenance of repeated gene families without altering the substitution rate at other loci that experience the same amount of conversion bias, provided that the number of genes in the family is sufficiently large. Bias can, therefore, be important at the genomic level and yet be unimportant at the populational level. Finally, we discuss the role of biased gene conversion in speciation events, concluding that this type of molecular turnover acting independently at many individual loci is very unlikely to decrease the time required for two allopatric populations to speciate.     

Allelic divergence precedes and promotes gene duplication

One of the striking observations from recent whole-genome comparisons is that changes in the number of specialized genes in existing gene families, as opposed to novel taxon-specific gene families, are responsible for the majority of the difference in genome composition between major taxa. Previous models of duplicate gene evolution focused primarily on the role that neutral processes can play in evolutionary divergence after the duplicates are already fixed in the population. By instead including the entire cycle of duplication and divergence, we show that specialized functions are most likely to evolve through strong selection acting on segregating alleles at a single locus, even before the duplicate arises. We show that the fitness relationships that allow divergent alleles to evolve at a single locus largely overlap with the conditions that allow divergence of previously duplicated genes. Thus, a solution to the paradox of the origin of organismal complexity via the expansion of gene families exists in the form of the deterministic spread of novel duplicates via natural selection.     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

Genetic Variation in a Heterogeneous Environment. II. Temporal Heterogeneity and Directional Selection

The maintenance of genetic variation is investigated in a finite population where selection at an autosomal locus with two alleles varies temporally between two environments and the heterozygote has an intermediate fitness value. When there is additive gene action and equal selection in both environments, the autocorrelation between subsequent environments must be negative for more maintenance of genetic variation than for neutrality. The maximum maintenance occurs when there is equal selection in the two environments and the autocorrelation approaches ^-1.0 (for a stochastic model), or when there is short repeating cycle such as one related to seasons. Also comparison of the effects of stochastic variation in selection in finite and infinite populations is made by using Monte Carlo simulation. One situation was found where temporal environmental variation maintains genetic variation very effectively even in a small population and that is when there is evolution of dominance, i.e., the heterozygote is closer in fitness to the favored homozygote than the other homozygote. An important conclusion is that in a finite population genetic tracing of environmental change, particularly when there is a positive autocorrelation between environments or a long environmental cycle, leads to an increased loss of genetic variation making such a response undesirable in the long term, a result different from that in infinite populations.     

Interchromosomal Biased Gene Conversion, Mutation and Selection in a Multigene Family

A mathematical model of the effects of interchromosomal biased gene conversion, mutation and natural selection on a multigene family is developed and analyzed. The model assumes two allelic states at each of n loci. The effects of genetic drift are ignored. The model is developed under the assumption of no recombination, but the analysis shows that, at equilibrium, there is no linkage disequilibrium, which implies that the conclusions are valid for arbitrary recombination among loci. At equilibrium, the balance between mutation, gene conversion and selection depends on the ratio of the mutation rates to the quantity [s + g(2α - 1)/ n], where s is the increment or decrement in relative fitness with each additional copy of one of the alleles, g is the conversion rate, and α is a measure of the bias in favor of one of the alleles. When this quantity is large relative to the mutation rates, the allele that has the net advantage, combining the effects of selection and conversion, will be nearly fixed in the multigene family. A comparison of these results with those from a comparable model of intrachromosomal biased conversion shows that biased interchromosomal conversion leads to approximately the same equilibrium copy number as does intrachromosomal conversion of the same strength. Interchromosomal conversion is much more effective in causing the substitution of one allele by another. The relative frequencies of interchromosomal and intrachromosomal conversion is indicated by the extent of the linkage disequilibrium among the loci in a multigene family.     

The effect of linkage and population size on inbreeding depression due to mutational load.

Using a stochastic model of a finite population in which there is mutation to partially recessive detrimental alleles at many loci, we study the effects of population size and linkage between the loci on the population mean fitness and inbreeding depression values. Although linkage between the selected loci decreases the amount of inbreeding depression, neither population size nor recombination rate have strong effects on these quantities, unless extremely small values are assumed. We also investigate how partial linkage between the loci that determine fitness affects the invasion of populations by alleles at a modifier locus that controls the selfing rate. In most of the cases studied, the direction of selection on modifiers was consistent with that found in our previous deterministic calculations. However, there was some evidence that linkage between the modifier locus and the selected loci makes outcrossing less likely to evolve; more losses of alleles promoting outcrossing occurred in runs with linkage than in runs with free recombination. We also studied the fate of neutral alleles introduced into populations carrying detrimental mutations. The times to loss of neutral alleles introduced at low frequency were shorter than those predicted for alleles in the absence of selected loci, taking into account the reduction of the effective population size due to inbreeding. Previous studies have been confined to outbreeding populations, and to alleles at frequencies close to one-half, and have found an effect in the opposite direction. It therefore appears that associations between neutral and selected loci may produce effects that differ according to the initial frequencies of the neutral alleles.     

Polymorphism and Balancing Selection at Major Histocompatibility Complex Loci

Amino acid replacements in the peptide-binding region (PBR) of the functional major histocompatibility complex (Mhc) genes appear to be driven by balancing selection. Of the various types of balancing selection, we have examined a model equivalent to overdominance that confers heterozygote advantage. As discussed by A. Robertson, overdominance selection tends to maintain alleles that have more or less the same degree of heterozygote advantage. Because of this symmetry, the model makes various testable predictions about the genealogical relationships among different alleles and provides ways of analyzing DNA sequences of Mhc alleles. In this paper, we analyze DNA sequences of 85 alleles at the HLA-A, -B, -C, -DRB1 and -DQB1 loci with respect to the number of alleles and extent of nucleotide differences at the PBR, as well as at the synonymous (presumably neutral) sites. Theory suggests that the number of alleles that differ at the sites targeted by selection (presumably the nonsynonymous sites in the PBR) should be equal to the mean number of nucleotide substitutions among pairs of alleles. We also demonstrate that the nucleotide substitution rate at the targeted sites relative to that of neutral sites may be much larger than 1. The predictions of the presented model are in surprisingly good agreement with the actual data and thus provide means for inferring certain population parameters. For overdominance selection in a finite population at equilibrium, the product of selection intensity (s) against homozygotes and the effective population size (N) is estimated to be 350-3000, being largest at the B locus and smallest at the C locus. We argue that N is of the order of 10(5) and s is several percent at most, if the mutation rate per site per generation is 10(-8).     

Evolution of the segregation ratio: modification of gene conversion and meiotic drive

We compare the evolutionary pressures that direct the modification of gene conversion and meiotic drive at loci subject to purifying and overdominant viability selection. Gene conversion differs from meiotic drive in that modifers do not affect their own segregation ratios, even when linked to the viability locus. Segregation distortion generates gametic level disequilibria between alleles at the viability locus and modifiers of gene conversion and meiotic drive: enhancers of segregation distortion become positively associated with driven alleles. Suppression of gene conversion evolves if the driven allele is marginally disadvantageous (overdominant viability selection), and higher rates evolve if the driven alleles are relatively advantageous (purifying viability selection). Gametic disequilibria permit enhancers of meiotic drive that are linked to the driven locus to promote their own segregation. We attribute the failure of genetic modifiers of gene conversion and meiotic drive to maximinize mean fitness to the generation of such associations.     

The evolutionary rate of duplicated genes under concerted evolution

The effect of directional selection on the fixation process of a single mutation that spreads in a multigene family by gene conversion is investigated. A simple two-locus model with two alleles, A and a, is first considered in a random-mating diploid population with size N. There are four haplotypes, AA, Aa, aA, and aa, and selection works on the number of alleles A in a diplod (i = 0, 1, 2, 3, 4). Because gene conversion is allowed between the two loci, when the mutation rate is very low, either AA or aa will fix in the population eventually. We consider a situation where a single mutant, A, arises in one locus when a is fixed in both loci. Then, we derive the fixation probability analytically, and the fixation time is investigated by simulations. It is found that gene conversion has an effect to increase the "effective" population size, so that weak selection works more efficiently in a multigene family. With these results, we discuss the effect of gene conversion on the rate of molecular evolution in a multigene family undergoing concerted evolution. We also argue about the applicability of the theoretical results to models of multigene families with more than two loci.     

Recombination,selection, and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.     

Genetic Differentiation of Transposable Elements under Mutation and Unbiased Gene Conversion

A model is developed to predict the extent of genetic differentiation in a family of transposable elements under the combined effects of genetic drift, transposition, mutation and unbiased gene conversion. The model is based on simplifying assumptions that are valid when transposition is always to new sites and copy number per site is low. In the absence of gene conversion, the degree of differentiation as measured by the probability of identity of different elements is the same as at a single locus with the same mutation rate but in a population of effective size Nc/2, where N is the population size and c is the number of copies per individual. The inclusion of unbiased gene conversion does not significantly change this result. If, as seems to be the case, families of transposable elements are relatively homogeneous, then the model implies either that mutation rates for transposable elements are much lower than at comparable single-copy loci or that some other force, such as natural selection or biased gene conversion, is at work. Transposition is a very ineffective force for homogenizing a family of transposable elements.     

Inclusive fitness in finite populations—effects of heterogeneity and synergy

I review recent results concerning the relationship between the inclusive fitness (IF) effect and standard measures of allele fitness in a finite‐population, with attention to the effect of heterogeneity in population structure and nonadditive fitness effects. In both cases, existing theoretical work is somewhat technical and I try to provide a more transparent account. In a heterogeneous population it is known that inclusive fitness will generally fail to incorporate the effects of selection on the distribution of alleles among states unless a reproductive‐value weighting is used. But even given that, recent work shows that under certain updating rules, the IF effect can fail to be equivalent to standard measures such as fixation probability. In terms of synergistic fitness effects, I review the result that in the finite population model, the IF effect can be calculated using only “additive” relatedness coefficients so that computational difficulties found in the infinite‐population model do not arise. In my own work, there is an interaction here in that my 2012 work on synergy with Maciejewski made an assumption about inclusive fitness that my 2014 work on heterogeneity with Tarnita showed to be wrong. I include (Appendix C) a corrected argument for the 2012 result.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

A gene's eye view of epistasis,selection and speciation

In this mini‐review, I discuss the effects of gene interaction or epistasis from a `gene's eye view.' By a `gene's eye view' of epistasis, I mean that I will consider a single, bi‐allelic locus, A , whose effects on fitness result only from its interactions with alleles of another, unknown locus, X . I will show how changes in the frequencies of alleles at the background locus affect the relationship of alleles at the A ‐locus to fitness. Changing the genetic background changes the fundamental characteristics of the A ‐locus, such as the magnitude and sign of allelic effects on fitness, and, consequently, it changes the strength and pattern of selection. I consider each of the four kinds of pair–wise interactions between two loci and show that some kinds of epistasis are more sensitive than others to population genetic subdivision. Lastly, I show that some kinds of epistasis are more likely than others to affect the process of speciation and contribute to or be responsible for general genetic features of interspecific hybrids, such as Haldane's rule.     

ANTAGONISTIC VERSUS NONANTAGONISTIC MODELS OF BALANCING SELECTION: CHARACTERIZING THE RELATIVE TIMESCALES AND HITCHHIKING EFFECTS OF PARTIAL SELECTIVE SWEEPS

Antagonistically selected alleles‐–those with opposing fitness effects between sexes, environments, or fitness components‐–represent an important component of additive genetic variance in fitness‐related traits, with stably balanced polymorphisms often hypothesized to contribute to observed quantitative genetic variation. Balancing selection hypotheses imply that intermediate‐frequency alleles disproportionately contribute to genetic variance of life‐history traits and fitness. Such alleles may also associate with population genetic footprints of recent selection, including reduced genetic diversity and inflated linkage disequilibrium at linked, neutral sites. Here, we compare the evolutionary dynamics of different balancing selection models, and characterize the evolutionary timescale and hitchhiking effects of partial selective sweeps generated under antagonistic versus nonantagonistic (e.g., overdominant and frequency‐dependent selection) processes. We show that the evolutionary timescales of partial sweeps tend to be much longer, and hitchhiking effects are drastically weaker, under scenarios of antagonistic selection. These results predict an interesting mismatch between molecular population genetic and quantitative genetic patterns of variation. Balanced, antagonistically selected alleles are expected to contribute more to additive genetic variance for fitness than alleles maintained by classic, nonantagonistic mechanisms. Nevertheless, classical mechanisms of balancing selection are much more likely to generate strong population genetic signatures of recent balancing selection.     

PHYSIOLOGICAL AND BEHAVIORAL ADAPTATION TO VARYING ENVIRONMENTS: A MATHEMATICAL MODEL

We develop a mathematical model to explore the evolution of habitat selection and physiological adaptation in a heterogeneous environment. The model assumes the following conditions: 1) a panmictic population of infinite size; 2) prereproductive individuals mobile enough to move between patches; 3) alleles at one locus code for absence or presence of adaptation to detrimental patches; 4) alleles at a second locus code for absence or presence of behavior(s) that cause avoidance of the detrimental patches; 5) additive effects of alleles controlling physiology and behavior; 6) frequency-independent fitness. Results of the model indicate that nontrivial, polymorphic equilibria do not exist. The pattern of genotypic fitnesses and the initial allelic frequencies can influence whether the population adapts by physiological or behavioral mechanisms, or by both. Linkage between the two loci can alter the outcome of evolution, given specified genotypic fitness values and initial allelic frequencies.     

Muller''s Ratchet under Epistatic Selection

In a finite asexual population mean fitness may decrease by a process known as Muller's ratchet, which proceeds if all individuals with the minimum number of deleterious alleles are randomly lost. If these alleles have independent effects on fitness, previous analysis suggested that the rate of this decrease either remains constant or, if accumulation of mutations leads to the decline of the population size, grows. Here I show that this conclusion is quite sensitive to the assumption of independence. If deleterious alleles have synergistic fitness effects, then, as the ratchet advances, the frequency of the best available genotype will necessarily increase, making its loss less and less probable. As a result, sufficiently strong synergistic epistasis can effectively halt the action of Muller's ratchet. Instead of being driven extinct, a finite asexual population could then survive practically indefinitely, although with lower mean fitness than without random drift.     

Selection and segregation distortion in a sex-differentiated population

We extend the classical model for selection at an autosomal locus in a sex-differentiated population to include segregation distortion. The equations remain the same, but the fitness parameters are interpreted differently and refer to alleles instead of genotypes. We derive conditions for internal and external stability of the equilibria, i.e., stability with respect to perturbations of alleles that are already present at equilibrium and stability with respect to invasion attempts by newly arising alleles. We show that, in a sex-differentiated population, external stability of an equilibrium can be judged on the basis of Shaw--Mohler criteria. Throughout, we compare the situation in populations with and without sex differentiation. Interestingly, internal stability is more difficult to achieve in a population without sex differentiation than in a population in which selection and segregation distortion are restricted to one sex. In a companion paper we show how the general results of the present paper can lead to new insights into specific systems such as the t complex of the house mouse.     

Correlation of heterozygosity and the number of alleles in different frequency classes

Analytic expressions for the expectations and variance of the number of alleles with gene frequencies in a specified class are derived in the entire population as well as in a random sample of genes drawn from the population. The correlation of this quantity with heterozygosity at the locus is also obtained. The derivations are given in details for a steady state population of finite size under the infinite allele model of selectively neutral alleles. The results are extended to include weak selection pressures and non-stationarity of the population. The relevance of the correlation of heterozygosity and the number of rare alleles in connection with the neutralist-selectionist controversy is also discussed.     

Reduced fertility in patients' families is consistent with the sexual selection model of schizophrenia and schizotypy

Background

Schizophrenia is a mental disorder marked by an evolutionarily puzzling combination of high heritability, reduced reproductive success, and a remarkably stable prevalence. Recently, it has been proposed that sexual selection may be crucially involved in the evolution of schizophrenia. In the sexual selection model (SSM) of schizophrenia and schizotypy, schizophrenia represents the negative extreme of a sexually selected indicator of genetic fitness and condition. Schizotypal personality traits are hypothesized to increase the sensitivity of the fitness indicator, thus conferring mating advantages on high-fitness individuals but increasing the risk of schizophrenia in low-fitness individuals; the advantages of successful schzotypy would be mediated by enhanced courtship-related traits such as verbal creativity. Thus, schizotypy-increasing alleles would be maintained by sexual selection, and could be selectively neutral or even beneficial, at least in some populations. However, most empirical studies find that the reduction in fertility experienced by schizophrenic patients is not compensated for by increased fertility in their unaffected relatives. This finding has been interpreted as indicating strong negative selection on schizotypy-increasing alleles, and providing evidence against sexual selection on schizotypy.

Methodology

A simple mathematical model is presented, showing that reduced fertility in the families of schizophrenic patients can coexist with selective neutrality of schizotypy-increasing alleles, or even with positive selection on schizotypy in the general population. If the SSM is correct, studies of patients'' families can be expected to underestimate the true fertility associated with schizotypy.

Significance

This paper formally demonstrates that reduced fertility in the families of schizophrenic patients does not constitute evidence against sexual selection on schizotypy-increasing alleles. Futhermore, it suggests that the fertility estimates derived from extant studies may be biased to an unknown extent. These results have important implications for the evolutionary genetics of psychosis.