Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

The mutation rate and cancer

The selection of advantageous mutations underlies tumorigenesis. The growth of a tumor is therefore a form of evolution at the somatic level, in which the population is comprised of individual cells within the tumor. Models of tumorigenesis have considered the relative importance of mutation and selection. We show that selection is more important than an increased mutation rate in the growth of a tumor. Some cancers may acquire a “mutator phenotype,” probably leading to faster growth, but mutator phenotypes are not necessary for carcinogenesis.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

The high spontaneous mutation rate: Is it a health risk?*

The human mutation rate for base substitutions is much higher in males than in females and increases with paternal age. This effect is mainly, if not entirely, due to the large number of cell divisions in the male germ line. The mutation-rate increase is considerably greater than expected if the mutation rate were simply proportional to the number of cell divisions. In contrast, those mutations that are small deletions or rearrangements do not show the paternal age effect. The observed increase with the age of the father in the incidence of children with different dominant mutations is variable, presumably the result of different mixtures of base substitutions and deletions. In Drosophila, the rate of mutations causing minor deleterious effects is estimated to be about one new mutation per zygote. Because of a larger number of genes and a much larger amount of DNA, the human rate is presumably higher. Recently, the Drosophila data have been reanalyzed and the mutation-rate estimate questioned, but I believe that the totality of evidence supports the original conclusion. The most reasonable way in which a species can cope with a high mutation rate is by quasi-truncation selection, whereby a number of mutant genes are eliminated by one “genetic death.”

Multilocus DNA fingerprinting reveals high rate of heritable genetic mutation in herring gulls nesting in an industrialized urban site.

