The impact of selection, gene flow and demographic history on heterogeneous genomic divergence: three-spine sticklebacks in divergent environments

Repeated evolution on oceanic islands: comparative genomics reveals species-specific processes in birds

Understanding the interplay between genetic drift, natural selection, gene flow, and demographic history in driving phenotypic and genomic differentiation of insular populations can help us gain insight into the speciation process. Comparing patterns across different insular taxa subjected to similar selective pressures upon colonizing oceanic islands provides the opportunity to study repeated evolution and identify shared patterns in their genomic landscapes of differentiation. We selected four species of passerine birds (Common Chaffinch Fringilla coelebs/canariensis, Red-billed Chough Pyrrhocorax pyrrhocorax, House Finch  Haemorhous mexicanus and Dark-eyed/island Junco Junco hyemalis/insularis) that have both mainland and insular populations. Changes in body size between island and mainland populations were consistent with the island rule. For each species, we sequenced whole genomes from mainland and insular individuals to infer their demographic history, characterize their genomic differentiation, and identify the factors shaping them. We estimated the relative (Fst) and absolute (dxy) differentiation, nucleotide diversity (π), Tajima's D, gene density and recombination rate. We also searched for selective sweeps and chromosomal inversions along the genome. All species shared a marked reduction in effective population size (Ne) upon island colonization. We found diverse patterns of differentiated genomic regions relative to the genome average in all four species, suggesting the role of selection in island-mainland differentiation, yet the lack of congruence in the location of these regions indicates that each species evolved differently in insular environments. Our results suggest that the genomic mechanisms involved in the divergence upon island colonization-such as chromosomal inversions, and historical factors like recurrent selection-differ in each species, despite the highly conserved structure of avian genomes and the similar selective factors involved. These differences are likely influenced by factors such as genetic drift, the polygenic nature of fitness traits and the action of case-specific selective pressures.

Population genomic structure of the sea urchin Diadema africanum, a relevant species in the rocky reef systems across the Macaronesian archipelagos

The sea urchin Diadema africanum is a macro-herbivore found in the rocky reef systems of the West African region and Macaronesian archipelagos. Over several decades, high densities of this species have generated marine barrens in certain areas at the Canary Islands. In contrast, more recently, during the last few years, the species has suffered mass mortality events that continue to the present day. In this study, we used 9,109 Single Nucleotide Polymorphisms (SNPs) and a fragment of a mitochondrial gene to evaluate the species' population structure, effects of mass mortalities on its diversity, and potential local adaptation across the Canary Islands and Cabo Verde. Our research provides compelling evidence of low genomic diversity and homogeneity across the studied area for neutral markers, along with recent demographic fluctuations. The high connectivity among distant areas likely allows a rapid recovering of the populations from local mortality events. Interestingly, we also observed genomic sub-structure from 405 SNPs identified as candidate loci under selection for seawater temperature. The lack of divergence among distant sites and the low diversity found can be attributed to the species' divergence from a small ancestral genomic pool, followed by a contemporary demographic expansion, and ongoinggene flow.

Microevolutionary change in wild stickleback: Using integrative time-series data to infer responses to selection

A central goal in evolutionary biology is to understand how different evolutionary processes cause trait change in wild populations. However, quantifying evolutionary change in the wild requires linking trait change to shifts in allele frequencies at causal loci. Nevertheless, datasets that allow for such tests are extremely rare and existing theoretical approaches poorly account for the evolutionary dynamics that likely occur in ecological settings. Using a decade-long integrative phenome-to-genome time-series dataset on wild threespine stickleback (Gasterosteus aculeatus), we identified how different modes of selection (directional, episodic, and balancing) drive microevolutionary change in correlated traits over time. Most strikingly, we show that feeding traits changed by as much 25% across 10 generations which was driven by changes in the genetic architecture (i.e., in both genomic breeding values and allele frequencies at genetic loci for feeding traits). Importantly, allele frequencies at genetic loci related to feeding traits changed at a rate greater than expected under drift, suggesting that the observed change was a result of directional selection. Allele frequency dynamics of loci related to swimming traits appeared to be under fluctuating selection evident in periodic population crashes in this system. Our results show that microevolutionary change in a wild population is characterized by different modes of selection acting simultaneously on different traits, which likely has important consequences for the evolution of correlated traits. Our study provides one of the most thorough descriptions to date of how microevolutionary processes result in trait change in a natural population.

