Demographic history of speciation in a Senecio altitudinal hybrid zone on Mt. Etna

Weak Genetic Isolation and Putative Phenotypic Selection in the Wild Carnation Dianthus virgineus (Caryophyllaceae)

By relating genetic divergence at neutral loci, phenotypic variation, and geographic and environmental distances, it is possible to dissect micro-evolutionary scenarios involving natural selection and neutral evolution. In this work, we tested the patterns of intraspecific genetic and phenotypic variation along an elevational gradient, using Dianthus virgineus as study system. We genotyped genome-wide SNPs through ddRAD sequencing and quantified phenotypic variation through multivariate morphological variation. We assessed patterns of variation by testing the statistical association between genetic, phenotypic, geographic, and elevational distances and explored the role of genetic drift and selection by comparing the Fst and Pst of morphometric traits. We revealed a weak genetic structure related to geographic distance among populations, but we excluded the predominant role of genetic drift acting on phenotypic traits. A high degree of phenotypic differentiation with respect to genetic divergence at neutral loci allowed us to hypothesize the effect of selection, putatively fuelled by changing conditions at different sites, on morphological traits. Thus, natural selection acting despite low genetic divergence at neutral loci can be hypothesized as a putative driver explaining the observed patterns of variation.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Rapid evolution of hybrid breakdown following recent divergence with gene flow in Senecio species on Mount Etna, Sicily

How do nascent species evolve reproductive isolation during speciation with on-going gene flow? How do hybrid lineages become stabilised hybrid species? While commonly used genomic approaches provide an indirect way to identify species incompatibility factors, synthetic hybrids generated from interspecific crosses allow direct pinpointing of phenotypic traits involved in incompatibilities and the traits that are potentially adaptive in hybrid species. Here we report the analysis of phenotypic variation and hybrid breakdown in crosses between closely-related Senecio aethnensis and S. chrysanthemifolius, and their homoploid hybrid species, S. squalidus. The two former species represent a likely case of recent (<200 ky) speciation with gene flow driven by adaptation to contrasting conditions of high- and low-elevations on Mount Etna, Sicily. As these species form viable and fertile hybrids, it remains unclear whether they have started to evolve reproductive incompatibility. Our analysis represents the first study of phenotypic variation and hybrid breakdown involving multiple Senecio hybrid families. It revealed wide range of variation in multiple traits, including the traits previously unrecorded in synthetic hybrids. Leaf shape, highly distinct between S. aethnensis and S. chrysanthemifolius, was extremely variable in F₂ hybrids, but more consistent in S. squalidus. Our study demonstrates that interspecific incompatibilities can evolve rapidly despite on-going gene flow between the species. Further work is necessary to understand the genetic bases of these incompatibilities and their role in speciation with gene flow.

Whole-genome analysis of multiple wood ant population pairs supports similar speciation histories, but different degrees of gene flow, across their European ranges

The application of demographic history modelling and inference to the study of divergence between species has become a cornerstone of speciation genomics. Speciation histories are usually reconstructed by analysing single populations from each species, assuming that the inferred population history represents the actual speciation history. However, this assumption may not be met when species diverge with gene flow, for example, when secondary contact may be confined to specific geographic regions. Here, we tested whether divergence histories inferred from heterospecific populations may vary depending on their geographic locations, using the two wood ant species Formica polyctena and F. aquilonia. We performed whole‐genome resequencing of 20 individuals sampled in multiple locations across the European ranges of both species. Then, we reconstructed the histories of distinct heterospecific population pairs using a coalescent‐based approach. Our analyses always supported a scenario of divergence with gene flow, suggesting that divergence started in the Pleistocene (c. 500 kya) and occurred with continuous asymmetrical gene flow from F. aquilonia to F. polyctena until a recent time, when migration became negligible (2–19 kya). However, we found support for contemporary gene flow in a sympatric pair from Finland, where the species hybridise, but no signature of recent bidirectional gene flow elsewhere. Overall, our results suggest that divergence histories reconstructed from a few individuals may be applicable at the species level. Nonetheless, the geographical context of populations chosen to represent their species should be taken into account, as it may affect estimates of migration rates between species when gene flow is spatially heterogeneous.

