Evidence for past and present hybridization in three Antarctic icefish species provides new perspectives on an evolutionary radiation

Determining the timing, extent and underlying causes of interspecific gene exchange during or following speciation is central to understanding species' evolution. Antarctic notothenioid fish, thanks to the acquisition of antifreeze glycoproteins during Oligocene transition to polar conditions, experienced a spectacular radiation to >100 species during Late Miocene cooling events. The impact of recent glacial cycles on this group is poorly known, but alternating warming and cooling periods may have affected species' distributions, promoted ecological divergence into recurrently opening niches and/or possibly brought allopatric species into contact. Using microsatellite markers and statistical methods including Approximate Bayesian Computation, we investigated genetic differentiation, hybridization and the possible influence of the last glaciation/deglaciation events in three icefish species of the genus Chionodraco. Our results provide strong evidence of contemporary and past introgression by showing that: (i) a substantial fraction of contemporary individuals in each species has mixed ancestry, (ii) evolutionary scenarios excluding hybridization or including it only in ancient times have small or zero posterior probabilities, (iii) the data support a scenario of interspecific gene flow associated with the two most recent interglacial periods. Glacial cycles might therefore have had a profound impact on the genetic composition of Antarctic fauna, as newly available shelf areas during the warmer intervals might have favoured secondary contacts and hybridization between diversified groups. If our findings are confirmed in other notothenioids, they offer new perspectives for understanding evolutionary dynamics of Antarctic fish and suggest a need for new predictions on the effects of global warming in this group.     

Population Genetic Studies Revealed Local Adaptation in a High Gene-Flow Marine Fish,the Small Yellow Croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.     

Phylogeography of the Chionodraco genus (Perciformes, Channichthydae) in the Southern Ocean

Antarctic fish of the suborder Notothenioidei represent one of the most notable examples of adaptive radiation in the marine environment. The evolutionary relationships between and within the eight families of this suborder have been well established by numerous studies, whereas the microevolutionary processes of notothenioid species remain largely unexplored. In the present paper we investigated the evolutionary relationships between three closely related species of the genus Chionodraco (family Channichthyidae), namely Chionodraco hamatus, Chionodraco rastrospinosus, and Chionodraco myersi by analysing portions of the mitochondrial genome (D-loop and 16S rRNA). The taxonomic status of C. hamatus and C. rastrospinosus as separate species has been questioned because of the limited number of key morphological characters that distinguish these two taxa. Our results, based on the analysis of several specimens belonging to both morphological groups revealed a small genetic differentiation among haplotypes, however, a clear separation between the two nominal species emerged since all individuals of each of the two taxa clustered together in distinct monophyletic groups. C. myersi appeared more distantly related in the phylogenetic analysis. For one species, C. hamatus, sampling was carried out at three different geographic locations in the area of the Ross Sea and Weddell Sea. The results showed that the partition of the genetic variation within this species is not compatible with the hypothesis of panmixia as gene flow between populations was significantly reduced.     

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (F_ST ~0.1), over a short coastline (<60 km) in the North Sea‐Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post‐glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Detecting population structure in a high gene-flow species,Atlantic herring (Clupea harengus): direct,simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.     

Spatial and ecological population genetic structures within two island‐endemic Aeonium species of different niche width

The Crassulacean genus Aeonium is a well‐known example for plant species radiation on oceanic archipelagos. However, while allopatric speciation among islands is documented for this genus, the role of intra‐island speciation due to population divergence by topographical isolation or ecological heterogeneity has not yet been addressed. The aim of this study was to investigate intraspecific genetic structures and to identify spatial and ecological drivers of genetic population differentiation on the island scale. We analyzed inter simple sequence repeat variation within two island‐endemic Aeonium species of La Palma: one widespread generalist that covers a large variety of different habitat types (Ae. davidbramwellii) and one narrow ecological specialist (Ae. nobile), in order to assess evolutionary potentials on this island. Gene pool differentiation and genetic diversity patterns were associated with major landscape structures in both species, with phylogeographic implications. However, overall levels of genetic differentiation were low. For the generalist species, outlier loci detection and loci–environment correlation approaches indicated moderate signatures of divergent selection pressures linked to temperature and precipitation variables, while the specialist species missed such patterns. Our data point to incipient differentiation among populations, emphasizing that ecological heterogeneity and topographical structuring within the small scales of an island can foster evolutionary processes. Very likely, such processes have contributed to the radiation of Aeonium on the Canary Islands. There is also support for different evolutionary mechanisms between generalist and specialist species.     

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late‐successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.     

Signatures of adaptation and genetic structure among the mainland populations of Pinus radiata (D. Don) inferred from SNP loci

