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

Nomenclature and placental mammal phylogeny

An issue arising from recent progress in establishing the placental mammal Tree of Life concerns the nomenclature of high-level clades. Fortunately, there are now several well-supported clades among extant mammals that require unambiguous, stable names. Although the International Code of Zoological Nomenclature does not apply above the Linnean rank of family, and while consensus on the adoption of competing systems of nomenclature does not yet exist, there is a clear, historical basis upon which to arbitrate among competing names for high-level mammalian clades. Here, we recommend application of the principles of priority and stability, as laid down by G.G. Simpson in 1945, to discriminate among proposed names for high-level taxa. We apply these principles to specific cases among placental mammals with broad relevance for taxonomy, and close with particular emphasis on the Afrotherian family Tenrecidae. We conclude that no matter how reconstructions of the Tree of Life change in years to come, systematists should apply new names reluctantly, deferring to those already published and maximizing consistency with existing nomenclature.

The new framework for understanding placental mammal evolution

An unprecedented level of confidence has recently crystallized around a new hypothesis of how living placental mammals share a pattern of common descent. The major groups are afrotheres (e.g., aardvarks, elephants), xenarthrans (e.g., anteaters, sloths), laurasiatheres (e.g., horses, shrews), and euarchontoglires (e.g., humans, rodents). Compared with previous hypotheses this tree is remarkably stable; however, some uncertainty persists about the location of the placental root, and (for example) the position of bats within laurasiatheres, of sea cows and aardvarks within afrotheres, and of dermopterans within euarchontoglires. A variety of names for sub-clades within the new placental mammal tree have been proposed, not all of which follow conventions regarding priority and stability. More importantly, the new phylogenetic framework enables the formulation of new hypotheses and testing thereof, for example regarding the possible developmental dichotomy that seems to distinguish members of the newly identified southern and northern radiations of living placental mammals.

Pattern and timing of diversification of the mammalian order Carnivora inferred from multiple nuclear gene sequences

The mammalian order Carnivora has attracted the attention of scientists of various disciplines for decades, leading to intense interest in defining its supra-familial relationships. In the last few years, major changes to the topological structure of the carnivoran tree have been proposed and supported by various molecular data sets, radically changing the traditional view of family composition in this order. Although a sequence of molecular studies have established a growing consensus with respect to most inter-familial relationships, no analysis so far has included all carnivoran lineages (both feliform and caniform) in an integrated data set, so as to determine comparative patterns of diversification. Moreover, no study conducted thus far has estimated divergence dates among all carnivoran families, which is an important requirement in the attempt to understand the patterns and tempo of diversification in this group. In this study, we have investigated the phylogenetic relationships among carnivoran families, and performed molecular dating analyses of the inferred nodes. We assembled a molecular supermatrix containing 14 genes (7765 bp), most of which have not been previously used in supra-familial carnivoran phylogenetics, for 50 different genera representing all carnivoran families. Analysis of this data set led to consistent and robust resolution of all supra-familial nodes in the carnivoran tree, and allowed the construction of a molecular timescale for the evolution of this mammalian order.

Estimating the phylogeny and divergence times of primates using a supermatrix approach

Background

The primates are among the most broadly studied mammalian orders, with the published literature containing extensive analyses of their behavior, physiology, genetics and ecology. The importance of this group in medical and biological research is well appreciated, and explains the numerous molecular phylogenies that have been proposed for most primate families and genera. Composite estimates for the entire order have been infrequently attempted, with the last phylogenetic reconstruction spanning the full range of primate evolutionary relationships having been conducted over a decade ago.

Results

To estimate the structure and tempo of primate evolutionary history, we employed Bayesian phylogenetic methods to analyze data supermatrices comprising 7 mitochondrial genes (6,138 nucleotides) from 219 species across 67 genera and 3 nuclear genes (2,157 nucleotides) from 26 genera. Many taxa were only partially represented, with an average of 3.95 and 5.43 mitochondrial genes per species and per genus, respectively, and 2.23 nuclear genes per genus. Our analyses of mitochondrial DNA place Tarsiiformes as the sister group of Strepsirrhini. Within Haplorrhini, we find support for the primary divergence of Pitheciidae in Platyrrhini, and our results suggest a sister grouping of African and non-African colobines within Colobinae and of Cercopithecini and Papionini within Cercopthecinae. Date estimates for nodes within each family and genus are presented, with estimates for key splits including: Strepsirrhini-Haplorrhini 64 million years ago (MYA), Lemuriformes-Lorisiformes 52 MYA, Platyrrhini-Catarrhini 43 MYA and Cercopithecoidea-Hominoidea 29 MYA.

Conclusion

We present an up-to-date, comprehensive estimate of the structure and tempo of primate evolutionary history. Although considerable gaps remain in our knowledge of the primate phylogeny, increased data sampling, particularly from nuclear loci, will be able to provide further resolution.

Evolution of Picornaviridae: an examination of phylogenetic relationships and cophylogeny

Picornaviruses are responsible for some of the most common and debilitating illnesses affecting humans and animals worldwide. To extend our knowledge of the evolution of picornaviruses and their molecular epidemiology, phylogenetic relationships among 11 genera and the unassigned seal picornavirus type 1 were estimated from the conserved proteins 2C, 3C(pro), and 3D(pol). Each gene was analyzed separately and as a combined dataset. Different tree topologies were recovered from each gene. However, their sequences were determined to be combinable based on our finding of no recombination among genera and failing to reject the hypothesis of homogeneity among datasets using ILD tests. The combined data tree topology was identical to the 3D(pol) gene tree; a topology largely consistent with previous phylogenetic hypotheses based on 3D(pol) and the coding genome. Phylogenetic trees estimated from six phenotypic characters were not congruent with those recovered from molecular datasets; further supporting the hypothesis that viral phenotypes are highly plastic. Finally, we tested the hypothesis of host-virus cophylogeny. Both global and individual tests of the relationships between host and virus trees failed to detect a significant association. These results emphasize the importance of horizontal transmission among host species for picornavirus diversification rather than vertical transmission accompanying speciation.

