Additional support for Afrotheria and Paenungulata, the performance of mitochondrial versus nuclear genes, and the impact of data partitions with heterogeneous base composition

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

An Evaluation of Different Partitioning Strategies for Bayesian Estimation of Species Divergence Times

The explosive growth of molecular sequence data has made it possible to estimate species divergence times under relaxed-clock models using genome-scale data sets with many gene loci. In order to improve both model realism and to best extract information about relative divergence times in the sequence data, it is important to account for the heterogeneity in the evolutionary process across genes or genomic regions. Partitioning is a commonly used approach to achieve those goals. We group sites that have similar evolutionary characteristics into the same partition and those with different characteristics into different partitions, and then use different models or different values of model parameters for different partitions to account for the among-partition heterogeneity. However, how to partition data in practical phylogenetic analysis, and in particular in relaxed-clock dating analysis, is more art than science. Here, we use computer simulation and real data analysis to study the impact of the partition scheme on divergence time estimation. The partition schemes had relatively minor effects on the accuracy of posterior time estimates when the prior assumptions were correct and the clock was not seriously violated, but showed large differences when the clock was seriously violated, when the fossil calibrations were in conflict or incorrect, or when the rate prior was mis-specified. Concatenation produced the widest posterior intervals with the least precision. Use of many partitions increased the precision, as predicted by the infinite-sites theory, but the posterior intervals might fail to include the true ages because of the conflicting fossil calibrations or mis-specified rate priors. We analyzed a data set of 78 plastid genes from 15 plant species with serious clock violation and showed that time estimates differed significantly among partition schemes, irrespective of the rate drift model used. Multiple and precise fossil calibrations reduced the differences among partition schemes and were important to improving the precision of divergence time estimates. While the use of many partitions is an important approach to reducing the uncertainty in posterior time estimates, we do not recommend its general use for the present, given the limitations of current models of rate drift for partitioned data and the challenges of interpreting the fossil evidence to construct accurate and informative calibrations.

Life on the rocks: Multilocus phylogeography of rock hyrax (Procavia capensis) from southern Africa

Understanding the role of geography and climatic cycles in determining patterns of biodiversity is important in comparative and evolutionary biology and conservation. We studied the phylogeographic pattern and historical demography of a rock-dwelling small mammal species from southern Africa, the rock hyrax Procavia capensis capensis. Using a multilocus coalescent approach, we assessed the influence of strong habitat dependence and fluctuating regional climates on genetic diversity. We sequenced a mitochondrial gene (cytochrome b) and two nuclear introns (AP5, PRKC1) supplemented with microsatellite genotyping, in order to assess evolutionary processes over multiple temporal scales. In addition, distribution modelling was used to investigate the current and predicted distribution of the species under different climatic scenarios. Collectively, the data reveal a complex history of isolation followed by secondary contact shaping the current intraspecific diversity. The cyt b sequences confirmed the presence of two previously proposed geographically and genetically distinct lineages distributed across the southern African Great Escarpment and north-western mountain ranges. Molecular dating suggests Miocene divergence of the lineages, yet there are no discernible extrinsic barriers to gene flow. The nuclear markers reveal incomplete lineage sorting or ongoing mixing of the two lineages. Although the microsatellite data lend some support to the presence of two subpopulations, there is weak structuring within and between lineages. These data indicate the presence of gene flow from the northern into the southern parts of the southern African sub-region likely following the secondary contact. The distribution modelling predictably reveal the species' preference for rocky areas, with stable refugia through time in the northern mountain ranges, the Great Escarpment, as well as restricted areas of the Northern Cape Province and the Cape Fold Mountains of South Africa. Different microclimatic variables appear to determine the distributional range of the species. Despite strong habitat preference, the micro-habitat offered by rocky crevices and unique life history traits likely promoted the adaptability of P. capensis, resulting in the widespread distribution and persistence of the species over a long evolutionary period. Spatio-temporal comparison of the evolutionary histories of other co-distributed species across the rocky landscapes of southern Africa will improve our understanding of the regional patterns of biodiversity and local endemism.

Automatic selection of partitioning schemes for phylogenetic analyses using iterative k-means clustering of site rates

BACKGROUND

Model selection is a vital part of most phylogenetic analyses, and accounting for the heterogeneity in evolutionary patterns across sites is particularly important. Mixture models and partitioning are commonly used to account for this variation, and partitioning is the most popular approach. Most current partitioning methods require some a priori partitioning scheme to be defined, typically guided by known structural features of the sequences, such as gene boundaries or codon positions. Recent evidence suggests that these a priori boundaries often fail to adequately account for variation in rates and patterns of evolution among sites. Furthermore, new phylogenomic datasets such as those assembled from ultra-conserved elements lack obvious structural features on which to define a priori partitioning schemes. The upshot is that, for many phylogenetic datasets, partitioned models of molecular evolution may be inadequate, thus limiting the accuracy of downstream phylogenetic analyses.

RESULTS

We present a new algorithm that automatically selects a partitioning scheme via the iterative division of the alignment into subsets of similar sites based on their rates of evolution. We compare this method to existing approaches using a wide range of empirical datasets, and show that it consistently leads to large increases in the fit of partitioned models of molecular evolution when measured using AICc and BIC scores. In doing so, we demonstrate that some related approaches to solving this problem may have been associated with a small but important bias.

CONCLUSIONS

Our method provides an alternative to traditional approaches to partitioning, such as dividing alignments by gene and codon position. Because our method is data-driven, it can be used to estimate partitioned models for all types of alignments, including those that are not amenable to traditional approaches to partitioning.

Summary of Laurasiatheria (mammalia) phylogeny

Laurasiatheria is one of the richest and most diverse superorders of placental mammals. Because this group had a rapid evolutionary radiation, the phylogenetic relationships among the six orders of Laurasiatheria remain a subject of heated debate and several issues related to its phylogeny remain open. Reconstructing the true phylogenetic relationships of Laurasiatheria is a significant case study in evolutionary biology due to the diversity of this suborder and such research will have significant implications for biodiversity conservation. We review the higher-level (inter-ordinal) phylogenies of Laurasiatheria based on previous cytogenetic, morphological and molecular data, and discuss the controversies of its phylogenetic relationship. This review aims to outline future researches on Laurasiatheria phylogeny and adaptive evolution.

Comparative Anatomy of the Bony Labyrinth (Inner Ear) of Placental Mammals

BACKGROUND

Variation is a naturally occurring phenomenon that is observable at all levels of morphology, from anatomical variations of DNA molecules to gross variations between whole organisms. The structure of the otic region is no exception. The present paper documents the broad morphological diversity exhibited by the inner ear region of placental mammals using digital endocasts constructed from high-resolution X-ray computed tomography (CT). Descriptions cover the major placental clades, and linear, angular, and volumetric dimensions are reported.

