DIVERGENCE IS FOCUSED ON FEW GENOMIC REGIONS EARLY IN SPECIATION: INCIPIENT SPECIATION OF SUNFLOWER ECOTYPES

Computing the effects of temperature and osmotic stress on the seed germination of Helianthus annuus L. by using a mathematical model

An extremely important oil crop in the world, Helianthus annuus L. is one of the world's most significant members of the Asteraceae family. The rate and extent of seed germination and agronomic features are consistently affecting  by temperature (T) and changes in water potential (ψ). A broad hydrothermal time model with T and ψ components could explain sunflower responses over suboptimal T and ψ. A lab experiment was performed using the HTT model to discover both T and ψ and their interactive effects on sunflower germination and also to figure  out the cardinal Ts values. The sunflower seeds were germinated at temperatures (15 °C, 20 °C, 25 °C and 30 °C); each Ts had five constant ψs of 0, 0.3, 0.6, 0.9, and 1.2 MPa via PEG 6000 as osmotic stress inducer. The results revealed that highest germination index was found in seed grown at 20 °C in distilled water (0 MPa) and the lowest at 30 °C with osmotic stress of (- 1.2 MPa). The highest value of germination rate index was found in seed grown at 20 °C in distilled water (0 MPa) and the lowest at 15 °C with an osmotic stress of (- 1.2 MPa). In conclusion, water potential, temperature, and their interactions have a considerable impact on seed germination rate, and other metrics (GI, SVI-I, GRI, GE, SVI-II, and MGT). Seeds sown  at 20 °C with zero water potential showed high germination metrics such as GE, GP, GRI, and T50%. The maximum value to TTsub noted at 30 °C in - 0.9 MPa osmotic stress and the minimum value was calculated at 15 °C in - 1.2 MPa osmotic stress. The result of TTsupra recorded highest at 15 °C in  controlled group (0 MPa). Moreover, θH was  highest at 30 °C in controlled condition (0 MPa) and minimum value was observed at  20 °C under - 1.2 MPa osmotic stress. The value of θHTT were  maximum at  30 °C in controlled group (0 MPa) and minimum value was  recorded at 15 °C under - 1.2 MPa osmotic potential. The base, optimum and ceiling temperatures for sunflower germination metrics in this experiment were noted  6.8, 20 and 30 °C respectively.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Adaptation to distinct habitats is maintained by contrasting selection at different life stages in sunflower ecotypes

Conspecific populations living in adjacent but contrasting microenvironments represent excellent systems for studying natural selection. These systems are valuable because gene flow is expected to force genetic homogeneity except at loci experiencing divergent selection. A history of reciprocal transplant and common garden studies in such systems, and a growing number of genomic studies, have contributed to understanding how selection operates in natural populations. While selection can vary across different fitness components and life stages, few studies have investigated how this ultimately affects allele frequencies and the maintenance of divergence between populations. Here, we study two sunflower ecotypes in distinct, adjacent habitats by combining demographic models with genome-wide sequence data to estimate fitness and allele frequency change at multiple life stages. This framework allows us to estimate that only local ecotypes are likely to experience positive population growth (λ > 1) and that the maintenance of divergent adaptation appears to be mediated via habitat- and life stage-specific selection. We identify genetic variation, significantly driven by loci in chromosomal inversions, associated with different life history strategies in neighbouring ecotypes that optimize different fitness components and may contribute to the maintenance of distinct ecotypes.

Peripatric speciation within Torreya fargesii (Taxaceae) in the Hengduan Mountains inferred from multi-loci phylogeography

BACKGROUND

The Hengduan Mountains (HDM) are one of the major global biodiversity hotspots in the world. Several evolutionary scenarios, especially in-situ diversification, have been proposed to account for the high species richness of temperate plants. However, peripatric speciation, an important mode of allopatric speciation, has seldom been reported in this region.

RESULTS

