Find research datasets worth reusing

Data from: Genome-wide patterns of divergence and gene flow across a butterfly radiation

The Heliconius butterflies are a diverse recent radiation comprising multiple levels of divergence with on-going gene flow between species. The recently sequenced genome of Heliconius melpomene allowed us to investigate the genomic evolution of this group using dense RAD marker sequencing. Phylogenetic analysis of 54 individuals robustly supported reciprocal monophyly of H. melpomene and H. cydno and refuted previous phylogenetic hypotheses that H. melpomene may be paraphylectic with respect to H. cydno. H. timareta also formed a monophyletic clade closely related but distinct from H. cydno with H. heurippa falling within this clade. We find evidence for pervasive gene flow between sympatric populations of the sister clades H. melpomene and H. cydno/timareta, particularly between H. cydno and H. melpomene from Central America and between H. timareta and H. melpomene from the eastern slopes of the Andes. Between races, divergence is primarily explained by isolation-by-distance; and there is little, if any, genetic population structure between parapatric races, suggesting that hybrid zones between races are not zones of secondary contact. Our results support previous findings that colour pattern loci are shared between populations and species with similar colour pattern elements. Further this pattern is almost unique to these genomic regions with only a very small number of other loci showing significant similarity between populations and species with similar colour patterns.

opencc-zeroDec 2015View details →

Data from: Genome-wide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River

The lower Congo River (LCR) is a freshwater biodiversity hotspot in Africa characterized by some of the world's largest rapids. However, little is known about the evolutionary forces shaping this diversity, which include numerous endemic fishes. We investigated phylogeographic relationships in Teleogramma, a small clade of rheophilic cichlids, in the context of regional geography and hydrology. Previous studies have been unable to resolve phylogenetic relationships within Teleogramma due to lack of variation in nuclear genes and discrete morphological characters among putative species. To sample more broadly across the genome we analyzed double-digest restriction-associated sequencing (ddRAD) data from 53 individuals across all described species in the genus. We also assessed body shape and mitochondrial variation within and between taxa. Phylogenetic analyses reveal previously unrecognized lineages and instances of microallopatric divergence across as little as ~1.5 km. Species ranges appear to correspond to geographic regions broadly separated by major hydrological and topographic barriers, indicating these features are likely important drivers of diversification. Mitonuclear discordance indicates one or more introgressive hybridization events, but no clear evidence of admixture is present in nuclear genomes, suggesting these events were likely ancient. A survey of female fin patterns hints that previously undetected lineage-specific patterning may be acting to reinforce species cohesion. These analyses highlight the importance of hydrological complexity in generating diversity in certain freshwater systems, as well as the utility of ddRAD-Seq data in understanding diversification processes operating both below and above the species level.

Data from: Genomic tests of the species-pump hypothesis: recent island connectivity cycles drive population divergence but not speciation in Caribbean crickets across the Virgin Islands

Harnessing the power of genomic scans, we test the debated 'species pump' hypothesis that implicates repeated cycles of island connectivity and isolation as drivers of divergence. This question has gone understudied given the limited resolution of past molecular markers for studying such dynamic phenomena. With an average of 32000 SNPs from the genome of 136 individuals from ten populations of a Caribbean flightless ground cricket species (Amphiacusta sanctaecrucis) and a complementary set of statistical approaches, we infer a stepping-stone colonization model and high levels of genetic differentiation across the Virgin Islands, which have been periodically inter-connected until 8 ka. Estimates of divergence times from models based on the site frequency spectrum coincide with a period of repeated connection and fragmentation of the islands at 75–130 ka. These results are consistent with a role of island connectivity cycles in promoting genomic divergence and indicate that the genetic distinctiveness of island populations has persisted despite subsequent and extended interisland connections identified from bathymetric data. We discuss these findings in the broader context of Caribbean biogeography, and more specifically why high levels of genomic divergence across the Virgin Islands associated with repeated connectivity cycles do not actually translate into species diversification.

opencc-zeroDec 2015View details →

Data from: The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear-cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200,000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high-density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F(ST)), but not absolute (d(XY)), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to 'cold-spots' of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.

opencc-zeroDec 2017View details →

Data from: Rapid divergence of genome architectures following the origin of an ectomycorrhizal symbiosis in the genus Amanita

