Mesquite: a modular system for evolutionary analysis. Version 2.6

A new statistical method for estimating divergence dates of species from DNA sequence data by a molecular clock approach is developed. This method takes into account effectively the information contained in a set of DNA sequence data. The molecular clock of mitochondrial DNA (mtDNA) was calibrated by setting the date of divergence between primates and ungulates at the Cretaceous-Tertiary boundary (65 million years ago), when the extinction of dinosaurs occurred. A generalized least-squares method was applied in fitting a model to mtDNA sequence data, and the clock gave dates of 92.3 +/- 11.7, 13.3 +/- 1.5, 10.9 +/- 1.2, 3.7 +/- 0.6, and 2.7 +/- 0.6 million years ago (where the second of each pair of numbers is the standard deviation) for the separation of mouse, gibbon, orangutan, gorilla, and chimpanzee, respectively, from the line leading to humans. Although there is some uncertainty in the clock, this dating may pose a problem for the widely believed hypothesis that the pipedal creature Australopithecus afarensis, which lived some 3.7 million years ago at Laetoli in Tanzania and at Hadar in Ethiopia, was ancestral to man and evolved after the human-ape splitting. Another likelier possibility is that mtDNA was transferred through hybridization between a proto-human and a proto-chimpanzee after the former had developed bipedalism.

Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.

Two simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions are presented. Although they give no weights to different types of codon substitutions, these methods give essentially the same results as those obtained by Miyata and Yasunaga's and by Li et al.'s methods. Computer simulation indicates that estimates of synonymous substitutions obtained by the two methods are quite accurate unless the number of nucleotide substitutions per site is very large. It is shown that all available methods tend to give an underestimate of the number of nonsynonymous substitutions when the number is large.

/pdf/simple-methods-for-estimating-the-numbers-of-synonymous-and-32dmft0e36.pdf

Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions.

A new measure of the extent of population subdivision as inferred from allele frequencies at microsatellite loci is proposed and tested with computer simulations. This measure, called R(ST), is analogous to Wright's F(ST) in representing the proportion of variation between populations. It differs in taking explicit account of the mutation process at microsatellite loci, for which a generalized stepwise mutation model appears appropriate. Simulations of subdivided populations were carried out to test the performance of R(ST) and F(ST). It was found that, under the generalized stepwise mutation model, R(ST) provides relatively unbiased estimates of migration rates and times of population divergence while F(ST) tends to show too much population similarity, particularly when migration rates are low or divergence times are long [corrected].

https://www.genetics.org/content/genetics/139/1/457.full.pdf

A measure of population subdivision based on microsatellite allele frequencies.

With the aim of analyzing and interpreting data on DNA polymorphism obtained by DNA sequencing or restriction enzyme technique, a mathematical theory on the expected evolutionary relationship among DNA sequences (nucleons) sampled is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. The statistical property of the number of nucleotide differences between randomly chosen nucleons and that of heterozygosity or nucleon diversity is investigated using this theory. These studies indicate that the estimates of the average number of nucleotide differences and nucleon diversity have a large variance, and a large part of this variance is due to stochastic factors. Therefore, increasing sample size does not help reduce the variance significantly The distribution of sample allele (nucleomorph) frequencies is also studied, and it is shown that a small number of samples are sufficient in order to know the distribution pattern.

https://www.genetics.org/content/genetics/105/2/437.full.pdf

Evolutionary relationship of dna sequences in finite populations

We analyzed the complete mitochondrial DNA (mtDNA) sequences of three humans (African, European, and Japanese), three African apes (common and pygmy chimpanzees, and gorilla), and one orangutan in an attempt to estimate most accurately the substitution rates and divergence times of hominoid mtDNAs. Nonsynonymous substitutions and substitutions in RNA genes have accumulated with an approximately clock-like regularity. From these substitutions and under the assumption that the orangutan and African apes diverged 13 million years ago, we obtained a divergence time for humans and chimpanzees of 4.9 million years. This divergence time permitted calibration of the synonymous substitution rate (3.89 x 10(-8)/site per year). To obtain the substitution rate in the displacement (D)-loop region, we compared the three human mtDNAs and measured the relative abundance of substitutions in the D-loop region and at synonymous sites. The estimated substitution rate in the D-loop region was 7.00 x 10(-8)/site per year. Using both synonymous and D-loop substitutions, we inferred the age of the last common ancestor of the human mtDNAs as 143,000 +/- 18,000 years. The shallow ancestry of human mtDNAs, together with the observation that the African sequence is the most diverged among humans, strongly supports the recent African origin of modern humans, Homo sapiens sapiens.

