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 (alpha(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 alpha(m) and c. Alpha(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 alpha(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 alpha(m), which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F(13,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 alpha(m).

Testing for adaptive evolution of the female reproductive protein ZPC in mammals, birds and fishes reveals problems with the M7-M8 likelihood ratio test.

BACKGROUND: Adaptive evolution appears to be a common feature of reproductive proteins across a very wide range of organisms. A promising way of addressing the evolutionary forces responsible for this general phenomenon is to test for adaptive evolution in the same gene but among groups of species, which differ in their reproductive biology. One can then test evolutionary hypotheses by asking whether the variation in adaptive evolution is consistent with the variation in reproductive biology. We have attempted to apply this approach to the study of a female reproductive protein, zona pellucida C (ZPC), which has been previously shown by the use of likelihood ratio tests (LRTs) to be under positive selection in mammals. RESULTS: We tested for evidence of adaptive evolution of ZPC in 15 mammalian species, in 11 avian species and in six fish species using three different LRTs (M1a-M2a, M7-M8, and M8a-M8). The only significant findings of adaptive evolution came from the M7-M8 test in mammals and fishes. Since LRTs of adaptive evolution may yield false positives in some situations, we examined the properties of the LRTs by several different simulation methods. When we simulated data to test the robustness of the LRTs, we found that the pattern of evolution in ZPC generates an excess of false positives for the M7-M8 LRT but not for the M1a-M2a or M8a-M8 LRTs. This bias is strong enough to have generated the significant M7-M8 results for mammals and fishes. CONCLUSION: We conclude that there is no strong evidence for adaptive evolution of ZPC in any of the vertebrate groups we studied, and that the M7-M8 LRT can be biased towards false inference of adaptive evolution by certain patterns of non-adaptive evolution.

A novel method with improved power to detect recombination hotspots from polymorphism data reveals multiple hotspots in human genes.

We introduce a new method for detection of recombination hotspots from population genetic data. This method is based on (a) defining an (approximate) penalized likelihood for how recombination rate varies with physical position and (b) maximizing this penalized likelihood over possible sets of recombination hotspots. Simulation results suggest that this is a more powerful method for detection of hotspots than are existing methods. We apply the method to data from 89 genes sequenced in African American and European American populations. We find many genes with multiple hotspots, and some hotspots show evidence of being population-specific. Our results suggest that hotspots are randomly positioned within genes and could be as frequent as one per 30 kb.

Analysis of recombination in Campylobacter jejuni from MLST population data.

We analyze recombination in C. jejuni using MLST data from isolates taken from wild birds, cattle, wild rabbits, and water in a 100-km2 study region in Cheshire, UK. We use a recent approximate likelihood method for inference, based on combining likelihood information from all pairs of segregating (polymorphic) sites in the data. We find substantial evidence for recombination, but only for recombination with short tract lengths, of around 225-750 bp. We estimate that the rate of recombination is of a similar magnitude to the rate of mutation.

A comparison of three estimators of the population-scaled recombination rate: accuracy and robustness.

We have performed simulations to assess the performance of three population genetics approximate-likelihood methods in estimating the population-scaled recombination rate from sequence data. We measured performance in two ways: accuracy when the sequence data were simulated according to the (simplistic) standard model underlying the methods and robustness to violations of many different aspects of the standard model. Although we found some differences between the methods, performance tended to be similar for all three methods. Despite the fact that the methods are not robust to violations of the underlying model, our simulations indicate that patterns of relative recombination rates should be inferred reasonably well even if the standard model does not hold. In addition, we assess various techniques for improving the performance of approximate-likelihood methods. In particular we find that the composite-likelihood method of Hudson (2001) can be improved by including log-likelihood contributions only for pairs of sites that are separated by some prespecified distance.

Male-driven biased gene conversion governs the evolution of base composition in human alu repeats.

Regional biases in substitution pattern are likely to be responsible for the large-scale variation in base composition observed in vertebrate genomes. However, the evolutionary forces responsible for these biases are still not clearly defined. In order to study the processes of mutation and fixation across the entire human genome, we analyzed patterns of substitution in Alu repeats since their insertion. We also studied patterns of human polymorphism within the repeats. There is a highly significant effect of recombination rate on the pattern of substitution, whereas no such effect is seen on the pattern of polymorphism. These results suggest that regional biases in substitution are caused by biased gene conversion, a process that increases the probability of fixation of mutations that increase GC content. Furthermore, the strongest correlate of substitution patterns is found to be male recombination rates rather than female or sex-averaged recombination rates. This indicates that in addition to sexual dimorphism in recombination rates, the sexes also differ in the relative rates of crossover and gene conversion.

Comparison of the chicken and turkey genomes reveals a higher rate of nucleotide divergence on microchromosomes than macrochromosomes.

A distinctive feature of the avian genome is the large heterogeneity in the size of chromosomes, which are usually classified into a small number of macrochromosomes and numerous microchromosomes. These chromosome classes show characteristic differences in a number of interrelated features that could potentially affect the rate of sequence evolution, such as GC content, gene density, and recombination rate. We studied the effects of these factors by analyzing patterns of nucleotide substitution in two sets of chicken-turkey sequence alignments. First, in a set of 67 orthologous introns, divergence was significantly higher in microchromosomes (chromosomes 11-38; 11.7% divergence) than in both macrochromosomes (chromosomes 1-5; 9.9% divergence; P = 0.016) and intermediate-sized chromosomes (chromosomes 6-10; 9.5% divergence; P = 0.026). At least part of this difference was due to the higher incidence of CpG sites on microchromosomes. Second, using 155 orthologous coding sequences we noted a similar pattern, in which synonymous substitution rates on microchromosomes (13.1%) were significantly higher than were rates on macrochromosomes (10.3%; P = 0.024). Broadly assuming neutrality of introns and synonymous sites, or constraints on such sequences do not differ between chromosomal classes, these observations imply that microchromosomal genes are exposed to more germ line mutations than those on other chromosomes. We also find that dN/dS ratios for genes located on microchromosomes (average, 0.094) are significantly lower than those of macrochromosomes (average, 0.185; P = 0.025), suggesting that the proteins of genes on microchromosomes are under greater evolutionary constraint.

Evidence for turnover of functional noncoding DNA in mammalian genome evolution.

The vast majority of the mammalian genome does not code for proteins, and a fundamental question in genomics is: What proportion of the noncoding mammalian genome is functional? Most attempts to address this issue use sequence comparisons between highly diverged mammals such as human and mouse to identify conservation due to negative selection. But such comparisons will underestimate the true proportion of functional noncoding DNA if there is turnover, if patterns of negative selection change over time. Here we test whether the inferred level of negative selection differs between different pairwise species comparisons. Using a multiple alignment of more than a megabase of contiguous sequence from eight mammalian species, we find a strong negative relationship between inferred levels of negative selection and pairwise divergence using 21 pairwise comparisons. This result suggests that there is a high rate of turnover of functional noncoding elements in the mammalian genome, so measures of functional constraint based on human-mouse comparisons may seriously underestimate the true value.

Gene expression, synteny, and local similarity in human noncoding mutation rates.

The human genome is organized with regard to many features such as isochores, Giemsa bands, clusters of genes with similar expression patterns, and contiguous regions with shared evolutionary histories (synteny blocks). In addition to these genomic features, it is clear that mutation rates also vary across the human genome. To address how mutation rates and genomic features are related, we analyzed substitution rates at three classes of putatively neutral noncoding sites (nongenic, intronic, and ancestral repeats) in approximately 14 Mb of human-chimpanzee alignments covering human chromosome 7. Patterns of mutation rate variation inferred from substitution rate variation differ among the three site classes. In particular, we find that intronic mutation rates are strongly affected by the breadth of expression of the genes in which they reside, with broadly expressed genes exhibiting low mutation rates, probably as a consequence of the transcription-coupled repair process acting in the germ line. All site classes show significant local similarities in mutation rate at the megabase scale, and regional similarities in nongenic mutation rate covary with blocks of synteny between the human and mouse genomes, indicating that the evolutionary history of a genomic region is an important determinant of mutation rate.

Male-biased mutation rate and divergence in autosomal, z-linked and w-linked introns of chicken and Turkey.

To investigate mutation-rate variation between autosomes and sex chromosomes in the avian genome, we have analyzed divergence between chicken (Gallus gallus) and turkey (Meleagris galopavo) sequences from 33 autosomal, 28 Z-linked, and 14 W-linked introns with a total ungapped alignment length of approximately 43,000 bp. There are pronounced differences in the mean divergence among autosomes and sex chromosomes (autosomes [A] = 10.08%, Z chromosome = 10.99%, and W chromosome = 5.74%), and we use these data to estimate the male-to-female mutation-rate ratio (alpha(m)) from Z/A, Z/W, and A/W comparisons at 1.71, 2.37, and 2.52, respectively. Because the alpha(m) estimates of the three comparisons do not differ significantly, we find no statistical support for a specific reduction in the Z chromosome mutation rate (Z reduction estimated at 4.89%, P = 0.286). The idea of mutation-rate reduction in the sex chromosome hemizygous in one sex (i.e., X in mammals, Z in birds) has been suggested on the basis of theory on adaptive mutation-rate evolution. If it exists in birds, the effect would, thus, seem to be weak; a preliminary power analysis suggests that it is significantly less than 18%. Because divergence may vary within chromosomal classes as a result of variation in mutation and/or selection, we developed a novel double-bootstrapping method, bootstrapping both by introns and sites from concatenated alignments, to estimate confidence intervals for chromosomal class rates and for alpha(m). The narrowest interval for the alpha(m) estimate is 1.88 to 2.97 from the Z/W comparison. We also estimated alpha(m) using maximum likelihood on data from all three chromosome classes; this method yielded alpha(m) = 2.47 and approximate 95% confidence intervals of 2.27 to 2.68. Our data are broadly consistent with the idea that mutation-rate differences between chromosomal classes can be explained by the male mutation bias alone.

Do avian mitochondria recombine?

The dogma of strict maternal inheritance of mitochondria is now being tested with population genetics methods on sequence data from many species. In this study we investigated whether recombination occurs in the mitochondria of the blue tit ( Parus caeruleus) by studying polymorphisms in the mitochondrial control region and in a recently identified (A)(n) microsatellite on the W chromosome. The female heterogamety of avian sex chromosomes allows a test of whether mitochondrial recombination affects genealogical inference by comparison of mitochondrial and W-linked sequence variation. There is no discrepancy between mitochondrial and W-linked genealogies in blue tits, consistent with no recombination. We also analyzed mitochondrial sequence variation in both blue tits and peregrine falcons ( Falco peregrinus) using a coalescent-based approach which accounts for recurrent mutation; in neither bird species did we find evidence of recombination. We conclude that it is unlikely that mitochondrial recombination has large effects on mitochondrial genetic variability in birds.

Fixation biases affecting human SNPs.

Under neutrality all classes of mutation have an equal probability of becoming fixed in a population. In this article, we describe our analysis of the frequency distributions of >5000 human SNPs and provide evident of biases in the process of fixation of certain classes of point mutation that are most likely to be attributable to biased gene conversion. The results indicate an increased fixation probability of mutations that result in the incorporation of a GC base pair. Furthermore, in transcribed regions this process exhibits strand asymmetry, and is biased towards preserving a G base on the coding strand. Biased gene conversion has the potential to explain both existence of isochores and the compositional asymmetry in mammalian transcribed regions.

Life history and the male mutation bias.

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 (alpha(m)). 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 alpha(m), which was estimated by comparing Z-linked and W-linked substitution rates in gametologous introns. Alpha(m) 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 alpha(m) 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 alpha(m) between clades were due to differences in W-linked or Z-linked substitution rates. In one case, a significant difference in alpha(m) 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.

Mutation rate variation in the mammalian genome.

Recent advances in the large-scale sequencing of mammalian genomes have provided a means to study divergence in not only genic sequences but also in the non-coding bulk of DNA. There is evidence of significant variation in the levels of divergence between presumably neutral regions, pointing at an underlying variation in the rate of mutation across the genome. Apparently, such variation occurs on different scales, including sequence context effects (the influence of neighboring nucleotides on the rate of mutation at individual sites), variation within chromosomes (on the scales of kilobases as well as megabases), and between chromosomes (among autosomes as well as between autosomes and sex chromosomes). An important aspect for further research in this area is to study whether there is an ultimate evolutionary explanation for mutation rate variation within mammalian genomes.

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.

Compositional evolution of noncoding DNA in the human and chimpanzee genomes.

We have examined the compositional evolution of noncoding DNA in the primate genome by comparison of lineage-specific substitutions observed in 1.8 Mb of genomic alignments of human, chimpanzee, and baboon with 6542 human single-nucleotide polymorphisms (SNPs) rooted using chimpanzee sequence. The pattern of compositional evolution, measured in terms of the numbers of GC-->AT and AT-->GC changes, differs significantly between fixed and polymorphic sites, and indicates that there is a bias toward fixation of AT-->GC mutations, which could result from weak directional selection or biased gene conversion in favor of high GC content. Comparison of the frequency distributions of a subset of the SNPs revealed no significant difference between GC-->AT and AT-->GC polymorphisms, although AT-->GC polymorphisms in regions of high GC segregate at slightly higher frequencies on average than GC-->AT polymorphisms, which is consistent with a fixation bias favoring high GC in these regions. However, the substitution data suggest that this fixation bias is relatively weak, because the compositional structure of the human and chimpanzee genomes is becoming homogenized, with regions of high GC decreasing in GC content and regions of low GC increasing in GC content. The rate and pattern of nucleotide substitution in 333 Alu repeats within the human-chimpanzee-baboon alignments are not significantly affected by the GC content of the region in which they are inserted, providing further evidence that, since the time of the human-chimpanzee ancestor, there has been little or no regional variation in mutation bias.

Partitioning the variation in mammalian substitution rates.

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.

A low rate of simultaneous double-nucleotide mutations in primates.

The occurrence of double-nucleotide (doublet) mutations is contrary to the normal assumption that point mutations affect single nucleotides. Here we develop a new method for estimating the doublet mutation rate and apply it to more than a megabase of human-chimpanzee-baboon genomic DNA alignments and more than a million human single-nucleotide polymorphisms. The new method accounts for the effect of regional variation in evolutionary rates, which may be a confounding factor in previous estimates of the doublet mutation rate. Furthermore we determine sequence context effects by using sequence comparisons over a variety of lineage lengths. This approach yields a new estimate of the doublet mutation rate of 0.3% of the singleton rate, indicating that doublet mutations are far rarer than previously thought. Our results suggest that doublet mutations are unlikely to have caused the correlation between synonymous and nonsynonymous substitution rates in mammals, and also show that regional variation and sequence context effects play an important role in primate DNA sequence evolution.

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.

Quantifying the slightly deleterious mutation model of molecular evolution.

We have attempted to quantify the frequency and effects of slightly deleterious mutations (SDMs), those that have selective effects close to the reciprocal of the effective population size of a species, by comparing the level of selective constraint in protein-coding genes of related species that have different present-day effective population sizes. In our two comparisons, the species with the smaller effective population size showed lower constraint, implying that SDMs had become fixed. The fixation of SDMs was supported by the observation of a higher fraction of radical to conservative amino acid substitutions in species with smaller effective population sizes. The fraction of strongly deleterious mutations (which rarely become fixed) is >70% in most species. Only approximately 10% or fewer of mutations seem to behave as SDMs, but SDMs could comprise a substantial fraction of mutations in protein-coding genes that have a chance of becoming fixed between species.