Here, two chloroplast DNA regions and 14 nuclear loci were sequenced for 112 individuals from 10 populations of Torreya fargesii var. fargesii and 63 individuals from 6 populations of T. fargesii var. yunnanensis. Population genetic analyses revealed that the two varieties are well differentiated genetically (FST, 0.5765) and have uneven genetic diversity (π, 0.00221 vs. 0.00073 on an average of nuclear loci). The gene genealogical relationship showed that T. fargesii var. yunnanensis is inferred as derived from T. fargesii var. fargesii, which was further supported by the coalescent simulations (DIYABC, fastsimcoal2 and IMa2). By the coalescent simulations, the divergence time (~ 2.50-3.65 Ma) and the weak gene flow between the two varieties were detected. The gene flow was asymmetrical and only occurred in later stages of divergence, which is caused by second contact due to the population expansion (~ 0.61 Ma) in T. fargesii var. fargesii. In addition, niche modeling indicated that the two varieties are differentiated geographically and ecologically and have unbalanced distribution range.

CONCLUSIONS

Overall, T. fargesii var. fargesii is always parapatric with respect to T. fargesii var. yunnanensis, and the latter derived from the former in peripatry of the HDM following a colonization from central China during the late Pliocene. Our findings demonstrate that peripatric speciation following dispersal events may be an important evolutionary scenario for the formation of biodiversity hotspot of the HDM.

Where to draw the boundaries? Using landscape genomics to disentangle the scribbly gum species complex

PREMISE

Species delimitation is an integral part of evolution and ecology and is vital in conservation science. However, in some groups, species delimitation is difficult, especially where ancestral relationships inferred from morphological or genetic characters are discordant, possibly due to a complicated demographic history (e.g., recent divergences between lineages). Modern genetic techniques can take into account complex histories to distinguish species at a reasonable cost and are increasingly used in numerous applications. We focus on the scribbly gums, a group of up to five closely related and morphologically similar "species" within the eucalypts.

METHODS

Multiple populations of each recognized scribbly gum species were sampled over a wide region across climates, and genomewide scans were used to resolve species boundaries.

RESULTS

None of the taxa were completely divergent, and there were two genetically distinct entities: the inland distributed Eucalyptus rossii and a coastal conglomerate consisting of four species forming three discernible, but highly admixed groups. Divergence among taxa was likely driven by temporal vicariant processes resulting in partial separation across biogeographic barriers. High interspecific gene flow indicated separated taxa reconnected at different points in time, blurring species boundaries.

CONCLUSIONS

Our results highlight the need for genetic screening when dealing with closely related taxonomic entities, particularly those with modest morphological differences. We show that high-throughput sequencing can be effective at identifying species groupings and processes driving divergence, even in the most taxonomically complex groups, and be used as a standard practice for disentangling species complexes.

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.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Genomic divergence between two sister Ostrya species through linked selection and recombination

Studying the evolution of genomic divergence between lineages is a topical issue in evolutionary biology. However, the evolutionary forces that shape the heterogeneous divergence of the genomic landscape are still poorly understood. Here, two wind-pollinated sister-species (Ostrya japonica and O. chinensis) are used to explore what these potential forces might be. A total of 40 individuals from 16 populations across their main distribution areas in China were sampled for genome-wide resequencing. Population demography analyses revealed that these two sister-species diverged at 3.06-4.43 Mya. Both population contraction and increased gene flow were detected during glacial periods, suggesting secondary contact at those times. All three parameters (D XY, π, and ρ) decreased in those regions showing high levels of differentiation (F ST). These findings indicate that linked selection and recombination played a key role in the genomic heterogeneous differentiation between the two Ostrya species. Genotype-environment association analyses showed that precipitation was the most important ecological factor for speciation. Such environmentally related genes and positive selection genes may have contributed to local adaptation and the maintenance of species boundaries.

Eco-Geography and Phenology Are the Major Drivers of Reproductive Isolation in the Royal Irises, a Species Complex in the Course of Speciation

The continuous nature of speciation implies that different species are found at different stages of divergence, from no- to complete reproductive isolation. This process and its underlying mechanisms are best viewed in incipient species. Moreover, the species complex can offer unique insight into how reproductive isolation (RI) has evolved. The royal irises (Iris section Oncocyclus) are a young group of species in the course of speciation, providing an ideal system for speciation study. We quantified pre- and post-zygotic reproductive barriers between the eight Israeli species of this complex and estimated the total RI among them. We tested for both pre-pollination and post-pollination reproductive barriers. Pre-pollination barriers, i.e., eco-geographic divergence and phenological differentiation were the major contributors to RI among the Iris species. On the other hand, post-pollination barriers, namely pollen-stigma interactions, fruit set, and seed viability had negligible contributions to total RI. The strength of RI was not uniform across the species complex, suggesting that species may have diverged at different rates. Overall, this study in a young, recently diverged group of species provides insight into the first steps of speciation, suggesting a crucial role of the pre-zygotic barriers.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae)

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats, and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographic barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site associated DNA sequencing, (ii) genetic structure analyses, (iii) demographic modeling, (iv) morphometrics, and (v) climatic niche modeling and niche overlap analysis. Results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread I. linifolia during Late Pliocene. Meanwhile, geographic isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographic barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

