*Primate Genetics, German Primate Center, Göttingen;
Department of Biology, Faculty of Mathematics and Natural Sciences, Bogor Agricultural University, Bogor, Indonesia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The phylogenetic analysis of the complete mitochondrial (mt) genome sequence of Cynocephalus variegatus presented in this study first appears to corroborate interpretations of primates as a paraphyletic group. However, more detailed analyses disclosed that mt nucleotide composition and consequently amino acid (AA) composition varied considerably among the species analyzed. This led us to assume that the flying lemur may be incorrectly grouped with anthropoids on the basis of similar mt nucleotide and AA compositions, rather than reflecting the true evolutionary relationship.
To reanalyze the flying lemur's evolutionary association with other eutherian orders from a completely different molecular perspective, a molecular cladistic approach was applied. To this end, we determined the presence/absence pattern of transposable elements that provide a nearly homoplasy-free and copious source of molecular evolutionary markers, with well-defined character polarity. We could identify transposable elements, both on a multilocus and single-locus level, being present in all extant primate infraorders but absent in the flying lemur, thus clearly supporting the monophyly of primates by retropositional evidence.
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Using molecular evidence, several authors have recently split the eutherian orders into four major groups, displaying primates in a cluster with rodents, lagomorphs, tree shrews, and flying lemurs (Madsen et al. 2001
; Murphy et al. 2001a
, 2001b
). Although the positioning of the major groups is largely confirmed, both by nuclear and mitochondrial (mt) DNA data sets, the phylogenetic affiliations among the members of the major groups remain, in part, contradictory. This is illustrated by the varying affiliations of colugos to other eutherian taxa, as described in Murphy et al. (2001a
, 2001b)
. Although it is difficult to compare distinct topologies because of the different taxonomic sampling, it is clear that primate affiliations to other eutherians have not been resolved to a degree that would permit the launch of large-scale comparative projects, e.g., on molecular character evolution.
Another difficulty in establishing a firm phylogenetic framework for linking primates to other eutherians is represented by the scanty early primate fossil record, so that the discovery of new fossils often requires major revisions of phylogenetic hypotheses (Beard 1990
, 1993
; Kay, Thorington, and Houde 1990
; Martin 1993
).
The concept of the superorder Archonta, first proposed by Gregory (1910)
and revived in the last decades (Novacek 1992
), indicates a close relationship between primates and dermopterans. This superorder comprises primates, colugos, bats, and tree shrews. Although the monophyly of the Archonta is partially refuted both by molecular (Schmitz, Ohme, and Zischler 2000
; Teeling et al. 2000
; Lin and Penny 2001
) and by palaeontological (Kay, Thorington, and Houde 1990
) evidence, it is commonly agreed that the primate origin lies somewhere among archontan representatives (for an overview see Fleagle and MacPhee 1993
, pp. 1383).
To complement early primate fossil records with molecular evidence, we analyzed the complete mtDNA of Cynocephalus variegatus. Results were combined with retropositional data for testing different phylogenetic interpretations in primates and other eutherians (for GenBank accession numbers see Supplementary Material).
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Sequence alignments were carried out by CLUSTAL X (Thompson et al. 1997
). Phylogenetic analyses of both DNA and amino acids (AA) were performed using maximum parsimony (MP, heuristic search), maximum likelihood (ML, with eight categories of the discrete, gamma-distribution model of rate heterogeneity over sites, mtREV model of sequence evolution), and distance-based (neighbor-joining, LogDet distance measure, either including or excluding constant sites, or Dayhoff PAM matrix) methods implemented in PAUP* 4.0b8 (Swofford 2000
), PHYLIP 3.6 (Felsenstein 1995
), and TREE-PUZZLE 5.0 (Strimmer and von Haeseler 1996
). All bootstrap and puzzling supports are based on 1,000 replicates.
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Because this result is at odds with previous hypotheses, we used an entirely different class of phylogenetic markers, performing single-locus and multilocus analyses of primate-specific short interspersed nuclear element (SINE) transpositions (for review see Shedlock and Okada 2000
). This approach makes use of a peculiar feature of the primate genome, which is the high abundance of dimerized, seven SL-RNAderived SINEs or Alu-sequences. Because of a wavelike spreading of Alu-SINEs during primate evolution, Alu-SINEs cover more than 10% of the entire human genome (Li et al. 2001
) and apparently also considerable parts of the genome of other primates (Hwu et al. 1986
). Given the irreversible integration of individual Alu-SINEs at mainly random genomic target sites, they are ideally suited as molecular cladistic markers and landmarks of primate evolution (Schmitz, Ohme, and Zischler 2001
; Singer et al. 2003
).
In Southern blot experiments using primate Alu elements as a probe, we challenged HinfI-restricted total DNA of humans and lemursthe two primates representing the descendants of the deepest split in the primate treealong with tree shrew, colugo, bat, rabbit, hedgehog, and mouse (fig. 2A ). Hybridization conditions were chosen to detect hybrids with at least 70% sequence similarity between the genomic DNAs and an Alu SINE probe. Only the two primate representatives revealed a signal indicating high copy numbers of the Alu-element. No signal was detected in the colugo, as would be expected with colugos as a sister group of the Anthropoidea, as indicated by the mtDNA-based phylogenetic reconstructions.
|
Finally, we analyzed three haplorhine-specific Alu SINE markers published earlier (Schmitz, Ohme, and Zischler 2001
), which unequivocally cluster together tarsiers and anthropoids as sister groups. All were absent in C. variegatus (as an example see fig. 2C
). Cynocephalus variegatus thus displays the ancestral character state at all analyzed orthologous loci that became targets for retroposon integration at various times during primate evolution.
The independent retropositional evidence clearly demonstrates that the apparently significant affiliation of colugos to higher primates obtained in mtDNA-based phylogenetic reconstructions is due to a misleading mtDNA signal. We consider this effect sufficiently strong to also suggest a close relationship between colugos and higher primates in data sets composed of both nuclear and mt sequences (see fig. 1
of Murphy et al. 2001a
).
We recently presented evidence that mt-specific, directional mutational pressure might be responsible for the artifact of tarsiers and strepsirhines appearing as sister groups on account of high bootstrap values. Depending on the algorithm of reconstruction, they even fall together with nonprimate representatives (Schmitz, Ohme, and Zischler 2002
). In this study we propose that the same mechanism inversely places the colugo amidst the primates. Directional mutation pressure leads to an alteration of the nucleotide composition; hence, in extreme cases, phylogenetic reconstructions might include an artificial clustering of phylogenetically unrelated taxa due to similar base composition (Tarrio, Rodriguez-Trelles, and Ayala 2001
). This is obviously most pronounced where two successive splits occur within a short period of time, leaving behind only weak phylogenetic signals. In an analysis of mt genes and their respective nuclear pseudogenes, we were able to show that in primates the plasticity of nucleotide composition is specific for mtDNA (Schmitz, Ohme, and Zischler 2002
).
We compared the nucleotide composition for several mammals (fig. 3A ). Cynocephalus variegatus clearly represents the nucleotide composition exclusively found in higher primates. Using correspondence analysis to examine the AA composition of C. variegatus mtDNA, a multivariate method in which observations (AA composition) and variables (taxa) can be jointly displayed in a two-dimensional space, a clustering of C. variegatus and higher primates was revealed, disclosing a similar AA preference in these taxa. In contrast, the prosimians (tarsiers and stepsirhines) group with nonprimate eutherians, exhibiting an affinity for AAs encoded by AT-rich codons (circled in fig. 3B ).
|
To circumvent the problem of variable base composition among sequences, we applied the LogDet transformation by both including and excluding invariant sites in separate analyses. However, this method did not resolve the conflicting tree shown in figure 1 . One potential reason for this could be that LogDet does not compensate for rate heterogeneity across sites (see also Waddell et al. 1999
).
Thus, from an ad hoc analysis of the mtDNA of mammalian representatives, we cannot gain conclusive evidence for a sister group relationship of dermopterans and primates. To avoid an inadvertent effect of the nucleotide composition on phylogenetic reconstructions, we excluded all higher primates from further computations, thus taking into account only taxa with similar mt nucleotide composition. As a result, independent of the method of phylogenetic reconstruction applied, the mtDNA sequence of colugo supports a sister taxon relationship to tarsiers and strepsirhines as representatives of the primate order (tree not shown).
In summary, the SINE results clearly refute both the mtDNA-based phylogenetic position of the colugo presented in the study by Arnason et al. (2002)
and the interpretations obtained from the composite mtDNA and nuclear data of Murphy et al. (2001a)
. Thus, the newly proposed clade "Dermosimii" joining together Dermoptera and higher or anthropoid primates (Arnason et al. 2002
) is not supported by SINE-evidence. To tease apart the influence of mtDNA and nuclear DNA information in the composite nuclear and mt data set of Murphy et al. (2001a)
, we reanalyzed this data set by including or excluding mt sequence information. Whereas a separate analysis of nuclear sequence information generated a largely unresolved tree that involved among others the dermopteran representative, the composite data set gives rise to a Dermoptera-Anthropoidea sister group with high support. It therefore seems that the clustering of Dermoptera and higher primates in Murphy et al. (2001a)
is largely driven by mt sequence information.
In addition, our mt data do not agree with the modified phylogenetic position of dermopterans as presented in Murphy et al. (2001b)
. In this article the colugo shows a strong affiliation to Tupaia (up to 100% Bayesian posterior probability). No supporting signals for this affiliation could be found in the mt genome.
Using the presence or absence of retrotransposable elements as evidence, we were essentially able to demonstrate that the nucleotide compositional plasticity of mammalian mt genomes strongly interferes with phylogenetic reconstructions involving primate and primate-related taxa (figs. 1 and 2D ). Sampling only mammalian species with similar nucleotide composition gives rise to a phylogenetic constellation in which colugos and primates share a common ancestor that is exclusive to all other extant eutherian taxa. Given this unsatisfactory incomplete taxonomic sampling, future comparative sequencing projects that reveal the presence of retrotransposable elements shared by dermopterans and primates but not by other living eutherians are urgently required. For the part of nuclear DNA-based evidence, this could provide a final settlement for the debate of primate origins.
![]() |
Supplementary Material |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: Cynocephalus variegatus
colugo
mitochondrial genome
mammalian phylogeny
Address for correspondence and reprints: Jürgen Schmitz, Primate Genetics, German Primate Center, Kellnerweg 4, 37077 Göttingen, Germany. E-mail: jschmitz{at}z.gwdg.de
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Arnason U., J. A. Adegoke, K. Bodin, E. W. Born, Y. B. Esa, A. Gullberg, M. Nilsson, R. V. Short, X. F. Xu, A. Janke, 2002 Mammalian mitogenomic relationships and the root of the eutherian tree Proc. Natl. Acad. Sci 99:8151-8156
Beard K. C., 1990 Gliding behaviour and palaeoecology of the alleged primate family Paromomyidae (Mammalia, Dermoptera) Nature 345:340-341[ISI]
. 1993 Origin and evolution of gliding in early Cenozoic Dermoptera (Mammalia, Primatomorpha) Pp. 6390 in J. G. Fleagle and R. D. E. Macphee, eds. Primates and their relatives in phylogenetic perspective. Plenum Press, N.Y
Felsenstein J., 1995 PHYLIP (phylogeny inference package). Version 3.6 Distributed by the author, Department of Genetics, University of Washington, Seattle
Fleagle J. C., R. D. E. Macphee, 1993 Primates and their relatives in phylogenetic perspective Plenum Press, N.Y
Foster P. G., L. S. Jermiin, D. A. Hickey, 1997 Nucleotide composition bias affects amino acid content in proteins coded by animal mitochondria J. Mol. Evol 44:282-288[ISI][Medline]
Graur D., W. A. Hide, W.-H. Li, 1991 Is the guinea-pig a rodent? Nature 351:649-652[ISI][Medline]
Gregory W. K., 1910 The orders of mammals Bull. Am. Mus. Nat. Hist 27:1-524
Hwu H. R., J. W. Roberts, E. H. Davidson, R. J. Britten, 1986 Insertion and/or deletion of many repeated DNA sequences in human and higher ape evolution Proc. Natl. Acad. Sci. USA 83:3875-3879[Abstract]
Kay R. F., R. W. Thorington, P. Houde, 1990 Eocene plesiadapiform shows affinities with flying lemurs not primates Nature 345:342-344[ISI]
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]
Li W.-H., Z. Gu, H. Wang, A. Nekrutenko, 2001 Evolutionary analyses of the human genome Nature 409:847-849[ISI][Medline]
Lin Y.-H., D. Penny, 2001 Implications for bat evolution from two new complete mitochondrial genomes Mol. Biol. Evol 18:684-688
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 adaptive radiations in two major clades of placental mammals Nature 409:610-614[ISI][Medline]
Martin R. D., 1993 Primate origins: plugging the gaps Nature 363:223-234[ISI][Medline]
Murphy W. J., E. Eizirik, W. E. Johnson, 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
Novacek M. J., 1992 Mammalian phylogeny: shaking the tree Nature 356:121-125[ISI][Medline]
Schmitz J., M. Ohme, H. Zischler, 2000 The complete mitochondrial genome of Tupaia belangeri and the phylogenetic affiliation of Scandentia to other eutherian orders Mol. Biol. Evol 17:1334-1343
. 2001 SINE insertions in cladistic analyses and the phylogenetic affiliations of Tarsius bancanus to other primates Genetics 157:777-784
. 2002 The complete mitochondrial sequence of Tarsius bancanus: evidence for an extensive nucleotide compositional plasticity of primate mitochondrial DNA Mol. Biol. Evol 19:544-553
Shedlock A. M., N. Okada, 2000 SINE insertions: powerful tools for molecular systematics Bioessays 22:148-160[ISI][Medline]
Singer S. S., J. Schmitz, C. Schwiegk, H. Zischler, 2003 Molecular cladistic markers in New World monkey phylogeny (Platyrrhini, Primates) Mol. Phylogenet. Evol. (in press)
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., 2000 PAUP*: phylogenetic analysis using parsimony (*and other methods). Version 4 Sinauer, Sunderland, Mass
Tarrio R., F. Rodriguez-Trelles, F. J. Ayala, 2001 Shared nucleotide composition biases among species and their impact on phylogenetic reconstructions of the Drosophilidae Mol. Biol. Evol 18:1464-1473
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]
Thompson J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, D. G. Higgins, 1997 The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools Nucleic Acids Res 25:4876-4882
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]