Sex and conflict.

Evolutionary conflict occurs when the deterministic spread of an allele lowers the fitness either of its bearer or of other individuals in the population, leading to selection for suppressors. Sex promotes conflict because associations between alleles are temporary. Differing selection on males and females, sexual selection, and differences in transmission patterns between classes of nuclear and cytoplasmic genes can all give rise to conflict. Inert Y chromosomes, uniparental inheritance of cytoplasmic genes, mating strains and sexes, and many features of sexual behavior may have evolved in part as a result of evolutionary conflict. Estimates of its quantitative importance, however, are still needed.

Do we understand the evolution of genomic imprinting?

The conflict theory is the only hypothesis to have attracted any critical attention for the evolution of genomic imprinting. Although the earliest data appeared supportive, recent systematic analyses have not confirmed the model's predictions. The status of theory remains undecided, however, as post-hoc explanations can be provided as to why these predictions are not borne out.

Evolutionary genetics. The silence of the genes.

A low mutation rate is required for the evolution of large genomes. But is the repression of inadvertent gene expression also important, and is the evolution of complexity limited by the efficiency of noise reduction?

Do essential genes evolve slowly?

Approximately two thirds of all knockouts of individual mouse genes give rise to viable fertile mice. These genes have thus been termed 'non-essential' in contrast to 'essential' genes, the knockouts of which result in death or infertility. Although non-essential genes are likely to be under selection that favours sequence conservation [1], it is predicted that they are less subject to such stabilising selection than essential genes, and hence evolve faster [2]. We have addressed this issue by analysing the molecular evolution of 108 non-essential and 67 essential genes that have been sequenced in both mouse and rat. On preliminary analysis, the non-essential genes appeared to be faster evolving than the essential ones. We found, however, that the non-essential class contains a disproportionate number of immune-system genes that may be under directional selection (that is, selection favouring change) because of host-parasite coevolution. After correction for this bias, we found that the rate at which genes evolve does not correlate with the severity of the knockout phenotype. This was corroborated by the finding that, whereas neuron-specific genes have significantly lower rates of change than other genes, essential and non-essential neuronal genes have comparable rates of evolution. Our findings most probably reflect strong selection acting against even very subtle deleterious phenotypes, and indicate that the putative involvement of directional selection in host-parasite coevolution and gene expression within the nervous system explains much more of the variance in rates of gene evolution than does the knockout phenotype.

A difficult phase for introns-early. Molecular evolution.

Close analysis of intron phase - the position of introns within codons - is claimed to provide novel evidence supporting the view that introns predate the divergence of bacteria and eukaryotes and, via 'exon shuffling', played a crucial role in protein evolution. But just how compelling is this evidence?

Growth effects of uniparental disomies and the conflict theory of genomic imprinting.

While numerous theories have been proposed for the evolution of genomic imprinting, few have been tested. The conflict theory proposes that imprinting is an intra-individual manifestation of classical parent-offspring conflict. This theory is unique in predicting that imprinted genes expressed from the paternally derived genome should be enhancers of pre- and post-natal growth, while those expressed from the maternally derived genome should be growth suppressors. We examine this prediction by reviewing the literature on growth of human and mouse progeny that have inherited both copies (or part thereof) of a particular chromosome from only one parent. Perhaps surprisingly, we find that much of the data do not support the hypothesis.

Sex biases in the mutation rate.

Men have more germ-line cell divisions than women. Does this lead to a higher mutation rate in males? Most estimates of the proportion of mutations originating in men come either from direct observation of disease-inducing mutations or from analysis of the relative rate of evolution of sex-linked and autosomal genes in primates. The latter mode of analysis has also been applied to other mammals, birds and files. For unknown reasons, this method produces contradictory results. A majority of estimates using the best direct methods in humans indicate a male bias for point mutations, but the variance in estimates is high. It is unclear how the evolutionary and direct data correspond and a consensus as to the extent of any male bias is not presently possible. While the number of germ-line cell divisions might contribute to differences, this by no means accounts for all of the data.

Evidence for purifying selection acting on silent sites in BRCA1.

In mammals, it is usually assumed that selection cannot be strong enough to act on nucleotide mutations that do not cause a change at the protein level (i.e. 'silent' or 'synonymous' mutations). Here we report the results of a molecular evolutionary analysis of BRCA1. We find a repeatable pronounced peak in the ratio of nonsynonymous to synonymous substitutions between codons 200-300. Unusually, this peak is caused by a plummet in the silent-site rate of evolution. The most parsimonious interpretation of these data is that purifying selection is acting on silent sites.