PRINCIPAL FINDINGS

The size of the labyrinth is correlated to the overall body mass of individuals, such that large bodied mammals have absolutely larger labyrinths. The ratio between the average arc radius of curvature of the three semicircular canals and body mass of aquatic species is substantially lower than the ratios of related terrestrial taxa, and the volume percentage of the vestibular apparatus of aquatic mammals tends to be less than that calculated for terrestrial species. Aspects of the bony labyrinth are phylogenetically informative, including vestibular reduction in Cetacea, a tall cochlear spiral in caviomorph rodents, a low position of the plane of the lateral semicircular canal compared to the posterior canal in Cetacea and Carnivora, and a low cochlear aspect ratio in Primatomorpha.

SIGNIFICANCE

The morphological descriptions that are presented add a broad baseline of anatomy of the inner ear across many placental mammal clades, for many of which the structure of the bony labyrinth is largely unknown. The data included here complement the growing body of literature on the physiological and phylogenetic significance of bony labyrinth structures in mammals, and they serve as a source of data for future studies on the evolution and function of the vertebrate ear.

Phylogeny and taxonomy of the round-eared sengis or elephant-shrews, genus Macroscelides (Mammalia, Afrotheria, Macroscelidea)

The round-eared sengis or elephant-shrews (genus Macroscelides) exhibit striking pelage variation throughout their ranges. Over ten taxonomic names have been proposed to describe this variation, but currently only two taxa are recognized (M. proboscideus proboscideus and M. p. flavicaudatus). Here, we review the taxonomic history of Macroscelides, and we use data on the geographic distribution, morphology, and mitochondrial DNA sequence to evaluate the current taxonomy. Our data support only two taxa that correspond to the currently recognized subspecies M. p. proboscideus and M. p. flavicaudatus. Mitochondrial haplotypes of these two taxa are reciprocally monophyletic with over 13% uncorrected sequence divergence between them. PCA analysis of 14 morphological characters (mostly cranial) grouped the two taxa into non-overlapping clusters, and body mass alone is a relatively reliable distinguishing character throughout much of Macroscelides range. Although fieldworkers were unable to find sympatric populations, the two taxa were found within 50 km of each other, and genetic analysis showed no evidence of gene flow. Based upon corroborating genetic data, morphological data, near sympatry with no evidence of gene flow, and differences in habitat use, we elevate these two forms to full species.

Sources of signal in 62 protein-coding nuclear genes for higher-level phylogenetics of arthropods

BACKGROUND

This study aims to investigate the strength of various sources of phylogenetic information that led to recent seemingly robust conclusions about higher-level arthropod phylogeny and to assess the role of excluding or downweighting synonymous change for arriving at those conclusions.

METHODOLOGY/PRINCIPAL FINDINGS

The current study analyzes DNA sequences from 68 gene segments of 62 distinct protein-coding nuclear genes for 80 species. Gene segments analyzed individually support numerous nodes recovered in combined-gene analyses, but few of the higher-level nodes of greatest current interest. However, neither is there support for conflicting alternatives to these higher-level nodes. Gene segments with higher rates of nonsynonymous change tend to be more informative overall, but those with lower rates tend to provide stronger support for deeper nodes. Higher-level nodes with bootstrap values in the 80% - 99% range for the complete data matrix are markedly more sensitive to substantial drops in their bootstrap percentages after character subsampling than those with 100% bootstrap, suggesting that these nodes are likely not to have been strongly supported with many fewer data than in the full matrix. Data set partitioning of total data by (mostly) synonymous and (mostly) nonsynonymous change improves overall node support, but the result remains much inferior to analysis of (unpartitioned) nonsynonymous change alone. Clusters of genes with similar nonsynonymous rate properties (e.g., faster vs. slower) show some distinct patterns of node support but few conflicts. Synonymous change is shown to contribute little, if any, phylogenetic signal to the support of higher-level nodes, but it does contribute nonphylogenetic signal, probably through its underlying heterogeneous nucleotide composition. Analysis of seemingly conservative indels does not prove useful.

CONCLUSIONS

Generating a robust molecular higher-level phylogeny of Arthropoda is currently possible with large amounts of data and an exclusive reliance on nonsynonymous change.

Phylogenomic analysis resolves the formerly intractable adaptive diversification of the endemic clade of east Asian Cyprinidae (Cypriniformes)

Despite their great diversity and biological importance, evolutionary relationships among the endemic clade of East Asian Cyprinidae remain ambiguous. Understanding the phylogenetic history of this group involves many challenges. For instance, ecomorphological convergence may confound morphology-based phylogenetic inferences, and previous molecular phylogenetic studies based on single genes have often yielded contradictory and poorly supported trees. We assembled a comprehensive data matrix of 100 nuclear gene segments (∼ 71132 base pairs) for representative species of the endemic East Asian cyprinid fauna and recovered a robust phylogeny from this genome-wide signal supported by multiple analytical methods, including maximum parsimony, maximum likelihood and Bayesian inference. Relaxed molecular clock analyses indicated species radiations of this clade concentrated at approximately 1.9–7.6 MYA. We provide evidence that the bursts of diversification in this fauna are directly linked to major paleoenvironmental events associated with monsoon evolution occurring from late Miocene to Pliocene. Ancestral state reconstruction reveals convergent morphological characters are hypothesized to be independent products of similar selective pressures in ecosystems. Our study is the first comprehensive phylogenetic study of the enigmatic East-Asian cyprinids. The explicit molecular phylogeny provides a valuable framework for future research in genome evolution, adaptation and speciation of cyprinids.

Base composition, selection, and phylogenetic significance of indels in the recombination activating gene-1 in vertebrates

BACKGROUND

The Recombination Activating Proteins, RAG1 and RAG2, play a crucial role in the immune response in vertebrates. Among the nuclear markers currently used for phylogenetic purposes, Rag1 has especially enjoyed enormous popularity, since it successfully contributed to elucidating the relationships among and within a large variety of vertebrate lineages. We here report on a comparative investigation of the genetic variation, base composition, presence of indels, and selection in Rag1 in different vertebrate lineages (Actinopterygii, Amphibia, Aves, Chondrichthyes, Crocodylia, Lepidosauria, Mammalia, and Testudines) through the analysis of 582 sequences obtained from Genbank. We also analyze possible differences between distinct parts of the gene with different type of protein functions.

RESULTS