Independent Evolutionary Lineages in a Globular Cactus Species Complex Reveals Hidden Diversity in a Central Chile Biodiversity Hotspot

Unraveling the processes involved in the origin of a substantial fraction of biodiversity can be a particularly difficult task in groups of similar, and often convergent, morphologies. The genus Eriosyce (Cactaceae) might present a greater specific diversity since much of its species richness might be hidden in morphological species complexes. The aim of this study was to investigate species delimitation using the molecular data of the globose cacti “E. curvispina”, which harbor several populations of unclear evolutionary relationships. We ran phylogenetic inferences on 87 taxa of Eriosyce, including nine E. curvispina populations, and by analyzing three plastid noncoding introns, one plastid and one nuclear gene. Additionally, we developed 12 new pairs of nuclear microsatellites to evaluate the population-level genetic structure. We identified four groups that originated in independent cladogenetic events occurring at different temporal depths; these groups presented high genetic diversity, and their populations were genetically structured. These results suggest a complex evolutionary history in the origin of globular cacti, with independent speciation events occurring at different time spans. This cryptic richness is underestimated in the Mediterranean flora of central Chile, and thus unique evolutionary diversity could be overlooked in conservation and management actions.

Modulating response of sunflower (Hellianthus annuus) to induced salinity stress through application of engineered urea functionalized hydroxyapatite nanoparticles

Agro-nanotechnology aims to improve the quality and quantity of plants and plant products while preserving environmental health. Contemporary anecdotal studies that provide representation of the use of nanostructures as fertilizers, pesticides, and plant growth regulators have highlighted the need to determine the effect of such modified nanofertilizers on transforming plant yield under abiotic stress. Present study was performed to modulate the physiological response of Hellianthus annuus through the application of Urea capped hydroxyapatite nanoparticles (Urea-HANPs) in stressed environment. Hydroxyapatite nanoparticles were synthesized via co-precipitation method, functionalized with urea and characterized through a series of contemporary techniques of transmission electron microscopy (TEM), energy dispersive X-ray spectroscopy (EDX), X-ray diffraction (XRD), and Fourier transform infrared (FTIR) spectroscopy. We observed that Urea-HANPs significantly (p < .05) ameliorated resistivity in plant to osmotic stress by enhancing agronomic and physiobiochemical attributes. Elevated chlorophyll contents were reported from tested leaves treated with Urea-HANPs in T6 (0.05 M NaCl + 10 μg/ml Urea-HANP) under induced salinity stress. Data revealed significant decrease in osmolytes at T3 (0.1 M NaCl), and T4 (0.2 M NaCl) which was significantly ameliorated in T9 (0.1 M NaCl + 10 μg/ml Urea-HANPs) and T12 (0.2 M NaCl + 10 μg/ml Urea-HANPs). The antioxidant system was boosted up by the application of Urea-HANPs preventing the plant from oxidative stress by scavenging reactive oxygen species (ROS). It has been concluded that alleviation in impact of osmotic stress on plant through the use of Urea-HANPs was coupled with elevation in photosynthetic performance, stimulation of osmolytes and boosting antioxidant system of plants.

Standing variation rather than recent adaptive introgression probably underlies differentiation of the texanus subspecies of Helianthus annuus

The origins of geographic races in wide-ranging species are poorly understood. In Texas, the texanus subspecies of Helianthus annuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies.

Microsatellites as Agents of Adaptive Change: An RNA-Seq-Based Comparative Study of Transcriptomes from Five Helianthus Species

Mutations that provide environment-dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates achieved through symmetric or asymmetric mutation processes and are abundant in both coding and non-coding regions in positions that could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from five closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5′UTRs and coding regions compared to 3′UTRs and non-coding transcripts among Helianthus spp. Furthermore, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

Population genomic structure of Eurasian and African foxtail millet landrace accessions inferred from genotyping-by-sequencing

