Hervé Philippe
email: herve.philippe@umontreal.ca

Frédéric Delsuc

Henner Brinkmann

Nicolas Lartillot

PHYLOGENOMICS

Keywords: inconsistency, molecular phylogeny, systematic bias, taxon sampling, tree of life

The continuous flow of genomic data is creating unprecedented opportunities for the reconstruction of molecular phylogenies. Access to wholegenome data means that phylogenetic analysis can now be performed at different genomic levels, such as primary sequences and gene order, allowing for reciprocal corroboration of the results. We critically review the different kinds of phylogenomic methods currently available, paying particular attention to method reliability. Our emphasis is on methods for the analysis of primary sequences because these are the most advanced. We discuss the important issue of statistical inconsistency and show how failing to fully capture the process of sequence evolution in the underlying models leads to tree reconstruction artifacts. We suggest strategies for detecting and potentially overcoming these problems. These strategies involve the development of better models, the use of an improved taxon sampling and the exclusion of phylogenetically misleading data.

INTRODUCTION

The newly arising discipline of phylogenomics owes its existence to the revolutionizing progress in DNA sequencing technology. The number of complete genome sequences is already high and increases at an ever-accelerating pace.

The newly coined term phylogenomics [START_REF] Eisen | Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis[END_REF][START_REF] O'brien | Phylogenomics. Ancestral primate viewed[END_REF] comprises several areas of research at the interplay between molecular biology and evolution. The main issues are: (1) using molecular data to infer species relationships, and (2) using information on species evolutionary history to gain insights into the mechanisms of molecular evolution. The majority of publications on phylogenomics deal with the second aspect (see [START_REF] Sjolander | Phylogenomic inference of protein molecular function: advances and challenges[END_REF] for review). However, our concern here is the use of data at the genomic scale to reconstruct the phylogeny of organisms.

A novel and interesting aspect of phylogenomics lies in the possibility of using molecular information above the primary sequence level. In particular, trees can be inferred from whole genome features such as gene content [START_REF] Fitz-Gibbon | Whole genome-based phylogenetic analysis of free-living microorganisms[END_REF][START_REF] Snel | Genome phylogeny based on gene content[END_REF][START_REF] Tekaia | The genomic tree as revealed from whole proteome comparisons[END_REF], gene-order [START_REF] Korbel | SHOT: a web server for the construction of genome phylogenies[END_REF][START_REF] Sankoff | Gene order comparisons for phylogenetic inference: evolution of the mitochondrial genome[END_REF], intron positions [START_REF] Roy | Resolution of a deep animal divergence by the pattern of intron conservation[END_REF], or protein domain structure [START_REF] Lin | Whole-genome trees based on the occurrence of folds and orthologs: implications for comparing genomes on different levels[END_REF][START_REF] Yang | Phylogeny determined by protein domain content[END_REF]. A potential advantage of these methods is that the complexity of some of these characters (e.g. gene order) renders the character-state space very large, reducing the risk of homoplasy by convergence and reversal, thus rendering the inferred phylogenies more reliable. However, these integrated approaches imply the use of a reduced number of characters relative to primary sequence based approaches. There are about 300 times fewer genes than amino acid positions (assuming a mean protein length of about 300 amino acids), thus increasing the risk of stochastic error.

In fact, stochastic or sampling error constitutes one of the major limitations of standard phylogenetics based on single genes. Because the number of positions of a single gene is small, random noise influences the inference of numerous nodes, leading generally to poorly resolved phylogenetic trees. The idea of using large amounts of genomic data as a way to address this problem is not new. For example, to resolve the original tritomy between chimpanzee, human, and gorilla, about ten kilobases were sequenced as early as the late 1980's [START_REF] Miyamoto | Molecular systematics of higher primates: genealogical relations and classification[END_REF]. However, technical and financial limitations have often confined molecular systematics to the use of a few markers (e.g. ribosomal RNA (rRNA)), for which a large diversity of organisms have been sequenced.

Numerous important phylogenetic questions remained unsolved and great hope was placed into the wealth of genomic data soon to be available.

Sampling (or stochastic) error should vanish as the number of genes added to the analysis gets large enough. In practice, this means that statistical support (e.g. bootstrap support) will eventually rise to 100% as more genes are considered. The use of tens of thousands, or millions, of aligned positions that provide a great deal of phylogenetic information should ultimately lead to fully resolved trees. Indeed, several empirical studies confirmed this premise [START_REF] Bapteste | The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba[END_REF][START_REF] Madsen | Parallel adaptive radiations in two major clades of placental mammals[END_REF][START_REF] Murphy | Molecular phylogenetics and the origins of placental mammals[END_REF][START_REF] Qiu | The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes[END_REF][START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF], Soltis et al 1999). This increased resolution leads to the optimistic view that phylogenomics would "end the incongruencies" observed in single gene phylogenies [START_REF] Gee | Evolution: ending incongruence[END_REF]. However, whether the resulting, highly supported, phylogenetic trees are the true ones is not certain.

Systematic error and consistency

The most important challenge of phylogenomics is to verify that tree reconstruction methods are consistent, i.e. converge towards the correct answer as more and more characters are considered [START_REF] Felsenstein | Cases in which parsimony or compatibility methods will be positively misleading[END_REF][START_REF] Felsenstein | Phylogenies from molecular sequences: inference and reliability[END_REF]. In principle, at least in a probabilistic framework, a lack of consistency can always be traced back to some violation of model assumptions by the data analyzed. Note that methods that are not explicitly model based, such as Maximum Parsimony (MP), are equivalent to a statistical analysis under an implicit model [START_REF] Steel | Parsimony, likelihood, and the role of models in molecular phylogenetics[END_REF]. The best understood causes of method inconsistency stem from models that do not properly account for: (1) variable evolutionary rates, leading to the long branch attraction (LBA) artifact [START_REF] Felsenstein | Cases in which parsimony or compatibility methods will be positively misleading[END_REF], (2) heterogeneous nucleotide/amino acid compositions, resulting in the artificial grouping of species that share the same bias [START_REF] Lockhart | Recovering evolutionary trees under a more realistic model of sequence evolution[END_REF], and (3) heterotachy, i.e. shift of position-specific evolutionary rates [START_REF] Kolaczkowski | Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous[END_REF][START_REF] Lockhart | Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis[END_REF], Philippe & Germot 2000).

These systematic biases could be interpreted respectively as rate signal, compositional signal and heterotachous signal, which we will collectively refer to as non-phylogenetic signals [START_REF] Ho | Tracing the decay of the historical signal in biological sequence data[END_REF]. In other words, nonphylogenetic signals are due to substitutions that occurred along the true phylogeny, but that are misinterpreted by tree reconstruction methods as supporting an alternative topology.