Insights into the relative contributions of locus specific and genome-wide effects on population genetic diversity can be gained through separation of their resulting genetic signals. Here we explore patterns of adaptive and neutral genetic diversity in the disjunct natural populations of Pinus radiata (D. Don) from mainland California. A first-generation common garden of 447 individuals revealed significant differentiation of wood phenotypes among populations (P _ST), possibly reflecting local adaptation in response to environment. We subsequently screened all trees for genetic diversity at 149 candidate gene single nucleotide polymorphism (SNP) loci for signatures of adaptation. Ten loci were identified as being possible targets of diversifying selection following F _ST outlier tests. Multivariate canonical correlation performed on a data set of 444 individuals identified significant covariance between environment, adaptive phenotypes and outlier SNP diversity, lending support to the case for local adaptation suggested from F _ST and P _ST tests. Covariation among discrete sets of outlier SNPs and adaptive phenotypes (inferred from multivariate loadings) with environment are supported by existing studies of candidate gene function and genotype–phenotype association. Canonical analyses failed to detect significant correlations between environment and 139 non-outlier SNP loci, which were applied to estimate neutral patterns of genetic differentiation among populations (F _ST 4.3 %). Using this data set, significant hierarchical structure was detected, indicating three populations on the mainland. The hierarchical relationships based on neutral SNP markers (and SSR) were in contrast with those inferred from putatively adaptive loci, potentially highlighting the independent action of selection and demography in shaping genetic structure in this species.     

Molecular signatures of lineage‐specific adaptive evolution in a unique sea basin: the example of an anadromous goby Leucopsarion petersii

Climate changes on various time scales often shape genetic novelty and adaptive variation in many biotas. We explored molecular signatures of directional selection in populations of the ice goby Leucopsarion petersii inhabiting a unique sea basin, the Sea of Japan, where a wide variety of environments existed in the Pleistocene in relation to shifts in sea level by repeated glaciations. This species consisted of two historically allopatric lineages, the Japan Sea (JS) and Pacific Ocean (PO) lineages, and these have lived under contrasting marine environments that are expected to have imposed different selection regimes caused by past climatic and current oceanographic factors. We applied a limited genome‐scan approach using seven candidate genes for phenotypic differences between two lineages in combination with 100 anonymous microsatellite loci. Neuropeptide Y (NPY) gene, which is an important regulator of food intake and potent orexigenic agent, and three anonymous microsatellites were identified as robust outliers, that is, candidate loci potentially under directional selection, by multiple divergence‐ and diversity‐based outlier tests in comparisons focused on multiple populations of the JS vs. PO lineages. For these outlier loci, populations of the JS lineage had putative signals of selective sweeps. Additionally, real‐time quantitative PCR analysis using fish reared in a common environment showed a higher expression level for NPY gene in the JS lineage. Thus, this study succeeded in identifying candidate genomic regions under selection across populations of the JS lineage and provided evidence for lineage‐specific adaptive evolution in this unique sea basin.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

An outlier locus relevant in habitat-mediated selection in an alpine plant across independent regional replicates

Habitat types can induce genetic responses in species and may drive adaptive differentiation and evolutionary divergence of populations. In this study, we aimed at detecting loci indicative of adaptation for different habitat types in the alpine plant Arabis alpina. We used a dataset consisting of A. alpina plants collected in scree, nutrient-rich and moist habitat types in two independent regional replicates of the European Alps (the Swiss and French Alps). Genome scans resulting in 825 amplified fragment length polymorphisms (AFLPs) followed by outlier analysis, i.e. looking for excessive differentiation between habitat types, after accounting for heterozygosity and population structure, was used to detect loci under divergent selection for habitat type within and across the alpine regions. The outlier analyses resulted in the detection of a consistent single outlier locus, which showed a higher fragment frequency in moist compared to the other habitat types in both alpine regions. In addition, a posteriori tests for hierarchical population structuring in the dataset did not detect signals confounding selection at this locus (i.e. signals of regional population structure). Thus, we consider this locus indicative of habitat-mediated selection, and we subsequently sequence-characterized and compared it to the Arabidopsis genome. The sequence was found to be a putative homologue to the SIT4 phosphatase-associated family protein. The detection of this locus in two alpine regions and the availability of its genome sequence make this locus a strong candidate worth further exploration in the habitat-mediated selection and genetic adaptation of natural populations in the alpine plant A. alpina.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

Genetic Divergence and Signatures of Natural Selection in Marginal Populations of a Keystone,Long-Lived Conifer,Eastern White Pine (Pinus strobus) from Northern Ontario

Marginal populations are expected to provide the frontiers for adaptation, evolution and range shifts of plant species under the anticipated climate change conditions. Marginal populations are predicted to show genetic divergence from central populations due to their isolation, and divergent natural selection and genetic drift operating therein. Marginal populations are also expected to have lower genetic diversity and effective population size (N _e) and higher genetic differentiation than central populations. We tested these hypotheses using eastern white pine (Pinus strobus) as a model for keystone, long-lived widely-distributed plants. All 614 eastern white pine trees, in a complete census of two populations each of marginal old-growth, central old-growth, and central second-growth, were genotyped at 11 microsatellite loci. The central populations had significantly higher allelic and genotypic diversity, latent genetic potential (LGP) and N _e than the marginal populations. However, heterozygosity and fixation index were similar between them. The marginal populations were genetically diverged from the central populations. Model testing suggested predominant north to south gene flow in the study area with curtailed gene flow to northern marginal populations. Signatures of natural selection were detected at three loci in the marginal populations; two showing divergent selection with directional change in allele frequencies, and one balancing selection. Contrary to the general belief, no significant differences were observed in genetic diversity, differentiation, LGP, and N _e between old-growth and second-growth populations. Our study provides information on the dynamics of migration, genetic drift and selection in central versus marginal populations of a keystone long-lived plant species and has broad evolutionary, conservation and adaptation significance.     

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site‐associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.