Foxtail millet [Setaria italica (L.) P. Beauv.] is the second most important millet species globally and is adapted to cultivation in diverse environments. Like its wild progenitor, green foxtail [S. viridis (L.) P. Beauv.], it is a model species for C4 photosynthetic pathways and stress tolerance genes in related bioenergy crops. We addressed questions regarding the evolution and spread of foxtail millet through a population genomic study of landraces from across its cultivated range in Europe, Asia, and Africa. We sought to determine population genomic structure and the relationship of domesticated lineages relative to green foxtail. Further, we aimed to identify genes involved in environmental stress tolerance that have undergone differential selection between geographical and genetic groups. Foxtail millet landrace accessions (n = 328) and green foxtail accessions (n = 12) were sequenced by genotyping-by-sequencing (GBS). After filtering, 5,677 single nucleotide polymorphisms (SNPs) were retained for the combined foxtail millet-green foxtail dataset and 5,020 for the foxtail millet dataset. We extended geographic coverage of green foxtail by including previously published GBS sequence tags, yielding a 4,515-SNP dataset for phylogenetic reconstruction. All foxtail millet samples were monophyletic relative to green foxtail, suggesting a single origin of foxtail millet, although no group of foxtail millet was clearly the most ancestral. Four genetic clusters were found within foxtail millet, each with a distinctive geographical distribution. These results, together with archaeobotanical evidence, suggest plausible routes of spread of foxtail millet. Selection scans identified nine candidate genes potentially involved in environmental adaptations, particularly to novel climates encountered, as domesticated foxtail millet spread to new altitudes and latitudes.

Hybridization and introgression of native and foreign Sorbus tree species in unique environments of protected mountainous areas

Hybridization and introgression are important processes influencing the genetic diversity and evolution of species. These processes are of particular importance in protected areas, where they can lead to the formation of hybrids between native and foreign species and may ultimately result in the loss of parental species from their natural range. Despite their importance, the contribution of hybridization and introgression to genetic diversity in Sorbus genus remains not fully recognized. We analysed the genetic and morphological variability of several Sorbus species including native (Sorbus aria), foreign (S. intermedia) and potentially hybrid (S. carpatica) individuals from the Polish Carpathian range. Patterns of variation at 13 nuclear microsatellite loci show hybridization between the tested species and confirm the existence of the hybrid form S. carpatica. Biometric analysis on leaves, based of 10 metric features and three parameters, identified several characters for preliminary taxonomic classification; however, none of them could be used as a fully diagnostic marker for faultless annotation of S. intermedia and S. carpatica. The genetic structure analysis indicated complex patterns of population differentiation and its diverse origin. The results allow assessment of genetic variation and identification of parental species participating in hybridization. This knowledge will advance the management of genetic diversity and development of conservation strategies for efficient maintenance of the unique protected ecosystem.

Detecting Selection from Linked Sites Using an F-Model

Allele frequencies vary across populations and loci, even in the presence of migration. While most differences may be due to genetic drift, divergent selection will further increase differentiation at some loci. Identifying those is key in studying local adaptation, but remains statistically challenging. A particularly elegant way to describe allele frequency differences among populations connected by migration is the F-model, which measures differences in allele frequencies by population specific FST coefficients. This model readily accounts for multiple evolutionary forces by partitioning FST coefficients into locus- and population-specific components reflecting selection and drift, respectively. Here we present an extension of this model to linked loci by means of a hidden Markov model (HMM), which characterizes the effect of selection on linked markers through correlations in the locus specific component along the genome. Using extensive simulations, we show that the statistical power of our method is up to twofold higher than that of previous implementations that assume sites to be independent. We finally evidence selection in the human genome by applying our method to data from the Human Genome Diversity Project (HGDP).

The Genomic Landscape of Divergence Across the Speciation Continuum in Island-Colonising Silvereyes (Zosterops lateralis)