Compared to single gene studies, inconsistency is more pronounced in phylogenomic analyses. For example, a reanalysis of the large dataset of Rokas et al. [START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF] demonstrates that, depending on the method used, mutually incongruent, yet 100% supported, trees could be obtained [START_REF] Phillips | Genome-scale phylogeny and the detection of systematic biases[END_REF]. It is well-accepted that the analysis of phylogenomic datasets will necessarily increase the resolution of the trees through the "increase of the signal-to-noise ratio" [START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF]. Indeed, the signal-to-random-noise ratio increases but the phylogenetic-to-non-phylogenetic signal ratio remains constant whatever the number of genes considered (assuming that the gene sampling is not biased). In this review, we will focus on best practices for enhancing the phylogenetic signal in genomic data, while reducing the impact of erroneous signals, in order to obtain accurate and robust trees.

ASSEMBLY OF PHYLOGENOMIC DATASETS

The reliability of a phylogenetic tree depends on the quality of the data and the accuracy of the reconstruction method. In 1988, Felsenstein noted that "molecular evolutionists who use methods for inferring phylogenies do not engage in much discussion of the properties of the methods they use since they focus on the difficult task of collecting the data" p. 523 [START_REF] Felsenstein | Phylogenies from molecular sequences: inference and reliability[END_REF]. Almost 20 years later, molecular systematists still spend much of their time assembling larger and larger datasets, and the crucial discussion about inference methods remains neglected. In phylogenomics, the reliability of the inference is often simply justified by the large number of characters used. Nevertheless, the problem of data acquisition deserves further discussion, as it can heavily compromise the subsequent analysis.

Importance of a rich taxon sampling

A long-standing debate in phylogenetics concerns the relative importance of improving taxon versus gene sampling [START_REF] Graybeal | Is it better to add taxa or characters to a difficult phylogenetic problem?[END_REF][START_REF] Hillis | Is sparse taxon sampling a problem for phylogenetic inference?[END_REF][START_REF] Rosenberg | Taxon sampling, bioinformatics, and phylogenomics[END_REF]. In the genomic age, gene sampling would seem not to be an issue. However the limited resources devoted to systematics often prevent sequencing the genomes of all relevant species. Two strategies can be used, depending on the importance accorded to taxon sampling: (1) gathering complete genome sequences from a few key organisms, or (2) gathering incomplete, yet large, genome sequences from a great diversity of organisms.

The first approach is supported by some computer simulation studies [START_REF] Rosenberg | Taxon sampling, bioinformatics, and phylogenomics[END_REF] and is the most frequently used in phylogenomic analyzes [START_REF] Blair | The evolutionary position of nematodes[END_REF][START_REF] Goremykin | The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm[END_REF][START_REF] Misawa | Revisiting the Glires concept--phylogenetic analysis of nuclear sequences[END_REF][START_REF] Philip | The Opisthokonta and the Ecdysozoa may not be clades: Stronger support for the grouping of plant and animal than for animal and fungi and stronger support for the Coelomata than Ecdysozoa[END_REF][START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF][START_REF] Wolf | Coelomata and not ecdysozoa: evidence from genome-wide phylogenetic analysis[END_REF]. However, the design of computer simulations and the interpretation of their results, make it difficult to draw firm conclusions from this approach [START_REF] Hillis | Is sparse taxon sampling a problem for phylogenetic inference?[END_REF][START_REF] Rosenberg | Taxon sampling, bioinformatics, and phylogenomics[END_REF].

Empirical evidence seems nevertheless to argue against the taxon-poor approach, as illustrated by the phylogeny of metazoans. Two new clades (Ecdysozoa, the moulting animals, including among others arthropods and nematodes, and Lophotrochozoa including among others annelids, molluscs and platyhelminthes) were proposed from rRNA analyses [START_REF] Aguinaldo | Evidence for a clade of nematodes, arthropods and other moulting animals[END_REF].

However, several phylogenomic studies strongly supported the paraphyly of Ecdysozoa, when considering a few model organisms and using a distant outgroup [START_REF] Blair | The evolutionary position of nematodes[END_REF][START_REF] Dopazo | Phylogenomics and the number of characters required for obtaining an accurate phylogeny of eukaryote model species[END_REF][START_REF] Philip | The Opisthokonta and the Ecdysozoa may not be clades: Stronger support for the grouping of plant and animal than for animal and fungi and stronger support for the Coelomata than Ecdysozoa[END_REF][START_REF] Wolf | Coelomata and not ecdysozoa: evidence from genome-wide phylogenetic analysis[END_REF].

In contrast, when 49 species and 71 genes (20,705 positions) are used [START_REF] Philippe | Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia[END_REF], the monophyly of Ecdysozoa and Lophotrochozoa is recovered with strong support (Figure 1a). However, the removal of close outgroups leads to drastic changes (Figure 1b): the fast evolving lineages emerge paraphyletically at the base of the tree. Such an asymmetrical tree-shape is expected to result from LBA artifact, when the outgroup is distantly related [START_REF] Philippe | How good are deep phylogenetic trees?[END_REF], as fungi are. Contrary to the situation with a few species discussed above, the statistical support for these incorrect placements is weak (bootstrap values between 32 and 63), demonstrating that an increased taxon sampling (from 4-10 to 45) has reduced, but not eliminated, the impact of LBA. In fact, the low bootstrap support (Figure 1b) demonstrates that non-phylogenetic signal becomes equivalent to phylogenetic signal when species sampling is impoverished.

A rich-taxon sampling is not the panacea however. First, computation time increases rapidly with the number of species, rendering exhaustive searches impossible with more than 20 species and most heuristic searches with probabilistic methods intractable with more than ~200 species. Second, adding taxa can sometimes degrade the phylogenetic inference [START_REF] Kim | General inconsistency conditions for maximum parsimony: effects of branch lengths and increasing numbers of taxa[END_REF]. Third, the number of extant species can be naturally sparse, forever preventing the assembly of a rich and balanced taxon sample. For example, Amborella is proposed to constitute the first, or one of the first, emerging angiosperm lineages [START_REF] Qiu | The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes[END_REF], but is the only extant representative of an ancient group. Even if the heated debate about its placement [START_REF] Goremykin | The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm[END_REF][START_REF] Soltis | Genome-scale data, angiosperm relationships, and "ending incongruence": a cautionary tale in phylogenetics[END_REF] could be solved by improved taxon sampling, the assumed basal position of Amborella might prove difficult to attest in the absence of closely related extant taxa [START_REF] Stefanovic | Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?[END_REF]. In conclusion, an adequate taxon sampling, as balanced as possible, is important to increase the accuracy of phylogenomic trees.

Missing data

