Genomic parallelism and lack thereof in contrasting systems of three-spined sticklebacks

Genomic insights on conservation priorities for North Sea houting and European lake whitefish (Coregonus spp.)

Population genomics analysis holds great potential for informing conservation of endangered populations. We focused on a controversial case of European whitefish (Coregonus spp.) populations. The endangered North Sea houting is the only coregonid fish that tolerates oceanic salinities and was previously considered a species (C. oxyrhinchus) distinct from European lake whitefish (C. lavaretus). However, no firm evidence for genetic-based salinity adaptation has been available. Also, studies based on microsatellite and mitogenome data suggested surprisingly recent divergence (c. 2500 years bp) between houting and lake whitefish. These data types furthermore have provided no evidence for possible inbreeding. Finally, a controversial taxonomic revision recently classified all whitefish in the region as C. maraena, calling conservation priorities of houting into question. We used whole-genome and ddRAD sequencing to analyse six lake whitefish populations and the only extant indigenous houting population. Demographic inference indicated post-glacial expansion and divergence between lake whitefish and houting occurring not long after the Last Glaciation, implying deeper population histories than previous analyses. Runs of homozygosity analysis suggested not only high inbreeding (FROH up to 30.6%) in some freshwater populations but also FROH up to 10.6% in the houting prompting conservation concerns. Finally, outlier scans provided evidence for adaptation to high salinities in the houting. Applying a framework for defining conservation units based on current and historical reproductive isolation and adaptive divergence led us to recommend that the houting be treated as a separate conservation unit regardless of species status. In total, the results underscore the potential of genomics to inform conservation practices, in this case clarifying conservation units and highlighting populations of concern.

Reproductive isolation in a three-way contact zone

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

Integrating Pool-seq uncertainties into demographic inference

Next-generation sequencing of pooled samples (Pool-seq) is a popular method to assess genome-wide diversity patterns in natural and experimental populations. However, Pool-seq is associated with specific sources of noise, such as unequal individual contributions. Consequently, using Pool-seq for the reconstruction of evolutionary history has remained underexplored. Here we describe a novel Approximate Bayesian Computation (ABC) method to infer demographic history, explicitly modelling Pool-seq sources of error. By jointly modelling Pool-seq data, demographic history and the effects of selection due to barrier loci, we obtain estimates of demographic history parameters accounting for technical errors associated with Pool-seq. Our ABC approach is computationally efficient as it relies on simulating subsets of loci (rather than the whole-genome) and on using relative summary statistics and relative model parameters. Our simulation study results indicate Pool-seq data allows distinction between general scenarios of ecotype formation (single versus parallel origin) and to infer relevant demographic parameters (e.g. effective sizes and split times). We exemplify the application of our method to Pool-seq data from the rocky-shore gastropod Littorina saxatilis, sampled on a narrow geographical scale at two Swedish locations where two ecotypes (Wave and Crab) are found. Our model choice and parameter estimates show that ecotypes formed before colonization of the two locations (i.e. single origin) and are maintained despite gene flow. These results indicate that demographic modelling and inference can be successful based on pool-sequencing using ABC, contributing to the development of suitable null models that allow for a better understanding of the genetic basis of divergent adaptation.

Ghost introgression in ricefishes of the genus Adrianichthys in an ancient Wallacean lake

Because speciation might have been promoted by ancient introgression from an extinct lineage, it is important to detect the existence of 'ghost introgression' in focal taxa and examine its contribution to their diversification. In this study, we examined possible ghost introgression and its contributions to the diversification of ricefishes of the genus Adrianichthys in Lake Poso, an ancient lake on Sulawesi Island, in which some extinctions are known to have occurred. Population-genomic analysis revealed that two extant Adrianichthys species, A. oophorus and A. poptae are reproductively isolated from each other. Comparisons of demographic models demonstrated that introgression from a ghost population, which diverged from the common ancestor of A. oophorus and A. poptae, is essential for reconstructing the demographic history of Adrianichthys. The best model estimated that the divergence of the ghost population greatly predated the divergence between A. oophorus and A. poptae, and that the ghost population secondarily contacted the two extant species within Lake Poso more recently. Genome scans and simulations detected a greatly divergent locus, which cannot be explained without ghost introgression. This locus was also completely segregated between A. oophorus and A. poptae. These findings suggest that variants that came from a ghost population have contributed to the divergence between A. oophorus and A. poptae, but the large time-lag between their divergence and ghost introgression indicates that the contribution of introgression may be restricted.

