A rapid loss of stripes: the evolutionary history of the extinct quagga

Twenty years ago, the field of ancient DNA was launched with the publication of two short mitochondrial (mt) DNA sequences from a single quagga (Equus quagga) museum skin, an extinct South African equid (Higuchi et al. 1984 Nature312, 282-284). This was the first extinct species from which genetic information was retrieved. The DNA sequences of the quagga showed that it was more closely related to zebras than to horses. However, quagga evolutionary history is far from clear. We have isolated DNA from eight quaggas and a plains zebra (subspecies or phenotype Equus burchelli burchelli). We show that the quagga displayed little genetic diversity and very recently diverged from the plains zebra, probably during the penultimate glacial maximum. This emphasizes the importance of Pleistocene climate changes for phylogeographic patterns in African as well as Holarctic fauna.     

High variation and very low differentiation in wide ranging plains zebra (Equus quagga): insights from mtDNA and microsatellites

Patterns of genetic differentiation in the plains zebra ( Equus quagga ) were analysed using mitochondrial DNA control region variation and seven microsatellites. The six morphologically defined subspecies of plains zebra lacked the population genetic structure indicative of distinct evolutionary units. Both marker sets showed high levels of genetic variation and very low levels of differentiation. There was no geographical structuring of mitochondrial DNA haplotypes in the phylogenetic tree, and the plains zebra showed the lowest overall differentiation recorded in any African ungulate studied so far. Arid-adapted African ungulates have shown significant regional genetic structuring in support of the Pleistocene refuge theory. This was not the case in the zebra, and the data are discussed in relation to the impact of Pleistocene climate change on a nonbovid member of the savannah ungulate community. The only other species showing a similar absence of genetic structuring is the African buffalo ( Syncerus caffer ), but this taxon lacks the high levels of morphological variation present in the plains zebra.     

A survey of equid mitochondrial DNA: Implications for the evolution, genetic diversity and conservation of Equus

The evolution, taxonomy and conservation of the genus Equuswere investigated by examining the mitochondrial DNA sequences of thecontrol region and 12S rRNA gene. The phylogenetic analysis of thesesequences provides further evidence that the deepest node in thephylogeny of the extant species is a divergence between twolineages; one leading to the ancestor of modern horses (E.ferus, domestic and przewalskii) and the other to thezebra and ass ancestor, with the later speciation events of the zebrasand asses occurring either as one or more rapid radiations, or withextensive secondary contact after speciation. Examination of the geneticdiversity within species suggested that two of the E. hemionussubspecies (E. h. onager and E. h. kulan) onlyrecently diverged, and perhaps, are insufficiently different to beclassified as separate subspecies. The genetic divergence betweendomestic and wild forms of E. ferus (horse) and E.africanus (African ass) was no greater than expected within anequid species. In E. burchelli (plains zebra) there was anindication of mtDNA divergence between populations increasing withdistance. The implications of these results for equid conservation arediscussed and recommendations are made for conservation action.     

Variation in salmonid mitochondrial DNA: Evolutioinary constraints and mechanisms of substitution

Summary Sequence comparisons were made from 2214 bp of mitochondrial DNA cloned from six Pacific salmonid species. These sequences include the genes for ATPase subunit 6, cytochrome oxidase subunit 3, NADH dehydrogenase subunit 3, NADH dehydrogenase subunit 4L, tRNA^GLY, and tRNA^ARG. Variation is found at 338 silent and 12 nonsilent positions of protein coding genes and 10 positions in the two tRNA sequences. A single 3-bp length difference was also detected. In all pairwise comparisons the sequence divergence observed in the fragment was higher than that previously predicted by restriction enzyme analysis of the entire molecule. The inferred evolutionary relationship of these species is consistent between methods. The distribution of silent variation shows a complex pattern with greatly reduced variation at the junctions of genes. The variation in the tRNA sequences is concentrated in the DHU loop. The close relationship of these species and extensive sequence analyzed allows for an analysis of the spectrum of substitutions that includes the frequencies of all 12 possible substitutions. The observed spectrum of substitutions is related to potential pathways of spontaneous substitution. The salmonid sequences show an extremely high ratio of silent to replacement substitutions. In addition the amino acid sequences of the four proteins coded in this fragment show a consistently high level of identity with theXenopus sequences. Taken together these data are consistent with a slower rate of amino acid substitution among the cold-blooded vertebrates when compared to mammals.     

Protamine P1 sequences in equids: comparison with even-toed animals