An integrative genomic and phenomic analysis to investigate the nature of plant species in Escallonia (Escalloniaceae)

What we mean by species and whether they have any biological reality has been debated since the early days of evolutionary biology. Some biologists even suggest that plant species are created by taxonomists as a subjective, artificial division of nature. However, the nature of plant species has been rarely tested critically with data while ignoring taxonomy. We integrate phenomic and genomic data collected across hundreds of individuals at a continental scale to investigate this question in Escallonia (Escalloniaceae), a group of plants which includes 40 taxonomic species (the species proposed by taxonomists). We first show that taxonomic species may be questionable as they match poorly to patterns of phenotypic and genetic variation displayed by individuals collected in nature. We then use explicit statistical methods for species delimitation designed for phenotypic and genomic data, and show that plant species do exist in Escallonia as an objective, discrete property of nature independent of taxonomy. We show that such species correspond poorly to current taxonomic species (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%) and that phenomic and genomic data seldom delimit congruent entities (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%). These discrepancies suggest that evolutionary forces additional to gene flow can maintain the cohesion of species. We propose that phenomic and genomic data analyzed on an equal footing build a broader perspective on the nature of plant species by helping delineate different ‘types of species’. Our results caution studies which take the accuracy of taxonomic species for granted and challenge the notion of plant species without empirical evidence. Note: A version of the complete manuscript in Spanish is available in the Supplemental Materials.

Evidence of altitudinal gradient modifying genomic and chemical diversity in populations of Lychnophora pinaster Mart

Lychnophora pinaster Mart. (Asteraceae) is endemic to the Brazilian Cerrado. It is distributed along the altitudinal gradient of the mountainous ranges of the state of Minas Gerais. This study aimed to evaluate the influence of altitude on the genetic diversity of L. pinaster populations and the effects of altitude and climatic factors on essential oil chemical composition. Essential oils from L. pinaster populations from the north (North 01, North 02, and North 03, 700-859 m) and the Metropolitan region of Belo Horizonte (MhBH 01 and MrBH 02, 1366-1498 m) were analyzed. SNP markers from L. pinaster in these regions and Campos das Vertentes (CV 01, CV 02, and CV 03, 1055-1292 m) were also analyzed. The main compounds in essential oils were 14-hydroxy-α-humulene (North 01 and North 03), cedr-8(15)-en-9-α-ol (North 02), 14-acetoxy-α-humulene (MrBH 01), and 4-oxo-15-nor-eudesman-11-ene (MrBH 02). Hierarchical cluster and heatmap analyses showed that the North and MrBH populations included five different groups, indicating the chemical composition of essential oils is distinct in each population. Furthermore, principal component analysis showed that higher altitudes (1366 m and 1498 m) in the MrBH influence the chemical composition of essential oils, and climatic factors determine the chemical composition in North region. The genetic diversity showed that most alleles are in Hardy-Weinberg equilibrium and imply high genetic variation and genetic polymorphisms between populations. Furthermore, the results of Mantel tests (R = 0.3861517; p = 0.04709529; R = 0.9423121; p = 0.02739726) also showed that higher altitude (>1360 m) shapes the genetic diversity at the MrBH. The genetic structure showed that higher altitudes (>1360 m) contribute to the structure of the MrBH populations, but not to North and CV populations. Therefore, the altitudinal ranges of Minas Gerais mountainous ranges determine the higher genetic and chemical diversity of L. pinaster populations.

Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: Current knowledge and future perspectives

