Molecular Phylogenetic Inference of the Woolly Mammoth Mammuthus primigenius, Based on Complete Sequences of Mitochondrial Cytochrome b and 12S Ribosomal RNA Genes

Complete sequences of cytochrome b (1,137 bases) and 12S ribosomal RNA (961 bases) genes in mitochondrial DNA were successfully determined from the woolly mammoth (Mammuthus primigenius), African elephant (Loxodonta africana), and Asian elephant (Elephas maximus). From these sequence data, phylogenetic relationships among three genera were examined. Molecular phylogenetic trees reconstructed by the neighbor-joining and the maximum parsimony methods provided an identical topology both for cytochrome b and 12S rRNA genes. These results support the ``Mammuthus-Loxodonta' clade, which is contrary to some previous morphological reports that Mammuthus is more closely related to Elephas than to Loxodonta. Received: 8 April 1997 / Accepted: 23 July 1997     

Molecular and morphological evidence on the phylogeny of the Elephantidae

The African and Asian elephants and the mammoth diverged ca. 4-6 million years ago and their phylogenetic relationship has been controversial. Morphological studies have suggested a mammoth Asian elephant relationship, while molecular studies have produced conflicting results. We obtained cytochrome b sequences of up to 545 base pairs from five mammoths, 14 Asian and eight African elephants. A high degree of polymorphism is detected within species. With a dugong sequence used as the outgroup, parsimony and maximum-likelihood analyses support a mammoth African elephant clade. As the dugong is a very distant outgroup, we employ likelihood analysis to root the tree with a molecular clock, and use bootstrap and Bayesian analyses to quantify the relative support for different topologies. The analyses support the mammoth African elephant relationship, although other trees cannot be rejected. Ancestral polymorphisms may have resulted in gene trees differing from the species phylogeny Examination of morphological data, especially from primitive fossil members, indicates that some supposed synapomorphies between the mammoth and Asian elephant are variable, others convergent or autapomorphous. A mammoth African elephant relationship is not excluded. Our results highlight the need, in both morphological and molecular phylogenetics, for multiple markers and close attention to within-taxon variation and outgroup selection.     

Sequential Loss of Two Neighboring Exons of the Tropoelastin Gene During Primate Evolution

Previous evidence has demonstrated the absence of exons 34 and 35 within the 3′ end of the human tropoelastin (ELN) gene. These exons encode conserved polypeptide domains within tropoelastin and are found in the ELN gene in vertebrate species ranging from chickens to rats to cows. We have analyzed the ELN gene in a variety of primate species to determine whether the absence of exons 34 and 35 in humans either is due to allelic variation within the human population or is a general characteristic of the Primates order. An analysis of the 3′ end of the ELN gene in several nonhuman primates and in 546 chromosomes from humans of varying ethnic background demonstrated a sequential loss of exons 34 and 35 during primate evolution. The loss of exon 35 occurred at least 35–45 million years ago, when Catarrhines diverged from Platyrrhines (New World monkeys). Exon 34 loss, in contrast, occurred only about 6–8 million years ago, when Homo separated from the common ancestor shared with chimpanzees and gorillas. Loss of both exons was probably facilitated by Alu-mediated recombination events and possibly conferred a functional evolutionary advantage in elastic tissue. Received: 6 July 1998 / Accepted: 18 February 1999     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996     

The complete mitochondrial DNA (mtDNA) of the donkey and mtDNA comparisons among four closely related mammalian species-pairs