Protamine P1 amino acid sequences were determined from semen samples of the Przewalski horse, donkey, Somali wild ass, Grevy's zebra, and Grant's zebra (odd-toed perissodactyls), and compared with those of the domestic horse. Although the rate of amino acid variation of protamine P1 is known to be among the most rapidly diverging polypeptides, the equid sequences revealed only little variation. The sequence from the Przewalski horse was identical with that from the domestic horse. The other sequences differed from the corresponding sequences of the domestic and Przewalski horses in two positions-Ser29 was replaced by Cys and Gln32 was replaced by Arg. The presence of the Cys residue at position 29 in the protamine P1 from the zebras, the donkey, and the Somali wild ass may allow formation of one extra protamine disulfide bridge during chromosome condensation in these species. Comparison with protamines from various even-toed animals (artiodactyls) indicated amino acid changes specific for those but different from the equid sequences.     

Development of a molecular diagnostic system to discriminate Dreissena polymorpha (zebra mussel) and Dreissena bugensis (quagga mussel)

A 3-primer PCR system was developed to discriminate invasive zebra (Dreissena polymorpha) and quagga (Dreissena bugensis) mussel. The system is based on: 1) universal primers that amplifies a region of the nuclear 28s rDNA gene from both species and 2) a species-specific primer complementary to either zebra or quagga mussel. The species-specific primers bind to sequences between the binding sites for the universal primers resulting in the amplification of two products from the target species and one product from the nontarget species. Therefore, nontarget products are positive amplification controls. The 3-primer system accurately discriminated zebra and quagga mussels from seven geographically distinct populations.     

Modern Northern Domestic Horses Carry Mitochondrial DNA Similar to Przewalski’s Horse

Several recent studies have suggested past gene flow between the Przewalski’s horse and modern domestic horse and questioned the wild origin of the Przewalski’s horse. Mitochondrial DNA has placed representatives of the Przewalski’s horse into three among the eighteen haplogroups detected from the modern horse. Of these, two haplogroups have so far been found exclusively in the Przewalski’s horse, while the one shared with the domestic horse includes captive individuals that have uncertain pedigrees. We recently found five domestic horse individuals of North European horse breeds to carry a mitochondrial haplogroup that was previously confined only to the Przewalski’s horse. These individuals were sequenced for 6039 bp of mitochondrial DNA and used, together with domestic and Przewalski’s horse sequences presenting all horse haplogroups, to examine the phylogenetic relationships and to date the divergence time between Przewalski’s and domestic horse clusters within this haplogroup. The divergence was dated to have likely occurred about 13,300–11,400 years ago, which coincides with the time of the Younger Dryas.

     

Mitochondrial DNA evolution in the genus Equus

Employing mitochondrial DNA (mtDNA) restriction-endonuclease maps as the basis of comparison, we have investigated the evolutionary affinities of the seven species generally recognized as the genus Equus. Individual species' cleavage maps contained an average of 60 cleavage sites for 16 enzymes, of which 29 were invariant for all species. Based on an average divergence rate of 2%/Myr, the variation between species supports a divergence of extant lineages from a common ancestor approximately 3.9 Myr before the present. Comparisons of cleavage maps between Equus przewalskii (Mongolian wild horse) and E. caballus (domestic horse) yielded estimates of nucleotide sequence divergence ranging from 0.27% to 0.41%. This range was due to intraspecific variation, which was noted only for E. caballus. For pairwise comparisons within this family, estimates of sequence divergence ranged from 0% (E. hemionus onager vs. E. h. kulan) to 7.8% (E. przewalskii vs. E. h. onager). Trees constructed according to the parsimony principle, on the basis of 31 phylogenetically informative restriction sites, indicate that the three extant zebra species represent a monophyletic group with E. grevyi and E. burchelli antiquorum diverging most recently. The phylogenetic relationships of E. africanus and E. hemionus remain enigmatic on the basis of the mtDNA analysis, although a recent divergence is unsupported.      

Palaeogenomics

In many ways, palaeogenomics began when the first ancient DNA sequence was reported. This first sequence was derived from a stuffed museum specimen of the quagga, an extinct mammal related to the zebra. Unspecified and unselected DNA was extracted from the quagga specimen, cloned into a bacterial library, and then sequenced. It took another 17 years and the development of PCR before two independent groups successfully sequenced the complete mitochondrial genomes from several extinct moa species. Only 4 years later, using the original approach of cloning nonspecific ancient DNA extract and shotgun sequencing, the first ancient nuclear DNA sequences were determined, this time from the extinct cave bear. Since these early successes, palaeogenomics has rapidly expanded, because of both technological development and increasing interest in ancient DNA research. New methods, developed since the cave bear sequence was reported, have produced nuclear DNA on a megabase scale from two extinct species, the mammoth and the Neanderthal, our closest relative. For both species, low-coverage genome-sequencing projects have been proposed. It is likely that these will be successful, given the rapid technical development in sequencing techniques. This review carefully examines both the promise and the current limitations of palaeogenomic analyses for both mitochondrial and nuclear DNA.     