Inferring the evolutionary dynamics at play during the process of speciation by analyzing the genomic landscape of divergence is a major pursuit in population genomics. However, empirical assessments of genomic landscapes under varying evolutionary scenarios that are known a priori are few, thereby limiting our ability to achieve this goal. Here we combine RAD-sequencing and individual-based simulations to evaluate the genomic landscape of divergence in the silvereye (Zosterops lateralis). Using pairwise comparisons that differ in divergence timeframe and the presence or absence of gene flow, we document how genomic patterns accumulate along the speciation continuum. In contrast to previous predictions, our results provide limited support for the idea that divergence accumulates around loci under divergent selection or that genomic islands widen with time. While a small number of genomic islands were found in populations diverging with and without gene flow, in few cases were SNPs putatively under selection tightly associated with genomic islands. The transition from localized to genome-wide levels of divergence was captured using individual-based simulations that considered only neutral processes. Our results challenge the ubiquity of existing verbal models that explain the accumulation of genomic differences across the speciation continuum and instead support the idea that divergence both within and outside of genomic islands is important during the speciation process.

Adaptive zones shape the magnitude of premating reproductive isolation in Timema stick insects

Simpson's fossil-record inspired model of 'adaptive zones' proposes that evolution is dominated by small fluctuations within adaptive zones, occasionally punctuated by larger shifts between zones. This model can help explain why the process of population divergence often results in weak or moderate reproductive isolation (RI), rather than strong RI and distinct species. Applied to the speciation process, the adaptive zones hypothesis makes two inter-related predictions: (i) large shifts between zones are relatively rare, (ii) when large shifts do occur they generate stronger RI than shifts within zones. Here, we use ecological, phylogenetic and behavioural data to test these predictions in Timema stick insects. We show that host use in Timema is dominated by moderate shifts within the systematic divisions of flowering plants and conifers, with only a few extreme shifts between these divisions. However, when extreme shifts occur, they generate greater RI than do more moderate shifts. Our results support the adaptive zones model, and suggest that the net contribution of ecological shifts to diversification is dependent on both their magnitude and frequency. We discuss the generality of our findings in the light of emerging evidence from diverse taxa that the evolution of RI is not always the only factor determining the origin of species diversity. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Massive haplotypes underlie ecotypic differentiation in sunflowers

Species often include multiple ecotypes that are adapted to different environments¹. However, it is unclear how ecotypes arise and how their distinctive combinations of adaptive alleles are maintained despite hybridization with non-adapted populations^2-4. Here, by resequencing 1,506 wild sunflowers from 3 species (Helianthus annuus, Helianthus petiolaris and Helianthus argophyllus), we identify 37 large (1-100 Mbp in size), non-recombining haplotype blocks that are associated with numerous ecologically relevant traits, as well as soil and climate characteristics. Limited recombination in these haplotype blocks keeps adaptive alleles together, and these regions differentiate sunflower ecotypes. For example, haplotype blocks control a 77-day difference in flowering between ecotypes of the silverleaf sunflower H. argophyllus (probably through deletion of a homologue of FLOWERING LOCUS T (FT)), and are associated with seed size, flowering time and soil fertility in dune-adapted sunflowers. These haplotypes are highly divergent, frequently associated with structural variants and often appear to represent introgressions from other-possibly now-extinct-congeners. These results highlight a pervasive role of structural variation in ecotypic adaptation.

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time-consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome-environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described "islands of differentiation," and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.

Population Genomics in Rhamdia quelen (Heptapteridae, Siluriformes) Reveals Deep Divergence and Adaptation in the Neotropical Region

Rhamdia quelen, a Neotropical fish with hybridization between highly divergent mitochondrial DNA (mtDNA) lineages, represents an interesting evolutionary model. Previous studies suggested that there might be demographic differences between coastal lagoons and riverine environments, as well as divergent populations that could be reproductively isolated. Here, we investigated the genetic diversity pattern of this taxon in the Southern Neotropical Basin system that includes the La Plata Basin, Patos-Merin lagoon basin and the coastal lagoons draining to the SW Atlantic Ocean, through a population genomics approach using 2b-RAD-sequencing-derived single nucleotide polymorphisms (SNPs). The genomic scan identified selection footprints associated with divergence and suggested local adaptation environmental drivers. Two major genomic clusters latitudinally distributed in the Northern and Southern basins were identified, along with consistent signatures of divergent selection between them. Population structure based on the whole set of loci and on the presumptive neutral vs. adaptive loci showed deep genomic divergence between the two major clusters. Annotation of the most consistent SNPs under divergent selection revealed some interesting candidate genes for further functional studies. Moreover, signals of adaptation to a coastal lagoon environment mediated by purifying selection were found. These new insights provide a better understanding of the complex evolutionary history of R. quelen in the southernmost basin of the Neotropical region.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA-sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL-charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL-charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI-charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.

