Nick G.C. Smith

Adaptive protein evolution in Drosophila

For over 30 years a central question in molecular evolution has been whether natural selection plays a substantial role in evolution at the DNA sequence level. Evidence has accumulated over the last decade that adaptive evolution does occur at the protein level, but it has remained unclear how prevalent adaptive evolution is. Here we present a simple method by which the number of adaptive substitutions can be estimated and apply it to data from Drosophila simulans and D. yakuba. We estimate that 45% of all amino-acid substitutions have been fixed by natural selection, and that on average one adaptive substitution occurs every 45 years in these species.

Microsatellite evolution inferred from human– chimpanzee genomic sequence alignments

Most studies of microsatellite evolution utilize long, highly mutable loci, which are unrepresentative of the majority of simple repeats in the human genome. Here we use an unbiased sample of 2,467 microsatellite loci derived from alignments of 5.1 Mb of genomic sequence from human and chimpanzee to investigate the mutation process of tandemly repetitive DNA. The results indicate that the process of microsatellite evolution is highly heterogeneous, exhibiting differences between loci of different lengths and motif sizes and between species. We find a highly significant tendency for human dinucleotide repeats to be longer than their orthologues in chimpanzees, whereas the opposite trend is observed in mononucleotide repeat arrays. Furthermore, the rate of divergence between orthologues is significantly higher at longer loci, which also show significantly greater mutability per repeat number. These observations have important consequences for understanding the molecular mechanisms of microsatellite mutation and for the development of improved measures of genetic distance.

Human disease genes: patterns and predictions

We compared genes at which mutations are known to cause human disease (disease genes) with other human genes (nondisease genes) using a large set of human-rodent alignments to infer evolutionary patterns. Such comparisons may be of use both in predicting disease genes and in understanding the general evolution of human genes. Four features were found to differ significantly between disease and nondisease genes, with disease genes (i) evolving with higher nonsynonymous/synonymous substitution rate ratios (Ka/Ks), (ii) evolving at higher synonymous substitution rates, (iii) with longer protein-coding sequences, and (iv) expressed in a narrower range of tissues. Discriminant analysis showed that these differences may help to predict human disease genes. We also investigated other factors affecting the mode of evolution in the disease genes: Ka/Ks is significantly affected by protein function, mode of inheritance, and the reduction of life expectancy caused by disease.

LIFE HISTORY AND THE MALE MUTATION BIAS

Abstract If DNA replication is a major cause of mutation, then those life‐history characters, which are expected to affect the number of male germline cell divisions, should also affect the male to female mutation bias (am). We tested this hypothesis by comparing several clades of bird species, which show variation both in suitable life‐history characters (generation time as measured by age at first breeding and sexual selection as measured by frequency of extrapair paternity) and in am, which was estimated by comparing Z‐linked and W‐linked substitution rates in gametologous introns. am differences between clades were found to positively covary with both generation time and sexual selection, as expected if DNA replication causes mutation. The effects of extrapair paternity frequency on am suggests that increased levels of sexual selection cause higher mutation rates, which offers an interesting solution to the paradox of the loss of genetic variance associated with strong directional sexual selection. We also used relative rate tests to examine whether the observed differences in am between clades were due to differences in W‐linked or Z‐linked substitution rates. In one case, a significant difference in am between two clades was shown to be due to W‐linked rates and not Z‐linked rates, a result that suggests that mutation rates are not determined by replication alone.

The evolutionary dynamics of male‐killers and their hosts

Male-killing bacteria are cytoplasmic sex-ratio distorters that are transmitted vertically through females of their insect hosts. The killing of male hosts by their bacteria is thought to be an adaptive bacterial trait because it augments the fitness of female hosts carrying clonal relatives of those bacteria. Here we attempt to explain observations of multiple male-killers in natural host populations. First we show that such male-killer polymorphism cannot be explained by a classical model of male-killing. We then show that more complicated models incorporating the evolution of resistance in hosts can explain male-killer polymorphism. However, this is only likely if resistance genes are very costly. We also consider the long-term evolutionary dynamics of male-killers, and show that evolution towards progressively more ‘efficient’ male-killers can be thwarted by the appearance of host resistance. The presence of a resistance gene can allow a less efficient male-killer to outcompete its rival and hence reverse the trend towards more efficient transmission and reduced metabolic load on the host.

Vertebrate genome evolution: a slow shuffle or a big bang?

In vertebrates it is often found that if one considers a group of genes clustered on a certain chromosome, then the homologues of those genes often form another cluster on a different chromosome. There are four explanations, not necessarily mutually exclusive, to explain how such homologous clusters appeared. Homologous clusters are expected at a low probability even if genes are distributed at random. The duplication of a subset of the genome might create homologous clusters, as would a duplication of the entire genome. Alternatively, it may be adaptive for certain combinations of genes to cluster, although clearly the genes must have duplicated prior to rearrangement into clusters. Molecular phylogenetics provides a means to examine the origins of homologous clusters, although it is difficult to discriminate between the different explanations using current data. However, with more extensive sequencing and mapping of vertebrate genomes, especially those of the early diverging chordates, it should soon become possible to resolve the origins of homologous clusters.

Analysis of the phylogenetic distribution of isochores in vertebrates and a test of the thermal stability hypothesis

Warm-blooded vertebrates show large-scale variation in G + C content along their chromosomes, a pattern which appears to be largely absent from cold-blooded vertebrates. However, compositional variation in poikilotherms has generally been studied by ultracentrifugation rather than sequence analysis. In this paper, we investigate the compositional properties of coding sequences from a broad range of vertebrate poikilotherms using DNA sequence analysis. We find that on average poikilotherms have lower third-codon position GC contents (GC3) than homeotherms but that some poikilotherms have higher mean GC3 values. We find that most poikilotherms have lower variation in GC3 than homeotherms but that there is a correlation between GC12 and GC3 for some species, indicating that there is systematic variation in base composition across their genomes. We also demonstrate that the GC3 of genes in the zebrafish, Danio rerio, is correlated with that in humans, suggesting that vertebrates share a basic isochore structure. However, we find no correlation between either the mean GC3 or the standard deviation in GC3 and body temperature.

Are Radical and Conservative Substitution Rates Useful Statistics in Molecular Evolution

A DNA mutation in a protein coding gene which causes an amino acid change can be classified as “conservative” or “radical” depending on the magnitude of the physicochemical difference between the two amino acids: radical mutations involve larger changes than conservative mutations. Here, I examine two key issues in determining whether radical and conservative substitution rates are useful statistics in molecular evolution. The first issue is whether such rates can be estimated reliably, and for this purpose I demonstrate considerable improvements achieved by simple modifications to an existing method. The second issue is whether conservative and radical substitution rates can tell us something about selection on protein function. I address this problem by estimating positive and negative selection on conservative and radical mutations using polymorphism and divergence data from Drosophila. These analyses show that negative selection, but not positive selection, differs significantly between conservative and radical mutations. The power of conservative and radical substitution rates in testing the nearly neutral theory of molecular evolution is illustrated by the analysis of two mammalian datasets.

Small introns tend to occur in GC-rich regions in some but not all vertebrates

We thank two anonymous referees for their comments and L. Duret for access to unpublished data.

Substitution Rate Heterogeneity and the Male Mutation Bias

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (αm) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of αm and c. αm is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of αm based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G15 = 37.6, p = 0.0010) as well as within W chromosome introns (G15 = 44.0, p = 0.0001). This heterogeneity also affected the estimates of αm, which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F13,14 =2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of αm.