Genomics reveals repeated landlocking of diadromous fish on an isolated island

Landlocking of diadromous fish in freshwater systems can have significant genomic consequences. For instance, the loss of the migratory life stage can dramatically reduce gene flow across populations, leading to increased genetic structuring and stronger effects of local adaptation. These genomic consequences have been well-studied in some mainland systems, but the evolutionary impacts of landlocking in island ecosystems are largely unknown. In this study, we used a genotyping-by-sequencing (GBS) approach to examine the evolutionary history of landlocking in common smelt (Retropinna retropinna) on Chatham Island, a small isolated oceanic island 800 kilometres east of mainland New Zealand. We examined the relationship between Chatham Island and mainland smelt and used coalescent analyses to test the number and timing of landlocking events on Chatham Island. Our genomic analysis, based on 21,135 SNPs across 169 individuals, revealed that the Chatham Island smelt was genomically distinct from the mainland New Zealand fish, consistent with a single ancestral colonisation event of Chatham Island in the Pleistocene. Significant genetic structure was also evident within the Chatham Island smelt, with a diadromous Chatham Island smelt group, along with three geographically structured landlocked groups. Coalescent demographic analysis supported three independent landlocking events, with this loss of diadromy significantly pre-dating human colonisation. Our results illustrate how landlocking of diadromous fish can occur repeatedly across a narrow spatial scale, and highlight a unique system to study the genomic basis of repeated adaptation.

Genome architecture underlying salinity adaptation in the invasive copepod Eurytemora affinis species complex: A review

With climate change, habitat salinity is shifting rapidly throughout the globe. In addition, many destructive freshwater invaders are recent immigrants from saline habitats. Recently, populations of the copepod Eurytemora affinis species complex have invaded freshwater habitats multiple times independently from saline estuaries on three continents. This review discusses features of this species complex that could enhance their evolutionary potential during rapid environmental change. Remarkably, across independent freshwater invasions, natural selection has repeatedly favored the same alleles far more than expected. This high degree of parallelism is surprising, given the expectation of nonparallel evolution for polygenic adaptation. Factors such as population structure and the genome architecture underlying critical traits under selection might help drive rapid adaptation and parallel evolution. Given the preponderance of saline-to-freshwater invasions and climate-induced salinity change, the principles found here could provide invaluable insights into mechanisms operating in other systems and the potential for adaptation in a changing planet.

Range-wide and temporal genomic analyses reveal the consequences of near-extinction in Swedish moose

Ungulate species have experienced severe declines over the past centuries through overharvesting and habitat loss. Even if many game species have recovered thanks to strict hunting regulation, the genome-wide impacts of overharvesting are still unclear. Here, we examine the temporal and geographical differences in genome-wide diversity in moose (Alces alces) over its whole range in Sweden by sequencing 87 modern and historical genomes. We found limited impact of the 1900s near-extinction event but local variation in inbreeding and load in modern populations, as well as suggestion of a risk of future reduction in genetic diversity and gene flow. Furthermore, we found candidate genes for local adaptation, and rapid temporal allele frequency shifts involving coding genes since the 1980s, possibly due to selective harvesting. Our results highlight that genomic changes potentially impacting fitness can occur over short time scales and underline the need to track both deleterious and selectively advantageous genomic variation.

Genomic Signatures of Freshwater Adaptation in Pacific Herring (Clupea pallasii)