Fungi are evolutionary shape shifters and adapt quickly to new environments. Ectomycorrhizal (EM) symbioses are mutualistic associations between fungi and plants and have evolved repeatedly and independently across the fungal tree of life, suggesting lineages frequently reconfigure genome content to take advantage of open ecological niches. To date analyses of genomic mechanisms facilitating EM symbioses have involved comparisons of distantly related species, but here, we use the genomes of three EM and two asymbiotic (AS) fungi from the genus Amanita as well as an AS outgroup to study genome evolution following a single origin of symbiosis. Our aim was to identify the defining features of EM genomes, but our analyses suggest no clear differentiation of genome size, gene repertoire size or transposable element content between EM and AS species. Phylogenetic inference of gene gains and losses suggests the transition to symbiosis was dominated by the loss of plant cell wall decomposition genes, a confirmation of previous findings. However, the same dynamic defines the AS species A. inopinata, suggesting loss is not strictly associated with origin of symbiosis. Gene expansions in the common ancestor of EM Amanita were modest, but lineage specific and large gene family expansions are found in two of the three EM extant species. Even closely related EM genomes appear to share few common features. The genetic toolkit required for symbiosis appears already encoded in the genomes of saprotrophic species, and this dynamic may explain the pervasive, recurrent evolution of ectomycorrhizal associations.

opencc-zeroDec 2013View details →

Data from: Range-wide multilocus phylogeography of the red fox reveals ancient continental divergence, minimal genomic exchange, and distinct demographic histories

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization timeframe suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y-chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long-term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.

Data from: SNP-array reveals genome wide patterns of geographical and potential adaptive divergence across the natural range of Atlantic salmon (Salmo salar)

Atlantic salmon (Salmo salar) is one of the most extensively studied fish species in the world due to its significance in aquaculture, fisheries and ongoing conservation efforts to protect declining populations. Yet, limited genomic resources have hampered our understanding of genetic architecture in the species and the genetic basis of adaptation to the wide range of natural and artificial environments it occupies. In this paper, we describe the development of a medium density Atlantic salmon SNP-array based on Expressed Sequence Tags (ESTs) and genomic sequencing. The array was used in the most extensive assessment of population genetic structure performed to date in this species. A total of 6176 informative SNPs were successfully genotyped in 38 anadromous and freshwater wild populations distributed across the species natural range. Principal component analysis clearly differentiated European and North American populations, and within Europe, three major regional genetic groups were identified for the first time in a single analysis. We assessed the potential for the array to disentangle neutral and putative adaptive divergence of SNP allele frequencies across populations and among regional groups. In Europe, secondary contact zones were identified between major clusters where endogenous and exogenous barriers could be associated, rendering the interpretation of environmental influence on potentially adaptive divergence equivocal. A small number of markers highly divergent in allele frequencies (outliers) were observed between (multiple) freshwater and anadromous populations, between northern and southern latitudes, and when comparing Baltic populations to all others. We also discuss the potential future applications of the SNP-array for conservation, management and aquaculture.

opencc-zeroDec 2013View details →

Data from: Genomic variation in a widespread Neotropical bird (Xenops minutus) reveals divergence, population expansion, and gene flow

The demographic and phylogeographic histories of species provide insight into the processes responsible for generating biological diversity, and genomic datasets are now permitting the estimation of species histories with unprecedented accuracy. We used a genomic single nucleotide polymorphism (SNP) dataset generated using a RAD-Seq method to investigate the historical demography and phylogeography of a widespread lowland Neotropical bird (Xenops minutus). As expected, we found that prominent landscape features that act as dispersal barriers, such as Amazonian rivers and the Andes Mountains, are associated with the deepest phylogeographic breaks, and also that isolation by distance is limited in areas between these barriers. In addition, we inferred positive population growth for most populations and detected evidence of historical gene flow between populations that are now physically isolated. Although we were able to reconstruct the history of Xenops minutus with unprecedented resolution, we had difficulty conclusively relating this history to the landscape events implicated in many Neotropical diversification hypotheses. We suggest that even if many traditional diversification hypotheses remain untestable, investigations using genomic datasets will provide greater resolution of species histories in the Neotropics and elsewhere.

opencc-zeroDec 2018View details →

Data from: Parallel signatures of selection at genomic islands of divergence and the major histocompatibility complex in ecotypes of sockeye salmon across Alaska

