Uncovering the genomic signature of ancient introgression between white oak lineages (Quercus).

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.

Genomic landscape of the global oak phylogeny.

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.

Sympatric parallel diversification of major oak clades in the Americas and the origins of Mexican species diversity.

Oaks (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 time-calibrated 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.

Phylogenomics reveals a complex evolutionary history of lobed-leaf white oaks in western North America.

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.

A genetic legacy of introgression confounds phylogeny and biogeography in oaks.

Oaks (Quercus 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 ca 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.

The evolution and diversification of the red oaks of the California Floristic Province (Quercus section Lobatae, series Agrifoliae).

PREMISE OF THE STUDY: The California Floristic Province (CA-FP) is a unique and diverse region of floral endemism, yet the timing and nature of divergence and diversification of many lineages remain underexplored. We seek to elucidate the evolutionary history of the red oaks of the CA-FP, the Agrifoliae. METHODS: We collected PstI-associated RAD-seq data as well as morphometrics from individuals of the four species across their ranges, including varieties and hybrids. Phylogeny and divergence times were estimated. We analyzed morphological differentiation in over 70 plants using PCA and assessed species delimitation and admixture using genotype clustering analysis in over 40 plants. KEY RESULTS: We find that the Agrifoliae are monophyletic and sister to all other red oak species. Within the Agrifoliae, all species are supported, with Quercus kelloggii sister to a clade of subevergreen taxa: (Quercus agrifolia - (Q. parvula + Q. wislizeni)). Molecular and morphometric analyses are equivocal for named varieties. Notably, Q. parvula var. tamalpaisensis appears to be part of a hybrid swarm between Q. parvula and Q. wislizeni. Dating estimates were concordant with previous hypotheses and geological evidence, with diversification occurring between 10 and 20 million years ago. CONCLUSIONS: The Agrifoliae represent a geographically discrete, early-diverging red oak lineage that diversified during the period of drying and warming associated with Sierran uplift during the middle Miocene. Molecular differentiation within the clade supports the current taxonomy, including an east-west species level pattern (Q. parvula and Q. wislizeni) and north-south intraspecific patterns to some degree, although the latter require additional study.

Poor fit to the multispecies coalescent is widely detectable in empirical data.

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.

Evidence of multiple paternity in Morelet's Crocodile (Crocodylus moreletii) in Belize, CA, inferred from microsatellite markers.

Microsatellite data were generated from hatchlings collected from ten nests of Morelet's Crocodile (Crocodylus moreletii) from New River Lagoon and Gold Button Lagoon in Belize to test for evidence of multiple paternity. Nine microsatellite loci were genotyped for 188 individuals from the 10 nests, alongside 42 nonhatchlings from Gold Button Lagoon. Then mitochondrial control region sequences were generated for the nonhatchlings and for one individual from each nest to test for presence of C. acutus-like haplotypes. Analyses of five of the nine microsatellite loci revealed evidence that progeny from five of the ten nests were sired by at least two males. These data suggest the presence of multiple paternity as a mating strategy in the true crocodiles. This information may be useful in the application of conservation and management techniques to the 12 species in this genus, most of which are threatened or endangered.