https://www.pnas.org/content/pnas/92/2/532.full.pdf

Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs.

In the last few years, more than 500 primate major histocompatibility complex (Mhc) genes or parts thereof have been sequenced. The extraordinary sequence information is used here to draw conclusions about the manner of Mhc evolution. The Mhc genes are found to evolve at a relatively slow rate with the regularity of a clock. It takes from 1 to 6 million years for a new mutation to be incorporated into an Mhc allele, and the mutation rate is comparable to that of most other primate genes. The nonsynonymous sites coding for the peptide-binding region (PBR) are under relatively weak positive selection pressure (selection coefficient of a few percent only); the nonsynonymous non-PBR sites are under moderate negative selection pressure. The positive pressure is probably provided by parasites and is responsible for the trans-species persistence of allelic lineages at functional Mhc loci for more than 40 million years.

The Molecular Descent of the Major Histocompatibility Complex

A genealogical relationship among genes at a locus (gene tree) sampled from three related populations was examined with special reference to population relatedness (population tree). A phylogenetically informative event in a gene tree constructed from nucleotide differences consists of interspecific coalescences of genes in each of which two genes sampled from different populations are descended from a common ancestor. The consistency probability between gene and population trees in which they are topologically identical was formulated in terms of interspecific coalescences. It was found that the consistency probability thus derived substantially increases as the sample size of genes increases, unless the divergence time of populations is very long compared to population sizes. Hence, there are cases where large samples at a locus are very useful in inferring a population tree.

https://www.genetics.org/content/genetics/122/4/957.full.pdf

Gene Genealogy in Three Related Populations: Consistency Probability between Gene and Population Trees

Genetic variation at most loci examined in human populations indicates that the (effective) population size has been approximately 10(4) for the past 1 Myr and that individuals have been genetically united rather tightly. Also suggested is that the population size has never dropped to a few individuals, even in a single generation. These impose important requirements for the hypotheses for the origin of modern humans: a relatively large population size and frequent migration if populations were geographically subdivided. Any hypothesis that assumes a small number of founding individuals throughout the late Pleistocene can be rejected. Extraordinary polymorphism at some loci of the major histocompatibility complex (Mhc) rules out the past action of severe bottlenecks, or the so-called founder principle, which invokes only a small number of founding individuals when a new species emerges. This conclusion may be extended to the 35-Myr-old history of the human lineage, because some polymorphism at Mhc loci seems to have lasted that long. Furthermore, although the population structure prior to the late Pleistocene is less clear, owing to the insensitivity of Mhc alleles, even to low levels of migration, the nature of Mhc polymorphism suggests that the effective size of populations leading to humans was as large as 10(5). Hence, the effective population size of humans might have become somewhat smaller in most of the late Pleistocene. The reduction could be due either to the then adverse environment in the Old World and/or to the increased migration rate. It is also argued that population explosion fostered by the agriculture revolution has had significant effects on incorporating new alleles into human populations.

Allelic genealogy and human evolution.

Different alleles undergoing strong symmetric balancing selection show a simple genealogical structure (allelic genealogy), similar to the gene genealogy described by the coalescence process for a sample of neutral genes randomly drawn from a panmictic population at equilibrium. The only difference between the two genealogies lies in the different time scales. An approximate scaling factor for allelic genealogy relative to that of neutral gene genealogy is [square root of S/(2M)].[In[S/(16 pi M2)]]-3/2, where M = Nu and S = 2Ns (N, effective population size; u, mutation rate to selected alleles per locus per generation; s, selection coefficient). The larger the value of square root of S/M (greater than or equal to 100), the larger the scaling factor. These findings, supported by simulation results, allow one to apply the theoretical results of the coalescence process directly to the allelic genealogy. Combined with the trans-species evolution of the major histocompatibility complex polymorphism for which balancing selection is believed to be responsible, allelic genealogy predicts that the number of breeding individuals in the human population could not be as small as 50-100 at any time of its evolutionary history. The analysis appears to contradict the founder principle as being important in recent mammalian evolution.

https://www.pnas.org/content/pnas/87/7/2419.full.pdf

Naoyuki Takahata

Papers

Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs.

The Molecular Descent of the Major Histocompatibility Complex

Gene Genealogy in Three Related Populations: Consistency Probability between Gene and Population Trees

Allelic genealogy and human evolution.

A simple genealogical structure of strongly balanced allelic lines and trans-species evolution of polymorphism.