Mountain ecosystems support a significant one-third of all terrestrial biodiversity, but our understanding of the spatiotemporal maintenance of this high biodiversity remains poor, or at best controversial. The Himalaya hosts a complex mountain ecosystem with high topographic and climatic heterogeneity and harbors one of the world's richest floras. The high species endemism, together with increasing anthropogenic threats, has qualified the Himalaya as one of the most significant global biodiversity hotspots. The topographic and climatic complexity of the Himalaya makes it an ideal natural laboratory for studying the mechanisms of floral exchange, diversification, and spatiotemporal distributions. Here, we review literature pertaining to the Himalaya in order to generate a concise synthesis of the origin, distribution, and climate change responses of the Himalayan flora. We found that the Himalaya supports a rich biodiversity and that the Hengduan Mountains supplied the majority of the Himalayan floral elements, which subsequently diversified from the late Miocene onward, to create today's relatively high endemicity in the Himalaya. Further, we uncover links between this Miocene diversification and the joint effect of geological and climatic upheavals in the Himalaya. There is marked variance regarding species dispersal, elevational gradients, and impact of climate change among plant species in the Himalaya, and our review highlights some of the general trends and recent advances on these aspects. Finally, we provide some recommendations for conservation planning and future research. Our work could be useful in guiding future research in this important ecosystem and will also provide new insights into the maintenance mechanisms underpinning other mountain systems.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Molecular phylogeny, biogeography and character evolution of the montane genus Incarvillea Juss. (Bignoniaceae)

The complex orogeny of the Himalaya and the Qinghai-Tibet Plateau (QTP) fosters habitat fragmentation that drives morphological differentiation of mountain plant species. Consequently, determining phylogenetic relationships between plant subgenera using morphological characters is unreliable. Therefore, we used both molecular phylogeny and historical biogeographic analysis to infer the ancestral states of several vegetative and reproductive characters of the montane genus Incarvillea. We determined the taxonomic position of the genus Incarvillea within its family and inferred the biogeographical origin of taxa through Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP) analyses using three molecular data sets (trnL-trnF sequences, nr ITS sequences, and a data set of combined sequences) derived from 81% of the total species of the genus Incarvillea. Within the genus-level phylogenetic framework, we examined the character evolution of 10 key morphological characters, and inferred the ancestral area and biogeographical history of the genus. Our analyses revealed that the genus Incarvillea is monophyletic and originated in Central Asia during mid-Oligocene ca. 29.42 Ma. The earliest diverging lineages were subsequently split into the Western Himalaya and Sino-Himalaya during the early Miocene ca. 21.12 Ma. These lineages resulted in five re-circumscribed subgenera (Amphicome, Olgaea, Niedzwedzkia, Incarvillea, and Pteroscleris). Moreover, character mapping revealed the ancestral character states of the genus Incarvillea (e.g., suffruticose habit, cylindrical capsule shape, subligneous capsule texture, absence of capsule wing, and loculicidal capsule dehiscence) that are retained at the earliest diverging ancestral nodes across the genus. Our phylogenetic tree of the genus Incarvillea differs from previously proposed phylogenies, thereby recommending the placement of the subgenus Niedzwedzkia close to the subgenus Incarvillea and maintaining two main divergent lineages.

Hybridization and cryptic speciation in the Iberian endemic plant genus Phalacrocarpum (Asteraceae-Anthemideae)

Understanding the role and impact of reticulation in phylogenetic inquiry has improved with extended use of high throughput sequencing data. Yet, due to the dynamism of genomes over evolutionary time, disentangling old hybridization events remains a serious challenge. Phalacrocarpum (DC.) Willk. is one of the 27 Iberian endemic plant genera, currently considered monotypic but including three subspecies. Its uncertain phylogenetic relationships within tribe Anthemideae (Asteraceae) point to an Early Miocene divergence from its sister group, and its persistent taxonomic instability has been proposed to be due to hybridization. We aim at understanding the evolutionary history of this genus using SNPs called from a genotyping-by-sequencing (GBS) analysis, Sanger sequences-from three plastid DNA regions (psbJ-petA, petB-petD, trnH-psbA) and the nuclear ribosomal ITS regions (cloned)-as well as leaf morphometric multivariate analysis. SNP data and Sanger sequences strongly support the unforeseen existence of a cryptic species in the eastern populations of P. oppositifolium subsp. anomalum. Broad molecular and morphometric patterns of variation found in conflictive populations from the Sanabria Valley region convincingly identify a recent previously undocumented hybrid zone. By contrast, evidence is less conclusive on relationships between subspecies hoffmannseggii, oppositifolium and a second conflictive group distributed along the Galician-Portuguese border (Orense massifs). Although genetic clustering analysis of SNP data suggests that the former subspecies was the maternal progenitor in hybridization events that gave rise to the other two groups, we found considerable uniqueness of ITS ribotypes and plastid haplotypes in them. This result, in the context of Pleistocene climatically-driven range shifts in NW Iberian Peninsula, can be due to periods of isolation, genetic bottlenecks and drift superimposed on old hybridization events. Our study confirms the idea that unravelling old hybridization events may be compromised by the suite of evolutionary processes accumulated subsequently, particularly in areas with a history of climatic instability.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

