Gustavo A Bravo
email: gustavo_bravo@fas.harvard.edu

Alexandre Antonelli

Christine D Bacon

Krzysztof Bartoszek

Mozes P K Blom

Stella Huynh

Graham Jones

L Lacey Knowles

Sangeet Lamichhaney

Thomas Marcussen

Hélène Morlon

Luay K Nakhleh

Bengt Oxelman

Bernard Pfeil

Alexander Schliep

Niklas Wahlberg

Fernanda P Werneck

John Wiedenhoeft

Sandi Willows-Munro

Scott V Edwards

Collection: "Endless forms: Advances in evolutionary analyses of biodiversity

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to 2 benefit by embracing the many heterogeneous genomic signals emerging from the first decade of largescale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Introduction

Phylogenomics has been greatly enriched with the introduction of high-throughput sequencing (HTS) and increased breadth of phylogenomic sampling, which have allowed researchers interested in the ToL to scale up in several dimensions and placing both fields squarely in the era of 'big data'. Additionally, conceptual advances and improvements of statistical models used to analyze these data are helping bridge what some have perceived as a gap between phylogenetics and phylogeography (e.g., [START_REF] Felsenstein | Phylogenies from molecular sequences: inference and reliability[END_REF][START_REF] Huson | Application of Phylogenetic Networks in Evolutionary Studies[END_REF]Edwards et al., 2016a). However, as datasets become larger, researchers are inevitably faced with a plethora of heterogeneous signals across different genomic regions that often depart from a dichotomously-branching phylogeny [START_REF] Kunin | The net of life: reconstructing the microbial phylogenetic network[END_REF][START_REF] Jeffroy | Phylogenomics: the beginning of incongruence?[END_REF][START_REF] Mallet | How reticulated are species?[END_REF]. These signals cover an increasingly large array of biological processes at the level of genes and genomes, as well as individual organisms and populations, including processes such as recombination, hybridization, gene flow, and polyploidization; they can be thought of as conflicting, but in truth, they are simply a record of the singular history that we commonly refer to as the ToL [START_REF] King | Embracing Uncertainty in Reconstructing Early Animal Evolution[END_REF]. One of the grand challenges of evolutionary biology is deciphering this history, whether at the level of genes, populations, species, or genomes, however, currently phylogenomicists have not yet determined how to fully exploit the diverse signals in molecular data.