In the vertebrate lineages studied, Rag1 is over 3 kb long. We observed a high level of heterogeneity in base composition at the 3(rd )codon position in some of the studied vertebrate lineages and in some specific taxa. This result is also paralleled by taxonomic differences in the GC content at the same codon position. Moreover, positive selection occurs at some sites in Aves, Lepidosauria and Testudines. Indels, which are often used as phylogenetic characters, are more informative across vertebrates in the 5' than in the 3'-end of the gene. When the entire gene is considered, the use of indels as phylogenetic character only recovers one major vertebrate clade, the Actinopterygii. However, in numerous cases insertions or deletions are specific to a monophyletic group.

CONCLUSIONS

Rag1 is a phylogenetic marker of undoubted quality. Our study points to the need of carrying out a preliminary investigation on the base composition and the possible existence of sites under selection of this gene within the groups studied to avoid misleading resolution. The gene shows highly heterogeneous base composition, which affects some taxa in particular and contains sites under positive selection in some vertebrate lineages in the 5'-end. The first part of the gene (5'-end) is more variable than the second (3'-end), and less affected by a heterogeneous base composition. However, in some vertebrate lineages the 5'-end of the gene is not yet widely used for phylogenetic studies.

The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals

BACKGROUND

During the last ten years, major advances have been made in characterizing and understanding the evolution of mitochondrial DNA, the most popular marker of molecular biodiversity. Several important results were recently reported using mammals as model organisms, including (i) the absence of relationship between mitochondrial DNA diversity and life-history or ecological variables, (ii) the absence of prominent adaptive selection, contrary to what was found in invertebrates, and (iii) the unexpectedly large variation in neutral substitution rate among lineages, revealing a possible link with species maximal longevity. We propose to challenge these results thanks to the bird/mammal comparison. Direct estimates of population size are available in birds, and this group presents striking life-history trait differences with mammals (higher mass-specific metabolic rate and longevity). These properties make birds the ideal model to directly test for population size effects, and to discriminate between competing hypotheses about the causes of substitution rate variation.

RESULTS

A phylogenetic analysis of cytochrome b third-codon position confirms that the mitochondrial DNA mutation rate is quite variable in birds, passerines being the fastest evolving order. On average, mitochondrial DNA evolves slower in birds than in mammals of similar body size. This result is in agreement with the longevity hypothesis, and contradicts the hypothesis of a metabolic rate-dependent mutation rate. Birds show no footprint of adaptive selection on cytochrome b evolutionary patterns, but no link between direct estimates of population size and cytochrome b diversity. The mutation rate is the best predictor we have of within-species mitochondrial diversity in birds. It partly explains the differences in mitochondrial DNA diversity patterns observed between mammals and birds, previously interpreted as reflecting Hill-Robertson interferences with the W chromosome.

CONCLUSION

Mitochondrial DNA diversity patterns in birds are strongly influenced by the wide, unexpected variation of mutation rate across species. From a fundamental point of view, these results are strongly consistent with a relationship between species maximal longevity and mitochondrial mutation rate, in agreement with the mitochondrial theory of ageing. Form an applied point of view, this study reinforces and extends the message of caution previously expressed for mammals: mitochondrial data tell nothing about species population sizes, and strongly depart the molecular clock assumption.

A new perspective on phylogeny and evolution of tetraodontiform fishes (Pisces: Acanthopterygii) based on whole mitochondrial genome sequences: basal ecological diversification?

Background

The order Tetraodontiformes consists of approximately 429 species of fishes in nine families. Members of the order exhibit striking morphological diversity and radiated into various habitats such as freshwater, brackish and coastal waters, open seas, and deep waters along continental shelves and slopes. Despite extensive studies based on both morphology and molecules, there has been no clear resolution except for monophyly of each family and sister-group relationships of Diodontidae + Tetraodontidae and Balistidae + Monacanthidae. To address phylogenetic questions of tetraodontiform fishes, we used whole mitochondrial genome (mitogenome) sequences from 27 selected species (data for 11 species were newly determined during this study) that fully represent all families and subfamilies of Tetraodontiformes (except for Hollardinae of the Triacanthodidae). Partitioned maximum likelihood (ML) and Bayesian analyses were performed on two data sets comprising concatenated nucleotide sequences from 13 protein-coding genes (all positions included; third codon positions converted into purine [R] and pyrimidine [Y]), 22 transfer RNA and two ribosomal RNA genes (total positions = 15,084).

Results

The resultant tree topologies from the two data sets were congruent, with many internal branches showing high support values. The mitogenomic data strongly supported monophyly of all families and subfamilies (except the Tetraodontinae) and sister-group relationships of Balistidae + Monacanthidae and Tetraodontidae + Diodontidae, confirming the results of previous studies. However, we also found two unexpected basal splits into Tetraodontoidei (Triacanthidae + Balistidae + Monacanthidae + Tetraodontidae + Diodontidae + Molidae) and Triacanthodoidei (Ostraciidae + Triodontidae + Triacanthodidae).

Conclusion

This basal split into the two clades has never been reported and challenges previously proposed hypotheses based on both morphology and nuclear gene sequences. It is likely that the basal split had involved ecological diversification, because most members of Tetraodontoidei exclusively occur in shallow waters (freshwater, brackish and coastal waters, and open seas), while those of Triacanthodoidei occur mainly in relatively deep waters along continental shelves and slopes except for more derived ostraciids. This suggests that the basal split between the two clades led to subsequent radiation into the two different habitats.

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei: Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A>T, C>G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA.

Primate Origins: Implications of a Cretaceous Ancestry

It has long been accepted that the adaptive radiation of modern placental mammals, like that of modern birds, did not begin until after the Cretaceous/Tertiary (K/T) boundary 65 million years (Ma) ago, following the extinction of the dinosaurs. The first undoubted fossil relatives of modern primates appear in the record 55 Ma ago. However, in agreement with evidence from molecular phylogenies calibrated with dates from denser parts of the fossil record, a statistical analysis of the primate record allowing for major gaps now indicates a Cretaceous origin of euprimates 80-90 Ma ago. If this interpretation is correct, primates overlapped with dinosaurs by some 20 Ma prior to the K/T boundary, and the initial radiation of primates was probably truncated as part of the major extinction event that occurred at the end of the Cretaceous. Following a review of evidence for an early origin of primates, implications of this are discussed with respect to the likely ancestral condition for primates, including a southern continental area of origin and moderately large body size. The known early Tertiary primates are re-interpreted as northern continental offshoots of a 'second wave' of primate evolution.

Indel evolution of mammalian introns and the utility of non-coding nuclear markers in eutherian phylogenetics

Nuclear DNA intron sequences are increasingly used to investigate evolutionary relationships among closely related organisms. The phylogenetic usefulness of intron sequences at higher taxonomic levels has, however, not been firmly established and very few studies have used these markers to address evolutionary questions above the family level. In addition, the mechanisms driving intron evolution are not well understood. We compared DNA sequence data derived from three presumably independently segregating introns (THY, PRKC I and MGF) across 158 mammalian species. All currently recognized extant eutherian mammalian orders were included with the exception of Cingulata, Dermoptera and Scandentia. The total aligned length of the data was 6366 base pairs (bp); after the exclusion of autapomorphic insertions, 1511 bp were analyzed. In many instances the Bayesian and parsimony analyses were complementary and gave significant posterior probability and bootstrap support (>80) for the monophyly of Afrotheria, Euarchontoglires, Laurasiatheria and Boreoeutheria. Apart from finding congruent support when using these methods, the intron data also provided several indels longer than 3 bp that support, among others, the monophyly of Afrotheria, Paenungulata, Ferae and Boreoeutheria. A quantitative analysis of insertions and deletions suggested that there was a 75% bias towards deletions. The average insertion size in the mammalian data set was 16.49 bp +/- 57.70 while the average deletion was much smaller (4.47 bp +/- 14.17). The tendency towards large insertions and small deletions is highlighted by the observation that out of a total of 17 indels larger than 100 bp, 15 were insertions. The majority of indels (>60% of all events) were 1 or 2 bp changes. Although the average overall indel substitution rate of 0.00559 per site is comparable to that previously reported for rodents and primates, individual analyses among different evolutionary lineages provide evidence for differences in the formation rate of indels among the different mammalian groups.

Chromosome painting in the manatee supports Afrotheria and Paenungulata

BACKGROUND

Sirenia (manatees, dugongs and Stellar's sea cow) have no evolutionary relationship with other marine mammals, despite similarities in adaptations and body shape. Recent phylogenomic results place Sirenia in Afrotheria and with elephants and rock hyraxes in Paenungulata. Sirenia and Hyracoidea are the two afrotherian orders as yet unstudied by comparative molecular cytogenetics. Here we report on the chromosome painting of the Florida manatee.

RESULTS

The human autosomal and X chromosome paints delimited a total of 44 homologous segments in the manatee genome. The synteny of nine of the 22 human autosomal chromosomes (4, 5, 6, 9, 11, 14, 17, 18 and 20) and the X chromosome were found intact in the manatee. The syntenies of other human chromosomes were disrupted in the manatee genome into two to five segments. The hybridization pattern revealed that 20 (15 unique) associations of human chromosome segments are found in the manatee genome: 1/15, 1/19, 2/3 (twice), 3/7 (twice), 3/13, 3/21, 5/21, 7/16, 8/22, 10/12 (twice), 11/20, 12/22 (three times), 14/15, 16/19 and 18/19.

CONCLUSION

There are five derived chromosome traits that strongly link elephants with manatees in Tethytheria and give implicit support to Paenungulata: the associations 2/3, 3/13, 8/22, 18/19 and the loss of the ancestral eutherian 4/8 association. It would be useful to test these conclusions with chromosome painting in hyraxes. The manatee chromosome painting data confirm that the associations 1/19 and 5/21 phylogenetically link afrotherian species and show that Afrotheria is a natural clade. The association 10/12/22 is also ubiquitous in Afrotheria (clade I), present in Laurasiatheria (clade IV), only partially present in Xenarthra (10/12, clade II) and absent in Euarchontoglires (clade III). If Afrotheria is basal to eutherians, this association could be part of the ancestral eutherian karyotype. If afrotherians are not at the root of the eutherian tree, then the 10/12/22 association could be one of a suite of derived associations linking afrotherian taxa.

The evolutionary radiation of Arvicolinae rodents (voles and lemmings): relative contribution of nuclear and mitochondrial DNA phylogenies

BACKGROUND

Mitochondrial and nuclear genes have generally been employed for different purposes in molecular systematics, the former to resolve relationships within recently evolved groups and the latter to investigate phylogenies at a deeper level. In the case of rapid and recent evolutionary radiations, mitochondrial genes like cytochrome b (CYB) are often inefficient for resolving phylogenetic relationships. One of the best examples is illustrated by Arvicolinae rodents (Rodentia; Muridae), the most impressive mammalian radiation of the Northern Hemisphere which produced voles, lemmings and muskrats. Here, we compare the relative contribution of a nuclear marker--the exon 10 of the growth hormone receptor (GHR) gene--to the one of the mitochondrial CYB for inferring phylogenetic relationships among the major lineages of arvicoline rodents.

RESULTS

The analysis of GHR sequences improves the overall resolution of the Arvicolinae phylogeny. Our results show that the Caucasian long-clawed vole (Prometheomys schaposnikowi) is one of the basalmost arvicolines, and confirm that true lemmings (Lemmus) and collared lemmings (Dicrostonyx) are not closely related as suggested by morphology. Red-backed voles (Myodini) are found as the sister-group of a clade encompassing water vole (Arvicola), snow vole (Chionomys), and meadow voles (Microtus and allies). Within the latter, no support is recovered for the generic recognition of Blanfordimys, Lasiopodomys, Neodon, and Phaiomys as suggested by morphology. Comparisons of parameter estimates for branch lengths, base composition, among sites rate heterogeneity, and GTR relative substitution rates indicate that CYB sequences consistently exhibit more heterogeneity among codon positions than GHR. By analyzing the contribution of each codon position to node resolution, we show that the apparent higher efficiency of GHR is due to their third positions. Although we focus on speciation events spanning the last 10 million years (Myr), CYB sequences display highly saturated codon positions contrary to the nuclear exon. Lastly, variable length bootstrap predicts a significant increase in resolution of arvicoline phylogeny through the sequencing of nuclear data in an order of magnitude three to five times greater than the size of GHR exon 10.

CONCLUSION

Our survey provides a first resolved gene tree for Arvicolinae. The comparison of CYB and GHR phylogenetic efficiency supports recent assertions that nuclear genes are useful for resolving relationships of recently evolved animals. The superiority of nuclear exons may reside both in (i) less heterogeneity among sites, and (ii) the presence of highly informative sites in third codon positions, that evolve rapidly enough to accumulate synapomorphies, but slow enough to avoid substitutional saturation.

The pattern and timing of diversification of Philippine endemic rodents: evidence from mitochondrial and nuclear gene sequences

The 22 genera and 64 species of rodents (Muridae: Murinae) distributed in the Philippine Islands provide a unique opportunity to study patterns and processes of diversification in island systems. Over 90% of these rodent species are endemic to the archipelago, but the relative importance of dispersal from the mainland, dispersal within the archipelago, and in situ differentiation as explanations of this diversity remains unclear, as no phylogenetic hypothesis for these species and relevant mainland forms is currently available. Here we report the results of phylogenetic analyses of the endemic Philippine murines and a wide sampling of murine diversity from outside the archipelago, based on the mitochondrial cytochrome b gene and the nuclear-encoded IRBP exon 1. Analysis of our combined gene data set consistently identified five clades comprising endemic Philippine genera, suggesting multiple invasions of the archipelago. Molecular dating analyses using parametric and semiparametric methods suggest that colonization occurred in at least two stages, one ca. 15 Mya, and another 8 to 12 million years later, consistent with the previous recognition of "Old" and "New" endemic rodent faunas. Ancestral area analysis suggests that the Old Endemics invaded landmasses that are now part of the island of Luzon, whereas the three New Endemic clades may have colonized through either Mindanao, Luzon, or both. Further, our results suggest that most of the diversification of Philippine murines took place within the archipelago. Despite heterogeneity between nuclear and mitochondrial genes in most model parameters, combined analysis of the two data sets using both parsimony and likelihood increased phylogenetic resolution; however, the effect of data combination on support for resolved nodes was method dependent. In contrast, our results suggest that combination of mitochondrial and nuclear data to estimate relatively ancient divergence times can severely compromise those estimates, even when specific methods that account for rate heterogeneity among genes are employed. [Biogeography; divergence date estimation; mitochondrial DNA; molecular systematics; Murinae; nuclear exon; Philippines; phylogeny.].

Preliminary genetic characterization of two lineages of black rats (Rattus rattus sensu lato) in Japan, with evidence for introgression at several localities

We conducted a pilot survey of genetic diversity among 37 karyotyped individuals of the black rat Rattus rattus (sensu lato) from six localities on the Japanese Islands, using complete gene sequences of mitochondrial cytochrome b (cyt b) and nuclear interphotoreceptor retinoid binding protein (IRBP). Our sampling included two previously documented karyotypic groups: 'Oceanian' with 2n = 38 and 'Asian' with 2n = 42. Cyt b sequences for most individuals clustered according to their karyotypic groups, with an average between-group divergence of 3.8%. One exception was that individuals from Kagoshima (Kyushu Island) showed 'Asian' karyotypes combined with a cyt b haplotype that differed by a single nucleotide substitution from the haplotype of the 'Oceanian' karyotypic group. Six IRBP haplotypes were identified. They belonged to three distinct IRBP lineages (I-III), with an average inter-lineage divergence of 1%. Among homozygous individuals, these lineages showed good association with the karyotypic groups: IRBP lineage I occurred only with 'Oceanian' karyotypes, while IRBP lineages II and III both occurred with 'Asian' karyotypes. Individuals from Kagoshima all possessed IRBP of 'Asian' lineages, despite the presence of an 'Oceanian' mitochondrial type. The Chichijima population (Ogasawara Islands) featured exclusively 'Asian' karyotypes and cyt b sequences, but various combinations of all three IRBP lineages. The Kagoshima and Chichijima populations thus provide strong evidence of viable hybridization and genetic introgression between the two karyotypic groups, but with variable genetic outcomes. Our results demonstrate the potential of combined analysis of karyotypes and mitochondrial and nuclear gene sequences to elucidate the complex dispersal and population history of the black rat.

On the phylogenetic position of a rare Iberian endemic mammal, the Pyrenean desman (Galemys pyrenaicus)

The nucleotide sequences of the complete mitochondrial genome and nine partial nuclear genes of the Pyrenean desman (Galemys pyrenaicus) were determined in order to establish the relative phylogenetic position of this species at different taxonomic levels within the placental tree. Phylogenetic relationships of desman within the family Talpidae were inferred based on complete mitochondrial cytochrome b gene nucleotide sequence data. The Pyrenean desman was unambiguously recovered as sister group of the Russian desman (Desmana moschata) confirming the monophyly of the subfamily Desmaninae. However, phylogenetic relationships among major lineages within the Talpidae could not be confidently resolved. Phylogenetic analyses based on mitochondrial (at the amino acid level) and nuclear (at the nucleotide level) sequence data sets confidently placed desman within the Eulipotyphla (that also included moles, shrews, and hedgehogs), and partially recovered placental interordinal relationships. The monophyly of Laurasiatheria (including Eulipotyphla, Chiroptera, Carnivora, Pholidota, Perissodactyla, and Cetartiodactyla) was strongly supported. Mitochondrial amino acid sequences of Pholidota (pangolins) were found to bias phylogenetic inferences due to long-branch attraction effects. A Bayesian inference based on a combined mitochondrial and nuclear data set without Pholidota arrived at an almost fully resolved tree that supported the basal position of Eulipotyphla within Laurasiatheria.

Combined Mitochondrial and Nuclear DNA Sequences Resolve the Interrelations of the Major Australasian Marsupial Radiations

Australasian marsupials include three major radiations, the insectivorous/carnivorous Dasyuromorphia, the omnivorous bandicoots (Peramelemorphia), and the largely herbivorous diprotodontians. Morphologists have generally considered the bandicoots and diprotodontians to be closely related, most prominently because they are both syndactylous (with the 2nd and 3rd pedal digits being fused). Molecular studies have been unable to confirm or reject this Syndactyla hypothesis. Here we present new mitochondrial (mt) genomes from a spiny bandicoot (Echymipera rufescens) and two dasyurids, a fat-tailed dunnart (Sminthopsis crassicaudata) and a northern quoll (Dasyurus hallucatus). By comparing trees derived from pairwise base-frequency differences between taxa with standard (absolute, uncorrected) distance trees, we infer that composition bias among mt protein-coding and RNA sequences is sufficient to mislead tree reconstruction. This can explain incongruence between trees obtained from mt and nuclear data sets. However, after excluding major sources of compositional heterogeneity, both the "reduced-bias" mt and nuclear data sets clearly favor a bandicoot plus dasyuromorphian association, as well as a grouping of kangaroos and possums (Phalangeriformes) among diprotodontians. Notably, alternatives to these groupings could only be confidently rejected by combining the mt and nuclear data. Elsewhere on the tree, Dromiciops appears to be sister to the monophyletic Australasian marsupials, whereas the placement of the marsupial mole (Notoryctes) remains problematic. More generally, we contend that it is desirable to combine mt genome and nuclear sequences for inferring vertebrate phylogeny, but as separately modeled process partitions. This strategy depends on detecting and excluding (or accounting for) major sources of non-historical signal, such as from compositional non-stationarity. [Base composition; combined data; marsupial; mitochondrial genome; phylogeny.].

A retroposon analysis of Afrotherian phylogeny

Recent comprehensive studies of DNA sequences support the monophyly of Afrotheria, comprising elephants, sirenians (dugongs and manatees), hyraxes, tenrecs, golden moles, aardvarks, and elephant shrews, as well as that of Paenungulata, comprising elephants, sirenians, and hyraxes. However, phylogenetic relationships among paenungulates, as well as among nonpaenungulates, have remained ambiguous. Here we applied an extensive retroposon analysis to these problems to support the monophyly of aardvarks, tenrecs, and golden moles, with elephant shrews as their sister group. Regarding phylogenetic relationships in Paenungulata, we could characterize only one informative locus, although we could isolate many insertions specific to each of three lineages, namely, Proboscidea, Sirenia, and Hyracoidea. These data prompted us to reexamine phylogenetic relationships among Paenungulata using 19 nuclear gene sequences resulting in three different analyses, namely, short interspersed element (SINE) insertions, nuclear sequence analyses, and morphological cladistics, supporting different respective phylogenies. We concluded that these three lineages diverged very rapidly in a very short evolutionary period, with the consequence that ancestral polymorphism present in the last common ancestor of Paenungulata results in such incongruence. Our results suggest the rapid fixation of many large-scale morphological synapomorphies for Tethytheria; implications of this in relation to the morphological evolution in Paenungulata are discussed.

Phylogenomics and the reconstruction of the tree of life

As more complete genomes are sequenced, phylogenetic analysis is entering a new era - that of phylogenomics. One branch of this expanding field aims to reconstruct the evolutionary history of organisms on the basis of the analysis of their genomes. Recent studies have demonstrated the power of this approach, which has the potential to provide answers to several fundamental evolutionary questions. However, challenges for the future have also been revealed. The very nature of the evolutionary history of organisms and the limitations of current phylogenetic reconstruction methods mean that part of the tree of life might prove difficult, if not impossible, to resolve with confidence.

Empirical problems of the hierarchical likelihood ratio test for model selection

Advocates of maximum likelihood (ML) approaches to phylogenetics commonly cite as one of their primary advantages the use of objective statistical criteria for model selection. Currently, a particular implementation of the likelihood ratio test (LRT) is the most commonly used model-selection criterion in phylogenetics. This approach requires the choice of a starting point and a parameter addition (or removal) sequence that can affect all ML inferences (i.e., topology, model, and all evolutionary parameters). Here, several alternative starting points and parameter sequences are tested in empirical data sets to assess their influence on model selection and optimal topology. In the studied data sets, varying model-selection protocols leads to selection of different models that, in some cases, lead to different ML trees. Given the sensitivity of the LRT, some possible solutions to model selection (within the hypothesis testing approach) are outlined, and alternative model-selection criteria are discussed. Some of the suggested alternatives seem to lack these problems, although their behavior and adequacy for phylogenetics needs to be further explored.

Phylogenetic relationships and divergence times of the genusTokudaiawithin Murinae (Muridae; Rodentia) inferred from the nucleotide sequences encoding the Cytbgene, RAG 1, and IRBP

The phylogenetic relationships of six genera of Murinae (Apodemus, Diplothrix, Micromys, Mus, Rattus, and Tokudaia) were examined using the nucleotide sequences for the mitochondrial cytochrome b (Cytb), as well as the nuclear recombination activating gene 1 (RAG1) and interphotoreceptor retinoid-binding protein (IRBP), with special emphasis on the position of the genus Tokudaia, which is endemic to the Ryukyu Islands. Compared with Cytb at all codon positions, the first and second codon positions of Cytb, RAG1 (1002 base pairs (bp)), and IRBP (1586 bp) sequences were less prone to saturation. Close affinity between the genera Tokudaia and Apodemus was observed in the analyses using the IRBP (1586 bp) and combined nuclear (2588 bp; RAG1 + IRBP) sequences. The divergence time for the Tokudaia–Apodemus clade was estimated at approximately 6.5–8.0 Ma, which is more recent than previously reported, thereby indicating the recent colonization of the Ryukyu Islands by the genus Tokudaia. The other relationships among the main genera were highly ambiguous, owing either to saturation or insufficient phylogenetic information. The radiation of the main genera within a relatively short period of evolutionary time may explain the unresolved topologies, although molecular sources that are less subject to saturation are required to resolve the outstanding issues.

Phylogeny of Neotropical oryzomyine rodents (Muridae: Sigmodontinae) based on the nuclear IRBP exon

Sigmodontine rodents are the most diverse family-level mammalian clade in the Neotropical region, with about 70 genera and 320 recognized species. Partial sequences (1266 bp) from the first exon of the nuclear gene encoding the Interphotoreceptor Retinoid Binding Protein (IRBP) were used to infer the phylogenetic relationships among 44 species representing all 16 currently recognized genera of the largest sigmodontine tribe, the Oryzomyini. Monophyly of the tribe was assessed relative to 15 non-oryzomyine sigmodontine taxa representing all major sigmodontine lineages. Twelve taxa from seven muroid subfamilies were used as outgroups. The resulting matrix included 71 taxa and 386 parsimony-informative characters. Phylogenetic analysis of this matrix resulted in 16 equally parsimonious cladograms, which contained the following well-supported groups: (i). a monophyletic Oryzomyini, (ii). a clade containing all oryzomyines except Scolomys and Zygodontomys, (iii). a clade containing Oecomys, Handleyomys, and several species of forest-dwelling Oryzomys, and (iv). a clade containing the remaining oryzomyine taxa. The last clade is composed of two large subclades, each with lower nodal support, containing the following taxa: (i). Microryzomys, Oligoryzomys, Neacomys, and Oryzomys balneator; (ii). Holochilus, Lundomys, Pseudoryzomys, Nectomys, Amphinectomys, Sigmodontomys, and several species of open-vegetation or semiaquatic Oryzomys. Regarding relationships among non-oryzomyine taxa, sigmodontines, neotomines, and tylomyines do not form a monophyletic group; a clade containing Rheomys and Sigmodon is basal relative to all other sigmodontines; and the remaining sigmodontines are grouped in three clades: the first containing Thomasomyini, Akodontini, and Reithrodon; the second containing Abrothrichini, and Phyllotini, plus Wiedomys, Juliomys, Irenomys, and Delomys; and the third containing the oryzomyines. No conflict is observed between IRBP results and previous robust hypotheses from mitochondrial data, while a single case of incongruence is present between the IRBP topology and robust hypothesis from morphological studies.

Genome-scale phylogeny and the detection of systematic biases

Phylogenetic inference from sequences can be misled by both sampling (stochastic) error and systematic error (nonhistorical signals where reality differs from our simplified models). A recent study of eight yeast species using 106 concatenated genes from complete genomes showed that even small internal edges of a tree received 100% bootstrap support. This effective negation of stochastic error from large data sets is important, but longer sequences exacerbate the potential for biases (systematic error) to be positively misleading. Indeed, when we analyzed the same data set using minimum evolution optimality criteria, an alternative tree received 100% bootstrap support. We identified a compositional bias as responsible for this inconsistency and showed that it is reduced effectively by coding the nucleotides as purines and pyrimidines (RY-coding), reinforcing the original tree. Thus, a comprehensive exploration of potential systematic biases is still required, even though genome-scale data sets greatly reduce sampling error.

Phylogeny of muroid rodents: relationships within and among major lineages as determined by IRBP gene sequences

The rodent family Muridae is the single most diverse family of mammals with over 1300 recognized species. We used DNA sequences from the first exon ( approximately 1200bp) of the IRBP gene to infer phylogenetic relationships within and among the major lineages of muroid rodents. We included sequences from every recognized muroid subfamily except Platacanthomyinae and from all genera within the endemic Malagasy subfamily Nesomyinae, all recognized tribes of Sigmodontinae, and a broad sample of genera in Murinae. Phylogenetic analysis of the IRBP data suggest that muroid rodents can be sorted into five major lineages: (1) a basal clade containing the fossorial rodents in the subfamilies Spalacinae, Myospalacinae, and Rhizomyinae, (2) a clade of African and Malagasy genera comprising the subfamilies Petromyscinae, Mystromyinae, Cricetomyinae, Nesomyinae, and core dendromurines, (3) a clade of Old World taxa belonging to Murinae, Otomyinae, Gerbillinae, Acomyinae, and Lophiomyinae, (4) a clade uniting the subfamilies Sigmodontinae, Arvicolinae, and Cricetinae, and (5) a unique lineage containing the monotypic Calomyscinae. Although relationships among the latter four clades cannot be resolved, several well-supported supergeneric groupings within each are identified. A preliminary examination of molar tooth morphology on the resulting phylogeny suggests the triserial murid molar pattern as conceived by evolved at least three times during the course of muroid evolution.

“Lipotyphlan” phylogeny based on the growth hormone receptor gene: a reanalysis

From an evolutionary perspective, "insectivores" have been one of the most important mammalian groups for over a century. Morphologists have successively pruned flying lemurs, elephant shrews, and tree shrews from Insectivora, but have retained chrysochlorids, tenrecs, erinaceids, soricids, talpids, and solenodontids in crown-group Lipotyphla. With the appearance of large molecular data sets, the monophyly of Lipotyphla has proved untenable. Rather, an emerging consensus is that Lipotyphla is a diphyletic taxon comprised of two monophyletic groups, Afrosoricida and Eulipotyphla. A recent paper by Malia et al. [Mol. Phylogenet. Evol. 24 (2002) 91-101] challenged this view and argued that "While the data [Growth Hormone Receptor] were unable to support the orders Lipotyphla, Eulipotyphla, and Tenrecoidea [= Afrosoricida] this was most likely due to the polyphyly of these groups and not to problems associated with the gene itself such as saturation or highly divergent sequences em leader " (p. 100). We analyzed Malia et al.'s original GHR data set (at both nuclear and protein level), an expanded GHR data set that included 49 additional sequences, and a concatenated data set that included GHR, BRCA1, vWF, and A2AB for a diverse selection of lipotyphlan taxa. Although protein analyses proved inconclusive, all analyses at the DNA level clearly show that the statement of Malia et al. is erroneous. Indeed, likelihood analyses with GHR and with the concatenated data set provide more support for Eulipotyphla and Afrosoricida than for competing hypotheses. These results also highlight the potential pitfalls of single-gene and parsimony-only analyses.

Exhaustive sample set among Viverridae reveals the sister-group of felids: the linsangs as a case of extreme morphological convergence within Feliformia

Although molecular studies have helped to clarify the phylogeny of the problematic family Viverridae, a recent phylogenetic investigation based on cytochrome b (cyt b) has excluded the Asiatic linsangs (genus Prionodon) from the family. To assess the phylogenetic position of the Asiatic linsangs within the Feliformia, we analysed an exhaustive taxonomic sample set with cyt b and newly produced transthyretin intron I sequences (TR-I-I). TR-I-I alone and cyt b +TR-I-I combined (maximum-likelihood analysis) highly support the position of Asiatic linsangs as sister-group of the Felidae. The estimation of minimum divergence dates from molecular data suggests a splitting event ca. 33.3 million years (Myr) ago, which lends support to historical assertions that the Asiatic linsangs are "living fossils" that share a plesiomorphic morphotype with the Oligocene feliform Paleoprionodon. The African linsang is estimated to appear more than 20 Myr later and represents the sister-group of the genus Genetta. Our phylogenetic results illustrate numerous morphological convergences of "diagnostic" characters among Feliformia that might be problematic for the identification of fossil taxa. The morphotype reappearance from the Asiatic to the African linsangs suggests that the genome of the Feliformia conserved its potential ability of expression for a peculiar adaptive phenotype throughout evolution, in this case arboreality and hypercarnivory in tropical forest.

RNA-based phylogenetic methods: application to mammalian mitochondrial RNA sequences

The PHASE software package allows phylogenetic tree construction with a number of evolutionary models designed specifically for use with RNA sequences that have conserved secondary structure. Evolution in the paired regions of RNAs occurs via compensatory substitutions, hence changes on either side of a pair are correlated. Accounting for this correlation is important for phylogenetic inference because it affects the likelihood calculation. In the present study we use the complete set of tRNA and rRNA sequences from 69 complete mammalian mitochondrial genomes. The likelihood calculation uses two evolutionary models simultaneously for different parts of the sequence: a paired-site model for the paired sites and a single-site model for the unpaired sites. We use Bayesian phylogenetic methods and a Markov chain Monte Carlo algorithm is used to obtain the most probable trees and posterior probabilities of clades. The results are well resolved for almost all the important branches on the mammalian tree. They support the arrangement of mammalian orders within the four supra-ordinal clades that have been identified by studies of much larger data sets mainly comprising nuclear genes. Groups such as the hedgehogs and the murid rodents, which have been problematic in previous studies with mitochondrial proteins, appear in their expected position with the other members of their order. Our choice of genes and evolutionary model appears to be more reliable and less subject to biases caused by variation in base composition than previous studies with mitochondrial genomes.

A new phylogenetic marker, apolipoprotein B, provides compelling evidence for eutherian relationships

Higher-level relationships within, and the root of Placentalia, remain contentious issues. Resolution of the placental tree is important to the choice of mammalian genome projects and model organisms, as well as for understanding the biogeography of the eutherian radiation. We present phylogenetic analyses of 63 species representing all extant eutherian mammal orders for a new molecular phylogenetic marker, a 1.3kb portion of exon 26 of the apolipoprotein B (APOB) gene. In addition, we analyzed a multigene concatenation that included APOB sequences and a previously published data set (Murphy et al., 2001b) of three mitochondrial and 19 nuclear genes, resulting in an alignment of over 17kb for 42 placentals and two marsupials. Due to computational difficulties, previous maximum likelihood analyses of large, multigene concatenations for placental mammals have used quartet puzzling, less complex models of sequence evolution, or phylogenetic constraints to approximate a full maximum likelihood bootstrap. Here, we utilize a Unix load sharing facility to perform maximum likelihood bootstrap analyses for both the APOB and concatenated data sets with a GTR+Gamma+I model of sequence evolution, tree-bisection and reconnection branch-swapping, and no phylogenetic constraints. Maximum likelihood and Bayesian analyses of both data sets provide support for the superordinal clades Boreoeutheria, Euarchontoglires, Laurasiatheria, Xenarthra, Afrotheria, and Ostentoria (pangolins+carnivores), as well as for the monophyly of the orders Eulipotyphla, Primates, and Rodentia, all of which have recently been questioned. Both data sets recovered an association of Hippopotamidae and Cetacea within Cetartiodactyla, as well as hedgehog and shrew within Eulipotyphla. APOB showed strong support for an association of tarsier and Anthropoidea within Primates. Parsimony, maximum likelihood and Bayesian analyses with both data sets placed Afrotheria at the base of the placental radiation. Statistical tests that employed APOB to examine a priori hypotheses for the root of the placental tree rejected rooting on myomorphs and hedgehog, but did not discriminate between rooting at the base of Afrotheria, at the base of Xenarthra, or between Atlantogenata (Xenarthra+Afrotheria) and Boreoeutheria. An orthologous deletion of 363bp in the aligned APOB sequences proved phylogenetically informative for the grouping of the order Carnivora with the order Pholidota into the superordinal clade Ostentoria. A smaller deletion of 237-246bp was diagnostic of the superordinal clade Afrotheria.

Evaluating placental inter-ordinal phylogenies with novel sequences including RAG1, gamma-fibrinogen, ND6, and mt-tRNA, plus MCMC-driven nucleotide, amino acid, and codon models

It is essential to test a priori scientific hypotheses with independent data, not least to partly negate factors such as gene-specific base composition biases misleading our models. Seven new gene segments and sequences plus Bayesian likelihood phylogenetic methods were used to compare and test five recent placental phylogenies. These five phylogenies are similar to each other, yet quite different from Fthose of previously proposed trees, and span Waddell et al. [Syst. Biol. 48 (1999) 1] to Murphy et al. [Science 294 (2001b) 2348]. Trees for RAG1, gamma-fibrinogen, ND6, mt-tRNA, mt-RNA, c-MYC, epsilon -globin, and GHR are significantly congruent with the four main groups of mammals common to the five phylogenies, i.e., Afrotheria, Laurasiatheria, Euarchontoglires, Xenarthra plus Boreoeutheria (Laurasiatheria plus Euarchontoglires). Where these five a priori phylogenies differ, remain areas generally hard to resolve with the new sequences. The root remains ambiguous and does not reject a basal Afrotheria (the Exafroplacentalia hypothesis), Afrotheria plus Xenarthra together with basal (Atlantogenata), or Epitheria (Xenarthra basal) convincingly. Good evidence is found that Eulipotyphla is monophyletic and is located at the base of Laurasiatheria. The shrew mole, Uropsilus, is found to cluster consistently with other moles, while Solenodon may be the sister taxa to all other eulipotyphlans. Support is found for a probable sister pairing of just hedgehogs/gymnures and shrews. Relationships within Afrotheria, except the Paenungulata clade, remain hard to resolve, although there is congruent support for Afroinsectiphillia (aardvark, elephant shrews, golden moles, and tenrecs). A first-time use is made of MCMC enacted general time-reversible (GTR) amino acid and codon-based models for general tree selection. Even with ND6, a GTR amino acid model provided resolution of fine features, such as the sister group relationship of walrus to Otatriidae, and with BRCA a more reasonable rooting. An extensive analysis of GHR sequences reveals strong congruence with prior phylogenies, including strong support for Eulipotyphla, and good resolution within Rodentia. A codon model gives a worse likelihood than a nucleotide model and sometimes switches support, e.g., with RAG1+gamma-fibrinogen from a hyrax-sirenian association to support for Tethytheria. An analysis of the concatenated data is in accordance with well-resolved features of the gene trees. Taken all together, this work suggests that we are on the right path finding strong confirmation of prior phylogenies. However, with the use of robust criteria for assessing trees (i.e., not Bayesian posteriors), it is apparent parts of the tree remain hard to resolve. Since our current models are far from fitting the sequence data, we should continue with our exploratory analyses to arrive at a refined set of hypotheses for future testing using more model independent characters (e.g., rare indels, gene rearrangement, and SINE data).

Reciprocal chromosome painting among human, aardvark, and elephant (superorder Afrotheria) reveals the likely eutherian ancestral karyotype

The Afrotheria, a supraordinal grouping of mammals whose radiation is rooted in Africa, is strongly supported by DNA sequence data but not by their disparate anatomical features. We have used flow-sorted human, aardvark, and African elephant chromosome painting probes and applied reciprocal painting schemes to representatives of two of the Afrotherian orders, the Tubulidentata (aardvark) and Proboscidea (elephants), in an attempt to shed additional light on the evolutionary affinities of this enigmatic group of mammals. Although we have not yet found any unique cytogenetic signatures that support the monophyly of the Afrotheria, embedded within the aardvark genome we find the strongest evidence yet of a mammalian ancestral karyotype comprising 2n = 44. This karyotype includes nine chromosomes that show complete conserved synteny to those of man, six that show conservation as single chromosome arms or blocks in the human karyotype but that occur on two different chromosomes in the ancestor, and seven neighbor-joining combinations (i.e., the synteny is maintained in the majority of species of the orders studied so far, but which corresponds to two chromosomes in humans). The comparative chromosome maps presented between human and these Afrotherian species provide further insight into mammalian genome organization and comparative genomic data for the Afrotheria, one of the four major evolutionary clades postulated for the Eutheria.