Genetic differentiation in Sichuan jay (Perisoreus internigrans) and its sibling species Siberian jay (P. infaustus)

The Pleistocene ice age and recent forest fragmentation have both played a significant role in shaping the population genetic variation of boreal coniferous species in the Qinghai-Tibet Plateau (QTP) and Eurasian coniferous forests. The Sichuan jay is one of the least know endemic bird species in QTP coniferous forests in western China while its sibling species, the Siberian Jay, is widespread within the coniferous forests in northern Eurasia. Here we used 11 microsatellite markers to assay genetic diversity across 58 Sichuan jay samples from China and 205 Siberian jay samples from Sweden and Russia. Results showed three distinct genetic clusters from the Sichuan jay sampling. Furthermore, the pair-wise FST values indicated high genetic differentiation not only among the two species but also between Swedish and Russian Siberian jay populations. What is more, a pattern of isolation by distance was found among the analyzed populations. Our study suggests that targeted habitat restoration in fragmented forests and more genetic work is urgently needed for conservation of the Sichuan jay.

Demographic Inference of Divergence and Gene Exchange Between Castanopsis fabri and Castanopsis lamontii

The cytoplasmic genome of one species may be replaced by that of another species without leaving any trace of past hybridization in its nuclear genome, which can thus confuse the inference of genealogical relationship and evolutionary history of many congeneric species. In this study, we used sequence variations of chloroplast DNA and restriction site-associated DNA to investigate gene exchange between Castanopsis fabri and Castanopsis lamontii, and to infer the divergence history of the two species by comparing different divergence models based on the joint allele frequency spectrum. We evaluated climatic niche similarity of the two species using climatic variables across their entire distribution range in subtropical China. Clear genetic differentiation was revealed between C. fabri and C. lamontii, and gene exchange between the two species was discovered as a consequence of secondary contact. The gene exchange rates were variable across the genome. Gene exchange could allow C. fabri to widen its habitat through pollen swamping and broaden its climatic niche, and the chloroplast genome of C. lamontii is captured by C. fabri during this process. These results further our understanding of the timing and contribution of gene exchange to species divergence in forests.

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome-wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly-differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs' allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16-0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species - pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.

Incomplete reproductive isolation between Rhododendron taxa enables hybrid formation and persistence

The evolutionary consequences of hybridization ultimately depend on the magnitude of reproductive isolation between hybrids and their parents. We evaluated the relative contributions of pre- and post-zygotic barriers to reproduction for hybrid formation, hybrid persistence and potential for reproductive isolation of hybrids formed between two Rhododendron species, R. spiciferum and R. spinuliferum. Our study established that incomplete reproductive isolation promotes hybrid formation and persistence and delays hybrid speciation. All pre-zygotic barriers to reproduction leading to hybrid formation are incomplete: parental species have overlapping flowering; they share the same pollinators; reciprocal assessments of pollen tube germination and growth do not differ among parents. The absence of post-zygotic barriers between parental taxa indicates that the persistence of hybrids is likely. Reproductive isolation was incomplete between hybrids and parents in all cases studied, although asymmetric differences in reproductive fitness were prevalent and possibly explain the genetic structure of natural hybrid swarms where hybridization is known to be bidirectional but asymmetric. Introgression, rather than speciation, is a probable evolutionary outcome of hybridization between the two Rhododendron taxa. Our study provides insights into understanding the evolutionary implications of natural hybridization in woody plants.

Completing the hybridization triangle: the inheritance of genetic incompatibilities during homoploid hybrid speciation in ragworts (Senecio)

