Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA

More than an order of magnitude difference in substitution rate exists among sites within hypervariable region 1 of the control region of human mitochondrial DNA. A two-rate Poisson mixture and a negative binomial distribution are used to describe the distribution of the inferred number of changes per nucleotide site in this region. When three data sets are pooled, however, the two-rate model cannot explain the data. The negative binomial distribution always fits, suggesting that substitution rates are approximately gamma distributed among sites. Simulations presented here provide support for the use of a biased, yet commonly employed, method of examining rate variation. The use of parsimony in the method to infer the number of changes at each site introduces systematic errors into the analysis. These errors preclude an unbiased quantification of variation in substitution rate but make the method conservative overall. The method can be used to distinguish sites with highly elevated rates, and 29 such sites are identified in hypervariable region 1. Variation does not appear to be clustered within this region. Simulations show that biases in rates of substitution among nucleotides and non-uniform base composition can mimic the effects of variation in rate among sites. However, these factors contribute little to the levels of rate variation observed in hypervariable region 1.

[1]  C. I. Bliss,et al.  FITTING THE NEGATIVE BINOMIAL DISTRIBUTION TO BIOLOGICAL DATA AND NOTE ON THE EFFICIENT FITTING OF THE NEGATIVE BINOMIAL , 1953 .

[2]  S. Benzer,et al.  ON THE TOPOGRAPHY OF THE GENETIC FINE STRUCTURE. , 1961, Proceedings of the National Academy of Sciences of the United States of America.

[3]  SOME ELEMENTARY TESTS FOR MIXTURES OF DISCRETE DISTRIBUTION , 1963 .

[4]  A. Cohen.  Estimation in mixtures of poisson and mixtures of exponential distributions , 1965 .

[5]  L. Pauling,et al.  Evolutionary Divergence and Convergence in Proteins , 1965 .

[6]  J. L. King,et al.  Non-Darwinian evolution. , 1969, Science.

[7]  T. Jukes.  CHAPTER 24 – Evolution of Protein Molecules , 1969 .

[8]  Thomas Uzzell,et al.  Fitting Discrete Probability Distributions to Evolutionary Events , 1971, Science.

[9]  W. Fitch.  Toward Defining the Course of Evolution: Minimum Change for a Specific Tree Topology , 1971 .

[10]  T. Jukes,et al.  Estimation of evolutionary changes in certain homologous polypeptide chains. , 1972, Journal of molecular biology.

[11]  M. Kimura.  Model of effectively neutral mutations in which selective constraint is incorporated. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[12]  W. Brown,et al.  Rapid evolution of animal mitochondrial DNA. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[13]  G. Attardi,et al.  Mapping of nascent light and heavy strand transcripts on the physical map of HeLa cell mitochondrial DNA. , 1980, Nucleic acids research.

[14]  C. Reynolds,et al.  Bloom-forming cyanobacterium Microcystis aeruginosa overwinters on sediment surface , 1980, Nature.

[15]  Theoretical foundations for quantitative paleogenetics , 1980, Journal of molecular evolution.

[16]  M. Walberg,et al.  Sequence and properties of the human KB cell and mouse L cell D-loop regions of mitochondrial DNA. , 1981, Nucleic acids research.

[17]  F. Sanger,et al.  Sequence and organization of the human mitochondrial genome , 1981, Nature.

[18]  J. Cairns,et al.  Random components in mutagenesis , 1982, Nature.

[19]  W. Hauswirth,et al.  Nucleotide sequence evidence for rapid genotypic shifts in the bovine mitochondrial DNA D-loop , 1983, Nature.

[20]  R. Hudson.  Properties of a neutral allele model with intragenic recombination. , 1983, Theoretical population biology.

[21]  Richard R. Hudson,et al.  TESTING THE CONSTANT‐RATE NEUTRAL ALLELE MODEL WITH PROTEIN SEQUENCE DATA , 1983, Evolution; international journal of organic evolution.

[22]  B. Greenberg,et al.  Intraspecific nucleotide sequence variability surrounding the origin of replication in human mitochondrial DNA. , 1983, Gene.

[23]  C. Aquadro,et al.  Human mitochondrial DNA variation and evolution: analysis of nucleotide sequences from seven individuals. , 1983, Genetics.

[24]  G. B. Golding,et al.  Estimates of DNA and protein sequence divergence: an examination of some assumptions. , 1983, Molecular biology and evolution.

[25]  F. Tajima.  Evolutionary relationship of DNA sequences in finite populations. , 1983, Genetics.

[26]  D. Chang,et al.  Precise identification of individual promoters for transcription of each strand of human mitochondrial DNA , 1984, Cell.