Pacific herring (Clupea pallasii) is an essential target of commercial fishing in the North Pacific Ocean. Previous studies have suggested the existence of marine and lake ecological forms of this species within its range. The lake ecological form of herring has a shortened life cycle, spending the winter and spawning in brackish waters near the shoreline without long migrations for feeding; it also has a relatively smaller body size than the marine form. Genetic-based studies have shown that brackish water Pacific herring not only can be distinguished as a separate lake ecological form but possibly has its genetic legacy. Here, as part of an ongoing study, using ddRAD-sequencing data for marine and lake ecological forms from a total of 54 individuals and methods of comparative bioinformatics, we describe genomic signatures of freshwater adaptivity in Pacific herring. In total, 253 genes containing discriminating SNPs were found, and part of those genes was organized into genome clusters, also known as “genomic islands of divergence”. Moreover, the Tajima’s D test showed that these loci are under directional selection in the lake populations of the Pacific herring. Yet, most discriminating loci between the lake and marine ecological forms of Pacific herring do not intersect (by gene name) with those in other known marine fish species with known freshwater/brackish populations. However, some are associated with the same physiological trait—osmoregulation.

Sexually mediated phenotypic variation within and between sexes as a continuum structured by ecology: The mosaic nature of skeletal variation across body regions in Threespine stickleback (Gasterosteus aculeatus L.)

Ecological character displacement between the sexes, and sexual selection, integrate into a convergent set of factors that produce sexual variation. Ecologically modulated, sexually mediated variation within and between sexes may be a major contributor to the amount of total variation that selection can act on in species. Threespine stickleback (Gasterosteus aculeatus) display rapid adaptive responses and sexual variation in many phenotypic traits. We examined phenotypic variation in the skull, pectoral and pelvic girdles of threespine stickleback from two freshwater and two coastal marine sites on the Sunshine Coast of British Columbia, Canada, using an approach that avoids a priori assumptions about bimodal patterns of variation. We quantified shape and size of the cranial, pectoral and pelvic regions of sticklebacks in marine and freshwater habitats using 3D geometric morphometrics and an index of sexually mediated variation. We show that the expression of phenotypic variation is structured in part by the effects of both habitat marine vs freshwater and the effects of individual sites within each habitat. Relative size exerts variable influence, and patterns of phenotypic variation associated with sex vary among body regions. This fine-grained quantification of sexually mediated variation in the context of habitat difference and different anatomical structures indicates a complex relationship between genetically inferred sex and environmental factors, demonstrating that the interplay between shared genetic background and sexually mediated, ecologically based selective pressures structures the phenotypic expression of complex traits.

Adaptation insights from comparative transcriptome analysis of two Opisthopappus species in the Taihang mountains

Opisthopappus is a major wild source of Asteraceae with resistance to cold and drought. Two species of this genus (Opisthopappus taihangensis and O. longilobus) have been employed as model systems to address the evolutionary history of perennial herb biomes in the Taihang Mountains of China. However, further studies on the adaptive divergence processes of these two species are currently impeded by the lack of genomic resources. To elucidate the molecular mechanisms involved, a comparative analysis of these two species was conducted. Among the identified transcription factors, the bHLH members were most prevalent, which exhibited significantly different expression levels in the terpenoid metabolic pathway. O. longilobus showed higher level of expression than did O. taihangensis in terms of terpenes biosynthesis and metabolism, particularly monoterpenoids and diterpenoids. Analyses of the positive selection genes (PSGs) identified from O. taihangensis and O. longilobus revealed that 1203 genes were related to adaptative divergence, which were under rapid evolution and/or have signs of positive selection. Differential expressions of PSG occurred primarily in the mitochondrial electron transport, starch degradation, secondary metabolism, as well as nucleotide synthesis and S-metabolism pathway processes. Several PSGs were obviously differentially expressed in terpenes biosynthesis that might result in the fragrances divergence between O. longilobus and O. taihangensis, which would provide insights into adaptation of the two species to different environments that characterized by sub-humid warm temperate and temperate continental monsoon climates. The comparative analysis for these two species in Opisthopappus not only revealed how the divergence occurred from molecular perspective, but also provided novel insights into how differential adaptations occurred in Taihang Mountains.