Breakpoint graphs and ancestral genome reconstructions

Recently completed whole-genome sequencing projects marked the transition from gene-based phylogenetic studies to phylogenomics analysis of entire genomes. We developed an algorithm MGRA for reconstructing ancestral genomes and used it to study the rearrangement history of seven mammalian genomes: human, chimpanzee, macaque, mouse, rat, dog, and opossum. MGRA relies on the notion of the multiple breakpoint graphs to overcome some limitations of the existing approaches to ancestral genome reconstructions. MGRA also generates the rearrangement-based characters guiding the phylogenetic tree reconstruction when the phylogeny is unknown.

Retroposon analysis and recent geological data suggest near-simultaneous divergence of the three superorders of mammals

As a consequence of recent developments in molecular phylogenomics, all extant orders of placental mammals have been grouped into 3 lineages: Afrotheria, Xenarthra, and Boreotheria, which originated in Africa, South America, and Laurasia, respectively. Despite this advancement, the order of divergence of these 3 lineages remains unresolved. Here, we performed extensive retroposon analysis with mammalian genomic data. Surprisingly, we identified a similar number of informative retroposon loci that support each of 3 possible phylogenetic hypotheses: the basal position for Afrotheria (22 loci), Xenarthra (25 loci), and Boreotheria (21 loci). This result indicates that the divergence of the placental common ancestor into the 3 lineages occurred nearly simultaneously. Thus, we examined whether these molecular data could be integrated into the geological context by incorporating recent geological data. We obtained firm evidence that complete separation of Gondwana into Africa and South America occurred 120 +/- 10 Ma. Accordingly, the previous reported time frame (division of Pangea into Gondwana and Laurasia at 148-138 Ma and division of Gondwana at 105 Ma) cannot be used to validate mammalian divergence order. Instead, we use our retroposon results and the recent geological data to propose that near-simultaneous divisions of continents leading to isolated Africa, South America, and Laurasia caused nearly concomitant divergence of the ancient placental ancestor into 3 lineages, Afrotheria, Xenarthra, and Boreotheria, approximately 120 Ma.

Chromosomal rearrangements underlying karyotype differences between Chinese pangolin (Manis pentadactyla) and Malayan pangolin (Manis javanica) revealed by chromosome painting

The Chinese pangolin (Manis pentadactyla), a representative species of the order Pholidota, has been enlisted in the mammalian whole-genome sequencing project mainly because of its phylogenetic importance. Previous studies showed that the diploid number of M. pentadactyla could vary from 2n = 36 to 42. To further characterize the genome organization of M. pentadactyla and to elucidate chromosomal mechanism underlying the karyotype diversity of Pholidota, we flow-sorted the chromosomes of 2n = 40 M. pentadactyla, and generated a set of chromosome-specific probes by DOP-PCR amplification of flow-sorted chromosomes. A comparative chromosome map between M. pentadactyla and the Malayan pangolin (Manis javanica, 2n = 38), as well as between human and M. pentadactyla, was established by chromosome painting for the first time. Our results demonstrate that seven Robertsonian rearrangements, together with considerable variations in the quantity of heterochromatin and in the number of nucleolar organizer regions (NORs) differentiate the karyotypes of 2n = 38 M. javanica and 2n = 40 M. pentadactyla. Moreover, we confirm that the M. javanica Y chromosome bears one NOR. Comparison of human homologous segment associations found in the genomes of M. javanica and M. pentadactyla revealed seven shared associations (HSA 1q/11, 2p/5, 2q/10q, 4p+q/20, 5/13, 6/19p and 8q/10p) that could constitute the potential Pholidota-specific signature rearrangements.

Mosaic retroposon insertion patterns in placental mammals

One and a half centuries after Charles Darwin and Alfred Russel Wallace outlined our current understanding of evolution, a new scientific era is dawning that enables direct observations of genetic variation. However, pure sequence-based molecular attempts to resolve the basal origin of placental mammals have so far resulted only in apparently conflicting hypotheses. By contrast, in the mammalian genomes where they were highly active, the insertion of retroelements and their comparative insertion patterns constitute a neutral, virtually homoplasy-free archive of evolutionary histories. The "presence" of a retroelement at an orthologous genomic position in two species indicates their common ancestry in contrast to its "absence" in more distant species. To resolve the placental origin controversy we extracted approximately 2 million potentially phylogenetically informative, retroposon-containing loci from representatives of the major placental mammalian lineages and found highly significant evidence challenging all current single hypotheses of their basal origin. The Exafroplacentalia hypothesis (Afrotheria as the sister group to all remaining placentals) is significantly supported by five retroposon insertions, the Epitheria hypothesis (Xenarthra as the sister group to all remaining placentals) by nine insertion patterns, and the Atlantogenata hypothesis (a monophyletic clade comprising Xenarthra and Afrotheria as the sister group to Boreotheria comprising all remaining placentals) by eight insertion patterns. These findings provide significant support for a "soft" polytomy of the major mammalian clades. Ancestral successive hybridization events and/or incomplete lineage sorting associated with short speciation intervals are viable explanations for the mosaic retroposon insertion patterns of recent placental mammals and for the futile search for a clear root dichotomy.

Support patterns from different outgroups provide a strong phylogenetic signal

