Toward a semi-supervised learning approach to phylogenetic estimation.

Models have always been central to inferring molecular evolution and to reconstructing phylogenetic trees. Their use typically involves the development of a mechanistic framework reflecting our understanding of the underlying biological processes, such as nucleotide substitu- tions, and the estimation of model parameters by maximum likelihood or Bayesian inference. However, deriving and optimizing the likelihood of the data is not always possible under complex evolutionary scenarios or even tractable for large datasets, often leading to unrealistic simplifying assumptions in the fitted models. To overcome this issue, we coupled stochastic simulations of genome evolution with a new supervised deep learning model to infer key parameters of molecular evolution. Our model is designed to directly analyze multiple sequence alignments and estimate per-site evolutionary rates and divergence, without requiring a known phylogenetic tree. The accuracy of our predictions matched that of likelihood-based phylogenetic inference, when rate heterogeneity followed a simple gamma distribution, but it strongly exceeded it under more complex patterns of rate variation, such as codon models. Our approach is highly scalable and can be efficiently applied to genomic data, as we showed on a dataset of 26 million nucleotides from the clownfish clade. Our simulations also showed that the integration of per-site rates obtained by deep learning within a Bayesian framework led to significantly more accu- rate phylogenetic inference, particularly with respect to the estimated branch lengths. We thus propose that future advancements in phylogenetic analysis will benefit from a semi-supervised learning approach that combines deep-learning estimation of substitution rates, which allows for more flexible models of rate variation, and probabilistic inference of the phylogenetic tree, which guarantees interpretability and a rigorous assessment of statistical support.

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish.

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex.

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Hyper-Cryptic radiation of a tropical montane plant lineage.

Species are seen as the fundamental unit of biotic diversity, and thus their delimitation is crucial for defining measures for diversity assessments and studying evolution. Differences between species have traditionally been associated with variation in morphology. And yet, the discovery of cryptic diversity suggests that the evolution of distinct lineages does not necessarily involve morphological differences. Here, we analyze 1,684,987 variant sites and over 4,000 genes for more than 400 samples to show how a tropical montane plant lineage (Geonoma undata species complex) is composed of numerous unrecognized genetic groups that are not morphologically distinct. We find that 11 to 14 clades do not correspond to the three currently recognized species. Most clades are genetically different and geographic distance and topography are the most important factors determining this genetic divergence. The genetic structure of this lineage does not match its morphological variation. Instead, this species complex constitutes the first example of a hyper-cryptic plant radiation in tropical mountains.

Conceptual and empirical bridges between micro- and macroevolution.

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

Insights into the Genomics of Clownfish Adaptive Radiation: The Genomic Substrate of the Diversification.

Clownfishes are an iconic group of coral reef fishes that evolved a mutualistic interaction with sea anemones, which triggered the rapid diversification of the group. Following the emergence of this mutualism, clownfishes diversified into different ecological niches and developed convergent phenotypes associated with their host use. The genetic basis of the initial acquisition of the mutualism with host anemones has been described, but the genomic architecture underlying clownfish diversification once the mutualism was established and the extent to which clownfish phenotypic convergence originated through shared genetic mechanisms are still unknown. Here, we investigated these questions by performing comparative genomic analyses on the available genomic data of five pairs of closely related but ecologically divergent clownfish species. We found that clownfish diversification was characterized by bursts of transposable elements, an overall accelerated coding evolution, incomplete lineage sorting, and ancestral hybridization events. Additionally, we detected a signature of positive selection in 5.4% of the clownfish genes. Among them, five presented functions associated with social behavior and ecology, and they represent candidate genes involved in the evolution of the size-based hierarchical social structure so particular to clownfishes. Finally, we found genes with patterns of either relaxation or intensification of purifying selection and signals of positive selection linked with clownfish ecological divergence, suggesting some level of parallel evolution during the diversification of the group. Altogether, this work provides the first insights into the genomic substrate of clownfish adaptive radiation and integrates the growing collection of studies investigating the genomic mechanisms governing species diversification.

Large-scale whole-genome resequencing unravels the domestication history of Cannabis sativa.