De Novo Mutation Rates in Sticklebacks

Mutation rate is a fundamental parameter in population genetics. Apart from being an important scaling parameter for demographic and phylogenetic inference, it allows one to understand at what rate new genetic diversity is generated and what the expected level of genetic diversity is in a population at equilibrium. However, except for well-established model organisms, accurate estimates of de novo mutation rates are available for a very limited number of organisms from the wild. We estimated mutation rates (µ) in two marine populations of the nine-spined stickleback (Pungitius pungitius) with the aid of several 2- and 3-generational family pedigrees, deep (>50×) whole-genome resequences and a high-quality reference genome. After stringent filtering, we discovered 308 germline mutations in 106 offspring translating to µ = 4.83 × 10-9 and µ = 4.29 × 10-9 per base per generation in the two populations, respectively. Up to 20% of the mutations were shared by full-sibs showing that the level of parental mosaicism was relatively high. Since the estimated µ was 3.1 times smaller than the commonly used substitution rate, recalibration with µ led to substantial increase in estimated divergence times between different stickleback species. Our estimates of the de novo mutation rate should provide a useful resource for research focused on fish population genetics and that of sticklebacks in particular.

Plate reduction in southern Japanese freshwater populations of threespine stickleback (Gasterosteus aculeatus)

Adaptation to similar environments can lead to the evolution of similar phenotypes in phylogenetically independent lineages. However, the extent of parallel evolution often varies. Because such variations can be due to environmental heterogeneity among seemingly similar habitats, identification of the environmental factors that cause non-parallel patterns can provide valuable insight into the ecological factors associated with phenotypic diversification. Armor plate reduction in replicate freshwater populations of the threespine stickleback (Gasterosteus aculeatus) represents a well-known example of parallel evolution. Many freshwater populations in multiple regions of the Northern Hemisphere have reduced plate numbers, but not all freshwater populations exhibit plate reduction. In this study, we characterized plate number variation in Japanese freshwater populations and investigated the association between plate number and several abiotic environmental factors. We found that most freshwater populations have not reduced plate numbers in Japan. Plate reduction tends to occur in habitats with warmer winter temperatures at lower latitudes in Japan. In contrast, low dissolved calcium levels or water turbidity had no significant effects on plate reduction, although these were reported to be associated with plate reduction in Europe. Although our data are consistent with the hypothesis that winter temperatures are associated with plate reduction, further studies on the relationship between temperatures and fitness using sticklebacks with varying plate numbers are necessary to confirm this hypothesis and understand the factors causing variations in the extent of parallel evolution.

Multiple colonizations and hybridization of a freshwater fish group on a satellite island of Sulawesi

Repeated colonizations and resultant hybridization may increase lineage diversity on an island if introgression occurs only in a portion of the indigenous island lineage. Therefore, to precisely understand how island biodiversity was shaped, it is essential to reconstruct the history of secondary colonization and resultant hybridization both in time and space. In this study, we reconstructed the history of multiple colonizations of the Oryzias woworae species group, a freshwater fish group of the family Adrianichthyidae, from Sulawesi Island to its southeast satellite island, Muna Island. Phylogenetic and species tree analyses using genome-wide single-nucleotide polymorphisms revealed that all local populations on Muna Island were monophyletic, but that there were several genetically distinct lineages within the island. Population structure and phylogenetic network analyses demonstrated that colonization of this island occurred more than once, and that secondary colonization and resultant introgressive hybridization occurred only in one local population on the island. The spatially heterogeneous introgression induced by the multiple colonizations were also supported by differential admixture analyses. In addition, the differential admixture analyses detected reverse colonization from Muna Island to the Sulawesi mainland. Coalescence-based demographic inference estimated that these mutual colonizations occurred during the middle to late Quaternary period, during which sea level repeatedly declined; this indicates that the colonizations occurred via land bridges. We conclude that these mutual colonizations between Muna Island and the Sulawesi mainland, and the resultant spatially heterogeneous introgression shaped the current biodiversity of this species group in this area.

Genomic responses to parallel temperature gradients in the eelgrass Zostera marina in adjacent bays

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

Comparative phylogeography of two commensal rat species (Rattus tanezumi and Rattus norvegicus) in China: Insights from mitochondrial DNA, microsatellite, and 2b-RAD data

