Cytochromeb gene of marine mammals: Phylogeny and evolution

The DNA sequences of the mitochondrial cytochromeb gene of marine mammals (Cetacea, Pinnipedia, Sirenia) were compared with cytochromeb genes of terrestrial mammals including the semiaquatic hippopotamus. The comparison included 28 sequences, representing 22 families and 10 orders. The dugong (order Sirenia) sequence associated with that of the elephant, supporting the Tethytheria clade. The fin whale and dolphin (order Cetacea) sequences are more closely related to those of the artiodactyls, and the comparison suggests that the hippopotamus may be the extant artiodactyl species that is most closely related to the cetaceans. The seal sequence may be more closely related to those of artiodactyls, cetaceans, and perissodactyls than to tethytheres, rodents, lagomorphs, or primates. The cytochromeb proteins of mammals do not evolve at a uniform rate. Human and elephant cytochromeb amino acid sequences were found to evolve the most rapidly, while those of myomorph rodents evolved slowest. The cytochromeb of marine mammals evolves at an intermediate rate. The pattern of amino acid substitutions in marine mammals is similar to that of terrestrial mammals.     

Nucleotide sequence of the D-loop region of the sperm whale (Physeter macrocephalus) mitochondrial genome

We have amplified, by the polymerase chain reaction, and have sequenced the D-loop region of the mitochondrial DNA from the sperm whale (Physeter macrocephalus). The sperm whale D-loop was aligned with D- loop sequences from four other cetaceans (Commerson's dolphin, orca, fin whale, and minke whale) and an out-group (cow). This alignment showed the sperm whale sequence to be larger than that of other cetaceans. In addition, some sequence blocks were highly conserved among all six species, suggesting roles in the functioning of mitochondrial DNA. Other blocks that were previously reported to be well conserved among cetaceans showed little sequence conservation with the sperm whale D-loop, which argues against the functional importance of these sequence blocks in cetaceans.      

Evolution of the cytochromeb gene of mammals

Summary With the polymerase chain reaction (PCR) and versatile primers that amplify the whole cytochromeb gene (∼ 1140 bp), we obtained 17 complete gene sequences representing three orders of hoofed mammals (ungulates) and dolphins (cetaceans). The fossil record of some ungulate lineages allowed estimation of the evolutionary rates for various components of the cytochromeb DNA and amino acid sequences. The relative rates of substitution at first, second, and third positions within codons are in the ratio 10 to 1 to at least 33. For deep divergences (>5 million years) it appears that both replacements and silent transversions in this mitochondrial gene can be used for phylogenetic inference. Phylogenetic findings include the association of (1) cetaceans, artiodactyls, and perissodactyls to the exclusion of elephants and humans, (2) pronghorn and fallow deer to the exclusion of bovids (i. e., cow, sheep, and goat), (3) sheep and goat to the exclusion of other pecorans (i. e., cow, giraffe, deer, and pronghorn), and (4) advanced ruminants to the exclusion of the chevrotain and other artiodactyls. Comparisons of these cytochromeb sequences support current structure-function models for this membrane-spanning protein. That part of the outer surface which includes the Q_o redox center is more constrained than the remainder of the molecule, namely, the transmembrane segments and the surface that protrudes into the mitochondrial matrix. Many of the amino acid replacements within the transmembrane segments are exchanges between hydrophobic residues (especially leucine, isoleucine, and valine). Replacement changes at first and second positions of codons approximate a negative binomial distribution, similar to other protein-coding sequences. At four-fold degenerate positions of codons, the nucleotide substitutions approximate a Poisson distribution, implying that the underlying mutational spectrum is random with respect to position.     

The complete mitochondrial DNA sequence of the harbor seal,Phoca vitulina

Summary The nucleotide sequence of the mitochondrial DNA (mtDNA) of the harbor seal, Phoca vitulina, was determined. The total length of the molecule was 16,826 bp. The organization of the coding regions of the molecule conforms with that of other mammals, but the control region is unusually long. A considerable portion of the control region is made up of short repeats with the motif GTACAC particularly frequent. The two rRNA genes and the 13 peptide-coding genes of the harbor seal, fin whale, cow, human, mouse, and rat were compared and the relationships between the different species assessed. At ordinal level the 12S rRNA gene and 7 out of the 13 peptide-coding genes yielded a congruent topological tree of the mtDNA relationship between the seal, cow, whale, human, and the rodents. In this tree the whale and the cow join first, and this clade is most closely related to the seal.Offprint requests to: Ú. Árnason     

Comparison between the complete mtDNA sequences of the blue and the fin whale,two species that can hybridize in nature

The sequence of the mitochondrial DNA (mtDNA) molecule of the blue whale (Balaenoptera musculus) was determined. The molecule is 16,402 by long and its organization conforms with that of other eutherian mammals. The molecule was compared with the mtDNA of the congeneric fin whale (B. physalus). It was recently documented that the two species can hybridize and that male offspring are infertile whereas female offspring may be fertile. The present comparison made it possible to determine the degree of mtDNA difference that occurs between two species that are not completely separated by hybridization incompatibility. The difference between the complete mtDNA sequences was 7.4%. Lengths of peptide coding genes were the same in both species. Except for a small portion of the control region, disruption in alignment was usually limited to insertion/deletion of a single nucleotide. Nucleotide differences between peptide coding genes ranged from 7.1 to 10.5%, and difference at the inferred amino acid level was 0.0–7.9%. In the rRNA genes the mean transition difference was 3.8%. This figure is similar in degree to the difference (3.4%) between the 12S rRNA gene of humans and the chimpanzee. The mtDNA differences between the two whale species, involving both peptide coding and rRNA genes, suggest an evolutionary separation of 5 million years. Although hybridization between more distantly related mammalian species may not be excluded, it is probable that the blue and fin whales are nearly as different in their mtDNA sequences as hybridizing mammal species may be. Correspondence to: Ú. Árnason     

Evolutionary genetics of the suiformes as reconstructed using mtDNA sequencing

We have amplified and sequnced the entire mitochondrial DNA cytochromeb gene from four species of Suidae: babirusa, warthog, bearded pig, and some specimens belonging to different subspecies and populations of wild and domestic pigs (Sus scrofa). These sequences were aligned with additional mammalian sequences retrieved from the literature and were used to obtain phylogenetic trees of the Suiformes (Artiodactyla). Several species of Carnivora, Perissodactyla. Cetacea, and other Artiodactyla were used as outgroups. Molecular phylogenetic relationships among the Suiformes reflect their current taxonomy: Hippopotamidae, Tayassuidae, and Suidae are separated by deep genetic gaps, and the division of the Suidae into the subfamilies Babyrousinae., Phacochoerinae, and Suinae has strong genetic correlates. Cytochromeb sequences show differences among Asian and Western populations ofSus scrofa, agreeing with other genetic information (karyotypes blood groups, and protein variability). The two Italian subspecies of wild boar have unique mtDNA cytochromeb haplotypes. The evolutionary rates of cytochromeb sequences are different at transitions versus transversions as well as at first, second, and third positions of codons. Therefore, these classes of substitutions reached different levels of mutational saturation. Only transversions and the conservative first and second position substitutions are linearly related to genetic distances among the Suiformes. Therefore, divergence times were computed using unsaturated conserved nucleotide substitutions and calibrated using paleontological divergence times between some Artiodactyla. Transversions apparently evolve at remarkably regular rates in ungulate taxa which have accumulated less than 20% estimated sequence divergence, corresponding to about 40–45 million years of independent evolution. Molecular, information suggests that Hippopotamidae and Tayassuidae are not closely related (as stated by Pickford, 1986, 1989, 1993) and that the origin of babirusa and warthog (about 10–19 and 5–15 million years ago, respectively) is more recent than supported by current evolutionary reconstructions. The inferred origin of bearded pig is about 2.1 million years old, and genetic divergence among differentSus scrofa populations is probably a Pleistocene event. The addition of new sequences of Suiformes does not help in resolving the phylogenetic position ofHippopotamus amphibius, which shows weak but recurrent linkages with the cetacean evolutionary lineage.To whom correspondence should be addressed.     

The Mitochondrial Genome of the Sperm Whale and a New Molecular Reference for Estimating Eutherian Divergence Dates

Extant cetaceans are systematically divided into two suborders: Mysticeti (baleen whales) and Odontoceti (toothed whales). In this study, we have sequenced the complete mitochondrial (mt) genome of an odontocete, the sperm whale (Physeter macrocephalus), and included it in phylogenetic analyses together with the previously sequenced complete mtDNAs of two mysticetes (the fin and blue whales) and a number of other mammals, including five artiodactyls (the hippopotamus, cow, sheep, alpaca, and pig). The most strongly supported cetartiodactyl relationship was: outgroup,((pig, alpaca),((cow, sheep),(hippopotamus,(sperm whale,(baleen whales))))). As in previous analyses of complete mtDNAs, the sister-group relationship between the hippopotamus and the whales received strong support, making both Artiodactyla and Suiformes (pigs, peccaries, and hippopotamuses) paraphyletic. In addition, the analyses identified a sister-group relationship between Suina (the pig) and Tylopoda (the alpaca), although this relationship was not strongly supported. The paleontological records of both mysticetes and odontocetes extend into the Oligocene, suggesting that the mysticete and odontocete lineages diverged 32–34 million years before present (MYBP). Use of this divergence date and the complete mtDNAs of the sperm whale and the two baleen whales allowed the establishment of a new molecular reference, O/M-33, for dating other eutherian divergences. There was a general consistency between O/M-33 and the two previously established eutherian references, A/C-60 and E/R-50. Cetacean (whale) origin, i.e., the divergence between the hippopotamus and the cetaceans, was dated to ≈55 MYBP, while basal artiodactyl divergences were dated to ≥65 MYBP. Molecular estimates of Tertiary eutherian divergences were consistent with the fossil record. Received: 12 July 1999 / Accepted: 28 February 2000     

Sequence analysis of the complete mitochondrial DNA molecule of the hedgehog,Erinaceus europaeus,and the phylogenetic position of the Lipotyphla

The sequence of the mitochondrial DNA (mtDNA) molecule of the European hedgehog (Erinaceus europaeus) was determined. The length of the sequence presented is 17,442 nucleotides (nt). The molecule is thus the largest eutherian mtDNA molecule so far reported. The organization of the molecule conforms with that of other eutherians, but the control region of the molecule is exceptionally long, 1,988 nt, due to the presence of repeated motifs at two different positions in the 3 part of the control region. The length of the control region is not absolute due to pronounced heteroplasmy caused by variable numbers of the motif TACGCA in one of the repetitive regions. The sequence presented includes 46 repeats of this type. The other repeated region is composed of different AT-rich repeats. This region was identical among four clones studied. Comparison of mitochondrial peptide-coding genes identified a separate position of the hedgehog among several mammalian orders. The concatenated protein sequence of the 13 peptide-coding genes was used in a phylogenetic study using the opossum as outgroup. The position of the hedgehog sequence was basal among the other eutherian sequences included: human, rat, mouse, cow, blue whale, harbor seal, and horse. The analysis did not resolve the relationship among carnivores, perissodactyls, and artiodactyls/cetaceans, suggesting a closer relationship among these orders than acknowledged by classical approaches. Correspondence to: U. Arnason     

Multiple nuclear pseudogenes of mitochondrial DNA exist in the canine genome

Many copies of nuclear counterparts of mitochondrial DNA (mtDNA) were found in nuclear DNA from sperm heads of the domestic dog, Canis familiaris, by DNA-DNA hybridization and DNA sequencing. Nuclear counterparts homologous to the mtDNA D-loop region were cloned into lambda phage vectors (EMBL4 and lambda gt11), and nucleotide sequences of seven different mtDNA pseudogenes were then determined. The seven pseudogenes were E3 (474 bp; 82% homology with canine mtDNA), E13 (1867 bp; 67%), 8B (2375 bp; 78%), 12A (2650 bp; 79%), 33 (4131 bp; 86%), 47 (4251 bp; 86%), and E17 (5721 bp; 71%). These seven mtDNA pseudogenes corresponded to portions of cytoplasmic mtDNA containing the genes ile, ND1, leu, 16S rRNA, val, 12S rRNA, phe, D-loop, pro, thr, cytb, and glu. A neighbor-joining phylogenetic tree constructed from 12S rRNA sequences in mtDNA pseudogenes 8B, 33, 47, and E17 and in 10 mtDNA fragments from other species showed that these four pseudogenes form a monophyletic clade with canine mtDNA. A neighbor-joining phylogenetic tree based on the 318-bp cytb region showed that the canine pseudogenes existed before the divergence of 17 related canids, and their divergence dates were calculated at around 4.4 to 8.6 million years ago.     

comparison between the complete mitochondrial DNA sequences ofHomo and the common chimpanzee based on nonchimeric sequences

The complete mitochondrial DNA (mtDNA) molecules ofHomo and of the common chimpanzee were sequenced. Each sequence was established from tissue of one individual and thus nonchimeric. Both sequences were assembled in their entirety from natural (not PCR amplified) clones. Comparison with sequences in the literature identified the chimpanzee specimen asPan troglodytes verus, the West African variety of the species. The nucleotide difference between the complete human and chimpanzee sequences is 8.9%. The difference between the control regions of the two sequences is 13.9% and that between the remaining portions of the sequences 8.5%. The mean amino acid difference between the inferred products of the 13 peptide-coding genes is 4.4%. Sequences of the complete control regions, the complete 12S rRNA genes, the complete cytochromeb genes, and portions of the NADH4 and NADH5 genes of two other chimpanzee specimens showed that they were similar but strikingly different from the same regions of the completely sequenced molecule fromPan troglodytes verus. The two specimens were identified asPan troglodytes troglodytes, the Central African variety of the common chimpanzee.     

A molecular perspective on the phylogeny of placental mammals based on mitochondrial 12S rDNA sequences,with special reference to the problem of the Paenungulata

Part of the 12S rDNA gene was amplified and sequenced for 11 placental mammals, 3 marsupials, and 2 monotremes. Multiple alignments for these sequences and nine additional placental sequences taken from GenBank were obtained using CLUSTAL. Phylogenetic analyses were performed using standard parsimony, transversion parsimony, and Lake's method of invariants. All of our analyses uniteLoxodontia withDugong. Procavia, in turn, is a sister group to these taxa, thus supporting the monophyly of the Paenungulata. Perissodactyls are a sister group to paenungulates when transitions and transversions are both included but not when transitions are omitted. Likewise, cetaceans are a sister group to artiodactyls on minimum length trees under standard parsimony but not under transversion parsimony. Rodent monophyly and bat monophyly also receive mixed support, as does a putative alliance between primates and lagomorphs. Interestingly, the percentage divergence between the echidna and the platypus is less than for the rat and mouse.     

SPATIOTEMPORAL PREDICTION MODELS OF CETACEAN HABITATS IN THE MID-WESTERN NORTH ATLANTIC OCEAN (FROM CAPE HATTERAS, NORTH CAROLINA, U.S.A. TO NOVA SCOTIA, CANADA)

Habitat prediction models were developed for 13 cetacean species of the mid-western North Atlantic Ocean: beaked whale, fin whale, humpback whale, minke whale, pilot whale, sperm whale, bottlenose dolphin, common dolphin, Risso's dolphin, spotted dolphin, whitesided dolphin, and harbor porpoise. Using the multiple logistic regression, sightings of cetaceans during the 1990–1996 summer (June-September) surveys were modeled with oceanographic (sea surface temperature, monthly probability of front occurrence) and topographic (depth, slope) variables for the same period. Predicted habitat maps for June and August were created for each species using a Geographical Information System. The predicted habitat locations matched with current and historic cetacean sighting locations. The model also predicted habitat shifts for some species associated with oceanographic changes. The correct classification rate of the prediction models with 1997–1998 summer survey data ranged from 44% to 70%, of which most of the misclassifications were caused by false positives ( i.e. , absence of sightings at locations where the models predicted).     

More DNA support for a Cetacea/Hippopotamidae clade: the blood-clotting protein gene gamma-fibrinogen

Recent phylogenetic analyses of DNA sequences suggest that cetaceans (whales) and hippopotamid artiodactyls (hippos) are extant sister taxa. Consequently, the shared aquatic specializations of these taxa may be synapomorphies. This molecular view is contradicted by paleontological data that overwhelmingly support a monophyletic Artiodactyla (even-toed ungulates) and a close relationship between Cetacea and extinct mesonychian ungulates. According to the fossil evidence, molecular, behavioral, and anatomical resemblances between hippos and whales are interpreted as convergences or primitive retentions. In this report, competing interpretations of whale origins are tested through phylogenetic analyses of the blood-clotting protein gene gamma- fibrinogen from cetaceans, artiodactyls, perissodactyls (odd-toed ungulates), and carnivores (cats, dogs, and kin). In combination with published DNA sequences, the gamma-fibrinogen data unambiguously support a hippo/whale clade and are inconsistent with the paleontological perspective. If the phylogeny favored by fossil evidence is accepted, the convergence at the DNA level between Cetacea and Hippopotamidae is remarkable in its distribution across three genetic loci: gamma-fibrinogen, the linked milk casein genes, and mitochondrial cytochrome b.      

Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.     

Species identification and likely catch time period of whale bones from South Georgia

Skeletal remains of baleen whales killed during the onset of 20th century commercial whaling lie scattered across the shores and abandoned whaling stations of the subantarctic island of South Georgia. Here we report on genetic species identification of whale bones collected from South Georgia using standard historical DNA protocols. We amplified and sequenced short fragments of the mitochondrial DNA (mtDNA) control region from 281 available bone samples. Of these, 231 provided mtDNA sequences of sufficient quality and length (174–194 bp) for species identification: 158 bones were identified as humpback whale (Megaptera novaeangliae), 51 bones were identified as fin whale (Balaenoptera physalus), 18 bones were identified as blue whale (B. musculus), two bones were identified as sei whale (B. borealis), one bone was identified as a southern right whale (Eubalaena australis), and one bone was identified as a southern elephant seal (Mirounga leonina). The prominence of humpback, fin, and blue whale bones in the sample collection corresponds to the catch record of the early years of whaling on the island of South Georgia (pre‐1915), prior to the depletion of these populations.     

Speciation of two sympatric coastal fish species, Girella punctata and Girella leonina (Perciformes, Kyphosidae)

Girella punctata and Girella leonina are sympatric sister species showing extensive distributional overlap in shallow rocky reefs in the Pacific Ocean south of the Japanese Islands. Differences between the two species in external morphological characters, such as number of pored lateral line scales, colour of opercular flap and shape of caudal fin, are congruent with genetic divergence. Nucleotide identity between the two species in the 3.3 kbp region of partial mitochondrial DNA containing the D-loop region, in 12S and 16S ribosomal RNA (rRNA) and transfer RNA genes is 95%. To estimate divergence time, Bayesian analysis was conducted using a dataset comprising concatenated nucleotide sequences from the two rRNA genes of three girellid and nine other fish species. Using the Elopomorpha – Clupeocephala split (265 million years ago (mya)) as a calibration point, divergence between G. punctata and G. leonina is estimated as having occurred 6.0±1.4 mya. Speciation is suggested to have been caused by geographical isolation associated with formation of the Japanese Islands, which resulted in disjunction of Girella habitat.     

Molecular characterization of a cloned dolphin mitochondrial genome

Summary DNA clones have been isolated that span the complete mitochondrial (mt) genome of the dolphin,Cephalorhynchus commersonii. Hybridization experiments with purified primate mtDNA probes have established that there is close resemblance in the general organization of the dolphin mt genome and the terrestrial mammalian mt genomes. Sequences covering 2381 bp of the dolphin mt genome from the major noncoding region, three tRNA genes, and parts of the genes encoding cytochrome b, NADH dehydrogenase subunit 3 (ND3), and 16S rRNA have been compared with corresponding regions from other mammalian genomes. There is a general tendency throughout the sequenced regions for greater similarity between dolphin and bovine mt genomes than between dolphin and rodent or human mt genomes.     

The main regulatory region of mammalian mitochondrial DNA: Structure-function model and evolutionary pattern

Summary The evolution of the main regulatory region (D-loop) of the mammalian mitochondrial genome was analyzed by comparing the sequences of eight mammalian species: human, common chimpanzee, pygmy chimpanzee, dolphin, cow, rat, mouse, and rabbit. The best alignment of the sequences was obtained by optimization of the sequence similarities common to all these species.The two peripheral left and right D-loop domains, which contain the main regulatory elements so far discovered, evolved rapidly in a species-specific manner generating heterogeneity in both length and base composition. They are prone to the insertion and deletion of elements and to the generation of short repeats by replication slippage. However, the preservation of some sequence blocks and similar cloverleaf-like structures in these regions, indicates a basic similarity in the regulatory mechanisms of the mitochondrial genome in all mammalian species.We found, particularly in the right domain, significant similarities to the telomeric sequences of the mitochondrial (mt) and nuclear DNA ofTetrahymena thermophila. These sequences may be interpreted as relics of telomeres present in ancestral linear forms of mtDNA or may simply represent efficient templates of RNA primase-like enzymes.Due to their peculiar evolution, the two peripheral domains cannot be used to estimate in a quantitative way the genetic distances between mammalian species. On the other hand the central domain, highly conserved during evolution, behaves as a good molecular clock.Reliable estimates of the times of divergence between closely and distantly related species were obtained from the central domain using a Markov model and assuming nonhomogeneous evolution of nucleotide sites.     

Non-invasive genotyping of primates

Using DNA amplified from shed or plucked hair follicles it is now possible to genotype individual primates at many nuclear and mitochondrial gene loci. Sequence specific primers and the polymerase chain reaction permit the rapid production of sufficient DNA from a single hair for numerous analyses. The direct sequencing of relatively conservative mtDNA sequences like cytochromeb is proving useful in establishing species and subspecies-level relationships. More variable sequences (e.g. the mtDNA control region or D-loop) are useful at the population and social community levels. Paternity exclusion, pedigree relationships, and community structure can be determined using simple sequence length polymorphisms (SSLPs) of multiple hypervariable nuclear microsatellite or simple sequence repeat (SSR) loci. Studies involving captive and free-ranging chimpanzees, gibbons, and macaques illustrate the resolving power of these new non-invasive molecular genetic genotyping techniques.     

Structure and evolution of opossum,guinea pig,and porcupine cytochrome b genes

Summary We have sequenced the mitochondrial cytochrome b gene from the guinea pig, the African porcupine, and a South American opossum. A phylogenetic analysis, which includes 22 eutherian and four other vertebrate cytochrome b sequences, indicates that the guinea pig and the porcupine constitute a natural clade (Hystricomorpha) that is not a sister group to the clade of mice and rats (Myomorpha). Therefore, the hypothesis that the Rodentia is paraphyletic receives additional support. The artiodactyls, the perissodactyls, and the cetaceans form a group that is separated from the primates and the rodents. The 26 sequences are used to study the structure/function relationships in cytochrome b, whose function is electron transport. Most of the amino acid residues involved in the two reaction centers are well conserved in evolution. The four histidines that are believed to ligate the two hemes are invariant among the 26 sequences, but their nearby residues are not well conserved in evolution. The eight transmembrane domains represent some of the most divergent regions in the cytochrome b sequence. The rate of nonsynonymous substitution is considerably faster in the human and elephant lineages than in other eutherian lineages; the faster rate might be due to coevolution between cytochrome b and cytochrome c. Offprint requests to: W.-H. Li