BACKGROUND
recently the resolution of the mammalian tree made a quantum leap forward with the analysis of protein coding sequences of whole genome data  <cit> .
the genome data allowed collecting the protein coding sequences from some  <dig> genes, which is equivalent to  <dig>  million nucleotides , representing ≈10% of all coding sequences  <cit> .
during the last  <dig> years since the influential review of novacek  <cit>  on the mammalian evolution, many previously uncertain relationships are now becoming consistently resolved by different data sets and analytical approaches.
however, lately released mammal genome data have not been investigated by phylogenomic analyses.
presently the monophyly of each of the four major clades, euarchontoglires, laurasiatheria, xenarthra, and afrotheria is supported by most analyses based on sequence data  <cit>  and genome level data such as line and ltr elements  <cit> .
the smallest clade, xenarthra, consists of only two orders, cingulata  and pilosa   <cit> .
the afrotheria  <cit>  comprises the traditional orders sirenia, proboscidea, hyracoidea, and tubulidentata plus the members from the paraphyletic lipotyphla : the tenrecidae, macroscelidea and chrysochloridea.
the afrotheria and xenarthra are grouped in the clade named xenafrotheria  <cit> , a clade that is supported by phylogenomic analyses  <cit>  and retroposon and indel analysis  <cit> .
other studies, however, find support for a basal xenarthra clade from retroposon data  <cit> , or a basal afrotheria clade from some sequence data analyses  <cit> .
the two remaining clades, euarchontoglires  and laurasiatheria  together comprise the species rich taxon boreoplacentalia, which is supported by all current analyses.
the new name boreoplacentalia has been suggested for this clade for being more consistent and specific than the previous name boreoeutheria  <cit> .
many branches of the mammalian tree remained difficult to resolve, even with the analysis of several  <dig> thousand nt  of sequence data.
the most basal placental mammal divergences were inconclusive even with the analysis of some  <dig> knt of coding sequence data  <cit> .
the analysis of a ten-fold larger dataset,  <dig>  mnt of protein coding sequence data supported the grouping of xenarthra with the afrotheria in the higher order clade xenafrotheria  <cit> .
jackknife analysis showed that at least  <dig> knt sites of protein coding data were needed to significantly resolve this relationship  <cit> .
the notoriously difficult to resolve split between xenafrotheria and boreoplacentalia took place at ≈ <dig> million years ago .
however, the two clades diverged into xenarthra and afrotheria, and laurasiatheria and euarchontoglires respectively only 3– <dig> million years  later  <cit> .
the resolution of these narrow temporal occurrences required more data than available to the study of nikolaev et al.
<cit> , however the placement of xenarthra and afrotheria remains challenged  <cit> .
in comparison, more distant branching events such as that of monotreme, marsupial and placental mammals can be significantly resolved with as little as  <dig> aa sites of nuclear coding genes  <cit> .
molecular dating on phylogenomic data estimated that 30– <dig> myr separated the splits among the three mammalian infra-classes  <cit> .
this leaves enough time to accumulate phylogenetically relevant sequence differences and to reduce the effects of lineage sorting and prevents introgression.
while the monophyly of the clade laurasiatheria is well established, phylogenomic analysis did not yet resolve the relationships within this clade.
the internal branches of dog, cow, bat, and shrew received only 86% and 93% bootstrap support, respectively, in a recent phylogenomic analyses  <cit> .
this non-significant support may be due to the "insufficient" amount of  <dig> knt of coding sequence that was available at the time of the analysis.
however, the splits may also have occurred at a very narrow time window of a few millions years or less, and more data than expected may be required to resolve such tight divergences.
it is noteworthy that the few million years between mammal ordinal divergences are in the order of the average duration of mammalian species and their typical speciation times  <cit> .
the fact that the splits between the extant orders occurred several tens of million years ago does not diminish the problem that these splits represent divergences among species or even lineages that were at the brink of becoming genetically separated species.
therefore, speciation adds an additional level of complication for resolving such narrow divergences.
not only is the time for such divergences relatively short to accumulate sufficient sequence differences that are needed to resolve these splits, but lineage sorting  <cit> , species hybridization  <cit>  and hybrid speciation  <cit>  can make the resolution of significant parts of the mammalian tree all but impossible.
the continuously increasing amount of sequence data from protein coding genes now makes it possible to provide a phylogenomic analysis of placental mammals for  <dig> species on the basis of several million characters long alignments.
the quantity of new genome data also allows using a new and more stringent assembly of the data, which aims at avoiding data from gene families and thus paralogous sequences in the alignment.
we investigate if the resolution of the basal split among placental mammals  <cit>  prevails with the increased taxon sampling and if other difficult to resolve divergences such as those within the laurasiatheria and euarchontoglires can be resolved with statistical confidence with an increased amount of phylogenomic data and new analytical approaches  <cit> .
RESULTS
the first step in the data collection, searching for human protein-coding genes that do not have obviously closely related members of a gene family resulted in a set of  <dig>  human sequences, from the original  <dig>  unique genes.
these sequences were used for identifying related sequences from the remaining  <dig> species.
after selecting against genes from multi-gene families in any of the species  and using only alignments in which more than  <dig> species are represented,  <dig>  multiple sequence alignments remained.
from this data set sequence alignments that were shorter than  <dig> bp or observed aa distances  larger than 30% for any species pair were removed.
this left a data set of  <dig>  multiple sequence alignments.
the total size of this alignment was  <dig> , <dig> nucleotides for  <dig> species.
the average length of the individual sequences is  <dig> ±  <dig> nt and the average p distance between human and platypus is  <dig>  ±  <dig> %, with a maximum of 30% , which is expected from the filtering constraints.
the rigorous filtering by recursive blast search removed sequences that are part of gene families notably reduced the amount of data.
the filtering was rather effective, even though there is no formal proof for its efficiency yet.
already after the first filtering step  sequences of the hemoglobin, myosin, olfactory receptor, hox gene and other sequences that are part of known gene families were eliminated from the data set.
it is conceivable that by this approach data from other, less known gene families were also excluded from the analysis.
the recursive blast search with a cutoff value of 10- <dig> identified genes in gene families that had more than 75% sequence similarity and eliminated these.
it can be argued that more distant sequences will not be found when searching for such sequences in other species with the same cutoff value of 10- <dig>
when inspected by eye the base and amino acid frequencies appear to be very similar among the species.
however, due to the large size of the dataset the chi-square test may become overly strict.
compositional homogeneity was rejected for many species both for nucleotide and amino acid data .
the compositional heterogeneity remained even when the data were r/y-coded.
a species for which compositional homogeneity could not be rejected with 5% significance, and are thus assumed to be homogeneous.
the evolutionary model and the number of rate heterogeneity classes were estimated from previous phylogenomic analyses  <cit> .
exact tests such as implemented in modeltest  <cit>  could not be performed due to computational constraints and the inability of most programs to analyze a data set of this magnitude.
a manual analysis would have been prohibitively time consuming to perform and we suspect would have resulted in the same model selection.
the phylogenetic analyses of the data produced an unambiguous picture of most placental mammals relationships.
the ml tree in figure  <dig> was constructed from first and second codon positions  and analyzed under a gtr 8Γ+i model of sequence evolution for the  <dig> , <dig> nt long sequences in tf.
although in this and most other analyses the bat groups with the cow, the branch is shown as unresolved in figure  <dig>  because the placement of the bat on the mammalian tree is indefinite in some analyses.
bayesian analysis of the same data set, on first plus second codon positions or on aa sequence data reconstructed the same topology.
also r/y coded sequences that were analyzed using a two-state  <cit>  ml 4Γ+i model reconstructed the topology shown in figure  <dig>  ml analysis with a non-stationary model in paml  <cit>  confirmed the topology in figure  <dig>  finally, the same tree was also reconstructed after removing sites with more than three different amino acids across all species.
placing the treeshrew on the mammalian tree by sequence analysis was somewhat problematic.
ml reconstructed an alternative position of the treeshrew, outside the primates plus glires clade, when the nt data were not partitioned by codon position.
an extended ml analysis on aa and first and second codon position sequences  could, however, differentiate between alternative placements of the treeshrew on the tree  for the different datasets.
a sister group position to the primates is the best ml option for most analyses.
the treeshrew as sister group to primates plus glires can be rejected by other datasets with statistical confidence.
also the grouping of the treeshrew with the glires can be rejected by ml analysis of aa and nt <dig>  neither a codon analysis of the concatenated dataset or a separate analysis of  <dig> individual genes, a total of  <dig>  nt of sequence data, for which full taxon sampling was available, clarified the position of the treeshrew, because none of the topologies in table 2a received significant support.
alternative phylogenetic positions of the guinea pig among the rodents can clearly be rejected by an sh test on all analyzed data types .
the highly unconventional phylogenetic position of the bat as sister group to the cow and thus probably to all artiodactyla and cetacea is recovered by most analyses and datasets .
the only alternative position is as sister group to the lipotyphla  which cannot be rejected by the sh test on aa or codon sequence data.
it is the favored topology for separate aa analysis of individual genes.
albeit, the support is very limited.
varied results were obtained from ml analyses on data that were partitioned according to the average aa distance.
without following logic the bat grouped either with the cow or the shrew plus hedgehog.
analyses of partitions with over 10% aa distance grouped the bat indistinguishably by ml on any of the neighboring branches.
although most analyses favor a grouping of bat plus cow, the branch leading to the bat, cow and carnivores is shown as unresolved in figure  <dig>
ml analyses of all nt sequences , nt <dig>  codon model analysis and separate analysis of nt sequences favor a split between xenafrotheria and the remaining placentals .
however, analysis of aa sequences ever so slightly favor a split between afrotheria and the remaining placentals , but cannot reject alternative hypothesis.
partitioned ml analysis of codon positions and separate analysis of aa sequences support a split between xenarthra and the remaining placentals .
ml bootstrap support values for the branches that do not receive maximum support are shown in figure  <dig>  thus, bootstrap analysis unambiguously supports the relationships among xenafrotheria, glires and primates.
as expected, nt sequences find only limited support for the treeshrew as sister group to the primates, but even the analysis of aa sequences does not provide maximal bootstrap support for this relationship.
while nt sequences provide the strongest bootstrap support  for placing the cow and bat on a common branch, ml bootstrap analysis of aa sequences data support the bat plus cow grouping only by 31% and instead favors to group the bat with the shrew plus hedgehog on a common branch.
this relationship receives 64% bootstrap support.
divergence times were calculated for the tree topology as shown in figure  <dig>  estimated ages for all divergences are depicted in a chronogram in figure  <dig>  the exact values for different methods and their standard deviations for the r8s dates based on aa and nt sequences are shown in table  <dig>  virtually the same divergence times among boreoplacentalia were calculated for the two alternatively rooted trees .
for some dates the r8s program tend to provide 4– <dig> myr younger time estimates than tf.
however, the dates provided by tf generally conform within a few  <dig>  years to those previously published on phylogenomic analyses  <cit> , and most importantly, only marginal differences are observed for the relative dates.
the standard deviations appear to be unrealistically small, but reflect the amount of data that was used for the analysis.
most divergence time estimates are only negligibly different from analysis of much smaller nuclear or mitogenomic alignments  <cit> .
a notable exception is the age of the lagomorph and rodent divergence, which has previously been estimated 5– <dig> myr younger by phylogenomic analysis  <cit> .
it should be observed that many of the difficult to resolve divergences are splits that are separated by only a few myr.
these involve e.g.
the treeshrew, bat, primate, glires, and the boreoplacentalia and xenafrotheria divergence.
a refers to the cow-bat branching which is shown unresolved in figure  <dig>
bestimated date is on the edge of the calibration interval and thus no standard deviation could be calculated.
the branches c and e do not exist and were not determined .
a refers to the cow-bat branching
discussion
the study of placental mammal relationships has drawn the attention of generations of scientists.
however, the resolution of the mammalian tree has been problematic and many relationships remained unresolved or only poorly supported.
initially phylogenetic studies were made on morphological data from recent and fossil species  <cit> , immunological data  <cit> , sequence data of single genes  <cit> , sequence data from whole mitochondrial genomes  <cit> , concatenation of several nuclear encoded genes  <cit> , hundreds of nuclear genes from expressed sequence tag data  <cit> , and now phylogenomic analyses, which cover 10% of the coding region of the genome.
the use of rare genomic characters such as retroposed elements  <cit>  or indels  <cit>  as cladistic markers has become an important addition to the analysis of sequence data.
the steadily increasing amount of genomic sequence data from placental mammals does not only allow an ever-deeper insight into their evolution, but it also allows development of rigorous data handling procedures.
a central aim in all phylogenomic analyses is selecting orthologous sequences for data analysis.
the lack of complete chromosome maps for most genomes make it currently impractical to utilize synteny information and orthology can only be established when the species phylogeny is known  <cit> .
however, the mammalian phylogeny is often unclear and the very issue of such studies.
therefore, in most phylogenetic analyses orthology is determined by a sequence similarity criterion  <cit>  or choosing from a pre-determined database of orthologous sequences, such as cog  <cit>  and the orthologs matrix project, oma  <cit> .
these databases usually use similarity criteria too and are sometimes combined with manual curation for increasing the confidence in the orthology of the sequences.
in the mammalian genome, the ratio of gene to gene family is about two and the size of about half of the gene families have changed during mammalian evolution  <cit> .
the globin genes are a textbook example of such gene families  <cit> .
the human beta hemoglobin family has six members that have been duplicated at different times during evolution of mammals and they fulfill different functions.
other mammals, such as goat and sheep, which have different history of gene duplications, have as many as  <dig> members of the beta hemoglobin family  <cit> .
thus, there is a considerable risk that the orthology of sequences from gene families cannot be securely established, unless by laborious case-to-case studies for each gene and each species involved.
such an approach is prohibitive for large-scale phylogenetic analyses.
the rigorous approach of the data assembly, recursive blast, in our analysis excludes the vast majority of sequences from multi-gene families, which are the source of paralogous sequences, from the analysis.
already the first filtering step successfully removed sequences from typical gene families.
while this is no formal proof that all genes from gene families were removed, it is conceivable that other, less thoroughly studied, gene families were also eliminated from the analysis during this process.
thus, the risk for including paralogous sequences in the current analysis was greatly reduced by the recursive blast approach.
the data collection and filtering produced a  <dig>  mbp long dataset composed of  <dig> alignments of inferred transcripts from genes for  <dig> mammalian and one avian species with a moderate evolutionary rate.
this alignment is well suited for studying deep phylogenetic divergences in the placental tree.
the compositional bias of the data remained to be a major obstacle in the analysis, which could not be overcome for most sequences by recoding the sequences to r/y.
however, congruency of the reconstructed trees with analyses using a non-stationary model, indicate that compositional bias does not influence the tree topology to a greater extend.
further and less well-defined systematic errors are a general threat to phylogenetic analysis  <cit>  and are the source for strong support for erroneous relationships.
removing sites with more than three different aa, i.e relatively variable sites, or analyzing different data partitions according to evolutionary rate are a recommend cure to the problem.
however, these analyses did not change the topology.
the taxon sampling currently available for genome analysis does not yet allow for investigating if one or several rouge taxa influence the ml reconstruction.
the reconstructed phylogeny is largely congruent with that of previous studies.
most analyses favor a basal divergence of xenafrotheria and boreoplacentalia.
this split is supported by phylogenomic sequence analyses, two retroposon insertions and indels  <cit> .
the revised naming of the boreoeutheria clade is preferable to the previously used term, because of the improved consistency and logic.
recent naming of new clades were largely based on molecular phylogenetics, which by definition exclude the stem group members of a clade.
thus, the term "eutheria" in the formerly used boreoeutheria encompasses also stem group eutherians only known from the fossil record.
in comparison "placentalia" describes the respective clade more accurately, referring only to crown group members, which are accessible for molecular studies  <cit> .
the elephant and tenrec group, as expected, in the afrotheria clade  <cit> .
the boreoplacentalia are split into laurasiatheria and eurchontoglires, the latter with a well-supported glires clade.
within the rodents the guinea pig fell between myomorpha  and sciuromorpha , which is consistent with some previous findings of partially constraint tree analyses  <cit> , but inconsistent with some mitogenomic studies, which find the squirrel  closer to the muroid mouse and rat to the exclusion of the guinea pig  <cit> .
the phylogenetic position of the treeshrew  within the euarchontoglires was problematic even when analyzing nearly  <dig> mnt of protein coding data.
the uncertainty of placing the treeshrews among placental mammals is evident from the numerous positions it had been placed on the tree in previous studies.
the treeshrews have been placed as the first divergence among euarchontaglires  <cit> , as sister group to lagomorpha or rodents  <cit> , sister group to dermoptera  <cit> , sister group to primates  <cit>  or part of the traditional archonta clade outside a , chiroptera) clade  <cit> .
two recent studies on rare genomic changes, found support for a grouping the treeshrew with primates, from four indels  <cit>  and five retroposed elements  <cit>  respectively.
this conforms with our sequence analysis based findings.
the placement of the bat as sister group to the cow in most ml analyses is highly unexpected and has not been proposed in any other analyses.
currently the result can only be regarded as peculiar.
the only alternative place of the bat that cannot be rejected by ml analysis of aa sequences is on a common branch to the lipotyphla .
this grouping receives inconclusive 64% bootstrap support.
mitogenomic and nuclear gene data generally find the bat  as the sister group to cetferungulata  <cit> , which is the currently preferred hypothesis.
the notable exception to the general consensus came from a recent retroposon analysis that found support for a sister group relationship of chiroptera to perissodactyla plus carnivora  <cit> .
if this particular relationship was supported by our sequence data, the bat would have been placed on the carnivore branch, because there is no species from perissodactyla present in our data set.
however, ml analyses can exclude a relationship of the bat and carnivores and most alternative positions with high statistical confidence.
systematic errors can never be formally excluded as a cause for reconstructing inaccurate relationships.
therefore we have investigated the effect of two major sources of systematic error, namely highly variable sites and highly variable sequences  <cit> .
the exclusion or separate analyses of these had no clear consequence on the reconstructed topology.
until genome data from further species become available for analysis, we regard the phylogenetic position of the bat as unsettled by genomic sequence data.
the divergences among the placental mammal orders in this study occurred between ≈ <dig> and ≈ <dig> mya, thus well inside the cretaceous  and agrees with previous molecular based estimates on mitogenomic and phylogenomic data  <cit>  but is in steep contrast with some recent interpretations of the fossil record  <cit> .
many ordinal and higher level divergences occur within a few myr from each other.
tf estimates the xenafrotheria and boreoplacentalia divergence ≈ <dig> mya and the first split within these clades is estimated to ≈ <dig> mya and ≈ <dig> mya.
the r8s estimates are somewhat younger, but present the same picture.
similarly, the euarchontoglires clade diverges at ≈ <dig> mya into glires and euarchonta, and at ≈ <dig> mya into treeshrew  and the primate lineage.
these splits are remarkably close to each other and would prove problematic to resolve even from recent species.
one major problem of resolving short divergence time intervals is that an insufficient amount of phylogenetic information accumulates in these short intervals.
this obstacle can usually be overcome simply by increasing the amount of data, which is the current strategy in phylogenomic analyses.
the other, more severe, problem that cannot as easily be solved is connected to genomic processes around speciation.
before closely related species become completely genetically isolated they can hybridize with each other.
this poses a problem for phylogenetic analysis, because the species history can be obscured by introgression, the incorporation of genes from one species into the gene pool of another.
it has been estimated that such hybridization can occur up to 7– <dig> myr after separation in birds  <cit>  and thus probably also in mammals.
cattle-zebu  hybrids are common and fully fertile despite that the respective lineages have been separated for 1– <dig> myr  <cit> .
there are surprisingly many examples where viable, fertile and even well adapted offspring form by species hybridization in mammals  <cit> .
it has been estimated that 10% of all animals , actually hybridize with at least one other species  <cit> .
many more may have gone through a phase of hybridization that remains undetected, because the species that now are observed to hybridize are usually the evolutionarily youngest.
another process that can obscure phylogenetic analysis is lineage sorting.
lineage sorting is the result of an ancestral polymorphism that survives a speciation event, succeeded by fixation of different alleles in the descendant lineages.
this gives rise to a gene history that is incongruent with the species' history  <cit> .
the extent to which lineage sorting jumbles the genome of closely related species has recently been investigated for the anthropoid genomes, i.e.
the genomes of human, chimpanzee, gorilla, and orangutan.
for about one quarter of the human genome gorilla and not chimpanzee is the closest genetic relative  <cit> .
in  <dig> % of the studied genes the human sequence share closest ancestry even with the orangutan.
lineage sorting is currently the favored hypothesis to explain the mosaic of the primate genomes  <cit> .
however, other natural processes such as hybridization have been used to explain the fact that some parts of the genome are more similar between human and gorilla than between human and chimpanzee  <cit> .
finally, a process with a similar effect on the genome as species hybridization is hybrid speciation  <cit> , which may be a new and radical explanation for the concoction of the mammalian genome.
in this process two previously isolated sister-species hybridize and give offspring to a new and possibly better-adapted species  <cit> .
while such a process is common in plants, there is growing evidence that this also occurs in animals, including vertebrates  <cit> .
which of these processes, lineage sorting or hybridization, has the dominating effect on the mammalian genome cannot be determined yet and needs to be studied in more detail.
in any case, these processes seem to have a sizeable impact on the resolution of short internal branches of the mammalian tree and phylogenetic reconstruction in general.
the studies of the jumbled anthropoid genomes  <cit>  are the first to put some exact values on the amount of gene trees that are incongruent with the species tree.
it must be recognized that there is a non-negligible chance that a fair amount of single loci have a history that is not compatible with the species history.
therefore introgression or linage sorting are problematic especially to super-tree analyses  <cit>  and analyses from rare genomic events  <cit> .
in these approaches of phylogenetic reconstruction many branches are only supported by single or at most a few genes or loci.
incompatibilities in the support from the analysis of retroposon insertion for basal and ordinal mammalian divergences  <cit>  probably document introgression or linage sorting events that confound the genome history.
probably the only solution is the collection of large amounts of information as in phylogenomic analyses in the hope that the massive amount of data evens out events that lead to a genome history that is inconsistent with the species history.
large phylogenomic analyses produce a tree that, at least on average, reflects the species history.
completely or nearly completely sequenced genome data are becoming available for phylogenomic analysis of mammalian relationships at an increasing rate.
this allows ever more detailed reconstruction of their history and their major branches are becoming with time more consistently recovered.
however, some mammalian relationships will require more detailed studies of their history, taking into account that the genome is not a fixed entity but malleable by speciation events.
CONCLUSIONS
phylogenomic analysis of  <dig>  genes  from  <dig> placental mammal species could significantly resolve most relationships and date their basal divergences in the cretaceous at  <dig> –  <dig> mya.
however, the mammalian divergences that occurred in very narrow time windows of 2– <dig> myr remain surprisingly difficult to resolve, even by the huge amount of genomic data.
divergences of lineages that are now considered to be orders, started as speciation events even if they occurred millions of years ago.
recent species that have been separated for about the same amount of time can still hybridize and have fully fertile offspring, leading to introgression, or may be affected by lineage sorting that results in a gene history that is different from the species history.
thus, some deep mammalian divergences that are separated by only a few myr appear to be affected by speciation processes that obscure the phylogenetic reconstruction.

 recently the resolution of the mammalian tree made a quantum leap forward with the analysis of protein coding sequences of whole genome data. the genome data allowed collecting the protein coding sequences from some genes, which is equivalent to million nucleotides, representing ≈10% of all coding sequences. during the last years since the influential review of novacek on the mammalian evolution, many previously uncertain relationships are now becoming consistently resolved by different data sets and analytical approaches. however, lately released mammal genome data have not been investigated by phylogenomic analyses. presently the monophyly of each of the four major clades, euarchontoglires, laurasiatheria, xenarthra, and afrotheria is supported by most analyses based on sequence data and genome level data such as line and ltr elements. the smallest clade, xenarthra, consists of only two orders, cingulata and pilosa. the afrotheria comprises the traditional orders sirenia, proboscidea, hyracoidea, and tubulidentata plus the members from the paraphyletic lipotyphla: the tenrecidae, macroscelidea and chrysochloridea. the afrotheria and xenarthra are grouped in the clade named xenafrotheria, a clade that is supported by phylogenomic analyses and retroposon and indel analysis. other studies, however, find support for a basal xenarthra clade from retroposon data, or a basal afrotheria clade from some sequence data analyses. the two remaining clades, euarchontoglires and laurasiatheria together comprise the species rich taxon boreoplacentalia, which is supported by all current analyses. the new name boreoplacentalia has been suggested for this clade for being more consistent and specific than the previous name boreoeutheria. many branches of the mammalian tree remained difficult to resolve, even with the analysis of several thousand nt of sequence data. the most basal placental mammal divergences were inconclusive even with the analysis of some knt of coding sequence data. the analysis of a ten-fold larger dataset, mnt of protein coding sequence data supported the grouping of xenarthra with the afrotheria in the higher order clade xenafrotheria. jackknife analysis showed that at least knt sites of protein coding data were needed to significantly resolve this relationship. the notoriously difficult to resolve split between xenafrotheria and boreoplacentalia took place at ≈ million years ago. however, the two clades diverged into xenarthra and afrotheria, and laurasiatheria and euarchontoglires respectively only 3– million years later. the resolution of these narrow temporal occurrences required more data than available to the study of nikolaev et al., however the placement of xenarthra and afrotheria remains challenged. in comparison, more distant branching events such as that of monotreme, marsupial and placental mammals can be significantly resolved with as little as aa sites of nuclear coding genes. molecular dating on phylogenomic data estimated that 30– myr separated the splits among the three mammalian infra-classes. this leaves enough time to accumulate phylogenetically relevant sequence differences and to reduce the effects of lineage sorting and prevents introgression. while the monophyly of the clade laurasiatheria is well established, phylogenomic analysis did not yet resolve the relationships within this clade. the internal branches of dog, cow, bat, and shrew received only 86% and 93% bootstrap support, respectively, in a recent phylogenomic analyses. this non-significant support may be due to the "insufficient" amount of knt of coding sequence that was available at the time of the analysis. however, the splits may also have occurred at a very narrow time window of a few millions years or less, and more data than expected may be required to resolve such tight divergences. it is noteworthy that the few million years between mammal ordinal divergences are in the order of the average duration of mammalian species and their typical speciation times. the fact that the splits between the extant orders occurred several tens of million years ago does not diminish the problem that these splits represent divergences among species or even lineages that were at the brink of becoming genetically separated species. therefore, speciation adds an additional level of complication for resolving such narrow divergences. not only is the time for such divergences relatively short to accumulate sufficient sequence differences that are needed to resolve these splits, but lineage sorting, species hybridization and hybrid speciation can make the resolution of significant parts of the mammalian tree all but impossible. the continuously increasing amount of sequence data from protein coding genes now makes it possible to provide a phylogenomic analysis of placental mammals for species on the basis of several million characters long alignments. the quantity of new genome data also allows using a new and more stringent assembly of the data, which aims at avoiding data from gene families and thus paralogous sequences in the alignment. we investigate if the resolution of the basal split among placental mammals prevails with the increased taxon sampling and if other difficult to resolve divergences such as those within the laurasiatheria and euarchontoglires can be resolved with statistical confidence with an increased amount of phylogenomic data and new analytical approaches.  the first step in the data collection, searching for human protein-coding genes that do not have obviously closely related members of a gene family resulted in a set of human sequences, from the original unique genes. these sequences were used for identifying related sequences from the remaining species. after selecting against genes from multi-gene families in any of the species and using only alignments in which more than species are represented, multiple sequence alignments remained. from this data set sequence alignments that were shorter than bp or observed aa distances larger than 30% for any species pair were removed. this left a data set of multiple sequence alignments. the total size of this alignment was, nucleotides for species. the average length of the individual sequences is ± nt and the average p distance between human and platypus is ± %, with a maximum of 30%, which is expected from the filtering constraints. the rigorous filtering by recursive blast search removed sequences that are part of gene families notably reduced the amount of data. the filtering was rather effective, even though there is no formal proof for its efficiency yet. already after the first filtering step sequences of the hemoglobin, myosin, olfactory receptor, hox gene and other sequences that are part of known gene families were eliminated from the data set. it is conceivable that by this approach data from other, less known gene families were also excluded from the analysis. the recursive blast search with a cutoff value of 10- identified genes in gene families that had more than 75% sequence similarity and eliminated these. it can be argued that more distant sequences will not be found when searching for such sequences in other species with the same cutoff value of 10- when inspected by eye the base and amino acid frequencies appear to be very similar among the species. however, due to the large size of the dataset the chi-square test may become overly strict. compositional homogeneity was rejected for many species both for nucleotide and amino acid data. the compositional heterogeneity remained even when the data were r/y-coded. a species for which compositional homogeneity could not be rejected with 5% significance, and are thus assumed to be homogeneous. the evolutionary model and the number of rate heterogeneity classes were estimated from previous phylogenomic analyses. exact tests such as implemented in modeltest could not be performed due to computational constraints and the inability of most programs to analyze a data set of this magnitude. a manual analysis would have been prohibitively time consuming to perform and we suspect would have resulted in the same model selection. the phylogenetic analyses of the data produced an unambiguous picture of most placental mammals relationships. the ml tree in figure was constructed from first and second codon positions and analyzed under a gtr 8Γ+i model of sequence evolution for the, nt long sequences in tf. although in this and most other analyses the bat groups with the cow, the branch is shown as unresolved in figure because the placement of the bat on the mammalian tree is indefinite in some analyses. bayesian analysis of the same data set, on first plus second codon positions or on aa sequence data reconstructed the same topology. also r/y coded sequences that were analyzed using a two-state ml 4Γ+i model reconstructed the topology shown in figure ml analysis with a non-stationary model in paml confirmed the topology in figure finally, the same tree was also reconstructed after removing sites with more than three different amino acids across all species. placing the treeshrew on the mammalian tree by sequence analysis was somewhat problematic. ml reconstructed an alternative position of the treeshrew, outside the primates plus glires clade, when the nt data were not partitioned by codon position. an extended ml analysis on aa and first and second codon position sequences could, however, differentiate between alternative placements of the treeshrew on the tree for the different datasets. a sister group position to the primates is the best ml option for most analyses. the treeshrew as sister group to primates plus glires can be rejected by other datasets with statistical confidence. also the grouping of the treeshrew with the glires can be rejected by ml analysis of aa and nt neither a codon analysis of the concatenated dataset or a separate analysis of individual genes, a total of nt of sequence data, for which full taxon sampling was available, clarified the position of the treeshrew, because none of the topologies in table 2a received significant support. alternative phylogenetic positions of the guinea pig among the rodents can clearly be rejected by an sh test on all analyzed data types. the highly unconventional phylogenetic position of the bat as sister group to the cow and thus probably to all artiodactyla and cetacea is recovered by most analyses and datasets. the only alternative position is as sister group to the lipotyphla which cannot be rejected by the sh test on aa or codon sequence data. it is the favored topology for separate aa analysis of individual genes. albeit, the support is very limited. varied results were obtained from ml analyses on data that were partitioned according to the average aa distance. without following logic the bat grouped either with the cow or the shrew plus hedgehog. analyses of partitions with over 10% aa distance grouped the bat indistinguishably by ml on any of the neighboring branches. although most analyses favor a grouping of bat plus cow, the branch leading to the bat, cow and carnivores is shown as unresolved in figure ml analyses of all nt sequences, nt codon model analysis and separate analysis of nt sequences favor a split between xenafrotheria and the remaining placentals. however, analysis of aa sequences ever so slightly favor a split between afrotheria and the remaining placentals, but cannot reject alternative hypothesis. partitioned ml analysis of codon positions and separate analysis of aa sequences support a split between xenarthra and the remaining placentals. ml bootstrap support values for the branches that do not receive maximum support are shown in figure thus, bootstrap analysis unambiguously supports the relationships among xenafrotheria, glires and primates. as expected, nt sequences find only limited support for the treeshrew as sister group to the primates, but even the analysis of aa sequences does not provide maximal bootstrap support for this relationship. while nt sequences provide the strongest bootstrap support for placing the cow and bat on a common branch, ml bootstrap analysis of aa sequences data support the bat plus cow grouping only by 31% and instead favors to group the bat with the shrew plus hedgehog on a common branch. this relationship receives 64% bootstrap support. divergence times were calculated for the tree topology as shown in figure estimated ages for all divergences are depicted in a chronogram in figure the exact values for different methods and their standard deviations for the r8s dates based on aa and nt sequences are shown in table virtually the same divergence times among boreoplacentalia were calculated for the two alternatively rooted trees. for some dates the r8s program tend to provide 4– myr younger time estimates than tf. however, the dates provided by tf generally conform within a few years to those previously published on phylogenomic analyses, and most importantly, only marginal differences are observed for the relative dates. the standard deviations appear to be unrealistically small, but reflect the amount of data that was used for the analysis. most divergence time estimates are only negligibly different from analysis of much smaller nuclear or mitogenomic alignments. a notable exception is the age of the lagomorph and rodent divergence, which has previously been estimated 5– myr younger by phylogenomic analysis. it should be observed that many of the difficult to resolve divergences are splits that are separated by only a few myr. these involve e.g. the treeshrew, bat, primate, glires, and the boreoplacentalia and xenafrotheria divergence. a refers to the cow-bat branching which is shown unresolved in figure bestimated date is on the edge of the calibration interval and thus no standard deviation could be calculated. the branches c and e do not exist and were not determined. a refers to the cow-bat branching discussion the study of placental mammal relationships has drawn the attention of generations of scientists. however, the resolution of the mammalian tree has been problematic and many relationships remained unresolved or only poorly supported. initially phylogenetic studies were made on morphological data from recent and fossil species, immunological data, sequence data of single genes, sequence data from whole mitochondrial genomes, concatenation of several nuclear encoded genes, hundreds of nuclear genes from expressed sequence tag data, and now phylogenomic analyses, which cover 10% of the coding region of the genome. the use of rare genomic characters such as retroposed elements or indels as cladistic markers has become an important addition to the analysis of sequence data. the steadily increasing amount of genomic sequence data from placental mammals does not only allow an ever-deeper insight into their evolution, but it also allows development of rigorous data handling procedures. a central aim in all phylogenomic analyses is selecting orthologous sequences for data analysis. the lack of complete chromosome maps for most genomes make it currently impractical to utilize synteny information and orthology can only be established when the species phylogeny is known. however, the mammalian phylogeny is often unclear and the very issue of such studies. therefore, in most phylogenetic analyses orthology is determined by a sequence similarity criterion or choosing from a pre-determined database of orthologous sequences, such as cog and the orthologs matrix project, oma. these databases usually use similarity criteria too and are sometimes combined with manual curation for increasing the confidence in the orthology of the sequences. in the mammalian genome, the ratio of gene to gene family is about two and the size of about half of the gene families have changed during mammalian evolution. the globin genes are a textbook example of such gene families. the human beta hemoglobin family has six members that have been duplicated at different times during evolution of mammals and they fulfill different functions. other mammals, such as goat and sheep, which have different history of gene duplications, have as many as members of the beta hemoglobin family. thus, there is a considerable risk that the orthology of sequences from gene families cannot be securely established, unless by laborious case-to-case studies for each gene and each species involved. such an approach is prohibitive for large-scale phylogenetic analyses. the rigorous approach of the data assembly, recursive blast, in our analysis excludes the vast majority of sequences from multi-gene families, which are the source of paralogous sequences, from the analysis. already the first filtering step successfully removed sequences from typical gene families. while this is no formal proof that all genes from gene families were removed, it is conceivable that other, less thoroughly studied, gene families were also eliminated from the analysis during this process. thus, the risk for including paralogous sequences in the current analysis was greatly reduced by the recursive blast approach. the data collection and filtering produced a mbp long dataset composed of alignments of inferred transcripts from genes for mammalian and one avian species with a moderate evolutionary rate. this alignment is well suited for studying deep phylogenetic divergences in the placental tree. the compositional bias of the data remained to be a major obstacle in the analysis, which could not be overcome for most sequences by recoding the sequences to r/y. however, congruency of the reconstructed trees with analyses using a non-stationary model, indicate that compositional bias does not influence the tree topology to a greater extend. further and less well-defined systematic errors are a general threat to phylogenetic analysis and are the source for strong support for erroneous relationships. removing sites with more than three different aa, i.e relatively variable sites, or analyzing different data partitions according to evolutionary rate are a recommend cure to the problem. however, these analyses did not change the topology. the taxon sampling currently available for genome analysis does not yet allow for investigating if one or several rouge taxa influence the ml reconstruction. the reconstructed phylogeny is largely congruent with that of previous studies. most analyses favor a basal divergence of xenafrotheria and boreoplacentalia. this split is supported by phylogenomic sequence analyses, two retroposon insertions and indels. the revised naming of the boreoeutheria clade is preferable to the previously used term, because of the improved consistency and logic. recent naming of new clades were largely based on molecular phylogenetics, which by definition exclude the stem group members of a clade. thus, the term "eutheria" in the formerly used boreoeutheria encompasses also stem group eutherians only known from the fossil record. in comparison "placentalia" describes the respective clade more accurately, referring only to crown group members, which are accessible for molecular studies. the elephant and tenrec group, as expected, in the afrotheria clade. the boreoplacentalia are split into laurasiatheria and eurchontoglires, the latter with a well-supported glires clade. within the rodents the guinea pig fell between myomorpha and sciuromorpha, which is consistent with some previous findings of partially constraint tree analyses, but inconsistent with some mitogenomic studies, which find the squirrel closer to the muroid mouse and rat to the exclusion of the guinea pig. the phylogenetic position of the treeshrew within the euarchontoglires was problematic even when analyzing nearly mnt of protein coding data. the uncertainty of placing the treeshrews among placental mammals is evident from the numerous positions it had been placed on the tree in previous studies. the treeshrews have been placed as the first divergence among euarchontaglires, as sister group to lagomorpha or rodents, sister group to dermoptera, sister group to primates or part of the traditional archonta clade outside a, chiroptera) clade. two recent studies on rare genomic changes, found support for a grouping the treeshrew with primates, from four indels and five retroposed elements respectively. this conforms with our sequence analysis based findings. the placement of the bat as sister group to the cow in most ml analyses is highly unexpected and has not been proposed in any other analyses. currently the result can only be regarded as peculiar. the only alternative place of the bat that cannot be rejected by ml analysis of aa sequences is on a common branch to the lipotyphla. this grouping receives inconclusive 64% bootstrap support. mitogenomic and nuclear gene data generally find the bat as the sister group to cetferungulata, which is the currently preferred hypothesis. the notable exception to the general consensus came from a recent retroposon analysis that found support for a sister group relationship of chiroptera to perissodactyla plus carnivora. if this particular relationship was supported by our sequence data, the bat would have been placed on the carnivore branch, because there is no species from perissodactyla present in our data set. however, ml analyses can exclude a relationship of the bat and carnivores and most alternative positions with high statistical confidence. systematic errors can never be formally excluded as a cause for reconstructing inaccurate relationships. therefore we have investigated the effect of two major sources of systematic error, namely highly variable sites and highly variable sequences. the exclusion or separate analyses of these had no clear consequence on the reconstructed topology. until genome data from further species become available for analysis, we regard the phylogenetic position of the bat as unsettled by genomic sequence data. the divergences among the placental mammal orders in this study occurred between ≈ and ≈ mya, thus well inside the cretaceous and agrees with previous molecular based estimates on mitogenomic and phylogenomic data but is in steep contrast with some recent interpretations of the fossil record. many ordinal and higher level divergences occur within a few myr from each other. tf estimates the xenafrotheria and boreoplacentalia divergence ≈ mya and the first split within these clades is estimated to ≈ mya and ≈ mya. the r8s estimates are somewhat younger, but present the same picture. similarly, the euarchontoglires clade diverges at ≈ mya into glires and euarchonta, and at ≈ mya into treeshrew and the primate lineage. these splits are remarkably close to each other and would prove problematic to resolve even from recent species. one major problem of resolving short divergence time intervals is that an insufficient amount of phylogenetic information accumulates in these short intervals. this obstacle can usually be overcome simply by increasing the amount of data, which is the current strategy in phylogenomic analyses. the other, more severe, problem that cannot as easily be solved is connected to genomic processes around speciation. before closely related species become completely genetically isolated they can hybridize with each other. this poses a problem for phylogenetic analysis, because the species history can be obscured by introgression, the incorporation of genes from one species into the gene pool of another. it has been estimated that such hybridization can occur up to 7– myr after separation in birds and thus probably also in mammals. cattle-zebu hybrids are common and fully fertile despite that the respective lineages have been separated for 1– myr. there are surprisingly many examples where viable, fertile and even well adapted offspring form by species hybridization in mammals. it has been estimated that 10% of all animals, actually hybridize with at least one other species. many more may have gone through a phase of hybridization that remains undetected, because the species that now are observed to hybridize are usually the evolutionarily youngest. another process that can obscure phylogenetic analysis is lineage sorting. lineage sorting is the result of an ancestral polymorphism that survives a speciation event, succeeded by fixation of different alleles in the descendant lineages. this gives rise to a gene history that is incongruent with the species' history. the extent to which lineage sorting jumbles the genome of closely related species has recently been investigated for the anthropoid genomes, i.e. the genomes of human, chimpanzee, gorilla, and orangutan. for about one quarter of the human genome gorilla and not chimpanzee is the closest genetic relative. in % of the studied genes the human sequence share closest ancestry even with the orangutan. lineage sorting is currently the favored hypothesis to explain the mosaic of the primate genomes. however, other natural processes such as hybridization have been used to explain the fact that some parts of the genome are more similar between human and gorilla than between human and chimpanzee. finally, a process with a similar effect on the genome as species hybridization is hybrid speciation, which may be a new and radical explanation for the concoction of the mammalian genome. in this process two previously isolated sister-species hybridize and give offspring to a new and possibly better-adapted species. while such a process is common in plants, there is growing evidence that this also occurs in animals, including vertebrates. which of these processes, lineage sorting or hybridization, has the dominating effect on the mammalian genome cannot be determined yet and needs to be studied in more detail. in any case, these processes seem to have a sizeable impact on the resolution of short internal branches of the mammalian tree and phylogenetic reconstruction in general. the studies of the jumbled anthropoid genomes are the first to put some exact values on the amount of gene trees that are incongruent with the species tree. it must be recognized that there is a non-negligible chance that a fair amount of single loci have a history that is not compatible with the species history. therefore introgression or linage sorting are problematic especially to super-tree analyses and analyses from rare genomic events. in these approaches of phylogenetic reconstruction many branches are only supported by single or at most a few genes or loci. incompatibilities in the support from the analysis of retroposon insertion for basal and ordinal mammalian divergences probably document introgression or linage sorting events that confound the genome history. probably the only solution is the collection of large amounts of information as in phylogenomic analyses in the hope that the massive amount of data evens out events that lead to a genome history that is inconsistent with the species history. large phylogenomic analyses produce a tree that, at least on average, reflects the species history. completely or nearly completely sequenced genome data are becoming available for phylogenomic analysis of mammalian relationships at an increasing rate. this allows ever more detailed reconstruction of their history and their major branches are becoming with time more consistently recovered. however, some mammalian relationships will require more detailed studies of their history, taking into account that the genome is not a fixed entity but malleable by speciation events.  phylogenomic analysis of genes from placental mammal species could significantly resolve most relationships and date their basal divergences in the cretaceous at – mya. however, the mammalian divergences that occurred in very narrow time windows of 2– myr remain surprisingly difficult to resolve, even by the huge amount of genomic data. divergences of lineages that are now considered to be orders, started as speciation events even if they occurred millions of years ago. recent species that have been separated for about the same amount of time can still hybridize and have fully fertile offspring, leading to introgression, or may be affected by lineage sorting that results in a gene history that is different from the species history. thus, some deep mammalian divergences that are separated by only a few myr appear to be affected by speciation processes that obscure the phylogenetic reconstruction. 
