Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
Phylogenetic relationships of the order Carnivora have been extensively studied. However, depending on the type of data, species sampling, and method of analysis used, carnivores have been placed in nearly every possible position throughout the eutherian phylogenetic tree (for a review, see Novacek 1992
). The integrity of the order itself has, however, remained intact. Undisputedly monophyletic, the order Carnivora nonetheless constitutes a very adaptable and heterogeneous group (Wayne et al. 1989
) whose evolution has been marked by several events of parallel or convergent evolution and rapid radiation. Consequently, phylogenetic relationships between and within many carnivore families are also still largely unresolved. For example, the families Procyonidae, Ailuridae, Mustelidae, Pinnidedia, and Ursidae are joined at a polytomy. Within Pinnipedia, the closer affinity of Otariidae to either Phocidae or Obodenidae has long been debated. Additionally, the monophyly of the family Mustelidae has been challenged, with true mustelids and mephitids thought to belong to different branches within the superfamily Musteloidea, also including Procyonidae and Ailuridae. The phylogeny of the family Ursidae is another typical example. Six species of bears, grouped under the subfamily Ursinae, are thought to have diverged at the beginning of the Pliocene through such a rapid radiation event (Thenius 1990
) that the order of species divergence is difficult to determine. Hence, the speciation event that led to extant ursine bears (American black bear, Ursus americanus; brown bear, U. arctos; polar bear, U. maritimus; Asiatic black bear, U. thibetanus; sun bear, U. malayanus; and sloth bear, U. ursinus) is usually represented as a polytomy.
Rapid radiations remain unresolved largely because available data contain very few informative changes essential to infer the correct phylogenetic tree. More resolving power can be reached by using longer DNA sequences to provide a sufficient number of informative characters. Moreover, it has been shown that the combination of multiple genes in a single large data set has the potential to raise a weak phylogenetic signal above the noise level (Bull et al. 1993
). Combining these two advantages, complete mitochondrial (mtDNA) sequences (approximately 17 kbp and 13 protein-coding genes) offer great potential for recovering phylogenies.
With the specific goal of resolving the ursine polytomy, and keeping in mind the variety of questions related to carnivore phylogeny that are yet to be resolved, a new approach to complete mitochondrial genome sequencing was developed to obtain complete mitochondrial sequences from closely related carnivore species. A series of primers were designed, based on conserved regions identified from an alignment of published complete mitochondrial genomes from carnivores (harbor seal, Phoca vitulina, X63768, Arnason and Johnsson 1992
; gray seal, Halichoerus grypus, X72004, Arnason and Gullberg 1993
; domestic cat, Felis catus, U20753, Lopez, Cevario, and O'Brien 1996
; and dog, Canis familiaris, U96639, Kim et al. 1998
). Eleven primer pairs were designed for the amplification of fragments covering the entire mitochondrial genome. Each fragment is longer than 1.4 kbp and covers more than one gene, therefore, limiting the risk of amplifying nuclear copies. The primer design also ensures that there is sufficient overlap of the fragments, in order to obtain the sequence of the primer sites and their flanking nucleotides. Also, 312 internal primers were designed to complete the sequence of each fragment in both the directions. A list of all the primers used to sequence the complete mitochondrial genome from a bear species is presented in table 1
.
|
|
To confirm the authenticity and the position of our new sequences, the cytochrome b gene from our complete mtDNAs was analyzed along with previously published cytochrome b sequences from the eight extant species of bears. To facilitate comparisons with previous studies, the weighted parsimony method used by Talbot and Shields (1996)
as well as the maximum likelihood method using the HKY + gamma model used by Waits et al. (1999)
were applied to a data set comprising cytochrome b sequences from Talbot and Shields (1996)
, an additional Asiatic black bear sequence (Matsuhashi et al. 1999
), and our three new sequences. Two different trees that could, however, not be distinguished statistically were obtained. This is not surprising because previous studies have shown that such short sequences are insufficient to resolve the ursine polytomy. At this point in our study, our purpose is not to thoroughly revisit bear phylogeny but mainly to validate our approach. A few findings are relevant to this objective. The cytochrome b sequence from the black bear used in this study is characteristic of the western haplotype previously identified in several independent studies (Cronin et al. 1991
; Paetkau and Strobeck 1996
; Byun, Koop, and Reimchen 1997
; Wooding and Ward 1997
). With a "western" North American black bear and an Asiatic black bear from Japan (Matsuhashi et al. 1999
) being added to the original study of Talbot and Shields (1996)
, who used an "eastern" black bear and an Asiatic black bear, unfortunately of an unknown origin, the bootstrap support for a clade comprising the two black bear species was improved from 40 to 98. This analysis suggests that the use of highly divergent conspecifics, as we did in this case for the two black bear species, can substantially improve confidence in the tree topology. A few representatives from each ursid species, corresponding to subspecific clades, will likely be helpful in resolving all of the nodes of the ursine polytomy, thereby emphasizing the need for an efficient method to obtain sequence data from very closely related taxa.
The complete mitochondrial sequences obtained here were used in a preliminary phylogenetic analysis, along with previously published sequences from other carnivores (as listed previously). Phylogenetic analyses were performed using various data subsets from the complete genome nucleotide sequence. Artiodactyl and perissodactyl sequences (Bos taurus, V00654, Anderson et al. 1982
; and Equus caballus, X79547, Xu and Arnason 1994
) were used to root the trees. Maximum parsimony, maximum likelihood, and neighbor-joining methods for estimating trees from complete sequence data gave the same tree topology. This corroborates the conclusions of Russo, Takezaki, and Nei (1996)
showing that when many mitochondrial genes are used, all methods appear to converge toward a single tree. The tree topology obtained corresponds to what is currently accepted as the putative true tree (see Cao et al. 1994
; Kuma and Miyata 1994
; Janke, Xu, and Arnason 1997
; Cao et al. 1998
). This analysis verifies that extensive sequence from the mitochondrial genome provides reliable data and sufficient resolution for building mammalian phylogenies (Honeycutt and Adkins 1993
) and distinguishing nodes resulting from rapid radiation episodes such as the ursine speciation events.
Our strategy offers many advantages when the ultimate goal is to sequence numerous individuals of the same or closely related species. A complete mitochondrial sequence can be obtained in just 73 sequencing reactions, that is, at least 10 times less than procedures that involve cloning of random fragments. Moreover, the sequencing procedure is completely PCR-based, thus bypassing purification of mitochondria, DNA restriction digestion, and cloning steps common to other methods.
Finally, the 11 primer pairs used for the amplification of mitochondrial fragments are located in highly conserved regions so as to transfer easily to other carnivores. They have indeed been used to amplify mtDNA from three other ursids, nine phocids, one odobenid, two otariids, five true mustelids, a mephitid, one procyonid, three canids, and two felids. Although we expect the majority of the primers to be widely useful, some of the internal primers presented here will show unavoidable mismatches that may lead to poorer PCR performance. Consequently, only a few new primers may be required to complete the sequence of any carnivore mitochondrial genome. Moreover, as more species are sequenced, a list of alternative primers will become available to facilitate further sequencing. This opens the door to various complete mtDNA phylogenetic analyses at the order, family, subfamily, and species levels among carnivores and should certainly lead to a better understanding of carnivore phylogeny.
Acknowledgements
I.D. holds a scholarship from the NSERC. C.S. is supported by grants from the NSERC and from Parks Canada. We would also like to thank two anonymous reviewers of an early version of this manuscript for their helpful comments.
Footnotes
Keywords: direct sequencing
complete mitochondrial genome
mtDNA
phylogeny
Carnivora
Ursidae
Address for correspondence and reprints: Isabelle Delisle, Department of Biological Sciences, CW 405, Biological Sciences Building, University of Alberta, Edmonton, Alberta, Canada T6G 2E9. idelisle{at}ualberta.ca
.
References
Anderson S., M. H. de Bruijn, A. R. Coulson, I. C. Eperon, F. Sanger, I. G. Young, 1982 Complete sequence of bovine mitochondrial DNA. Conserved features of the mammalian mitochondrial genome J. Mol. Biol 156:683-717[ISI][Medline]
Arnason U., A. Gullberg, 1993 Comparison between the complete mtDNA sequences of the blue and the fin whale, two species that can hybridize in nature J. Mol. Evol 37:312-322[ISI][Medline]
Arnason U., E. Johnsson, 1992 The complete mitochondrial DNA sequence of the harbor seal, Phoca vitulina J. Mol. Evol 34:493-505[ISI][Medline]
Bull J. J., J. P. Huelsenbeck, C. W. Cunningham, D. L. Swofford, P. J. Waddell, 1993 Partitioning and combining data in phylogenetic analysis Syst. Biol 42:384-397[ISI]
Byun S. A., B. F. Koop, T. E. Reimchen, 1997 North American black bear mtDNA phylogeography: implications for morphology and the Haida Gwaii glacial refugium controversy Evolution 51:1647-1653[ISI]
Cao Y., J. Adachi, A. Janke, S. Pääbo, M. Hasewaga, 1994 Phylogenetic relationships among eutherian orders estimated from inferred sequences of mitochondrial proteins: instability of a tree based on a single gene J. Mol. Evol 39:519-527[ISI][Medline]
Cao Y., A. Janke, P. J. Waddell, M. Westerman, O. Takenaka, S. Murata, N. Okada, S. Pääbo, M. Hasegawa, 1998 Conflict among individual mitochondrial proteins in resolving the phylogeny of eutherian orders J. Mol. Evol 47:307-322[ISI][Medline]
Cronin M. A., S. C. Amstrup, G. W. Garner, E. R. Vyse, 1991 Interspecific and intraspecific mitochondrial DNA variation in North American bears (Ursus) Can. J. Zool 69:2985-2992[ISI]
Hoelzel A. R., J. V. Lopez, G. A. Dover, S. J. O'Brien, 1994 Rapid evolution of a heteroplasmic repetitive sequence in the mitochondrial DNA control region of Carnivores J. Mol. Evol 39:191-199[ISI][Medline]
Honeycutt R. L., R. M. Adkins, 1993 Higher level systematics of eutherian mammals: an assessment of molecular characters and phylogenetic hypotheses Annu. Rev. Ecol. Syst 24:279-305[ISI]
Janke A., X. Xu, U. Arnason, 1997 The complete mitochondrial genome of the wallaroo (Macropus robustus) and the phylogenetic relationship among Monotremata, Marsupialia and Eutheria Proc. Natl. Acad. Sci. USA 94:1276-1281
Kim K. S., S. E. Lee, H. W. Jeong, J. H. Ha, 1998 The complete nucleotide sequence of the domestic dog (Canis familiaris) mitochondrial genome Mol. Phylogenet. Evol 10:210-220[ISI][Medline]
Kuma K., T. Miyata, 1994 Mammalian phylogeny inferred from multiple protein data Jpn. J. Genet 69:555-566[Medline]
Lopez J. V., S. Cevario, S. J. O'Brien, 1996 Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (Numt) in the nuclear genome Genomics 33:229-246[ISI][Medline]
Matsuhashi T., R. Masuda, T. Mano, M. C. Yoshida, 1999 Microevolution of the mitochondrial DNA control region in the Japanese brown bear (Ursus arctos) population Mol. Biol. Evol 16:676-684[Abstract]
Novacek M. J., 1992 Mammalian phylogeny: shaking the tree Nature 356:121-125[ISI][Medline]
Paetkau D., C. Strobeck, 1996 Mitochondrial DNA and the phylogeography of Newfoundland black bears Can. J. Zool 74:192-196[ISI]
Russo C. A. M., N. Takezaki, M. Nei, 1996 Efficiencies of different genes and different tree-building methods in recovering a known vertebrate phylogeny Mol. Biol. Evol 13:525-536[Abstract]
Talbot S. L., G. F. Shields, 1996 A phylogeny of the Ursidae inferred from complete sequence of mtDNA genes Mol. Phylogenet. Evol 5:567-575[ISI][Medline]
Thenius E., 1990 Carnivores: phylogeny Pp. 370385 in B. Grzimek, ed. Encyclopedia of mammals, Vol. 3. McGraw-Hill, New York
Waits L., J. Sullivan, S. J. O'Brien, R. H. Ward, 1999 Rapid radiation event in the family Ursidae indicated by likelihood phylogenetic estimation from multiple fragments of mtDNA Mol. Phylogenet. Evol 13:82-92[ISI][Medline]
Wayne R. K., R. E. Benveniste, D. N. Janczewski, S. J. O'Brien, 1989 Molecular and biochemical evolution of the Carnivora Pp. 465488 in J. L. Gittleman, ed. Carnivore behavior, ecology and evolution. Cornell University Press, Ithaca, N.Y
Wooding S., R. Ward, 1997 Phylogeography and Pleistocene evolution in the North American black bear Mol. Biol. Evol 14:1096-1105[Abstract]
Xu X., U. Arnason, 1994 The complete mitochondrial DNA sequence of the horse, Equus caballus: extensive heteroplasmy of the control region Gene 148:357-362[ISI][Medline]