Nonneutral Mitochondrial DNA Variation in Humans and Chimpanzees

We sequenced the NADH dehydrogenase subunit 3 (ND3) gene from a sample of 61 humans, five common chimpanzees, and one gorilla to test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution. Within humans and within chimpanzees, the ratio of replacement to silent nucleotide substitutions was higher than observed in comparisons between species, contrary to neutral expectations. To test the generality of this result, we reanalyzed published human RFLP data from the entire mitochondrial genome. Gains of restriction sites relative to a known human mtDNA sequence were used to infer unambiguous nucleotide substitutions. We also compared the complete mtDNA sequences of three humans. Both the RFLP data and the sequence data reveal a higher ratio of replacement to silent nucleotide substitutions within humans than is seen between species. This pattern is observed at most or all human mitochondrial genes and is inconsistent with a strictly neutral model. These data suggest that many mitochondrial protein polymorphisms are slightly deleterious, consistent with studies of human mitochondrial diseases.     

Isolation of Y chromosome-specific microsatellites in the horse and cross-species amplification in the genus Equus

Y chromosome polymorphisms such as microsatellites or single nucleotide polymorphisms represent a paternal counterpart to mitochondrial DNA (mtDNA) for evolutionary and phylogeographic studies. The use of Y chromosome haplotyping in natural populations of species other than humans is still hindered by the lack of sequence information necessary for polymorphism screening. Here we used representational difference analysis (RDA) followed by a screen of a bacterial artificial chromosome (BAC) library for repetitive sequences to obtain polymorphic Y-chromosomal markers. The procedure was performed for the domestic horse (Equus caballus) and we report the first six Y-chromosomal microsatellite markers for this species. Three markers were also useful for haplotyping taxa of the zebra/ass lineage. Y-chromosomal microsatellite markers show a single haplotype in the domestic horse, whereas notable variation has been observed in the other members of the genus Equus.     

Fixed nucleotide differences on the Y chromosome indicate clear divergence between Equus przewalskii and Equus caballus

The phylogenetic relationship between Equus przewalskii and E. caballus is often a matter of debate. Although these taxa have different chromosome numbers, they do not form monophyletic clades in a phylogenetic tree based on mtDNA sequences. Here we report sequence variation from five newly identified Y chromosome regions of the horse. Two fixed nucleotide differences on the Y chromosome clearly display Przewalski's horse and domestic horse as sister taxa. At both positions the Przewalski's horse haplotype shows the ancestral state, in common with the members of the zebra/ass lineage. We discuss the factors that may have led to the differences in mtDNA and Y-chromosomal observations.     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary 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.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

Diagnostic genetic markers and evolutionary relationships among invasive dreissenoid and corbiculoid bivalves in North America: phylogenetic signal from mitochondrial 16S rDNA.

Diagnostic genetic markers from 486 aligned nucleotide sequences of mitochondrial 16S ribosomal DNA were developed for the four closely related species of dreissenoid and corbiculoid bivalves that have invaded North America; the zebra mussel Dreissena polymorpha, the quagga mussel D. bugensis, and the dark false mussel Mytilopsis leucophaeata of the superfamily Dreissenoidea, and the Asian clam Corbicula fluminea of the sister superfamily Corbiculoidea. Evolutionary relationships were examined among the four genera and comparisons were made with native Eurasian populations of D. polymorpha and D. bugensis. Tests were conducted for gender-specific mitochondrial lineages, which occur in some other bivalves. Genetic variability and divergence rates were tested between stem (paired) and loop (unpaired) regions of secondary structure. There were 251 variable nucleotide sites, of which 99 were phylogenetically informative. Overall transition to transversion ratio was 0.76:1.00 and both accumulated linearly in stem and loop regions, suggesting appropriate phylogenetic signal. Genetic distance calibration with the fossil record estimated the pairwise sequence divergence as 0. 0057 +/- 0.0004 per million years. Mytilopsis and Dreissena appear to have diverged about 20.7 +/- 2.7 million years ago. D. bugensis and D. polymorpha appear separated by about 13.2 +/- 2.2 million years. No intraspecific variation was found, including between Eurasian and North American populations, among shallow and deep morphotypes of D. bugensis and between the sexes. Restriction endonuclease markers were developed to distinguish among the species at all life history stages, allowing rapid identification in areas of sympatric distribution.     

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

The complete mitochondrial genome of dhole <Emphasis Type="Italic">Cuon alpinus</Emphasis>: phylogenetic analysis and dating evolutionary divergence within canidae