Causes of sex ratio bias may account for unisexual sterility in hybrids: a new explanation of Haldane's rule and related phenomena.

Unisexual hybrid disruption can be accounted for by interactions between sex ratio distorters which have diverged in the species of the hybrid cross. One class of unisexual hybrid disruption is described by Haldane's rule, namely that the sex which is absent, inviable or sterile is the heterogametic sex. This effect is mainly due to incompatibility between X and Y chromosomes. We propose that this incompatibility is due to a mutual imbalance between meiotic drive genes, which are more likely to evolve on sex chromosomes than autosomes. The incidences of taxa with sex chromosome drive closely matches those where Haldane's rule applies: Aves, Mammalia, Lepidoptera and Diptera. We predict that Haldane's rule is not universal but is correct for taxa with sex chromosome meiotic drive. A second class of hybrid disruption affects the male of the species regardless of which sex is heterogametic. Typically the genes responsible for this form of disruption are cytoplasmic. These instances are accounted for by the release from suppression of cytoplasmic sex ratio distorters when in a novel nuclear cytotype. Due to the exclusively maternal transmission of cytoplasm, cytoplasmic sex ratio distorters cause only female-biased sex ratios. This asymmetry explains why hybrid disruption is limited to the male.

The causes of synonymous rate variation in the rodent genome. Can substitution rates be used to estimate the sex bias in mutation rate?

Miyata et al. have suggested that the male-to-female mutation rate ratio (alpha) can be estimated by comparing the neutral substitution rates of X-linked (X), Y-linked (Y), and autosomal (A) genes. Rodent silent site X/A comparisons provide very different estimates from X/Y comparisons. We examine three explanations for this discrepancy: (1) statistical biases and artifacts, (2) nonneutral evolution, and (3) differences in mutation rate per germline replication. By estimating errors and using a variety of methodologies, we tentatively reject explanation 1. Our analyses of patterns of codon usage, synonymous rates, and nonsynonymous rates suggest that silent sites in rodents are evolving neutrally, and we can therefore reject explanation 2. We find both base composition and methylation differences between the different sets of chromosomes, a result consistent with explanation 3, but these differences do not appear to explain the observed discrepancies in estimates of alpha. Our finding of significantly low synonymous substitution rates in genomically imprinted genes suggests a link between hemizygous expression and an adaptive reduction in the mutation rate, which is consistent with explanation 3. Therefore our results provide circumstantial evidence in favor of the hypothesis that the discrepancies in estimates of alpha are due to differences in the mutation rate per germline replication between different parts of the genome. This explanation violates a critical assumption of the method of Miyata et al., and hence we suggest that estimates of alpha, obtained using this method, need to be treated with caution.

Molecular evolution of an imprinted gene: repeatability of patterns of evolution within the mammalian insulin-like growth factor type II receptor.

The repeatability of patterns of variation in Ka/Ks and Ks is expected if such patterns are the result of deterministic forces. We have contrasted the molecular evolution of the mammalian insulin-like growth factor type II receptor (Igf2r) in the mouse-rat comparison with that in the human-cow comparison. In so doing, we investigate explanations for both the evolution of genomic imprinting and for Ks variation (and hence putatively for mutation rate evolution). Previous analysis of Igf2r, in the mouse-rat comparison, found Ka/Ks patterns that were suggested to be contrary to those expected under the conflict theory of imprinting. We find that Ka/Ks variation is repeatable and hence confirm these patterns. However, we also find that the molecular evolution of Igf2r signal sequences suggests that positive selection, and hence conflict, may be affecting this region. The variation in Ks across Igf2r is also repeatable. To the best of our knowledge this is the first demonstration of such repeatability. We consider three explanations for the variation in Ks across the gene: (1) that it is the result of mutational biases, (2) that it is the result of selection on the mutation rate, and (3) that it is the product of selection on codon usage. Explanations 2 and 3 predict a Ka-Ks correlation, which is not found. Explanation 3 also predicts a negative correlation between codon bias and Ks, which is also not found. However, in support of explanation 1 we do find that in rodents the rate of silent C --> T mutations at CpG sites does covary with Ks, suggesting that methylation-induced mutational patterns can explain some of the variation in Ks. We find evidence to suggest that this CpG effect is due to both variation in CpG density, and to variation in the frequency with which CpGs mutate. Interestingly, however, a GC4 analysis shows no covariance with Ks, suggesting that to eliminate methyl-associated effects CpG rates themselves must be analyzed. These results suggest that, in contrast to previous studies of intragenic variation, Ks patterns are not simply caused by the same forces responsible for Ka/Ks correlations.

Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection.