Cannabis sativa has long been an important source of fiber extracted from hemp and both medicinal and recreational drugs based on cannabinoid compounds. Here, we investigated its poorly known domestication history using whole-genome resequencing of 110 accessions from worldwide origins. We show that C. sativa was first domesticated in early Neolithic times in East Asia and that all current hemp and drug cultivars diverged from an ancestral gene pool currently represented by feral plants and landraces in China. We identified candidate genes associated with traits differentiating hemp and drug cultivars, including branching pattern and cellulose/lignin biosynthesis. We also found evidence for loss of function of genes involved in the synthesis of the two major biochemically competing cannabinoids during selection for increased fiber production or psychoactive properties. Our results provide a unique global view of the domestication of C. sativa and offer valuable genomic resources for ongoing functional and molecular breeding research.

Evolution of life cycles and reproductive traits: Insights from the brown algae.

A vast diversity of types of life cycles exists in nature, and several theories have been advanced to explain how this diversity has evolved and how each type of life cycle is retained over evolutionary time. Here, we exploited the diversity of life cycles and reproductive traits of the brown algae (Phaeophyceae) to test several hypotheses on the evolution of life cycles. We investigated the evolutionary dynamics of four life-history traits: life cycle, sexual system, level of gamete dimorphism and gamete parthenogenetic capacity. We assigned states to up to 77 representative species of the taxonomic diversity of the brown algal group, in a multi-gene phylogeny. We used maximum likelihood and Bayesian analyses of correlated evolution, while taking the phylogeny into account, to test for correlations between traits and to investigate the chronological sequence of trait acquisition. Our analyses are consistent with the prediction that diploid growth evolves when sexual reproduction is preferred over asexual reproduction, possibly because it allows the complementation of deleterious mutations. We also found that haploid sex determination is ancestral in relation to diploid sex determination. However, our results could not address whether increased zygotic and diploid growth are associated with increased sexual dimorphism. Our analyses suggest that in the brown algae, isogamous species evolved from anisogamous ancestors, contrary to the commonly reported pattern where evolution proceeds from isogamy to anisogamy.

Fossil-Informed Models Reveal a Boreotropical Origin and Divergent Evolutionary Trajectories in the Walnut Family (Juglandaceae).

Temperate woody plants in the Northern Hemisphere have long been known to exhibit high species richness in East Asia and North America and significantly lower diversity in Europe, but the causes of this pattern remain debated. Here, we quantify the roles of dispersal, niche evolution, and extinction in shaping the geographic diversity of the temperate woody plant family Juglandaceae (walnuts and their relatives). Integrating evidence from molecular, morphological, fossil, and (paleo)environmental data, we find strong support for a Boreotropical origin of the family with contrasting evolutionary trajectories between the temperate subfamily Juglandoideae and the tropical subfamily Engelhardioideae. Juglandoideae rapidly evolved frost tolerance when the global climate shifted to ice-house conditions from the Oligocene, with diversification at high latitudes especially in Europe and Asia during the Miocene. Subsequent range contraction at high latitudes and high levels of extinction in Europe driven by global cooling led to the current regional disparity in species diversity. Engelhardioideae showed temperature conservatism while adapting to increased humidity, tracking tropical climates to low latitudes since the middle Eocene with comparatively little diversification, perhaps due to high competition in the tropical zone. The biogeographic history of Juglandaceae shows that the North Atlantic land bridge and Europe played more critical roles than previously thought in linking the floras of East Asia and North America, and showcases the complex interplay among climate change, niche evolution, dispersal, and extinction that shaped the modern disjunct pattern of species richness in temperate woody plants. [Boreotropical origin; climatic niche evolution; disjunct distribution; dispersal; diversity anomaly; extinction; Juglandaceae.].

Thyroid hormones regulate the formation and environmental plasticity of white bars in clownfishes.

Determining how plasticity of developmental traits responds to environmental conditions is a challenge that must combine evolutionary sciences, ecology, and developmental biology. During metamorphosis, fish alter their morphology and color pattern according to environmental cues. We observed that juvenile clownfish (Amphiprion percula) modulate the developmental timing of their adult white bar formation during metamorphosis depending on the sea anemone species in which they are recruited. We observed an earlier formation of white bars when clownfish developed with Stichodactyla gigantea (Sg) than with Heteractis magnifica (Hm). As these bars, composed of iridophores, form during metamorphosis, we hypothesized that timing of their development may be thyroid hormone (TH) dependent. We treated clownfish larvae with TH and found that white bars developed earlier than in control fish. We further observed higher TH levels, associated with rapid white bar formation, in juveniles recruited in Sg than in Hm, explaining the faster white bar formation. Transcriptomic analysis of Sg recruits revealed higher expression of duox, a dual oxidase implicated in TH production as compared to Hm recruits. Finally, we showed that duox is an essential regulator of iridophore pattern timing in zebrafish. Taken together, our results suggest that TH controls the timing of adult color pattern formation and that shifts in duox expression and TH levels are associated with ecological differences resulting in divergent ontogenetic trajectories in color pattern development.