The nucleotide sequence of the complete mitochondrial genome of the donkey, Equus asinus, was determined. The length of the molecule is 16,670 bp. The length, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of two types of repetitive motifs in the control region. The sequence of the repeats is (a) 5′-CACACCCA and (b) 5′-TGCGCGCA, respectively. The order of (a) and (b) can be expressed as {n[2(a)+(b)]+m(a)}. In 32 different clones analyzed the number of n and m ranged from 0 to 9 and 1 to 7. The two rRNA genes, the 13 peptide-coding genes, and the 22 tRNA genes of the donkey and the horse, Equus caballus, were compared in detail. Total nucleotide difference outside the control region was 6.9%. Nucleotide difference between peptide-coding genes ranged from 6.4% to 9.4% with a mean of 8.0%. In the inferred protein sequences of the 13 peptide-coding genes the amino acid difference was 0.2–8.8%, and the mean for the 13 concatenated amino acid sequences was 1.9%. In the 22 tRNA genes, the mean difference was 3.5%, and that in the two rRNA genes was 4.1%. The mtDNA differences between the donkey and the horse suggest that the evolutionary separation of the two species occurred ≈9 million years ago. Analyses of differences among the mtDNAs of three other species-pairs, harbor seal/grey seal, fin whale/blue whale, and Homo/common chimpanzee, showed that the relative evolutionary rate of individual peptide-coding genes varies among different species-pairs and modes of comparison. The findings show that the superimposition of sequence data of one lineage for resolving and dating evolutionary divergences of other lineages should be performed with caution unless based on comprehensive data. Received: 15 October 1995 / Accepted: 15 April 1996     

Molecular Clock Calibrations and Metazoan Divergence Dates

It has recently been argued that living metazoans diverged over 800 million years ago, based on evidence from 22 nuclear genes for such a deep divergence between vertebrates and arthropods (Gu 1998). Two ``internal' calibration points were used. However, only one fossil divergence date (the mammal–bird split) was directly used to calibrate the molecular clock. The second calibration point (the primate–rodent split) was based on molecular estimates that were ultimately also calibrated by the same mammal–bird split. However, the first tetrapods that can be assigned with confidence to either the mammal (synapsid) lineage or the bird (diapsid) lineage are approximately 288 million years old, while the first mammals that can be assigned with confidence to either the primate or the rodent lineages are 65 million years old, or 85 million years old if ferungulates are part of the primate lineage and zhelestids are accepted as ferungulate relatives. Recalibration of the protein data using these fossil dates indicates that metazoans diverged between 791 and 528 million years ago, a result broadly consistent with the palaeontological documentation of the ``Cambrian explosion.' The third, ``external' calibration point (the metazoan–fungal divergence) was similarly problematic, since it was based on a controversial molecular study (which in turn used fossil dates including the mammal–bird split); direct use of fossils for this calibration point gives the absurd dating of 455 million years for metazoan divergences. Similar calibration problems affect another recent study (Wang et al. 1999), which proposes divergences for metazoans of 1000 million years or more: recalibrations of their clock again yields much more recent dates, some consistent with a ``Cambrian explosion' scenario. Molecular clock studies have persuasively argued for the imperfection of the fossil record but have rarely acknowledged that their inferences are also directly based on this same record. Received: 26 January 1999 / Accepted: 14 April 1999     

Phylogenetic Relationships Within the Genus Equus and the Evolution of α and θ Globin Genes

Sequences of the α1, α2 and θ globin genes from six equid species have been determined to investigate relationships within the genus Equus. Analyses using standard phylogenetic methods, or an approach designed to account for the effects of gene conversion between the α genes, gave broadly similar results and show that the horses diverged from the zebra/ass ancestor ∼2.4 million years ago and that the zebra and ass species arose in a rapid radiation ∼0.9 million years ago. These results from the α genes are corroborated by θ gene data and are in contrast to mitochondrial DNA studies of the phylogeny of this genus, which suggest a more gradual set of speciation events. Received: 22 April 1997 / Accepted: 20 July 1998     

Phylogeographic Patterns and Evolution of the Mitochondrial DNA Control Region in Two Neotropical Cats (Mammalia, Felidae)

The ocelot (Leopardus pardalis) and margay (L. wiedii) are sister-species of Neotropical cats which evolved from a lineage that migrated into South America during the formation of the Panamanian land bridge 3–5 million years ago. Patterns of population genetic divergence of each species were studied by phylogenetic analyses of mitochondrial DNA (mtDNA) control region sequences in individuals sampled across the distribution of these taxa. Abundant genetic diversity and remarkably concordant phylogeographic partitions for both species were observed, identifying parallel geographic regions which likely reflect historical faunal barriers. Inferred aspects of phylogeography, population genetic structure, and demographic history were used to formulate conservation recommendations for these species. In addition, observed patterns of sequence variation provided insight into the molecular evolution of the mtDNA control region in closely related felids. Received: 26 January 1998 / Accepted: 14 May 1998     