Rattus norvegicus and Rattus tanezumi are dominant species of Chinese house rats, but the colonization and demographic history of two species in China have not been thoroughly explored. Phylogenetic analyses with mitochondrial DNA including 486 individuals from 31 localities revealed that R. norvegicus is widely distributed in China, R. tanezumi is mainly distributed in southern China with currently invading northward; northeast China was the natal region of R. norvegicus, while the spread of R. tanezumi in China most likely started from the southeast coast. A total of 123 individuals from 18 localities were subjected to 2b-RAD analyses. In neighbor-joining tree, individuals of R. tanezumi grouped into geographic-specific branches, and populations from southeast coast were ancestral groups, which confirmed the colonization route from southeast coast to central and western China. However, individuals of R. norvegicus were generally grouped into two clusters instead of geographic-specific branches. One cluster comprised inland populations, and another cluster included both southeast coast and inland populations, which indicated that spread history of R. norvegicus in China was complex; in addition to on-land colonization, shipping transportation also have played great roles. ADMIXTURE and principal component analyses provided further supports for the colonization history. Demographic analyses revealed that climate changes at ~40,000 to 18,000 years ago and ~4000 years ago had led to population declines of both species; the R. norvegicus declined rapidly while the population of R. tanezumi continuously expanded since ~1500 years ago, indicating the importance of interspecies' competition in their population size changes. Our study provided a valuable framework for further investigation on phylogeography of two species in China.

Compiling forty years of guppy research to investigate the factors contributing to (non)parallel evolution

Examples of parallel evolution have been crucial for our understanding of adaptation via natural selection. However, strong parallelism is not always observed even in seemingly similar environments where natural selection is expected to favour similar phenotypes. Leveraging this variation in parallelism within well-researched study systems can provide insight into the factors that contribute to variation in adaptive responses. Here we analyse the results of 36 studies reporting 446 average trait values in Trinidadian guppies, Poecilia reticulata, from different predation regimes. We examine how the extent of predator-driven phenotypic parallelism is influenced by six factors: sex, trait type, rearing environment, ecological complexity, evolutionary history, and time since colonization. Analyses show that parallel evolution in guppies is highly variable and weak on average, with only 24.7% of the variation among populations being explained by predation regime. Levels of parallelism appeared to be especially weak for colour traits, and parallelism decreased with increasing complexity of evolutionary history (i.e., when estimates of parallelism from populations within a single drainage were compared to estimates of parallelism from populations pooled between two major drainages). Suggestive - but not significant - trends that warrant further research include interactions between the sexes and different trait categories. Quantifying and accounting for these and other sources of variation among evolutionary 'replicates' can be leveraged to better understand the extent to which seemingly similar environments drive parallel and nonparallel aspects of phenotypic divergence.

Genotyping-by-Sequencing Reveals the Impact of Restocking on Wild Common Carp Populations of the Southern Caspian Basin

Understanding the population structure and level of genetic diversity of wild populations is fundamental for appropriate stock management and species conservation. The common carp (Cyprinus carpio) is one of the most important bony fish throughout the Southern coastline of the Caspian Sea, but captures of this species have seen a dramatic reduction during the last decade. As a consequence, a restocking program has been put in place to maintain C. carpio populations, but its impact is not clear. In the present study, the population structure and genetic diversity of C. carpio in the Southern Caspian basin was determined using 17,828 single-nucleotide polymorphism (SNP) markers. A total of 117 individuals collected from four different locations in the southern Caspian basin and a farm were genotyped by genotyping-by-sequencing. The overall Fst obtained was 0.04, indicating a low level of differentiation between populations, and most genetic diversity was attributed to within population variation (97%). The low Fst values suggest that frequent migration events between different locations occur, and three migration events were inferred in the present study. However, each population still showed a distinct genetic profile, which allowed distinguishing the origin of the fish. This indicates that the ongoing restocking program is maintaining the differences between populations to some extent. Nonetheless, high inbreeding and low heterozygosity were detected in all populations, suggesting that additional conservation efforts are required to protect C. carpio populations in the Southern coast of the Caspian Sea.

The first steps toward a global pandemic: Reconstructing the demographic history of parasite host switches in its native range

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real-time, and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after being brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploiting the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

Life histories as mosaics: plastic and genetic components differ among traits that underpin life-history strategies