Phylogenomics of Gesneriaceae using targeted capture of nuclear genes.

Gesneriaceae (ca. 3400 species) is a pantropical plant family with a wide range of growth form and floral morphology that are associated with repeated adaptations to different environments and pollinators. Although Gesneriaceae systematics has been largely improved by the use of Sanger sequencing data, our understanding of the evolutionary history of the group is still far from complete due to the limited number of informative characters provided by this type of data. To overcome this limitation, we developed here a Gesneriaceae-specific gene capture kit targeting 830 single-copy loci (776,754 bp in total), including 279 genes from the Universal Angiosperms-353 kit. With an average of 557,600 reads and 87.8% gene recovery, our target capture was successful across the family Gesneriaceae and also in other families of Lamiales. From our bait set, we selected the most informative 418 loci to resolve phylogenetic relationships across the entire Gesneriaceae family using maximum likelihood and coalescent-based methods. Upon testing the phylogenetic performance of our baits on 78 taxa representing 20 out of 24 subtribes within the family, we showed that our data provided high support for the phylogenetic relationships among the major lineages, and were able to provide high resolution within more recent radiations. Overall, the molecular resources we developed here open new perspectives for the study of Gesneriaceae phylogeny at different taxonomical levels and the identification of the factors underlying the diversification of this plant group.

On the Effect of Asymmetrical Trait Inheritance on Models of Trait Evolution.

Current phylogenetic comparative methods modeling quantitative trait evolution generally assume that, during speciation, phenotypes are inherited identically between the two daughter species. This, however, neglects the fact that species consist of a set of individuals, each bearing its own trait value. Indeed, because descendent populations after speciation are samples of a parent population, we can expect their mean phenotypes to randomly differ from one another potentially generating a "jump" of mean phenotypes due to asymmetrical trait inheritance at cladogenesis. Here, we aim to clarify the effect of asymmetrical trait inheritance at speciation on macroevolutionary analyses, focusing on model testing and parameter estimation using some of the most common models of quantitative trait evolution. We developed an individual-based simulation framework in which the evolution of phenotypes is determined by trait changes at the individual level accumulating across generations, and cladogenesis occurs then by separation of subsets of the individuals into new lineages. Through simulations, we assess the magnitude of phenotypic jumps at cladogenesis under different modes of trait inheritance at speciation. We show that even small jumps can strongly alter both the results of model selection and parameter estimations, potentially affecting the biological interpretation of the estimated mode of evolution of a trait. Our results call for caution when interpreting analyses of trait evolution, while highlighting the importance of testing a wide range of alternative models. In the light of our findings, we propose that future methodological advances in comparative methods should more explicitly model the intraspecific variability around species mean phenotypes and how it is inherited at speciation.

A Cautionary Note on the Use of Genotype Callers in Phylogenomics.

Next-generation-sequencing genotype callers are commonly used in studies to call variants from newly sequenced species. However, due to the current availability of genomic resources, it is still common practice to use only one reference genome for a given genus, or even one reference for an entire clade of a higher taxon. The problem with traditional genotype callers, such as the one from GATK, is that they are optimized for variant calling at the population level. However, when these callers are used at the phylogenetic level, the consequences for downstream analyses can be substantial. Here, we performed simulations to compare the performance between the genotype callers of GATK and ATLAS, and present their differences at various phylogenetic scales. We show that the genotype caller of GATK substantially underestimates the number of variants at the phylogenetic level, but not at the population level. We also found that the accuracy of heterozygote calls declines with increasing distance to the reference genome. We quantified this decline and found that it is very sharp in GATK, while ATLAS maintains high accuracy even at moderately divergent species from the reference. We further suggest that efforts should be taken towards acquiring more reference genomes per species, before pursuing high-scale phylogenomic studies. [ATLAS; efficiency of SNP calling; GATK; heterozygote calling; next-generation sequencing; reference genome; variant calling.].

Rapid climate change results in long-lasting spatial homogenization of phylogenetic diversity.