Phylogenetic Analysis of Vertebrate Lactate Dehydrogenase (LDH) Multigene Families

In this paper we analyzed 49 lactate dehydrogenase (LDH) sequences, mostly from vertebrates. The amino acid sequence differences were found to be larger for a human–killifish pair than a human–lamprey pair. This indicates that some protein sequence convergence may occur and reduce the sequence differences in distantly related species. We also examined transitions and transversions separately for several species pairs and found that the transitions tend to be saturated in the distantly related species pair, while transversions are increasing. We conclude that transversions maintain a conservative rate through the evolutionary time. Kimura's two-parameter model for multiple-hit correction on transversions only was used to derive a distance measure and then construct a neighbor-joining (NJ) tree. Three findings were revealed from the NJ tree: (i) the branching order of the tree is consistent with the common branch pattern of major vertebrates; (ii) Ldh-A and Ldh-B genes were duplicated near the origin of vertebrates; and (iii) Ldh-C and Ldh-A in mammals were produced by an independent gene duplication in early mammalian history. Furthermore, a relative rate test showed that mammalian Ldh-C evolved more rapidly than mammalian Ldh-A. Under a two-rate model, this duplication event was calibrated to be approximately 247 million years ago (mya), dating back to the Triassic period. Other gene duplication events were also discovered in Xenopus, the first duplication occurring approximately 60–70 mya in both Ldh-A and Ldh-B, followed by another recent gene duplication event, approximately 20 mya, in Ldh-B. Received: 5 October 2001 / Accepted: 24 October 2001     

Evolutionary History of B1 Retroposons in the Genus Mus

Short interspersed DNA elements (SINEs) amplify by retroposition either by (i) successive waves of amplification from one or a few evolving master genes or by (ii) the generation of new master genes that coexist with their progenitors. Individual, highly conserved, elements of the B1 SINE family were identified from the GenBank nucleotide database using various B1 subfamily consensus query sequences to determine their integration times into the mouse genome. A comparison of orthologous loci in various species of the genus Mus demonstrated that four subfamilies of B1 elements have been amplifying within the last 1–3 million years. Therefore, B1 sequences are generated by coexisting source genes. Additionally, three B1 subfamilies have been concurrently propagated during subspecies divergence and strain formation in Mus, indicating very recent activity of this retroposon family. The patterns of intra- and interspecies variations of orthologous loci demonstrate the usefulness of B1 integrations as a phylogenetic tool. A single inconsistency in the phylogenetic trends was depicted by the presence of a B1 insert in an orthologous locus exclusively in M. musculus and M. pahari. However, DNA sequence analysis revealed that these were independent integrations at the same genomic site. One highly conserved B1 element that integrated at least 4–6 million years ago suggests the possibility of occasional function for B1 integrations. Received: 25 February 2000 / Accepted: 5 June 2000     

Ascidian and Amphioxus Adh Genes Correlate Functional and Molecular Features of the ADH Family Expansion During Vertebrate Evolution

The alcohol dehydrogenase (ADH) family has evolved into at least eight ADH classes during vertebrate evolution. We have characterized three prevertebrate forms of the parent enzyme of this family, including one from an urochordate (Ciona intestinalis) and two from cephalochordates (Branchiostoma floridae and Branchiostoma lanceolatum). An evolutionary analysis of the family was performed gathering data from protein and gene structures, exon–intron distribution, and functional features through chordate lines. Our data strongly support that the ADH family expansion occurred 500 million years ago, after the cephalochordate/vertebrate split, probably in the gnathostome subphylum line of the vertebrates. Evolutionary rates differ between the ancestral, ADH3 (glutathione-dependent formaldehyde dehydrogenase), and the emerging forms, including the classical alcohol dehydrogenase, ADH1, which has an evolutionary rate 3.6-fold that of the ADH3 form. Phylogenetic analysis and chromosomal mapping of the vertebrate Adh gene cluster suggest that family expansion took place by tandem duplications, probably concurrent with the extensive isoform burst observed before the fish/tetrapode split, rather than through the large-scale genome duplications also postulated in early vertebrate evolution. The absence of multifunctionality in lower chordate ADHs and the structures compared argue in favor of the acquisition of new functions in vertebrate ADH classes. Finally, comparison between B. floridae and B. lanceolatum Adhs provides the first estimate for a cephalochordate speciation, 190 million years ago, probably concomitant with the beginning of the drifting of major land masses from the Pangea. Received: 10 April 2001 / Accepted: 23 May 2001     