the rigorous data handling, recursive blast, successfully removed the sequences from gene families, including those from well-known families hemoglobin, olfactory, myosin and hox genes, thus avoiding alignment of possibly paralogous sequences. a number of the deeper divergences in the placental mammal tree are still inconclusively resolved despite extensive phylogenomic analyses. relationships with narrow divergence time intervals received unexpectedly limited support even from the phylogenomic analyses. therefore these and other difficult splits were examined by phylogenomic analysis of the hitherto largest sequence alignment. the current phylogenomic analysis of genes from a total of species yielded statistically significant support for most relationships. phylogenomic analysis of divergence times dated the basal placental mammal splits at 95– million years ago. the increasingly complete genome data of placental mammals also allowed developing a novel and stringent data search method. these processes obscure phylogenetic analysis, making some parts of the tree difficult to resolve even with genome data. while some major clades were confirmed using genomic sequence data, the placement of the treeshrew, bat and the relationship between boreoeutheria, xenarthra and afrotheria remained problematic to resolve despite the size of the alignment. the narrow temporal window within which some placental divergences took place suggests that inconsistencies and limited resolution of the mammalian tree may have their natural explanation in speciation processes such as lineage sorting, introgression from species hybridization or hybrid speciation.
 a number of the deeper divergences in the placental mammal tree are still inconclusively resolved despite extensive phylogenomic analyses. a recent analysis of kbp of protein coding sequences yielded only limited support for the relationships among laurasiatheria, probably because the divergences occurred only within a few million years from each other. it is generally expected that increasing the amount of data and improving the taxon sampling enhance the resolution of narrow divergences. therefore these and other difficult splits were examined by phylogenomic analysis of the hitherto largest sequence alignment. the increasingly complete genome data of placental mammals also allowed developing a novel and stringent data search method.  the rigorous data handling, recursive blast, successfully removed the sequences from gene families, including those from well-known families hemoglobin, olfactory, myosin and hox genes, thus avoiding alignment of possibly paralogous sequences. the current phylogenomic analysis of genes from a total of species yielded statistically significant support for most relationships. while some major clades were confirmed using genomic sequence data, the placement of the treeshrew, bat and the relationship between boreoeutheria, xenarthra and afrotheria remained problematic to resolve despite the size of the alignment. phylogenomic analysis of divergence times dated the basal placental mammal splits at 95– million years ago. many of the following divergences occurred only a few million years later. relationships with narrow divergence time intervals received unexpectedly limited support even from the phylogenomic analyses.  the narrow temporal window within which some placental divergences took place suggests that inconsistencies and limited resolution of the mammalian tree may have their natural explanation in speciation processes such as lineage sorting, introgression from species hybridization or hybrid speciation. these processes obscure phylogenetic analysis, making some parts of the tree difficult to resolve even with genome data. 