Scientific understanding of biodiversity dynamics, resulting from past climate oscillations and projections of future changes in biodiversity, has advanced over the past decade. Little is known about how these responses, past or future, are spatially connected. Analyzing the spatial variability in biodiversity provides insight into how climate change affects the accumulation of diversity across space. Here, we evaluate the spatial variation of phylogenetic diversity of European seed plants among neighboring sites and assess the effects of past rapid climate changes during the Quaternary on these patterns. Our work shows a marked homogenization in phylogenetic diversity across Central and Northern Europe linked to high climate change velocity and large distances to refugia. Our results suggest that the future projected loss in evolutionary heritage may be even more dramatic, as homogenization in response to rapid climate change has occurred among sites across large landscapes, leaving a legacy that has lasted for millennia.

Population Genetic Structure and Demographic History of Primula fasciculata in Southwest China.

Understanding the factors that drive the genetic structure of a species and its responses to past climatic changes is an important first step in modern population management. The response to the last glacial maximum (LGM) has been well studied, however, the effect of previous glaciation periods on plant demographic history is still not well studied. Here we investigated the population structure and demographic history of Primula fasciculata that widely occurs in the Hengduan Mountains and Qinghai-Tibetan Plateau. We obtained genomic data for 234 samples of the species using restriction site-associated DNA (RAD) sequencing and combined approximate Bayesian computation (ABC) and species distribution modeling (SDM) to evaluate the effects of multiple glaciation periods by testing several population divergence models and demographic scenarios. The analyses of population structure showed that P. fasciculata displays a striking population structure with six groups that could be identified genetically. Our ABC modeling suggested that the current groups diverged from ancestral populations located in the eastern Hengduan Mountains after the largest glaciation occurred in the region (~ 0.8-0.5 million years ago), which is consistent with the result of SDMs. Each current group has survived in different glacial refugia during the LGM and experienced expansions and/or bottlenecks since their divergence during or across the following Quaternary glacial cycles. Our study demonstrates the usefulness of population genomics for evaluating the effects of past climatic changes in alpine plant species with shallow population structure.

Genomic footprints of repeated evolution of CAM photosynthesis in a Neotropical species radiation.

The adaptive radiation of Bromeliaceae (pineapple family) is one of the most diverse among Neotropical flowering plants. Diversification in this group was facilitated by shifts in several adaptive traits or "key innovations" including the transition from C3 to CAM photosynthesis associated with xeric (heat/drought) adaptation. We used phylogenomic approaches, complemented by differential gene expression (RNA-seq) and targeted metabolite profiling, to address the mechanisms of C3 /CAM evolution in the extremely species-rich bromeliad genus, Tillandsia, and related taxa. Evolutionary analyses of whole-genome sequencing and RNA-seq data suggest that evolution of CAM is associated with coincident changes to different pathways mediating xeric adaptation in this group. At the molecular level, C3 /CAM shifts were accompanied by gene expansion of XAP5 CIRCADIAN TIMEKEEPER homologs, a regulator involved in sugar- and light-dependent regulation of growth and development. Our analyses also support the re-programming of abscisic acid-related gene expression via differential expression of ABF2/ABF3 transcription factor homologs, and adaptive sequence evolution of an ENO2/LOS2 enolase homolog, effectively tying carbohydrate flux to abscisic acid-mediated abiotic stress response. By pinpointing different regulators of overlapping molecular responses, our results suggest plausible mechanistic explanations for the repeated evolution of correlated adaptive traits seen in a textbook example of an adaptive radiation.

Prototypic SNARE Proteins Are Encoded in the Genomes of Heimdallarchaeota, Potentially Bridging the Gap between the Prokaryotes and Eukaryotes.

A defining feature of eukaryotic cells is the presence of numerous membrane-bound organelles that subdivide the intracellular space into distinct compartments. How the eukaryotic cell acquired its internal complexity is still poorly understood. Material exchange among most organelles occurs via vesicles that bud off from a source and specifically fuse with a target compartment. Central players in the vesicle fusion process are the soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins. These small tail-anchored (TA) membrane proteins zipper into elongated four-helix bundles that pull membranes together. SNARE proteins are highly conserved among eukaryotes but are thought to be absent in prokaryotes. Here, we identified SNARE-like factors in the genomes of uncultured organisms of Asgard archaea of the Heimdallarchaeota clade, which are thought to be the closest living relatives of eukaryotes. Biochemical experiments show that the archaeal SNARE-like proteins can interact with eukaryotic SNARE proteins. We did not detect SNAREs in α-proteobacteria, the closest relatives of mitochondria, but identified several genes encoding for SNARE proteins in γ-proteobacteria of the order Legionellales, pathogens that live inside eukaryotic cells. Very probably, their SNAREs stem from lateral gene transfer from eukaryotes. Together, this suggests that the diverse set of eukaryotic SNAREs evolved from an archaeal precursor. However, whether Heimdallarchaeota actually have a simplified endomembrane system will only be seen when we succeed studying these organisms under the microscope.