Intraspecific genetic divergence within Helianthus niveus and the status of two new morphotypes from Mexico

PREMISE

Collecting and characterizing the genetic diversity of wild relatives of crops can contribute importantly to sustainable crop production and food security. Wild sunflower, Helianthus niveus, occurs in arid regions in western North America and is partially cross-compatible with the cultivated sunflower (H. annuus). We assessed phylogenetic relationships and patterns of genetic divergence among three previously described subspecies (subsp. niveus, subsp. canescens, and subsp. tephrodes) as well as two new morphotypes of H. niveus recently discovered in extreme drought and dune habitats in Baja California, Mexico.

METHODS

We measured 50 plants growing in a common garden for 27 morphological traits and conducted principal component analysis to assess patterns of phenotypic variation. Genome size of each accession was determined using flow cytometry. Pollen viability of first generation hybrids between taxa was tested to infer the strength of intrinsic postzygotic reproductive barriers. Finally, genotyping-by-sequencing data were used to investigate the genetic structure and phylogenetic relationships among the previously described subspecies and new morphotypes.

RESULTS

The intraspecific genetic and phenotypic divergence of H. niveus populations closely tracks their geographical distribution. Subspecies niveus is phenotypically, genetically, and reproductively distinct from the other two subspecies and has a larger genome. Therefore, H. niveus as currently circumscribed should be considered to contain two distinct species, H. niveus and H. tephrodes. ABBA-BABA tests revealed substantial introgression between subsp. canescens and its sympatric congener H. petiolaris, which might contribute to their morphological similarities. The two new morphotypes collected in Mexico represent local ecotypes of subsp. niveus that occur in extreme drought and dune environments. Mantel tests showed a strong positive correlation between genetic and geographic distances.

CONCLUSIONS

We conclude that geographic isolation is primarily responsible for intraspecific genomic divergence within H. niveus, while patterns of phenotypic variation appear to have been shaped by ecological selection and interspecific introgression.

Divergence with gene flow is driven by local adaptation to temperature and soil phosphorus concentration in teosinte subspecies (Zea mays parviglumis and Zea mays mexicana)

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations.

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.

The Rapid Evolution of an Ohnolog Contributes to the Ecological Specialization of Incipient Yeast Species

Identifying the molecular changes that lead to ecological specialization during speciation is one of the major goals of molecular evolution. One question that remains to be thoroughly investigated is whether ecological specialization derives strictly from adaptive changes and their associated trade-offs, or from conditionally neutral mutations that accumulate under relaxed selection. We used whole-genome sequencing, genome annotation and computational analyses to identify genes that have rapidly diverged between two incipient species of Saccharomyces paradoxus that occupy different climatic regions along a south-west to north-east gradient. As candidate loci for ecological specialization, we identified genes that show signatures of adaptation and accelerated rates of amino acid substitutions, causing asymmetric evolution between lineages. This set of genes includes a glycyl-tRNA-synthetase, GRS2, which is known to be transcriptionally induced under heat stress in the model and sister species S. cerevisiae. Molecular modelling, expression analysis and fitness assays suggest that the accelerated evolution of this gene in the Northern lineage may be caused by relaxed selection. GRS2 arose during the whole-genome duplication (WGD) that occurred 100 million years ago in the yeast lineage. While its ohnolog GRS1 has been preserved in all post-WGD species, GRS2 has frequently been lost and is evolving rapidly, suggesting that the fate of this ohnolog is still to be resolved. Our results suggest that the asymmetric evolution of GRS2 between the two incipient S. paradoxus species contributes to their restricted climatic distributions and thus that ecological specialization derives at least partly from relaxed selection rather than a molecular trade-off resulting from adaptive evolution.

Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.

A neutral view of the evolving genomic architecture of speciation

Analyses of genomewide polymorphism data have begun to shed light on speciation and adaptation. Genome scans to identify regions of the genome that are unusually different between populations or species, possibly due to divergent natural or sexual selection, are widespread in speciation genomics. Theoretical and empirical work suggests that such outlier regions may grow faster than linearly during speciation with gene flow due to a rapid transition between low and high reproductive isolation. We investigate whether this pattern could be attributed to neutral processes by simulating genomes under neutral evolution with varying amounts and timing of gene flow. Under both neutral evolution and divergent selection, simulations with little or no gene flow, or with a long allopatric period after its cessation, resulted in faster than linear growth of the proportion of the genome lying in outlier regions. Without selection, higher recent gene flow erased differentiation; with divergent selection, these same scenarios produced nonlinear growth to a plateau. Our results suggest that, given a history of gene flow, the growth of the divergent genome is informative about selection during divergence, but that in many scenarios, this pattern does not easily distinguish neutral and non-neutral processes during speciation with gene flow.