Understanding the genetic mechanisms that facilitate adaptive radiation is an important component of evolutionary biology. Here, we genotyped 82 neutral SNPs, seven SNPs in islands of divergence identified in a previous study (island SNPs), and a region of the major histocompatibility complex (MHC) in 32 populations of sockeye salmon to investigate whether conserved genes and genomic regions are involved in adaptive radiation. Populations representing three ecotypes were sampled from seven drainages with differing habitats and colonization histories spanning a range of 2,000 km. We found strong signatures of parallel selection across drainages at the island SNPs and MHC, suggesting that the same loci undergo divergent selection during adaptive radiation. However, patterns of differentiation at most island SNPs and the MHC were not associated with ecotypes, suggesting that these loci are responding differently to a mosaic of selective pressures. Our study provides some of the first evidence that conserved genomic islands may be involved in adaptive divergence of salmon populations. Additionally, our data provide further support for the hypothesis that sockeye salmon inhabiting rivers unconnected to lakes harbor similar genetic diversity across large distances, are likely the ancestral form of the species, and have repeatedly recolonized lake systems as they have become available after glacial recession. Finally, our results highlight the value and importance of validating outlier loci by screening additional populations and regions, a practice that will hopefully become more common in the future.

Data from: Evaluating genomic divergence and parallelism in replicate ecomorphs from young and old cichlid adaptive radiations

Comparative genomic studies of closely related species typically focus on single species pairs at one given stage of divergence. That makes it difficult to infer the continuum of evolutionary process during speciation and beyond. Here, we use whole-genome resequencing to examine genomic patterns of divergence in three sympatric cichlid species pairs with very similar functional and ecological differentiation, but different ages. We find a strong signature of increasing genomic divergence with time in both the mitochondrial genome and the nuclear genome. In contrast to many other systems, we find that in these cichlids regions of elevated relative differentiation also exhibit increased absolute differentiation. We detect a signature of convergent evolution in a comparison of outlier regions across all three species pair comparisons but the extent of it is modest, and regions that are strongly divergent in any one pair tend to be only slightly elevated in the other pairs, consistent with a repeatable but polygenic basis of traits that characterize the ecomorphs. Our results suggest that strong functional phenotypic differentiation, as seen in all three species pairs, is generally associated with a clear signature of genomic divergence, even in the youngest species pair.

Data from: Genomic architecture of habitat-related divergence and signature of directional selection in the body shapes of Gnathopogon fishes

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat-related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat-related divergence in the body shape of Gnathopogon fishes, a novel example of lake–stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream-dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F2 interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction-site associated DNA sequencing-derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape-related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape-related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.

opencc-zeroDec 2016View details →

Data from: Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome

A whole-genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid-origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid-origin species. We demonstrate a new method that enables genome-wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus-specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid-origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications.

Data from: Genome divergence and the genetic architecture of barriers to gene flow between Lycaeides idas and L. melissa

Genome divergence during speciation is a dynamic process that is affected by various factors, including the genetic architecture of barriers to gene flow. Herein we quantitatively describe aspects of the genetic architecture of two sets of traits, male genitalic morphology and oviposition preference, that putatively function as barriers to gene flow between the butterfly species Lycaeides idas and L. melissa. Our analyses are based on unmapped DNA sequence data and a recently developed Bayesian regression approach that includes variable selection and explicit parameters for the genetic architecture of traits. A modest number of nucleotide polymorphisms explained a small to large proportion of the variation in each trait, and average genetic variant effects were non-negligible. Several genetic regions were associated with variation in multiple traits or with trait variation within- and among-populations. In some instances genetic regions associated with trait variation also exhibited exceptional genetic differentiation between speices or exceptional introgression in hybrids. These results are consistent with the hypothesis that divergent selection on male genitalia has contributed to heterogeneous genetic differentiation, and that both sets of traits affect fitness in hybrids. Although these results are encouraging, we highlight several difficulties related to understanding the genetics of speciation.

opencc-zeroSep 2019View details →

Data from: Parallel genomic architecture underlies repeated sexual signal divergence in Hawaiian Laupala crickets