Life‐history phenotypes emerge from clusters of traits that are the product of genes and phenotypic plasticity. If the impact of the environment differs substantially between traits, then life histories might not evolve as a cohesive whole. We quantified the sensitivity of components of the life history to food availability, a key environmental difference in the habitat occupied by contrasting ecotypes, for 36 traits in fast‐ and slow‐reproducing Trinidadian guppies. Our dataset included six putatively independent origins of the slow‐reproducing, derived ecotype. Traits varied substantially in plastic and genetic control. Twelve traits were influenced only by food availability (body lengths, body weights), five only by genetic differentiation (interbirth intervals, offspring sizes), 10 by both (litter sizes, reproductive timing), and nine by neither (fat contents, reproductive allotment). Ecotype‐by‐food interactions were negligible. The response to low food was aligned with the genetic difference between high‐ and low‐food environments, suggesting that plasticity was adaptive. The heterogeneity among traits in environmental sensitivity and genetic differentiation reveals that the components of the life history may not evolve in concert. Ecotypes may instead represent mosaics of trait groups that differ in their rate of evolution.

Cast away in the Adriatic: Low degree of parallel genetic differentiation in three-spined sticklebacks

The three-spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b-RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average F_ST  ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine-freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates.

A transposon-mediated reciprocal translocation promotes environmental adaptation but compromises domesticability of wild soybeans

Large structural variations frequently occur in higher plants; however, the impact of such variations on plant diversification, adaptation and domestication remains elusive. Here, we mapped and characterised a reciprocal chromosomal translocation in soybeans and assessed its effects on diversification and adaptation of wild (Glycine soja) and semiwild (Glycine gracilis) soybeans, and domestication of cultivated soybean (Glycine max), by tracing the distribution of the translocation in the USDA Soybean Germplasm Collection and population genetics analysis. We demonstrate that the translocation occurred through CACTA transposon-mediated chromosomal breakage in wild soybean c. 0.34 Ma and is responsible for semisterility in translocation heterozygotes and reduces their reproductive fitness. The translocation has differentiated Continental (i.e. China and Russia) populations from Maritime (i.e. Korea and Japan) populations of G. soja and predominately adapted to cold and dry climates. Further analysis revealed that the divergence of G. max from G. soja predates the translocation event and that G. gracilis is an evolutionary intermediate between G. soja and G. max. Our results highlight the effects of a chromosome rearrangement on the processes leading to plant divergence and adaptation, and provides evidence that suggests G. gracilis, rather than G. soja, as the ancestor of cultivated soybean.

Highly Replicated Evolution of Parapatric Ecotypes

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

Species divergence and repeated ancient hybridization in a Sulawesian lake system

An increasing volume of empirical studies demonstrated that hybridization between distant lineages may have promoted speciation in various taxa. However, the timing, extent and direction of introgressive hybridization remain unknown in many cases. Here, we report a possible case in which repeated hybridization promoted divergence of Oryzias ricefishes (Adrianichthyidae) on Sulawesi, an island of Wallacea. Four Oryzias species are endemic to the Malili Lake system in central Sulawesi, which is composed of five tectonic lakes; of these, one lake is inhabited by two species. Morphological and population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that these two sympatric species are phylogenetically sister to but substantially reproductively isolated from each other. Analyses of admixture and comparison of demographic models revealed that the two sympatric species experienced several substantial introgressions from outgroup populations that probably occurred soon after they had secondary contact with each other in the lake. However, the ratio of migrants from the outgroups was estimated to be different between the two species, which is consistent with the hypothesis that these introgressions aided their divergence or prevented them from forming a hybrid swarm. Repeated lake fragmentations and fusions may have promoted diversification of this freshwater fish species complex that is endemic to this ancient lake system.

Genetic diversity and selection signatures in maize landraces compared across 50 years of in situ and ex situ conservation

Genomics-based, longitudinal comparisons between ex situ and in situ agrobiodiversity conservation strategies can contribute to a better understanding of their underlying effects. However, landrace designations, ambiguous common names, and gaps in sampling information complicate the identification of matching ex situ and in situ seed lots. Here we report a 50-year longitudinal comparison of the genetic diversity of a set of 13 accessions from the state of Morelos, Mexico, conserved ex situ since 1967 and retrieved in situ from the same donor families in 2017. We interviewed farmer families who donated in situ landraces to understand their germplasm selection criteria. Samples were genotyped by sequencing, producing 74,739 SNPs. Comparing the two sample groups, we show that ex situ and in situ genome-wide diversity was similar. In situ samples had 3.1% fewer SNPs and lower pairwise genetic distances (F_st 0.008-0.113) than ex situ samples (F_st 0.031-0.128), but displayed the same heterozygosity. Despite genome-wide similarities across samples, we could identify several loci under selection when comparing in situ and ex situ seed lots, suggesting ongoing evolution in farmer fields. Eight loci in chromosomes 3, 5, 6, and 10 showed evidence of selection in situ that could be related with farmers' selection criteria surveyed with focus groups and interviews at the sampling site in 2017, including wider kernels and larger ear size. Our results have implications for ex situ collection resampling strategies and the in situ conservation of threatened landraces.