A new homoploid hybrid lineage needs to establish a degree of reproductive isolation from its parent species if it is to persist as an independent entity, but the role hybridization plays in this process is known in only a handful of cases. The homoploid hybrid ragwort species, Senecio squalidus (Oxford ragwort), originated following the introduction of hybrid plants to the UK approximately 320 years ago. The source of the hybrid plants was from a naturally occurring hybrid zone between S. aethnensis and S. chrysanthemifolius on Mount Etna, Sicily. Previous studies of the parent species found evidence for multiple incompatibility loci causing transmission ratio distortion of genetic markers in their hybrid progeny. This study closes the hybridization triangle by reporting a genetic mapping analysis of the remaining two paired cross combinations between S. squalidus and its parents. Genetic maps produced from F₂ mapping families were generally collinear but with half of the linkage groups showing evidence of genomic reorganization between genetic maps. The new maps produced from crosses between S. squalidus and each parent showed multiple incompatibility loci distributed across the genome, some of which co-locate with previously reported incompatibility loci between the parents. These findings suggest that this young homoploid hybrid species has inherited a unique combination of genomic rearrangements and incompatibilities from its parents that contribute to its reproductive isolation.

Using Genomic Data to Infer Historic Population Dynamics of Nonmodel Organisms

Genome sequence data are now being routinely obtained from many nonmodel organisms. These data contain a wealth of information about the demographic history of the populations from which they originate. Many sophisticated statistical inference procedures have been developed to infer the demographic history of populations from this type of genomic data. In this review, we discuss the different statistical methods available for inference of demography, providing an overview of the underlying theory and logic behind each approach. We also discuss the types of data required and the pros and cons of each method. We then discuss how these methods have been applied to a variety of nonmodel organisms. We conclude by presenting some recommendations for researchers looking to use genomic data to infer demographic history.

Pleistocene glacial cycles drive isolation, gene flow and speciation in the high-elevation Andes

Mountain ranges are amongst the most species-rich habitats, with many large and rapid evolutionary radiations. The tempo and mode of diversification in these systems are key unanswered questions in evolutionary biology. Here we study the Andean Lupinus radiation to understand the processes driving very rapid diversification in montane systems. We use genomic and transcriptomic data of multiple species and populations, and apply phylogenomic and demographic analyses to test whether diversification proceeded without interspecific gene flow - as expected if Andean orogeny and geographic isolation were the main drivers of diversification - or if diversification was accompanied by gene flow, in which case other processes were probably involved. We uncover several episodes of gene flow between species, including very recent events likely to have been prompted by changes in habitat connectivity during Pleistocene glacial cycles. Furthermore, we find that gene flow between species was heterogeneously distributed across the genome. We argue that exceptionally fast diversification of Andean Lupinus was partly a result of Late Pleistocene glacial cycles, with associated cycles of expansion and contraction driving geographic isolation or secondary contact of species. Furthermore, heterogeneous gene flow across the genome suggests a role for selection and ecological speciation in rapid diversification in this system.

Genomic data reveal cryptic lineage diversification and introgression in Californian golden cup oaks (section Protobalanus)

Here we study hybridization, introgression and lineage diversification in the widely distributed canyon live oak (Quercus chrysolepis) and the relict island oak (Q. tomentella), two Californian golden cup oaks with an intriguing biogeographical history. We employed restriction-site-associated DNA sequencing and integrated phylogenomic and population genomic analyses to study hybridization and reconstruct the evolutionary past of these taxa. Our analyses revealed the presence of two cryptic lineages within Q. chrysolepis. One of these lineages shares its most recent common ancestor with Q. tomentella, supporting the paraphyly of Q. chrysolepis. The split of these lineages was estimated to take place during the late Pliocene or the early Pleistocene, a time corresponding well with the common presence of Q. tomentella in the fossil records of continental California. Analyses also revealed historical hybridization among lineages, high introgression from Q. tomentella into Q. chrysolepis in their current area of sympatry, and widespread admixture between the two lineages of Q. chrysolepis in contact zones. Our results support that the two lineages of Q. chrysolepis behave as a single functional species phenotypically and ecologically well differentiated from Q. tomentella, a situation that can be only accommodated considering hybridization and speciation as a continuum with diffuse limits.