Rapid Evolution of the Ribonuclease A Superfamily: Adaptive Expansion of Independent Gene Clusters in Rats and Mice

The two eosinophil ribonucleases, eosinophil-derived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3), are among the most rapidly evolving coding sequences known among primates. The eight mouse genes identified as orthologs of EDN and ECP form a highly divergent, species-limited cluster. We present here the rat ribonuclease cluster, a group of eight distinct ribonuclease A superfamily genes that are more closely related to one another than they are to their murine counterparts. The existence of independent gene clusters suggests that numerous duplications and diversification events have occurred at these loci recently, sometime after the divergence of these two rodent species (∼10–15 million years ago). Nonsynonymous substitutions per site (d _N) calculated for the 64 mouse/rat gene pairs indicate that these ribonucleases are incorporating nonsilent mutations at accelerated rates, and comparisons of nonsynonymous to synonymous substitution (d _N / d _S) suggest that diversity in the mouse ribonuclease cluster is promoted by positive (Darwinian) selection. Although the pressures promoting similar but clearly independent styles of rapid diversification among these primate and rodent genes remain uncertain, our recent findings regarding the function of human EDN suggest a role for these ribonucleases in antiviral host defense. Received: 8 April 1999 / Accepted: 22 June 1999     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

Secondary structure and patterns of evolution among mammalian mitochondrial 12S rRNA molecules

Forty-nine complete 12S ribosomal RNA (rRNA) gene sequences from a diverse assortment of mammals (one monotreme, 11 marsupials, 37 placentals), including 11 new sequences, were employed to establish a ``core' secondary structure model for mammalian 12S rRNA. Base-pairing interactions were assessed according to the criteria of potential base-pairing as well as evidence for base-pairing in the form of compensatory mutations. In cases where compensatory evidence was not available among mammalian sequences, we evaluated evidence among other vertebrate 12S rRNAs. Our results suggest a core model for secondary structure in mammalian 12S rRNAs with deletions as well as additions to the Gutell (1994: Nucleic Acids Res. 22) models for Bos and Homo. In all, we recognize 40 stems, 34 of which are supported by at least some compensatory evidence within Mammalia. We also investigated the occurrence and conservation in mammalian 12S rRNAs of nucleotide positions that are known to participate in the decoding site in E. coli. Twenty-four nucleotide positions known to participate in the decoding site in E. coli also occur among mammalian 12S rRNAs and 17 are invariant for the same base as in E. coli. Patterns of nucleotide substitution were assessed based on our secondary structure model. Transitions in loops become saturated by approximately 10–20 million years. Transitions in stems, in turn, show partial saturation at 20 million years but divergence continues to increase beyond 100 million years. Transversions accumulate linearly beyond 100 million years in both stems and loops although the rate of accumulation of transversions is three- to fourfold higher in loops. Presumably, this difference results from constraints to maintain pairing in stems. Received: 21 June 1995 / Accepted: 25 March 1996     