Ancient and modern stickleback genomes reveal the demographic constraints on adaptation

Adaptation is typically studied by comparing modern populations with contrasting environments. Individuals persisting in the ancestral habitat are typically used to represent the ancestral founding population; however, it has been questioned whether these individuals are good proxies for the actual ancestors.¹ To address this, we applied a paleogenomics approach² to directly access the ancestral genepool: partially sequencing the genomes of two 11- to 13,000-year-old stickleback recovered from the transitionary layer between marine and freshwater sediments of two Norwegian isolation lakes³ and comparing them with 30 modern stickleback genomes from the same lakes and adjacent marine fjord, in addition to a global dataset of 20 genomes.⁴ The ancient stickleback shared genome-wide ancestry with the modern fjord population, whereas modern lake populations have lost substantial ancestral variation following founder effects, and subsequent drift and selection. Freshwater-adaptive alleles found in one ancient stickleback genome have not risen to high frequency in the present-day population from the same lake. Comparison to the global dataset suggested incomplete adaptation to freshwater in our modern lake populations. Our findings reveal the impact of population bottlenecks in constraining adaptation due to reduced efficacy of selection on standing variation present in founder populations.

Population structure limits parallel evolution in sticklebacks

Population genetic theory predicts that small effective population sizes (N_e) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.

An island-hopping bird reveals how founder events shape genome-wide divergence

When populations colonize new areas, both strong selection and strong drift can be experienced due to novel environments and small founding populations, respectively. Empirical studies have predominantly focused on the phenotype when assessing the role of selection, and limited neutral-loci when assessing founder-induced loss of diversity. Consequently, the extent to which processes interact to influence evolutionary trajectories is difficult to assess. Genomic-level approaches provide the opportunity to simultaneously consider these processes. Here, we examine the roles of selection and drift in shaping genomic diversity and divergence in historically documented sequential island colonizations by the silvereye (Zosterops lateralis). We provide the first empirical demonstration of the rapid appearance of highly diverged genomic regions following population founding, the position of which are highly idiosyncratic. As these regions rarely contained loci putatively under selection, it is most likely that these differences arise via the stochastic nature of the founding process. However, selection is required to explain rapid evolution of larger body size in insular silvereyes. Reconciling our genomic data with these phenotypic patterns suggests there may be many genomic routes to the island phenotype, which vary across populations. Finally, we show that accelerated divergence associated with multiple founding steps is the product of genome-wide rather than localized differences, and that diversity erodes due to loss of rare alleles. However, even multiple founder events do not result in divergence and diversity levels seen in evolutionary older subspecies, and therefore do not provide a shortcut to speciation as proposed by founder-effect speciation models.

Heterogeneous genetic basis of age at maturity in salmonid fishes

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10^-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

Bridgehead effect and multiple introductions shape the global invasion history of a termite

Native to eastern Asia, the Formosan subterranean termite Coptotermes formosanus (Shiraki) is recognized as one of the 100 worst invasive pests in the world, with established populations in Japan, Hawaii and the southeastern United States. Despite its importance, the native source(s) of C. formosanus introductions and their invasive pathway out of Asia remain elusive. Using ~22,000 SNPs, we retraced the invasion history of this species through approximate Bayesian computation and assessed the consequences of the invasion on its genetic patterns and demography. We show a complex invasion history, where an initial introduction to Hawaii resulted from two distinct introduction events from eastern Asia and the Hong Kong region. The admixed Hawaiian population subsequently served as the source, through a bridgehead, for one introduction to the southeastern US. A separate introduction event from southcentral China subsequently occurred in Florida showing admixture with the first introduction. Overall, these findings further reinforce the pivotal role of bridgeheads in shaping species distributions in the Anthropocene and illustrate that the global distribution of C. formosanus has been shaped by multiple introductions out of China, which may have prevented and possibly reversed the loss of genetic diversity within its invasive range.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.