RADseq provides evidence for parallel ecotypic divergence in the autotetraploid Cochlearia officinalis in Northern Norway

Speciation encompasses a continuum over time from freely interbreeding populations to reproductively isolated species. Along this process, ecotypes - the result of local adaptation - may be on the road to new species. We investigated whether three autotetraploid Cochlearia officinalis ecotypes, adapted to different habitats (beach, estuary, spring), are genetically differentiated and result from parallel ecotypic divergence in two distinct geographical regions. We obtained genetic data from thousands of single nucleotide polymorphisms (SNPs) from restriction-site associated DNA sequencing (RADseq) and from six microsatellite markers for 12 populations to assess genetic divergence at ecotypic, geographic and population level. The genetic patterns support differentiation among ecotypes as suggested by morphology and ecology. The data fit a scenario where the ancestral beach ecotype has recurrently and polytopically given rise to the estuary and spring ecotypes. Several ecologically-relevant loci with consistent non-random segregating patterns are identified across the recurrent origins, in particular around genes related to salt stress. Despite being ecologically distinct, the Cochlearia ecotypes still represent an early stage in the process of speciation, as reproductive isolation has not (yet) developed. A sequenced annotated genome is needed to specifically target candidate genes underlying local adaptation.

Pollutants and Insecticides Drive Local Adaptation in African Malaria Mosquitoes

The Anopheles gambiae complex contains a number of highly anthropophilic mosquito species that have acquired exceptional ability to thrive in complex human habitats. Thus, examining the evolutionary history of this Afrotropical mosquito may yield vital information on the selective processes that occurred during the adaptation to human-dominated environments. We performed reduced representation sequencing on 941 mosquitoes of the Anopheles gambiae complex collected across four ecogeographic zones in Cameroon. We find evidence for genetic and geographic subdivision within An. coluzzii and An. gambiae sensu stricto-the two most significant malaria vectors in the region. Importantly, in both species, rural and urban populations are genetically differentiated. Genome scans reveal pervasive signatures of selection centered on genes involved in xenobiotic resistance. Notably, a selective sweep containing detoxification enzymes is prominent in urban mosquitoes that exploit polluted breeding sites. Overall, our study suggests that recent anthropogenic environmental modifications and widespread use of insecticides are driving population differentiation and local adaptation in vectors with potentially significant consequences for malaria epidemiology.

Incipient ecological speciation between successional varieties of a dominant tree involves intrinsic postzygotic isolating barriers

Whereas disruptive selection imposed by heterogeneous environments can lead to the evolution of extrinsic isolating barriers between diverging populations, the evolution of intrinsic postzygotic barriers through divergent selection is less certain. Long-lived species such as trees may be especially slow to evolve intrinsic isolating barriers. We examined postpollination reproductive isolating barriers below the species boundary, in an ephemeral hybrid zone between two successional varieties of the landscape-dominant Hawaiian tree, Metrosideros polymorpha, on volcanically active Hawai'i Island. These archipelago-wide sympatric varieties show the weakest neutral genetic divergence of any taxon pair on Hawai'i Island but significant morphological and ecological differentiation consistent with adaptation to new and old lava flows. Cross-fertility between varieties was high and included heterosis of F1 hybrids at the seed germination stage, consistent with a substantial genetic load apparent within varieties through low self-fertility and a lack of self-pollen discrimination. However, a partial, but significant, barrier was observed in the form of reduced female and male fertility of hybrids, especially backcross hybrids, consistent with the accumulation of genetic incompatibilities between varieties. These results suggest that partial intrinsic postzygotic barriers can arise through disruptive selection acting on large, hybridizing populations of a long-lived species.

A combination of sexual and ecological divergence contributes to rearrangement spread during initial stages of speciation