Repeated Genetic Targets of Natural Selection Underlying Adaptation of Fishes to Changing Salinity

Ecological transitions across salinity boundaries have led to some of the most important diversification events in the animal kingdom, especially among fishes. Adaptations accompanying such transitions include changes in morphology, diet, whole-organism performance, and osmoregulatory function, which may be particularly prominent since divergent salinity regimes make opposing demands on systems that maintain ion and water balance. Research in the last decade has focused on the genetic targets underlying such adaptations, most notably by comparing populations of species that are distributed across salinity boundaries. Here, we synthesize research on the targets of natural selection using whole-genome approaches, with a particular emphasis on the osmoregulatory system. Given the complex, integrated and polygenic nature of this system, we expected that signatures of natural selection would span numerous genes across functional levels of osmoregulation, especially salinity sensing, hormonal control, and cellular ion exchange mechanisms. We find support for this prediction: genes coding for V-type, Ca2+, and Na+/K+-ATPases, which are key cellular ion exchange enzymes, are especially common targets of selection in species from six orders of fishes. This indicates that while polygenic selection contributes to adaptation across salinity boundaries, changes in ATPase enzymes may be of particular importance in supporting such transitions.

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.

The genomics of phenotypically differentiated Asellus aquaticus cave, surface stream and lake ecotypes

Organisms well suited for the study of ecotype formation have wide distribution ranges, where they adapt to multiple drastically different habitats repeatedly over space and time. Here we study such ecotypes in a Crustacean model, Asellus aquaticus, a commonly occurring isopod found in freshwater habitats as diverse as streams, caves and lakes. Previous studies focusing on cave vs. surface ecotypes have attributed depigmentation, eye loss and prolonged antennae to several south European cave systems. Likewise, surveys across multiple Swedish lakes have identified the presence of dark-pigmented "reed" and light-pigmented "stonewort" ecotypes, which can be found within the same lake. In this study, we sequenced the first draft genome of A. aquaticus, and subsequently use this to map reads and call variants in surface stream, cave and two lake ecotypes. In addition, the draft genome was combined with a RADseq approach to perform a quantitative trait locus (QTL) mapping study using a laboratory bred F₂ and F₄ cave × surface intercross. We identified genomic regions associated with body pigmentation, antennae length and body size. Furthermore, we compared genome-wide differentiation between natural populations and found several genes potentially associated with these habitats. The assessment of the cave QTL regions in the light-dark comparison of lake populations suggests that the regions associated with cave adaptation are also involved with genomic differentiation in the lake ecotypes. These demonstrate how troglomorphic adaptations can be used as a model for related ecotype formation.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Local environment-driven adaptive evolution in a marine invasive ascidian (Molgula manhattensis)

Elucidating molecular mechanisms of environment-driven adaptive evolution in marine invaders is crucial for understanding invasion success and further predicting their future invasions. Although increasing evidence suggests that adaptive evolution could contribute to organisms' adaptation to varied environments, there remain knowledge gaps regarding how environments influence genomic variation in invaded habitats and genetic bases underlying local adaptation for most marine invaders. Here, we performed restriction-site-associated DNA sequencing (RADseq) to assess population genetic diversity and further investigate genomic signatures of local adaptation in the marine invasive ascidian, Molgula manhattensis. We revealed that most invasive populations exhibited significant genetic differentiation, low recent gene flow, and no signal of significant population bottleneck. Based on three genome scan approaches, we identified 109 candidate loci potentially under environmental selection. Redundancy analysis and variance partitioning analysis suggest that local environmental factors, particularly the salinity-related variables, represent crucial evolutionary forces in driving adaptive divergence. Using the newly developed transcriptome as a reference, 14 functional genes were finally obtained with potential roles in salinity adaptation, including SLC5A1 and SLC9C1 genes from the solute carrier gene (SLC) superfamily. Our findings confirm that differed local environments could rapidly drive adaptive divergence among invasive populations and leave detectable genomic signatures in marine invaders.