When the same phenotype evolves repeatedly, we can explore the predictability of genetic changes underlying phenotypic evolution. Theory suggests that genetic parallelism is less likely when phenotypic changes are governed by many small-effect loci compared to few of major effect, because different combinations of genetic changes can result in the same quantitative outcome. However, some genetic trajectories might be favoured over others, making a shared genetic basis to repeated polygenic evolution more likely. To examine this, we studied the genetics of parallel male mating song evolution in the Hawaiian cricket Laupala. We compared quantitative trait loci (QTL) underlying song divergence in three species pairs varying in phenotypic distance. We tested whether replicated song divergence between species involves the same QTL and the likelihood that sharing QTL is related to phenotypic effect sizes. Contrary to theoretical predictions, we find substantial parallelism in polygenic genetic architectures underlying repeated song divergence. QTL overlapped more than expected based on simulated QTL analyses. Interestingly, QTL effect size did not predict QTL sharing, but did correlate with magnitude of phenotypic divergence. We highlight potential mechanisms driving these constraints on cricket song evolution and discuss a scenario that consolidates empirical quantitative genetic observations with micro-mutational theory.

opencc-zeroFeb 2020View details →

Genomics confirms surprising ecological divergence and isolation in an endangered butterfly

Phylogeographic patterns in phytophagous organisms are often contextualized in light of geographic isolation and ecological (host, habitat) specialization. However, assessing the relative impact of these phenomena is not straightforward, even in areas where phylogeography is well-studied, such as the California Floristic Province. Here, we use genome-wide markers to elucidate population genomic and phylgeographic patterns for a group of monophytophagous butterflies in southern California. This group is of high conservation interest because it includes the El Segundo blue, Euphilotes battoides allyni, one of the first insects listed under the U.S. Endangered Species Act, and a newly discovered population putatively assigned to E. b. allyni. Despite using the same unique host and coastal habitat, our results indicate that the newly discovered populations are not E. b. allyni and are more closely related to geographically proximate populations of the E. battoides group using a different habitat host. Aside from E. b. allyni and the newly discovered populations, the rest of the group shows only fine-scale structure and apparently maintains genetic connectivity throughout southern California, across a vast range of habitats and climates, and on multiple hosts. Thus, habitat and host specialization did not elicit genetic isolation in neighboring populations suggesting that: 1) other phenomena are needed to explain the remarkable and idiosyncratic divergence of these highly restricted, proximate, taxa, and 2) fine-scale genomic markers suggest broader implications for understanding the mechanisms of speciation and reinvestigation of phylogeographic patterns in regions like the California Floristic Province.

opencc-zeroDec 2012View details →

Data from: Comparative species divergence across eight triplets of spiny lizards (Sceloporus) using genomic sequence data

Species divergence is typically thought to occur in the absence of gene flow, but many empirical studies are discovering that gene flow may be more pervasive during species formation. Although many examples of divergence with gene flow have been identified, only few clades have been investigated in a comparative manner, and fewer have been studied using genome-wide sequence data. We contrast species divergence genetic histories across eight triplets of North American Sceloporus lizards using a maximum likelihood implementation of the isolation–migration (IM) model. Gene flow at the time of species divergence is modeled indirectly as variation in species divergence time across the genome or explicitly using a migration rate parameter. Likelihood ratio tests (LRTs) are used to test the null model of no gene flow at speciation against these two alternative gene flow models. We also use the Akaike information criterion to rank the models. Hundreds of loci are needed for the LRTs to have statistical power, and we use genome sequencing of reduced representation libraries to obtain DNA sequence alignments at many loci (between 340 and 3,478; mean 1⁄4 1,678) for each triplet. We find that current species distributions are a poor predictor of whether a species pair diverged with gene flow. Interrogating the genome using the triplet method expedites the comparative study of species divergence history and the estimation of genetic parameters associated with speciation.

opencc-zeroDec 2015View details →

Data from: Adaptive evolution and segregating load contribute to the genomic landscape of divergence in two tree species connected by episodic gene flow

Speciation often involves repeated episodes of genetic contact between divergent populations before reproductive isolation (RI) is complete. Whole-genome sequencing (WGS) holds great promise for unravelling the genomic bases of speciation. We have studied two ecologically divergent, hybridizing species of the 'model tree' genus Populus (poplars, aspens, cottonwoods), Populus alba and P. tremula, using >8.6 million single nucleotide polymorphisms (SNPs) from WGS of population pools. We used the genomic data to (i) scan these species' genomes for regions of elevated and reduced divergence, (ii) assess key aspects of their joint demographic history based on genomewide site frequency spectra (SFS) and (iii) infer the potential roles of adaptive and deleterious coding mutations in shaping the genomic landscape of divergence. We identified numerous small, unevenly distributed genome regions without fixed polymorphisms despite high overall genomic differentiation. The joint SFS was best explained by ancient and repeated gene flow and allowed pinpointing candidate interspecific migrant tracts. The direction of selection (DoS) differed between genes in putative migrant tracts and the remainder of the genome, thus indicating the potential roles of adaptive divergence and segregating deleterious mutations on the evolution and breakdown of RI. Genes affected by positive selection during divergence were enriched for several functionally interesting groups, including well-known candidate 'speciation genes' involved in plant innate immunity. Our results suggest that adaptive divergence affects RI in these hybridizing species mainly through intrinsic and demographic processes. Integrating genomic with molecular data holds great promise for revealing the effects of particular genetic pathways on speciation.

