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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

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.

The evolution of isochores: evidence from SNP frequency distributions.

The large-scale systematic variation in nucleotide composition along mammalian and avian genomes has been a focus of the debate between neutralist and selectionist views of molecular evolution. Here we test whether the compositional variation is due to mutation bias using two new tests, which do not assume compositional equilibrium. In the first test we assume a standard population genetics model, but in the second we make no assumptions about the underlying population genetics.

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.

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.

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.

Deterministic mutation rate variation in the human genome.

Several studies of substitution rate variation have indicated that the local mutation rate varies over the mammalian genome. In the present study, we show significant variation in substitution rates within the noncoding part of the human genome using 4.7 Mb of human-chimpanzee pairwise comparisons.

The compositional evolution of the murid genome.

Murid rodents show much less variation in isochore base composition than do most other mammals, a difference which has been referred to as the murid shift. We have investigated the murid shift by asking (1) whether the murid shift is ongoing and (2) whether there is any evidence of selection or biased gene conversion affecting base composition in the present-day mouse genome. By estimating the ancestral base composition of protein-coding genes in murids we can confirm that the murid shift is ongoing.

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.

Regional similarities in polymorphism in the human genome extend over many megabases.

The human genome exhibits extensive regional variation both in base composition and in the synonymous and nonsynonymous substitution rates of protein-coding genes. If such regional variation is due to variation in mutation rates, then levels of polymorphism should also vary across the human genome. Building on recent advances in mapping the human genome, we demonstrate regional variation in single nucleotide polymorphism density extending over many megabases.