Evolution of a supergene that regulates a trans-species social polymorphism.

Supergenes are clusters of linked genetic loci that jointly affect the expression of complex phenotypes, such as social organization. Little is known about the origin and evolution of these intriguing genomic elements. Here we analyse whole-genome sequences of males from native populations of six fire ant species and show that variation in social organization is under the control of a novel supergene haplotype (termed Sb), which evolved by sequential incorporation of three inversions spanning half of a 'social chromosome'. Two of the inversions interrupt protein-coding genes, resulting in the increased expression of one gene and modest truncation in the primary protein structure of another. All six socially polymorphic species studied harbour the same three inversions, with the single origin of the supergene in their common ancestor inferred by phylogenomic analyses to have occurred half a million years ago. The persistence of Sb along with the ancestral SB haplotype through multiple speciation events provides a striking example of a functionally important trans-species social polymorphism presumably maintained by balancing selection. We found that while recombination between the Sb and SB haplotypes is severely restricted in all species, a low level of gene flux between the haplotypes has occurred following the appearance of the inversions, potentially mitigating the evolutionary degeneration expected at genomic regions that cannot freely recombine. These results provide a detailed picture of the structural genomic innovations involved in the formation of a supergene controlling a complex social phenotype.

Slowly but surely: gradual diversification and phenotypic evolution in the hyper-diverse tree fern family Cyatheaceae.

BACKGROUND AND AIMS: The tremendously unbalanced distribution of species richness across clades in the tree of life is often interpreted as the result of variation in the rates of diversification, which may themselves respond to trait evolution. Even though this is likely a widespread pattern, not all diverse groups of organisms exhibit heterogeneity in their dynamics of diversification. Testing and characterizing the processes driving the evolution of clades with steady rates of diversification over long periods of time are of importance in order to have a full understanding of the build-up of biodiversity through time. METHODS: We studied the macroevolutionary history of the species-rich tree fern family Cyatheaceae and inferred a time-calibrated phylogeny of the family including extinct and extant species using the recently developed fossilized birth-death method. We tested whether the high diversity of Cyatheaceae is the result of episodes of rapid diversification associated with phenotypic and ecological differentiation or driven by stable but low rates of diversification. We compared the rates of diversification across clades, modelled the evolution of body size and climatic preferences and tested for trait-dependent diversification. KEY RESULTS: This ancient group diversified at a low and constant rate during its long evolutionary history. Morphological and climatic niche evolution were found to be overall highly conserved, although we detected several shifts in the rates of evolution of climatic preferences, linked to changes in elevation. The diversification of the family occurred gradually, within limited phenotypic and ecological boundaries, and yet resulted in a remarkable species richness. CONCLUSIONS: Our study indicates that Cyatheaceae is a diverse clade which slowly accumulated morphological, ecological and taxonomic diversity over a long evolutionary period and provides a compelling example of the tropics as a museum of biodiversity.

Linking micro and macroevolution in the presence of migration.

Understanding macroevolutionary patterns is central to evolutionary biology. This involves the process of divergence within a species, which starts at the microevolutionary level, for instance, when two subpopulations evolve towards different phenotypic optima. The speed at which these optima are reached is controlled by the degree of stabilising selection, which pushes the mean trait towards different optima in the different subpopulations, and ongoing migration that pulls the mean phenotype away from that optimum. Traditionally, macro phenotypic evolution is modelled by directional selection processes, but these models usually ignore the role of migration within species. Here, our goal is to reconcile the processes of micro and macroevolution by modelling migration as part of the speciation process. More precisely, we introduce an Ornstein-Uhlenbeck (OU) model where migration happens between two subpopulations within a branch of a phylogeny and this migration decreases over time as it happens during speciation. We then use this model to study the evolution of trait means along a phylogeny, as well as the way phenotypic disparity between species changes with successive epochs. We show that ignoring the effect of migration in sampled time-series data biases significantly the estimation of the selective forces acting upon it. We also show that migration decreases the expected phenotypic disparity between species and we analyse the effect of migration in the particular case of niche filling. We further introduce a method to jointly estimate selection and migration from time-series data. Our model extends traditional quantitative genetics results of selection and migration from a microevolutionary time frame to multiple speciation events at a macroevolutionary scale. Our results further support that not accounting for gene flow has important consequences in inferences at both the micro and macroevolutionary scale.