Identification of putative key genes for thermal adaptation in the Japanese mantis shrimp (Oratosquilla oratoria) through population genomic analysis

Little is known about the mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different habitat temperature. Here, the genome-wide genetic information of three O. oratoria populations were obtained by IIB restriction site-associated DNA (2b-RAD) sequencing and 2403 single-nucleotide polymorphisms (SNPs) were identified. Based on the 2403 SNPs, we found a remarkable genetic differentiation between the Yellow Sea and the East China Sea groups of O. oratoria. Furthermore, 63 SNPs are thought to be associated with different sea temperatures. Based on the 63 SNPs, it is hypothesised that the long-term temperature differences may contribute to the variation of genes associated with multiple biological functions, such as material metabolism, cytoskeleton, cellular processes, inflammatory response and hormonal regulation. This study provides new information for elucidating the molecular mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different temperature.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

Landscape genetic structure of Scirpus mariqueter reveals a putatively adaptive differentiation under strong gene flow in estuaries

Estuarine organisms grow in highly heterogeneous habitats, and their genetic differentiation is driven by selective and neutral processes as well as population colonization history. However, the relative importance of the processes that underlie genetic structure is still puzzling. Scirpus mariqueter is a perennial grass almost limited in the Changjiang River estuary and its adjacent Qiantang River estuary. Here, using amplified fragment length polymorphism (AFLP), a moderate-high level of genetic differentiation among populations (range F ST: 0.0310-0.3325) was showed despite large ongoing dispersal. FLOCK assigned all individuals to 13 clusters and revealed a complex genetic structure. Some genetic clusters were limited in peripheries compared with very mixing constitution in center populations, suggesting local adaptation was more likely to occur in peripheral populations. 21 candidate outliers under positive selection were detected, and further, the differentiation patterns correlated with geographic distance, salinity difference, and colonization history were analyzed with or without the outliers. Combined results of AMOVA and IBD based on different dataset, it was found that the effects of geographic distance and population colonization history on isolation seemed to be promoted by divergent selection. However, none-liner IBE pattern indicates the effects of salinity were overwhelmed by spatial distance or other ecological processes in certain areas and also suggests that salinity was not the only selective factor driving population differentiation. These results together indicate that geographic distance, salinity difference, and colonization history co-contributed in shaping the genetic structure of S. mariqueter and that their relative importance was correlated with spatial scale and environment gradient.

Historical demography of four gecko species specializing in boulder cave habitat: Implications in the evolutionary dead end hypothesis and conservation

Specialization in narrow ecological niches may not only help species to survive in competitive or unique environments but also contribute to their extermination over evolutionary time. Although the "evolutionary dead end" hypothesis has long been debated, empirical evidence from species with detailed information on niche specialization and evolutionary history remains rare. Here we use a group of four closely related Cnemaspis gecko species that depend highly on granite boulder caves in the Mekong Delta to investigate the potential impact of ecological specialization on their evolution and population dynamics. Isolated by unsuitable floodplain habitats, these boulder-dwelling geckos are among the most narrowly distributed Squamata in the world. We applied several coalescence-based approaches combined with the RAD-seq technique to estimate their divergence times, gene flow and demographic fluctuations during the speciation and population differentiation processes. Our results reveal long-term population shrinkage in the four geckos and limited gene flow during their divergence. The results suggest that the erosion and fragmentation of the granite boulder hills have greatly impacted population divergence and declines. The habitat specialization of these geckos has led to fine-scaled speciation in these granite rocky hills; in contrast, specialization might also have pushed these species toward the edge of extinction. Our study also emphasizes the conservation urgency of these vulnerable, cave-dependent geckos.

Elevational divergence in the great tit complex revealed by major hemoglobin genes