[27]  K. Risko,et al.  An analysis of the dynamics of mammalian mitochondrial DNA sequence evolution. , 1984, Molecular biology and evolution.

[28]  M. Stoneking,et al.  Mitochondrial DNA and two perspectives on evolutionary genetics , 1985 .

[29]  J. Hixson,et al.  Initiation of transcription from each of the two human mitochondrial promoters requires unique nucleotides at the transcriptional start sites. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Evidence for local DNA influences on patterns of substitutions in the human alpha-interferon gene family. , 1986, Canadian journal of genetics and cytology. Journal canadien de genetique et de cytologie.

[31]  C. Saccone,et al.  Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA. , 1986, Journal of molecular biology.

[32]  B. Mignotte,et al.  Sequence deduced physical properties in the D-loop region common to five vertebrate mitochondrial DNAs. , 1987, Journal of theoretical biology.

[33]  J. Wrigglesworth,et al.  Achievements and perspectives of mitochondrial research , 1987 .

[34]  N. Saitou,et al.  The neighbor-joining method: a new method for reconstructing phylogenetic trees. , 1987, Molecular biology and evolution.

[35]  J. N. Topper,et al.  Promoter selection in human mitochondria involves binding of a transcription factor to orientation-independent upstream regulatory elements , 1987, Cell.

[36]  James E. Hixson,et al.  Comparisons of ape and human sequences that regulate mitochondrial DNA transcription and D-loop DNA synthesis , 1988, Nucleic Acids Res..

[37]  M. King,et al.  50 – DETERMINATION OF FAMILIAL RELATIONSHIPS , 1989 .

[38]  A. Wilson,et al.  Patterns of ribosomal RNA evolution in salamanders. , 1989, Molecular biology and evolution.

[39]  S. Horai,et al.  Intraspecific nucleotide sequence differences in the major noncoding region of human mitochondrial DNA. , 1990, American journal of human genetics.

[40]  J. Hein.  Reconstructing evolution of sequences subject to recombination using parsimony. , 1990, Mathematical biosciences.

[41]  A. Wilson,et al.  Sequence Evolution of Mitochondrial DNA in Humans and Chimpanzees: Control Region and a Protein-Coding Region , 1991 .

[42]  大沢 省三,et al.  Evolution of life : fossils, molecules, and culture , 1991 .

[43]  A. Di Rienzo,et al.  Branching pattern in the evolutionary tree for human mitochondrial DNA. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[44]  K. Hawkes,et al.  African populations and the evolution of human mitochondrial DNA. , 1991, Science.

[45]  M. Stoneking,et al.  Human origins and analysis of mitochondrial DNA sequences. , 1992, Science.

[46]  M. Stoneking.  Mitochondrial DNA and human evolution , 1994, Journal of bioenergetics and biomembranes.

[47]  E. Markowitz.  Estimation and testing goodness-of-fit for some models of codon fixation variability , 1970, Biochemical Genetics.

[48]  Walter M. Fitch,et al.  A method for estimating the number of invariant amino acid coding positions in a gene using cytochrome c as a model case , 1967, Biochemical Genetics.

[49]  Minimum Distance and Discrepancy Approach for Estimation and Testing Goodness of Fit , 2004 .

[50]  W. Fitch,et al.  An improved method for determining codon variability in a gene and its application to the rate of fixation of mutations in evolution , 1970, Biochemical Genetics.

[51]  T. Jukes.  Evolutionary pattern of specificity regions in light chains of immunoglobulins , 2004, Biochemical Genetics.

[52]  H. Kishino,et al.  Dating of the human-ape splitting by a molecular clock of mitochondrial DNA , 2005, Journal of Molecular Evolution.

[53]  Walter M. Fitch,et al.  The molecular evolution of cytochrome c in eukaryotes , 1976, Journal of Molecular Evolution.

[54]  S. Palumbi,et al.  Rates of molecular evolution and the fraction of nucleotide positions free to vary , 1989, Journal of Molecular Evolution.

[55]  M. Hasegawa,et al.  Time of the deepest root for polymorphism in human mitochondrial DNA , 2005, Journal of Molecular Evolution.

[56]  Richard E. Dickerson,et al.  The structure of cytochromec and the rates of molecular evolution , 2005, Journal of Molecular Evolution.

[57]  Morris Goodman,et al.  The spatial distribution of fixed mutations within genes coding for proteins , 1983, Journal of Molecular Evolution.

[58]  J. Otsuka,et al.  Accumulation pattern of amino acid substitutions in protein evolution , 2005, Journal of Molecular Evolution.

[59]  G. Pesole,et al.  The main regulatory region of mammalian mitochondrial DNA: Structure-function model and evolutionary pattern , 1991, Journal of Molecular Evolution.