Convergent and divergent evolution in carnivorous pitcher plant traps

Contents Summary 1035 I. Introduction 1035 II. Evolution of the pitcher 1036 III. Convergent evolution 1036 IV. Divergent evolution 1038 V. Adaptive radiation and speciation 1040 VI. Conclusions and perspectives 1040 Acknowledgements 1040 References 1040 SUMMARY: The pitcher trap is a striking example of convergent evolution across unrelated carnivorous plant lineages. Convergent traits that have evolved across pitcher plant lineages are essential for trap function, suggesting that key selective pressures are in action. Recent studies have also revealed patterns of divergent evolution in functional pitcher morphology within genera. Adaptations to differences in local prey assemblages may drive such divergence and, ultimately, speciation. Here, we review recent research on convergent and divergent evolution in pitcher plant traps, with a focus on the genus Nepenthes, which we propose as a new model for research into adaptive radiation and speciation.

Genome-wide analyses of the Bemisia tabaci species complex reveal contrasting patterns of admixture and complex demographic histories

Once considered a single species, the whitefly, Bemisia tabaci, is a complex of numerous morphologically indistinguishable species. Within the last three decades, two of its members (MED and MEAM1) have become some of the world's most damaging agricultural pests invading countries across Europe, Africa, Asia and the Americas and affecting a vast range of agriculturally important food and fiber crops through both feeding-related damage and the transmission of numerous plant viruses. For some time now, researchers have relied on a single mitochondrial gene and/or a handful of nuclear markers to study this species complex. Here, we move beyond this by using 38,041 genome-wide Single Nucleotide Polymorphisms, and show that the two invasive members of the complex are closely related species with signatures of introgression with a third species (IO). Gene flow patterns were traced between contemporary invasive populations within MED and MEAM1 species and these were best explained by recent international trade. These findings have profound implications for delineating the B. tabaci species status and will impact quarantine measures and future management strategies of this global pest.

Hybrids and horizontal transfer: introgression allows adaptive allele discovery

Evolution has devised countless remarkable solutions to diverse challenges. Understanding the mechanistic basis of these solutions provides insights into how biological systems can be subtly tweaked without maladaptive consequences. The knowledge gained from illuminating these mechanisms is equally important to our understanding of fundamental evolutionary mechanisms as it is to our hopes of developing truly rational plant breeding and synthetic biology. In particular, modern population genomic approaches are proving very powerful in the detection of candidate alleles for mediating consequential adaptations that can be tested functionally. Especially striking are signals gained from contexts involving genetic transfers between populations, closely related species, or indeed between kingdoms. Here we discuss two major classes of these scenarios, adaptive introgression and horizontal gene flow, illustrating discoveries made across kingdoms.

Adaptive introgression from distant Caribbean islands contributed to the diversification of a microendemic adaptive radiation of trophic specialist pupfishes

Rapid diversification often involves complex histories of gene flow that leave variable and conflicting signatures of evolutionary relatedness across the genome. Identifying the extent and source of variation in these evolutionary relationships can provide insight into the evolutionary mechanisms involved in rapid radiations. Here we compare the discordant evolutionary relationships associated with species phenotypes across 42 whole genomes from a sympatric adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas and several outgroup pupfish species in order to understand the rarity of these trophic specialists within the larger radiation of Cyprinodon. 82% of the genome depicts close evolutionary relationships among the San Salvador Island species reflecting their geographic proximity, but the vast majority of variants fixed between specialist species lie in regions with discordant topologies. Top candidate adaptive introgression regions include signatures of selective sweeps and adaptive introgression of genetic variation from a single population in the northwestern Bahamas into each of the specialist species. Hard selective sweeps of genetic variation on San Salvador Island contributed 5 times more to speciation of trophic specialists than adaptive introgression of Caribbean genetic variation; however, four of the 11 introgressed regions came from a single distant island and were associated with the primary axis of oral jaw divergence within the radiation. For example, standing variation in a proto-oncogene (ski) known to have effects on jaw size introgressed into one San Salvador Island specialist from an island 300 km away approximately 10 kya. The complex emerging picture of the origins of adaptive radiation on San Salvador Island indicates that multiple sources of genetic variation contributed to the adaptive phenotypes of novel trophic specialists on the island. Our findings suggest that a suite of factors, including rare adaptive introgression, may be necessary for adaptive radiation in addition to ecological opportunity.