Although genome sequencing has become ever easier, it seems unlikely that complete genomes will be soon available for a rich diversity of organisms. In addition, a bias in favor of sequencing small genomes leads to potential problems. Since small genomes are generally derived from larger genomes, whole-genome features, such as gene content or gene order, will evolve much faster, rendering tree reconstruction susceptible to artifacts such as LBA and compositional bias [START_REF] Copley | Systematic searches for molecular synapomorphies in model metazoan genomes give some support for Ecdysozoa after accounting for the idiosyncrasies of Caenorhabditis elegans[END_REF][START_REF] House | Using homolog groups to create a wholegenomic tree of free-living organisms: an update[END_REF][START_REF] Korbel | SHOT: a web server for the construction of genome phylogenies[END_REF][START_REF] Lake | Deriving the genomic tree of life in the presence of horizontal gene transfer: conditioned reconstruction[END_REF][START_REF] Wolf | Genome trees constructed using five different approaches suggest new major bacterial clades[END_REF]. Moreover, genome reduction is often associated with an accelerated rate of protein evolution (Brinkmann et al accepted, Dufresne et al 2005), or extremely biased nucleotide compositions [START_REF] Herbeck | Non-homogeneous model of sequence evolution indicates independent origins of primary endosymbionts within the Enterobacteriales ({gamma}-Proteobacteria)[END_REF]. The sampling of a fraction of the genome from species with huge genomes is therefore a necessity to represent some key taxa and/or to include less biased representatives.

Two low cost approaches can be used: (1) the selection of a limited set of genes potentially useful for the phylogenetic question of interest, followed by their targeted PCR amplification and sequencing [START_REF] Murphy | Molecular phylogenetics and the origins of placental mammals[END_REF]; (2) the sequencing of thousands of Expressed Sequence Tags (ESTs), which generally provides hundreds of relevant genes [START_REF] Bapteste | The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba[END_REF]. The first method is more adapted to phylogeny at small evolutionary scales and has the advantage that genes can be a priori selected to obtain an optimal phylogenetic signal. The second might be preferable at larger evolutionary scales (e.g. among protists) and allows the discovery of many other genes, which can shed light on the evolution of important features such as metabolic pathways.

Phylogenomic reconstruction methods based on gene content/order cannot be applied to incomplete genomic sampling, but those based on DNA strings and on primary sequences can. In the latter case, missing data will occur even when complete genomes are used, especially at a large evolutionary scale, since most, if not all, genes can be lost, duplicated or horizontally transferred in some organisms. To our knowledge, no theoretical reasons suggest that sequence-based approaches can not be used on incomplete alignments, i.e. containing cells coded as missing data. Nevertheless, "the problem of missing data is widely considered to be the most significant obstacle (…) in combining datasets (…) that do not include identical taxa", as suggested by empirical studies and computer simulations [START_REF] Wiens | Missing data, incomplete taxa, and phylogenetic accuracy[END_REF]. Two problems need to be distinguished: (1) the potential lack of resolution induced by the presence of taxa with too many missing cells, (2) the possible interaction between missing entries and artifact-inducing model violations.

Recent computer simulations using large number of characters [START_REF] Philippe | Phylogenomics of eukaryotes: impact of missing data on large alignments[END_REF][START_REF] Wiens | Missing data, incomplete taxa, and phylogenetic accuracy[END_REF] suggest that the inaccurate placement of incomplete taxa is not due to missing data but rather to the insufficient number of informative characters. As an extreme example, the tree reconstruction method remains accurate when positions have an average of four known and 32 unknown character states, because each species is nevertheless represented by about 3,000 amino acids [START_REF] Philippe | Phylogenomics of eukaryotes: impact of missing data on large alignments[END_REF]. However, the presence of missing cells unevenly distributed across the data matrix potentially affects estimates of model parameters. It is not yet clear how the induced model misspecifications in turn influence phylogenetic inference. Interestingly, it seems that the advantage of adding an incomplete taxon that breaks a long branch is greater than the disadvantage of the induced model misspecification (Wiens in press).

Few attempts at assessing the effect of missing data have been made with empirical data [START_REF] Bapteste | The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba[END_REF][START_REF] House | Using homolog groups to create a wholegenomic tree of free-living organisms: an update[END_REF], Philippe et al 2004). For instance, a bipartition of the supermatrix (25% of missing data) into the most complete genes and the less complete genes appears to be indistinguishable from random bipartitions of the same size [START_REF] Bapteste | The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba[END_REF], Philippe et al 2004). However, when the level of missing data is extreme (92%), the quality of the inference appears to be affected (e.g. strong support for the paraphyly of Glires and of Ecdysozoa) [START_REF] Driskell | Prospects for building the tree of life from large sequence databases[END_REF], despite the large size of the data set (70 taxa and 1131 genes). In summary, even if the problem generated by missing data has been overrated, additional work is needed to characterize its impact more precisely.

INFERENCE OF PHYLOGENOMIC TREES

The methods used in phylogenomic inference are of two kinds: (1) primary sequence based methods, which are very similar to the classical tree reconstruction, and for which several excellent reviews and textbooks are available [START_REF] Felsenstein | Inferring phylogenies[END_REF][START_REF] Holder | Phylogeny estimation: Traditional and Bayesian approaches[END_REF], Swofford et al 1996), and (2) methods above the sequence level [START_REF] Wolf | Genome trees and the tree of life[END_REF]. The two approaches are fundamentally similar, the main difference being the characters used. In both cases, we believe that probabilistic methods are more powerful and more reliable. First, they have a more robust theoretical justification in that they rely on an explicit account of their assumptions by using stochastic models describing the pattern of molecular evolution [START_REF] Felsenstein | Inferring phylogenies[END_REF]). Second, not only do they allow estimating the phylogeny and a confidence-level, but they also provide general methods to evaluate the fit of the model used [START_REF] Goldman | Statistical tests of models of DNA substitution[END_REF].

Approaches based on whole-genome features