The dhole (Cuon alpinus) is the only existent species in the genus Cuon (Carnivora: Canidae). In the present study, the complete mitochondrial genome of the dhole was sequenced. The total length is 16672 base pairs which is the shortest in Canidae. Sequence analysis revealed that most mitochondrial genomic functional regions were highly consistent among canid animals except the CSB domain of the control region. The difference in length among the Canidae mitochondrial genome sequences is mainly due to the number of short segments of tandem repeated in the CSB domain. Phylogenetic analysis was progressed based on the concatenated data set of 14 mitochondrial genes of 8 canid animals by using maximum parsimony (MP), maximum likelihood (ML) and Bayesian (BI) inference methods. The genera Vulpes and Nyctereutes formed a sister group and split first within Canidae, followed by that in the Cuon. The divergence in the genus Canis was the latest. The divarication of domestic dogs after that of the Canis lupus laniger is completely supported by all the three topologies. Pairwise sequence divergence data of different mitochondrial genes among canid animals were also determined. Except for the synonymous substitutions in protein-coding genes, the control region exhibits the highest sequence divergences. The synonymous rates are approximately two to six times higher than those of the non-synonymous sites except for a slightly higher rate in the non-synonymous substitution between Cuon alpinus and Vulpes vulpes. 16S rRNA genes have a slightly faster sequence divergence than 12S rRNA and tRNA genes. Based on nucleotide substitutions of tRNA genes and rRNA genes, the times since divergence between dhole and other canid animals, and between domestic dogs and three subspecies of wolves were evaluated. The result indicates that Vulpes and Nyctereutes have a close phylogenetic relationship and the divergence of Nyctereutes is a little earlier. The Tibetan wolf may be an archaic pedigree within wolf subspecies. The genetic distance between wolves and domestic dogs is less than that among different subspecies of wolves. The domestication of dogs was about 1.56–1.92 million years ago or even earlier.     

Mammalian genes as molecular clocks?

In analyzing the silent nucleotide substitutions in some mammalian mitochondrial mRNA coding genes, we had found that the frequency of each of the four nucleotides in rat, mouse, and cow, but not in humans, is the same in the silent third codon position (Lanave C, Preparata G, Saccone C, Serio G (1984) J Mol Evol 20:86-93). Because our findings for these three species were compatible with a stationary Markov process for the evolution of nucleotide sequences, we applied such a model to calculate the effective evolutionary silent substitution rate (vs) and the divergence times among the species. In this paper we have analyzed the first and second codon positions in the same mammalian mitochondrial genes. We found that in the first and second codon positions the human mitochondrial genes satisfy the stationarity conditions. This has allowed us to use the stochastic model mentioned above to calculate the divergence times among mouse, rat, cow, and human. Furthermore, we have analyzed the silent substitution rate in one nuclear gene for these four mammals. We found that in this gene the effective silent substitution rate is about 3 times lower than in mitochondrial genes, and that humans are in this case stationary with respect to the other three mammals in the third codon position as well. Application of our Markov model to this latter gene yields divergence times consistent with our previous determinations.     

Analysis of genetic diversity of domestic water buffaloes and anoas based on variations in the mitochondrial gene for cytochrome b

There are two major groups of domestic water buffaloes in East and Southeast Asia: swamp buffaloes and river buffaloes. Genetic diversity among swamp and river buffaloes was studied by DNA sequence analysis of the mitochondrial gene for cytochrome b. The results showed that each of the two groups has mitochondrial DNA (mtDNA) with a specific cytochrome b haplotype. The pairwise nucleotide sequence divergence was calculated to be 2·67% between swamp and river buffaloes, suggesting that they might have diverged from the ancestral populations of Asiatic domestic water buffaloes, approximately 1 million years ago. In addition, the sequences of the same gene from three subspecies of anoa (lowland, mountain and quarles anoa) were determined and compared with that of a domestic water buffalo. The sequence divergence was 1·2% for mountain anoa vs quarles anoa, 3·6% for mountain anoa vs lowland anoa and 3·3% for quarles anoa vs lowland anoa. Moreover, the sequence divergence between water buffaloes and anoas was found to be approximately 3·33%. Our results provide molecular evidence to support the taxonomic classification, namely, that Asiatic buffaloes may be classified into four lineages, swamp buffalo, river buffalo, lowland anoa and mountain plus quarles anoa. However, the sequence divergence values among these four groups were lower than the sequence divergence values found in the genus and subgenus levels within the subfamily Bovinae. In particular, in contrast to some proposed taxonomic classifications, our results indicated that mtDNA in the water buffaloes and anoas did not diverge at the genus level.