Population-genomic inference of the strength and timing of selection against gene flow

The interplay of divergent selection and gene flow is key to understanding how populations adapt to local environments and how new species form. Here, we use DNA polymorphism data and genome-wide variation in recombination rate to jointly infer the strength and timing of selection, as well as the baseline level of gene flow under various demographic scenarios. We model how divergent selection leads to a genome-wide negative correlation between recombination rate and genetic differentiation among populations. Our theory shows that the selection density (i.e., the selection coefficient per base pair) is a key parameter underlying this relationship. We then develop a procedure for parameter estimation that accounts for the confounding effect of background selection. Applying this method to two datasets from Mimulus guttatus, we infer a strong signal of adaptive divergence in the face of gene flow between populations growing on and off phytotoxic serpentine soils. However, the genome-wide intensity of this selection is not exceptional compared with what M. guttatus populations may typically experience when adapting to local conditions. We also find that selection against genome-wide introgression from the selfing sister species M. nasutus has acted to maintain a barrier between these two species over at least the last 250 ky. Our study provides a theoretical framework for linking genome-wide patterns of divergence and recombination with the underlying evolutionary mechanisms that drive this differentiation.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

Genetic consequences of Quaternary climatic oscillations in the Himalayas: Primula tibetica as a case study based on restriction site-associated DNA sequencing

The effects of Quaternary climatic oscillations on the demography of organisms vary across regions and continents. In taxa distributed in Europe and North America, several paradigms regarding the distribution of refugia have been identified. By contrast, less is known about the processes that shaped the species' spatial genetic structure in areas such as the Himalayas, which is considered a biodiversity hotspot. Here, we investigated the phylogeographic structure and population dynamics of Primula tibetica by combining genomic phylogeography and species distribution models (SDMs). Genomic data were obtained for 293 samples of P. tibetica using restriction site-associated DNA sequencing (RADseq). Ensemble SDMs were carried out to predict potential present and past distribution ranges. Four distinct lineages were identified. Approximate Bayesian computation analyses showed that each of them have experienced both expansions and bottlenecks since their divergence, which occurred during or across the Quaternary glacial cycles. The two lineages at both edges of the distribution were found to be more vulnerable and responded in different ways to past climatic changes. These results illustrate how past climatic changes affected the demographic history of Himalayan organisms. Our findings highlight the significance of combining genomic approaches with environmental data when evaluating the effects of past climatic changes.

Genomic architecture of phenotypic divergence between two hybridizing plant species along an elevational gradient

Knowledge of the genetic basis of phenotypic divergence between species and how such divergence is caused and maintained is crucial to an understanding of speciation and the generation of biodiversity. The hybrid zone between Senecio aethnensis and S. chrysanthemifolius on Mount Etna, Sicily, provides a well-studied example of species divergence in response to conditions at different elevations, despite hybridization and gene flow. Here, we investigate the genetic architecture of divergence between these two species using a combination of quantitative trait locus (QTL) mapping and genetic differentiation measures based on genetic marker analysis. A QTL architecture characterized by physical QTL clustering, epistatic interactions between QTLs, and pleiotropy was identified, and is consistent with the presence of divergent QTL complexes resistant to gene flow. A role for divergent selection between species was indicated by significant negative associations between levels of interspecific genetic differentiation at mapped marker gene loci and map distance from QTLs and hybrid incompatibility loci. Within-species selection contributing to interspecific differentiation was evidenced by negative associations between interspecific genetic differentiation and genetic diversity within species. These results show that the two Senecio species, while subject to gene flow, maintain divergent genomic regions consistent with local selection within species and selection against hybrids between species which, in turn, contribute to the maintenance of their distinct phenotypic differences.