Data from: Genomic regions with a history of divergent selection affect fitness of hybrids between two butterfly species

Speciation is the process by which reproductively isolated lineages arise, and is one of the fundamental means by which the diversity of life increases. Whereas numerous studies have documented an association between ecological divergence and reproductive isolation, relatively little is known about the role of natural selection in genome divergence during the process of speciation. Here we use genome-wide DNA sequences and Bayesian models to test the hypothesis that loci under divergent selection between two butterfly species (Lycaeides idas and L. melissa) also affect fitness in an admixed population. Locus-specific measures of genetic differentiation between L. idas and L. melissa and genomic introgression in hybrids varied across the genome. The most differentiated genetic regions were characterized by elevated L. idas ancestry in the admixed population, which occurs in L. idas-like habitat, consistent with the hypothesis that local adaptation contributes to speciation. Moreover, locus-specific measures of genetic differentiation (a metric of divergent selection) were positively associated with extreme genomic introgression (a metric of hybrid fitness). Interestingly, concordance of differentiation and introgression was only partial. We discuss multiple, complementary explanations for this partial concordance.

opencc-zeroDec 2012View details →

Data from: Divergence is focused on few genomic regions early in speciation: incipient speciation of sunflower ecotypes

Early in speciation, when populations are undergoing the transition from local adaptation to incipient species, is when a number of transient, but potentially important, processes appear to be most easily detected. These include signatures of selective sweeps that can point to asymmetry in selection between habitats, divergence hitchhiking and associations of adaptive genes with environments. In a genomic comparison of ecotypes of the prairie sunflower, Helianthus petiolaris, occurring at Great Sand Dunes National Park (Colorado), we found that selective sweeps were mainly restricted to the dune ecotype and that there was variation across the genome in whether proximity to the non-dune population constrained or promoted divergence. The major regions of divergence were few and large between ecotypes, in contrast with an interspecific comparison between H. petiolaris and a sympatric congener, H. annuus. In general, the large regions of divergence observed in the ecotypic comparison swamped locus-specific associations with environmental variables. In both comparisons, regions of high divergence occurred in portions of the genetic map with high marker density, probably reflecting regions of low recombination. The difference in genomic distributions of highly divergent regions between ecotypic and interspecific comparisons highlights the value of studies spanning the spectrum of speciation in related taxa.

Data from: Shared and non-shared genomic divergence in parallel ecotypes of Littorina saxatilis at a local scale

Parallel speciation occurs when selection drives repeated, independent adaptive divergence that reduces gene flow between ecotypes. Classical examples show parallel speciation originating from shared genomic variation, but this does not seem to be the case in the rough periwinkle (Littorina saxatilis) that has evolved considerable phenotypic diversity across Europe, including several distinct ecotypes. Small 'wave' ecotype snails inhabit exposed rocks and experience strong wave action, while thick-shelled, 'crab' ecotype snails are larger and experience crab predation on less exposed shores. Crab and wave ecotypes appear to have arisen in parallel and recent evidence suggests only marginal sharing of molecular variation linked to evolution of similar ecotypes in different parts of Europe. However the extent of genomic sharing is expected to increase with gene flow and more recent common ancestry. To test this we used de novo RAD-sequencing to quantify the extent of shared genomic divergence associated with phenotypic similarities among ecotype-pairs on three close islands (< 10 km distance) connected by weak gene flow (Nm~0.03) and with recent common ancestry (< 10,000 years). After accounting for technical issues, including a large proportion of null alleles due to a large effective population size, we found ~8-28% of positive outliers were shared between two islands and ~2-9% were shared among all three islands. This low level of sharing suggests that parallel phenotypic divergence in this system is not matched by shared genomic divergence despite a high probability of gene flow and standing genetic variation.