It is known that the accuracy of phylogenetic reconstruction decreases when more distant outgroups are used. We quantify this phenomenon with a novel scoring method, the outgroup score pOG. This score expresses if the support for a particular branch of a tree decreases with increasingly distant outgroups. Large-scale simulations confirmed that the outgroup support follows this expectation and that the pOG score captures this pattern. The score often identifies the correct topology even when the primary reconstruction methods fail, particularly in the presence of model violations. In simulations of problematic phylogenetic scenarios such as rate variation among lineages (which can lead to long-branch attraction artifacts) and quartet-based reconstruction, the pOG analysis outperformed the primary reconstruction methods. Because the pOG method does not make any assumptions about the evolutionary model (besides the decreasing support from increasingly distant outgroups), it can detect cases of violations not treated by a specific model or too strong to be fully corrected. When used as an optimization criterion in the construction of a tree of 23 mammals, the outgroup signal confirmed many well-accepted mammalian orders and superorders. It supports Atlantogenata, a clade of Afrotheria and Xenarthra, and suggests an Artiodactyla-Chiroptera clade.

Estimating maximum likelihood phylogenies with PhyML

Our understanding of the origins, the functions and/or the structures of biological sequences strongly depends on our ability to decipher the mechanisms of molecular evolution. These complex processes can be described through the comparison of homologous sequences in a phylogenetic framework. Moreover, phylogenetic inference provides sound statistical tools to exhibit the main features of molecular evolution from the analysis of actual sequences. This chapter focuses on phylogenetic tree estimation under the maximum likelihood (ML) principle. Phylogenies inferred under this probabilistic criterion are usually reliable and important biological hypotheses can be tested through the comparison of different models. Estimating ML phylogenies is computationally demanding, and careful examination of the results is warranted. This chapter focuses on PhyML, a software that implements recent ML phylogenetic methods and algorithms. We illustrate the strengths and pitfalls of this program through the analysis of a real data set. PhyML v3.0 is available from (http://atgc_montpellier.fr/phyml/).

Suprafamilial relationships among Rodentia and the phylogenetic effect of removing fast-evolving nucleotides in mitochondrial, exon and intron fragments

BACKGROUND

The number of rodent clades identified above the family level is contentious, and to date, no consensus has been reached on the basal evolutionary relationships among all rodent families. Rodent suprafamilial phylogenetic relationships are investigated in the present study using approximately 7600 nucleotide characters derived from two mitochondrial genes (Cytochrome b and 12S rRNA), two nuclear exons (IRBP and vWF) and four nuclear introns (MGF, PRKC, SPTBN, THY). Because increasing the number of nucleotides does not necessarily increase phylogenetic signal (especially if the data is saturated), we assess the potential impact of saturation for each dataset by removing the fastest-evolving positions that have been recognized as sources of inconsistencies in phylogenetics.

RESULTS

Taxonomic sampling included multiple representatives of all five rodent suborders described. Fast-evolving positions for each dataset were identified individually using a discrete gamma rate category and sites belonging to the most rapidly evolving eighth gamma category were removed. Phylogenetic tree reconstructions were performed on individual and combined datasets using Parsimony, Bayesian, and partitioned Maximum Likelihood criteria. Removal of fast-evolving positions enhanced the phylogenetic signal to noise ratio but the improvement in resolution was not consistent across different data types. The results suggested that elimination of fastest sites only improved the support for nodes moderately affected by homoplasy (the deepest nodes for introns and more recent nodes for exons and mitochondrial genes).

CONCLUSION

The present study based on eight DNA fragments supports a fully resolved higher level rodent phylogeny with moderate to significant nodal support. Two inter-suprafamilial associations emerged. The first comprised a monophyletic assemblage containing the Anomaluromorpha (Anomaluridae + Pedetidae) + Myomorpha (Muridae + Dipodidae) as sister clade to the Castorimorpha (Castoridae + Geomyoidea). The second suprafamilial clustering identified a novel association between the Sciuromorpha (Gliridae + (Sciuridae + Aplodontidae)) and the Hystricomorpha (Ctenodactylidae + Hystricognathi) which together represents the earliest dichotomy among Rodentia. Molecular time estimates using a relaxed Bayesian molecular clock dates the appearance of the five suborders nearly contemporaniously at the KT boundary and this is congruent with suggestions of an early explosion of rodent diversity. Based on these newly proposed phylogenetic relationships, the evolution of the zygomasseteric pattern that has been used for a long time in rodent systematics is evaluated.

Hemiplasy and homoplasy in the karyotypic phylogenies of mammals

Phylogenetic reconstructions are often plagued by difficulties in distinguishing phylogenetic signal (due to shared ancestry) from phylogenetic noise or homoplasy (due to character-state convergences or reversals). We use a new interpretive hypothesis, termed hemiplasy, to show how random lineage sorting might account for specific instances of seeming "phylogenetic discordance" among different chromosomal traits, or between karyotypic features and probable species phylogenies. We posit that hemiplasy is generally less likely for underdominant chromosomal polymorphisms (i.e., those with heterozygous disadvantage) than for neutral polymorphisms or especially for overdominant rearrangements (which should tend to be longer-lived), and we illustrate this concept by using examples from chiropterans and afrotherians. Chromosomal states are especially powerful in phylogenetic reconstructions because they offer strong signatures of common ancestry, but their evolutionary interpretations remain fully subject to the principles of cladistics and the potential complications of hemiplasy.

Confirming the phylogeny of mammals by use of large comparative sequence data sets

The ongoing generation of prodigious amounts of genomic sequence data from myriad vertebrates is providing unparalleled opportunities for establishing definitive phylogenetic relationships among species. The size and complexities of such comparative sequence data sets not only allow smaller and more difficult branches to be resolved but also present unique challenges, including large computational requirements and the negative consequences of systematic biases. To explore these issues and to clarify the phylogenetic relationships among mammals, we have analyzed a large data set of over 60 megabase pairs (Mb) of high-quality genomic sequence, which we generated from 41 mammals and 3 other vertebrates. All sequences are orthologous to a 1.9-Mb region of the human genome that encompasses the cystic fibrosis transmembrane conductance regulator gene (CFTR). To understand the characteristics and challenges associated with phylogenetic analyses of such a large data set, we partitioned the sequence data in several ways and utilized maximum likelihood, maximum parsimony, and Neighbor-Joining algorithms, implemented in parallel on Linux clusters. These studies yielded well-supported phylogenetic trees, largely confirming other recent molecular phylogenetic analyses. Our results provide support for rooting the placental mammal tree between Atlantogenata (Xenarthra and Afrotheria) and Boreoeutheria (Euarchontoglires and Laurasiatheria), illustrate the difficulty in resolving some branches even with large amounts of data (e.g., in the case of Laurasiatheria), and demonstrate the valuable role that very large comparative sequence data sets can play in refining our understanding of the evolutionary relationships of vertebrates.

The molecular evolution of acrosin in placental mammals

Acrosin is thought to fulfill several different roles in fertilization including that of a serine protease and in secondary zona pellucida (ZP) binding. However, acrosin's importance as a fertilization protein has been questioned. Especially since it was discovered that acrosin knockout mice are fertile. In this study, we explored the sites involved in serine protease activity and secondary binding. We also assessed conservation in functional sites across species and examined whether amino acid changes present in the human population have the potential to affect fertility. In addition, since many mammalian reproduction proteins have been found to evolve rapidly, we tested for positive selection. Sequences from 43 mammals from all 19 placental orders, which included a total of 828 nucleotides from acrosin exons 2, 3, 4, and a portion of exon 5, were obtained. We found that all sites of the serine catalytic triad as well as three other sites linked to catalytic activity were completely conserved. Five of six sites proposed to play a role in secondary binding were 100% conserved as basic residues. These results support an evolutionary conserved role for acrosin as a serine protease and secondary binding protein across placental mammals. We found statistically significant support for positive selection within acrosin, but no single amino acid site reached the significance level of P > 0.95 for inclusion within the category omega > 1. Based upon two amino acid mutation scoring systems, three out of seven human residue changing single nucleotide polymorphisms (SNPs) were found to be potentially protein-altering mutations.

Rooting the eutherian tree: the power and pitfalls of phylogenomics

In an attempt to root the eutherian tree using genome-scale data with the maximum likelihood method, a concatenate analysis supports a putatively wrong tree, whereas separate analyses of different genes reduced the bias.

Background

Ongoing genome sequencing projects have led to a phylogenetic approach based on genome-scale data (phylogenomics), which is beginning to shed light on longstanding unresolved phylogenetic issues. The use of large datasets in phylogenomic analysis results in a global increase in resolution due to a decrease in sampling error. However, a fully resolved tree can still be wrong if the phylogenetic inference is biased.

Results

Here, in an attempt to root the eutherian tree using genome-scale data with the maximum likelihood method, we demonstrate a case in which a concatenate analysis strongly supports a putatively wrong tree, whereas the total evaluation of separate analyses of different genes grossly reduced the bias of the phylogenetic inference. A conventional method of concatenate analysis of nucleotide sequences from our dataset, which includes a more than 1 megabase alignment of 2,789 nuclear genes, suggests a misled monophyly of Afrotheria (for example, elephant) and Xenarthra (for example, armadillo) with 100% bootstrap probability. However, this tree is not supported by our 'separate method', which takes into account the different tempos and modes of evolution among genes, and instead the basal Afrotheria tree is favored.

Conclusion

Our analysis demonstrates that in cases in which there is great variation in evolutionary features among different genes, the separate model, rather than the concatenate model, should be used for phylogenetic inference, especially in genome-scale data.

Stable isotope evidence for an amphibious phase in early proboscidean evolution

The order Proboscidea includes extant elephants and their extinct relatives and is closely related to the aquatic sirenians (manatees and dugongs) and terrestrial hyracoids (hyraxes). Some analyses of embryological, morphological, and paleontological data suggest that proboscideans and sirenians shared an aquatic or semiaquatic common ancestor, but independent tests of this hypothesis have proven elusive. Here we test the hypothesis of an aquatic ancestry for advanced proboscideans by measuring delta(18)O in tooth enamel of two late Eocene proboscidean genera, Barytherium and Moeritherium, which are sister taxa of Oligocene-to-Recent proboscideans. The combination of low delta(18)O values and low delta(18)O standard deviations in Barytherium and Moeritherium matches the isotopic pattern seen in aquatic and semiaquatic mammals, and differs from that of terrestrial mammals. delta(13)C values of these early proboscideans suggest that both genera are likely to have consumed freshwater plants, although a component of C(3) terrestrial vegetation cannot be ruled out. The simplest explanation for the combined evidence from isotopes, dental functional morphology, and depositional environments is that Barytherium and Moeritherium were at least semiaquatic and lived in freshwater swamp or riverine environments, where they grazed on freshwater vegetation. These results lend new support to the hypothesis that Oligocene-to-Recent proboscideans are derived from amphibious ancestors.

Molecular phylogeny and divergence times of Malagasy tenrecs: influence of data partitioning and taxon sampling on dating analyses

BACKGROUND

Malagasy tenrecs belong to the Afrotherian clade of placental mammals and comprise three subfamilies divided in eight genera (Tenrecinae: Tenrec, Echinops, Setifer and Hemicentetes; Oryzorictinae: Oryzorictes, Limnogale and Microgale; Geogalinae:Geogale). The diversity of their morphology and incomplete taxon sampling made it difficult until now to resolve phylogenies based on either morphology or molecular data for this group. Therefore, in order to delineate the evolutionary history of this family, phylogenetic and dating analyses were performed on a four nuclear genes dataset (ADRA2B, AR, GHR and vWF) including all Malagasy tenrec genera. Moreover, the influence of both taxon sampling and data partitioning on the accuracy of the estimated ages were assessed.

RESULTS

Within Afrotheria the vast majority of the nodes received a high support, including the grouping of hyrax with sea cow and the monophyly of both Afroinsectivora (Macroscelidea + Afrosoricida) and Afroinsectiphillia (Tubulidentata + Afroinsectivora). Strongly supported relationships were also recovered among all tenrec genera, allowing us to firmly establish the grouping of Geogale with Oryzorictinae, and to confirm the previously hypothesized nesting of Limnogale within the genus Microgale. The timeline of Malagasy tenrec diversification does not reflect a fast adaptive radiation after the arrival on Madagascar, indicating that morphological specializations have appeared over the whole evolutionary history of the family, and not just in a short period after colonization. In our analysis, age estimates at the root of a clade became older with increased taxon sampling of that clade. Moreover an augmentation of data partitions resulted in older age estimates as well, whereas standard deviations increased when more extreme partition schemes were used.

CONCLUSION

Our results provide as yet the best resolved gene tree comprising all Malagasy tenrec genera, and may lead to a revision of tenrec taxonomy. A timeframe of tenrec evolution built on the basis of this solid phylogenetic framework showed that morphological specializations of the tenrecs may have been affected by environmental changes caused by climatic and/or subsequent colonization events. Analyses including various taxon sampling and data partitions allow us to point out some possible pitfalls that may lead to biased results in molecular dating; however, further analyses are needed to corroborate these observations.

Multigene phylogeny of the Mustelidae: resolving relationships, tempo and biogeographic history of a mammalian adaptive radiation

BACKGROUND

Adaptive radiation, the evolution of ecological and phenotypic diversity from a common ancestor, is a central concept in evolutionary biology and characterizes the evolutionary histories of many groups of organisms. One such group is the Mustelidae, the most species-rich family within the mammalian order Carnivora, encompassing 59 species classified into 22 genera. Extant mustelids display extensive ecomorphological diversity, with different lineages having evolved into an array of adaptive zones, from fossorial badgers to semi-aquatic otters. Mustelids are also widely distributed, with multiple genera found on different continents. As with other groups that have undergone adaptive radiation, resolving the phylogenetic history of mustelids presents a number of challenges because ecomorphological convergence may potentially confound morphologically based phylogenetic inferences, and because adaptive radiations often include one or more periods of rapid cladogenesis that require a large amount of data to resolve.

RESULTS

We constructed a nearly complete generic-level phylogeny of the Mustelidae using a data matrix comprising 22 gene segments (approximately 12,000 base pairs) analyzed with maximum parsimony, maximum likelihood and Bayesian inference methods. We show that mustelids are consistently resolved with high nodal support into four major clades and three monotypic lineages. Using Bayesian dating techniques, we provide evidence that mustelids underwent two bursts of diversification that coincide with major paleoenvironmental and biotic changes that occurred during the Neogene and correspond with similar bursts of cladogenesis in other vertebrate groups. Biogeographical analyses indicate that most of the extant diversity of mustelids originated in Eurasia and mustelids have colonized Africa, North America and South America on multiple occasions.

CONCLUSION

Combined with information from the fossil record, our phylogenetic and dating analyses suggest that mustelid diversification may have been spurred by a combination of faunal turnover events and diversification at lower trophic levels, ultimately caused by climatically driven environmental changes. Our biogeographic analyses show Eurasia as the center of origin of mustelid diversity and that mustelids in Africa, North America and South America have been assembled over time largely via dispersal, which has important implications for understanding the ecology of mustelid communities.

OrthoMaM: a database of orthologous genomic markers for placental mammal phylogenetics

Background

Molecular sequence data have become the standard in modern day phylogenetics. In particular, several long-standing questions of mammalian evolutionary history have been recently resolved thanks to the use of molecular characters. Yet, most studies have focused on only a handful of standard markers. The availability of an ever increasing number of whole genome sequences is a golden mine for modern systematics. Genomic data now provide the opportunity to select new markers that are potentially relevant for further resolving branches of the mammalian phylogenetic tree at various taxonomic levels.

Description

The EnsEMBL database was used to determine a set of orthologous genes from 12 available complete mammalian genomes. As targets for possible amplification and sequencing in additional taxa, more than 3,000 exons of length > 400 bp have been selected, among which 118, 368, 608, and 674 are respectively retrieved for 12, 11, 10, and 9 species. A bioinformatic pipeline has been developed to provide evolutionary descriptors for these candidate markers in order to assess their potential phylogenetic utility. The resulting OrthoMaM (Orthologous Mammalian Markers) database can be queried and alignments can be downloaded through a dedicated web interface .

Conclusion

The importance of marker choice in phylogenetic studies has long been stressed. Our database centered on complete genome information now makes possible to select promising markers to a given phylogenetic question or a systematic framework by querying a number of evolutionary descriptors. The usefulness of the database is illustrated with two biological examples. First, two potentially useful markers were identified for rodent systematics based on relevant evolutionary parameters and sequenced in additional species. Second, a complete, gapless 94 kb supermatrix of 118 orthologous exons was assembled for 12 mammals. Phylogenetic analyses using probabilistic methods unambiguously supported the new placental phylogeny by retrieving the monophyly of Glires, Euarchontoglires, Laurasiatheria, and Boreoeutheria. Muroid rodents thus do not represent a basal placental lineage as it was mistakenly reasserted in some recent phylogenomic analyses based on fewer taxa. We expect the OrthoMaM database to be useful for further resolving the phylogenetic tree of placental mammals and for better understanding the evolutionary dynamics of their genomes, i.e., the forces that shaped coding sequences in terms of selective constraints.

A new estimate of afrotherian phylogeny based on simultaneous analysis of genomic, morphological, and fossil evidence

BACKGROUND

The placental mammalian clade Afrotheria is now supported by diverse forms of genomic data, but interordinal relationships within, and morphological support for, the group remains elusive. As a means for addressing these outstanding problems, competing hypotheses of afrotherian interordinal relationships were tested through simultaneous parsimony analysis of a large data set (> 4,590 parsimony informative characters) containing genomic data (> 17 kb of nucleotide data, chromosomal associations, and retroposons) and 400 morphological characters scored across 16 extant and 35 extinct afrotherians.

RESULTS

Parsimony analysis of extant taxa alone recovered the interordinal topology (Afrosoricida, ((Macroscelidea, Tubulidentata), (Hyracoidea, (Proboscidea, Sirenia)))). Analysis following addition of extinct taxa instead supported Afroinsectivora (Afrosoricida + Macroscelidea) and Pseudoungulata (Tubulidentata + Paenungulata), as well as Tethytheria (Proboscidea + Sirenia). This latter topology is, however, sensitive to taxon deletion and different placements of the placental root, and numerous alternative interordinal arrangements within Afrotheria could not be statistically rejected. Relationships among extinct stem members of each afrotherian clade were more stable, but one alleged stem macroscelidean (Herodotius) never grouped with that clade and instead consistently joined pseudoungulates or paenungulates. When character transformations were optimized onto a less resolved afrotherian tree that reflects uncertainty about the group's interordinal phylogeny, a total of 21 morphological features were identified as possible synapomorphies of crown Afrotheria, 9 of which optimized unambiguously across all character treatments and optimization methods.

CONCLUSION

Instability in afrotherian interordinal phylogeny presumably reflects rapid divergences during two pulses of cladogenesis - the first in the Late Cretaceous, at and just after the origin of crown Afrotheria, and the second in the early Cenozoic, with the origin of crown Paenungulata. Morphological evidence for divergences during these two pulses either never existed or has largely been "erased" by subsequent evolution along long ordinal branches. There may, nevertheless, be more morphological character support for crown Afrotheria than is currently assumed; the features identified here as possible afrotherian synapomorphies can be further scrutinized through future phylogenetic analyses with broader taxon sampling, as well as recovery of primitive fossil afrotherians from the Afro-Arabian landmass, where the group is likely to have first diversified.

Chromosomal instability in Afrotheria: fragile sites, evolutionary breakpoints and phylogenetic inference from genome sequence assemblies

BACKGROUND

Extant placental mammals are divided into four major clades (Laurasiatheria, Supraprimates, Xenarthra and Afrotheria). Given that Afrotheria is generally thought to root the eutherian tree in phylogenetic analysis of large nuclear gene data sets, the study of the organization of the genomes of afrotherian species provides new insights into the dynamics of mammalian chromosomal evolution. Here we test if there are chromosomal bands with a high tendency to break and reorganize in Afrotheria, and by analyzing the expression of aphidicolin-induced common fragile sites in three afrotherian species, whether these are coincidental with recognized evolutionary breakpoints.

RESULTS

We described 29 fragile sites in the aardvark (OAF) genome, 27 in the golden mole (CAS), and 35 in the elephant-shrew (EED) genome. We show that fragile sites are conserved among afrotherian species and these are correlated with evolutionary breakpoints when compared to the human (HSA) genome. Inddition, by computationally scanning the newly released opossum (Monodelphis domestica) and chicken sequence assemblies for use as outgroups to Placentalia, we validate the HSA 3/21/5 chromosomal synteny as a rare genomic change that defines the monophyly of this ancient African clade of mammals. On the other hand, support for HSA 1/19p, which is also thought to underpin Afrotheria, is currently ambiguous.

CONCLUSION

We provide evidence that (i) the evolutionary breakpoints that characterise human syntenies detected in the basal Afrotheria correspond at the chromosomal band level with fragile sites, (ii) that HSA 3p/21 was in the amniote ancestor (i.e., common to turtles, lepidosaurs, crocodilians, birds and mammals) and was subsequently disrupted in the lineage leading to marsupials. Its expansion to include HSA 5 in Afrotheria is unique and (iii) that its fragmentation to HSA 3p/21 + HSA 5/21 in elephant and manatee was due to a fission within HSA 21 that is probably shared by all Paenungulata.

Waves of genomic hitchhikers shed light on the evolution of gamebirds (Aves: Galliformes)

BACKGROUND

The phylogenetic tree of Galliformes (gamebirds, including megapodes, currassows, guinea fowl, New and Old World quails, chicken, pheasants, grouse, and turkeys) has been considerably remodeled over the last decades as new data and analytical methods became available. Analyzing presence/absence patterns of retroposed elements avoids the problems of homoplastic characters inherent in other methodologies. In gamebirds, chicken repeats 1 (CR1) are the most prevalent retroposed elements, but little is known about the activity of their various subtypes over time. Ascertaining the fixation patterns of CR1 elements would help unravel the phylogeny of gamebirds and other poorly resolved avian clades.

RESULTS

We analyzed 1,978 nested CR1 elements and developed a multidimensional approach taking advantage of their transposition in transposition character (TinT) to characterize the fixation patterns of all 22 known chicken CR1 subtypes. The presence/absence patterns of those elements that were active at different periods of gamebird evolution provided evidence for a clade (Cracidae + (Numididae + (Odontophoridae + Phasianidae))) not including Megapodiidae; and for Rollulus as the sister taxon of the other analyzed Phasianidae. Genomic trace sequences of the turkey genome further demonstrated that the endangered African Congo Peafowl (Afropavo congensis) is the sister taxon of the Asian Peafowl (Pavo), rejecting other predominantly morphology-based groupings, and that phasianids are monophyletic, including the sister taxa Tetraoninae and Meleagridinae.

CONCLUSION

The TinT information concerning relative fixation times of CR1 subtypes enabled us to efficiently investigate gamebird phylogeny and to reconstruct an unambiguous tree topology. This method should provide a useful tool for investigations in other taxonomic groups as well.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Calibration Choice, Rate Smoothing, and the Pattern of Tetrapod Diversification According to the Long Nuclear Gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.

Phylogenomic data analyses provide evidence that Xenarthra and Afrotheria are sister groups

The phylogenetic positions of the 4 clades, Euarchontoglires, Laurasiatheria, Afrotheria, and Xenarthra, have been major issues in the recent discussion of basal relationships among placental mammals. However, despite considerable efforts these relationships, crucial to the understanding of eutherian evolution and biogeography, have remained essentially unresolved. Euarchontoglires and Laurasiatheria are generally joined into a common clade (Boreoeutheria), whereas the position of Afrotheria and Xenarthra relative to Boreoeutheria has been equivocal in spite of the use of comprehensive amounts of nuclear encoded sequences or the application of genome-level characters such as retroposons. The probable reason for this uncertainty is that the divergences took place long time ago and within a narrow temporal window, leaving only short common branches. With the aim of further examining basal eutherian relationships, we have collected conserved protein-coding sequences from 11 placental mammals, a marsupial and a bird, whose nuclear genomes have been largely sequenced. The length of the alignment of homologous sequences representing each individual species is 2,168,859 nt. This number of sites, representing 2840 protein-coding genes, exceeds by a considerable margin that of any previous study. The phylogenetic analysis joined Xenarthra and Afrotheria on a common branch, Atlantogenata. This topology was found to fit the data significantly better than the alternative trees.

Placentation in dolphins from the Amazon River Basin: the Boto, Inia geoffrensis, and the Tucuxi, Sotalia fluviatilis

A recent reassessment of the phylogenetic affinities of cetaceans makes it timely to compare their placentation with that of the artiodactyls. We studied the placentae of two sympatric species of dolphin from the Amazon River Basin, representing two distinct families. The umbilical cord branched to supply a bilobed allantoic sac. Small blood vessels and smooth muscle bundles were found within the stroma of the cord. Foci of squamous metaplasia occurred in the allanto-amnion and allantochorion. The interhemal membrane of the placenta was of the epitheliochorial type. Two different types of trophoblastic epithelium were seen. Most was of the simple columnar type and indented by fetal capillaries. However, there were also areolar regions with tall columnar trophoblast and these were more sparsely supplied with capillaries. The endometrium was well vascularised and richly supplied with actively secreting glands. These findings are consistent with the current view that Cetacea are nested within Artiodactyla as sister group to the hippopotamids.

Chromosome painting among Proboscidea, Hyracoidea and Sirenia: support for Paenungulata (Afrotheria, Mammalia) but not Tethytheria

Despite marked improvements in the interpretation of systematic relationships within Eutheria, particular nodes, including Paenungulata (Hyracoidea, Sirenia and Proboscidea), remain ambiguous. The combination of a rapid radiation, a deep divergence and an extensive morphological diversification has resulted in a limited phylogenetic signal confounding resolution within this clade both at the morphological and nucleotide levels. Cross-species chromosome painting was used to delineate regions of homology between Loxodonta africana (2n=56), Procavia capensis (2n=54), Trichechus manatus latirostris (2n=48) and an outgroup taxon, the aardvark (Orycteropus afer, 2n=20). Changes specific to each lineage were identified and although the presence of a minimum of 11 synapomorphies confirmed the monophyly of Paenungulata, no change characterizing intrapaenungulate relationships was evident. The reconstruction of an ancestral paenungulate karyotype and the estimation of rates of chromosomal evolution indicate a reduced rate of genomic repatterning following the paenungulate radiation. In comparison to data available for other mammalian taxa, the paenungulate rate of chromosomal evolution is slow to moderate. As a consequence, the absence of a chromosomal character uniting two paenungulates (at the level of resolution characterized in this study) may be due to a reduced rate of chromosomal change relative to the length of time separating successive divergence events.

Early history of mammals is elucidated with the ENCODE multiple species sequencing data

Understanding the early evolution of placental mammals is one of the most challenging issues in mammalian phylogeny. Here, we addressed this question by using the sequence data of the ENCODE consortium, which include 1% of mammalian genomes in 18 species belonging to all main mammalian lineages. Phylogenetic reconstructions based on an unprecedented amount of coding sequences taken from 218 genes resulted in a highly supported tree placing the root of Placentalia between Afrotheria and Exafroplacentalia (Afrotheria hypothesis). This topology was validated by the phylogenetic analysis of a new class of genomic phylogenetic markers, the conserved noncoding sequences. Applying the tests of alternative topologies on the coding sequence dataset resulted in the rejection of the Atlantogenata hypothesis (Xenarthra grouping with Afrotheria), while this test rejected the second alternative scenario, the Epitheria hypothesis (Xenarthra at the base), when using the noncoding sequence dataset. Thus, the two datasets support the Afrotheria hypothesis; however, none can reject both of the remaining topological alternatives.

Application of molecular phylogenetic methods drastically changed the conception of relationships within mammals. Recent molecular phylogenetic studies have shown that living placental mammals belong to one of the three subgroups: Boreoeutheria, Afrotheria, or Xenarthra, but the relations between these are still unknown. In a previous analysis using 16 genes, Boreoeutheria and Xenarthra grouped together. However, a study based on LINE insertions supported the grouping of Boreoeutheria and Afrotheria. To resolve this discrepancy, we applied sequence data from 1% of a genome in a subset of 18 mammalian species. We used concatenated coding sequence data from 218 genes encompassing 205 kilobases of DNA sequence. Phylogenetic analyses have shown Afrotheria as a basal group of Placentalia with high statistical support. To further validate these results, we analyzed a new phylogenetic marker: conserved noncoding sequence alignments (430 kilobases), which resulted in the same position of the placental root. Topological tests rejected the possibility of Afrotheria-Xenarthra grouping with the coding sequence dataset and Boreoeutheria-Afrotheria grouping with the noncoding sequence dataset. Ascertaining the relationships between mammals is of great importance for the investigation of evolutionary behavior of the different functional genomic elements.

Phylogeny of the Procyonidae (Mammalia: Carnivora): Molecules, morphology and the Great American Interchange

The Procyonidae (Mammalia: Carnivora) have played a central role in resolving the controversial systematics of the giant and red pandas, but phylogenetic relationships of species within the family itself have received much less attention. Cladistic analyses of morphological characters conducted during the last two decades have resulted in topologies that group ecologically and morphologically similar taxa together. Specifically, the highly arboreal and frugivorous kinkajou (Potos flavus) and olingos (Bassaricyon) define one clade, whereas the more terrestrial and omnivorous coatis (Nasua), raccoons (Procyon), and ringtails (Bassariscus) define another clade, with the similar-sized Nasua and Procyon joined as sister taxa in this latter group. These relationships, however, have not been tested with molecular sequence data. We examined procyonid phylogenetics based on combined data from nine nuclear and two mitochondrial gene segments totaling 6534bp. We were able to fully resolve relationships within the family with strongly supported and congruent results from maximum parsimony, maximum likelihood, minimum evolution, and Bayesian analyses. We identified three distinct lineages within the family: a (Nasua, Bassaricyon) clade, a (Bassariscus, Procyon) clade, and a Potos lineage, the last of which is sister to the other two clades. These findings, which are in strong disagreement with prior fossil and morphology-based assessments of procyonid relationships, reemphasize the morphological and ecological flexibility of these taxa. In particular, morphological similarities between unrelated genera possibly reflect convergence associated with similar lifestyles and diets rather than ancestry. Furthermore, incongruence between the molecular supermatrix and a morphological character matrix comprised mostly of dental characters [Baskin, J.A., 2004. Bassariscus and Probassariscus (Mammalia, Carnivora, Procyonidae) from the early Barstovian (Middle Miocene). J. Vert. Paleo. 24, 709-720] may be due to non-independence among atomized dental characters that does not take into account the high developmental genetic correlation of these characters. Finally, molecular divergence dating analyses using a relaxed molecular clock approach suggest that intergeneric and intrageneric splits in the Procyonidae mostly occurred in the Miocene. The inferred divergence times for intrageneric splits for several genera whose ranges are bisected by the Panamanian Isthmus is significant because they suggest diversification well precedes the Great American Interchange, which has long been considered a primary underlying mechanism for procyonid evolution.

Using genomic data to unravel the root of the placental mammal phylogeny

The phylogeny of placental mammals is a critical framework for choosing future genome sequencing targets and for resolving the ancestral mammalian genome at the nucleotide level. Despite considerable recent progress defining superordinal relationships, several branches remain poorly resolved, including the root of the placental tree. Here we analyzed the genome sequence assemblies of human, armadillo, elephant, and opossum to identify informative coding indels that would serve as rare genomic changes to infer early events in placental mammal phylogeny. We also expanded our species sampling by including sequence data from >30 ongoing genome projects, followed by PCR and sequencing validation of each indel in additional taxa. Our data provide support for a sister-group relationship between Afrotheria and Xenarthra (the Atlantogenata hypothesis), which is in turn the sister-taxon to Boreoeutheria. We failed to recover any indels in support of a basal position for Xenarthra (Epitheria), which is suggested by morphology and a recent retroposon analysis, or a hypothesis with Afrotheria basal (Exafricoplacentalia), which is favored by phylogenetic analysis of large nuclear gene data sets. In addition, we identified two retroposon insertions that also support Atlantogenata and none for the alternative hypotheses. A revised molecular timescale based on these phylogenetic inferences suggests Afrotheria and Xenarthra diverged from other placental mammals approximately 103 (95-114) million years ago. We discuss the impacts of this topology on earlier phylogenetic reconstructions and repeat-based inferences of phylogeny.

Genome size: a novel genomic signature in support of Afrotheria

Molecular phylogenetic analyses suggest an emerging phylogeny for the extant Placentalia (eutherian) that radically departs from morphologically based constructions of the past. Placental mammals are partitioned into four supraordinal clades: Afrotheria, Xenarthra, Laurasiatheria, and Euarchontoglires. Afrotheria form an endemic African clade that includes elephant shrews, golden moles, tenrecs, aardvarks, hyraxes, elephants, dugongs, and manatees. Datamining databases of genome size (GS) shows that till today just one afrotherian GS has been evaluated, that of the aardvark Orycteropus afer. We show that the GSs of six selected representatives across the Afrotheria supraordinal group are among the highest for the extant Placentalia, providing a novel genomic signature of this enigmatic group. The mean GS value of Afrotheria, 5.3 +/- 0.7 pg, is the highest reported for the extant Placentalia. This should assist in planning new genome sequencing initiatives.

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.

Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent Anomalurus

Background

Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious.

Results

We have sequenced and performed extensive evolutionary analyses on the mitochondrial genome of the scaly-tailed flying squirrel Anomalurus sp., an enigmatic rodent whose phylogenetic affinities have been obscure and extensively debated. Our phylogenetic analyses of the coding regions of available complete mitochondrial genome sequences from Euarchontoglires suggest that Anomalurus is a sister taxon to the Hystricognathi, and that this clade represents the most basal divergence among sampled Rodentia. Bayesian dating methods incorporating a relaxed molecular clock provide divergence-time estimates which are consistently in agreement with the fossil record and which indicate a rapid radiation within Glires around 60 million years ago.

Conclusion

Taken together, the data presented provide a working hypothesis as to the phylogenetic placement of Anomalurus, underline the utility of mitochondrial sequences in the resolution of even relatively deep divergences and go some way to explaining the difficulty of conclusively resolving higher-level relationships within Glires with available data and methodologies.