Since genomes are the results of evolution, virtually any features comparable between organisms can be used to infer phylogenies, as evidenced by the plethora of new approaches recently published [START_REF] Fitz-Gibbon | Whole genome-based phylogenetic analysis of free-living microorganisms[END_REF][START_REF] Henz | Wholegenome prokaryotic phylogeny[END_REF][START_REF] House | Using homolog groups to create a wholegenomic tree of free-living organisms: an update[END_REF][START_REF] House | Geobiological analysis using whole genome-based tree building applied to the Bacteria, Archaea and Eukarya[END_REF][START_REF] Korbel | SHOT: a web server for the construction of genome phylogenies[END_REF][START_REF] Lin | Whole-genome trees based on the occurrence of folds and orthologs: implications for comparing genomes on different levels[END_REF][START_REF] Pride | Evolutionary implications of microbial genome tetranucleotide frequency biases[END_REF][START_REF] Qi | Whole proteome prokaryote phylogeny without sequence alignment: A K-string composition approach[END_REF][START_REF] Snel | Genome phylogeny based on gene content[END_REF][START_REF] Tekaia | The genomic tree as revealed from whole proteome comparisons[END_REF][START_REF] Wolf | Genome trees constructed using five different approaches suggest new major bacterial clades[END_REF][START_REF] Yang | Phylogeny determined by protein domain content[END_REF]. The justifications of whole genome tree approaches are generally that the phylogeny of organisms can not be equated to the phylogeny of single genes (such as rRNA) and that this classical phylogeny is sensitive to hidden paralogy, horizontal gene transfer (HGT) or tree reconstruction artifacts (Fitz-Gibbon & House 1999, [START_REF] Lin | Whole-genome trees based on the occurrence of folds and orthologs: implications for comparing genomes on different levels[END_REF][START_REF] Snel | Genome phylogeny based on gene content[END_REF][START_REF] Tekaia | The genomic tree as revealed from whole proteome comparisons[END_REF]. Since tree reconstruction artifacts can affect any approach and are difficult to detect, these new approaches based on various character types are of paramount importance to corroborate phylogenetic inference [START_REF] Miyamoto | Testing species phylogenies and phylogenetic methods with congruence[END_REF][START_REF] Swofford | When are phylogeny estimates from molecular and morphological data incongruent? In Phylogenetic analysis of DNA sequences[END_REF][START_REF] Wolf | Genome trees and the tree of life[END_REF].

Distribution of sequence strings

The frequencies of small oligonucleotides (up to 8) or oligopeptides (up to 6) observed in the genome or the proteome can be transformed into distances then used to construct phylogenies [START_REF] Blaisdell | A measure of the similarity of sets of sequences not requiring sequence alignment[END_REF]). Numerous well-accepted clades, including deep ones, were recovered using this approach, confirming that a phylogenetic signal is present in these characters [START_REF] Edwards | A genomic schism in birds revealed by phylogenetic analysis of DNA strings[END_REF][START_REF] Pride | Evolutionary implications of microbial genome tetranucleotide frequency biases[END_REF][START_REF] Qi | Whole proteome prokaryote phylogeny without sequence alignment: A K-string composition approach[END_REF]. However, several undisputed clades were significantly rejected, and the comparison of the branch lengths of the genome tree obtained using tetranucleotide frequencies with those of the tree obtained from the standard analysis of rRNA sequences [START_REF] Pride | Evolutionary implications of microbial genome tetranucleotide frequency biases[END_REF] suggests that the phylogenetic signal contained in sequence strings saturates rapidly. However, the methods proposed so far are extremely crude. Oligonucleotide or oligopeptide frequencies are transformed into distances without any underlying model of evolution. It is nevertheless remarkable that something considered as a bias in standard sequence-based methods [START_REF] Lockhart | Recovering evolutionary trees under a more realistic model of sequence evolution[END_REF]

Homology and orthology assessment

All other approaches require the definition of homologous, or more often, orthologous genes. By definition, the phylogenetic history displayed by orthologous genes is the organismal phylogeny [START_REF] Fitch | Distinguishing homologous from analogous proteins[END_REF]), whereas paralogous or xenologous genes display a combination of organismal and genespecific history. In practice, the identification of orthologous genes involves a certain amount of circularity because it requires an a priori knowledge of the organismal phylogeny. Indeed, finding orthologous genes is difficult since the organismal phylogeny is generally unknown (or at best partially known), and its reconstruction represents the goal. This problem is much akin to the alignment/phylogeny problem, in which alignment and phylogeny should be estimated simultaneously [START_REF] Sankoff | Evolution of 5S RNA and the nonrandomness of base replacement[END_REF][START_REF] Wheeler | Iterative pass optimization of sequence data[END_REF]), but the practical difficulties are such that the two steps are generally separated. Moreover, a careful phylogenetic reconstruction of all gene families is a Herculean labor most researchers want to avoid (but see [START_REF] Storm | Automated ortholog inference from phylogenetic trees and calculation of orthology reliability[END_REF]. Therefore, an operational, yet approximate, definition of orthology is used. Schematically, all genomes are compared to each other at the amino acid level; only the pairs of sequences that are best reciprocal hits are further considered. The clusters of orthologous groups are then constructed by a single linkage analysis of all orthologous pairs [START_REF] Sunderland | A genomic perspective on protein families[END_REF] or by Markov cluster algorithms [START_REF] Remm | Automatic clustering of orthologs and in-paralogs from pairwise species comparisons[END_REF]. These approaches, albeit reasonably effective in practice, do not guarantee the identification of orthologous genes, since, if a gene has been transferred from a distant organism and replaced the original copy, it will fulfill all the requirements while being xenologous. Information from synteny will probably improve the accuracy of orthology assignment [START_REF] Zheng | Using shared genomic synteny and shared protein functions to enhance the identification of orthologous gene pairs[END_REF], but much more work is needed before obtaining perfect assignment of orthologous genes.

Gene order methods

Since the character space of gene order data is huge, it probably constitutes the most promising genome feature based method. However, it is also technically and computationally the most difficult (see Moret et al 2005 for review). Briefly, distances can be computed by minimising the number of inversions, transpositions, insertions and deletions necessary to transform one unichromosomal genome into another [START_REF] Sankoff | Gene order comparisons for phylogenetic inference: evolution of the mitochondrial genome[END_REF], a method which is most often further simplified by considering only inversions [START_REF] Bourque | Genome-scale evolution: reconstructing gene orders in the ancestral species[END_REF]. Alternatively, the break-point distance between two genomes, defined as the minimum number of pairs of genes next to each other in one of the two genomes, but not in the other [START_REF] Nadeau | Lengths of Chromosomal Segments Conserved since Divergence of Man and Mouse[END_REF], can be used to infer phylogeny [START_REF] Blanchette | Breakpoint phylogenies[END_REF]. However, the evolution of gene order involves additional rearrangement mechanisms, not easily accountable for with these methods, such as translocation (i.e., transposition between chromosomes), and fusion or fission events of chromosomes.

More recently, methods based on maximum-likelihood distances [START_REF] Wang | Distance-based genome rearrangement phylogeny[END_REF], or on Bayesian inference [START_REF] Larget | A Bayesian analysis of metazoan mitochondrial genome arrangements[END_REF][START_REF] York | Bayesian estimation of the number of inversions in the history of two chromosomes[END_REF] have been proposed that take multiple changes in gene order into account. They revealed a non-negligible level of saturation [START_REF] York | Bayesian estimation of the number of inversions in the history of two chromosomes[END_REF] and a limited treeresolving power [START_REF] Larget | A Bayesian analysis of metazoan mitochondrial genome arrangements[END_REF]. The small size of mitochondrial genomes might explain the predominance of stochastic noise in the latter case. Prokaryotic genomes contain much more information, but are too large for analysis by current software. The drastic simplifying assumption that gene order can be reduced to the presence/absence of gene pairs allows the inference of prokaryotic phylogenies that are similar to the ones based on gene content [START_REF] Korbel | SHOT: a web server for the construction of genome phylogenies[END_REF][START_REF] Wolf | Genome trees constructed using five different approaches suggest new major bacterial clades[END_REF]. Further methodological and computational developments are needed to realize the full potential of gene-order methods.

Gene content methods

Phylogenomic analyses based on gene content generally use orthologs, but a few variants exist which use homologous instead of orthologous genes (Fitz-Gibbon & House 1999), protein domain content [START_REF] Yang | Phylogeny determined by protein domain content[END_REF], or fold occurrence [START_REF] Lin | Whole-genome trees based on the occurrence of folds and orthologs: implications for comparing genomes on different levels[END_REF][START_REF] Yang | Phylogeny determined by protein domain content[END_REF]. Interestingly, an orthologous gene present in all organisms (the Holy Grail of the sequence-based approach) has the same character state and is not informative for gene content methods. This is an important source of corroboration, since patterns informative for sequence-based approaches are not informative for gene content approaches.

Otherwise, the distribution of orthologs will be informative in the cases of (1) gene loss, which have a non-negligible probability of being convergent, (2) gene genesis, which is potentially the most informative, and (3) horizontal gene transfers (i.e. the acquisition of a new gene from a distantly related organism at the base of a clade will create a synapomorphy for this clade and a homoplasy for locating this clade). Gene content, albeit more integrated than primary sequences, is thus far from providing an unambiguous phylogenetic signal.

The binary matrices of presence/absence of homologous or orthologous genes can be analyzed by distance methods [START_REF] Lin | Whole-genome trees based on the occurrence of folds and orthologs: implications for comparing genomes on different levels[END_REF][START_REF] Snel | Genome phylogeny based on gene content[END_REF], parsimony (Fitz-Gibbon & House 1999) or Dollo parsimony [START_REF] Wolf | Genome trees constructed using five different approaches suggest new major bacterial clades[END_REF]. However, big/small genome attraction [START_REF] Lake | Deriving the genomic tree of life in the presence of horizontal gene transfer: conditioned reconstruction[END_REF] appears to affect all these methods, as demonstrated by the artificial grouping of unrelated species with small genomes (e.g. Mycoplasma, Buchnera, Chlamydia or Rickettsia). For instance, [START_REF] Copley | Systematic searches for molecular synapomorphies in model metazoan genomes give some support for Ecdysozoa after accounting for the idiosyncrasies of Caenorhabditis elegans[END_REF] demonstrated that phylogenies based on gene content and protein domain combinations support the paraphyly of Ecdysozoa, but are biased by a systematic high rate of character loss in nematodes. When this bias is accounted for by computing the number of losses expected randomly, a slight support for the monophyly of Ecdysozoa is recovered. Interestingly, this artifact yields the same inconsistency phenotype as the one displayed by sequence-based analyses [START_REF] Philippe | Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia[END_REF]: in both cases, arthropods are sister of vertebrates to the exclusion of nematodes. Such a convergence between the two methods could be explained by the fact that the fast evolving species are also those that have undergone the most extreme genome reduction. This observation weakens the strength of the corroboration between gene content and primary sequence approaches.

Altogether, gene-content phylogenies are not in excellent agreement with previous knowledge, even if we ignore the problem due to small genome attraction. The monophyly of the three domains (Archaea, Bacteria and Eukaryota) is always recovered, but the phylogeny within domains is much more problematic. For instance, Halobacterium never clusters with Methanosarcina, probably because of many HGTs from Bacteria which attract it towards the base of Archaea [START_REF] Korbel | SHOT: a web server for the construction of genome phylogenies[END_REF]; except with threshold parsimony [START_REF] House | Geobiological analysis using whole genome-based tree building applied to the Bacteria, Archaea and Eukarya[END_REF], the monophyly of Proteobacteria is never recovered [START_REF] Dutilh | The consistent phylogenetic signal in genome trees revealed by reducing the impact of noise[END_REF][START_REF] Gu | Genome phylogenetic analysis based on extended gene contents[END_REF][START_REF] Henz | Wholegenome prokaryotic phylogeny[END_REF][START_REF] Wolf | Genome trees constructed using five different approaches suggest new major bacterial clades[END_REF].

Several technical improvements have recently been proposed. Since big/small genome attraction is akin to the problem of compositional bias in sequence-based approaches, this non-phylogenetic signal can be reduced by the use of the logdet/paralinear transformation [START_REF] Lake | Deriving the genomic tree of life in the presence of horizontal gene transfer: conditioned reconstruction[END_REF]. A simple model of gene genesis and gene loss allows ML estimates of evolutionary distances [START_REF] Gu | Genome phylogenetic analysis based on extended gene contents[END_REF][START_REF] Huson | Phylogenetic trees based on gene content[END_REF]), but simulations showed that their performance appears to be slightly poorer than the performance of Dollo parsimony [START_REF] Huson | Phylogenetic trees based on gene content[END_REF]).

