SummaryOaks (Quercus, Fagaceae) are the dominant tree genus of North America in species number and biomass, and Mexico is a global center of oak diversity. Understanding the origins of oak diversity is key to understanding biodiversity of northern temperate forests.A phylogenetic study of biogeography, niche evolution and diversification patterns in Quercus was performed using 300 samples, 146 species. Next-generation sequencing data were generated using the restriction-site associated DNA (RAD-seq) method. A timecalibrated maximum likelihood phylogeny was inferred and analyzed with bioclimatic, soils, and leaf habit data to reconstruct the biogeographic and evolutionary history of the American oaks.Our highly resolved phylogeny demonstrates sympatric parallel diversification in climatic niche, leaf habit, and diversification rates. The two major American oak clades arose in what is now the boreal zone and radiated, in parallel, from eastern North America into Mexico and Central America.Oaks adapted rapidly to niche transitions. The Mexican oaks are particularly numerous, not because Mexico is a center of origin, but because of high rates of lineage diversification associated with high rates of evolution along moisture gradients and between the evergreen and deciduous leaf habits. Sympatric parallel diversification in the oaks has shaped the diversity of North American forests.

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny.We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny.Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity.The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

Oaks ( L.) have long been suspected to hybridize in nature, and widespread genetic exchange between morphologically defined species is well documented in two- to six-species systems, but the phylogenetic consequences of hybridization in oaks have never been demonstrated in a phylogenetically diverse sample. Here, we present phylogenomic analyses of a 30 Myr clade that strongly support morphologically defined species and the resolution of novel clades of white oaks; however, historical hybridization across clade boundaries is detectable and, undiagnosed, would obscure the imprint of biogeographic history in the phylogeny. Phylogenetic estimation from restriction-site-associated DNA sequencing data for 156 individuals representing 81 species supports two successive intercontinental disjunctions of white oaks: an early vicariance between the Eurasian and American white oaks, and a second, independent radiation represented by two relictual species. A suite of subsampled and partitioned analyses, however, supports a more recent divergence of the Eurasian white oaks from within the American white oaks and suggests that historic introgression between the Eurasian white oaks and a now-relictual lineage biases concatenated phylogenetic estimates. We demonstrate how divergence and reticulation both influence our understanding of the timing and nature of diversification and global colonization in these ecologically and economically important taxa.

41• The tree of life is highly reticulate, with the history of population divergence buried amongst 42 phylogenies deriving from introgression and lineage sorting. In this study, we test the 43 hypothesis that there are regions of the oak (Quercus, Fagaceae) genome that are broadly 44 informative about phylogeny and investigate global patterns of oak diversity. 45• We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 46 individuals representing ca. 250 oak species to infer a time-calibrated phylogeny of the world's 47• The support we find for the phylogeny contrasts with high genomic heterogeneity in 56 phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may 57 in fact depend on the gene flow that shapes the oak genome. 58 59

Model checking is a critical part of Bayesian data analysis, yet it remains largely unused in systematic studies. Phylogeny estimation has recently moved into an era of increasingly complex models that simultaneously account for multiple evolutionary processes, the statistical fit of these models to the data has rarely been tested. Here we develop a posterior predictive simulation-based model check for a commonly used multispecies coalescent model, implemented in *BEAST, and apply it to 25 published data sets. We show that poor model fit is detectable in the majority of data sets; that this poor fit can mislead phylogenetic estimation; and that in some cases it stems from processes of inherent interest to systematists. We suggest that as systematists scale up to phylogenomic data sets, which will be subject to a heterogeneous array of evolutionary processes, critically evaluating the fit of models to data is an analytical step that can no longer be ignored.

Species within the genus Quercus (oak) hybridize in complex patterns that have yet to be fully explored with phylogenomic data. Analyses to date have recovered reasonable divergent patterns, suggesting that the impact of introgression may not always be obvious in inferred oak phylogenies. We explore this phenomenon using RADseq data for 136 samples representing 54 oak species by conducting phylogenetic analyses designed to distinguish signals of lineage diversification and hybridization, focusing on the lobed-leaf species Quercus gambelii, Q. lobata, and Q. garryana in the context of a broad sampling of allied white oaks (Quercus section Quercus), and particularly the midwestern Q. macrocarpa. We demonstrate that historical introgressive hybridization between once sympatric species affects phylogeny estimation. Historical range expansion during periods of favorable climate likely explains our observations; analyses support genetic exchange between ancestral populations of Q. gambelii and Q. macrocarpa. We conclude that the genomic consequences of introgression caused the attraction of distant lineages in phylogenetic tree space, and that introgressive and divergent signals can be disentangled to produce a robust estimate of the phylogenetic history of the species.

Botanists have long recognised interspecific gene flow as a common occurrence within white oaks (Quercus section Quercus). Historical allele exchange, however, has not been fully characterised and the complex genomic signals resulting from the combination of vertical and horizontal gene transmission may confound phylogenetic inference and obscure our ability to accurately infer the deep evolutionary history of oaks.Using anchored enrichment, we obtained a phylogenomic dataset consisting of hundreds of single-copy nuclear loci. Concatenation, species-tree and network analyses were carried out in an attempt to uncover the genomic signal of ancient introgression and infer the divergent phylogenetic topology for the white oak clade. Locus and site-level likelihood comparisons were then conducted to further explore the introgressed signal within our dataset.Historical, intersectional gene flow is suggested to have occurred between an ancestor of the Eurasian Roburoid lineage and Quercus pontica and North American Dumosae and Prinoideae lineages.Despite extensive time past, our approach proved successful in detecting the genomic signature of ancient introgression. Our results, however, highlight the importance of sampling and the use of a plurality of analytical tools and methods to sufficiently explore genomic datasets, uncover this signal, and accurately infer evolutionary history.

The phylogeny of thamnophiine snakes has not been fully resolved, largely because previous phylogenetic estimates have been based on incomplete taxon sampling or relied solely on mitochondrial sequence data. To address this deficiency, we collected data from multiple autosomal loci collected from 50 taxa before estimating the most robust phylogeny of Thamnophiini to date. Our findings clarify the relationships of taxa not previously included in molecular analyses and also lend evidence to previously recommended taxonomic revisions. Differences in topological estimates between competing models of evolution were minimal and not strongly supported; however, a multispecies coalescent model of evolution was highly favoured over a concatenated model based on marginal likelihood estimates. Additionally, we estimated the timing of divergence among the three major lineages to have occurred during the Miocene period (approximately 11–14 Mya), followed by a decline in speciation rates in all major lineages.