Explore the vast range of titles available.

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 (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.

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome‐scale data from diverse lineages. The development of high‐throughput sequencing and long‐read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project‐dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up‐and‐coming technologies that may help propel the field even further.

Aim Many subtropical organisms exhibit an East Asian‐Tethyan disjunction, a distribution split between East Asia and the Mediterranean. The underlying mechanisms and timing have remained unclear to date. The evolutionary history of Quercus section Ilex Loudon, a representative East Asian‐Tethyan disjunct lineage with a rich and widespread fossil record, was investigated to understand the key drivers of this disjunction. Location Eurasia. Methods The phylogeny of Quercus section Ilex was reconstructed using RAD‐seq. Divergence times were estimated based on three fossil calibrations. Ancestral range and niche were reconstructed on the calibrated tree to infer the timing of transitions in geographic distributions and niche. Convergence in ecological space was estimated by fitting alternative multiple‐regime Ornstein‐Uhlenbeck models. Leaf shape affinities among extant and fossil taxa of section Ilex were assessed using geometric morphometric approaches. Results Six clades were well resolved in section Ilex. Ancestral range reconstruction and divergence time dating suggest a wide distribution along the East Tethys seaway, with initial divergence at the mid‐Eocene, and all six clades originating before the Miocene. The section dispersed from East Asia to the Mediterranean at the Eocene‐Oligocene boundary. A shift toward higher elevations was detected in the Himalayan clade during the middle or late Miocene. European fossil lineages during the early Miocene differ in leaf morphology from later lineages, which we infer to be a consequence of adaptive differentiation or species turnover. Main conclusions Quercus section Ilex was widespread along the East Tethys seaway from the middle Eocene onward. The European holly oaks originated from an East Asian ancestral lineage that dispersed to Europe via the Tibet‐Himalaya corridor in the Oligocene. Lowlands along the margins of the Himalayas and through an Oligocene Tibetan valley served as the dispersal route(s) for these species. Changing climates drove Miocene extinction and local adaptation of European lineages.

• Investigating the interaction between environmental heterogeneity and local adaptation is critical for understanding the evolutionary history of a species, providing the premise for studying the response of organisms to rapid climate change. However, for most species how exactly the spatial heterogeneity promotes population divergence and how genomic variations contribute to adaptive evolution remain poorly understood. • We examine the contributions of geographical and environmental variables to population divergence of the relictual, alpine herb Circaeaster agrestis, as well as the genetic basis of local adaptation using RAD-seq and plastome data. • We detected significant genetic structure with an extraordinary disequilibrium of genetic diversity among regions, and signals of isolation-by-distance along with isolation-by-resistance. The populations were estimated to begin diverging in the late Miocene, along with a possible ancestral distribution of the Hengduan Mountains and adjacent regions. Both environmental gradient and redundancy analyses revealed significant association between genetic variation and temperature variables. Genome–environment association analyses identified 16 putatively adaptive loci related mainly to biotic and abiotic stress resistance. • Our genome-wide data provide new insights into the important role of environmental heterogeneity in shaping genetic structure, and access the footprints of local adaptation in an ancient relictual species, informing future conservation efforts.

Powerful approaches to inferring recent or current population structure based on nearest neighbor haplotype “coancestry” have so far been inaccessible to users without high quality genome-wide haplotype data. With a boom in nonmodel organism genomics, there is a pressing need to bring these methods to communities without access to such data. Here, we present RADpainter, a new program designed to infer the coancestry matrix from restriction-site-associated DNA sequencing (RADseq) data. We combine this program together with a previously published MCMC clustering algorithm into fineRADstructure—a complete, easy to use, and fast population inference package for RADseq data (https://github.com/millanek/fineRADstructure; last accessed February 24, 2018). Finally, with two example data sets, we illustrate its use, benefits, and robustness to missing RAD alleles in double digest RAD sequencing.

During speciation-with-gene-flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome-wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of as isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome-wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation-with-gene-flow mechanism.

Adaptation in the wild often involves standing genetic variation (SGV), which allows rapid responses to selection on ecological timescales. However, we still know little about how the evolutionary histories and genomic distributions of SGV influence local adaptation in natural populations. Here, we address this knowledge gap using the threespine stickleback fish (Gasterosteus aculeatus) as a model. We extend restriction site‐associated DNA sequencing (RAD‐seq) to produce phased haplotypes approaching 700 base pairs (bp) in length at each of over 50,000 loci across the stickleback genome. Parallel adaptation in two geographically isolated freshwater pond populations consistently involved fixation of haplotypes that are identical‐by‐descent. In these same genomic regions, sequence divergence between marine and freshwater stickleback, as measured by dXY, reaches tenfold higher than background levels and genomic variation is structured into distinct marine and freshwater haplogroups. By combining this dataset with a de novo genome assembly of a related species, the ninespine stickleback (Pungitius pungitius), we find that this habitat‐associated divergent variation averages six million years old, nearly twice the genome‐wide average. The genomic variation that is involved in recent and rapid local adaptation in stickleback has therefore been evolving throughout the 15‐million‐year history since the two species lineages split. This long history of genomic divergence has maintained large genomic regions of ancient ancestry that include multiple chromosomal inversions and extensive linked variation. These discoveries of ancient genetic variation spread broadly across the genome in stickleback demonstrate how selection on ecological timescales is a result of genome evolution over geological timescales, and vice versa.

Harnessing the power of genomic scans, we test the debated \"species pump\" hypothesis that implicates repeated cycles of island connectivity and isolation as drivers of divergence. This question has gone understudied given the limited resolution of past molecular markers for studying such dynamic phenomena. With an average of 32,000 SNPs from the genome of 136 individuals from 10 populations of a Caribbean flightless ground cricket species (Amphiacusta sanctaecrucis) and a complementary set of statistical approaches, we infer a stepping-stone colonization model and high levels of genetic differentiation across the Virgin Islands, which have been periodically interconnected until 8 ka. Estimates of divergence times from models based on the site frequency spectrum coincide with a period of repeated connection and fragmentation of the islands at 75–130 ka. These results are consistent with a role of island connectivity cycles in promoting genomic divergence and indicate that the genetic distinctiveness of island populations has persisted despite subsequent and extended interisland connections identified from bathymetric data. We discuss these findings in the broader context of Caribbean biogeography, and more specifically why high levels of genomic divergence across the Virgin Islands associated with repeated connectivity cycles do not actually translate into species diversification.

Advances in sequencing technology provide special opportunities for genotyping individuals with speed and thrift, but the lack of software to automate the calling of tens of thousands of genotypes over hundreds of individuals has hindered progress. Stacks is a software system that uses short-read sequence data to identify and genotype loci in a set of individuals either de novo or by comparison to a reference genome. From reduced representation Illumina sequence data, such as RAD-tags, Stacks can recover thousands of single nucleotide polymorphism (SNP) markers useful for the genetic analysis of crosses or populations. Stacks can generate markers for ultra-dense genetic linkage maps, facilitate the examination of population phylogeography, and help in reference genome assembly. We report here the algorithms implemented in Stacks and demonstrate their efficacy by constructing loci from simulated RAD-tags taken from the stickleback reference genome and by recapitulating and improving a genetic map of the zebrafish, Danio rerio.