Sequence of PRAT Satellite DNA ``Frozen' in Some Coleopteran Species

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera (Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50–60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.     

An Estimate of Divergence Time of Parazoa and Eumetazoa and That of Cephalochordata and Vertebrata by Aldolase and Triose Phosphate Isomerase Clocks

Previously we suggested that four proteins including aldolase and triose phosphate isomerase (TPI) evolved with approximately constant rates over long periods covering the whole animal phyla. The constant rates of aldolase and TPI evolution were reexamined based on three different models for estimating evolutionary distances. It was shown that the evolutionary rates remain essentially unchanged in comparisons not only between different classes of vertebrates but also between vertebrates and arthropods and even between animals and plants, irrespective of the models used. Thus these enzymes might be useful molecular clocks for inferring divergence times of animal phyla. To know the divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata, the aldolase cDNAs from Ephydatia fluviatilis, a freshwater sponge, and the TPI cDNAs from Ephydatia fluviatilis and Branchiostoma belcheri, an amphioxus, have been cloned and sequenced. Comparisons of the deduced amino acid sequences of aldolase and TPI from the freshwater sponge with known sequences revealed that the Parazoa–Eumetazoa split occurred about 940 million years ago (Ma) as determined by the average of two proteins and three models. Similarly, the aldolase and TPI clocks suggest that vertebrates and amphioxus last shared a common ancestor around 700 Ma and they possibly diverged shortly after the divergence of deuterostomes and protostomes.     

Molecular phylogeny of the genus Ceanothus (Rhamnaceae) using rbc L and ndh F sequences

 Intrageneric phylogeny among ten representative Ceanothus species was investigated using DNA sequences of the chloroplast encoded ndhF and rbcL genes. Parsimony analysis of the ndhF sequences identified two main clades corresponding to two subgenera Ceanothus and Cerastes. The phylogenetic results suggest that three monophyletic clades within the subgenus Ceanothus can be delimited on the basis of (1) evergreen or (2) deciduous leaves and (3) thorn presence within the evergreen clade. The estimated divergence time based on rbcL sequences suggests that the two subgenera diverged 18–39 million years ago whereas species within each subgenus diverged more recently. Taken together, the results support the division of Ceanothus into two monophyletic subgenera and are consistent with the postulated recent divergence of many species within each subgenus. Received: 25 September 1996/Accepted: 8 November 1996     

Slow Evolutionary Rate of GB Virus C/Hepatitis G Virus

With the aim of elucidating evolutionary features of GB virus C/hepatitis G virus (GBV-C/HGV), molecular evolutionary analyses were conducted using the entire coding region of this virus. In particular, the rate of nucleotide substitution for this virus was estimated to be less than 9.0 × 10⁻⁶ per site per year, which was much slower than those for other RNA viruses. The phylogenetic tree reconstructed for GBV-C/HGV, by using GB virus A (GBV-A) as outgroup, indicated that there were three major clusters (the HG, GB, and Asian types) in GBV-C/HGV, and the divergence between the ancestor of GB- and Asian-type strains and that of HG-type strains first took place more than 7000–10,000 years ago. The slow evolutionary rate for GBV-C/HGV suggested that this virus cannot escape from the immune response of the host by means of producing escape mutants, implying that it may have evolved other systems for persistent infection. Received: 2 June 1998 / Accepted: 8 August 1998     

Extraordinary Conservation of Cysteines Among Homologous Chironomus Silk Proteins sp185 and sp220

Aquatic larvae of the midge, Chironomus tentans, synthesize a 185-kDa silk protein (sp185) with the cysteine-containing motif Cys-X-Cys-X-Cys (where X is any residue) every 20–28 residues. We report here the cloning and full-length sequence of cDNAs encoding homologous silk proteins from Chironomus pallidivittatus (sp185) and Chironomus thummi (sp220). Deduced amino acid sequences reveal proteins of nearly identical mass composed of 72 blocks of 20–28 residues, 61% of which can be described by the motif X_5–8-Cys-X₅-(Trp/Phe/Tyr)-X₄-Cys-X-Cys-X-Cys. Spatial arrangement of these residues is preserved more than surrounding sequences. cDNA clones enabled us to map the genes on polytene chromosomes and identify for the first time the homolog of the Camptochironomus Balbiani ring 3 locus in Chironomus thummi. The apparent molecular weight difference between these proteins (185 vs 220 kDa) is not attributable to primary structure and may be due to differential N-linked glycosylation. DNA distances and codon substitutions indicate that the C. tentans and C. pallidivittatus genes are more related to each other than either is to C. thummi; however, substitution rates for the 5′- and 3′-halves of these genes are different. Blockwise sequence comparisons suggest intragenic variation in that some regions evolved slower or faster than the mean and may have been subjected to different selective pressures. Received: 30 August 1996 / Accepted: 6 November 1996     

Mutation Pattern Variation Among Regions of the Primate Genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (ΨAS-A1, ΨAS-A3, ΨAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene ΨAS-A1, ΨAS-A3, and ΨAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of ΨAS-A1, ΨAS-A3, and ΨAS-3, respectively. These ``regional' variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a ``sequence effect' has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region—i.e., the pseudogene appears as mutation ``hot' spots embedded in ``cold' regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30–40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution. Received: 10 February 1997 / Accepted: 21 April 1997