Gene flow and demographic history can play important roles in the adaptive genetic differentiation of species, which is rarely understood in the high-altitude adaptive evolution of birds. To elucidate genetic divergence of populations in the great tit complex (Parus major, P. minor and P. cinereus) at different elevations, we compared the genetic structure and gene flow in hemoglobin genes with neutral loci. Our results revealed the elevationally divergent structure of α^A-globin gene, distinctive from that of the β^A-globin gene and neutral loci. We further investigated gene flow patterns among the populations in the central-northern (> 1,000 m a.s.l.), south-eastern (< 1,000 m a.s.l.) regions and the Southwest Mountains (> 2,000 m a.s.l.) in China. The high-altitude (> 1,000 m a.s.l.) diverged α^A-globin genetic structure coincided with higher α^A-globin gene flow between highland populations, in contrast to restricted neutral gene flow concordant with the phylogeny. The higher α^A-globin gene flow suggests the possibility of adaptive evolution during population divergence, contrary to the lower α^A-globin gene flow homogenized by neutral loci during population expansion. In concordance with patterns of historical gene flow, genotypic and allelic profiles provide distinctive patterns of fixation in different high-altitude populations. The fixation of alleles at contrasting elevations may primarily due to highland standing variants α^A49Asn/72Asn/108Ala originating from the south-western population. Our findings demonstrate a pattern of genetic divergence with gene flow in major hemoglobin genes depending on population demographic history.

Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Background

Populations that have repeatedly colonized novel environments are useful for studying the role of ecology in adaptive divergence – particularly if some individuals persist in the ancestral habitat. Such “contemporary ancestors” can be used to demonstrate the effects of selection by comparing phenotypic and genetic divergence between the derived population and their extant ancestors. However, evolution and demography in these “contemporary ancestors” can complicate inferences about the source (standing genetic variation, de novo mutation) and pace of adaptive divergence. Marine threespine stickleback (Gasterosteus aculeatus) have colonized freshwater environments along the Pacific coast of North America, but have also persisted in the marine environment. To what extent are marine stickleback good proxies of the ancestral condition?

Results

We sequenced > 5800 variant loci in over 250 marine stickleback from eight locations extending from Alaska to California, and phenotyped them for platedness and body shape. Pairwise F_ST varied from 0.02 to 0.18. Stickleback were divided into five genetic clusters, with a single cluster comprising stickleback from Washington to Alaska. Plate number, Eda, body shape, and candidate loci showed evidence of being under selection in the marine environment. Comparisons to a freshwater population demonstrated that candidate loci for freshwater adaptation varied depending on the choice of marine populations.

Conclusions

Marine stickleback are structured into phenotypically and genetically distinct populations that have been evolving as freshwater stickleback evolved. This variation complicates their usefulness as proxies of the ancestors of freshwater populations. Lessons from stickleback may be applied to other “contemporary ancestor”-derived population studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1228-8) contains supplementary material, which is available to authorized users.

Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.

Heterogeneous Patterns of Genetic Diversity and Differentiation in European and Siberian Chiffchaff (Phylloscopus collybita abietinus/P. tristis)

Identification of candidate genes for trait variation in diverging lineages and characterization of mechanistic underpinnings of genome differentiation are key steps toward understanding the processes underlying the formation of new species. Hybrid zones provide a valuable resource for such investigations, since they allow us to study how genomes evolve as species exchange genetic material and to associate particular genetic regions with phenotypic traits of interest. Here, we use whole-genome resequencing of both allopatric and hybridizing populations of the European (Phylloscopus collybita abietinus) and the Siberian chiffchaff (P. tristis)—two recently diverged species which differ in morphology, plumage, song, habitat, and migration—to quantify the regional variation in genome-wide genetic diversity and differentiation, and to identify candidate regions for trait variation. We find that the levels of diversity, differentiation, and divergence are highly heterogeneous, with significantly reduced global differentiation, and more pronounced differentiation peaks in sympatry than in allopatry. This pattern is consistent with regional differences in effective population size and recurrent background selection or selective sweeps reducing the genetic diversity in specific regions prior to lineage divergence, but the data also suggest that postdivergence selection has resulted in increased differentiation and fixed differences in specific regions. We find that hybridization and backcrossing is common in sympatry, and that phenotype is a poor predictor of the genomic composition of sympatric birds. The combination of a differentiation scan approach with identification of fixed differences pinpoint a handful of candidate regions that might be important for trait variation between the two species.