Chromosomal rearrangements between sympatric species often contain multiple loci contributing to assortative mating, local adaptation and hybrid sterility. When and how these associations arise during the process of speciation remains a subject of debate. Here, we address the relative roles of local adaptation and assortative mating on the dynamics of rearrangement evolution by studying how a rearrangement covaries with sexual and ecological trait divergence within a species. Previously, a chromosomal rearrangement that suppresses recombination on the Z (sex) chromosome was identified in European corn borer moths (Ostrinia nubilalis). We further characterize this recombination suppressor and explore its association with variation in sex pheromone communication and seasonal ecological adaptation in pairs of populations that are divergent in one or both of these characteristics. Direct estimates of recombination suppression in pedigree mapping families indicated that more than 39% of the Z chromosome (encompassing up to ~10 megabases and ~300 genes) resides within a nonrecombining unit, including pheromone olfactory receptor genes and a major quantitative trait locus that contributes to ecotype differences (Pdd). Combining direct and indirect estimates of recombination suppression, we found that the rearrangement was occasionally present between sexually isolated strains (E vs. Z) and between divergent ecotypes (univoltine vs. bivoltine). However, it was only consistently present when populations differed in both sexual and ecological traits. Our results suggest that independent of the forces that drove the initial establishment of the rearrangement, a combination of sexual and ecological divergence is required for rearrangement spread during speciation.

Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST  > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80-402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole-genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination-reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non-Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour-joining trees based on whole-genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.

Evolution of heterogeneous genome differentiation across multiple contact zones in a crow species complex

Uncovering the genetic basis of species diversification is a central goal in evolutionary biology. Yet, the link between the accumulation of genomic changes during population divergence and the evolutionary forces promoting reproductive isolation is poorly understood. Here, we analysed 124 genomes of crow populations with various degrees of genome-wide differentiation, with parallelism of a sexually selected plumage phenotype, and ongoing hybridization. Overall, heterogeneity in genetic differentiation along the genome was best explained by linked selection exposed on a shared genome architecture. Superimposed on this common background, we identified genomic regions with signatures of selection specific to independent phenotypic contact zones. Candidate pigmentation genes with evidence for divergent selection were only partly shared, suggesting context-dependent selection on a multigenic trait architecture and parallelism by pathway rather than by repeated single-gene effects. This study provides insight into how various forms of selection shape genome-wide patterns of genomic differentiation as populations diverge.

Multiple reproductive barriers separate recently diverged sunflower ecotypes

Measuring reproductive barriers between groups of organisms is an effective way to determine the traits and mechanisms that impede gene flow. However, to understand the ecological and evolutionary factors that drive speciation, it is important to distinguish between the barriers that arise early in the speciation process and those that arise after speciation is largely complete. In this article, we comprehensively test for reproductive isolation between recently diverged (<10,000 years bp) dune and nondune ecotypes of the prairie sunflower, Helianthus petiolaris. We find reproductive barriers acting at multiple stages of hybridization, including premating, postmating-prezygotic, and postzygotic barriers, despite the recent divergence. Barriers include extrinsic selection against immigrants and hybrids, a shift in pollinator assemblage, and postpollination assortative mating. Together, these data suggest that multiple barriers can be important for reducing gene flow in the earliest stages of speciation.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

A genomic perspective on hybridization and speciation

Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.

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

Hybrid zones typically form as a result of species coming into secondary contact, but can also be established in situ as an ecotonal hybrid zone, a situation which has been reported far less frequently. An altitudinal hybrid zone on Mount Etna between two ragwort species (the low elevation Senecio chrysanthemifolius and high elevation S. aethnensis) could potentially represent either of these possibilities. However, a scenario of secondary contact vs. speciation with gene flow has not been explicitly tested. Here, we test these alternatives and demonstrate that the data do not support secondary contact. Furthermore, we report that the previous analyses of speciation history of these species were based on admixed populations, which has led to inflated estimates of ongoing, interspecific gene flow. Our new analyses, based on 'pure' S. aethnensis and S. chrysanthemifolius populations, reveal gene exchange of less than one effective migrant per generation, a level low enough to allow the species to accumulate neutral, genomewide differences. Overall, our results are consistent with a scenario of speciation with gene flow and a divergence time which coincides with the rise of Mt. Etna to altitudes above 2000 m (~150 KY). Further work to quantify the role of adaptation to contrasting environments of high and low altitudes will be needed to support the scenario of recent ecological speciation in this system.