In summary, the methods for inferring trees based on whole genome features are at an early stage of their development, which might be comparable to that of sequence-based methods in the early 1970's. In particular, they generally lack a global probabilistic modeling. Numerous works are ongoing and it will be important to extensively evaluate the accuracy of present and future methods.

Approaches based on primary sequences

Supermatrix versus supertree

The question of how to analyze multiple datasets has been the subject of intense debate (for review see [START_REF] Bull | Partitioning and combining characters in phylogenetic analysis[END_REF][START_REF] De Queiroz | Separate versus combined analysis of phylogenetic evidence[END_REF]. In brief, three approaches are mainly used: (1) total evidence [START_REF] Kluge | A concern for evidence and a phylogenetic hypothesis of relationships among Epicrates (Boidae, Serpentes)[END_REF], in which all datasets are combined together, called hereafter the supermatrix approach, (2) separate analysis [START_REF] Miyamoto | Testing species phylogenies and phylogenetic methods with congruence[END_REF], in which the datasets are analyzed individually and resulting topologies are combined using consensus or supertree methods, called hereafter the supertree approach, (3) conditional combination [START_REF] Bull | Partitioning and combining characters in phylogenetic analysis[END_REF][START_REF] Lecointre | Total evidence requires exclusion of phylogenetically misleading data[END_REF], in which only the datasets considered as congruent are combined, called hereafter the conditional supermatrix approach. Generally, their respective advantages are: (1) minimizing stochastic error, (2) increasing the significance of corroboration, and (3) minimizing conflicting signal. In the practice of phylogenomics, the supermatrix is by far preferred [START_REF] Murphy | Molecular phylogenetics and the origins of placental mammals[END_REF], Philippe et al 2005[START_REF] Qiu | The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes[END_REF][START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF], followed by a few cases of the supertree method [START_REF] Daubin | A phylogenomic approach to bacterial phylogeny: evidence of a core of genes sharing a common history[END_REF][START_REF] Philip | The Opisthokonta and the Ecdysozoa may not be clades: Stronger support for the grouping of plant and animal than for animal and fungi and stronger support for the Coelomata than Ecdysozoa[END_REF].

Although [START_REF] Bull | Partitioning and combining characters in phylogenetic analysis[END_REF] state that "no rational systematist would suggest combining genes with different histories to produce a single reconstruction", the reliance in the power of the supermatrix approach is quite strong, and this methodology is generally applied. Against Bull et al., one may argue that discordant genes will each display a different discordant history, which will be averaged away through a combined analysis (but see [START_REF] Matte-Tailliez | Archaeal phylogeny based on ribosomal proteins[END_REF].

In most of the large multigene studies, the possibility of incongruencies is voluntarily minimised by selecting genes having a priori the same evolutionary history (e.g. single-copy genes [START_REF] Lerat | From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the gamma-Proteobacteria[END_REF][START_REF] Murphy | Molecular phylogenetics and the origins of placental mammals[END_REF][START_REF] Philip | The Opisthokonta and the Ecdysozoa may not be clades: Stronger support for the grouping of plant and animal than for animal and fungi and stronger support for the Coelomata than Ecdysozoa[END_REF], organellar genes [START_REF] Qiu | The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes[END_REF], Soltis et al 1999), or orthology assessment based on synteny [START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF]). But generally, the homogeneity of the datasets was not tested, indicating a strict application of the total evidence principle.

The problem of homogeneity is important since several recurrent processes (gene duplication, HGT, or lineage sorting) can lead to incongruent gene trees. In particular the very notion of a Tree of Life has been questioned because of rampant HGTs [START_REF] Doolittle | Phylogenetic classification and the universal tree[END_REF]. Several studies have nevertheless argued that a strong phylogenetic signal is present in the prokaryotic genomes (for review see [START_REF] Brown | Ancient horizontal gene transfer[END_REF][START_REF] Philippe | Horizontal gene transfer and phylogenetics[END_REF] and therefore that HGTs do not wipe out the notion of organismal phylogeny. It is clear however, that few, and probably none, of the genes have followed exactly the organismal phylogeny during the entire history of life on Earth. Thus, HGTs constitute a source of nuisance that should be addressed to improve the accuracy of phylogenetic reconstructions.

The infrequent use of the conditional supermatrix approach is likely not due to a philosophical rejection of its principle, but rather to the difficulty of detecting homogeneous datasets in practice. The Incongruence Length Difference (ILD) test [START_REF] Farris | Testing significance of incongruence[END_REF] was initially designed for parsimony and has been recently adapted to distance methods [START_REF] Zelwer | Detecting phylogenetic incongruence using BIONJ: an improvement of the ILD test[END_REF]).

However, the interpretation of this test is complicated by the fact that stochastic noise can generate by itself significant results [START_REF] Dolphin | Noise and incongruence: interpreting results of the incongruence length difference test[END_REF]. Its efficiency in detecting incongruence and determining data combinability has been repeatedly questioned [START_REF] Darlu | When does the incongruence length difference test fail?[END_REF][START_REF] Dowton | Increased congruence does not necessarily indicate increased phylogenetic accuracy--the behavior of the incongruence length difference test in mixed-model analyses[END_REF][START_REF] Yoder | Failure of the ILD to determine data combinability for slow loris phylogeny[END_REF]. Parametric incongruence tests have also been proposed for ML methods [START_REF] Huelsenbeck | A likelihood ratio test to detect conflicting phylogenetic signal[END_REF], or in a Bayesian framework [START_REF] Nylander | Bayesian phylogenetic analysis of combined data[END_REF], although they are not yet in widespread use. A promising model has been proposed, in which each gene can, with a certain prior probability, choose between either conforming to the common global topology or relying on its own topology [START_REF] Suchard | Hierarchical phylogenetic models for analyzing multipartite sequence data[END_REF].

Nevertheless, the conditional supermatrix approach is used, despite the fact that the criteria used to discard gene/sequence are not well validated (Brochier et al 2002[START_REF] Brown | Universal trees based on large combined protein sequence data sets[END_REF][START_REF] Lecointre | Total evidence requires exclusion of phylogenetically misleading data[END_REF][START_REF] Matte-Tailliez | Archaeal phylogeny based on ribosomal proteins[END_REF]. However, the nature of the test used to decide whether a gene significantly supports a different topology yields divergent interpretations of the same data with regards to the importance of HGTs [START_REF] Bapteste | Phylogenetic reconstruction and lateral gene transfer[END_REF][START_REF] Lerat | From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the gamma-Proteobacteria[END_REF][START_REF] Zhaxybayeva | Genome mosaicism and organismal lineages[END_REF]. An improvement of these incongruence tests [START_REF] Goldman | Likelihood-based tests of topologies in phylogenetics[END_REF] therefore constitutes an important avenue of future research.

The supertree approach is most often used to combine trees from the literature that were obtained from diverse sources of data [START_REF] Sanderson | Phylogenetic supertrees: assembling the trees of live[END_REF], the most popular method being Matrix Representation with Parsimony (MRP, [START_REF] Baum | Combining trees as a way of combining data sets for phylogenetic inference, and the desirability of combining gene trees[END_REF][START_REF] Ragan | Matrix representation in reconstructing phylogenetic relationships among the eukaryotes[END_REF]. The comparative efficiency of supermatrix and supertree approaches is poorly studied, especially in a phylogenomic context [START_REF] Gatesy | Inconsistencies in arguments for the supertree approach: supermatrices versus supertrees of Crocodylia[END_REF][START_REF] Philip | The Opisthokonta and the Ecdysozoa may not be clades: Stronger support for the grouping of plant and animal than for animal and fungi and stronger support for the Coelomata than Ecdysozoa[END_REF]. We have therefore analyzed the dataset of 71 genes (Figure 1a) using a supertree MRP approach. Interestingly, almost all nodes inferred from the supermatrix or using the supertree are identical, except the position of urochordates within deuterostomes (not shown) and the relationships among protostomes, indicating an excellent congruence of the supermatrix and supertree (only 3 differences for 46 bipartitions). When the supertree is reconstructed with MP (Figure 1c), platyhelminthes are grouped with tardigrads+nematods instead of with other lophotrochozoans (annelids and molluscs), disrupting the monophyly of both Ecdysozoa and Lophotrochozoa.

When the supertree is reconstructed using Bayesian inference (Figure 1d), the results were slightly more consistent, since the monophyly of Ecdysozoa, but not Lophotrochozoa, was recovered. The MRP supertree approach appears to have difficulty in placing the fast-evolving lineages (e.g. Platyhelminthes). The synergy among all positions in the supermatrix might explain the ability of the method to better deal with LBA artifacts. However, refined studies are urgently required to evaluate the relative accuracy and efficiency of the supermatrix and supertree methods.

Supermatrix: scaling up current methods

In contrast to the genome-feature approaches mentioned above, the main advantage of sequence-based methods is that their properties have been intensively explored, tested and validated, so that many of their strengths and weaknesses are known. Indeed, in most sequence-based phylogenomic analyses published to date, almost the same protocols as for single gene studies have been applied. The congruence among results obtained by different methods is high, with some notable exceptions [START_REF] Canback | A phylogenomic study of endosymbiotic bacteria[END_REF][START_REF] Goremykin | The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm[END_REF][START_REF] Soltis | Phylogeny of seed plants based on evidence from eight genes[END_REF][START_REF] Stefanovic | Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?[END_REF]. These incongruencies confirm the presence of non-phylogenetic signal, and suggest that the increase in resolution obtained by analyzing larger datasets is not in itself a guarantee of accuracy.

Conversely, the agreement between the methods does not mean that the obtained tree is correct (see Brinkmann et al accepted).

The dramatic change of scale of the data matrices implies the need for a corresponding increase in computational power, in particular for probabilistic methods. Two factors have to be considered. First, there is a simple scaling-up of both memory requirement and computational load. Second, the reliability of the heuristic search procedures underlying ML programs, or the Monte Carlo devices of the Bayesian samplers, is anything but guaranteed. A particular concern is that conflicting signals result in the presence of many secondary maxima in the space of tree topologies, separated by high potential barriers. Standard procedures are likely to get trapped in these local optima [START_REF] Salter | Complexity of the likelihood surface for a large DNA dataset[END_REF], and thus do not yield reliable phylogenetic estimates.

A number of algorithmic innovations have led to better heuristic searches in the space of topologies: genetic algorithms [START_REF] Brauer | Genetic algorithms and parallel processing in maximum-likelihood phylogeny inference[END_REF][START_REF] Lemmon | The metapopulation genetic algorithm: An efficient solution for the problem of large phylogeny estimation[END_REF], disk-covering methods [START_REF] Huson | Disk-covering, a fast-converging method for phylogenetic tree reconstruction[END_REF], or parallelized computing [START_REF] Keane | DPRml: distributed phylogeny reconstruction by maximum likelihood[END_REF]. In the case of the Bayesian methods, an interesting approach has been proposed consisting of using coupled 'heated' Monte Carlo Markov Chains (MCMC), which can be easily parallelized [START_REF] Altekar | Parallel Metropolis coupled Markov chain Monte Carlo for Bayesian phylogenetic inference[END_REF][START_REF] Feng | Parallel algorithms for Bayesian phylogenetic inference[END_REF]. Thanks to these advances, maximum likelihood and Bayesian phylogenetic reconstructions will soon be able to handle phylogenomic datasets. However, their overall reliability has been evaluated mainly on simulated data, which are probably much more 'funnel-shaped' towards the true phylogeny than are real sequences (Brinkmann et al accepted, Stamatakis et al 2005).

Another possible stance towards efficient tree space searches is to restrict the analysis, by constraining nodes that have been found with high support when each gene of the concatenation was analyzed separately [START_REF] Philippe | Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia[END_REF].

The number of trees compatible with these constraints is still large, but accessible to an exhaustive analysis. Such approaches may not be considered as a definitive method, but could provide a good proxy.

Toward more complex models

If the availability of large data matrices poses new computational challenges to probabilistic methods, it allows the development of more realistic, parameter-rich, models. As long as the number of parameters increases more slowly than the number of sites, a model does not fall into the infinitely many parameter trap [START_REF] Felsenstein | Inferring phylogenies[END_REF], and thus, has good consistency properties.

Of course, no probabilistic model will ever capture evolutionary patterns in their full complexity, but their most important aspects, at least the ones that cause inconsistency of the current methods, can be accounted for [START_REF] Steel | Should phylogenetic models be trying to 'fit an elephant[END_REF]. The main idea would not so much be an increased realism, but an improved flexibility accounting for the diverse kinds of heterogeneities and disparities of the substitution processes.

One possible research direction is to account for disparities in the evolutionary process across the genes that make up the concatenation. A simple solution is to constrain the model to have a global topology, but gene-specific branch lengths [START_REF] Yang | Maximum-likelihood models for combined analyses of multiple sequence data[END_REF]. More generally, any parameter other than the topology can be considered as gene-specific in such "separate" (or partitioned) models. Another possible avenue of research is to account for site-specific patterns of substitution, using mixture models [START_REF] Kolaczkowski | Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous[END_REF][START_REF] Lartillot | A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process[END_REF][START_REF] Pagel | A phylogenetic mixture model for detecting patternheterogeneity in gene sequence or character-state data[END_REF]. A mixture model combines several different classes to describe the substitution process, each of which is characterized by its own set of parameters (e.g. equilibrium frequencies or exchangeability probabilities).

Thus far, few studies have tried to address the relative performances of alternative probabilistic models on phylogenomic datasets. A recent study (Brinkmann et al accepted) has confirmed that accounting for site-specific rates, or having a good empirical substitution rate matrix are important factors yielding a higher phylogenetic accuracy. On the other hand, separate models, in spite of their overall better statistical fit, do not seem to fundamentally improve phylogenetic inference. Since separate models handle a substantial part of heterotachy, this suggests that heterotachy, despite recent interest, may not constitute a major source of systematic bias. In any case, a much wider analysis of the impact of model choice on the prevalence of artifacts in phylogenomic inference has to be carried out.

Reducing systematic errors through data exclusion

Phylogenomic datasets contain a large amount of genuine phylogenetic signal but they also contain non-phylogenetic signals that current methods of tree reconstruction are not able to perfectly handle. To avoid the perils of inconsistency, one can take advantage of the fact that the quantity of phylogenetic signal is no longer a serious limiting factor. More precisely, the part of the datasets that contains mainly non-phylogenetic signals can be excluded, allowing the concentration of phylogenetic signal in the remaining dataset. This increase of the phylogenetic to non-phylogenetic signal ratio reduces the probability of inconsistency, even without the use of improved tree reconstruction methods.

The rationale of most methods of data exclusion is straightforward: tree reconstruction artifacts are due to multiple substitutions that are not correctly identified as convergences or reversions by inference methods [START_REF] Olsen | Earliest phylogenetic branching: comparing rRNA-based evolutionary trees inferred with various techniques[END_REF].

The simplest possibility consists in removing the fast evolving species, which by definition accumulate multiple substitutions. This approach efficiently reduces the misleading effect of the rate signal [START_REF] Philippe | Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia[END_REF][START_REF] Stefanovic | Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?[END_REF]. In many cases, the exclusion of odd taxa [START_REF] Sanderson | Troubleshooting molecular phylogenetic analyses[END_REF] is implicit, since investigators never envision using them in their analyses (e.g. microsporidia as a fungal representative).

When all of the available species representing a clade of interest are fastevolving, the specific removal of the fastest evolving sequences of this clade from the supermatrix appears to be efficient. For example, when a complete supermatrix of 133 genes is used, all tree reconstruction methods strongly, albeit artifactually, locate microsporidia at the base of eukaryotes, but probabilistic methods with a complex model (WAG+F+Γ) avoid this LBA artifact when >70% of the microsporidial sequences are coded as missing data (Brinkmann et al accepted). It is interesting to note that a highly incomplete taxon (70% missing data) is more accurately located than a complete one. A similar approach, in which genes in their totality are discarded, was successful at avoiding the attraction between the fast evolving nematodes and platyhelminthes (Philippe et al 2005) or between nematodes and outgroup [START_REF] Dopazo | Genome-scale evidence of the nematode-arthropod clade[END_REF].

Interestingly, in the first case, the statistical support for the monophyly of Ecdysozoa and Lophotrochozoa increases when more and more genes are discarded (as long as more than 40 genes are considered).

These approaches are rather crude, since sequences from genes and/or species are discarded in their totality. Nevertheless, even if these sequences contain much non-phylogenetic signal, they likely also contain some phylogenetic signal. Refined methods have been proposed to selectively eliminate fast evolving characters in part [START_REF] Lopez | The root of the tree of life in the light of the covarion model[END_REF] or completely [START_REF] Brinkmann | Archaea sister group of Bacteria? Indications from tree reconstruction artifacts in ancient phylogenies[END_REF][START_REF] Burleigh | Phylogenetic signal in nucleotide data from seed plants: Implications for resolving the seed plant tree of life[END_REF][START_REF] Dutilh | The consistent phylogenetic signal in genome trees revealed by reducing the impact of noise[END_REF][START_REF] Pisani | Identifying and removing fast-evolving sites using compatibility analysis: An example from the arthropoda[END_REF][START_REF] Ruiz-Trillo | Acoel flatworms: earliest extant bilaterian Metazoans, not members of Platyhelminthes[END_REF]. In several cases, taxa that emerged at the base of the tree when all the characters are used are relocated later in the tree when fast evolving positions are removed, strongly suggesting that their observed basal position is due to an LBA artifact (Brochier & Philippe 2002, Philippe et al 2000[START_REF] Pisani | Identifying and removing fast-evolving sites using compatibility analysis: An example from the arthropoda[END_REF]). However, the number of remaining slowly evolving positions is often too small to yield high statistical support for most of the clades in singlegene analyses, but not in phylogenomic analyses [START_REF] Burleigh | Phylogenetic signal in nucleotide data from seed plants: Implications for resolving the seed plant tree of life[END_REF][START_REF] Delsuc | Phylogenomics and the reconstruction of the tree of life[END_REF].

The RY coding strategy [START_REF] Woese | Archaeal phylogeny: reexamination of the phylogenetic position of Archaeoglobus fulgidus in light of certain composition-induced artifacts[END_REF] discards all fast-evolving transitions and improves inference without drastically compromising the resolution [START_REF] Phillips | Genome-scale phylogeny and the detection of systematic biases[END_REF]. Importantly, this coding not only addresses the problem of rate signal but also of compositional signal. Indeed, the GC content can be extremely variable among homologous sequences from various organisms whereas the frequency of purines is remarkably homogeneous [START_REF] Woese | Archaeal phylogeny: reexamination of the phylogenetic position of Archaeoglobus fulgidus in light of certain composition-induced artifacts[END_REF]. This constitutes a method of choice to avoid inconsistency due to compositional bias.

Finally, it is possible to remove characters whose evolutionary history violates most the assumptions of the underlying model of sequence evolution instead of the fastest evolving characters/species. In fact, the RY coding eliminates transitions that are mainly responsible for the non-stationarity of the nucleotide composition (see Hrdy et al 2004 for a similar approach in the case of proteins). Similarly, the constant sites violate the assumptions about the distribution of site rates [START_REF] Lockhart | Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis[END_REF] and their elimination constitutes an efficient way of improving inference [START_REF] Hirt | Microsporidia are related to fungi: evidence from the largest subunit of RNA polymerase II and other proteins[END_REF][START_REF] Phillips | Genome-scale phylogeny and the detection of systematic biases[END_REF].

Heterotachous sites violate the assumption that the evolutionary rate of a position is constant through time, made by all current models except the covarion model [START_REF] Fitch | An improved method for determining codon variability in a gene and its application to the rate of fixation of mutations in evolution[END_REF]. As expected, the elimination of these sites had reduced LBA artifacts in the case of the eukaryotic phylogeny [START_REF] Inagaki | Covarion Shifts Cause a Long-Branch Attraction Artifact That Unites Microsporidia and Archaebacteria in EF-1{alpha} Phylogenies[END_REF], Philippe & Germot 2000).

We believe that these data removal approaches are complementary to the improvement of tree reconstruction methods through the implementation of more realistic models of sequence evolution. So far, they are also more readily accessible in practice, since they are less demanding in terms of the complexity of bioinformatic methods and computational time.

CONCLUDING REMARKS

In this review, we have emphasized that inconsistency of tree reconstruction methods constitutes the major limitation of phylogenomics. We have explored ways to reduce its impact, whether at the level of data assembly or at the level of tree building per se. Adequate taxon sampling, probabilistic methods and the exclusion of phylogenetically misleading data constitute the three most important criteria to obtain reliable phylogenomic trees. However, this might not be sufficient to ensure that the inferred tree is the correct one. We believe that an important test of the reliability of tree reconstruction methods is their robustness with respect to species sampling. In particular, a reliable method should be able to recover exactly the same topology with a taxon-poor and a taxon-rich sampling, which is far from being the case at present (Figure 1). This is particularly important to be confident in the phylogenetic location of poorly diversified clades (e.g. Amborella, or monotremes). The stability of phylogenies in face of variation of species sampling will constitute one of the best guarantees that the non-phylogenetic signal has been correctly handled.

The correctness of inferences should also be verified via corroboration from independent sources [START_REF] Miyamoto | Testing species phylogenies and phylogenetic methods with congruence[END_REF][START_REF] Swofford | When are phylogeny estimates from molecular and morphological data incongruent? In Phylogenetic analysis of DNA sequences[END_REF]. When complete genomes are used, this may appear hopeless (even if internal verifications of homogeneity are possible). However, genomes can be conveniently subdivided into various character types that can be considered as more or less independent: oligonucleotide composition, sequences of orthologous gene, gene content, and gene order. If inferences based on these "independent" sets of characters converge to the same results, an increased confidence can be placed in the corresponding phylogeny, although the same bias can theoretically affect several approaches. We therefore encourage the developments of sophisticated probabilistic methods for all types of data, and not only for primary sequences.

In the long term, one might envision that the Tree of Life, or at least its global scaffold, will be established within the next 10 years. Then, to numerous molecular systematists the important question would be: what next?

Figure 1

 1 Figure1Super-matrix, taxon sampling and super-tree in animal phylogenomics.

Figure

  Figure 1: Philippe et al.

  contains a phylogenetic signal, but it is not yet clear whether accurate methods can be developed to extract it.

Acknowledgements

We wish to thank David Bryant, Emmanuel Douzery, David Moreira, Davide Pisani, Nicolas Rodrigue, Naiara Rodríguez-Ezpeleta, Béatrice Roure and Brad Shaffer for critical readings of the manuscript. This work was supported by operating funds from Genome Québec. H.P. is a member of the Program in Evolutionary Biology of the Canadian Institute for Advanced Research (CIAR), which is acknowledged for salary and interaction support. H.P. is also grateful to the Canada Research Chairs Program and the Canadian Foundation for Innovation for salary and equipment support.