Genotoxins, such as polycyclic aromatic compounds, are ubiquitous in urban and industrial environments. Our understanding of the role that these chemicals play in generating DNA sequence mutations is predominantly derived from laboratory studies with specific genotoxins or extracts of contaminants from environmental media. Most assays are not indicative of the germinal effects of exposure in situ to complex mixtures of common environmental mutagens. Using multilocus DNA fingerprinting, we found the mutation rate in herring gulls inhabiting a heavily industrialized urban harbor (Hamilton Harbour, Ontario) to be more than twice as high as three rural sites: Kent Island, Bay of Fundy; Chantry Island, Lake Huron; and Presqu'ile Provincial Park in Lake Ontario. Overall we found a mutation rate of 0.017 +/- 0.004 per offspring band in Hamilton, 0.006 +/- 0.002 at Kent Island, 0.002 +/- 0.002 from Chantry Island, and 0.004 +/- 0.002 from Presqu'ile Provincial Park. The mutation rate from the rural sites (pooled) was significantly lower than the rate observed in Hamilton Harbour (Fisher's exact test, two-tailed; P = 0.0006). These minisatellite DNA mutations may be important biomarkers for heritable genetic changes resulting from in situ exposure to environmental genotoxins in a free-living vertebrate species.

Mutation rates in the dihydrofolate reductase gene of Plasmodium falciparum

A new method has been established to define the limits on a spontaneous mutation rate for a gene in Plasmodium falciparum. The method combines mathematical modelling and large-scale in vitro culturing and calculates the difference in mutant frequencies at 2 separate time-points. We measured the mutation rate at 2 positions in the dihydrofolate reductase (DHFR) gene of 3D7, a pyrimethamine-sensitive line of P. fulciparum. This line was re-cloned and an effectively large population was treated with a selective pyrimethamine concentration of 40 nM. We detected point mutations at codon-46 (TTA to TCA) and codon-108 (ACC to AAC), resulting in serine replacing leucine and asparagine replacing serine respectively in the corresponding gene product. The substitutions caused a decrease in pyrimethamine sensitivity. By mathematical modelling we determined that the mutation rate at a given position in DHFR was low and occurred at less than 2(.)5 x 10(-9) mutations/DHFR gene/replication. This result has important implications for Plasmodium genetic diversity and antimalarial drug therapy by demonstrating that even with lon mutation rates anti-malarial resistance will inevitably arise when mutant alleles are selected under drug pressure.

Adaptive mutations in bacteria: high rate and small effects

Evolution by natural selection is driven by the continuous generation of adaptive mutations. We measured the genomic mutation rate that generates beneficial mutations and their effects on fitness in Escherichia coli under conditions in which the effect of competition between lineages carrying different beneficial mutations is minimized. We found a rate on the order of 10–5 per genome per generation, which is 1000 times as high as previous estimates, and a mean selective advantage of 1%. Such a high rate of adaptive evolution has implications for the evolution of antibiotic resistance and pathogenicity.

Stochastic demography and the neutral substitution rate in class-structured populations.

The neutral rate of allelic substitution is analyzed for a class-structured population subject to a stationary stochastic demographic process. The substitution rate is shown to be generally equal to the effective mutation rate, and under overlapping generations it can be expressed as the effective mutation rate in newborns when measured in units of average generation time. With uniform mutation rate across classes the substitution rate reduces to the mutation rate.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

Male mutation bias and possible long-term effects of human activities.

The ability of a population to adapt to changing environments depends critically on the amount and kind of genetic variability it possesses. Mutations are an important source of new genetic variability and may lead to new adaptations, especially if the population size is large. Mutation rates are extremely variable between and within species, and males usually have higher mutation rates as a result of elevated rates of male germ cell division. This male bias affects the overall mutation rate. We examined the factors that influence male mutation bias, and focused on the effects of classical life-history parameters, such as the average age at reproduction and elevated rates of sperm production in response to sexual selection and sperm competition. We argue that human-induced changes in age at reproduction or in sexual selection will affect male mutation biases and hence overall mutation rates. Depending on the effective population size, these changes are likely to influence the long-term persistence of a population.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Epstein-Barr virus (EBV) nuclear antigen (EBNA)-4 mutation in EBV-associated malignancies in three different populations

Different ethnic groups with a high human leukocyte antigen (HLA)-A11 prevalence have been shown to experience a high rate of Epstein-Barr virus (EBV) infection, EBV-associated malignancies, and Epstein-Barr nuclear antigen (EBNA)-4 mutations. The epitopes 393-408 and 416-424 of EBNA-4 are major antigenic epitopes that elicit an HLA-A11 cytotoxic T lymphocyte (CTL) response to EBV infection. Mutations selectively involving one or more nucleotide residues in these epitopes affect the antigenicity of EBNA-4, because the mutant EBV strains are not recognized by the HLA-A11-restricted CTLs. To investigate these mutations in common EBV-associated malignancies occurring in different populations, we studied the mutation rate of epitopes 393-408 and 416-424 of EBNA-4 in 25 cases of EBV-associated Hodgkin's disease (HD), nine cases of AIDS-related non-Hodgkin's lymphoma, and 37 cases of EBV-associated gastric carcinoma (GC) from the United States, Brazil, and Japan. We found one or more mutations in these two epitopes in 50% (6/12) of United States HD, 15% (2/13) of Brazilian HD, 50% (6/12) United States GC and 28% (7/25) Japanese GC, and 22% (2/9) of United States AIDS-lymphoma. Similar mutations were found in 30% (3/10) of United States reactive, 0% (0/6) of Brazilian reactive, and 25% (2/8) Japanese reactive tissues. The most frequent amino acid substitutions were virtually identical to those seen in previously reported isolates from EBV-associated nasopharyngeal carcinomas and Burkitt's lymphomas occurring in high prevalence HLA-A11 regions. However, only 2/28 (7%) mutations occurred in HLA-A11-positive patients. Our studies suggest that: 1) EBNA-4 mutations are a common phenomenon in EBV-associated HD, GC, and AIDS-lymphoma; 2) the mutation rate does not vary in these geographic areas and ethnic groups; 3) EBNA-4 mutations in EBV-associated United States and Brazilian HD, United States and Japanese GC, and United States AIDS lymphomas are not related to patients' HLA-A11 status.

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^