In numerous species, from bacteria to Drosophila, evidence suggests that selection acts even on synonymous codon usage: codon bias is greater in more abundantly expressed genes, the rate of synonymous evolution is lower in genes with greater codon bias, and there is consistency between genes in the same species in which codons are preferred. In contrast, in mammals, while nonequal use of alternative codons is observed, the bias is attributed to the background variance in nucleotide concentrations, reflected in the similar nucleotide composition of flanking noncoding and exonic third sites. However, a systematic examination of the covariants of codon usage controlling for background nucleotide content has yet to be performed. Here we present a new method to measure codon bias that corrects for background nucleotide content and apply this to 2396 human genes. Nearly all (99%) exhibit a higher amount of codon bias than expected by chance. The patterns associated with selectively driven codon bias are weakly recovered: Broadly expressed genes have a higher level of bias than do tissue-specific genes, the bias is higher for genes with lower rates of synonymous substitutions, and certain codons are repeatedly preferred. However, while these patterns are suggestive, the first two patterns appear to be methodological artifacts. The last pattern reflects in part biases in usage of nucleotide pairs. We conclude that we find no evidence for selection on codon usage in humans.

The effect of tandem substitutions on the correlation between synonymous and nonsynonymous rates in rodents.

Nonsynonymous substitutions in DNA cause amino acid substitutions while synonymous substitutions in DNA leave amino acids unchanged. The cause of the correlation between the substitution rates at nonsynonymous (K(A)) and synonymous (K(S)) sites in mammals is a contentious issue, and one that impacts on many aspects of molecular evolution. Here we use a large set of orthologous mammalian genes to investigate the causes of the K(A)-K(S) correlation in rodents. The strength of the K(A)-K(S) correlation exceeds the neutral theory expectation when substitution rates are estimated using algorithmic methods, but not when substitution rates are estimated by maximum likelihood. Irrespective of this methodological uncertainty the strength of the K(A)-K(S) correlation appears mostly due to tandem substitutions, an excess of which is generated by substitutional nonindependence. Doublet mutations cannot explain the excess of tandem synonymous-nonsynonymous substitutions, and substitution patterns indicate that selection on silent sites is the likely cause. We find no evidence for selection on codon usage. The nature of the relationship between synonymous divergence and base composition is unclear because we find a significant correlation if we use maximum-likelihood methods but not if we use algorithmic methods. Finally, we find that K(S) is reduced at the start of genes, which suggests that selection for RNA structure may affect silent sites in mammalian protein-coding genes.

Covariation of GC content and the silent site substitution rate in rodents: implications for methodology and for the evolution of isochores.

Many attempts to test selectionist and neutralist models employ estimates of synonymous (Ks) and non-synonymous (Ka) substitution rates of orthologous genes. For example, a stronger Ka-Ks correlation than expected under neutrality has been argued to indicate a role for selection and the absence of a Ks-GC4 correlation has been argued to be inconsistent with neutral models for isochore evolution. However, both of these results, we have shown previously, are sensitive to the method by which Ka and Ks are estimated. Using a maximum likelihood (ML) estimator (GY94) we found a positive correlation between Ks and GC4 and only a weak correlation between Ka and Ks, lower than expected under neutral expectations. This ML method is computationally slow. Recently, a new ad hoc approximation of this ML method has been provided (YN00). This is effectively an extension of Li's protocol but that also allows for codon usage bias. This method is computationally near-instantaneous and therefore potentially of great utility for analysis of large datasets. Here we ask whether this method might have such applicability. To this end we ask whether it too recovers the two unusual results. We report that when the ML and earlier ad hoc methods disagree, YN00 recovers the results described by the ML methods, i.e. a positive correlation between GC4 and Ks and only a weak correlation between Ks and Ka. If the ML method can be trusted, then YN00 can also be considered an adequately reliable method for analysis of large datasets. Assuming this to be so we also analyze further the patterns. We show, for example, that the positive correlation between GC4 and Ks is probably in part a mutational bias, there being more methyl induced CpG-->TpG mutations in GC rich regions. As regards the evolution of isochores, it seems inappropriate to use the claimed lack of a correlation between GC and Ks as definitive evidence either against or for any model. If the positive correlation is real then, we argue, this is hard to reconcile with the biased gene conversion model for isochore formation as this predicts a negative correlation.