In this perspective piece, our goal is to highlight opportunities for the field of phylogenomics to overcome conceptual and practical challenges as we navigate our way through the era of big data. We argue that conceptually and analytically embracing heterogeneity generated by population-level processes and associated with different histories across the genome will lead to increased insight into the ToL and its underlying processes. Because major theoretical advances in phylogenomics and population genetics resulting from the advent of genome-scale data have been reviewed elsewhere (e.g., [START_REF] Delsuc | Phylogenomics and the reconstruction of the Tree of Life[END_REF][START_REF] Degnan | Gene tree discordance, phylogenetic inference and the multispecies coalescent[END_REF][START_REF] Ellegren | Determinants of genetic diversity[END_REF][START_REF] Payseur | A genomic perspective on hybridization and speciation[END_REF][START_REF] Hahn | Molecular Population Genomics[END_REF], here we focus on advances and challenges to phylogenomics specifically brought about by population genetic processes, which inevitably leads us to focus on the major conceptual framework dealing with such processes, the multispecies coalescent model (MSC). Likewise, advances in other methods such as supertrees (reviewed by [START_REF] Warnow | Supertree Construction: Opportunities and Challenges[END_REF], supermatrices [START_REF] De Queiroz | The supermatrix approach to systematics[END_REF][START_REF] De Queiroz | The supermatrix approach to systematics[END_REF][START_REF] Philippe | Pitfalls in supermatrix phylogenomics[END_REF], and partitioned analyses (e.g., [START_REF] Lanfear | Selecting optimal partitioning schemes for phylogenomic datasets[END_REF]Kainier & Lanfear, 2015) are not central to the objectives of this perspective piece and further details on those topics can be found elsewhere.

A key concept introduced by the scaling up from phylogeography to phylogenomics is the continuum of processes and analytical methods, -the so-called phylogeography-phylogenetics continuum (Edwards et al., 2016a). We argue here that bridging this continuum is critical for advancing phylogenetics. This can be done by either developing phylogenomic approaches that acknowledge and explicitly account for phylogeographic processes, or by determining the regions of parameter space (e.g., branch lengths in tree, level of gene flow) if any, where such within-species processes are not relevant.

For example, the choice of markers in a given phylogenomics project is currently guided more by convenience and cost than by evaluating the biological properties and phylogenomic signals in those data; but comparisons of signals across various types of markers (e.g., transcriptomes, noncoding regions) reveal that marker choice is a critical step toward shedding light on the history of populations and unraveling potential processes underlying such history [START_REF] Rokas | Genome-scale approaches to resolving incongruence in molecular phylogenies[END_REF][START_REF] Cutter | Integrating phylogenetics, phylogeography and population genetics through genomes and evolutionary theory[END_REF]Jarvis et al., 2014;[START_REF] Reddy | Why Do Phylogenomic Data Sets Yield Conflicting Trees? Data Type Influences the Avian Tree of Life more than Taxon Sampling[END_REF]. On the analysis side, we are in desperate need of methods that can handle the increasingly large data sets being produced by empiricists, but at the same time there is a desire to include increasingly diverse sources of signal in estimates of divergence times, biogeographic history, and models of diversification [START_REF] Delsuc | Phylogenomics and the reconstruction of the Tree of Life[END_REF][START_REF] Jeffroy | Phylogenomics: the beginning of incongruence?[END_REF][START_REF] Kumar | Statistics and Truth in Phylogenomics[END_REF]. Finding the balance between breadth, depth, and computational feasibility in project design and statistical analysis is crucial for the field today.

Additionally, we discuss ways in which data archiving and integration can benefit access to phylogenomics data and the contributions of phylogenomics to society. Are the priorities that society places on the many disciplines feeding into scientific efforts toward the ToL -fieldwork, museum collections, databases -appropriate for this grand mission? Although we cannot possibly answer all of these questions within the scope of this perspective, we hope to at least spur discussion on the wide range of field, laboratory, conceptual, and societal issues that allow phylogenomics to move forward.

We first describe the types of genomic data that researchers can generate to perform phylogenomic analyses and how those are more or less suitable for phylogenomic and phylogeographic analyses. We then discuss key concepts around the MSC and highlight the need to expand this model beyond its current limitations. We then discuss how the interplay between phylogenomics and macroevolution might strengthen our understanding of diversity patterns and offer suggestions as to how the community can overcome limitations posed by current methods and models in both fields. Finally, we discuss desired practices that, as a community, phylogenomicists, museum scientists, field biologists, bioinformaticians, and other scientists can embrace toward the goal of assembling the ToL.

Survey methodology

During the 'Origins of Biodiversity Workshop' organized during May 15-19, 2017 by Chalmers University of Technology and the University of Gothenburg, Sweden, under the auspices of the Gothenburg Centre for Advanced Studies (GoCAS), we gathered scholars and students from several countries and scientific backgrounds to discuss future perspectives in the fields of phylogenomics and phylogeography. We spent one week sharing our recent experiences in these fields and outlining the topics presented here and continued remotely to complete this review. Our goal is not to provide a complete overview of phylogenomic and phylogeographic research, but rather present a number of conceptual and practical aspects that we feel are essential to keep the momentum that these fields have gained during recent times.

Data generation and data types in phylogenomics

The need for large-scale phylogenomic data

One of the fundamental challenges in evolutionary biology is to estimate a ToL for all species. The potential impact of such large phylogenies is reflected in their publication in the highest impact journals, but also in their broad contribution, which extends beyond big data, to methodological innovations, and downstream understanding of macroevolutionary processes (e.g., coalescent methods of species tree inference; accounting for hybridization and unsampled species or localities in datasets; understanding community or genome evolution through large-scale phylogenetics). Hence, the phylogenomics community is now placing a high priority on very large-scale trees, whether in terms of number of taxa, number of genes, or both. The current need for large phylogenies and the high priority placed on them by high-impact journals can also result in short-cuts, wherein extant species lacking any molecular data are often placed in trees based on current taxonomy [START_REF] Jetz | The global diversity of birds in space and time[END_REF][START_REF] Zanne | Three keys to the radiation of angiosperms into freezing environments[END_REF][START_REF] Faurby | A species-level phylogeny of all extant and late Quaternary extinct mammals using a novel heuristic-hierarchical Bayesian approach[END_REF], which can result in conflicts with more robust estimates based on actual data [START_REF] Brown | The development of scientific consensus: Analyzing conflict and concordance among avian phylogenies[END_REF]. At the same time, however, hypothesis-testing in areas such as macroevolution, macroecology, biodiversity, and systematics require these large-scale trees, even as they present challenges being built on high quality data. The phylogenetic knowledge on which we lay a foundation for downstream analyses must be robust, and therefore it is essential that the input phylogenetic hypotheses themselves are robust [START_REF] Pyron | Post-molecular systematics and the future of phylogenetics[END_REF]. Indeed, the current bottlenecks in large-scale phylogenomic data do not appear to be the sequencing, but rather the compilation of high quality, wellcurated genomic resources and the availability of adequate software and methods that can fuel phylogenomics for the next century (e.g., Global Genome Initiative, www.mnh.si.edu/ggi/).

Data quality

Genome-scale data in the form of multiple alignments and other homology statements are the foundation of phylogenomics. A major challenge is the difficulty of comprehensive quality checks of data, given that HTS datasets are so large. As researchers collect datasets consisting of thousands of alignments across scores of species and data quality is a serious concern that is left for detection and handling primarily by computer algorithms. In addition to inherent systematic errors in the data [START_REF] Kocot | Phylogenomics of Lophotrochozoa with Consideration of Systematic Error[END_REF], several examples of errors in phylogenomic data sets have been reported in the literature, including the use of unintended paralogous sequences in alignments (e.g., [START_REF] Struck | The Impact of Paralogy on Phylogenomic Studies -A Case Study on Annelid Relationships[END_REF]; mistaking the genome sequence of one species for another [START_REF] Philippe | Resolving Difficult Phylogenetic Questions: Why More Sequences Are Not Enough[END_REF]; and inclusion of genome sequence from parasites into the genome of the host [START_REF] Kumar | Blobology: exploring raw genome data for contaminants, symbionts and parasites using taxon-annotated GC-coverage plots[END_REF]. However, the incidence of smaller errors in alignments that are not easily discerned from natural allelic variation, such as base miscalls or misplaced indels, are probably much more widespread than has been reported in the literature. Combined with the sensitivity of some phylogenomic datasets to individual loci or Single Nucleotide Polymorphisms (SNP) within loci [START_REF] Shen | Contentious relationships in phylogenomic studies can be driven by a handful of genes[END_REF], such errors could have damaging consequences for phylogenomic studies, for the inference of topology and branch lengths of both gene phylogenies and species phylogenies [START_REF] Marcussen | Ancient hybridizations among the ancestral genomes of bread wheat[END_REF][START_REF] Bleidorn | Sources of Error and Incongruence in Phylogenomic Analyses[END_REF]. Furthermore, as phylogenomic datasets increase in size, it is likely that the accumulation of errors due to the use of different sequencing chemistries and sequencing depths [START_REF] Quail | A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers[END_REF][START_REF] Goodwin | Coming of age: ten years of next-generation sequencing technologies[END_REF] will ultimately influence phylogenetic inference. We predict that the impact of these errors will largely depend on the sampling breadth and taxonomic scale of each study, and whether the species phylogeny is a tree or a network.

Sequencing high-quality samples from well-archived voucher specimens is a good first step to increase reproducibility and alleviate issues related to sample identity [START_REF] Peterson | The need for proper vouchering in phylogenetic studies of birds[END_REF][START_REF] Pleijel | Phylogenies without roots? A plea for the use of vouchers in molecular phylogenetic studies[END_REF][START_REF] Chakrabarty | Genetypes: a concept to help integrate molecular phylogenetics and taxonomy[END_REF][START_REF] Turney | Non-repeatable science: assessing the frequency of voucher specimen deposition reveals that most arthropod research cannot be verified[END_REF][START_REF] Troudet | The Increasing Disconnection of Primary Biodiversity Data from Specimens: How Does It Happen and How to Handle It?[END_REF]. For individual phylogenomic studies, wholesale manual inspection of every locus is unsustainable [START_REF] Irisarri | Phylotranscriptomic consolidation of the jawed vertebrate timetree[END_REF]), but spot checks of a subset of the data (e.g., 5-10% of the alignments) is a recommended best practice (Phillipe et al., 2011) that is beginning to be encouraged in peer review and in published papers [START_REF] Montague | Comparative analysis of the domestic cat genome reveals genetic signatures underlying feline biology and domestication[END_REF][START_REF] Liu | Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary[END_REF]. Such checking is important not only for new data generated by a given study, but also for data downloaded from public repositories such as NCBI and Orthomam, which are well known to contain errors [START_REF] Wesche | DNA Sequence Error Rates in Genbank Records Estimated Using the Mouse Genome as a Reference[END_REF][START_REF] Ranwez | OrthoMaM: A database of orthologous genomic markers for placental mammal phylogenetics[END_REF][START_REF] Douzery | OrthoMaM v8: A Database of Orthologous Exons and Coding Sequences for Comparative Genomics in Mammals[END_REF]. Because several databases do not include the raw sequence data it is often impossible to evaluate whether oddities may derive from poor sequencing. Robust pipelines for flagging poorly aligned sites or non-homologous sequences, based on existing tools or novel scripts such as Gblocks [START_REF] Castresana | Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis[END_REF][START_REF] Talavera | Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments[END_REF] or TrimAl (Capella-Gutiérrez, Silla-Martínez & Gabaldon, 2009) are gradually being put into practice [START_REF] Marcussen | Ancient hybridizations among the ancestral genomes of bread wheat[END_REF][START_REF] He | Reducing long-branch effects in multi-protein data uncovers a close relationship between Alveolata and Rhizaria[END_REF][START_REF] Irisarri | Phylotranscriptomic consolidation of the jawed vertebrate timetree[END_REF].

Coding regions, whether derived from transcriptomes or whole-genome data, are particularly amenable to spot checking of alignments and to filtering out of low-quality data with bioinformatic pipelines (e.g., [START_REF] Dunn | Agalma: an automated phylogenomics workflow[END_REF][START_REF] Blom | EAPhy: A Flexible Tool for High-throughput Quality Filtering of Exon-alignments and Data Processing for Phylogenetic Methods[END_REF]. Coding regions have the advantage of allowing amino acids to guide alignments, which is particularly useful for highly divergent sequences. Stop codons can help flag errors or genuine pseudogenes. Examining gene tree topologies is also widely used to detect paralogs in phylogenomic data (e.g., [START_REF] Betancur | Conserved genes, sampling error, and phylogenomic inference[END_REF]. Examining gene trees for aberrantly long branch lengths can also reveal misalignments (e.g., [START_REF] He | Reducing long-branch effects in multi-protein data uncovers a close relationship between Alveolata and Rhizaria[END_REF]; sensitivity analyses of methods for indirectly detecting errors in alignments are sorely needed.

Data generation and marker development

Genome reduction methods

A growing number of genome reduction methods are now providing empiricists with the means to generate genomic subsets suitable for phylogenetic and phylogeographic inference (reviewed by [START_REF] Mccormack | Applications of next-generation sequencing to phylogeography and phylogenetics[END_REF][START_REF] Leaché | The Utility of Single Nucleotide Polymorphism (SNP) Data in Phylogenetics[END_REF][START_REF] Lemmon | High-Throughput Genomic Data in Systematics and Phylogenetics[END_REF]. For phylogenomics, most prominently featured are sequence-capture, focusing on highly conserved regions (e.g., [START_REF] Faircloth | Ultraconserved Elements Anchor Thousands of Genetic Markers Spanning Multiple Evolutionary Timescales[END_REF][START_REF] Lemmon | Anchored Hybrid Enrichment for Massively High-Throughput Phylogenomics[END_REF]reviewed by Jones & Good, 2016) and transcriptomes (e.g., Misof et al., 2014;[START_REF] Cohen | Comparative transcriptomics across the prokaryotic Tree of Life[END_REF][START_REF] Fernández | Evaluating Topological Conflict in Centipede Phylogeny Using Transcriptomic Data Sets[END_REF][START_REF] Park | Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle[END_REF][START_REF] Simion | A Large and Consistent Phylogenomic Dataset Supports Sponges as the Sister Group to All Other Animals[END_REF][START_REF] Irisarri | Phylotranscriptomic consolidation of the jawed vertebrate timetree[END_REF], but phylogenomic trees have also been constructed based on restriction-digest methods that focus on single nucleotide polymorphisms or SNPs [START_REF] Leaché | Phylogenomics of phrynosomatid lizards: conflicting signals from sequence capture versus restriction site associated DNA sequencing[END_REF][START_REF] Harvey | Sequence Capture versus Restriction Site Associated DNA Sequencing for Shallow Systematics[END_REF] and analysis of transposable elements (e.g., [START_REF] Suh | The Dynamics of Incomplete Lineage Sorting across the Ancient Adaptive Radiation of Neoavian Birds[END_REF]. This diversity of marker types for phylogenetics should be celebrated, but each marker type brings with it a list of pros and cons. For example, many questions in the higher level phylogenetics of animals and plants have so far relied almost exclusively on transcriptome data. However, the uncritical use of transcriptomes in phylogenetics is not without caveats. At high taxonomic levels, coding regions can exhibit extreme levels of among-taxon base composition, sometimes resulting in strong violations of phylogenetic models [START_REF] Romiguier | Phylogenomics Controlling for Base Compositional Bias Reveals a Single Origin of Eusociality in Corbiculate Bees[END_REF][START_REF] Romiguier | Analytical Biases Associated with GC-Content in Molecular Evolution[END_REF]. Coding regions can exhibit reduced levels of incomplete lineage sorting (ILS) compared to noncoding regions [START_REF] Scally | Insights into hominid evolution from the gorilla genome sequence[END_REF]. Such reduced ILS could in fact be helpful in building complex phylogenies with rapid radiations (Edwards, 2009a), but it will certainly distort estimated branch lengths when coalescent methods, which assume neutrality, are used. In addition to yielding abundant sequence variations and SNPs, transcriptome data also yields information on the levels of expression of various genes, and in which tissue-specific genes are expressed. However, using these aspects of transcriptome data are less likely to be of use to phylogeneticists, precisely because specific genes are often tissue-specific and because expression data can exhibit high levels of variation among individuals, populations and species in space and time [START_REF] Todd | The power and promise of RNA-seq in ecology and evolution[END_REF]. Such variation will certainly pose limitations for long term phylogenomic endeavors that will likely combine data collected originally for different purposes.

Although non-coding portions of the genome have been largely neglected in phylogenomics because they are difficult to align and analyze, we are now making progress in understanding their sequence evolution and how it might be informative for comparative purposes [START_REF] Ulitsky | Evolution to the rescue: using comparative genomics to understand long non-coding RNAs[END_REF][START_REF] Edwards | Conserved Nonexonic Elements: A Novel Class of Marker for Phylogenomics[END_REF]. For instance, studies on enhancer and promoter evolution in mammals have shown that despite low levels of sequence conservation, there is conservation of regulatory function and 3D structure across species that carries information for comparative purposes [START_REF] Villar | Enhancer evolution across 20 mammalian species[END_REF][START_REF] Berthelot | Complexity and conservation of regulatory landscapes underlie evolutionary resilience of mammalian gene expression[END_REF]. The development of methods to infer and interpret the evolutionary history and phylogenetic signal of non-coding elements and 3D genome structure is a critical priority.

SNPs have been advocated by some authors for higher level phylogenetics [START_REF] Leaché | The Utility of Single Nucleotide Polymorphism (SNP) Data in Phylogenetics[END_REF], but the available methods for analyzing such data are still extremely limited. For example, concatenation and two coalescent methods (SNAPP and SVD quartets: Bryant et al., 2012;[START_REF] Chifman | Quartet inference from SNP data under the coalescent model[END_REF] have recently been highlighted as the main methods available for phylogenomic analysis of SNPs [START_REF] Leaché | The Utility of Single Nucleotide Polymorphism (SNP) Data in Phylogenetics[END_REF]. But each of these methods has its shortcomings. It is likely that concatenation of SNPs will be misleading for many of the same reasons that concatenating genes can be misleading, due to different gene histories (see section 'Concepts and models in phylogenomics' for further details; [START_REF] Kubatko | Inconsistency of phylogenetic estimates from concatenated data under coalescence[END_REF]. SNAPP, a coalescent method suitable for analysis of SNPs [START_REF] Bryant | Inferring species trees directly from biallelic genetic markers: bypassing gene trees in a full coalescent analysis[END_REF], works well only on relatively small data sets, and it is unclear how well SVD quartets performs on some data sets [START_REF] Shi | Coalescent-Based Analyses of Genomic Sequence Data Provide a Robust Resolution of Phylogenetic Relationships among Major Groups of Gibbons[END_REF]. Although SNPs do provide a helpful route around the often-violated assumption in coalescent models of no recombination within loci [START_REF] Bryant | Inferring species trees directly from biallelic genetic markers: bypassing gene trees in a full coalescent analysis[END_REF], and are widely seen as excellent markers for phylogeography and population genetics, it remains to be seen how powerful they are at higher phylogenetic levels.

Despite the diversity of marker types for phylogenomics, it remains unclear whether features specific to each marker type can ultimately result in phylogenomic datasets that can strongly mislead. For example, incongruence in the phylogeny of modern birds developed by Jarvis et al. (2014;48 whole genomes) and [START_REF] Prum | A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing[END_REF]259 anchored phylogenomics loci, 198 species) has recently been attributed to differences in marker type rather than number of taxa [START_REF] Reddy | Why Do Phylogenomic Data Sets Yield Conflicting Trees? Data Type Influences the Avian Tree of Life more than Taxon Sampling[END_REF]. Whereas Jarvis et al. (2014) used primarily noncoding loci because they observed gross incongruence when using coding regions, the loci used by [START_REF] Prum | A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing[END_REF], although nominally focused on broadly "anchored" conserved regions, came primarily from coding regions. Thus, at least one marker type is likely inappropriate when applied across modern birds (ca. 100 Ma). These data type effects can stem from multiple sources. Selection on exons might lead to localized differences in effective population size across the genome, as previous studies have highlighted issues with base composition heterogeneity within exons across taxa [START_REF] Figuet | Biased Gene Conversion and GC-Content Evolution in the Coding Sequences of Reptiles and Vertebrates[END_REF][START_REF] Scally | Insights into hominid evolution from the gorilla genome sequence[END_REF]. On the other hand, alignment quality of introns and ultraconserved elements (UCEs) can sometimes be less than desired [START_REF] Edwards | Conserved Nonexonic Elements: A Novel Class of Marker for Phylogenomics[END_REF]. Clearly marker effects can potentially have substantial consequences on species tree estimates and need to be further evaluated and compared side-by-side by the phylogenetic community [START_REF] Shen | Contentious relationships in phylogenomic studies can be driven by a handful of genes[END_REF].

A priori versus a posteriori selection of loci for phylogenomics

In an ideal world, phylogeneticists would have whole and fully annotated genomes of all taxa available, allowing them to select loci for phylogenomics based on the relative merits of different loci. This a posteriori (after data generation) selection would be carried out after assessing desired phylogenetic and biological properties of a wide array of markers for which the data are already in hand. A posteriori selection of loci for phylogenomics is clearly a long-term goal that will yield greater choice and justification for specific loci. Today, the loci for phylogenomics are selected a priori (i.e., before data generation) based primarily on cost and ease of collection and alignment, disregarding potentially useful regions of the genome. Thus, an attractive aspect of whole-genome sequencing (WGS) for phylogenomics is to have the opportunity to select markers a posteriori once genomes are in hand (e.g., [START_REF] Edwards | Conserved Nonexonic Elements: A Novel Class of Marker for Phylogenomics[END_REF]Fig. 1). WGS is less prone to sequencing biases and also allows for further expansion into different research fields and questions based on the same initial data [START_REF] Lelieveld | Comparison of Exome and Genome Sequencing Technologies for the Complete Capture of Protein-Coding Regions[END_REF]. In contrast, a priori marker selection often limits the kinds of questions and methods that researchers can apply and represents a real constraint for phylogenomics and other disciplines.

An important constraint for using WGS for downstream phylogenomic analyses is genome quality. Obtaining high-coverage well-assembled and thoroughly annotated genomes is still very expensive and time-consuming, and even low coverage genomes are still outside reach for large portions of the scientific community. However, even low coverage genomes, which should be cautiously used for maker selection due to potential problems caused by poor annotation and coverage, can sometimes yield a modest number of markers for phylogenomics, and in the short term might even yield data sets allowing a broader diversity of markers for analysis. Although we are fully aware of its constraints, we are particularly excited about the potential that we see in routinely using WGS to produce phylogenomic data sets.

More taxa versus more loci

The question of whether to add more genes or more taxa was a dominant theme in phylogenetics in the 1990s and early 2000s (e.g., [START_REF] Hillis | Inferring complex phylogenies[END_REF][START_REF] Kim | General inconsistency conditions for maximum parsimony: effects of branch lengths and increasing numbers of taxa[END_REF], and remains a persistent question in phylogenomics today. After much debate in the literature (e.g., [START_REF] Hillis | Inferring complex phylogenies[END_REF][START_REF] Graybeal | Is it better to add taxa or characters to a difficult phylogenetic problem?[END_REF][START_REF] Hillis | Taxonomic sampling, phylogenetic accuracy, and investigator bias[END_REF][START_REF] Poe | Sensitivity of phylogeny estimation to taxonomic sampling[END_REF][START_REF] Mitchell | More taxa or more characters revisited: Combining data from nuclear protein-encoding genes for phylogenetic analyses of Noctuoidea (Insecta: Lepidoptera)[END_REF], the initial consensus view from the Sanger sequencing era of phylogenetics, is that adding more taxa generally improves phylogenetic analysis more so than more markers (e.g., [START_REF] Hillis | Inferring complex phylogenies[END_REF][START_REF] Graybeal | Is it better to add taxa or characters to a difficult phylogenetic problem?[END_REF][START_REF] Poe | Sensitivity of phylogeny estimation to taxonomic sampling[END_REF]. However, phylogenomics is adding a new twist to this consensus, both from the standpoint of data acquisition and from theory (e.g., [START_REF] Rokas | More Genes or More Taxa? The Relative Contribution of Gene Number and Taxon Number to Phylogenetic Accuracy[END_REF][START_REF] Nabhan | The impact of taxon sampling on phylogenetic inference: a review of two decades of controversy[END_REF][START_REF] Xi | Phylogenomics and a posteriori data partitioning resolve the Cretaceous angiosperm radiation Malpighiales[END_REF][START_REF] Patel | Error in Phylogenetic Estimation for Bushes in the Tree of Life[END_REF]. Amassing large data sets, both in terms of more taxa and more loci, is still a guiding principle of phylogenomics. But with the ability now to bring together many different types of markers in a single analysis, and to analyze them in ways that were not previously available, the "more taxa vs. more genes" debate is becoming more nuanced [START_REF] Nabhan | The impact of taxon sampling on phylogenetic inference: a review of two decades of controversy[END_REF]. For example, recent work shows that this debate can be highly context-specific and modeldependent (e.g., [START_REF] Baurain | Lack of Resolution in the Animal Phylogeny: Closely Spaced Cladogeneses or Undetected Systematic Errors?[END_REF][START_REF] Dell Ampio | Decisive Data Sets in Phylogenomics: Lessons from Studies on the Phylogenetic Relationships of Primarily Wingless Insects[END_REF]Edwards et al., 2016b). Also, some phylogenetic methods, such as coalescent methods, appear to be more robust to limited taxon sampling than traditional methods like concatenation [START_REF] Song | Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model[END_REF]Liu, Xi & Davis, 2015). Some researchers favor "horizontal" data matrices, wherein the number of loci far exceeds the number of taxa, whereas other researchers favor "vertical" matrices, where many taxa are analyzed at just a few (1-5) loci.

Whereas the PCR era of phylogenetics was often dominated by vertical matrices, HTS is allowing data matrices to become more horizontal (Fig. 2). It is important to note that as these more horizontal data become more prevalent, they increase the amount of missing data and aligning problems that can contribute to misleading or low phylogenetic resolution. At least in a coalescent framework, scaling up in both dimensions will be crucial for improved phylogenies and the number of loci required to resolve a given phylogenetic problem is often a function of the coalescent branch lengths in the phylogenetic tree being resolved, with longer branches requiring fewer loci [START_REF] Edwards | High-resolution species trees without concatenation[END_REF][START_REF] Huang | Sources of error inherent in species-tree estimation: impact of mutational and coalescent effects on accuracy and implications for choosing among different methods[END_REF]. At deeper time scales, increasing the number of loci have not proven particularly useful to resolve problematic areas in the ToL (e.g.sister group to animals; [START_REF] King | Embracing Uncertainty in Reconstructing Early Animal Evolution[END_REF]. Despite methodological complications in the accurate estimation of population genetic parameters and computational limitations for MCMC convergence at deep times in coalescent-based analyses, recalcitrant nodes likely represent true intrisicancies of the diversification history of these groups and not necessarily reflect failures of coalescent-based phylogenetics [START_REF] Lanier | Applying species-tree analyses to deep phylogenetic histories: Challenges and potential suggested from a survey of empirical phylogenetic studies[END_REF].

To study how researchers have resolved challenges of balancing numbers of taxa versus numbers of loci, we quantified trends in phylogenomic data set size and structure over the past 13 years, drawing data from 164 data sets across diverse taxa (Supplementary Table S1). We found that, whereas the number of species per paper has not increased significantly over time (Fig. 2A), there were significant increases with time in number of loci (Fig. 2B), total length of sequence analyzed (Fig. 2C), as well as total data set size, as measured by the product of species times locus number (Fig. 2E) or species times total alignment length (Fig. 2F). These results mirror similar trends evaluated for the size of data sets in phylogeography [START_REF] Garrick | The evolution of phylogeographic data sets[END_REF]. Surprisingly, we found no evidence for a tradeoff between the number of species investigated and the number of loci analyzed (Fig. 2D); perhaps HTS data sets have plateaued somewhat in terms of number of loci, whereas the number of species analyzed is more a function of the questions being asked and the clade being investigated. Regardless, we suspect that, in general, the number of loci and total alignment lengths in phylogenomic data sets are likely a function of resources and sequencing effort. The era of whole genome sequencing in phylogenomics is still dawning, given that most studies thus far have used targeted approaches for sampling loci (Supplementary Table S1). We suspect that once whole genome sequencing on a clade-wide basis become routine (e.g., Genome 10K Community of Scientists, 2009;[START_REF] Grigoriev | MycoCosm portal: gearing up for 1000 fungal genomes[END_REF][START_REF] Cheng | 10KP: A phylodiverse genome sequencing plan[END_REF], we will witness yet another jump in the sizes of phylogenomic data sets.

Filtering heterogeneous phylogenomic data sets

Recent studies show that the addition of more loci and more taxa can result in higher levels of gene-tree discordance (e.g., [START_REF] Smith | Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants[END_REF][START_REF] Shen | Contentious relationships in phylogenomic studies can be driven by a handful of genes[END_REF]. This is not unexpected -as the number of taxa and loci increase, the greater the likelihood that the dataset will capture the heterogeneous evolutionary history (e.g., incomplete lineage sorting [ILS], lateral gene transfer [LGT], hybridization, gene duplication and loss [GDL]; See Table 1 for a definition of these concepts), misidentification of orthologs from paralogs, and patterns of molecular evolution (e.g., noise/lack of signal in the sequences, and nonstationarity in base composition) that can contribute to gene tree discord. At the same time, the variance in gene tree topologies could also have been caused by errors in gene tree estimation. Such observations have been used to argue that the accuracy of gene tree inference should be maximized or at least evaluated, but it is not clear what criteria should be used to filter sets of gene trees. For example, filters can be based on rates of molecular evolution [START_REF] Klopfstein | More on the Best Evolutionary Rate for Phylogenetic Analysis[END_REF], levels of phylogenetic informativeness [START_REF] Fong | A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia[END_REF], or on the cause of gene-tree discord itself, if known [START_REF] Huang | Sources of error inherent in species-tree estimation: impact of mutational and coalescent effects on accuracy and implications for choosing among different methods[END_REF]. [START_REF] Chen | Selecting Question-Specific Genes to Reduce Incongruence in Phylogenomics: A Case Study of Jawed Vertebrate Backbone Phylogeny[END_REF] found that selecting genes whose trees contained a well-known uncontested long branch in a given species phylogeny was a better way to improve phylogenomic signal than selecting genes based on characteristics of sequence evolution. However, the effects of such culling on the distribution of gene trees, and whether it could distort the distribution so that it no longer conforms to models like the multispecies coalescent, are unknown, and potentially of concern (but see [START_REF] Huang | Cause of gene tree discord? Distinguishing incomplete lineage sorting and lateral gene transfer in phylogenetics[END_REF]. Aside from the use of some explicit methods for detection of outlier genes (e.g., de Vienne, Ollier & Aguileta, 2012), rogue taxa (e.g., [START_REF] Aberer | Pruning Rogue Taxa Improves Phylogenetic Accuracy: An Efficient Algorithm and Webservice[END_REF], outlier long branches (e.g., [START_REF] Mai | TreeShrink: fast and accurate detection of outlier long branches in collections of phylogenetic trees[END_REF], and tree space visualization (e.g., [START_REF] Huang | TreeScaper: Visualizing and Extracting Phylogenetic Signal from Sets of Trees[END_REF][START_REF] Jombart | TREESPACE: Statistical exploration of landscapes of phylogenetic trees[END_REF], an obvious way to alleviate potential effects of gene tree outliers is a more balanced taxon sampling [START_REF] Hedtke | Resolution of Phylogenetic Conflict in Large Data Sets by Increased Taxon Sampling[END_REF]. Nonetheless, we need further studies on the effects of different types of phylogenomic filters on the properties of large-scale phylogenomic datasets.

High throughput sequencing opens possibilities for new information and marker types

Heterozygosity and Intra-Individual Site Polymorphisms

Organisms with multiple ploidy levels are intra-individual polymorphisms and heterozygosity is, of course, common in diploid and polyploid organisms. However, due to issues such as lack of sufficient read coverage and connectivity, confident identification of such polymorphism continues to be challenging (Garrick, Sunnnucks & Dyer, 2010;[START_REF] Lischer | Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles[END_REF]Schrempff et al., 2016) and many data sets do not permit statistical approaches, such as PHASE [START_REF] Stephens | A new statistical method for haplotype reconstruction from population data[END_REF] to robustly determine haplotypes of different alleles (Garrick, Sunnnucks & Dyer, 2010).

Consequently, in phylogenetics, heterozygosity and intra-specific polymorphic sites are either accommodated using UIPAC ambiguity codes or ignored entirely or by randomly selecting alleles [START_REF] Iqbal | De novo assembly and genotyping of variants using colored de Bruijn graphs[END_REF]. In fact, most "one allele per individual/species" phylogenomic data sets consists of haplotypes that might not occur in nature because many methods yield single haplotypes consisting of consensus or other haplotype summaries from diploid organisms. The fact that HTS produces several reads of the same region allows for the identification of heterozygosity and intra-specific polymorphic sites represents an untapped opportunity to incorporate intra-individual variation in our phylogenetic estimates [START_REF] Lischer | Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles[END_REF]Schrempff et al., 2016;[START_REF] Andermann | Allele Phasing Greatly Improves the Phylogenetic Utility of Ultraconserved Elements[END_REF]. Recent models have been proposed to improve calling and sorting such polymorphisms [START_REF] De Maio | Linking Great Apes Genome Evolution across Time Scales Using Polymorphism-Aware Phylogenetic Models[END_REF][START_REF] Lischer | Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles[END_REF][START_REF] Potts | Constructing Phylogenies in the Presence Of Intra-Individual Site Polymorphisms (2ISPs) with a Focus on the Nuclear Ribosomal Cistron[END_REF]Schrempff et al., 2016) and, although results of different studies vary [START_REF] Kubatko | Inferring species-level phylogenies and taxonomic distinctiveness using multilocus data in Sistrurus rattlesnakes[END_REF][START_REF] Lischer | Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles[END_REF], estimation of individual, naturally occurring haplotypes has been shown to improve phylogenomic reconstructions based on genome-scale data [START_REF] Andermann | Allele Phasing Greatly Improves the Phylogenetic Utility of Ultraconserved Elements[END_REF].

Rare genomic changes

As noted above, molecular phylogenetics has primarily used alignments of sequence-level data for phylogenetic inference. This bias is perhaps driven by the notion that genome evolution occurs by aggregating small changes, such as point substitutions, over time, but more likely it is due to the challenges of characterizing rare genomic changes, such as indels, transpositions, inversions, and other large-scale genomic events [START_REF] Rokas | Rare genomic changes as a tool for phylogenetics[END_REF][START_REF] Boore | The use of genome-level characters for phylogenetic reconstruction[END_REF][START_REF] Bleidorn | Sources of Error and Incongruence in Phylogenomic Analyses[END_REF]. This emphasis on sequence data has produced a vast ecosystem of algorithms tailored to analyze such data, but most phylogeneticists would agree that rare genomic changes would be a welcome addition to the toolkit of phylogenomics, because they are generally regarded as highly informative markers, providing strong evidence of homology and monophyly [START_REF] Boore | The use of genome-level characters for phylogenetic reconstruction[END_REF][START_REF] Rogozin | Homoplasy in genome-wide analysis of rare amino acid replacements: the molecular-evolutionary basis for Vavilov's law of homologous series[END_REF]. With the increased availability and affordability of WGS, our view of genome plasticity has changed drastically in recent years and we are now capable of exploring other genomic features beyond the signals encapsulated in DNA or amino acid sequences (e.g., [START_REF] Ryan | The genome of the ctenophore Mnemiopsis leidyi and its implications for cell type evolution[END_REF]. The question then arises of how to identify and utilize these rare genomic markers, as well and assess their phylogenetic informativeness [START_REF] Rokas | Rare genomic changes as a tool for phylogenetics[END_REF][START_REF] King | Embracing Uncertainty in Reconstructing Early Animal Evolution[END_REF]. Genome-level characters will likely have different evolutionary properties than sequencebased markers, suggesting that one of the biggest challenges we face for incorporating genomic changes into phylogenetic analyses is to find informative evolutionary models and tools suited for these kinds of data and assess how congruent or discordant they are with respect to other markers (e.g., [START_REF] Rota-Stabelli | A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata[END_REF]. This will not only shed light on how phylogenetically informative different genomic changes are, but also will broaden our understanding of the evolutionary intricacies across different genomic regions [START_REF] Rokas | Phylogenetic Analysis of Protein Sequence Data Using the Randomized Axelerated Maximum Likelihood (RAXML) Program[END_REF][START_REF] Leigh | Evaluating Phylogenetic Congruence in the Post-Genomic Era[END_REF].

Gene order and synteny

Computational algorithms to use gene order and rearrangements as markers in phylogenetics [START_REF] Tang | Phylogenetic reconstruction from arbitrary gene-order data[END_REF][START_REF] Ghiurcuta | Evaluating synteny for improved comparative studies[END_REF][START_REF] Kowada | New Genome Similarity Measures Based on Conserved Gene Adjacencies[END_REF] were spurred in part by the seminal paper by [START_REF] Boore | Complete sequence, gene arrangement, and genetic code of mitochondrial DNA of the cephalochordate Branchiostoma floridae (Amphioxus)[END_REF] using mitochondrial gene rearrangements to understand the phylogeny of arthropods. Initially, algorithms for making use of gene order and synteny were applied primarily to microbial genomes, but recent efforts have extended such methods to the analysis of eukaryotes as well (see [START_REF] Lin | Maximum likelihood phylogenetic reconstruction from highresolution whole-genome data and a tree of 68 eukaryotes[END_REF]. Gene order and synteny appear most promising at high phylogenetic levels, although we still do not know how informative gene order will be at many levels, such as within mammals. Chromosomal rearrangements appear highly dynamic in some groups, such as mammals, and further study of their use in phylogenomics is warranted [START_REF] Murphy | Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps[END_REF].

Indels and transpositions

Indels and transpositions are two types of molecular characters that are underutilized in phylogenomics, the former perhaps because standard methods of analysis often treat indels as missing data and the latter because they are technically challenging to collect without whole genome data. Indels have been used sporadically in phylogenomics and several have argued for their utility and informativeness, given appropriate analytical tools (Jarvis et al., 2014;[START_REF] Ashkenazy | Indel Reliability in Indel-Based Phylogenetic Inference[END_REF][START_REF] Roncal | Active transposable elements recover species boundaries and geographic structure in Madagascan coffee species[END_REF]. [START_REF] Murphy | Using genomic data to unravel the root of the placental mammal phylogeny[END_REF] used indels in protein-coding regions to bolster estimates of mammalian phylogeny and found that the Atlantogenata hypothesis was supported after scrutinizing proteome-wide indels for spurious alignments and orthology. The Avian Phylogenomics Project found that indels had less homoplasy than SNPs and, despite showing high levels of ILS, was largely congruent with other markers across the avian tree. Transposable elements arguably are even more highly favored by phylogenomics researchers, but are much more difficult to isolate and analyze and have been used principally across various studies in mammals and birds [START_REF] Kaiser | Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds[END_REF][START_REF] Churakov | Rodent Evolution: Back to the Root[END_REF][START_REF] Kriegs | Retroposon Insertions Provide Insights into Deep Lagomorph Evolution[END_REF][START_REF] Suh | Mesozoic retroposons reveal parrots as the closest living relatives of passerine birds[END_REF][START_REF] Baker | Genomic Support for a Moa-Tinamou Clade and Adaptive Morphological Convergence in Flightless Ratites[END_REF]. Whereas they are generally considered to have a low rate of homoplasy, most researchers agree that they can in some circumstances exhibit insertional homoplasy. Moreover, no marker is immune to the challenges of ILS, and transposable elements and indels are no exception [START_REF] Matzke | Retroposon Insertion Patterns of Neoavian Birds: Strong Evidence for an Extensive Incomplete Lineage Sorting Era[END_REF][START_REF] Suh | The Dynamics of Incomplete Lineage Sorting across the Ancient Adaptive Radiation of Neoavian Birds[END_REF]. Still, the exceptional resolution afforded by some studies employing transposable elements is exciting, and we expect this marker type to increase in use as whole genomes are collected with higher frequency.

Copy Number Variations (CNV)

The 1000 Genomes Project estimates that in humans about 20 million base pairs are affected by structural variants, including copy number variations (CNV) and large deletions (1000 Genomes Project Consortium, 2015), suggesting that these types of mutations encompass a higher fraction of the genome than do SNPs in humans. A CNV is a DNA segment of at least one kilobase (kb) that varies in copy number compared with a reference genome [START_REF] Redon | Global variation in copy number in the human genome[END_REF]. CNVs appear as deletions, insertions, duplications, and complex multi-site variants [START_REF] Fredman | Complex SNPrelated sequence variation in segmental genome duplications[END_REF]. Such a profusion of CNVs across human genomes has proven useful in tracking population structure [START_REF] Sjödin | Population genetic nature of copy number variation. Population Genetic Nature of Copy Number Variation[END_REF], but still remains underappreciated in phylogenetics.

Newly available methods allow inference of CNV at high resolution with great accuracy [START_REF] Wiedenhoeft | Fast Bayesian Inference of Copy Number Variants using Hidden Markov Models with Wavelet Compression[END_REF]. The frequency with which CNVs occur in animal and plant populations raises the question of how informative they would be at higher phylogenetic levels, and whether they would incur unwanted homoplasy that would obscure homology and phylogenetic relationships. For example, some CNVs evolve so quickly that they can be used with success at the subindividual level, for example, in tracking clonal evolution of cancer cells using CNV-specific phylogenetic methods (e.g., [START_REF] Schwartz | The evolution of tumour phylogenetics: Principles and practice[END_REF][START_REF] Liang | An Improved Binary Differential Evolution Algorithm to Infer Tumor Phylogenetic Trees[END_REF][START_REF] Ricketts | Using LICHeE and BAMSE for Reconstructing Cancer Phylogenetic Trees[END_REF][START_REF] Urrutia | Integrative pipeline for profiling DNA copy number and inferring tumor phylogeny[END_REF]. Moreover, their interspecific variation has been shown to correlate with the phylogeny of some groups, such as the highly pathogenic and rapidly-evolving barley powdery mildews [START_REF] Frantzeskakis | Signatures of host specialization and a recent transposable element burst in the dynamic one-speed genome of the fungal barley powdery 52 mildew pathogen[END_REF]. However, such fast evolution may mean that these markers might be less useful at higher levels of biological organization. Additionally, the adaptive nature of CNVs may or may not facilitate clear phylogenetic signals. For example, a study in Arabidopsis thaliana [START_REF] Debolt | Copy Number Variation Shapes Genome Diversity in Arabidopsis Over Immediate Family Generational Scales[END_REF] showed that adaptation to novel environments, or to varying temperatures, is associated with mutations in CNVs. If CNVs are to become a useful tool in phylogenomics or phylogeography, we must understand their microevolutionary properties in greater detail. For example, the pattern of evolution of CNVs, wherein deletions of genetic material may not easily revert, resulting in a type of Dollo evolution, might help clarify the overall structure of the models applied to them [START_REF] Rogozin | Dollo parsimony and the reconstruction of genome evolution[END_REF][START_REF] Gusfield | Persistent phylogeny[END_REF].

Recent advances in the generation of high throughput sequence data and their impact on the reconstruction of the Tree of Life

As sequencing technology rapidly moves forward (reviewed by [START_REF] Goodwin | Coming of age: ten years of next-generation sequencing technologies[END_REF], our ability to accurately identify the aforementioned marker types increases considerably. For instance, the low per-base error rate of short-read sequencing technologies, such as the Illumina HiSeq X Ten and NovaSeq (Illumina, San Diego, CA), allow for a significant reduction in the cost of sequencing which can result in data for more taxa at a higher coverage. This is certainly beneficial for the accurate identification of SNPs, heterozygosity, and intra-individual site polymorphisms [START_REF] Goodwin | Coming of age: ten years of next-generation sequencing technologies[END_REF] and Oxford Nanopore (Oxford Nanopore Technologies, ONT) produce reads that exceed 10 kb [START_REF] Rhoads | PacBio Sequencing and Its Applications[END_REF][START_REF] Lu | Oxford Nanopore MinION Sequencing and Genome Assembly[END_REF]. The advent of these technologies has led to improved and more efficient assembly methods that allow accurate identification of structural changes (e.g., [START_REF] Khost | Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster[END_REF][START_REF] Merker | Long-read genome sequencing identifies causal structural variation in a Mendelian disease[END_REF]. When combined with data resulting from short-read sequencing, they represent a powerful tool to correctly identify and use a wide array of genomic markers for numerous purposes. These technical advances, which include portable devices that can be carried into the field (i.e., Oxford Nanopore; [START_REF] Johnson | Real-Time DNA Sequencing in the Antarctic Dry Valleys Using the Oxford Nanopore Sequencer[END_REF], will certainly yield an increase in the genomic loci and taxa available for genomic and phylogenomic studies.

Concepts and models in phylogenomics

For decades, phylogenetics has struggled with how best to translate evolutionary changes in DNA sequences and other characters into phylogenies, and genomic data are no exception to this trend.

Phylogenomics is still in a developing stage of formulating models that effectively represent the underlying mechanisms for genome-scale variation while remaining efficient and within reasonable analytical and bioinformatic capacities. The current focus on models and evolutionary forces generating the patterns that we recover as branching and reticulation events in our phylogenetic reconstructions is a healthy one and can be extended to other important topics in phylogenomics, such as species delimitation, character mapping, and trait evolution (e.g., [START_REF] Yang | Unguided species delimitation using DNA sequence data from multiple loci[END_REF]. All of these areas are developing rapidly and are in need of updated models and bioinformatics applications to cope with the heterogeneity brought by genome-scale data.

The multispecies coalescent model

One of the key practical advances in molecular phylogenetics has been the incorporation of gene tree stochasticity into the inference of species phylogenies, via the multispecies coalescent model (MSC: [START_REF] Rannala | Bayes estimations of species divergence times and ancestral population sizes using DNA sequences from multiple loci[END_REF]Liu & Pearl, 2007;[START_REF] Heled | Bayesian inference of species trees from multilocus data[END_REF]. The MSC allows gene trees to be inferred with their own histories, including coalescent-appropriate branching models, but contained within independent yet connected lineages within a species phylogeny, with speciation-appropriate branching models [START_REF] Degnan | Gene tree discordance, phylogenetic inference and the multispecies coalescent[END_REF]. The main conceptual advance has been to understand and separately manage the variation at different levels of biological organization -an advance that began years ago [START_REF] Doyle | Gene trees and species trees: molecular systematics as one-character taxonomy[END_REF][START_REF] Maddison | Gene trees in species trees[END_REF][START_REF] Pamilo | Relationships between gene trees and species trees[END_REF]), but has only recently been widely embraced and put into practice (Edwards, 2009b). Given its ability to accommodate heterogeneous histories across loci scattered throughout the genome, the MSC lays at the core of the conceptual framework to deal with genome-scale data (e.g., [START_REF] Rannala | Phylogenetic inference using whole genomes[END_REF]Liu et al., 2015). In the few instances in which model comparison and fit has been evaluated (Liu & Pearl 2007;[START_REF] Edwards | High-resolution species trees without concatenation[END_REF]), the MSC vastly outperforms concatenation. This of course does not mean that the MSC is the correct, or even an adequate, model for phylogenomic data [START_REF] Reid | Poor fit to the multispecies coalescent is widely detectable in empirical data[END_REF]. Despite concerns regarding some of its implementations when dealing with genomic data (e.g., [START_REF] Springer | The gene tree delusion[END_REF], the MSC is a powerful theoretical model for phylogenomics and there is room for refinement and improvement for its applications (e.g., Edwards et al., 2016b[START_REF] Xu | Challenges in Species Tree Estimation Under the Multispecies Coalescent Model[END_REF].

Bypassing full likelihood models by relying on summaries of the coalescent process

Given the huge computational resources required for modelling all the complexities of evolutionary processes in a statistical framework, there is interest in methods that will accommodate genome-scale data for large numbers of species within a coalescent framework. The utility of such methods cannot be overstated: the rapid rise of large-scale genomic data sets has clearly outstripped theoretical and computational methods required to analyze them. For example, although progress is being made regarding scalability of full Bayesian methods of species phylogeny inference (e.g., [START_REF] Ogilvie | StarBEAST2 Brings Faster Species Tree Inference and Accurate Estimates of Substitution Rates[END_REF], they are still unable to accommodate large phylogenomic datasets, which often consist of hundreds of species for thousands of loci (Supplementary Table S1). A common approach to speeding up species phylogeny inference consists of 'two-step' methods, wherein gene trees are estimated first and separately from the species phylogeny; then, using various summaries of the coalescent process for collections of gene trees, a species phylogeny is estimated. Many useful methods for estimating species phylogenies in this way have been proposed (see [START_REF] Marcussen | Ancient hybridizations among the ancestral genomes of bread wheat[END_REF][START_REF] Liu | Coalescent methods for estimating species trees from phylogenomic data[END_REF][START_REF] Mirarab | ASTRAL-II: coalescent-based species tree estimation with many hundreds of taxa and thousands of genes[END_REF][START_REF] Mirarab | Evaluating summary methods for multilocus species tree estimation in the presence of incomplete lineage sorting[END_REF], taking advantage of summary statistics of the coalescent process, such as the average ranks of pairs of species in the collection of gene trees (e.g., STAR: [START_REF] Liu | Estimating species phylogenies using coalescence times among sequences[END_REF]ASTRAL-II: Mirarab & Warnow, 2015) or the distribution of gene trees containing triplets of species (e.g., MP-EST; [START_REF] Liu | A maximum pseudo-likelihood approach for estimating species trees under the coalescent model[END_REF]. Some of these two-step methods, while approximate, nonetheless allow for statistical testing in a likelihood framework. For example, MP-EST can evaluate the (pseudo)likelihood of two proposed species phylogenies given a collection of gene trees and the difference in likelihood can be used to evaluate two proposed species phylogenies against each other. However, such statistical approaches have rarely been used thus far, and bootstrapping or approximate posterior probabilities on branches are by far the most common statistics applied to species phylogenies [START_REF] Sayyari | Fast Coalescent-Based Computation of Local Branch Support from Quartet Frequencies[END_REF]. Speeding up the estimation process using twostep methods can be effective, but it can also accumulate errors or misallocate sources of variance which cannot be corrected at later stages [START_REF] Xu | Challenges in Species Tree Estimation Under the Multispecies Coalescent Model[END_REF]. If gene trees are biased or uninformative, then downstream analyses for species phylogeny estimation or species delimitation may similarly be compromised (e.g., [START_REF] Olave | Upstream analyses create problems with DNA-based species delimitation[END_REF]. For example, MP-EST can sometimes perform poorly when PhyML [START_REF] Guindon | New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0[END_REF] is used to build low-information gene trees because PhyML may produce biased gene trees when the alignments contain very similar sequences (Xi, Liu & Davis, 2015). This may account for the lower performance of MP-EST compared to ASTRAL in some simulation conditions, because ASTRAL resolves input polytomies and zero-length branches in gene trees more appropriately. This difference between MP-EST and ASTRAL is eliminated when RAxML [START_REF] Stamatakis | RAxML Version 8: A tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies[END_REF] is used to build gene trees (Xi, Liu & Davis, 2015).

Beyond the multispecies coalescent model

Reticulation at multiple levels challenges the standard multispecies coalescent model

The phylogenetic processes of branching and reticulation can operate at several levels of organization, including within genes, within genomes, and within populations or species (Figs. 3, 4 and5). For example, recombination can cause reticulations within genes, allopolyploidization can cause reticulations at the level of whole genomes, and introgression and hybridization can cause reticulations at the level of populations. These levels are nested so that branching processes (and in part reticulations) acting at a higher level will cause correlated branching patterns at lower levels. At the same time, reticulations at lower levels, such as recombination acting within genes, will cause inference problems at higher levels, such as estimating population histories. Crucially, however, it is only recombination that will break one key element driving many recent models of phylogenetics and population histories, namely dichotomous gene trees. Reticulations at levels of organization higher than the genome, such as the fusing of populations, as well as gene duplication, will still yield collections of dichotomous gene trees, even if the higher-level history is reticulated. Ultimately, the additive effects of these reticulate processes result in our observed phylogenetic reconstructions, and we expect all of these scenarios to produce bifurcating, dichotomous gene trees. From a modelling point of view, another key distinction is whether at the species level, we still have a phylogeny that is tree-like, or whether a network is needed. The process whereby two populations jointly produce a third requires a network to model properly. Allopolyploidy is another situation requiring a network. There are several statistical methods for inferring homoploid networks [START_REF] Yu | Maximum likelihood inference of reticulate evolutionary histories[END_REF][START_REF] Solís-Lemus | Inferring Phylogenetic Networks with Maximum Pseudolikelihood under Incomplete Lineage Sorting[END_REF][START_REF] Wen | Reticulate evolutionary history and extensive introgression in mosquito species revealed by phylogenetic network analysis[END_REF][START_REF] Wen | Co-estimating Reticulate Phylogenies and Gene Trees from Multi-locus Sequence Data[END_REF], species histories under allopolyploidy [START_REF] Jones | Statistical Inference of Allopolyploid Species Networks in the Presence of Incomplete Lineage Sorting[END_REF], and some two-step methods such as PADRE [START_REF] Huber | Reconstructing the evolutionary history of polyploids from multilabeled trees[END_REF][START_REF] Lott | PADRE: A package for analyzing and displaying reticulate evolution[END_REF]. In general, dealing with multiple simultaneous violations of the MSC, such as introgression and allopolyploidy, remains challenging [START_REF] Degnan | Modeling Hybridization Under the Network Multispecies Coalescent[END_REF]. It is likely that the history of many radiations involves parts of the genome with a dichotomous history and parts that exhibit reticulation, demanding methods that accommodate both scenarios. Alternatively, rather than trying to accommodate multiple processes in our methods for phylogenetic inference, we might instead focus our attention on subsets of loci that would not violate the MSC (e.g., [START_REF] Knowles | A matter of phylogenetic scale: distinguishing incomplete lineage sorting from lateral gene transfer as the cause of gene tree discord in recent versus deep diversification histories[END_REF]. In cases where processes other than incomplete lineage sorting contribute to gene tree discord (i.e., the distribution of trees is statistically inconsistent with expectations under the MSC; see [START_REF] Smith | Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants[END_REF], loci consistent with the MSC model can be identified (e.g., separated from loci with horizontal gene transfer) using CLASSIPHY [START_REF] Huang | Cause of gene tree discord? Distinguishing incomplete lineage sorting and lateral gene transfer in phylogenetics[END_REF]. It has also been suggested that in order to distinguish violations of the MSC, the MSC can be used as a null model to be compared with increasingly complex models that would invoke processes such as hybridization and recombination using networks [START_REF] Degnan | Modeling Hybridization Under the Network Multispecies Coalescent[END_REF]. To follow this promising approach, further research must be conducted to not only model specific processes, but also distinguish them.

Models accommodating dichotomous divergence with gene flow are somewhat limited. For example, in IMa2 [START_REF] Hey | Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis[END_REF][START_REF] Hey | Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics[END_REF][START_REF] Hey | Isolation with Migration Models for More Than Two Populations[END_REF] the species phylogeny must be known and fairly small; in the method of [START_REF] Dalquen | Maximum Likelihood Implementation of an Isolation-with-Migration Model for Three Species[END_REF], both the species phylogeny and gene trees are restricted to three tips. Looking forward, it may be useful to deal with two subproblems: The first is estimating the species phylogeny despite migration, for example by identifying which loci are interfering with the species phylogeny inference or causing reticulations in the form of gene flow. The second sub-problem is to incorporate a gradual speciation process (Fig. 6), where gene flow after speciation slowly declines, perhaps according to some simple function like an exponential. Such a model would capture what is thought to be a more common speciation process than the instantaneous process modelled by the MSC [START_REF] Jones | Divergence estimation in the presence of incomplete lineage sorting and migration[END_REF].

In some cases, it is possible to model one situation with either a species network or a tree with gene flow. [START_REF] Long | The genetic structure of admixed populations[END_REF] discussed two models of admixture: Intermixture and gene flow (Fig. 7). The phylogenetics community has mainly focused on methods for inference under the intermixture model (e.g., the multispecies network coalescent; [START_REF] Yu | Maximum likelihood inference of reticulate evolutionary histories[END_REF], whereas the population genetics community has focused more on models including gene flow (e.g., IM, admixture graphs, G-PhoCS, Phrapl). While some initial work to test inference based on one of these models on data generated by the other has recently appeared [START_REF] Wen | Co-estimating Reticulate Phylogenies and Gene Trees from Multi-locus Sequence Data[END_REF][START_REF] Solís-Lemus | PhyloNetworks: A Package for Phylogenetic Networks[END_REF]Blischak et al., 2018;[START_REF] Zhang | Bayesian Inference of Species Networks from Multilocus Sequence Data[END_REF], much more work is needed to bring together these two lines of work. Simulations and comparisons of observed and expected summary statistics, such as the site-frequency spectrum [START_REF] Excoffier | Robust demographic inference from genomic and SNP data[END_REF], have proven especially useful in distinguishing such scenarios (Fig. 7).

Reticulation in the form of gene flow or introgression is probably the most difficult violation of the MSC to address (but see [START_REF] Hibbins | The timing and direction of introgression under the multispecies network coalescent[END_REF] for a model that estimates the timing and direction of introgression based on the multispecies network coalescent), in part because the number of potential trees accommodating a reticulating network is even higher than the already high number of trees for a given number of taxa. There is at least one issue where reticulation presents an opportunity as well as a challenge. Any kind of gene flow/hybridization means that there is the possibility of inferring the existence of extinct species, because extinct species contribute novel alleles that exceed the coalescence time of most alleles in the focal species under study [START_REF] Hammer | Genetic evidence for archaic admixture in Africa[END_REF]. Well-known examples are the documented presence of Neanderthal genes in most human genomes due to introgression (e.g., [START_REF] Meyer | A high-coverage genome sequence from an archaic Denisovan individual[END_REF] and the presence of genomes derived from now-extinct diploids in extant allopolyploids (i.e. meso-allopolyploids; e.g., [START_REF] Mandáková | Fast Diploidization in Close Mesopolyploid Relatives of Arabidopsis[END_REF][START_REF] Marcussen | From gene trees to a dated allopolyploid network: Insights from the angiosperm genus Viola (Violaceae)[END_REF]. Some current models can explain the data as containing genetic information from extinct species, but they do not model the full species phylogeny: such a generalized approach seems a promising avenue to explore.

Polyploidy and the challenges of analyzing gene duplication and loss

The MSC model describes well allelic lineages and the mutations they accumulate (Fig. 3; [START_REF] Degnan | Gene tree discordance, phylogenetic inference and the multispecies coalescent[END_REF]Liu, Xi & Davis, 2015). The simple MSC model is challenging to apply to evolutionary events in which the evolving entities (genes or paralogs) duplicate and occasionally go extinct during the evolutionary history of the populations/species and thus cannot be sampled in contemporary population or species. Estimating the existence and number of these "ghost" lineages remains challenging. For example, how can we detect duplication events if one of the duplicated loci is lost in descendant lineages? In the case of polyploidy, two (or more) genomes having separate evolutionary histories end up together in a single individual. What consequences for evolutionary history do genomic conflicts and dosage variation in gene expression impose? Polyploidy also raises technical issues, such as whether or not homoeologous sequences are recovered in standard genomic surveys.

The complication that gene duplication and loss (GDL) brings to the inference of species phylogenies has long been recognized [START_REF] Fitch | Distinguishing homologous from analogous proteins[END_REF]. It is therefore surprising that practical solutions to the problem of GDL are almost non-existent, with empirical examples usually based on ad hoc methods and deductions. Ancient duplications where most additional copies are retained in descendent species can be fairly easy to diagnose based on phylogeny [START_REF] Oxelman | RPB2 Gene Phylogeny in Flowering Plants, with Particular Emphasis on Asterids[END_REF][START_REF] Pfeil | Paralogy and orthology in the Malvaceae rpb2 gene family: Investigation of gene duplication in Hibiscus[END_REF]. However, resolving duplications becomes more difficult when copy number changes quickly [START_REF] Ashfield | Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine[END_REF], or when duplications are recent and copy loss is complete or nearly so, thus returning the locus to a single-copy state [START_REF] Ramadugu | Coalescence simulation testing of hybridization versus lineage sorting in Citrus (Rutaceae) using six nuclear genes[END_REF]. In the latter case, the phylogenetic pattern can mimic that of ILS and become indistinguishable from it [START_REF] Sousa | Using genomic location and coalescent simulation to investigate gene tree discordance in Medicago L[END_REF], generally leaving no trace at all of the loss.

Why is GDL so challenging to implement in theory? The topological and coalescent-time similarities between ILS and GDL complicates extending the MSC to include both processes, unless copy number exceeds one in at least some samples (Fig. 4). Assuming that allelic and homoeologous variation is not confused with the copy number of independently duplicated genes, at the very least, duplicated genes could be handled as independent loci with missing data for some samples with MSC inference.

When copy loss is complete, or when the duplication is so recent so as to conflate allelic versus copy variation, these GDL loci have little effect on species phylogeny inference and divergence times, especially if the algorithms used employ averages over coalescence times or other parameters across many gene trees [START_REF] Liu | Estimating species phylogenies using coalescence times among sequences[END_REF][START_REF] Sousa | Using genomic location and coalescent simulation to investigate gene tree discordance in Medicago L[END_REF]. At high proportions, though, they may cause serious issues for phylogenetic reconstruction, because the unexpected positions of gene duplications in a species phylogeny, coupled with random copy loss, means that no specific pattern is expected among the affected gene trees (Fig. 4). This scenario contrasts with the retention of ancestral polymorphisms, where we know that branches in short species phylogenies (in coalescent units) are the cause [START_REF] Rosenberg | Genealogical trees, coalescent theory and the analysis of genetic polymorphisms[END_REF]. Thus, we expect deeply coalescing lineages to occur in specific parts of a species phylogeny with a limited number of topological outcomes and branch lengths limited by effective population size, which is not the case for duplicated genes. A recent approach to identifying genes that are single copy, but have nonetheless been affected by GDL, was made using the genomic location of the loci [START_REF] Sousa | Using genomic location and coalescent simulation to investigate gene tree discordance in Medicago L[END_REF], and could prove useful for distinguishing GDL and ILS.

Recombination

All existing methods for coalescent estimation of species trees and networks make two important assumptions, namely that (1) there is free recombination between loci, and (2) there is no recombination within a locus. These two assumptions address a key concept distinguishing MSC models from concatenation or supermatrix models: it is the conditional independence of loci, mediated by recombination between loci, and not the ability to address ILS or discordance among genes per se.

Moving forward, three important questions to address are: (1) How robust are methods to the presence of recombination within loci and/or to the violation of independence among loci? (2) How should we model recombination within the species phylogeny inference framework? and (3) How do we detect it and differentiate recombination-free loci?

Researchers have started to examine the first question and found a detectable effect of recombination only under extreme levels of ILS and gene tree heterogeneity (e.g., [START_REF] Lanier | Is Recombination a Problem for Species-Tree Analyses?[END_REF]. However, more analyses and studies are still needed to explore a wider range of factors and parameters that could affect phylogenetic inference when the assumption of recombination-free loci is violated (e.g., [START_REF] Li | Recombination-aware phylogenomics unravels the complex divergence of hybridizing species[END_REF]. For answering the second question, one approach involves combining the multispecies coalescent with hidden Markov models (e.g., [START_REF] Hobolth | Genomic relationships and speciation times of human, chimpanzee, and gorilla inferred from a coalescent hidden Markov model[END_REF]. These methods suffer from the "state explosion problem", where individual states are needed for the different coalescent histories, and they increase rapidly with the number of taxa in the dataset, making them infeasible except for very small (~4 taxa) datasets. New methods that scale to larger datasets are needed if such approaches are to be useful in practice. A different direction is to devise novel methods for inferring species phylogeny while assuming that the genealogies of the individual loci could take the form of an ancestral recombination graph (ARG: [START_REF] Griffiths | An ancestral recombination graph[END_REF][START_REF] Siepel | Phylogenomics of primates and their ancestral populations[END_REF][START_REF] Rasmussen | Genome-Wide Inference of Ancestral Recombination Graphs[END_REF].

Extending these approaches to address recombination would require the development of new models that significantly extend the multispecies coalescent to account for ARGs within the branches of a species phylogeny. For two-step species tree methods, this entails developing new methods that infer ARGs for the individual loci and methods that infer species phylogenies from collections of ARGs. For single-step Bayesian methods, novel developments are needed to sample species phylogenies, locusspecific ARGs, and their related parameters. It will also be important to better understand the conditions under which ignoring recombination will still yield reasonable estimates of phylogeny. Extending the theory to accommodate ARGs may be of intrinsic interest, but if the parameter space in which recombination is relevant is very small, then practitioners may be able to ignore recombination.

Species concepts and delimitation

Coalescent methods have played an important role in the development and critical evaluation of species delimitation methods because they provide hypotheses for species boundaries based on genetic data and be integrated with phenotypic data (e.g., [START_REF] Solís-Lemus | Bayesian species delimitation combining multiple genes and traits in a unified framework[END_REF]. Irrespective of traditional species concepts, it is essential that the entities at the tips of the species tree do not violate the assumptions of the MSC, wherein species are defined mathematically (e.g., [START_REF] Rannala | Bayes estimations of species divergence times and ancestral population sizes using DNA sequences from multiple loci[END_REF][START_REF] Degnan | Gene tree discordance, phylogenetic inference and the multispecies coalescent[END_REF]: the branches of the species tree constitute species or populations that do not exchange genes.

However, the MSC model also carries strict assumptions about the divergence process if the delimited units are to be interpreted as species. Specifically, it is important to emphasize that in the "standard" MSC model, these species represent populations that, immediately after divergence, no longer experience gene flow. Therefore, the species of the MSC model do not necessarily correspond with species as a taxonomic rank, defined by traditional species concepts [START_REF] Heled | Bayesian inference of species trees from multilocus data[END_REF]: "MSC" species could simply be populations by other criteria, so long as they have ceased to exchange genes, even for a short period of time. In other words, a species tree built under the MSC might then be interpreted as a depiction of the history of the barriers to gene flow among diverging structured populations and this may be particularly true when there is dense spatial and genomic sampling across individuals (Sukumaran & Knowles, 2017).

Therefore, in those species tree methods requiring a priori assignments of individuals to species, such assignments may strongly influence the inferred species phylogeny, in the same way that hybridization will have serious consequences on an estimated species phylogeny [START_REF] Leaché | The influence of gene flow on species tree estimation: A simulation study[END_REF].

Recently, several MSC-based methods that have the ability to simultaneously perform species delimitation and estimate the species phylogenies have been developed and implemented (e.g., [START_REF] Yang | Unguided species delimitation using DNA sequence data from multiple loci[END_REF][START_REF] Jones | DISSECT: an assignment-free Bayesian discovery method for species delimitation under the multispecies coalescent[END_REF][START_REF] Jones | Algorithmic improvements to species delimitation and phylogeny estimation under the multispecies coalescent[END_REF]. These methods seem to consistently recover the correct number of "MSC species" given the assumptions of the model. However, it is probable that the assumption of no gene flow between the descendant populations is often violated and that most reproductive isolation processes are gradual or episodic rather than sudden and permanent (e.g., [START_REF] Rosindell | Protracted speciation revitalizes the neutral theory of biodiversity[END_REF]. There is thus need for methods that perform simultaneous species phylogeny estimation and assignment of individuals to species while considering the limitations of the MSC [START_REF] Jones | DISSECT: an assignment-free Bayesian discovery method for species delimitation under the multispecies coalescent[END_REF].

If one prefers a species concept that affirms that most recently diverged populations are necessarily reproductively isolated, current methods will overestimate the number of species as defined by traditional species concepts and will likely reveal instead intraspecific population structure (Sukumaran & Knowles, 2017). [START_REF] Toprak | Species Delimitation Without Prior knowledge: DISSECT Reveals Extensive Cryptic Speciation in the Silene aegyptiaca complex (Caryophyllaceae)[END_REF] used DISSECT [START_REF] Jones | DISSECT: an assignment-free Bayesian discovery method for species delimitation under the multispecies coalescent[END_REF] but also employed checks as to the integrity of various hypotheses of species boundaries suggested by the data. From a computational point of view, any species delimitation method will need an operational definition of species. Therefore, a possible development of MSC-based species delimitation methods could be allowing migration and assuming that speciation is complete when a certain proportion of the migrations is reached or when the migration rate is sufficiently low. However, this solution will not be suited for the protracted speciation model because other kinds of information besides the movement of genes will still be needed to identify when a clade becomes reproductively isolated. Possibly the best way to avoid confusion is to restrict the word "species" to taxonomy and base it on multiple sources of information which are synthesized in an integrative fashion [START_REF] Dayrat | Towards integrative taxonomy[END_REF][START_REF] Will | The perils of DNA Barcoding and the need for integrative taxonomy[END_REF][START_REF] Bacon | Evaluating multiple criteria for species delimitation: an empirical example using Hawaiian palms (Arecaceae: Pritchardia)[END_REF][START_REF] Solís-Lemus | Bayesian species delimitation combining multiple genes and traits in a unified framework[END_REF], and refer to the reproductively isolated units of MSC analysis as "MSC units" or "MSC taxa".

Models at the intersection of phylogenomics and macroevolution

At the intersection of phylogenomics and macroevolution, phylogeneticists aim at shedding light on how patterns of organismal diversity have been generated and maintained through time. The former focuses primarily on building phylogenies, whereas the latter uses them to study the tempo and mode of diversification over time. In many important respects, these two sub-disciplines have remained distinct and non-communicative. On the one hand, phylogenomics and phylogeography have not exhaustively aimed to address the type of questions -related to diversification and trait evolution -that macroevolution focuses on. On the other hand, macroevolution ignores many kinds of complexities inherent to the phylogeny building process that phylogenomics has recently begun to address.

Macroevolutionary models focus on long-term processes, in terms of both species richness and phenotypic diversity. They rely on two types of models: birth-death models of diversification aimed at understanding how and why speciation and extinction rates vary through time and across lineages [START_REF] Hey | Using phylogenetic trees to study speciation and extinction[END_REF][START_REF] Nee | Tempo and mode of evolution revealed from molecular phylogenies[END_REF]; see [START_REF] Stadler | Recovering speciation and extinction dynamics based on phylogenies[END_REF][START_REF] Morlon | Phylogenetic approaches for studying diversification[END_REF] for review) and models of trait evolution aimed at understanding the mode and tempo of phenotypic evolution [START_REF] Felsenstein | Maximum-likelihood estimation of evolutionary trees from continuous characters[END_REF]see Pennell &Harmon, 2013 andManceau, Lambert &Morlon, 2017 for reviews). These models are typically constructed at the level of species, ignoring the populations or individuals that constitute these species (but see Manceau, Lambert &[START_REF] Manceau | Phylogenies support out-of-equilibrium models of biodiversity[END_REF]Rosindell, Harmon &Etienne, 2015 for exceptions). As a consequence, microevolutionary processes such as coalescence have informed phylogenetic methods for building phylogenies more so than have macroevolutionary methods that use them. For example, the most widely used phylogenetic dating methods generally do not acknowledge the critical distinction between speciation times, which are usually of primary interest, and coalescence times, which are often assumed to represent speciation times but in fact represent events older than the divergence of the species concerned [START_REF] Edwards | Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies[END_REF][START_REF] Dos Reis | Bayesian molecular clock dating of species divergences in the genomics era[END_REF][START_REF] Angelis | The impact of ancestral population size and incomplete lineage sorting on Bayesian estimation of species divergence times[END_REF]. In addition, macroevolutionary models are fit to species phylogenies (diversification models) or a combination of species phylogenies and phenotypic data (trait evolution models), most often assuming that evolution is best represented by a species tree, not a network (but see [START_REF] Jhwueng | Trait evolution on phylogenetic networks[END_REF][START_REF] Bastide | Phylogenetic Comparative Methods on Phylogenetic Networks with Reticulations[END_REF][START_REF] Solís-Lemus | PhyloNetworks: A Package for Phylogenetic Networks[END_REF] for models of trait evolution on networks), and that the species phylogeny is known. Nearly all models that use phylogenies to study character evolution assume a single underlying species phylogeny on which characters evolve. But it has become evident recently that different characters in principle have different phylogenies, for the same reason that genes themselves might have different phylogenies (Hahn & Nahkleh, 2016). Analyzing incongruences between character evolution inferred from the species tree versus from gene trees that are more directly linked to the character under study would provide a refined understanding of character evolution. Recent work on the phylogeny of quantitative characters may be helpful in this endeavor [START_REF] Felsenstein | A comparative method for both discrete and continuous characters using the threshold model[END_REF].

Developing research projects that integrate the heterogeneity inherent in phylogenomics and macroevolution will bring important new insights into the evolutionary process. For example, developing diversification and phenotypic evolution models to be fit to networks rather than dichotomous trees will allow estimates of rates of hybrid speciation and phenotypic evolution as well as a better understanding of factors influencing such rates [START_REF] Bastide | Phylogenetic Comparative Methods on Phylogenetic Networks with Reticulations[END_REF]. Embracing genetic heterogeneity and the incongruence between gene trees and species phylogenies when applying macroevolutionary models could help us to better understand how speciation proceeds, and also to analyze the coupling between genetic and phenotypic evolution (e.g., is phenotypic convergence coupled or not with genetic convergence in relevant genes?). Developing macroevolutionary models accounting for within-species heterogeneity linked to biogeography could help us understand how biogeographic structuring influences speciation, extinction, and phenotypic evolution.

More generally, evolutionary biologists have yet to thoroughly explore the type of new questions that we are going to be able to address if we are given genomic data at the tips of all species from a phylogeny. Such data could allow us to gain an integrative understanding of three fundamental aspects of evolution: evolution at the molecular level, at the phenotypic level, and at the clade level, as well as the links among them. Are rates of evolution at these three levels correlated? If so, how? Do features of genomes or of genome evolution, such as quantity of transposable elements, substitution rates, number of gene duplications, influence rates of diversification and phenotypic evolution? Clearly, we are only at the beginning of exploring these new possibilities.

Mapping trait evolution on heterogeneous genomic datasets

Mapping the genomic basis of phenotypic traits is a major trend in evolutionary biology today [START_REF] Elmer | Adaptation in the age of ecological genomics: insights from parallelism and convergence[END_REF][START_REF] Hoban | Finding the Genomic Basis of Local Adaptation: Pitfalls, Practical Solutions, and Future Directions[END_REF]. Such mapping can be conducted in the context of populations of a single species or, increasingly, via comparisons of species on a phylogeny (e.g., [START_REF] Hiller | A "forward genomics" approach links genotype to phenotype using independent phenotypic losses among related species[END_REF][START_REF] Marcovitz | Reverse Genomics" Predicts Function of Human Conserved Noncoding Elements[END_REF]. Phylogenetic genome-wide association studies ("PhyloGWAS") methods identify genomic features in coding or non-coding DNA that exhibit unusual patterns of evolution on branches concerned with repeated evolution of phenotypes, thereby drawing connections between the genomic and phenotypic levels [START_REF] Pease | Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation[END_REF]. Such phylogenomic mapping usually assumes a single phylogeny, the species phylogeny, as a framework for analysis, and therefore ignores genomic heterogeneity. To make phyloGWAS mapping most efficient it might be more appropriate to use the local topology in the genome for inference and estimation of ancestral states. Estimating genotypephenotype associations solely on the species phylogeny might yield misleading results regarding the origin and evolution of phenotypic traits (Hahn & Nakleh, 2016;[START_REF] Mendes | Gene tree discordance can generate patterns of diminishing convergence over time[END_REF][START_REF] Guerrero | Quantifying the risk of hemiplasy in phylogenetic inference[END_REF]. Heterogeneity across gene histories has been traditionally considered as "biological noise" when using comparative genomics to map traits, but of course such heterogeneity is the focus of gene mapping efforts at lower taxonomic levels. A recently proposed application of the MSC for quantitative traits, accounting for genealogical heterogeneity improves downstream estimates of mean trait values, phylogenetic signal, and evolutionary rates of traits [START_REF] Mendes | A multispecies coalescent model for quantitative traits[END_REF]. Furthermore, genome-wide or gene-specific selective sweeps associated with the evolution of a particular phenotypic trait are a major source of genetic heterogeneity among closely related populations or species, and can be captured using outlier statistics, such as Fst or Dxy [START_REF] Pease | Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation[END_REF]. Such selective sweep mapping of genes with large phenotypic effect can now be accomplished with high resolution and precision in genomically poorly studied organisms [START_REF] Lamichhaney | Evolution of Darwin's finches and their beaks revealed by genome sequencing[END_REF]. Apart from providing valuable knowledge on the genetic basis of trait diversification, such data are providing increasing support to the fact that cases of genetic heterogeneity can be profitably used in the effort to understand and resolve evolutionary history, rather than considering it "biological noise." Such thinking needs to be incorporated into comparative genomics more frequently (e.g., [START_REF] Mendes | A multispecies coalescent model for quantitative traits[END_REF].

Tree-free methods of character evolution

We have seen that incorporating phylogenetic heterogeneity is a challenge for macroevolutionary models of character evolution. At the other end of the spectrum are a class of methods (so called "tree-free methods") that attempt to draw inferences and principles about trait evolution without assuming a particular phylogeny. The common situation when analyzing character or trait data correlated by a phylogeny is to assume a stochastic process for the trait, commonly a variation of the Brownian motion (BM; [START_REF] Felsenstein | Phylogenies and the comparative method[END_REF] or Ornstein-Uhlenbeck (OU; [START_REF] Hansen | Stabilizing selection and the comparative analysis of adaptation[END_REF] processes. Then, using the estimated phylogeny and measured trait data for each species, the parameters of various evolutionary processestrait variation, patterns and rates of change, etc. -are estimated, often using maximum-likelihood or Bayesian approaches (see Pennell &Harmon, 2013 andManceau, Lambert &[START_REF] Manceau | A unifying comparative phylogenetic framework including 63 traits coevolving across interacting lineages[END_REF] for reviews).

However, given the various logistical and technical challenges of inferring robust phylogenies, exploring tree-free methods might represent a useful mechanism for guiding the study of character evolution for certain groups.

Tree-free comparative methods work by integrating over the space of trees (under a given branching process model). For example, under a pure birth model and with enough tip measurements, the optimum value of the OU process can be estimated as the sample average (Bartoszek & Sagitov, 2015a). Similar results have now been derived for other models of tree growth that include extinction (Adamczak and Miłoś, 2014;2015;[START_REF] Ané | Phase transition on the convergence rate of parameter estimation under an Ornstein-Uhlenbeck diffusion on a tree[END_REF]. Similarly, the rate of adaptation under the OU process, often modeled as the stationary variance -the ratio of the squared "rate of evolution" (sigma parameter in the OU model) and twice the "rate of adaptation" (the alpha parameter) can be estimated as the sample variance (Bartoszek & Sagitov, 2015a). Teasing sigma and alpha apart, however, requires a tree. The key parameter of the BM model, the rate of evolution, is similarly estimable directly from the trait sample (Bartoszek & Sagitov, 2015b;[START_REF] Crawford | Diversity, disparity, and evolutionary rate estimation for unresolved Yule trees[END_REF], whereas the root state cannot be consistently estimated without a tree [START_REF] Ané | Analysis of comparative data with hierarchical autocorrelation[END_REF][START_REF] Sagitov | Interspecies correlation for neutrally evolving traits[END_REF]. In addition to providing tree-free estimators of some model parameters, the studies mentioned above also derived Central Limit Theorems that allow computing confidence intervals around these point estimates as well as the sample sizes needed to obtain reliable estimates.

Extinct and unsampled species

A notable case when phylogenomics and macroevolution do meet is in the treatment of extinct or unsampled species in phylogenetic reconstruction and dating. Despite the avalanche of genomic data for an increasing number of species, we still lack sequence data for most species, making it difficult to place them in a phylogeny. Some researchers (e.g., [START_REF] Jetz | The global diversity of birds in space and time[END_REF][START_REF] Tonini | Fully-sampled phylogenies of squamates reveal evolutionary patterns in threat status[END_REF] have used a combination of consensus trees and current taxonomic classifications to impute the phylogenetic relationships of unsampled species. In this case, polytomies are often resolved by using distributions of branching times obtained from macroevolutionary birth-death models [START_REF] Kuhn | A Simple Polytomy Resolver for Dated Phylogenies[END_REF].

While we appreciate the value of these approaches given the real logistic difficulties researchers face as they attempt to obtain samples from around the globe and that methods are now easily applicable, such comparisons (reviewed by [START_REF] Ksepka | Synthesizing and databasing fossil calibrations: divergence dating and beyond[END_REF] or estimated fossil ages have been used to post-hoc-scale branch lengths [START_REF] Sanderson | Estimating Absolute Rates of Molecular Evolution and Divergence Times: A Penalized Likelihood Approach[END_REF][START_REF] Smith | treePL: divergence time estimation using penalized likelihood for large phylogenies[END_REF][START_REF] Bapst | Assessing the effect of time-scaling methods on phylogeny-based analyses in the fossil record[END_REF]. However, there is an increasing interest in simultaneously performing phylogenetic inference and time-scaling phylogenies [START_REF] Pyron | Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia[END_REF][START_REF] Ronquist | A total-evidence approach to dating with fossils, applied to the early radiation of the Hymenoptera[END_REF][START_REF] Heath | The fossilized birth-death process for coherent calibration of divergence-time estimates[END_REF]. These methods, known as tip-dating, have the advantage of accommodating phylogenetic uncertainty implemented in a Bayesian model-based framework [START_REF] Hunt | Integrating Paleontological and Phylogenetic Approaches to Macroevolution[END_REF] and have been implemented in several groups spanning a wide temporal scale (e.g., [START_REF] Pyron | Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia[END_REF][START_REF] Pyron | Trait-dependent diversification and the impact of paleontological data on evolutionary hypothesis testing in New World ratsnakes (tribe Lampropeltini)[END_REF][START_REF] Ronquist | A total-evidence approach to dating with fossils, applied to the early radiation of the Hymenoptera[END_REF][START_REF] Heath | The fossilized birth-death process for coherent calibration of divergence-time estimates[END_REF][START_REF] Slater | Iterative adaptive radiations of fossil canids show no evidence for diversity-dependent trait evolution[END_REF][START_REF] Larson-Johnson | Phylogenetic investigation of the complex evolutionary history of dispersal mode and diversification rates across living and fossil Fagales[END_REF][START_REF] Matzke | Inferring node dates from tip dates in fossil Canidae: the importance of tree priors[END_REF]. Lastly, these phylogenies containing both extinct and extant diversity serve as framework to conduct robust diversification and trait evolution studies [START_REF] Hunt | Testing the link between phenotypic evolution and speciation: an integrated palaeontological and phylogenetic analysis[END_REF][START_REF] Pennell | An integrative view of phylogenetic comparative methods: Connections to population genetics, community ecology, and paleobiology[END_REF][START_REF] Pyron | Post-molecular systematics and the future of phylogenetics[END_REF].

Building, updating and sustaining the Tree of Life

Scalability challenges

Inferring the phylogeny of all living organisms represents a different challenge than inferring the relationships of just a few terminals; often the scale at which new methods are developed and tested is on this latter scale. For instance, for eukaryotes alone, recent conservative estimates indicate that there are ~8.7 million species on Earth and only 9-14% of them have been formally described [START_REF] Mora | How Many Species Are There on Earth and in the Ocean?[END_REF]. Furthermore, out of 2.6 million taxa currently represented in the Open Tree of Life (https://tree.opentreeoflife.org; [START_REF] Hinchliff | Synthesis of phylogeny and taxonomy into a comprehensive Tree of Life[END_REF], only ~55,000 were gathered from hard-data phylogenies, whereas phylogenetic affinities of the rest were inferred from current taxonomic classifications [START_REF] Mctavish | Phylesystem: A git-based data store for community-curated phylogenetic estimates[END_REF]2017). These observations suggest that the vast majority of taxa on Earth still await formal taxonomic description and placement in the ToL [START_REF] Mora | How Many Species Are There on Earth and in the Ocean?[END_REF][START_REF] Mctavish | How and why to build a unified Tree of Life[END_REF]. As mentioned above, one common challenge that phylogeneticists encounter is the difficulty in accessing samples from rare, endangered, or extinct taxa, particularly in countries where collecting and exporting is difficult. Recent genomic techniques now allow successful results in obtaining valuable ancient DNA data from museum specimens (e.g., [START_REF] Staats | Genomic treasure troves: complete genome sequencing of herbarium and insect museum specimens[END_REF][START_REF] Hykin | Fixing Formalin: A Method to Recover Genomic-Scale DNA Sequence Data from Formalin-Fixed Museum Specimens Using High-Throughput Sequencing[END_REF][START_REF] Mccormack | Sequence capture of ultraconserved elements from bird museum specimens[END_REF][START_REF] Mccormack | Transforming Museum Specimens into Genomic Resources[END_REF][START_REF] Ruane | Phylogenomics using formalin-fixed and 100+ year-old intractable natural history specimens[END_REF], and here, we advocate for routine use of these resources to improve taxon sampling, enhance research in phylogenomics and phylogeography, and increase awareness and usefulness of natural history collections.

Despite the great increase in the generation of genomic data across organisms, we are often limited to fast and efficient but simpler, less realistic phylogenetic methods and assumptions, such as IQ-Tree [START_REF] Nguyen | IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies[END_REF], PhyloBayes [START_REF] Lartillot | PhyloBayes MPI: Phylogenetic Reconstruction with Infinite Mixtures of Profiles in a Parallel Environment[END_REF], and ExaML [START_REF] Kozlov | ExaML version 3: a tool for phylogenomic analyses on supercomputers[END_REF], to deal with large, heterogeneous datasets. The main downside of these options is that they are not a full coalescent framework and thus far rely on data concatenation. Fully coalescent methods, such as the popular phylogenetic software program *BEAST [START_REF] Heled | Bayesian inference of species trees from multilocus data[END_REF], are not capable of dealing with more than a few hundred taxa and some dozen loci at a time for a common analysis, and only recently the release of StarBeast2 allows for the use of thousands of loci for tens of taxa [START_REF] Ogilvie | StarBEAST2 Brings Faster Species Tree Inference and Accurate Estimates of Substitution Rates[END_REF]. To tackle this problem, we encourage the continuing development of methods that are fully scalable and ideally only increase analytical time linearly rather than exponentially with the number of taxa and loci. Phylogenetic methods should also be fully parallelizable (in order to run natively in computer clusters) and contain checkpoints, i.e., be able to resume the analyses from the latest logged file in case an analysis crashes or the user wishes to evaluate partial results. Another point of possible improvement is in dealing with new sequences to be added to a previously large dataset: should the analysis start from scratch, or could there be substantial time gains by letting those sequences find their placement in the phylogeny 'on the fly' (e.g., [START_REF] Siu-Ting | Evolutionary relationships of the Critically Endangered frog Ericabatrachus baleensis Largen, 1991 with notes on incorporating previously unsampled taxa into large-scale phylogenetic analyses[END_REF].

Large scale phylogenies should ideally be based on the best (or most comprehensive) available datasets in terms of taxonomic and molecular sampling and be constructed from the data itself. However, even supermatrix inference conducted under a single analysis can add bias on tree heights and coalescence times when performed across unbalanced sampled clades (a very common case for speciesrich clades or understudied taxa), and therefore affect downstream analyses that rely on these parameters (e.g., biogeography, trait evolution, diversification rates). Computing optimally populated datasets that combine the largest number of taxa and loci simultaneously is a complex mathematical problem, but recent approaches (e.g., SUPERSMART -Antonelli et al., 2017;PyPHLAWD -Smith & Brown, 2018) attempt to overcome it objectively, such as applying the knapsack problem to phylogenetics by packing the optimal choice of species and suitable alignments into a minimally sparse supermatrix.

Community initiatives

Building the ToL is a grand challenge in molecular phylogenetics, and one that cannot be accomplished by a single person or institution's efforts. Several initiatives have been developed in recent years to coordinate efforts and provide the research community with synthetic information. A prominent project is the Open Tree of Life (https://tree.opentreeoflife.org/; [START_REF] Hinchliff | Synthesis of phylogeny and taxonomy into a comprehensive Tree of Life[END_REF]. This project provides a synthesis of previously published phylogenies merged through supertree and other grafting methods. One issue faced by the initiative is that it relies on authors uploading their phylogenetic trees to open data repositories, such as Dryad Data Repository (http://datadryad.org/pages/organization; [START_REF] Vision | Open data and the social contract of scientific publishing[END_REF] or TreeBase [START_REF] Sanderson | TreeBASE: a prototype database of phylogenetic analyses and an interactive tool for browsing the phylogeny of life[END_REF][START_REF] Piel | Treebase v. 2: A Database of Phylogenetic Knowledge[END_REF], which at least until recently only occurred in about 17% of cases [START_REF] Drew | Data deposition: Missing data mean holes in Tree of Life[END_REF]. Substantial curatorial efforts are also critical to facilitate reusability of deposited trees [START_REF] Mctavish | Phylesystem: A git-based data store for community-curated phylogenetic estimates[END_REF]. A different approach was taken by [START_REF] Antonelli | Towards a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa[END_REF], who developed a framework for continuously inferring time-calibrated large phylogenies from raw sequence data deposited in GenBank (Clark et al., 2016) in a multi-step method. Similarly, various tools have been developed to make information contained in the ToL available for the general public (e.g., [START_REF] Rosindell | OneZoom: A Fractal Explorer for the Tree of Life[END_REF][START_REF] Harmon | Arbor: Comparative Analysis Workflows for the Tree of Life[END_REF].

Mapping the Tree of Life

While progress has been made in mapping species distributions at the large scale aiming for improved conservation practices (e.g., the Map of Life collaborative project; https://mol.org/), most initiatives do not map the tips of phylogenetic trees directly onto the geographic space, and therefore are limited by current taxonomic knowledge. As spatial variation in biodiversity results from interactions between evolutionary history and environmental factors, explicit connections between the tips of the ToL and geographic ranges will greatly improve biogeographic inferences [START_REF] Quintero | Historical Biogeography Using Species Geographical Ranges[END_REF] and our understanding of biodiversity patterns and future trends. Advances in mapping the ToL through earth history using genomic-based phylogenetic inferences over broad scales and explicit spatial models (e.g., geophylogenies and continuous diffusion models: [START_REF] Kidd | Geophylogenies and the Map of Life[END_REF]) depend directly on locality data that should be made available in raw and ready-to-use formats. Data sharing policies for associated data, such as geographic coordinates and voucher information, is not well established among journals. We argue that editorial boards should try as best as possible to establish data policies that value and encourage the deposit of geographic data associated to vouchered specimens and other associated information available for future reference.

Best practices for building the Tree of Life

Data must be well curated in databases and publicly available

As we are now in the era of big data in biological sciences, adequate reproducibility must be a fundamental endeavor of biodiversity research. Therefore, data publication in open access repositories represents a powerful tool that not only ensures long-term storage and public availability for future research, but also serves as a vehicle for clarifying intellectual rights and scientific merits [START_REF] Costello | Best practice for biodiversity data management and publication[END_REF]. Biocuration, the activity of organizing, representing, and making biological information accessible to both biologists and bioinformaticians, has now become an important consideration in building, updating, and sustaining the ToL [START_REF] Mctavish | How and why to build a unified Tree of Life[END_REF]. The exponential growth in the amount of genomic scale data and the increased dependence on the availability of each other's' data to answer complex biological questions means that there is a need for improved data management, analysis, and accessibility. GenBank has been the main open access repository for annotated collections of publicly available molecular data. Although the data stored in this database usually lists information such as organism of origin and publication details, the utility of molecular data in this database to answer multiple biological questions, such as biogeographic patterns of biodiversity, is often hampered by lack of associated information such as collection locality or [START_REF] Scotch | Enhancing phylogeography by improving geographical information from GenBank[END_REF]Gratton et al., 2017) attachment to a specific voucher specimen. Moreover, recent surveys have shown that fewer than 20% of phylogenetic studies provide access to phylogenetic data (i.e., alignments and phylogenies) and when they do, critical biological information such as complete taxon names is missing [START_REF] Drew | Lost Branches on the Tree of Life[END_REF][START_REF] Magee | The Dawn of Open Access to Phylogenetic Data[END_REF].

We propose two urgent actions to advance this key field. First, authors should be encouraged to submit molecular data that is linked to voucher specimens deposited in recognized scientific collection and museums. Second, authors, journals, and curators should encourage all molecular data submitted to include information such as collection locality and details of voucher specimens. In this regard, other global initiatives such as the International Barcode of Life Project (iBOL; http://www.ibolproject.org) have had great success linking molecular data with morphological and distributional data. When all the data produced or published are curated to high standards and made accessible as soon as available, biological research will be able to process massive amounts of complex data much more quickly.

Submitting sequence and tree data during publication is now routine. However, making available all analytical methods such as software and code used to process and analyze data is less widely employed by the phylogenetic community. Facilities such as TreeBase, Dryad Digital Repository, and Github (https://github.com/) provide a platform for the curated storage of the data and bioinformatic pipelines underlying the scientific literature (see [START_REF] Mctavish | Phylesystem: A git-based data store for community-curated phylogenetic estimates[END_REF]2017). Authors and journals should require all published research to include links to raw data, processed data, and all analytical methods used to produce the results presented. In general, we advocate for following best practices of data management and publication to ensure the quality and utility of phylogenomic data and their associated biological information [START_REF] Stoltzfus | Sharing and re-use of phylogenetic trees (and associated data) to facilitate synthesis[END_REF][START_REF] Drew | Lost Branches on the Tree of Life[END_REF][START_REF] Costello | Best practice for biodiversity data management and publication[END_REF].

In putting together Figure 2, for example, we found that basic information on a given phylogenomic study, such as the number of species or sequences analyzed, or the total number of base pairs in an alignment, were often not reported or difficult to recover; including such information in easy-to-access tables prior to article acceptance would greatly facilitate meta-analyses and syntheses as the number of studies grows (Supplementary Table S1).

The need of adequate curation of analytical tools

In the same way that data must be adequately stored and curated, analytical tools must be available for future use and should guarantee proper reproducibility [START_REF] Wilson | Best Practices for Scientific Computing[END_REF][START_REF] Darriba | The state of software for evolutionary biology[END_REF]. One of the reasons behind the dramatic increase in the number of phylogeographic and phylogenetic studies during the last 20 years is the proliferation of software and bioinformatic tools to process and analyze these data. Thanks to these new methods, it is now possible to implement a wide array of theoretical models that sustain the fields of phylogenomics and phylogeography. As stated above, genome-wide data have notoriously increased the necessity to expand our analytical models, ultimately leading to a stronger demand for computing resources [START_REF] Darriba | The state of software for evolutionary biology[END_REF]. Given their key role in phylogenomic research, it is advisable that software development, documentation, and availability follow the best possible practices (e.g., [START_REF] Leprevost | On best practices in the development of bioinformatics software[END_REF][START_REF] Wilson | Best Practices for Scientific Computing[END_REF][START_REF] Guang | An Integrated Perspective on Phylogenetic Workflows[END_REF][START_REF] Darriba | The state of software for evolutionary biology[END_REF]. Having both data and analytical tools adequately stored and accessible to the public not only will ensure high reproducibility of previous studies, but, more importantly, will facilitate continuing the construction of the ToL [START_REF] Mctavish | How and why to build a unified Tree of Life[END_REF].

All contributions toward building the Tree of Life must be properly recognized

Some current publishing practices in the scientific community may unintentionally represent hurdles toward the ultimate end of collecting and disseminating phylogenetic data on which to build a ToL. For instance, the increasing need in many countries and communities for publishing high-impact papers understandably often discourages researchers from releasing their data until their studies are complete and have passed the peer-review process. This is partially explained by the heavy emphasis of top journals on unusually novel and flashy findings as compared to those studies that represent more modest, but just as critical, advances in the understanding of the phylogenetic relationships of the groups. Similarly, this urge to publish high-impact papers often impedes adequate long-term studies that could potentially generate a wider variety of basic data. With the cultural emphasis on impact and numbers and rates of publication, in practice there is often a penalty for long-term studies. Our current climate often values novel results produced in the short term. Consequently, as a community, we must reach an equilibrium between shortand long-term scientific production in a way that values both, encouraging high impact studies bringing radical reorganizations of the ToL, without hurting lower impact research and the ongoing search for innovation.

Moreover, because building the ToL is a slow and daunting task, it is important that, as a scientific community, all contributors to the process receive proper recognition for their contributions, thereby keeping motivation high and retaining our best talent. Unfortunately, some contributors, both institutions and roles within them, receive less recognition in this grand task than others. For example, field biologists that obtain basic natural history information and specimens used for building the ToL (Suarez & Tsutsui, 2004), and the natural history museums that house those specimens, are often not recognized sufficiently. As a community, we have been following a trend in which, perhaps inadvertently, we do not value as much the production of basic biological and natural history data. This can certainly be recognized in our national funding practices, which often do not support basic taxonomic or natural history fieldwork at the expense of flashier end-uses of biological specimens. Specimens are the foundation of most phylogenomic and phylogeographic studies, and we should find standard mechanisms not only to acknowledge, but also to encourage the production of these data in an integrative framework.

It is time to strengthen those initiatives aimed at recognizing scientific production beyond citations of peer-reviewed literature (e.g., ORCID; https://orcid.org) by giving also credit to the production and impact of basic biology datasets and collected specimens. Providing credit for depositing and generating data by tracking, for example, number of access and downloads or number of studies using genetic data associated to specimens, could represent a formal recognition of the importance of producing and sharing basic biological data could help bridge the gap between naturalists, taxonomists, empiricists, and mathematicians invested on the study of life history.

It will be exciting to have objective estimates that allow tracking the direct and indirect impact of how these data and samples are being used. We are confident that such initiatives will highlight the importance of continuing field-and museum-based research in various fields of biological research [START_REF] Buerki | Collections-based research in the genomics era[END_REF]. Furthermore, such cultural shifts will undoubtedly encourage discerning young minds to embrace basic biological research in their academic endeavors, rather than embracing more lucrative and societally appreciated applied fields.

Conclusions

In this perspective, we have attempted to cover ground in the vast arena of issues facing modern phylogenomics today. In Table 2, we summarize some of the most pressing challenges that the field of phylogenomics is experiencing as we use coalescent-based methods toward building the ToL. We have seen how genome-scale phylogenomics, currently on a strong footing as a result of the multispecies coalescent model, is increasingly improved by models that recognize reticulate processes, such as recombination and introgression. In contrast, macroevolutionary models that use phylogenies have yet to embrace the heterogeneity that currently drives many theoretical innovations in phylogenetic reconstruction itself. We have emphasized the need for the phylogenomics community to embrace high standards of data quality, curation and accessibility in its long-term pursuit of the ToL. Such a grand mission requires value and recognition placed not only on the end products of the process, such as publications and trees, but also on the natural history specimens on which phylogenies are based and which are cared for by the community of natural history museums. Building the ToL will require contributions from all sectors of biological and related sciences -from field biology to theory and everything in between -and robust cyberinfrastructures to integrate these diverse and increasingly massive data streams.

Figures

Figure 1. A posteriori marker selection from whole genome alignments for phylogenomics and phylogeography. Whole genome analysis (A) permits researchers to choose different markers for specific purposes (B -D). By contrast, subsampling methods such as Rad-seq or hybrid capture, which dominate phylogenomics today, usually yield a specific set of markers that the researcher has chosen a priori. The generation of WGA thus greatly increases the use of genomic data in biological research, beyond the initial goals of the researcher producing those data. Here, we show how a hypothetical WGA that includes seven different loci (different colors) for four individuals allows extracting sequence data to generate gene trees (B), identifying SNPs to genotype individuals (C), and measuring copy depth to infer CNVs across genomic regions (D). Ultimately, these different kinds of data can be translated into species tree inferences (B -D). In the case of CNVs, only locus number 3 (orange) shows significant CNV. Because CNVs are measured as continuous characters (i.e., copy depth), the orange shading represents a hypothetical evolutionary scenario of copy number variation of genomic region number 3 within the inferred species tree, which is incongruent with those based on sequence and SNP data from other loci in the genome. S1), we observed no increase in the number of species per data set over time (A). On the other hand, there is a significant increase in the number of loci (B), total alignment length (C), and total data set size, as measured by the product of species times locus number (Data set size 1, E) and species times total alignment length (Data set size 2, F). Moreover, the advent of HTS does not support the notion of a tradeoff between the number of species and the number of loci in phylogenomic studies.

Figure 3. Some examples of violations of the multispecies coalescent. In event A, there is gene flow; in event B there is homoploid hybridization; in event C, there is a gene duplication; and in event D, incomplete lineage sorting. All of these processes contribute to gene tree heterogeneity but fall outside the standard multispecies coalescent model. Importantly, all of these processes also yield strictly dichotomous gene trees, whereas recombination (not illustrated here) does not.

Figure 4. Gene duplication and loss creates patterns that can mimic incomplete lineage sorting and can lead to spurious inferences of the species history. Genes and genomes of three species A, B, and C. Multicolored bars show (parts of) their genomes with a number of loci indicated in different colors. The orange gene is duplicated in species A and it was lost in species B. The blue gene was duplicated before the divergence between species A and the ancestor of species B and C. However, one of these copies was lost in species A, whereas both copies were maintained in species B and C. Reconstruction of the orange gene tree based on extant diversity will yield a wrong inference of its history due to the absence of data for species B. On the other hand, a phylogenetic reconstruction of the blue gene is difficult to predict. Depending on which of the duplicates are sampled for species B and C, different outcomes can be expected regarding the relationship among the three species. The duplication and loss history of these two genes may cause serious issues for phylogenetic reconstruction because no specific pattern can be expected between them. Figure 7. Two possible species phylogenies producing similar observations at present time. On the left (A), there is a species tree with gene flow. On the right (B), there is a species network with homoploid hybridization. Distinguishing two such scenarios usually requires simulations and comparison of observed and expected summary statistics.

  and their use in a phylogenomic context. Moreover, single molecule real-time sequencing technologies, such as Pacific Biosciences (Pacific Biosciences of California, Menlo Park, CA)

Figure 2 .

 2 Figure 2. Trends in phylogenomic data sets since the emergence of HTS. Based on a sample of 164 phylogenomic papers published since 2004 (see Supplementary TableS1), we observed no increase in the number of species per data set over time (A). On the other hand, there is a significant increase in the number of loci (B), total alignment length (C), and total data set size, as measured by the product of species times locus number (Data set size 1, E) and species times total alignment length (Data set size 2, F). Moreover, the advent of HTS does not support the notion of a tradeoff between the number of species and the number of loci in phylogenomic studies.

Figure 5 .

 5 Figure 5. Complex patterns of gene lineages with polyploidization and interspecific gene flow. Genes and genomes of four species A, B, C and D. Multi-colored bars show (parts of) genomes with a number of loci indicated in different colors. Two gene trees, one orange and one blue, evolve within the species network. Species B is an allopolyploid containing two genomes.

Figure 6 .

 6 Figure6. Gradual speciation or isolation-with migration. After starting to split, gene flow between species decreases gradually. Such a gradual decrease in the extent of gene flow between species might present an especially useful extension of the standard multispecies coalescent model. Colors depict different gene pools and their gradual change along branches describes how species gradually differentiate despite the existence of migration over time. Thickness and color intensity of arrows show that gene flow becomes weaker as species gradually isolate.

ACKNOWLEDGEMENTS

This paper is a product of the 'Origin of Biodiversity Workshop' organized by Chalmers University of Technology and the University of Gothenburg, under the auspices of the Gothenburg Centre for Advanced Studies (GoCAS). We are particularly grateful to the GoCAS organizers and facilitators, in particular Karin Hårding, Mattias Marklund, Bernt Wennberg, Sandra Johansson, and Lotta Fernström.

We thank Johnathan Clark, Alison Cloutier, Phil Grayson, Kathrin Näpflin, Flavia Termignoni, Jonathan Schmitt, Simon Sin, João Tonini, and Pengcheng Wang for help compiling Supplementary Table 1.

Thomas Couvreur, Tobias Andermann, Prosanta Chakrabarty, Alex Pyron, Claire Morgan, Chris Creevey, and one anonymous reviewer provided useful comments that improved the contents of this manuscript.

approaches have generated well-founded concerns about biases in our inferences [START_REF] Davies | Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism[END_REF][START_REF] Rabosky | No substitute for real data: a cautionary note on the use of phylogenies from birthdeath polytomy resolvers for downstream comparative analyses[END_REF] and the extent of these biases remain largely unknown [START_REF] Pennell | A simple approach for maximizing the overlap of phylogenetic and comparative data[END_REF]. Nonetheless, there are methods that can objectively alleviate issues arising from using incomplete phylogenies in downstream phylogenetic comparative analyses. For example, recent results using conditioned birth-death processes (e.g., Gernhard, 2008a;b;[START_REF] Sagitov | Interspecies correlation for neutrally evolving traits[END_REF] show that under constant rate processes the size of the clade contributes information on the height of the tree and also on the coalescence times. Such results can be used to improve the calibration and node dating of the phylogeny when some species are not sampled. One would expect that ignoring the non-sequenced species would incur a bias resulting in shorter tree heights, because less time is usually required to generate fewer tips. Conditioned branching process models can help alleviate this bias. Also, macroevolutionary birth-death models are used as branching process priors in Bayesian molecular dating.

The availability of likelihood expressions for incompletely sampled phylogenies [START_REF] Stadler | On incomplete sampling under birth-death models and connections to the samplingbased coalescent[END_REF][START_REF] Stadler | Distribution of branch lengths and phylogenetic diversity under homogeneous speciation models[END_REF][START_REF] Morlon | Reconciling molecular phylogenies with the fossil record[END_REF] thus allow to date phylogenies while accounting for the fact that we have observed only a certain fraction of unsampled species.

An important source of genetic data for extinct and unsampled species comes from ancient DNA.

Current methods now allow sequencing ancestral DNA and incorporating it into phylogenomic, phylogeographic, and population genetic analyses (reviewed by [START_REF] Leonardi | Evolutionary Patterns and Processes: Lessons from Ancient DNA[END_REF], and even obtaining whole-genomes from various samples dating back up to a few thousand years (e.g., [START_REF] Lynch | Elephantid Genomes Reveal the Molecular Bases of Woolly Mammoth Adaptations to the Arctic[END_REF][START_REF] Cloutier | First nuclear genome assembly of an extinct moa species, the little bush moa (Anomalopteryx didiformis)[END_REF]. The paleontological and paleobiological records represent another source of extinct species for which genomic data is out of reach. Having the opportunity of integrating the temporal and phylogenetic information of extinct data provided by fossils into phylogenomic estimates of extant diversity is a crucial task toward a complete ToL (reviewed by [START_REF] Hunt | Integrating Paleontological and Phylogenetic Approaches to Macroevolution[END_REF]. Besides the possibility of placing fossil data into phylogenetic trees, a growing field of development is that of timescaling phylogenies using fossil data. Traditionally, node-dating methods have used fossils to bound the minimum age of nodes at the base of extant clades in which fossils likely belong based on morphological