Genetic diversity and structure of wild and cultivated Amorphophallus paeoniifolius populations in southwestern China as revealed by RAD-seq

Amorphophallus paeoniifolius, is a commercially important vegetable crop because of its high production potential. In this study, we generated a total of 166 Gb of genomic data from 16 wild and 20 cultivated A. paeoniifolius individuals in southwestern China using restriction site associated DNA sequencing (RAD-seq). We compared the genome-wide variations between the wild and cultivated populations. Wild populations exhibited higher genetic diversity than did cultivated populations based on private allele number, expected heterozygosity, observed heterozygosity and nucleotide diversity. STRUCTURE analysis, principal component analysis (PCA) and a maximum likelihood (ML) tree indicated that A. paeoniifolius populations could be divided into three groups (a cultivated group and two wild groups) with significant genetic differentiation. The low genetic diversity and shallow genetic differentiation found within cultivated populations are likely caused by continuous selection and the clonal propagation methods used during domestication. The significant differentiation between the wild populations may suggest strong genetic drift due to small populations and human disturbance. The genome-wide single nucleotide polymorphisms (SNPs) identified in our study will provide a valuable resource for further breeding improvement and effective use of the germplasm.

Non-parallel divergence across freshwater and marine three-spined stickleback Gasterosteus aculeatus populations

This work investigated whether multiple freshwater populations of three-spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine-freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater-marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype-genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low-plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low-plated allele. Re-examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype-genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used.

Low genetic and phenotypic divergence in a contact zone between freshwater and marine sticklebacks: gene flow constrains adaptation

BACKGROUND

Distinct hybrid zones and phenotypic and genomic divergence is often observed between marine and freshwater threespine sticklebacks (Gasterosteus aculeatus). Nevertheless, cases also exist where marine-freshwater divergence is diffuse despite seemingly similar environmental settings. In order to assess what characterizes these highly different outcomes, we focused on the latter kind of system in the Odder River, Denmark. Here, a previous study based on RAD (Restriction site Associated DNA) sequencing found non-significant genome-wide differentiation between marine and freshwater sticklebacks. In the present study, we analyzed samples on a finer geographical scale. We assessed if the system should be regarded as panmictic, or if fine-scale genetic structure and local selection was present but dominated by strong migration. We also asked if specific population components, that is the two sexes and different lateral plate morphs, contributed disproportionally more to dispersal.

RESULTS

We assessed variation at 96 SNPs and the Eda gene that affects lateral plate number, conducted molecular sex identification, and analyzed morphological traits. Genetic differentiation estimated by F_ST was non-significant throughout the system. Nevertheless, spatial autocorrelation analysis suggested fine scale genetic structure with a genetic patch size of 770 m. There was no evidence for sex-biased dispersal, but full-plated individuals showed higher dispersal than low- and partial-plated individuals. The system was dominated by full-plated morphs characteristic of marine sticklebacks, but in the upstream part of the river body shape and frequency of low-plated morphs changed in the direction expected for freshwater sticklebacks. Five markers including Eda were under possible diversifying selection. However, only subtle clinal patterns were observed for traits and markers.

CONCLUSIONS

We suggest that gene flow from marine sticklebacks overwhelms adaptation to freshwater conditions, but the short genetic patch size means that the effect of gene flow on the most upstream region must be indirect and occurs over generations. The occurrence of both weak unimodal and strong bimodal hybrid zones within the same species is striking. We suggest environmental and demographic factors that could determine these outcomes, but also highlight the possibility that long-term population history and the presence or absence of genomic incompatibilities could be a contributing factor.

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.