The molecular evolution of signal peptides.

Signal peptides direct mature peptides to their appropriate cellular location, after which they are cleaved off. Very many random alternatives can serve the same function. Of all coding sequences, therefore, signal peptides might come closest to being neutrally evolving. Here we consider this issue by examining the molecular evolution of 76 mouse-rat orthologues, each with defined signal peptides. Although they do evolve rapidly, they evolve about half as fast as neutral sequences. This indicates that a substantial proportion of mutations must be under stabilizing selection. A few putative signal sequences lack a hydrophobic core and these tend to be more slowly evolving than others, indicating even stronger stabilizing selection. However, closer scrutiny suggests that some of these represent mis-annotations in GenBank. It is also likely that some of the substitutions are not neutral. We find, for example, that the rate of protein evolution correlates with that of the mature peptide. This may be a result of compensatory evolution. We also find that signal peptides of immune genes tend to be faster evolving than the average, which suggests an association with antagonistic co-evolution. Previous reports also indicated that the signal peptide of the imprinted gene, Igf2r, is also unusually fast evolving. This, it was hypothesized, might also be indicative of antagonistic co-evolution. Comparison of Igf2r's signal peptide evolution shows that, although it is not an outlier, its rate of evolution is comparable to that of many of the faster evolving immune system signal sequences and 5/6 of the amino acid changes do not conserve hydrophobicity. This is at least suggestive that there is something unusual about Igf2r's signal sequence.

A model for the mechanism of transmission ratio distortion and for t-associated hybrid sterility.

A mechanistic model is presented to account for the action of t-complex of mice. This model takes account of recent evidence suggesting that t-complex distorters are amorphs or hypomorphs. Following Lyon's (Genet. Res. 59, 27 (1992) scheme, the model proposes that the t-complex distorter (tcd+) loci for normal function than does the wild-type form of tcr. However, a tradeoff against this ability to drive is a reduced efficiency of the haploid specific product of tcrt in the absence of drive. Regulation of tcr could be achieved by differential splicing or post-translational modification under the control of the t-complex distorters. It is shown that the model is consistent with known fertility and distortion data, as well as with the finding that the mechanism of drive is intimately connected with the mechanism of intraspecific homozygous sterility. Importantly, the model predicts that the mechanism of hybrid sterility associated with the t-complex is the same as the mechanism of intraspecific homozygous sterility. If accepted then this will be, to the best of the author's knowledge, the first description and characterization of a Haldane rule sterility gene. The new understanding of the mechanisms of t-complex shows its mode of operation to be fundamentally different to the only other well-described autosomal meiotic driver, Segregation Distorter (SD) of Drosophila melanogaster.

Molecular evolution of imprinted genes: no evidence for antagonistic coevolution.

Genomically imprinted genes are those for which expression is dependent on the sex of the parent from which they are derived. Numerous theories have been proposed for the evolution of genomic imprinting: one theory is that it is an intra-individual manifestation of classical parent -offspring conflict. This theory is unique in predicting that an arms race may develop between maternally and paternally derived genes for the control of foetal growth demands. Such antagonistic coevolution may be mediated through changes in the structure of the proteins concerned. Comparable coevolution is the most likely explanation for the rapid changes seen in antigenic components of parasites and antigen recognition components of immune systems. We have examined the evolution of insulin-like growth factor Igf2, and its antagonistic receptor Igf2r) and find that in contrast to immune genes, at the sites of mutual binding they are highly conserved. In addition, we have analysed the rate of molecular evolution of seven imprinted genes including Igf2 and Igf2r), sequenced in both mouse and rat, and had that this is the same as that of nonimprinted receptors and significantly lower than that of immune genes controlling for differences in mutation rates. Contrary to the expectations of the conflict hypothesis, we hence find no evidence for antagonistic coevolution of imprinted genes mediated by changes in sequence.