*Laboratoire de Paléontologie, Paléobiologie et Phylogénie, Institut des Sciences de l'Evolution, Université Montpellier II, Montpellier, France;
Biology and Biochemistry, Queen's University of Belfast;
Department of Biology, University of California;
Department of Biochemistry, University of Nijmegen, Nijmegen, The Netherlands;
||Bioinformatics, GlaxoSmithKline, Pennsylvania;
¶Institute for Biodiversity and Ecosystems Dynamics, University of Amsterdam
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Solving the phylogenetic position of the order Xenarthra within Mammalia is of primary importance to understand the morphological and biogeographical processes that shaped the early stages of placental evolution. Indeed, despite their highly specialized morphology, xenarthrans retain anatomical and physiological characters thought to be primitive for placental mammals (McKenna 1975
). This composite morphology, mixing ancestral and derived characters, has made the order's position within placentals very difficult to assess (Engelmann 1985
; Gaudin et al. 1996
). On the basis of the retention of numerous archaic features, morphologists (Gregory 1910
; McKenna 1975
; Novacek 1992
; Shoshani and McKenna 1998
) have long proposed that Xenarthra represents the sister-group to all other eutherians thereby named Epitheria. However, morphological synapomorphies defining epitherians are weak, and their phylogenetic distribution is equivocal (Gaudin et al. 1996
). Early studies of complete mitochondrial genomes did not support a basal position for Xenarthra (Arnason, Gullberg, and Janke 1997
). More recent analysesincluding a larger taxon samplingsuggested a sister-group relationship between the nine-banded armadillo (Dasypus novemcinctus) and representatives of the African clade (Waddell et al. 1999
; Cao et al. 2000
; Mouchaty et al. 2000a
). However, complete mitochondrial genome analyses appear to be affected by insufficient taxon samplingespecially within Xenarthralikely responsible for long-branch attraction and rooting artifacts (Waddell et al. 1999
) and suffer from saturation in the deepest parts of the tree (Cao et al. 2000
; Springer et al. 2001
).
By contrast two recent independent analyses based on the large concatenation of mainly nuclear genes for a broad taxon sampling of eutherian mammals have shown that the order Xenarthra on its own represents one of the major clades of placentals (Madsen et al. 2001
; Murphy et al. 2001a,
2001b
). These studies provided convincing evidence for an arrangement of placental orders into four major clades: (I) Afrotheria, (II) Xenarthra, (III) Euarchontoglires and (IV) Laurasiatheria, emphasizing the crucial influence of tectonic events in their early differentiation. A close relationship between the Northern Hemisphere clades III and IV has been proposed (i.e., "Boreoeutheria"; Springer and de Jong 2001
) with either Afrotheria or Xenarthra as the most basal clade, suggesting a Southern Hemisphere origin for eutherian mammals (Eizirik, Murphy, and O'Brien 2001
; Madsen et al. 2001
). In fact, the relationships between these four clades directly depend on the unstable position of the root. Recent application of the Bayesian approach (Yang and Rannala 1997
; Huelsenbeck et al. 2001
) to this problem using a large character sampling supported a basal position of Afrotheria (Murphy et al. 2001b
). However, in this study, Xenarthra suffers from poor taxonomic representation (only three taxa) relative to the other three major clades that include 8 (Afrotheria), 11 (Euarchontoglires), and 20 (Laurasiatheria) taxa.
Here we present a study including the broadest taxonomic representation so far considered in a molecular approach to xenarthran phylogeny. We constructed a supermatrix of 47 placental taxa and three marsupial outgroups, including 12 of the 13 living xenarthran genera, for three genetically independent protein-coding nuclear genes: 2B Adrenergic receptor (ADRA2B), Breast Cancer Susceptibility exon 11 (BRCA1), and von Willebrand Factor exon 28 (VWF) representing a total of 5,130 aligned nucleotide sites. Choice of these nuclear markers was guided by their wide use for inferring the phylogeny of placental mammals (Springer et al. 1997
; Stanhope et al. 1998a,
1998b
; Madsen et al. 2001
), demonstrating that their resolving power was higher than for mitochondrial markers at the mammalian interordinal level (Springer et al. 1999
; Springer et al. 2001
). Furthermore, their coding status provides the opportunity to compare the phylogenetic signal contained in both nucleotides and amino acids. The extensive sampling realized within Xenarthra permits us to resolve intraordinal relationships, and to investigate the impact of an increased taxon sampling on the root position of the placental tree. Increased nucleotide sampling has a dominant impact on phylogenetic accuracy (Poe and Swofford 1999
; Rosenberg and Kumar 2001
), but increased taxon sampling has also been shown to facilitate phylogenetic inference (Lecointre et al. 1993
; Hillis 1996
; Rannala et al. 1998
). Indeed, breaking potentially long branches by adding taxa within one of the two most basal clades might help to stabilize the placental ingroup topology in its deepest parts. By using the maximum likelihood (ML) and the Bayesian framework this study aims to (1) compare the evolutionary dynamics and phylogenetic content of these three protein-coding nuclear markers evolving under different selective pressures, (2) resolve the phylogeny of living xenarthrans with special reference to armadillos, and (3) evaluate the effect of an increased taxon sampling within Xenarthra for finding the root of the placental tree.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Data Acquisition
Xenarthran samples preserved in 95% ethanol were stored in the mammalian tissue collection of the Institut des Sciences de l'Evolution de Montpellier (Catzeflis 1991
). Total DNAs were extracted for D. novemcinctus (nine-banded armadillo), Dasypus kappleri (great long-nosed armadillo), Chaetophractus villosus (larger hairy armadillo), Euphractus sexcinctus (six-banded armadillo), Zaedyus pichiy (pichi), Tolypeutes matacus (southern three-banded armadillo), Cabassous unicinctus (southern naked-tailed armadillo), P. maximus (giant armadillo), Cyclopes didactylus (pygmy anteater), Tamandua tetradactyla (collared anteater), Myrmecophaga tridactyla (giant anteater), Bradypus tridactylus (pale-throated three-toed sloth), and Choloepus didactylus (southern two-toed sloth).
The single exon gene ADRA2B was amplified and sequenced using the primers designed by Springer et al. (1997)
and the following additional pairs of forward (A) and reverse (B) primers: A4 (5'-GCCATCGCGGCNGYCRYCACCTTCCTCATC-3'), B4 (5'-GCTGCGYTTGGCAATCAGGTAGAGTCG-3'), and B5 (5'-GCGCCCAGGCTGTAGCTGAAGAAGAA-3'). The exon 28 of VWF was amplified according to Delsuc et al. (2001)
. PCR products for ADRA2B and VWF were purified from 1% agarose gels using Amicon Ultrafree-DA columns (Millipore) and sequenced on both strands using automatic sequencing (Big Dye Terminator cycle kit) on an ABI 310 (PE Applied Biosystems). BRCA1 gene sequences were obtained as described elsewhere (Teeling et al. 2000
; Madsen et al. 2001
). The 2,900-bp region of exon 11 was amplified in three overlapping segments, as described in Teeling et al. (2000)
and Madsen et al. (2001)
. PCR primers for each segment were as follows: [1] Forward UF1 (5'-GTTTCAAACTTGCATGTGGAGCC-3'), Reverse XR11 (5'-GCAGATTCTTTTTCCAATGATTCTG-3'); [2] Forward GF8 (5'-GGCCTTCATCCTGAGGATTTTATCAA-3'), Reverse R19 (5'-TGYAAATACTGAGTATCAAGTTCACT-3'); [3] Forward XBF17 (5'-TATGGCACTCARGAYAGTATCTCATT-3'), Reverse BRCA1B (5'-GTTGGAAGCAGGGAAGCTCTTCATC-3'). Sequencing was performed using PCR primers and additional internal primers in segment [1] Forward F4 (5'-GAAAGTTAATGAGTGGTTTTCCAGAA-3'), Reverse UR7 (5'-CTTCCTCCGATAGGTTTTCCCAA-3') and in segment [3] Forward F25 (5'-AACTAGGTAGAAACAGAGGRCCTA-3'), Reverse R26 (5'-TTAGGYCCTCTGTTTCTACCTAGTT-3').
The 26 xenarthran sequences new to this study have been deposited in the EMBL data bank. Taxonomy and accession numbers referring to all sequences used in this study are provided as Supplementary Material.
Sequence Alignment
Sequences were manually aligned with the ED editor of the MUST package (Philippe 1993
). We excluded a glutamic acid repeat region of ADRA2B and also a 21-bp region of BRCA1 that is repeated up to four times according to Madsen et al. (2001)
. The concatenated placental data set is 5,130 bp long with ADRA2B (1,152 bp), BRCA1 (2,788 bp), and VWF (1,190 bp). All introduced gaps were treated as missing data in subsequent analyses. Alignments are available upon request.
Phylogenetic Analyses
All phylogenetic analyses were conducted under the ML and Bayesian approaches. ML was chosen because it is known to be less sensitive to potential long-branch attraction artifacts, and it takes into account the underlying molecular evolutionary process (Swofford et al. 2001
). The Bayesian approach allows the analysis of large phylogenetic data sets under complex evolutionary models (Huelsenbeck et al. 2001
). Using a likelihood framework also permits the statistical comparison of competing hypotheses and topologies obtained from different partition and rooting schemes (Huelsenbeck and Rannala 1997
; Whelan, Lio, and Goldman 2001
).
Results from Modeltest 3.06 (Posada and Crandall 1998
) based on the Akaike Information Criterion (AIC) indicate that the General Time Reversible model (GTR, Yang 1994
) plus a gamma (
) distribution of parameter
(Yang 1996a
) and a proportion I of invariable sites was the best model for each of the three data sets. However, to render ML analyses comparable between PAUP* version 4.0b8 (Swofford 1998
) and PAML version 3.0c (Yang 1997
) which does not allow the use of invariable sites, we choose the GTR +
8 model with eight discrete rate categories to represent the continuous
distribution.
In order to avoid excessive computation time, ML analyses with PAUP* version 4.0b8 (Swofford 1998
) were conducted using a loop approach to estimate the best tree and the optimal likelihood parameters. First, substitution rate matrix and among-site substitution rate heterogeneity parameters were optimized on a neighbor-joining (NJ) topology derived from ML distances obtained using a GTR +
8 model (
= 0.5). Second, a ML heuristic search was conducted by Tree Bisection Reconnection (TBR) branch swapping to identify the optimal tree under these GTR +
8 parameter estimates. Third, likelihood parameters were reestimated on this new topology. Four, a new heuristic tree search was run under the reestimated GTR +
8 parameters. This loop procedure was performed until stabilization of both topology and parameters was attained (after three cycles). The Bayesian approach to phylogenetic reconstruction (Yang and Rannala 1997
; Huelsenbeck et al. 2001
) was implemented using MrBayes 2.01 (Huelsenbeck and Ronquist 2001
). Metropolis-coupled Markov chain Monte Carlo (MCMCMC) sampling was performed with four chains that were run for 500,000 generations, using default model parameters as starting values.
The robustness of nucleotide-derived trees was estimated by Bootstrap Percentages (BP) (Felsenstein 1985
) computed by PAUP* using the optimal ML estimated parameters, with NJ starting trees and TBR branch swapping. The number of TBR rearrangements was unlimited for the Xenarthra data set (1,000 bootstrap replications) and set to 1,000 per replicate for the Placental data set (100 bootstrap replications). Bayesian posterior probabilities were picked from the 50% majority rule consensus of trees sampled every 20 generations after removing trees obtained before chains reached apparent stationarity (burn-in determined by empirical checking of likelihood values).
Evaluation of the Saturation of Nucleotide Substitutions
The nucleotide substitution saturation of the phylogenetic markers was evaluated with the graphical method used by Philippe and Forterre (1999)
. The inferred number of substitutions between each pair of sequences was estimated from the ML tree as the sum of the lengths of all branches linking these two sequences. The saturation level was estimated by plotting the number of observed differences as a function of the ML inferred number of substitutions for all 1,225 pairwise comparisons for 50 sequences. In these graphics, the straight line (Y = X) represents the situation for which the number of inferred substitutions equals the number of observed differences, i.e., there is no detected homoplasy in the data. The nucleotide substitution saturation is evidenced when the number of inferred substitutions increased, whereas the number of observed differences remained constant (plateau shape).
Statistical Tests of the Position of the Root: Impact of Character Partitions and Taxon Sampling
Extensive studies (Madsen et al. 2001
; Murphy et al. 2001a,
2001b;
this article) identified four major clades of placentals: Afrotheria, Xenarthra, Euarchontoglires, and Laurasiatheria. Assuming the respective monophyly of these cladesas suggested by the ML analyses of the datathere are 15 possible ways to connect the four clades into a bifurcating topology rooted by marsupials. These 15 alternative topologies were compared by the ML test of Kishino and Hasegawa (1989)
, with introduction of the Shimodaira and Hasegawa (1999)
correction for multiple comparisons of topologies defined a posteriori (SH test).
To evaluate the impact of character sampling on the location of the root, we performed SH tests using different character partitions of the combined data set. As ADRA2B, BRCA1, and VWF locate to different chromosomes and their protein products markedly differ in function, these nuclear genes likely evolve under differential selective pressures. To accommodate expected differences in molecular evolution, we followed the approach of Yang (1996b)
for analyzing multiple genes. Thus, each codon position of the three genes has been attributed one independent GTR +
8 model for which the different sets of likelihood parameters (base frequencies, substitution rate parameters, the
parameter of the
distribution, and branch lengths) were estimated. This allows us to take into account both differences in substitution pattern and rate heterogeneity between both the three codon positions and the three different genes. This partitioned ML model was used in PAML to compute log-likelihoods and confidence probability values for the 15 topologies corresponding to the possible locations of the root between the four major clades of placentals. The following four sets of characters were analyzed: (1) first and second codon positions of the three genes (ABV 1 + 2; 3,421 sites; 6 partitions), (2) first and second codon positions of ADRA2B and VWF plus all codon positions of BRCA1 (A2BV2; 4,350 sites; 7 partitions), (3) all codon positions of the three genes (ABV 1 + 2 + 3; 5,130 sites; 9 partitions), and (4) concatenated amino acids of the three proteins (ABV AA; 1,709 sites). Because of computational time constraints, a single JTT (Jones, Taylor, and Thornton 1992
) plus
8 model was assumed for the amino acid concatenated data set.
To assess the impact of increased taxon sampling within Xenarthra on the location of the root, we conducted SH tests after resampling of xenarthran species belonging to its three subgroups: sloths (two taxa), anteaters (three), and armadillos (eight). To do this, we built all taxon combinations with each one of the 13 available xenarthrans (13 combinations), with one from all three xenarthran subclades (48 combinations of three taxa), with two xenarthrans per subclade (84 combinations of six taxa), with three (except for sloths) xenarthrans per subclade (56 combinations of eight taxa), and with all xenarthrans (the original combination including 13 taxa). Because of computational time limitations, PAUP* was used to perform SH tests assuming a single GTR + 8 model on the different concatenated data sets.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Evolutionary Dynamics
Comparisons of ML GTR + 8 parameters estimated for each codon position reveal a particular behavior of BRCA1, relative to ADRA2B and VWF (fig. 1
). Indeed, for the two latter genes, we observed almost the same pattern of substitutions at first and second codon positions. These positions are characterized by an excess of transitions over transversions especially marked at second codon positions but seem to exhibit comparable relative evolutionary rates. Furthermore, these positions share similarly small values of the gamma shape parameter (
< 0.65) indicating strong rate heterogeneity at these positions within each gene. As expected from genetic code redundancies, third codon positions of ADRA2B and VWF evolve about three times faster than their respective first and second codon positions. High gamma shape parameter values (
> 2.92) indicate that these fast rates have a rather homogenous distribution along third codon positions of these exons. In striking contrast to this typical coding gene pattern, all codon positions of BRCA1 seem to behave similarly and evolve at about the same rate (fig. 1
). Indeed, there is almost no difference between the three codon positions in terms of the transition-transversion ratio, the gamma shape parameter, and the relative evolutionary rate. They are all characterized by high values of the gamma shape parameter (
> 2.09) and an evolutionary rate about three times higher than the slowest evolving character partition, i.e. the second codon positions of ADRA2B.
|
Nucleotide Substitution Saturation Analyses
Saturation plots of the pairwise observed differences between the 50 sequences as a function of the pairwise number of substitutions inferred on the ML tree (Philippe and Forterre 1999
) are presented in figure 2
for third codon positions of each gene. The slope of the regression line between the numbers of observed differences and inferred substitutions is an indication of the relative saturation level of the characters considered: the greater the number of inferred substitutions (abscissa) is for a given number of differences (ordinate), the higher the level of saturation. Under this graphical representation, third codon positions of BRCA1 (slope = 0.44) appear less saturated than third codon positions of ADRA2B (slope = 0.28) and VWF (slope = 0.19) for which there is considerable dispersion of pairwise comparison points. This pattern indicates that multiple substitutions are more frequent at these positions and that a beginning of saturation is reached even for comparisons between quite closely related taxa. For the 141 pairwise comparisons involving the three marsupials against the 47 placentals (triangles in fig. 2
), the plots reach a plateau for ADRA2B and VWF (slopes of 0.004 and 0.05, respectively) but not for BRCA1 (slope = 0.31). So, third codon positions of BRCA1 show a weak tendency to form a plateau, being almost linearly distributed just below the Y = X straight line. These saturation analyses strengthen the peculiarity of BRCA1 third codon positions, which appear to be less saturated than those of ADRA2B and VWF. By contrast, saturation analyses for first and second codon positions of the three genes did not reveal significant saturation, graphical patterns being very similar to the one of BRCA1 third codon positions (data not shown). Therefore, third codon positions of BRCA1 do not seem to exhibit strong saturation in nucleotide substitutions and are more likely to have retained some deep phylogenetic signal than those of ADRA2B and VWF.
|
|
|
|
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Placental Mammal Phylogeny
In many aspects, the phylogenetic picture of placental mammal relationships that we obtained is fully compatible with the challenging studies of Madsen et al. (2001)
and Murphy et al. (2001a,
2001b)
. Our study based on the detailed ML analysis at the DNA and amino acid levels of 5,130 bp of three nuclear protein-coding genes yielded conclusions consistent with those obtained on 16,397 bp by Murphy et al. (2001b)
. They strongly confirmed the existence of the four major placental clades Xenarthra, Afrotheria, Euarchontoglires, and Laurasiatheria and add support to the grouping of the Northern Hemisphere clades Euarchontoglires and Laurasiatheria in the so-called Boreoeutheria (Springer and de Jong 2001
).
Xenarthra
We present here the first comprehensive study of xenarthran molecular phylogeny based on the analysis of three nuclear markers for 12 of the 13 living genera. Analyses of the placental data sets (fig. 4
) unambiguously supported the monophyly of the order Xenarthra and its division in the two suborders Cingulata (armadillos) and Pilosa (anteaters and sloths). This result confirms those obtained on both morphological (Engelmann 1985
; Patterson et al. 1992
; Gaudin 1999
) and molecular (van Dijk et al. 1999
; Delsuc et al. 2001
; Madsen et al. 2001
; Murphy et al. 2001a
) characters. Our results also led to a clear picture of Xenarthra interrelationships with almost all nodes robustly evidenced (fig. 3
).
Relationships within anteaters and sloths (Pilosa) are fully congruent with previous analyses (Delsuc et al. 2001
) and confirm both the respective monophyly of anteaters (Vermilingua) and sloths (Folivora), and the early emergence of the pygmy anteater (Cyclopes) within Vermilingua (fig. 3
). These results emphasize the very deep split between larger body sized anteaters (Tamandua and Myrmecophaga) and the pygmy Cyclopes reflected by numerous morphological peculiarities related to its strictly arboreal way of life (Gaudin and Branham 1998
).
Regarding armadillos (Cingulata, Dasypodidae), our results clearly identify three distinct lineages corresponding to the subfamilies defined by McKenna and Bell (1997, pp. 82102)
: Dasypodinae (D. kappleri and D. novemcinctus), Euphractinae (Chaetophractus, Euphractus, and Zaedyus) and Tolypeutinae (Priodontes, Cabassous, and Tolypeutes). The evidence strongly supports the early emergence of Dasypodinae, whereas Tolypeutinae and Euphractinae unequivocally cluster together (fig. 3
). Such a relationship has already been proposed based on the study of spermatozoa (Cetica et al. 1998
) and contradicts the early morphological assessment of Engelmann (1985)
. However, relationships within the subfamilies Euphractinae and Tolypeutinae remain unclear. Our results suggest the grouping of Euphractus and Zaedyus excluding Chaetophractus within Euphractinae. Such a grouping is consistently retrieved in most of the analyses but appears with only poor support in the total evidence ML tree (fig. 3
). This lack of resolution is not unexpected because the three genera are morphologically so similar that their interrelationships have always been considered as unresolved (Engelmann 1985
; Patterson, Segall, and Turnbull 1989
). More surprising is the trifurcation obtained within Tolypeutinae. Indeed, a close relationship between Cabassous and Priodontes was expected on the basis of their comparable external morphologies (Engelmann 1985
; Wetzel 1985
) and spermatozoa (Cetica et al. 1998
). By contrast, in our ML phylogeny, Priodontes emerges in a basal position relative to Cabassous and Tolypeutes but with low support (fig. 3
). Thus, relationships between these three genera are still unclear. However, Tolypeutes possesses a highly distinct morphology shaped by the anatomical constraints induced by its capacity to entirely roll into a ball, hence rendering the interpretation of homology in postcranial characters difficult. In this context, it is interesting to note that characters of the ear region tend to support a sister-group relationship between Tolypeutes and Cabassous as proposed by our data (Patterson, Segall, and Turnbull 1989
). The apparent lack of resolution observed within these two clades despite the high number of nucleotides analyzed raises the question of the rapidity by which the splitting events occurred. The shortness of internal nodes associated with low-support values suggest that rapid cladogenesis left only short time intervals for molecular synapomorphies to accumulate in these two groups. It is hoped that considering more rapidly evolving molecules, such as nuclear introns or mitochondrial genes, will help to resolve the only two remaining uncertainties in Xenarthra phylogeny.
Afrotheria
Despite its highly provocative nature for morphologists (Novacek 2001
), there is at present little doubt about the naturalness of the African clade for which support consistently increases as molecular evidences accumulate (de Jong, Leunissen, and Wistow 1993
; Douzery and Catzeflis 1995
; Lavergne et al. 1996
; Springer et al. 1997
, 1999
; Stanhope et al. 1998a,
1998b;
Mouchaty et al. 2000b;
Madsen et al. 2001
; Murphy et al. 2001a,
2001b;
van Dijk et al. 2001
). Our results provide strong ML bootstrap support for nodes within Afrotheria that were previously difficult to resolve (Madsen et al. 2001
; Murphy et al. 2001a
). Indeed, whereas the retrieval of the well-defined Paenungulata (elephants, sirenians, and hyraxes) was expected, a second major afrotherian clade regrouping Tubulidentata (aardvark), Macroscelidea (elephant shrews), Chrysochloridae (golden moles), and Tenrecidae (tenrecs) is robustly confirmed. Within this clade, the grouping of the two exinsectivoran families Chrysochloridae and Tenrecidae is also well supported. Moreover, in agreement with Murphy et al. (2001b)
, our results support the early emergence of Tubulidentata leading to the monophyly of an African insectivore-like clade (Macroscelidea + Afrosoricida). It is finally worth noting that, as in previous analyses of concatenated nuclear genes (Madsen et al. 2001
; Murphy et al. 2001a,
2001b
), relationships among paenungulates remain unsettled despite the high numbers of nucleotides analyzed. The resolution of this tricky question and the verification of the proposed sister-group relation between Macroscelidea and Afrosoricida might benefit from an expanding taxon sampling within Afrotheria.
Boreoeutheria
The association of Euarchontoglires and Laurasiatheria appears strongly supported in all ML analyses (fig. 4
), and its monophyly is suggested by most ML statistical tests (table 2
). Such an arrangement, fundamentally opposing Northern and Southern placental groups, argues for a more crucial importance of plate tectonic events in the early diversification of placental mammals than previously recognized (Eizirik, Murphy, and O'Brien 2001
; Madsen et al. 2001
; Murphy et al. 2001a,
2001b
). It also reveals the occurrence of parallel adaptive radiations on separated continental masses that are likely responsible for the long-standing difficulties encountered by systematic morphologists to reconstruct placental mammal evolution (Madsen et al. 2001
; Scally et al. 2001
).
The Euarchontoglires clade containing the orders Primates, Scandentia (tree shrews), Dermoptera (flying lemurs), Lagomorpha (rabbits, hares and pikas) and Rodentia (rodents) is consistently retrieved in all analyses with high support values (fig. 4
). The monophyly of rodents is unambiguously supported and their association with lagomorphs into Glires received moderate support. This arrangement, which has for a long time been favored by morphological studies (Gregory 1910
; Simpson 1945
), has proved difficult to retrieve in molecular studies, but was finally confirmed by larger taxon sampling (Murphy et al. 2001a;
Huchon et al. 2002
). However, it has been questioned whether the Murphy et al. (2001a)
data set is adequate to resolve the relationships of rabbits, rodents, and primates or needs more genes with longer sequences (Rosenberg and Kumar 2001
).
Similar to Afrotheria and Xenarthra, the clade Laurasiatheria (Waddell, Okada, and Hasegawa 1999
) is always supported by 100% bootstrap support values or Bayesian posterior probabilities of 1.00 in our analyses. The existence of this clade is corroborated by analyses of complete mitochondrial genomes, in which the basal position of the hedgehog is likely considered to be an artifact because of its peculiar nucleotide composition (Cao et al. 2000
; Mouchaty et al. 2000b;
Nikaido et al. 2001
). Within this clade, the proposal of an early divergence of Eulipotyphla from the remaining laurasiatherians (Waddell, Okada, and Hasegawa 1999
; Murphy et al. 2001b
) finds support in our data. The monophyly of Chiroptera is unambiguously supported by our data as is the paraphyly of Microchiroptera (Hipposideros, Megaderma, Tonatia, Myotis, Tadarida) and its subsequent implications for the evolution of echolocation and flight in bats (Hutcheon, Kirsch, and Pettigrew, 1998
; Teeling et al. 2000
, 2002
; Springer et al. 2001
). An interesting input provided by our data resides in the reliable support obtained for the sister-group relationship of Pholidota (pangolins) and Carnivora (fig. 4
). This grouping, called Ferae by Shoshani and McKenna (1998)
, was suggested by early molecular studies (de Jong et al. 1985
; Shoshani 1986
; de Jong, Leunissen, and Wistow 1993
). In the light of our results and those of Murphy et al. (2001b)
, it seems likely that the presence of a well-developed ossified tentorium in both carnivores and pangolins represents a true synapomorphy uniting these two orders (Shoshani and McKenna 1998
). Cetartiodactyla is indubitably monophyletic based on our data (fig. 4
) and other multiple sets of molecular evidence (Montgelard, Catzeflis, and Douzery 1997
; Shimamura et al. 1997
; Gatesy et al. 1999
). In addition, the artiodactyl affinities of whales have been confirmed by recent fossil discoveries (Gingerich et al. 2001
; Thewissen et al. 2001
). Among Cetartiodactyla, the relationships are entirely congruent with those obtained from studies of other genes (Montgelard, Catzeflis, and Douzery 1997
; Ursing and Arnason 1998
; Gatesy et al. 1999
; Madsen et al. 2001
; Murphy et al. 2001a,
2001b
) and SINE insertions (Shimamura et al. 1997
; Nikaido, Rooney, and Okada 1999
). Nevertheless, relationships among these five major Laurasiatherian clades remain unclear.
Character and Taxon Sampling and the Root of the Placental Tree
Despite extensive taxon sampling, our data did not allow the exact position of the root of the placental tree to be determined. Indeed, like Madsen et al. (2001)
and Murphy et al. (2001a)
, we identified three almost equally likely possible root locations corresponding to the following topological arrangements: basal Afrotheria, basal Xenarthra + Afrotheria, and basal Xenarthra (as sister-group to Epitheria). Recently, Murphy et al. (2001b)
proposed the root along the branch leading to Afrotheria. One of our goals was to evaluate to what extent increasing the taxon representation of a previously poorly sampled basal clade (Xenarthra) might help to stabilize the position of the root on the placental tree.
Although not statistically conclusive, a detailed examination of the impact of both character and taxon sampling on root location tests (SH tests) allows for the identification of clear trends. Thus, it is worth noting that Epitheria never appears as the most likely hypothesis, suggesting that Xenarthra does not represent the earliest offshoot of the placental tree as previously thought (McKenna 1975
). The location of the root for the two remaining hypotheses is sensitive to the molecular characters sampled. Indeed, the basal Afrotheria + Xenarthra hypothesis tends to be slightly favored when only first and second codon positions are used, whereas the basal Afrotheria solution is preferred when adding BRCA1 third codon positions or third codon positions of the three genes. It is however not clear whether this results from signals caused by the growing number of sampled characters (3,4215,130) or to noise coming from some saturated third codon positions. Moreover, it is striking to note that in amino acid analyses the three competing hypotheses for the position of the root (basal Afrotheria, Afrotheria + Xenarthra, and Epitheria) only differ by no more than 0.75 units of log-likelihood (table 2 ). Increasing the number of genes that are analyzed, in association with a critical examination of different codon or amino acid partitions, will be important in further investigations of the rooting problem.
It is also important to maintain adequate taxon sampling. Indeed, our results show that including only a single xenarthran representative could be misleading in terms of SH test results, especially when the number of analyzed characters increased. Actually, considering at least one xenarthran per major lineage (armadillos, anteaters, and sloths) might help to stabilize their position in the tree. It is likely that increasing the sampling within Xenarthra contributed toward reducing the effect of potential base composition biases or long branch attraction artifacts associated with the distantly related marsupial outgroups. We suggest here that additional sampling among Marsupials and the addition of a representative of Monotremes might also contribute toward attaining more reliability for assessing the root in placentals. Larger taxon sampling for different genes might help to find a reliable rooting of the ingroup by mitigating the peculiarities of certain taxa and thus increasing the signal-noise ratio (Hillis 1996
; Rannala et al. 1998
). It is hoped that the future consideration of selected positions of additional genes for a balanced taxon sampling of major clades will resolve the remaining uncertainties about the root position of the placental tree. This question is of primary interest for understanding and reconstructing the early history of placental mammals, in which tectonic and biogeographic events seem to have played a crucial role (Murphy et al. 2001b
).
![]() |
Supplementary Material |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: xenarthran phylogeny
taxon and character sampling
placentals
evolutionary dynamics
nuclear markers
maximum likelihood and Bayesian phylogenetics
Address for correspondence and reprints: Emmanuel J. P. Douzery, Laboratoire de Paléontologie, Paléobiologie et PhylogénieCC064, Institut des Sciences de l'Evolution UMR 5554/CNRS, Université Montpellier II, Place Eugène Bataillon, 34 095 Montpellier Cedex 05, France. douzery{at}isem.univ-montp2.fr
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Adkins R. M., E. L. Gelke, D. Rowe, R. L. Honeycutt, 2001 Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes Mol. Biol. Evol 18:777-791
Arnason U., A. Gullberg, A. Janke, 1997 Phylogenetic analyses of mitochondrial DNA suggest a sister group relationship between Xenarthra (Edentata) and Ferungulates Mol. Biol. Evol 14:762-768[Abstract]
Bernardi G., 2001 Misunderstandings about isochores. Part 1 Gene 276:3-13[ISI][Medline]
Cao Y., M. Fujiwara, M. Nikaido, N. Okada, M. Hasegawa, 2000 Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data Gene 259:149-158[ISI][Medline]
Catzeflis F., 1991 Animal tissue collections for molecular genetics and systematics Trends Ecol. Evol 6:168.[ISI]
Cetica P. D., A. J. Solari, M. S. Merani, J. C. De Rosas, M. H. Burgos, 1998 Evolutionary sperm morphology and morphometry in armadillos J. Submicroscop. Cytol. Pathol 30:309-314[ISI][Medline]
de Jong W. W., J. A. M. Leunissen, G. J. Wistow, 1993 Eye lens crystallins and the phylogeny of placental orders: evidence for a Macroscelid-Paenungulate clade? Pp. 512 in M. J. N. a. M. C. M. E. F. S. Sazlay, ed. Mammal phylogeny (placentals). Springer-Verlag, New York
de Jong W. W., A. Zweers, K. A. Joysey, J. T. Gleaves, D. Boulter, 1985 Protein sequence analysis applied to Xenarthran and Pholidote phylogeny Pp. 6576 in G. G. Montgomery, ed. The evolution and ecology of armadillos, sloths and vermilinguas. Smithsonian Institution, Washington, DC
Delsuc F., F. M. Catzeflis, M. J. Stanhope, E. J. Douzery, 2001 The evolution of armadillos, anteaters and sloths depicted by nuclear and mitochondrial phylogenies: implications for the status of the enigmatic fossil Eurotamandua Proc. R. Soc. Lond. B Biol. Sci 268:1605-1615.[ISI][Medline]
Deng C.-X., S. G. Brodie, 2000 Roles of BRCA1 and its interacting proteins Bioessays 22:728-737[ISI][Medline]
Douzery E., F. M. Catzeflis, 1995 Molecular evolution of the mitochondrial 12S rRNA in Ungulata (Mammalia) J. Mol. Evol 41:622-636[ISI][Medline]
Eizirik E., W. J. Murphy, S. J. O'Brien, 2001 Molecular dating and biogeography of the early placental mammal radiation J. Hered 92:212-219
Engelmann G. F., 1985 The phylogeny of the Xenarthra Pp. 5163 in G. G. Montgomery, ed. The evolution and ecology of armadillos, sloths and vermilinguas. Smithsonian Institution, Washington, DC
Felsenstein J., 1985 Confidence limits on phylogenies: an approach using the bootstrap Evolution 39:783-791[ISI]
Gatesy J., M. C. Milinkovitch, V. G. Waddell, M. S. Stanhope, 1999 Stability of cladistic relationships between Cetacea and higher-level Artiodactyl taxa Syst. Biol 48:6-20[ISI][Medline]
Gaudin T. J., 1999 The morphology of xenarthrous vertebrae (Mammalia: Xenarthra) Fieldiana: Geology, New Ser 41:1-38.
Gaudin T. J., D. G. Branham, 1998 The phylogeny of the Myrmecophagidae (Mammalia, Xenarthra, Vermilingua) and the relationship of Eurotamandua to the Vermilingua J. Mamm. Evol 5:237-265
Gaudin T. J., J. R. Wible, J. A. Hopson, W. D. Turnbull, 1996 Reexamination of the morphological evidence for the Cohort Epitheria (Mammalia, Eutheria) J. Mamm. Evol 3:31-78
Gingerich P. D., M. Haq, I. S. Zalmout, I. H. Khan, M. S. Malkani, 2001 Origin of whales from early artiodactyls: hands and feet of Eocene Protocetidae from Pakistan Science 293:2239-2242
Greenwood A. D., J. Castresana, G. Feldmaeir-Fuchs, S. Pääbo, 2001 A molecular phylogeny of two extinct sloths Mol. Phylogenet. Evol 18:94-103[ISI][Medline]
Gregory W. K., 1910 The orders of mammals Bull. Am. Mus. Nat. Hist 27:1-524
Guth C., 1961 La région temporale des Edentés Doctoral dissertation, University of Paris, Paris, France
Hillis D. M., 1996 Inferring complex phylogenies Nature 383:130-131[ISI][Medline]
Huchon D., O. Madsen, M. J. J. B. Sibbald, K. Ament, M. J. Stanhope, F. Catzeflis, W. W. de Jong, E. J. P. Douzery, 2002 Rodent phylogeny and a timescale for the evolution of Glires: evidence from an extensive taxon sampling using three nuclear genes Mol. Biol. Evol. 19:1053-1065
Huelsenbeck J. P., B. Rannala, 1997 Phylogenetic methods come of age: testing hypotheses in an evolutionary context Science 276:227-32
Huelsenbeck J. P., F. Ronquist, 2001 MRBAYES: Bayesian inference of phylogenetic trees Bioinformatics 17:754-755
Huelsenbeck J. P., F. Ronquist, R. Nielsen, J. P. Bollback, 2001 Bayesian inference of phylogeny and its impact on evolutionary biology Science 294:2310-2314
Hutcheon J. M., J. A. Kirsch, J. D. Pettigrew, 1998 Base-compositional biases and the bat problem. III. The questions of microchiropteran monophyly Philos. Trans. R. Soc. Lond. B Biol. Sci 353:607-617.[ISI][Medline]
Jones D. T., W. R. Taylor, J. M. Thornton, 1992 The rapid generation of mutation data matrices from protein sequences Comput. Appl. Biosci 8:275-282[Abstract]
Kishino H., M. Hasegawa, 1989 Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in Hominoidea J. Mol. Evol 29:170-179[ISI][Medline]
Lavergne A., E. Douzery, T. Stichler, F. M. Catzeflis, M. S. Springer, 1996 Interordinal mammalian relationships: evidence for Paenungulate monophyly is provided by complete mitochondrial 12S rRNA sequences Mol. Phylogenet. Evol 6:245-258[ISI][Medline]
Lecointre G., H. Philippe, H. L. vân Lê, H. Le Guyader, 1993 Species sampling has a major impact on phylogenetic inference Mol. Phylogenet. Evol 2:205-224[Medline]
Loughry W. J., P. A. Prodöhl, C. M. McDonough, J. C. Avise, 1998 Polyembryony in armadillos Am. Sci 86:274-279[ISI]
Madsen O., M. Scally, C. J. Douady, D. J. Kao, R. W. DeBry, R. Adkins, H. M. Amrine, M. J. Stanhope, W. W. de Jong, M. S. Springer, 2001 Parallel adaptative radiations in two major clades of placental mammals Nature 409:610-614[ISI][Medline]
McKenna M. C., 1975 Toward a phylogenetic classification of the Mammalia Pp. 2146 in W. P. Luckett and F. S. Szalay, eds. Phylogeny of the primates. Plenum Press, New York
McKenna M. C., S. K. Bell, 1997 Classification of mammals above the species level Columbia University Press, New York
Montgelard C., F. M. Catzeflis, E. Douzery, 1997 Phylogenetic relationships of artiodactyls and Ceteacans as deduced from the comparison of cytochrome b and 12S rRNA mitochondrial sequences Mol. Biol. Evol 14:550-559[Abstract]
Mouchaty S. K., A. Gullberg, A. Janke, U. Arnason, 2000a. Phylogenetic position of the Tenrecs (Mammalia: Tenrecidae) of Madagascar based on analysis of the complete mitochondrial genome sequence of Echinops telfairi Zool. Scr 29:307-317[ISI]
. 2000b. The phylogenetic position of the Talpidae within eutheria based on analysis of complete mitochondrial sequences Mol. Biol. Evol 17:60-67
Murphy W. J., E. Eizirik, W. E. Jonhson, Y. P. Zhang, O. A. Ryder, S. J. O'Brien, 2001a. Molecular phylogenetics and the origins of placental mammals Nature 409:614-618[ISI][Medline]
Murphy W. J., E. Eizirik, S. J. O'Brien, et al. (11 co-authors) 2001b. Resolution of the early placental mammal radiation using Bayesian phylogenetics Science 294:2348-2351
Nikaido M., K. Kawai, Y. Cao, M. Harada, S. Tomita, N. Okada, M. Hasegawa, 2001 Maximum likelihood analysis of the complete mitochondrial genomes of eutherians and a reevaluation of the phylogeny of bats and insectivores J. Mol. Evol 53:508-516[ISI][Medline]
Nikaido M., A. P. Rooney, N. Okada, 1999 Phylogenetic relationships among cetartiodactyls based on insertions of short and long interpersed elements: hippopotamuses are the closest extant relatives of whales Proc. Natl. Acad. Sci. USA 96:10261-10266
Novacek M. J., 1992 Mammalian phylogeny: shaking the tree Nature 356:121-125[ISI][Medline]
. 2001 Mammalian phylogeny: genes and supertrees Curr. Biol 11:R573-R575[ISI][Medline]
Patterson B., R. Pascual, 1972 The fossil mammal fauna of South-America Pp. 247309 in A. Keast, F. C. Erk, and B. P. Glass, eds. Evolution, mammals and southern continents. State University of New-York Press, Albany, NY
Patterson B., W. Segall, W. D. Turnbull, 1989 The ear region in xenarthrans (= Edentata, Mammalia). Part I Cingulates. Fieldiana: Geology, New Ser 18:1-46.
Patterson B., W. Segall, W. D. Turnbull, T. J. Gaudin, 1992 The ear region in xenarthrans (= Edentata, Mammalia). Part II. Sloths, Anteaters, Palaeanodonts, and a Miscellany Fieldiana: Geology, New Ser 24:1-79.
Philippe H., 1993 MUST: a computer package of management utilities for sequences and trees Nucleic Acids Res 21:5264-5272[Abstract]
Philippe H., P. Forterre, 1999 The rooting of the universal tree of life is not reliable J. Mol. Evol 49:509-523[ISI][Medline]
Poe S., D. L. Swofford, 1999 Taxon sampling revisited Nature 398:299-300[Medline]
Posada D., K. A. Crandall, 1998 MODELTEST: testing the model of DNA substitution Bioinformatics 19:817-818
Rannala B., J. P. Huelsenbeck, Z. H. Yang, R. Nielsen, 1998 Taxon sampling and the accuracy of large phylogenies Syst. Biol 47:702-710[ISI][Medline]
Rose K. D., R. J. Emry, 1993 Relationships of Xenarthra, Pholidota, and Fossil "Edentates": the morphological evidence Pp. 81102 in F. S. Sazlay, M. J. Novacek, and M. C. McKenna, eds. Mammal phylogeny (placentals). Springer-Verlag, New York
Rosenberg M. S., S. Kumar, 2001 Incomplete taxon sampling is not a problem for phylogenetic inference Proc. Natl. Acad. Sci. USA 98:10751-6
Scally M., O. Madsen, C. J. Douady, W. W. de Jong, M. J. Stanhope, M. S. Springer, 2001 Molecular evidence for the major clades of placental mammals J. Mamm. Evol 8:239-277
Shimamura M., H. Yasue, K. Ohshima, H. Abe, H. Kato, T. Kishiro, M. Goto, I. Munechika, N. Okada, 1997 Molecular evidence from retroposons that whales form a clade within even-toed ungulates Nature 388:666-670[ISI][Medline]
Shimodaira H., M. Hasegawa, 1999 Multiple comparisons of log-likelihoods with applications to phylogenetic inference Mol. Biol. Evol 16:1114-1116
Shoshani J., 1986 Mammalian phylogeny: comparison of morphological and molecular results Mol. Biol. Evol 3:222-242[Abstract]
Shoshani J., M. C. McKenna, 1998 Higher taxonomic relationships among extant Mammals based on morphology, with selected comparisons of results from molecular data Mol. Phylogenet. Evol 9:572-584[ISI][Medline]
Simpson G. G., 1945 The principles of classification and a classification of mammals Bull. Am. Mus. Nat. Hist 85:1-350[ISI]
Springer M. S., H. M. Amrine, A. Burk, M. J. Stanhope, 1999 Additional support for Afrotheria and Paenungulata, the performance of mitochondrial versus nuclear genes, and the impact of data partitions with heterogenous base composition Syst. Biol 48:65-75[ISI][Medline]
Springer M. S., G. C. Cleven, O. Madsen, W. W. de Jong, V. G. Waddell, H. M. Amrine, M. J. Stanhope, 1997 Endemic African mammals shake the phylogenetic tree Nature 388:61-64[ISI][Medline]
Springer M. S., R. W. DeBry, C. Douady, H. M. Amrine, O. Madsen, W. W. de Jong, M. J. Stanhope, 2001 Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction Mol. Biol. Evol 18:132-143
Springer M. S., W. W. de Jong, 2001 Which mammalian supertree to bark up? Science 291:1709-1711
Springer M. S., E. C. Teeling, O. Madsen, M. J. Stanhope, W. W. de Jong, 2001 Integrated fossil and molecular data reconstruct bat echolocation Proc. Natl. Acad. Sci. USA 98:6241-6246
Stanhope M. J., O. Madsen, V. G. Waddell, G. C. Cleven, W. W. de Jong, M. S. Springer, 1998b. Highly congruent molecular support for a diverse superordinal clade of endemic African mammals Mol. Phylogenet. Evol 9:501-508[ISI][Medline]
Stanhope M. J., V. G. Waddell, O. Madsen, W. W. de Jong, S. Blair Hedges, G. C. Cleven, D. Kao, M. S. Springer, 1998a. Molecular evidence for multiple origins of Insectivora and for a new order of endemic African insectivore mammals Proc. Natl. Acad. Sci. USA 95:9967-9972
Strimmer K., A. von Haeseler, 1996 Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies Mol. Biol. Evol 13:964-969
Swofford D. L., 1998 PAUP*. Phylogenetic analysis using parsimony (* and other methods). Version 4.0b8 Sinauer Associates, Sunderland, Mass
Swofford D. L., P. J. Waddell, J. P. Huelsenbeck, P. G. Foster, P. O. Lewis, J. S. Rogers, 2001 Bias in phylogenetic estimation and its relevance to the choice between parsimony and likelihood methods Syst. Biol 50:525-539[ISI][Medline]
Teeling E. C., O. Madsen, R. A. Van Den Bussche, W. W. de Jong, M. J. Stanhope, M. S. Springer, 2002 Microbat paraphyly and the convergent evolution of a key innovation in Old World rhinolophoid microbats Proc. Natl. Acad. Sci. USA 99:1431-1436
Teeling E. C., M. Scally, D. J. Kao, M. L. Romagnoli, M. S. Springer, M. J. Stanhope, 2000 Molecular evidence regarding the origin of echolocation and flight in bats Nature 403:188-192[ISI][Medline]
Thewissen J. G., E. M. Williams, L. J. Roe, S. T. Hussain, 2001 Skeletons of terrestrial cetaceans and the relationship of whales to artiodactyls Nature 413:277-281[ISI][Medline]
Ursing B. M., U. Arnason, 1998 Analyses of mitochondrial genomes strongly support a hippopotamus-whale clade Proc. R. Soc. Lond. B 265:2251-2255[ISI][Medline]
van Dijk M. A., O. Madsen, F. Catzeflis, M. J. Stanhope, W. W. de Jong, M. Pagel, 2001 Protein sequence signatures support the African clade of mammals Proc. Natl. Acad. Sci. USA 98:188-193
van Dijk M. A. M., E. Paradis, F. Catzeflis, W. W. de Jong, 1999 The virtues of gaps: Xenarthran (Edentate) monophyly supported by a unique deletion in alphaA-crystallin Syst. Biol 48:94-106[ISI][Medline]
Vizcaíno S. F., 1995 Identificación específica de las "mulitas", género Dasypus L. (Mammalia, Dasypodidae), del noroeste argentino Descripción de una nueva espécie. Mastozool. Neotrop 2:5-13
Waddell P. J., Y. Cao, J. Hauf, M. Hasegawa, 1999 Using novel phylogenetic methods to evaluate mammalian mtDNA, including amino acid-invariant sites-logDet plus site stripping, to detect internal conflicts in the data, with special reference to the positions of hedgehog, armadillo, and elephant Syst. Biol 48:31-53[ISI][Medline]
Waddell P. J., N. Okada, M. Hasegawa, 1999 Towards resolving the interordinal relationships of placental mammals Syst. Biol 48:1-5[ISI][Medline]
Wetzel R. M., 1985 The identification and distribution of recent Xenarthra (= Edentata) Pp. 521 in G. G. Montgomery, ed. The evolution and ecology of armadillos, sloths and vermilinguas. Smithsonian Institution, Washington, DC
Whelan S., P. Lio, N. Goldman, 2001 Molecular phylogenetics: state-of-the-art methods for looking into the past Trends Genet 17:262-272[ISI][Medline]
Yang Z., 1994 Estimating the pattern of nucleotide substitution J. Mol. Evol 39:105-111[ISI][Medline]
. 1996a. Among-site rate variation and its impact on phylogenetic analyses Trends Ecol. Evol 11:367-372[ISI]
. 1996b. Maximum-likelihood models for combined analyses of multiple sequence data J. Mol. Evol 42:587-596[ISI][Medline]
. 1997 PAML: phylogenetic analysis by maximum likelihood University College London, London
Yang Z., B. Rannala, 1997 Bayesian phylogenetic inference using DNA sequences: a Markov chain Monte Carlo method Mol. Biol. Evol 14:717-724[Abstract]