We present here a framework for the study of molecular variation within a single species. Information on DNA haplotype divergence is incorporated into an analysis of variance format, derived from a matrix of squared-distances among all pairs of haplotypes. This analysis of molecular variance (AMOVA) produces estimates of variance components and F-statistic analogs, designated here as phi-statistics, reflecting the correlation of haplotypic diversity at different levels of hierarchical subdivision. The method is flexible enough to accommodate several alternative input matrices, corresponding to different types of molecular data, as well as different types of evolutionary assumptions, without modifying the basic structure of the analysis. The significance of the variance components and phi-statistics is tested using a permutational approach, eliminating the normality assumption that is conventional for analysis of variance but inappropriate for molecular data. Application of AMOVA to human mitochondrial DNA haplotype data shows that population subdivisions are better resolved when some measure of molecular differences among haplotypes is introduced into the analysis. At the intraspecific level, however, the additional information provided by knowing the exact phylogenetic relations among haplotypes or by a nonlinear translation of restriction-site change into nucleotide diversity does not significantly modify the inferred population genetic structure. Monte Carlo studies show that site sampling does not fundamentally affect the significance of the molecular variance components. The AMOVA treatment is easily extended in several different directions and it constitutes a coherent and flexible framework for the statistical analysis of molecular data.

/pdf/analysis-of-molecular-variance-inferred-from-metric-89hl32sg7j.pdf

Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data.

Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

PAML, currently in version 4, is a package of programs for phylogenetic analyses of DNA and protein sequences using maximum likelihood (ML). The programs may be used to compare and test phylogenetic trees, but their main strengths lie in the rich repertoire of evolutionary models implemented, which can be used to estimate parameters in models of sequence evolution and to test interesting biological hypotheses. Uses of the programs include estimation of synonymous and nonsynonymous rates (d(N) and d(S)) between two protein-coding DNA sequences, inference of positive Darwinian selection through phylogenetic comparison of protein-coding genes, reconstruction of ancestral genes and proteins for molecular restoration studies of extinct life forms, combined analysis of heterogeneous data sets from multiple gene loci, and estimation of species divergence times incorporating uncertainties in fossil calibrations. This note discusses some of the major applications of the package, which includes example data sets to demonstrate their use. The package is written in ANSI C, and runs under Windows, Mac OSX, and UNIX systems. It is available at -- (http://abacus.gene.ucl.ac.uk/software/paml.html).

/pdf/paml-4-phylogenetic-analysis-by-maximum-likelihood-2gg19rsgos.pdf

PAML 4: Phylogenetic Analysis by Maximum Likelihood

Reconstructing phylogenies from intraspecific data (such as human mitochondrial DNA variation) is often a challenging task because of large sample sizes and small genetic distances between individuals. The resulting multitude of plausible trees is best expressed by a network which displays alternative potential evolutionary paths in the form of cycles. We present a method ("median joining" [MJ]) for constructing networks from recombination-free population data that combines features of Kruskal's algorithm for finding minimum spanning trees by favoring short connections, and Farris's maximum-parsimony (MP) heuristic algorithm, which sequentially adds new vertices called "median vectors", except that our MJ method does not resolve ties. The MJ method is hence closely related to the earlier approach of Foulds, Hendy, and Penny for estimating MP trees but can be adjusted to the level of homoplasy by setting a parameter epsilon. Unlike our earlier reduced median (RM) network method, MJ is applicable to multistate characters (e.g., amino acid sequences). An additional feature is the speed of the implemented algorithm: a sample of 800 worldwide mtDNA hypervariable segment I sequences requires less than 3 h on a Pentium 120 PC. The MJ method is demonstrated on a Tibetan mitochondrial DNA RFLP data set.

/pdf/median-joining-networks-for-inferring-intraspecific-19tuocf48l.pdf

Median-joining networks for inferring intraspecific phylogenies.

Mitochondrial DNA was purified from four species of higher primates (Guinea baboon, rhesus macaque, guenon, and human) and digested with 11 restriction endonucleases. A cleavage map was constructed for the mitochondrial DNA of each species. Comparison of the maps, aligned with respect to the origin and direction of DNA replication, revealed that the species differ from one another at most of the cleavage sites. The degree of divergence in nucleotide sequence at these sites was calculated from the fraction of cleavage sites shared by each pair of species. By plotting the degree of divergence in mitochondrial DNA against time of divergence, the rate of base substitution could be calculated from the initial slope of the curve. The value obtained, 0.02 substitutions per base pair per million years, was compared with the value for single-copy nuclear DNA. The rate of evolution of the mitochondrial genome appears to exceed that of the single-copy fraction of the nuclear genome by a factor of about 10. This high rate may be due, in part, to an elevated rate of mutation in mitochondrial DNA. Because of the high rate of evolution, mitochondrial DNA is likely to be an extremely useful molecule to employ for high-resolution analysis of the evolutionary process.

https://www.pnas.org/content/pnas/76/4/1967.full.pdf

Rapid evolution of animal mitochondrial DNA.

In the past decade, mitochondrial DNA (mtDNA) analysis has become established as a powerful tool for evolutionary studies of animals. These studies, recently reviewed by Avise (8) and Wilson et al (149), have used mtDNA analyses to provide insights into population structure and gene flow, hybridization, biogeography, and phylogenetic relationships. Major advances have been made in our understanding of the molecular biology of animal mtDNA (reviewed in 7, 28). A powerful synergism exists between the two fields of research: Evolutionary studies provide comparative data on mtDNA organization and function, and molecular investigations can, and should, improve the level of sophistication of evolutionary studies that use mtDNA. Here we focus on molecular aspects of animal mtDNA that are especially relevant to its use in evolutionary studies. Thus, we consider the form and frequency of three types of change in mtDNA: base substitution, length variation, and sequence rearrangement. Attention is also given to the nature and possible consequences of interactions between nuclear and mitochondrial (mt) genes. In the concluding sections, we discuss the way in which the knowledge of molecular processes allows informed use of mtDNA variation in evolutionary studies. Only animal mtDNAs are considered. The molecular biology and evolution of chloroplast DNA and nonmetazoan mtDNAs have been reviewed elsewhere (41, 66, 111, 121, 123).

https://www.jstor.org/stable/pdfplus/2097133.pdf

Evolution of animal mitochondrial dna: relevance for population biology and systematics

We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.

https://deepblue.lib.umich.edu/bitstream/handle/2027.42/48036/239_2005_Article_BF01734101.pdf?sequence=1&isAllowed=y

Mitochondrial DNA sequences of primates: tempo and mode of evolution.

The evolutionary relationships among the four major lineages of arthropods remain controversial, despite extensive study. We report here a derived gene rearrangement common to insects and crustaceans but absent in the other arthropod groups. This finding strongly supports an insect-crustacean evolutionary lineage that is separate from those leading to myriapods and chelicerates.

/pdf/gene-translocation-links-insects-and-crustaceans-4em96e5hfn.pdf

Wesley M. Brown

Papers

Rapid evolution of animal mitochondrial DNA.

Evolution of animal mitochondrial dna: relevance for population biology and systematics

Mitochondrial DNA sequences of primates: tempo and mode of evolution.

Gene translocation links insects and crustaceans

Big trees from little genomes: mitochondrial gene order as a phylogenetic tool.