Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers

Sympatric Speciation in the Genomic Era

Sympatric speciation has been of key interest to biologists investigating how natural and sexual selection drive speciation without the confounding variable of geographic isolation. The advent of the genomic era has provided a more nuanced and quantitative understanding of the different and often complex modes of speciation by which sympatric sister taxa arose, and a reassessment of some of the most compelling empirical case studies of sympatric speciation. However, I argue that genomic studies based on contemporary populations may never be able to provide unequivocal evidence of true primary sympatric speciation, and there is a need to incorporate palaeogenomic studies into this field. This inability to robustly distinguish cases of primary and secondary 'divergence with gene flow' may be inconsequential, as both are useful for understanding the role of large effect barrier loci in the progression from localised genic isolation to genome-wide reproductive isolation. I argue that they can be of equivalent interest due to shared underlying mechanisms driving divergence and potentially leaving similar patterns of coalescence.

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.

Genomics of adaptive divergence with chromosome-scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome-scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low-crossover centres relative to the high-crossover peripheries of chromosomes ("Chromosome Centre-Biased Differentiation", CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome-scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome-scale heterogeneity in crossover rate in evolutionary genomics.

Direct estimate of the spontaneous germ line mutation rate in African green monkeys

Here, I provide the first direct estimate of the spontaneous mutation rate in an Old World monkey, using a seven individual, three-generation pedigree of African green monkeys. Eight de novo mutations were identified within ∼1.5 Gbp of accessible genome, corresponding to an estimated point mutation rate of 0.94 × 10^-8 per site per generation, suggesting an effective population size of ∼12000 for the species. This estimation represents a significant improvement in our knowledge of the population genetics of the African green monkey, one of the most important nonhuman primate models in biomedical research. Furthermore, by comparing mutation rates in Old World monkeys with the only other direct estimates in primates to date-humans and chimpanzees-it is possible to uniquely address how mutation rates have evolved over longer time scales. While the estimated spontaneous mutation rate for African green monkeys is slightly lower than the rate of 1.2 × 10^-8 per base pair per generation reported in chimpanzees, it is similar to the lower range of rates of 0.96 × 10^-8 -1.28 × 10^-8 per base pair per generation recently estimated from whole genome pedigrees in humans. This result suggests a long-term constraint on mutation rate that is quite different from similar evidence pertaining to recombination rate evolution in primates.

Genomic differentiation and patterns of gene flow between two long-tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions ("genomic differentiation islands") spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene-flow patterns for autosomes using RAD-seq data between two closely related species of long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.

Mosaic genome evolution in a recent and rapid avian radiation

Recent genomic analyses of evolutionary radiations suggest that ancestral or standing genetic variation may facilitate rapid diversification, particularly in cases involving convergence in ecological traits. Likewise, lateral transfer of alleles via hybridization may also facilitate adaptive convergence, but little is known about the role of ancestral variation in examples of explosive diversification that primarily involve the evolution of species recognition traits. Here, we show that genomic regions distinguishing sympatric species in an extraordinary radiation of small finches called munias (genus Lonchura) have phylogenetic histories that are discordant with each other, with the overall pattern of autosomal differentiation among species, and with sex-linked and mitochondrial components of the genome. Genome-wide data for 11 species sampled in Australia and Papua New Guinea indicate substantial autosomal introgression between sympatric species, but also identify a limited number of divergent autosomal regions, several of which overlap known colour genes (ASIP, EDN3, IGSF11, KITLG, MC1R and SOX10). Phylogenetic analysis of these outlier regions shows that different munia species have acquired unique combinations of alleles across a relatively small set of phenotypically relevant genes. Our results demonstrate that the recombination of ancestral genetic variation across multiple loci may be an important mechanism for generating phenotypic novelty and diversity.

Genomic insights into adaptive divergence and speciation among malaria vectors of the Anopheles nili group

Ongoing speciation in the most important African malaria vectors gives rise to cryptic populations, which differ remarkably in their behavior, ecology, and capacity to vector malaria parasites. Understanding the population structure and the drivers of genetic differentiation among mosquitoes is crucial for effective disease control because heterogeneity within vector species contributes to variability in malaria cases and allow fractions of populations to escape control efforts. To examine population structure and the potential impacts of recent large-scale control interventions, we have investigated the genomic patterns of differentiation in mosquitoes belonging to the Anopheles nili group-a large taxonomic group that diverged ~3 Myr ago. Using 4,343 single nucleotide polymorphisms (SNPs), we detected strong population structure characterized by high-FST values between multiple divergent populations adapted to different habitats within the Central African rainforest. Delineating the cryptic species within the Anopheles nili group is challenging due to incongruence between morphology, ribosomal DNA, and SNP markers consistent with incomplete lineage sorting and/or interspecific gene flow. A very high proportion of loci are fixed (FST = 1) within the genome of putative species, which suggests that ecological and/or reproductive barriers are maintained by strong selection on a substantial number of genes.

Assembly and RNA-free annotation of highly heterozygous genomes: The case of the thick-billed murre (Uria lomvia)

Thanks to a dramatic reduction in sequencing costs followed by a rapid development of bioinformatics tools, genome assembly and annotation have become accessible to many researchers in recent years. Among tetrapods, birds have genomes that display many features that facilitate their assembly and annotation, such as small genome size, low number of repeats and highly conserved genomic structure. However, we found that high genomic heterozygosity could have a great impact on the quality of the genome assembly of the thick-billed murre (Uria lomvia), an arctic colonial seabird. In this study, we tested the performance of three genome assemblers, ray/sscape, soapdenovo2 and platanus, in assembling the highly heterozygous genome of the thick-billed murre. Our results show that platanus, an assembler specifically designed for heterozygous genomes, outperforms the other two approaches and produces a highly contiguous (N50 = 15.8 Mb) and complete genome assembly (93% presence of genes from the Benchmarking Universal Single Copy Ortholog [BUSCO] gene set). Additionally, we annotated the thick-billed murre genome using a homology-based approach that takes advantage of the genomic resources available for birds and other taxa. Our study will be useful for those researchers who are approaching assembly and annotation of highly heterozygous genomes, or genomes of species of conservation concern, and/or who have limited financial resources.

Interpreting the genomic landscape of introgression

Introgression, the transfer of genetic material between species through hybridisation, occurs in many taxa and has important consequences. Genomic studies allow us to characterise the landscape of introgression across the genome, shedding light on both its adaptive benefits and the incompatibilities that help to maintain species barriers. Studies taking a genome-wide view suggest that adaptive introgression may be common, but that introgressed variation between many species is selected against throughout much of the genome. Confounding factors can complicate interpretations from these data, and computational simulations have proved vital to illustrate expected patterns under different scenarios. Future developments will move beyond correlative evidence to explicit models that account for how selection and genetic drift influence introgressed variation.

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.

Genomewide patterns of variation in genetic diversity are shared among populations, species and higher-order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST , PBS, d_xy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.

Genetic differences between willow warbler migratory phenotypes are few and cluster in large haplotype blocks

It is well established that differences in migratory behavior between populations of songbirds have a genetic basis but the actual genes underlying these traits remains largely unknown. In an attempt to identify such candidate genes we de novo assembled the genome of the willow warbler Phylloscopus trochilus, and used whole‐genome resequencing and a SNP array to associate genomic variation with migratory phenotypes across two migratory divides around the Baltic Sea that separate SW migrating P. t. trochilus wintering in western Africa and SSE migrating P. t. acredula wintering in eastern and southern Africa. We found that the genomes of the two migratory phenotypes lack clear differences except for three highly differentiated regions located on chromosomes 1, 3, and 5 (containing 146, 135, and 53 genes, respectively). Within each migratory phenotype we found virtually no differences in allele frequencies for thousands of SNPs, even when comparing geographically distant populations breeding in Scandinavia and Far East Russia (>6000 km). In each of the three differentiated regions, multidimensional scaling‐based clustering of SNP genotypes from more than 1100 individuals demonstrates the presence of distinct haplotype clusters that are associated with each migratory phenotype. In turn, this suggests that recombination is absent or rare between haplotypes, which could be explained by inversion polymorphisms. Whereas SNP alleles on chromosome 3 correlate with breeding altitude and latitude, the allele distribution within the regions on chromosomes 1 and 5 perfectly matches the geographical distribution of the migratory phenotypes. The most differentiated 10 kb windows and missense mutations within these differentiated regions are associated with genes involved in fatty acid synthesis, possibly representing physiological adaptations to the different migratory strategies. The ∼200 genes in these regions, of which several lack described function, will direct future experimental and comparative studies in the search for genes that underlie important migratory traits.

The genomic mosaicism of hybrid speciation

Hybridization is widespread in nature and, in some instances, can result in the formation of a new hybrid species. We investigate the genetic foundation of this poorly understood process through whole-genome analysis of the hybrid Italian sparrow and its progenitors. We find overall balanced yet heterogeneous levels of contribution from each parent species throughout the hybrid genome and identify areas of novel divergence in the hybrid species exhibiting signals consistent with balancing selection. High-divergence areas are disproportionately located on the Z chromosome and overrepresented in gene networks relating to key traits separating the focal species, which are likely involved in reproductive barriers and/or species-specific adaptations. Of special interest are genes and functional groups known to affect body patterning, beak morphology, and the immune system, which are important features of diversification and fitness. We show that a combination of mosaic parental inheritance and novel divergence within the hybrid lineage has facilitated the origin and maintenance of an avian hybrid species.

Intrinsic incompatibilities evolving as a by-product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by-product. A crucial question is to what extent such by-product barriers contribute to divergence and speciation-that is whether they stably reduce gene flow. We discuss under which conditions by-product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by-product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

Gene flow, ancient polymorphism, and ecological adaptation shape the genomic landscape of divergence among Darwin's finches

Genomic comparisons of closely related species have identified “islands” of locally elevated sequence divergence. Genomic islands may contain functional variants involved in local adaptation or reproductive isolation and may therefore play an important role in the speciation process. However, genomic islands can also arise through evolutionary processes unrelated to speciation, and examination of their properties can illuminate how new species evolve. Here, we performed scans for regions of high relative divergence (F_ST) in 12 species pairs of Darwin's finches at different genetic distances. In each pair, we identify genomic islands that are, on average, elevated in both relative divergence (F_ST) and absolute divergence (d_XY). This signal indicates that haplotypes within these genomic regions became isolated from each other earlier than the rest of the genome. Interestingly, similar numbers of genomic islands of elevated d_XY are observed in sympatric and allopatric species pairs, suggesting that recent gene flow is not a major factor in their formation. We find that two of the most pronounced genomic islands contain the ALX1 and HMGA2 loci, which are associated with variation in beak shape and size, respectively, suggesting that they are involved in ecological adaptation. A subset of genomic island regions, including these loci, appears to represent anciently diverged haplotypes that evolved early during the radiation of Darwin's finches. Comparative genomics data indicate that these loci, and genomic islands in general, have exceptionally low recombination rates, which may play a role in their establishment.

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.

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole-genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.

Covariation in levels of nucleotide diversity in homologous regions of the avian genome long after completion of lineage sorting

Closely related species may show similar levels of genetic diversity in homologous regions of the genome owing to shared ancestral variation still segregating in the extant species. However, after completion of lineage sorting, such covariation is not necessarily expected. On the other hand, if the processes that govern genetic diversity are conserved, diversity may potentially covary even among distantly related species. We mapped regions of conserved synteny between the genomes of two divergent bird species-collared flycatcher and hooded crow-and identified more than 600 Mb of homologous regions (66% of the genome). From analyses of whole-genome resequencing data in large population samples of both species we found nucleotide diversity in 200 kb windows to be well correlated (Spearman's ρ = 0.407). The correlation remained highly similar after excluding coding sequences. To explain this covariation, we suggest that a stable avian karyotype and a conserved landscape of recombination rate variation render the diversity-reducing effects of linked selection similar in divergent bird lineages. Principal component regression analysis of several potential explanatory variables driving heterogeneity in flycatcher diversity levels revealed the strongest effects from recombination rate variation and density of coding sequence targets for selection, consistent with linked selection. It is also possible that a stable karyotype is associated with a conserved genomic mutation environment contributing to covariation in diversity levels between lineages. Our observations imply that genetic diversity is to some extent predictable.

Transitions between phases of genomic differentiation during stick-insect speciation

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Inferring Individual Inbreeding and Demographic History from Segments of Identity by Descent in Ficedula Flycatcher Genome Sequences

Individual inbreeding and historical demography can be estimated by analyzing runs of homozygosity (ROH), which are indicative of chromosomal segments of identity by descent (IBD). Such analyses have so far been rare in natural populations due to limited genomic resources. We analyzed ROH in whole genome sequences from 287 Ficedula flycatchers representing four species, with the objectives of evaluating the causes of genome-wide variation in the abundance of ROH and inferring historical demography. ROH were clearly more abundant in genomic regions with low recombination rate. However, this pattern was substantially weaker when ROH were mapped using genetic rather than physical single nucleotide polymorphism (SNP) coordinates in the genome. Empirical results and simulations suggest that high ROH abundance in regions of low recombination was partly caused by increased power to detect the very long IBD segments typical of regions with a low recombination rate. Simulations also showed that hard selective sweeps (but not soft sweeps or background selection) likely contributed to variation in the abundance of ROH across the genome. Comparisons of the abundance of ROH among several study populations indicated that the Spanish pied flycatcher population had the smallest historical effective population size (N_e) for this species, and that a putatively recently founded island (Baltic) population had the smallest historical N_e among the collared flycatchers. Analysis of pairwise IBD in Baltic collared flycatchers indicated that this population was founded <60 generations ago. This study provides a rare genomic glimpse into demographic history and the mechanisms underlying the genome-wide distribution of ROH.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

The Genome of the Trinidadian Guppy, Poecilia reticulata, and Variation in the Guanapo Population

For over a century, the live bearing guppy, Poecilia reticulata, has been used to study sexual selection as well as local adaptation. Natural guppy populations differ in many traits that are of intuitively adaptive significance such as ornamentation, age at maturity, brood size and body shape. Water depth, light supply, food resources and predation regime shape these traits, and barrier waterfalls often separate contrasting environments in the same river. We have assembled and annotated the genome of an inbred single female from a high-predation site in the Guanapo drainage. The final assembly comprises 731.6 Mb with a scaffold N50 of 5.3 MB. Scaffolds were mapped to linkage groups, placing 95% of the genome assembly on the 22 autosomes and the X-chromosome. To investigate genetic variation in the population used for the genome assembly, we sequenced 10 wild caught male individuals. The identified 5 million SNPs correspond to an average nucleotide diversity (π) of 0.0025. The genome assembly and SNP map provide a rich resource for investigating adaptation to different predation regimes. In addition, comparisons with the genomes of other Poeciliid species, which differ greatly in mechanisms of sex determination and maternal resource allocation, as well as comparisons to other teleost genera can begin to reveal how live bearing evolved in teleost fish.

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.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

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.

A Half-Century of Studies on a Chromosomal Hybrid Zone of the House Mouse

The first natural chromosomal variation in the house mouse was described nearly 50 years ago in Val Poschiavo on the Swiss side of the Swiss-Italian border in the Central Eastern Alps. Studies have extended into neighboring Valtellina, and the house mice of the Poschiavo-Valtellina area have been subject to detailed analysis, reviewed here. The maximum extent of this area is 70 km, yet it has 4 metacentric races and the standard 40-chromosome telocentric race distributed in a patchwork fashion. The metacentric races are characterized by highly reduced diploid numbers (2n = 22-26) resulting from Robertsonian fusions, perhaps modified by whole-arm reciprocal translocations. The races hybridize and the whole Poschiavo-Valtellina area can be considered a "hybrid zone." The studies of this area have provided insights into origin of races within hybrid zones, gene flow within hybrid zones and the possibility of speciation in hybrid zones. This provides a case study of how chromosomal rearrangements may impact the genetic structure of populations and their diversification.

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within-group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within-population selection-induced reductions in variation. An important implication of this 'sweep-before-differentiation' model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.