The Atlantic herring is one of the most common fish in the world and has been a crucial food resource in northern Europe. One school of herring may comprise billions of fish, but previous studies had only revealed very few genetic differences in herring from different geographic regions. This was unexpected since Atlantic herring is one of the few marine species that can reproduce throughout the brackish Baltic Sea, which can be about a tenth as salty as the Atlantic Ocean.

This unexpected finding could be explained in at least two different ways. Firstly, perhaps Atlantic herring are flexible enough to adapt to very different environments (i.e. high or low salinity) without much genetic change. Secondly, the previous studies only looked at a handful of sites in the Atlantic herring’s genome and so it is possible that genetic differences at other genes control this fish’s adaptation instead.

Now, Martinez Barrio, Lamichhaney, Fan, Rafati et al. have sequenced entire genomes from groups of Atlantic herring and revealed hundreds of sites that are associated with adaptation to the Baltic Sea. The analysis also identified a number of genes that control when these fish reproduce by comparing herring that spawn in the autumn with those that spawn in spring. This is important because natural populations must carefully time when they reproduce to maximize the survival of their young.

These new findings provide compelling evidence that changes in protein-coding genes and stretches of DNA that regulate the expression of other genes both contribute to adaptation in herrings. The analysis also clearly shows that variants of genes that contribute to adaptation were likely to evolve over time by accumulating multiple sequence changes affecting the same gene. Furthermore, these gene variants essentially form a rich “tool-box” that underlies the Atlantic herring’s adaptation to its environment, and different subpopulations of herring were found to have their own optimal sets of gene variants. For instance, autumn-spawning herring and spring-spawning herring from the Baltic Sea both have gene variants that favor adaptation to low salinity. However, autumn-spawning Baltic herring also share gene variants that favor spawning in the autumn with autumn-spawning herring from the North Sea, but not with spring-spawning Baltic herring.

The next step will be to study how the 500 or so genes identified affect adaptation at the molecular level. This will likely involve experiments with other model fish such as zebrafish and sticklebacks. Finally, these new findings can be directly applied to monitor stocks of herring to make herring fisheries more sustainable.

The Atlantic herring (Clupea harengus) is a pelagic fish that occurs in huge schools, up to billions of individuals. The herring fishery has been crucial for food security and economic development in Northern Europe and currently ranks among the five largest fisheries in the world with nearly 2 million tons fish landed annually (FAO, 2014). The herring is one of few marine fishes that reproduce throughout the Baltic Sea where the salinity drops to 2–3‰ in the Bothnian Bay, compared with 35‰ in the Atlantic Ocean (Figure 1A). This ecological adaptation must be recent because the brackish Baltic Sea has only existed for 10,000 years following the last glaciation (Andrén et al., 2011). Fishery biologists have for more than a century recognized stocks of herring defined by spawning location, spawning time, morphological characters and life history parameters (Iles and Sinclair, 1982; McQuinn, 1997). Several decades of genetic studies based on limited numbers of genetic markers (allozymes, microsatellites or SNPs) have not been able to verify this divergence; extremely low levels of differentiation even between geographically distant populations as well as between spring- and autumn-spawning herring have been observed (Andersson et al., 1981; Ryman et al., 1984; Larsson et al., 2007; 2010, Limborg et al., 2012). It has been proposed that lack of precision in homing behaviour of herring causes sufficient gene flow between stocks to counteract genetic differentiation (McQuinn, 1997). However, in a recent study we constructed an exome assembly and used this in combination with whole genome sequencing of eight population samples and found more than 400,000 SNPs (Lamichhaney et al., 2012). We confirmed lack of differentiation at most loci, whereas a small percentage showed highly significant differentiation. Simulations demonstrated that the distribution of fixation index (F_ST)-values among herring populations deviated significantly from expectation for selectively neutral loci.

Figure 1 with 1 supplement see all

Download asset Open asset

Genetic studies of ecological adaptation in natural populations is challenging because genetic differentiation caused by natural selection is often confounded with genetic differences due to genetic drift caused by restricted effective population sizes. An ideal species for studying the genetic basis of ecological adaptation should comprise subpopulations of infinite size and exposed to different ecological conditions. In such a species there is minute genetic drift and genetic differentiation is caused by selection resulting in local adaptation. The herring is close to being such an ideal subject for studies of ecological adaptation due to the extremely low levels of genetic differentiation at most loci as documented in previous studies (Andersson et al., 1981; Ryman et al., 1984; Larsson et al., 2007; 2010; Limborg et al., 2012; Lamichhaney et al., 2012). This unique opportunity together with herring being such a valuable natural resource prompted us to generate a genome assembly and perform genome sequencing of populations adapted to different ecological conditions.

Here we present a high-quality genome assembly for the Atlantic herring, and results of whole genome sequencing of 20 population samples using pooled DNA. The results were verified by individual genotyping using a custom-made 70k SNP array. Our study addresses two fundamentally different types of adaptations; one example of niche expansion (adaptation to low salinity), and one example of sympatric balancing selection (variation in the timing of reproduction). The results provide a comprehensive list of hundreds of independent loci underlying ecological adaptation and shed light on the relative importance of coding and non-coding variation. The results have important implications for sustainable fishery management, and provide a road map for cost effective high-resolution characterization of genetic diversity in natural populations.

Results

Genome assembly and annotation

Clupeiformes represents an early diverging clade of the otomorpha (Near et al., 2012) (Figure 2A). The genome size for herring has been estimated at ~850 Mb (Hinegardner and Rosen, 1972; Ida et al., 1991; Ohno et al., 1969) with no recent whole genome duplications reported. We performed whole genome assembly based on short read sequencing of libraries ranging from 170 bp to 20 kb insert sizes (Supplementary file 1A). The 808 Mb assembly had a scaffold N50 of 1.84 Mb with 23,336 predicted coding gene models. It showed a high degree of completeness based on RNAseq alignments, core gene analyses and comparisons to other fish gene sets (Table 1, Supplementary files 2, 3A–D, Figure 2B, Figure 2—figure supplements 1–2). The GC content was 44%, and repetitive elements made up 31% of the assembly (Table 1). Alignments of synthetic long reads (SLRs; Illumina) failed to significantly improve the assembly due to coincidental gaps between the assembly and the SLRs, but proved useful in phasing parental alleles (Materials and methods; Figure 2—figure supplements 3–4) and dramatically improved the discovery of indels larger than 30 bp compared to short Illumina reads (Supplementary file 1F). We identified 150 endogenous retroviruses (ERVs) constituting ~0.14% of the genomic sequence but none included open reading frames in all gag, pol and env genes (Supplementary file 1, Figure 2—figure supplement 5).

Figure 2 with 5 supplements see all

Download asset Open asset

Table 1

Species	Herring (Clupea harengus)	Zebrafish (Danio rerio)	Cod (Gadus morhua)	Coelacanth (Latimeria chalumnae)	Stickleback (Gasteosteus aculeatus)
Estimated genome size (Mb)	850	1,454a	830b	3,530c	530d
Assembly size (Mb)	808	1,412	753b	2,861e	463f
Contig N50 (kb)	21.3	25.0	2.8	12.7	83.2
Scaffold N50 (Mb)	1.84	1.55	0.69	0.92	10.8
Sequencing technologyg	I	S+I	R+I	I	S
Repeat content	30.9	52.2	25.4	27.7	25.2
%GC content	44.1	36.7	45.4	43.0	44.6
Heterozygosity	1/309	n.a.	1/500	1/435	1/700
Protein-coding gene count	23,336	26,459	22,154	19,033	20,787

1. ^a(Freeman et al., 2007; Vinogradov, 1998; Howe et al., 2013)

2. ^b(Star et al., 2011)

3. ^cGenome size calculated as pg x 0.978 × 10⁹ bp/pg; picogram values taken from Cimino and Bahr (1974)

4. ^d(Vinogradov, 1998; Jones et al., 2012)

5. ^e(Amemiya et al., 2013)

6. ^f(Jones et al., 2012)

7. ^gI=Illumina sequencing; S=Sanger sequencing; R=Roche 454 n.a.=not available

Population genetics and demographic history

Whole genome pooled sequencing was done using 20 population samples of herring from the Baltic Sea, Skagerrak, Kattegat, North Sea, Atlantic Ocean and Pacific Ocean (Figure 1A; Table 2); the latter sample represents the closely related Pacific herring (Clupea pallasii). Each pool comprised 47–100 fish and was sequenced to ~30x coverage. Furthermore, 16 fish, eight Baltic and eight Atlantic herring (Table 2), were sequenced individually to ~10x coverage. All data were aligned to the reference assembly and SNPs were called after rigorous quality filtering. We found 8.83 million SNPs when Pacific herring was included and 6.04 million among Atlantic and Baltic herring.

Table 2

Localitya	Sample	n	Position		Salinity (‰)	Date (yy/mm/dd)	Spawning season
Baltic Sea
Gulf of Bothnia (Kalix)b	BK	47	N 65°52’	E 22°43’	3	800629	spring
Bothnian Sea (Hudiksvall)	BU	100	N 61°45’	E 17°30’	6	120419	spring
Bothnian Sea (Gävle)	BÄV	100	N 60°43’	E 17°18’	6	120507	spring
Bothnian Sea (Gävle)	BÄS	100	N 60°43’	E 17°18’	6	120718	summer
Bothnian Sea (Gävle)	BÄH	100	N 60°44’	E 17°35’	6	120904	autumn
Bothnian Sea (Hästskär)c	BH	50	N 60°35’	E 17°48’	6	130522	spring
Central Baltic Sea (Vaxholm)b	BV	50	N 59°26’	E 18°18’	6	790827	spring
Central Baltic Sea (Gamleby)b	BG	49	N 57°50’	E 16°27’	7	790820	spring
Central Baltic Sea (Kalmar)	BR	100	N 57°39’	E 17°07’	7	120509	spring
Central Baltic Sea (Karlskrona)	BA	100	N 56°10’	E 15°33’	7	120530	spring
Central Baltic Sea	BC	100	N 55°24’	E 15°51’	8	111018	unknown
Southern Baltic Sea (Fehmarn)b	BF	50	N 54°50’	E 11°30’	12	790923	autumn
Kattegat, Skagerrak, North Sea, Atlantic Ocean
Kattegat (Träslövsläge)b	KT	50	N 57°03’	E 12°11’	20	781023	unknown
Kattegat (Björköfjorden)	KB	100	N 57°43’	E 11°42’	23	120312	spring
Skagerrak (Brofjorden)	SB	100	N 58°19’	E 11°21’	25	120320	spring
Skagerrak (Hamburgsund)b	SH	49	N 58°30’	E 11°13’	25	790319	spring
North Seab	NS	49	N 58°06’	E 06°10’	35	790805	autumn
Atlantic Ocean (Bergen)b	AB1	49	N 64°52’	E 10°15’	35	800207	spring
Atlantic Ocean (Bergen)c	AB2	8	N 60°35’	E 05°00’	33	130522	spring
Atlantic Ocean (Höfn)	AI	100	N 65°49’	W 12°58’	35	110915	spring
Pacific Ocean
Strait of Georgia (Vancouver)	PH	50	-	-	35	121124	-

1. ^aPlaces where the sample was landed (if known) are given in parenthesis

2. ^bSamples from previous study (Lamichhaney et al., 2012)

3. ^cEight Baltic herring from the BH sample and eight Atlantic herring from the AB2 sample were used for individual sequencing n=number of fish

Average nucleotide diversity was estimated by counting the frequency of heterozygous sites in the reference individual after stringent filtering for sequence quality and coverage (within one standard deviation of mean coverage). The estimate was one heterozygous site per 309 bp, giving a nucleotide diversity of 0.32%; no estimate based on the 16 herring sequenced individually deviated significantly from this value and there was no significant difference between Atlantic and Baltic herring. The average decay of linkage disequilibrium between loci was very steep, with average r² falling to 0.1 at a distance of 100 base pairs (Figure 1—figure supplement 1A).

The allele frequency distribution deviated significantly from the one expected for selectively neutral alleles at genetic equilibrium (p<2x10^-16, Kolmogorov-Smirnov test), due to an excess of rare alleles (Figure 1—figure supplement 1B) consistent with population expansion. The result is supported by the genome-wide distribution of Tajima’s D, which shows a global shift towards negative values (mean=−0.57 ± 0.01; Figure 1—figure supplement 1C). A demographic analysis using the diCal software (Sheehan et al., 2013) confirmed that herring have experienced an expansion in effective population size, roughly five- to ten-fold, and that the current N_e is on the order of 10⁶ individuals (Figure 1B); the results for Baltic and Atlantic herring were essentially identical. The result indicates that the effective population size minimum occurred at around one to two MYA, after the onset of the Quaternary ice age.

Phylogeny

The neighbor-joining phylogenetic tree including Atlantic, Baltic and Pacific herring shows a large phylogenetic distance between Pacific and Atlantic herring, as compared with the tiny genetic divergence among samples of Atlantic and Baltic herring (Figure 1C). We estimated the split between Atlantic and Pacific herring to ~2.2 million years ago based on mtDNA cytochrome B sequence divergence. The phylogenetic tree is consistent with minute differentiation at selectively neutral loci in Atlantic herring (Ryman et al., 1984; Lamichhaney et al., 2012); all subpopulations in the Eastern North Atlantic may have expanded from a common ancestral population after the last glaciation as indicated by demographic analysis (Figure 1B).

A closer examination of the tight cluster of Atlantic and Baltic herring populations reveals some structure consistent with geographic origin (Figure 1C). Samples from the Baltic Sea cluster on one half while samples from marine waters cluster on the other half of the tree. Only three populations are located at intermediate positions. Two of these are autumn-spawners from the Baltic Sea (BÄH and BF), indicating that autumn-spawning herring are genetically distinct from spring- and summer-spawning herring. The third sample (KT) at an intermediate position was sampled outside the spawning season and at the border between Kattegat and Baltic Sea and may represent a mixed sample of local Kattegat population and fish that spawn in the Baltic Sea but migrate into Kattegat for feeding.

Genetic adaptation to a new niche environment

The Atlantic (Clupea harengus harengus) and Baltic herring (Clupea harengus membras) were classified as subspecies by Linnaeus (1761) in the 18^th century. They are adapted to strikingly different environments, in particular regarding salinity that ranges from 2–3‰ in the Gulf of Bothnia to 12‰ in Southern Baltic Sea, whereas salinity in Kattegat, Skagerrak, North Sea and Atlantic Ocean is in the range 20‰–35‰ (Figure 1A; Table 2). To reveal loci underlying genetic adaptation associated with the recent niche expansion into brackish waters after the last glaciation we compared allele frequencies, SNP by SNP, in two superpools: one Atlantic including all populations from Atlantic Ocean, Skagerrak and Kattegat and a pool comprising all samples collected in Baltic Sea; this is justified by low differentiation at neutral loci as documented by the low F_ST-values when comparing all samples of Atlantic and Baltic herring (Figure 1D). Samples of autumn-spawning herring, a possible confounding factor, were excluded from the analysis. We used a stringent significance threshold of p<1x10^-10 (Bonferroni correction, p=8.2x10^-9).

We identified 46,045 SNPs that showed an allele frequency difference with p<1x10^-10 in the χ² test (Figure 3A; Supplementary file 3A). An important question is how many independent loci these represent. A conservative estimate of 472 independent loci was obtained (i) by only using SNPs with p<1x10^-20, (ii) by taking into account gaps in the assembly and (iii) by using the Comb-P software (Pedersen et al., 2012) to combine strongly correlated SNPs from the same genomic region (see Materials and methods). Figure 3A (lower panel) illustrates one of the most striking associations. For a large part of scaffold 218 there are no significant differences among Atlantic and Baltic samples whereas there are striking allele frequency differences over a 119.4 kb region; this is a characteristic pattern for differentiated regions, indicating that genetic adaptation typically occur as large haplotype blocks, often including multiple genes. A phylogenetic tree based on SNPs showing genetic differentiation between Atlantic and Baltic (Figure 3B) differs profoundly from the tree based on all SNPs (Figure 1C). With the exception of the two autumn-spawning populations BF and BÄH from the Baltic Sea, the position of all other populations match the variation in salinity perfectly with the population samples from the North Sea and Atlantic Ocean (35‰) at one end of the tree and samples from the brackish Baltic Sea (3‰–12‰) at the other end and with samples from Skagerrak (25‰) and Kattegat (20‰) at intermediate positions. The low genetic differentiation among Baltic samples, excluding the two autumn-spawning populations BF and BÄH, suggests that adaptation to brackish waters is a derived state.

Figure 3 with 2 supplements see all

Download asset Open asset

Figure 3C (upper panel) shows estimated allele frequencies for highly differentiated SNPs from five genomic regions in six population samples, each region showing an underlying genetic architecture with large and distinctly defined haplotype blocks. The Atlantic Ocean and North Sea samples are both nearly fixed for the reference allele at these SNPs. In contrast, the samples of Baltic herring were close to fixation for the alternate alleles. Interestingly, the sample (SB) collected in Skagerrak (salinity ~25‰) is most similar to the Atlantic Ocean and North Sea samples, but consistently shows a trend towards more intermediate allele frequencies at these loci.

We developed a 70k custom SNP chip to study differentiated regions in more detail and to use data from individual fish to confirm associations detected by pooled sequencing. The chip included 13,355 neutral SNPs evenly distributed across the genome and 59,205 SNPs showing genetic differentiation between subpopulations. Thirty fish each from 12 populations were used in the SNP screen. There was an excellent correlation between allele frequencies estimated with pooled sequencing and with the SNP chip (Figure 3—figure supplement 1). We constructed a phylogenetic tree (Figure 3C, lower panel) for haplotypes of highly differentiated SNPs from scaffold 218 present among individual fish from six representative populations, after phasing haplotypes using BEAGLE (Browning and Browning, 2007). As expected all fish from Atlantic Ocean and North Sea carried closely related “Atlantic” haplotypes. Two major haplotype groups were present among Baltic herring and with few exceptions Baltic herring carried only “Baltic” haplotypes. Fish from Skagerrak predominantly carried Atlantic haplotypes, but with a considerable proportion of Baltic haplotypes. Phylogenetic trees for other top scaffolds are presented in Figure 3—figure supplement 2.

There are many environmental and ecological differences between Atlantic Ocean and Baltic Sea e.g. temperature variability, eutrophication of the Baltic Sea, zooplankton and predator populations), but the most obvious difference concerns salinity. We used the Bayenv 2.0 (Günther and Coop, 2013) software to reveal which of the 472 independent loci detected with the χ² test showed the most consistent correlation with salinity. This analysis identified 3,335 SNPs from 122 independent regions with highly significant association to salinity (Supplementary file 3A). Twenty-one of the genes in these regions have previously been associated with hypertension in human and 36 of these genes showed differential expression in sticklebacks kept in freshwater or sea water (Supplementary file 3A).

Here we present three loci with striking association to salinity. Firstly, the 11 kb region in scaffold 899 (Figure 3C) contains a single gene, prolactin receptor (PRLR), that is essential for mammalian reproduction but has a central role for osmoregulation in fish (Manzon, 2002), and possibly in mammals (Schennink et al., 2015). Secondly, strong genetic differentiation was also observed at scaffold 346 (Figure 3A; p<1x10^-39). This signal overlaps HCE encoding high choriolytic enzyme. This locus was also identified as one of the most differentiated region in our screen for structural changes (Supplementary file 3B). A 4 kb region including part of the coding sequence showed a massive copy number amplification that had a strong negative correlation with salinity (Figure 3D). The outgroup, Pacific herring, showed an intermediate copy number. Interestingly, the Pacific herring spawns exclusively in shallow nearshore waters (Hay et al., 2009) often in estuaries and tidal zones where salinity varies, in contrast to deeper-spawning Atlantic herring. HCE is a protease, also denoted hatching enzyme, that solubilizes the inner layer of the egg envelope during hatching and adaptive evolution of this protein in relation to salinity has been reported (Kawaguchi et al., 2013). In herring, we found no coding changes implying altered transcriptional regulation. In fact, massive amplification of the promoter region is expected to alter gene expression. Hatching of the egg is probably a particularly challenging stage of development for a marine fish adapting to brackish conditions. Thirdly, a ~65 kb region downstream of solute carrier family 12 (sodium/chloride transporter) member 3 (SLC12A3) shows strong correlation with salinity (Figure 3E, Supplementary file 3A). SLC12A3, which has an established role in regulating osmotic balance, is associated with hypertension in human and shows differential expression in kidney tissue between sticklebacks kept in freshwater or sea water (Wang et al., 2014).

Genetic basis underlying timing of reproduction

Herring spawn from early spring to late fall. Prior to this study it was unknown if spawning time is entirely due to phenotypic plasticity, set by nutritional status and environmental conditions, or if genetic factors contribute (McQuinn, 1997). For example, it has been hypothesized that spawning time in the Baltic Sea is regulated by productivity of the system affecting maturation of fish prior to spawning (Aneer, 1985). To study this important question we collected spawning herring from the same geographic area, close to Gävle (Sweden), in May, July and September (Table 2). Our sampling included two other autumn-spawning populations collected in 1979, one from North Sea and the other from Southern Baltic Sea. We formed two superpools including three autumn-spawning and 10 spring-spawning population samples, respectively; the summer-spawners and one population of non-spawning herring (KT in Table 2) were excluded from the initial analysis. We identified 10,195 SNPs with significant allele frequency differences between pools (p<1x10^-10) and 69 regions with copy number variation (p<0.001) (Figure 4A); the highly differentiated SNPs represented at least 125 independent loci based on our strict criteria (see Materials and methods). The result demonstrates for the first time that autumn- and spring-spawning herring are genetically distinct and indicates that genetic factors affect spawning time. In a phylogenetic tree based on these 10,195 SNPs the autumn-spawning populations from the Baltic Sea and North Sea tended to cluster with spring-spawning herring from the Atlantic Ocean (Figure 4B).

Figure 4 with 2 supplements see all

Download asset Open asset

A general linear mixed model was used to identify which of the 125 independent loci showed the most consistent allele frequency differences between spring and autumn spawners. This analysis revealed 17 independent genomic regions that passed the stringent significance threshold of p<10^–10 (Bonferroni correction, p=4.9x10^-6) (Supplementary file 3C). We then illustrate the striking allele frequency differences at the four most significant regions using data from six different populations. As observed for the genetic adaptation to declined salinity (above), the most significant regions underlying seasonal reproductive timing typically consists of large haplotype blocks often containing multiple genes. Spring-spawning Atlantic and Baltic herring showed nearly identical allele frequencies at these loci while autumn-spawning herring from Baltic Sea and North Sea showed high frequencies of the alternate alleles (Figure 4C). Remarkably, summer-spawning herring showed a clear trend towards intermediate allele frequencies at all loci, most pronounced for scaffold 481 (Figure 4C). This may either reflect that this sample is an admixture of spring- and autumn-spawning herring or that it represents a distinct population. To explore this we investigated deviations from Hardy-Weinberg equilibrium using the F_IT statistics because we expect a heterozygote deficiency if this is a mix of two populations. The results, based on 1,500 SNPs all showing strong genetic differentiation between spring- and autumn-spawners and genotyped individually using the SNP chip, showed that the summer spawners (BÄS) did not deviate markedly from F_IT = 0 and in fact to a lesser extent than the spring-spawning population (BÄV) sampled at the same locality (Figure 4—figure supplement 1). For instance, individual genotyping of the highly differentiated SNPs from scaffold 481 (Figure 4C) resulted in mean F_IT = −0.10 (excess of heterozygotes) for the summer spawners (BÄS) whereas if the sample had constituted an equal mix of spring- and autumn spawners from the same locality (BÄV and BÄH) the expected F_IT-value would have been 0.46 (strong heterozygote deficiency). Thus, the data strongly suggest that these summer spawners represent a distinct population rather than admixture. Spawning time may be fine-tuned by the dosage of alleles affecting spawning time. The three populations from Gävle showed nearly identical allele frequencies at loci with strong genetic differentiation between Atlantic Ocean and Baltic Sea (Figure 3C), whereas they showed dramatic allele frequency differences at loci associated with spawning time (Figure 4C).

We used SNP-chip data to construct a haplotype tree based on highly differentiated SNPs in scaffold 190/1420. Two haplotype groups were strongly associated with spring- and autumn spawning (Figure 4D); haplotype trees for other top scaffolds are in Figure 4—figure supplement 2. The estimated average heterozygosity per polymorphic site across scaffold 1420 indicated a selective sweep among spring-spawning herring but not in autumn-spawning populations (Figure 4E). However, the nucleotide diversity did not show a significant difference between groups (spring: 0.24% ± 0.004%; autumn: 0.27% ± 0.003%). Thus, the number of variable sites are higher among spring-spawning herring, but the average heterozygosity per site is lower. One possible explanation for this observation is that a selective sweep happened at this locus in the past in spring-spawning herring, which was then followed by a population expansion allowing the accumulation of new mutations. This interpretation is supported by strong negative Tajima’s D-values in this region among spring-spawning Atlantic and Baltic herring (Figure 1—figure supplement 1E).

Genetic differences in spawning time are expected to involve photoperiodic regulation of reproduction. Interestingly, our strongest signals (p<1x10^-120) in this contrast is located within and up to 25 kb upstream of TSHR encoding thyroid-stimulating hormone receptor, which has a central role in this pathway in birds and mammals (Nakao et al., 2008; Ono et al., 2008; Hanon et al., 2008). Further, a second gene in the same scaffold (190/1420), calmodulin has a role in initiating reproduction following secretion of gonadotropin-releasing hormone (GnRH) (Melamed et al., 2012) downstream of TSHR signalling in photoperiodic regulation of reproduction. SOX11, one of the genes in the associated region in scaffold 1440 (Figure 4C), encodes a transcription factor that controls GnRH expression in GnRH-secreting neurons (Kim et al., 2011). Finally, ESR2a, in scaffold 312, encodes estrogen receptor beta that has a well established function in reproductive biology (Bondesson et al., 2015). Interestingly, a previous experimental study in sticklebacks also indicate that estrogen receptor signaling is involved in photoperiodic regulation of reproduction since treatment with aromatase inhibitors, which leads to an inhibition of the conversion of androgens to estrogens, altered photoperiodic regulation of male sexual maturation (Bornestaf et al., 1997). Also, the expression of ESR2 but not ESR1 is regulated by circadian factors in mice (Cai et al., 2008), consistent with our data suggesting that estrogen receptor beta (encoded by ESR2) is more important than estrogen receptor alpha (encoded by ESR1) for photoperiodic regulation of reproduction.

Adaptive haplotype blocks are maintained by selection

A common feature for the signatures of selection for adaptation to low salinity and for seasonal reproduction in herring is the presence of haplotype blocks (10–200 kb in size) showing strong differentiation (Figures 3C, 4C), despite the rapid decay of linkage disequilibrium at selectively neutral sites (Figure 1—figure supplement 1A). A possible explanation for the pattern is the presence of inversions suppressing recombination as previously shown in three-spined stickleback (Jones et al., 2012). We constructed 3.3 kb Nextera mate pair libraries for two Atlantic and two Baltic herring individuals to scan for inversions with a particular focus on regions under selection. However, few convincing inversion candidates were detected and none coincided with the regions highlighted in Figures 3C, 4C. Thus, inversions do not appear to be an important explanation for the presence of haplotype blocks.

Having excluded inversions as a major explanation for the long haplotype blocks, two other possible explanations were considered. Haplotype blocks may occur as a consequence of recent fast selective sweeps that leads to hitchhiking of neutral polymorphism in close genetic linkage with causal variants (Maynard-Smith and Haigh, 1974; Charlesworth et al., 1997). Alternatively, haplotype blocks involving multiple causal mutations may be maintained by natural selection. These two models give entirely different predictions as regards nucleotide diversity in the differentiated regions of the genome. The hitchhiking model predicts reduced levels of genetic diversity in the differentiated region whereas the haplotype evolution model implies that nucleotide diversity in the differentiated regions, even within populations, may be as high or even higher than in neutral regions because the haplotypes are expected to have been maintained during an evolutionary process. We decided to test this by comparing nucleotide diversity for the 30 most differentiated regions in the contrast Atlantic vs. Baltic within and between one population of Atlantic herring (Bergen) and one population of Baltic herring (Kalix). The nucleotide diversity turned out to be significantly higher in the differentiated regions than in random regions of the genome both within and between populations (Figure 5A). The same conclusion emerged from the analysis of the 30 most differentiated regions between autumn- and spring-spawning herring using the samples collected at the same locality (Gävle) in May and September (Figure 5B). Thus, we conclude that our data on genetic differentiation in herring is consistent with the evolution of haplotype blocks harbouring multiple causal variants. The model also implies that the presence of multiple alleles containing different combinations of causal variants is expected.

Figure 5

Download asset Open asset

Genomic distribution of causal variants

Genome-wide analysis combined with strong signatures of selection enabled us to explore the genomic distribution of sequence polymorphisms underlying ecological adaptation. We carried out an enrichment analysis as previously used to identify categories of SNPs showing differentiation between domestic and wild rabbits (Carneiro et al., 2014). We calculated the absolute allele frequency difference (dAF) for different categories of SNPs in the two contrasts Atlantic vs. Baltic and spring- vs. autumn spawning herring and sorted these into bins (dAF 0–0.05, etc.) for different categories of SNPs. In both contrasts the great majority of SNPs (>90%) showed a dAF lower than 0.10 (Figure 6, Supplementary file 3E).

Figure 6

Download asset Open asset

Non-synonymous substitutions showed the most striking enrichment in both contrasts and showed a steady increase above dAF=0.15 reaching a two-fold enrichment at dAF>0.50 (Figure 6, Supplementary file 3E). This enrichment must reflect natural selection acting on the protein sequence because synonymous substitutions did not show a similar strong enrichment at high dAF. All non-synonymous substitutions showing dAF>0.50 in any of the two contrasts are compiled in Supplementary file 3F. A striking feature of this list is the common occurrence of multiple high dAF SNPs in the same gene. The 74 non-synonymous changes with dAF>0.50 in the contrast Atlantic vs. Baltic occur in only 29 different genes and the corresponding figure for the contrast spring- vs. autumn-spawning is 21 non-synonymous changes in 9 genes. We excluded the possibility that the presence of multiple non-synonymous changes in many of the genes was explained by errors in gene models (non-coding sequences annotated as exons) by a comparative analysis with other teleosts. We identified the orthologous position for about two thirds of the positions listed in Supplementary file 3F, the great majority of these (58/62) were annotated as coding sequence also in other species (Supplementary file 3F).

SNPs located in the 5’untranslated and 3’untranslated regions (UTRs) showed a more consistent enrichment compared to synonymous changes implying that this enrichment is unlikely to be caused entirely by close linkage to coding sequences under selection. Thus, changes in UTRs have contributed to ecological adaptation in the herring, most likely due to their role in regulating mRNA stability and translation efficiency. In this analysis we combined 5’UTR and 3’UTR SNPs to avoid too small classes for the extremely high dAF. However, an analysis based on all SNPs showing a dAF > 0.1 in the Atlantic vs. Baltic contrast and all SNPs showing a dAF > 0.2 for the spring- vs. autumn-spawning contrast demonstrated that both 5’UTR and 3’UTR SNPs are overrepresented at high dAF and the trend is particularly strong for 5’UTR SNPs (Supplementary file 3G).

The importance of regulatory changes underlying ecological adaptation is evident from the highly significant enrichment of SNPs within 5 kb upstream and downstream of coding sequences (Figure 6, Supplementary file 3E). Further, the excess is particularly pronounced within 1 kb upstream of the coding sequence where the promoter is expected to be located (Supplementary file 3H). The enrichment is not as high as for non-synonymous changes but this does not mean that regulatory changes are less important than coding changes because a much higher proportion of SNPs within the 5 kb region flanking coding sequences are expected to be selectively neutral compared with those causing non-synonymous changes. Thus, it is possible that the enrichment of non-coding SNPs would be much higher if there was a better annotation of the functional significance of non-coding sequences in Atlantic herring.

Intergenic and intronic SNPs were in general underrepresented among SNPs showing high dAF (Figure 6). For the most differentiated SNPs (dAF > 0.50) the intergenic SNPs showed a marked underrepresentation in the Atlantic – Baltic contrast (M=-0.64; p=5.1 x 10^-25; Supplementary file 3E) while intronic SNPs were most underrepresented in the spring- vs. autumn-spawning contrast (M=-0.55; p=6.7 x 10^-7; Supplementary file 3E).

We also explored the possibility that loss of function-mutations have contributed to ecological adaptation. We identified a total of 469 nonsense mutations but expect that many of these will be false predictions due to errors in the gene model. Eight predicted nonsense mutations had a dAF higher than 0.20 in one of the contrasts and were further examined. Seven of these were unlikely to be correct annotations since the positions were not annotated as coding in zebrafish, and the remaining one had a dAF of 0.21 but was far from statistical significance. Thus, we conclude that gene inactivation is not a common mechanism for ecological adaptation.

1. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Herring genome

   Publicly available at NCBI Assembly (accession no. GCA_000966335.1).

   http://www.ncbi.nlm.nih.gov/assembly/318681
2. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Population sequencing

   Publicly available at the NCBI Sequence Read Archive (accession no. SRP056617).

   http://www.ncbi.nlm.nih.gov/sra/?term=SRP056617
3. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Data from: The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing

   Available at Dryad Digital Repository under a CC0 Public Domain Dedication.

   http://dx.doi.org/10.5061/dryad.5r774
4. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Population sequencing

   Publicly available at the NCBI Sequence Read Archive (accession no. SRP017094).

   http://www.ncbi.nlm.nih.gov/sra/?term=SRP017094
5. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Population sequencing

   Publicly available at the NCBI Sequence Read Archive (accession no. SRP017095).

   http://www.ncbi.nlm.nih.gov/sra/?term=SRP017095
6. 1. Martinez Barrio A
   2. Lamichhaney S
   3. Fan G
   4. Rafati N
   5. Pettersson M
   6. Zhang H
   7. Dainat J
   8. Ekman D
   9. Höppner M
   10. Jern P
   11. Martin M
   12. Nystedt B
   13. Liu X
   14. Chen W
   15. Liang X
   16. Shi C
   17. Fu Y
   18. Ma K
   19. Zhan X
   20. Feng C
   21. Gustafson U
   22. Rubin C
   23. Sällman Almén M
   24. Blass M
   25. Casini M
   26. Folkvord A
   27. Laikre L
   28. Ryman N
   29. Ming-Yuen Lee S
   30. Xu X
   31. Andersson L

   (2016) Population sequencing

   Publicly available at the NCBI Sequence Read Archive (accession no. SRP056617).

   http://www.ncbi.nlm.nih.gov/sra/?term=SRP056617

1. 1. Amemiya CT
   2. Alföldi J
   3. Lee AP
   4. Fan S
   5. Philippe H
   6. Maccallum I
   7. Braasch I
   8. Manousaki T
   9. Schneider I
   10. Rohner N
   11. Organ C
   12. Chalopin D
   13. Smith JJ
   14. Robinson M
   15. Dorrington RA
   16. Gerdol M
   17. Aken B
   18. Biscotti MA
   19. Barucca M
   20. Baurain D
   21. Berlin AM
   22. Blatch GL
   23. Buonocore F
   24. Burmester T
   25. Campbell MS
   26. Canapa A
   27. Cannon JP
   28. Christoffels A
   29. De Moro G
   30. Edkins AL
   31. Fan L
   32. Fausto AM
   33. Feiner N
   34. Forconi M
   35. Gamieldien J
   36. Gnerre S
   37. Gnirke A
   38. Goldstone JV
   39. Haerty W
   40. Hahn ME
   41. Hesse U
   42. Hoffmann S
   43. Johnson J
   44. Karchner SI
   45. Kuraku S
   46. Lara M
   47. Levin JZ
   48. Litman GW
   49. Mauceli E
   50. Miyake T
   51. Mueller MG
   52. Nelson DR
   53. Nitsche A
   54. Olmo E
   55. Ota T
   56. Pallavicini A
   57. Panji S
   58. Picone B
   59. Ponting CP
   60. Prohaska SJ
   61. Przybylski D
   62. Saha NR
   63. Ravi V
   64. Ribeiro FJ
   65. Sauka-Spengler T
   66. Scapigliati G
   67. Searle SM
   68. Sharpe T
   69. Simakov O
   70. Stadler PF
   71. Stegeman JJ
   72. Sumiyama K
   73. Tabbaa D
   74. Tafer H
   75. Turner-Maier J
   76. van Heusden P
   77. White S
   78. Williams L
   79. Yandell M
   80. Brinkmann H
   81. Volff JN
   82. Tabin CJ
   83. Shubin N
   84. Schartl M
   85. Jaffe DB
   86. Postlethwait JH
   87. Venkatesh B
   88. Di Palma F
   89. Lander ES
   90. Meyer A
   91. Lindblad-Toh K

   (2013) The african coelacanth genome provides insights into tetrapod evolution

   Nature 496:311–316.

   https://doi.org/10.1038/nature12027

   • Google Scholar
2. 1. Andersson L
   2. Ryman N
   3. Rosenberg R
   4. Ståhl G

   (1981) Genetic variability in atlantic herring (Clupea harengus harengus): Description of protein loci and population data

   Hereditas 95:69–78.

   https://doi.org/10.1111/j.1601-5223.1981.tb01330.x

   • Google Scholar
3. 1. Andersson L

   (2013) Molecular consequences of animal breeding

   Current Opinion in Genetics & Development 23:295–301.

   https://doi.org/10.1016/j.gde.2013.02.014

   • Google Scholar
4. Book

   1. Andrén T
   2. Björck S
   3. Andrén E
   4. Conley D
   5. Zillén L
   6. Anjar J

   (2011) The Development of the Baltic Sea Basin During the Last 130 ka

   In: Harff J, Bjorck S, Hoth P, editors. The Baltic Sea Basin. Berlin, Heidelberg: Springer Berlin Heidelberg. pp. 75–97.

   https://doi.org/10.1007/978-3-642-17220-5_4

   • Google Scholar
5. 1. Aneer G

   (1985) Some speculations about the Baltic herring (Clupea harengus membras) in connection with the eutrophication of the Baltic Sea

   Canadian Journal of Fisheries and Aquatic Sciences 42:s83–s90.

   https://doi.org/10.1139/f85-264

   • Google Scholar
6. 1. Ashburner M
   2. Ball CA
   3. Blake JA
   4. Botstein D
   5. Butler H
   6. Cherry JM
   7. Davis AP
   8. Dolinski K
   9. Dwight SS
   10. Eppig JT
   11. Harris MA
   12. Hill DP
   13. Issel-Tarver L
   14. Kasarskis A
   15. Lewis S
   16. Matese JC
   17. Richardson JE
   18. Ringwald M
   19. Rubin GM
   20. Sherlock G

   (2000) Gene ontology: Tool for the unification of biology. the gene ontology consortium

   Nature Genetics 25:25–29.

   https://doi.org/10.1038/75556

   • Google Scholar
7. 1. Bates D
   2. Maechler M
   3. Bolker BM
   4. Walker S

   (2014)

   lme4: Linear mixed-effects models using eigen and S4 .

   R package version.

   • Google Scholar
8. 1. Bondesson M
   2. Hao R
   3. Lin CY
   4. Williams C
   5. Gustafsson JÅ

   (2015) Estrogen receptor signaling during vertebrate development

   Biochimica Et Biophysica Acta 1849:142–151.

   https://doi.org/10.1016/j.bbagrm.2014.06.005

   • Google Scholar
9. 1. Bornestaf C
   2. Antonopoulou E
   3. Mayer I
   4. Borg B

   (1997) Effects of aromatase inhibitors on reproduction in male three-spined sticklebacks, gasterosteus aculeatus, exposed to long and short photoperiods

   Fish Physiology and Biochemistry 16:419–423.

   https://doi.org/10.1023/A:1007776517447

   • Google Scholar
10. 1. Browning SR
    2. Browning BL

    (2007) Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering

    American Journal of Human Genetics 81:1084–1097.

    https://doi.org/10.1086/521987

    • Google Scholar
11. 1. Burridge CP
    2. Craw D
    3. Fletcher D
    4. Waters JM

    (2008) Geological dates and molecular rates: Fish DNA sheds light on time dependency

    Molecular Biology and Evolution 25:624–633.

    https://doi.org/10.1093/molbev/msm271

    • Google Scholar
12. 1. Cai W
    2. Rambaud J
    3. Teboul M
    4. Masse I
    5. Benoit G
    6. Gustafsson JA
    7. Delaunay F
    8. Laudet V
    9. Pongratz I

    (2008) Expression levels of estrogen receptor beta are modulated by components of the molecular clock

    Molecular and Cellular Biology 28:784–793.

    https://doi.org/10.1128/MCB.00233-07

    • Google Scholar
13. 1. Cantarel BL
    2. Korf I
    3. Robb SM
    4. Parra G
    5. Ross E
    6. Moore B
    7. Holt C
    8. Sánchez Alvarado A
    9. Yandell M

    (2008) MAKER: An easy-to-use annotation pipeline designed for emerging model organism genomes

    Genome Research 18:188–196.

    https://doi.org/10.1101/gr.6743907

    • Google Scholar
14. 1. Carneiro M
    2. Rubin CJ
    3. Di Palma F
    4. Albert FW
    5. Alföldi J
    6. Barrio AM
    7. Pielberg G
    8. Rafati N
    9. Sayyab S
    10. Turner-Maier J
    11. Younis S
    12. Afonso S
    13. Aken B
    14. Alves JM
    15. Barrell D
    16. Bolet G
    17. Boucher S
    18. Burbano HA
    19. Campos R
    20. Chang JL
    21. Duranthon V
    22. Fontanesi L
    23. Garreau H
    24. Heiman D
    25. Johnson J
    26. Mage RG
    27. Peng Z
    28. Queney G
    29. Rogel-Gaillard C
    30. Ruffier M
    31. Searle S
    32. Villafuerte R
    33. Xiong A
    34. Young S
    35. Forsberg-Nilsson K
    36. Good JM
    37. Lander ES
    38. Ferrand N
    39. Lindblad-Toh K
    40. Andersson L

    (2014) Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication

    Science 345:1074–1079.

    https://doi.org/10.1126/science.1253714

    • Google Scholar
15. 1. Charlesworth B
    2. Nordborg M
    3. Charlesworth D

    (1997) The effects of local selection, balanced polymorphism and background selection on equilibrium patterns of genetic diversity in subdivided populations

    Genetical Research 70:155–174.

    https://doi.org/10.1017/S0016672397002954

    • Google Scholar
16. 1. Chen K
    2. Wallis JW
    3. McLellan MD
    4. Larson DE
    5. Kalicki JM
    6. Pohl CS
    7. McGrath SD
    8. Wendl MC
    9. Zhang Q
    10. Locke DP
    11. Shi X
    12. Fulton RS
    13. Ley TJ
    14. Wilson RK
    15. Ding L
    16. Mardis ER

    (2009) BreakDancer: An algorithm for high-resolution mapping of genomic structural variation

    Nature Methods 6:677–681.

    https://doi.org/10.1038/nmeth.1363

    • Google Scholar
17. 1. Cimino MC
    2. Bahr GF

    (1974) The nuclear DNA content and chromatin ultrastructure of the coelacanth Latimeria chalumnae

    Experimental Cell Research 88:263–272.

    https://doi.org/10.1016/0014-4827(74)90240-7

    • Google Scholar
18. 1. Cingolani P
    2. Platts A
    3. Wang le L
    4. Coon M
    5. Nguyen T
    6. Wang L
    7. Land SJ
    8. Lu X
    9. Ruden DM

    (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: Snps in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

    Fly 6:80–92.

    https://doi.org/10.4161/fly.19695

    • Google Scholar
19. 1. Danecek P
    2. Auton A
    3. Abecasis G
    4. Albers CA
    5. Banks E
    6. DePristo MA
    7. Handsaker RE
    8. Lunter G
    9. Marth GT
    10. Sherry ST
    11. McVean G
    12. Durbin R
    13. 1000 Genomes Project Analysis Group

    (2011) The variant call format and vcftools

    Bioinformatics 27:2156–2158.

    https://doi.org/10.1093/bioinformatics/btr330

    • Google Scholar
20. 1. Dickey-Collas M
    2. Nash RDM
    3. Brunel T
    4. van Damme CJG
    5. Marshall CT
    6. Payne MR
    7. Corten A
    8. Geffen AJ
    9. Peck MA
    10. Hatfield EMC
    11. Hintzen NT
    12. Enberg K
    13. Kell LT
    14. Simmonds EJ

    (2010) Lessons learned from stock collapse and recovery of North Sea herring: A review

    ICES Journal of Marine Science 67:1875–1886.

    https://doi.org/10.1093/icesjms/fsq033

    • Google Scholar
21. 1. Doležel J
    2. Bartoš J
    3. Voglmayr H
    4. Greilhuber J

    (2003)

    Cytometry Part A,51A

    127–128, Letter to the editor, Cytometry Part A,51A.

    • Google Scholar
22. 1. Edwards M
    2. Richardson AJ

    (2004) Impact of climate change on marine pelagic phenology and trophic mismatch

    Nature 430:881–884.

    https://doi.org/10.1038/nature02808

    • Google Scholar
23. 1. Ewing G
    2. Hermisson J

    (2010) Msms: A coalescent simulation program including recombination, demographic structure and selection at a single locus

    Bioinformatics 26:2064–2065.

    https://doi.org/10.1093/bioinformatics/btq322

    • Google Scholar
24. Book

    1. FAO

    (2014)

    Yearbook Fishery and Aquaculture Statistics

    Rome: Fao.

    • Google Scholar
25. 1. Felsenstein J

    (1989)

    Phylip - phylogeny inference package (version 3.2)

    Cladistics 5:164–166.

    • Google Scholar
26. 1. Finn RD
    2. Bateman A
    3. Clements J
    4. Coggill P
    5. Eberhardt RY
    6. Eddy SR
    7. Heger A
    8. Hetherington K
    9. Holm L
    10. Mistry J
    11. Sonnhammer EL
    12. Tate J
    13. Punta M

    (2014) Pfam: The protein families database

    Nucleic Acids Research 42:D222–D230.

    https://doi.org/10.1093/nar/gkt1223

    • Google Scholar
27. 1. Freeman JL
    2. Adeniyi A
    3. Banerjee R
    4. Dallaire S
    5. Maguire SF
    6. Chi J
    7. Ng BL
    8. Zepeda C
    9. Scott CE
    10. Humphray S
    11. Rogers J
    12. Zhou Y
    13. Zon LI
    14. Carter NP
    15. Yang F
    16. Lee C

    (2007) Definition of the zebrafish genome using flow cytometry and cytogenetic mapping

    BMC Genomics 8:.

    https://doi.org/10.1186/1471-2164-8-195

    • Google Scholar
28. 1. Glasauer SM
    2. Neuhauss SC

    (2014) Whole-genome duplication in teleost fishes and its evolutionary consequences

    Molecular Genetics and Genomics : MGG 289:1045–1060.

    https://doi.org/10.1007/s00438-014-0889-2

    • Google Scholar
29. 1. Grabherr MG
    2. Haas BJ
    3. Yassour M
    4. Levin JZ
    5. Thompson DA
    6. Amit I
    7. Adiconis X
    8. Fan L
    9. Raychowdhury R
    10. Zeng Q
    11. Chen Z
    12. Mauceli E
    13. Hacohen N
    14. Gnirke A
    15. Rhind N
    16. di Palma F
    17. Birren BW
    18. Nusbaum C
    19. Lindblad-Toh K
    20. Friedman N
    21. Regev A

    (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome

    Nature Biotechnology 29:644–652.

    https://doi.org/10.1038/nbt.1883

    • Google Scholar
30. 1. Günther T
    2. Coop G

    (2013) Robust identification of local adaptation from allele frequencies

    Genetics 195:205–220.

    https://doi.org/10.1534/genetics.113.152462

    • Google Scholar
31. 1. Haas BJ
    2. Delcher AL
    3. Mount SM
    4. Wortman JR
    5. Smith RK
    6. Hannick LI
    7. Maiti R
    8. Ronning CM
    9. Rusch DB
    10. Town CD
    11. Salzberg SL
    12. White O

    (2003) Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies

    Nucleic Acids Research 31:5654–5666.

    https://doi.org/10.1093/nar/gkg770

    • Google Scholar
32. Website

    1. Hall B
    2. Derego T
    3. Geib S

    (2014) Gag: The genome annotation generator [online]. available

    http://genomeannotation.github.io/Gag
33. 1. Hanon EA
    2. Lincoln GA
    3. Fustin JM
    4. Dardente H
    5. Masson-Pévet M
    6. Morgan PJ
    7. Hazlerigg DG

    (2008) Ancestral TSH mechanism signals summer in a photoperiodic mammal

    Current Biology : CB 18:1147–1152.

    https://doi.org/10.1016/j.cub.2008.06.076

    • Google Scholar
34. 1. Hay DE
    2. McCarter PB
    3. Daniel KS
    4. Schweigert JF

    (2009) Spatial diversity of Pacific herring (Clupea pallasi) spawning areas

    ICES Journal of Marine Science 66:1662–1666.

    https://doi.org/10.1093/icesjms/fsp139

    • Google Scholar
35. 1. Hayward A
    2. Cornwallis CK
    3. Jern P

    (2015) Pan-vertebrate comparative genomics unmasks retrovirus macroevolution

    Proceedings of the National Academy of Sciences of the United States of America 112:464–469.

    https://doi.org/10.1073/pnas.1414980112

    • Google Scholar
36. 1. Hinegardner R
    2. Rosen DE

    (1972) Cellular DNA content and the evolution of teleostean fishes

    American Naturalist 106:621–644.

    https://doi.org/10.1086/282801

    • Google Scholar
37. 1. Holt C
    2. Yandell M

    (2011) Maker2: An annotation pipeline and genome-database management tool for second-generation genome projects

    BMC Bioinformatics 12:.

    https://doi.org/10.1186/1471-2105-12-491

    • Google Scholar
38. 1. Howe K
    2. Clark MD
    3. Torroja CF
    4. Torrance J
    5. Berthelot C
    6. Muffato M
    7. Collins JE
    8. Humphray S
    9. McLaren K
    10. Matthews L
    11. McLaren S
    12. Sealy I
    13. Caccamo M
    14. Churcher C
    15. Scott C
    16. Barrett JC
    17. Koch R
    18. Rauch GJ
    19. White S
    20. Chow W
    21. Kilian B
    22. Quintais LT
    23. Guerra-Assunção JA
    24. Zhou Y
    25. Gu Y
    26. Yen J
    27. Vogel JH
    28. Eyre T
    29. Redmond S
    30. Banerjee R
    31. Chi J
    32. Fu B
    33. Langley E
    34. Maguire SF
    35. Laird GK
    36. Lloyd D
    37. Kenyon E
    38. Donaldson S
    39. Sehra H
    40. Almeida-King J
    41. Loveland J
    42. Trevanion S
    43. Jones M
    44. Quail M
    45. Willey D
    46. Hunt A
    47. Burton J
    48. Sims S
    49. McLay K
    50. Plumb B
    51. Davis J
    52. Clee C
    53. Oliver K
    54. Clark R
    55. Riddle C
    56. Elliot D
    57. Eliott D
    58. Threadgold G
    59. Harden G
    60. Ware D
    61. Begum S
    62. Mortimore B
    63. Mortimer B
    64. Kerry G
    65. Heath P
    66. Phillimore B
    67. Tracey A
    68. Corby N
    69. Dunn M
    70. Johnson C
    71. Wood J
    72. Clark S
    73. Pelan S
    74. Griffiths G
    75. Smith M
    76. Glithero R
    77. Howden P
    78. Barker N
    79. Lloyd C
    80. Stevens C
    81. Harley J
    82. Holt K
    83. Panagiotidis G
    84. Lovell J
    85. Beasley H
    86. Henderson C
    87. Gordon D
    88. Auger K
    89. Wright D
    90. Collins J
    91. Raisen C
    92. Dyer L
    93. Leung K
    94. Robertson L
    95. Ambridge K
    96. Leongamornlert D
    97. McGuire S
    98. Gilderthorp R
    99. Griffiths C
    100. Manthravadi D
    101. Nichol S
    102. Barker G
    103. Whitehead S
    104. Kay M
    105. Brown J
    106. Murnane C
    107. Gray E
    108. Humphries M
    109. Sycamore N
    110. Barker D
    111. Saunders D
    112. Wallis J
    113. Babbage A
    114. Hammond S
    115. Mashreghi-Mohammadi M
    116. Barr L
    117. Martin S
    118. Wray P
    119. Ellington A
    120. Matthews N
    121. Ellwood M
    122. Woodmansey R
    123. Clark G
    124. Cooper J
    125. Cooper J
    126. Tromans A
    127. Grafham D
    128. Skuce C
    129. Pandian R
    130. Andrews R
    131. Harrison E
    132. Kimberley A
    133. Garnett J
    134. Fosker N
    135. Hall R
    136. Garner P
    137. Kelly D
    138. Bird C
    139. Palmer S
    140. Gehring I
    141. Berger A
    142. Dooley CM
    143. Ersan-Ürün Z
    144. Eser C
    145. Geiger H
    146. Geisler M
    147. Karotki L
    148. Kirn A
    149. Konantz J
    150. Konantz M
    151. Oberländer M
    152. Rudolph-Geiger S
    153. Teucke M
    154. Lanz C
    155. Raddatz G
    156. Osoegawa K
    157. Zhu B
    158. Rapp A
    159. Widaa S
    160. Langford C
    161. Yang F
    162. Schuster SC
    163. Carter NP
    164. Harrow J
    165. Ning Z
    166. Herrero J
    167. Searle SM
    168. Enright A
    169. Geisler R
    170. Plasterk RH
    171. Lee C
    172. Westerfield M
    173. de Jong PJ
    174. Zon LI
    175. Postlethwait JH
    176. Nüsslein-Volhard C
    177. Hubbard TJ
    178. Roest Crollius H
    179. Rogers J
    180. Stemple DL

    (2013) The zebrafish reference genome sequence and its relationship to the human genome

    Nature 496:498–503.

    https://doi.org/10.1038/nature12111

    • Google Scholar
39. 1. Hunter S
    2. Jones P
    3. Mitchell A
    4. Apweiler R
    5. Attwood TK
    6. Bateman A
    7. Bernard T
    8. Binns D
    9. Bork P
    10. Burge S
    11. de Castro E
    12. Coggill P
    13. Corbett M
    14. Das U
    15. Daugherty L
    16. Duquenne L
    17. Finn RD
    18. Fraser M
    19. Gough J
    20. Haft D
    21. Hulo N
    22. Kahn D
    23. Kelly E
    24. Letunic I
    25. Lonsdale D
    26. Lopez R
    27. Madera M
    28. Maslen J
    29. McAnulla C
    30. McDowall J
    31. McMenamin C
    32. Mi H
    33. Mutowo-Muellenet P
    34. Mulder N
    35. Natale D
    36. Orengo C
    37. Pesseat S
    38. Punta M
    39. Quinn AF
    40. Rivoire C
    41. Sangrador-Vegas A
    42. Selengut JD
    43. Sigrist CJ
    44. Scheremetjew M
    45. Tate J
    46. Thimmajanarthanan M
    47. Thomas PD
    48. Wu CH
    49. Yeats C
    50. Yong SY

    (2012) InterPro in 2011: New developments in the family and domain prediction database

    Nucleic Acids Research 40:D306–D312.

    https://doi.org/10.1093/nar/gkr948

    • Google Scholar
40. 1. ICES

    (2014)

    919, Report of the Baltic Fisheries Assessment Working Group (WGBFAS)., ICES HQ, Copenhagen, Denmark

    919, Report of the Baltic Fisheries Assessment Working Group (WGBFAS)., ICES HQ, Copenhagen, Denmark.

    • Google Scholar
41. 1. Ida H
    2. Oka N
    3. Hayashigaki K

    (1991)

    Karyotypes and cellular DNA contents of three species of the subfamily Clupeinae

    Japanese Journal of Ichthyology 38:289–294.

    • Google Scholar
42. 1. Iles TD
    2. Sinclair M

    (1982) Atlantic herring: Stock discreteness and abundance

    Science 215:627–633.

    https://doi.org/10.1126/science.215.4533.627

    • Google Scholar
43. 1. Jones FC
    2. Grabherr MG
    3. Chan YF
    4. Russell P
    5. Mauceli E
    6. Johnson J
    7. Swofford R
    8. Pirun M
    9. Zody MC
    10. White S
    11. Birney E
    12. Searle S
    13. Schmutz J
    14. Grimwood J
    15. Dickson MC
    16. Myers RM
    17. Miller CT
    18. Summers BR
    19. Knecht AK
    20. Brady SD
    21. Zhang H
    22. Pollen AA
    23. Howes T
    24. Amemiya C
    25. Baldwin J
    26. Bloom T
    27. Jaffe DB
    28. Nicol R
    29. Wilkinson J
    30. Lander ES
    31. Di Palma F
    32. Lindblad-Toh K
    33. Kingsley DM
    34. Broad Institute Genome Sequencing Platform & Whole Genome Assembly Team

    (2012) The genomic basis of adaptive evolution in threespine sticklebacks

    Nature 484:55–61.

    https://doi.org/10.1038/nature10944

    • Google Scholar
44. 1. Jones P
    2. Binns D
    3. Chang HY
    4. Fraser M
    5. Li W
    6. McAnulla C
    7. McWilliam H
    8. Maslen J
    9. Mitchell A
    10. Nuka G
    11. Pesseat S
    12. Quinn AF
    13. Sangrador-Vegas A
    14. Scheremetjew M
    15. Yong SY
    16. Lopez R
    17. Hunter S

    (2014) InterProScan 5: Genome-scale protein function classification

    Bioinformatics 30:1236–1240.

    https://doi.org/10.1093/bioinformatics/btu031

    • Google Scholar
45. 1. Karasov T
    2. Messer PW
    3. Petrov DA

    (2010) Evidence that adaptation in Drosophila is not limited by mutation at single sites

    PLoS Genetics 6:e1000924.

    https://doi.org/10.1371/journal.pgen.1000924

    • Google Scholar
46. 1. Kawaguchi M
    2. Yasumasu S
    3. Shimizu A
    4. Kudo N
    5. Sano K
    6. Iuchi I
    7. Nishida M

    (2013) Adaptive evolution of fish hatching enzyme: One amino acid substitution results in differential salt dependency of the enzyme

    Journal of Experimental Biology 216:1609–1615.

    https://doi.org/10.1242/jeb.069716

    • Google Scholar
47. 1. Kijas JW
    2. Lenstra JA
    3. Hayes B
    4. Boitard S
    5. Porto Neto LR
    6. San Cristobal M
    7. Servin B
    8. McCulloch R
    9. Whan V
    10. Gietzen K
    11. Paiva S
    12. Barendse W
    13. Ciani E
    14. Raadsma H
    15. McEwan J
    16. Dalrymple B
    17. International Sheep Genomics Consortium Members

    (2012) Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection

    PLoS Biology 10:e1001258.

    https://doi.org/10.1371/journal.pbio.1001258

    • Google Scholar
48. 1. Kim HD
    2. Choe HK
    3. Chung S
    4. Kim M
    5. Seong JY
    6. Son GH
    7. Kim K

    (2011) Class-C SOX transcription factors control gnrh gene expression via the intronic transcriptional enhancer

    Molecular Endocrinology 25:1184–1196.

    https://doi.org/10.1210/me.2010-0332

    • Google Scholar
49. 1. Kim D
    2. Pertea G
    3. Trapnell C
    4. Pimentel H
    5. Kelley R
    6. Salzberg SL

    (2013) TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions

    Genome Biology 14:R36.

    https://doi.org/10.1186/gb-2013-14-4-r36

    • Google Scholar
50. 1. King MC
    2. Wilson AC

    (1975) Evolution at two levels in humans and chimpanzees

    Science 188:107–116.

    https://doi.org/10.1126/science.1090005

    • Google Scholar
51. 1. Kuleshov V
    2. Xie D
    3. Chen R
    4. Pushkarev D
    5. Ma Z
    6. Blauwkamp T
    7. Kertesz M
    8. Snyder M

    (2014) Whole-genome haplotyping using long reads and statistical methods

    Nature Biotechnology 32:261–266.

    https://doi.org/10.1038/nbt.2833

    • Google Scholar
52. 1. Lamichhaney S
    2. Martinez Barrio A
    3. Rafati N
    4. Sundström G
    5. Rubin CJ
    6. Gilbert ER
    7. Berglund J
    8. Wetterbom A
    9. Laikre L
    10. Webster MT
    11. Grabherr M
    12. Ryman N
    13. Andersson L

    (2012) Population-scale sequencing reveals genetic differentiation due to local adaptation in Atlantic herring

    Proceedings of the National Academy of Sciences of the United States of America 109:19345–19350.

    https://doi.org/10.1073/pnas.1216128109

    • Google Scholar
53. 1. Lamichhaney S
    2. Berglund J
    3. Almén MS
    4. Maqbool K
    5. Grabherr M
    6. Martinez-Barrio A
    7. Promerová M
    8. Rubin CJ
    9. Wang C
    10. Zamani N
    11. Grant BR
    12. Grant PR
    13. Webster MT
    14. Andersson L

    (2015) Evolution of Darwin's finches and their beaks revealed by genome sequencing

    Nature 518:371–375.

    https://doi.org/10.1038/nature14181

    • Google Scholar
54. 1. Lander ES
    2. Waterman MS

    (1988) Genomic mapping by fingerprinting random clones: A mathematical analysis

    Genomics 2:231–239.

    https://doi.org/10.1016/0888-7543(88)90007-9

    • Google Scholar
55. 1. Larsson LC
    2. Laikre L
    3. Palm S
    4. André C
    5. Carvalho GR
    6. Ryman N

    (2007) Concordance of allozyme and microsatellite differentiation in a marine fish, but evidence of selection at a microsatellite locus

    Molecular Ecology 16:1135–1147.

    https://doi.org/10.1111/j.1365-294X.2006.03217.x

    • Google Scholar
56. 1. Larsson LC
    2. Laikre L
    3. André C
    4. Dahlgren TG
    5. Ryman N

    (2010) Temporally stable genetic structure of heavily exploited Atlantic herring (Clupea harengus) in swedish waters

    Heredity 104:40–51.

    https://doi.org/10.1038/hdy.2009.98

    • Google Scholar
57. 1. Li L
    2. Stoeckert CJ
    3. Roos DS

    (2003) OrthoMCL: Identification of ortholog groups for eukaryotic genomes

    Genome Research 13:2178–2189.

    https://doi.org/10.1101/gr.1224503

    • Google Scholar
58. 1. Li H
    2. Durbin R

    (2009) Fast and accurate short read alignment with burrows-wheeler transform

    Bioinformatics 25:1754–1760.

    https://doi.org/10.1093/bioinformatics/btp324

    • Google Scholar
59. 1. Li R
    2. Zhu H
    3. Ruan J
    4. Qian W
    5. Fang X
    6. Shi Z
    7. Li Y
    8. Li S
    9. Shan G
    10. Kristiansen K
    11. Li S
    12. Yang H
    13. Wang J
    14. Wang J

    (2010) De novo assembly of human genomes with massively parallel short read sequencing

    Genome Research 20:265–272.

    https://doi.org/10.1101/gr.097261.109

    • Google Scholar
60. 1. Limborg MT
    2. Helyar SJ
    3. De Bruyn M
    4. Taylor MI
    5. Nielsen EE
    6. Ogden R
    7. Carvalho GR
    8. Bekkevold D
    9. FPT Consortium

    (2012) Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus)

    Molecular Ecology 21:3686–3703.

    https://doi.org/10.1111/j.1365-294X.2012.05639.x

    • Google Scholar
61. 1. Linnaeus C

    (1761)

    Fauna Suecica

    Fauna Suecica, Stockholm.

    • Google Scholar
62. 1. Lowe TM
    2. Eddy SR

    (1997) tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence

    Nucleic Acids Research 25:955–964.

    https://doi.org/10.1093/nar/25.5.0955

    • Google Scholar
63. 1. Magrane M
    2. Consortium U

    (2011) Uniprot knowledgebase: A hub of integrated protein data

    Database 2011:.

    https://doi.org/10.1093/database/bar009

    • Google Scholar
64. 1. Manzon LA

    (2002) The role of prolactin in fish osmoregulation: A review

    General and Comparative Endocrinology 125:291–310.

    https://doi.org/10.1006/gcen.2001.7746

    • Google Scholar
65. 1. Marco-Sola S
    2. Sammeth M
    3. Guigó R
    4. Ribeca P

    (2012) The GEM mapper: Fast, accurate and versatile alignment by filtration

    Nature Methods 9:1185–1188.

    https://doi.org/10.1038/nmeth.2221

    • Google Scholar
66. 1. Smith JM
    2. Haigh J

    (1974) The hitch-hiking effect of a favourable gene

    Genetical Research 23:23–35.

    https://doi.org/10.1017/s0016672300014634

    • Google Scholar
67. 1. McCoy RC
    2. Taylor RW
    3. Blauwkamp TA
    4. Kelley JL
    5. Kertesz M
    6. Pushkarev D
    7. Petrov DA
    8. Fiston-Lavier AS

    (2014) Illumina truseq synthetic long-reads empower de novo assembly and resolve complex, highly-repetitive transposable elements

    PLoS One 9:e106689.

    https://doi.org/10.1371/journal.pone.0106689

    • Google Scholar
68. 1. McKenna A
    2. Hanna M
    3. Banks E
    4. Sivachenko A
    5. Cibulskis K
    6. Kernytsky A
    7. Garimella K
    8. Altshuler D
    9. Gabriel S
    10. Daly M
    11. DePristo MA

    (2010) The genome analysis toolkit: A mapreduce framework for analyzing next-generation DNA sequencing data

    Genome Research 20:1297–1303.

    https://doi.org/10.1101/gr.107524.110

    • Google Scholar
69. 1. McQuinn I

    (1997) Metapopulations and the Atlantic herring

    Reviews in Fish Biology and Fisheries 7:297–329.

    https://doi.org/10.1023/A:1018491828875

    • Google Scholar
70. 1. Melamed P
    2. Savulescu D
    3. Lim S
    4. Wijeweera A
    5. Luo Z
    6. Luo M
    7. Pnueli L

    (2012) Gonadotrophin-releasing hormone signalling downstream of calmodulin

    Journal of Neuroendocrinology 24:1463–1475.

    https://doi.org/10.1111/j.1365-2826.2012.02359.x

    • Google Scholar
71. 1. Meuwissen T
    2. Hayes B
    3. Goddard M

    (2013) Accelerating improvement of livestock with genomic selection

    Annual Review of Animal Biosciences 1:221–237.

    https://doi.org/10.1146/annurev-animal-031412-103705

    • Google Scholar
72. 1. Nakane Y
    2. Ikegami K
    3. Iigo M
    4. Ono H
    5. Takeda K
    6. Takahashi D
    7. Uesaka M
    8. Kimijima M
    9. Hashimoto R
    10. Arai N
    11. Suga T
    12. Kosuge K
    13. Abe T
    14. Maeda R
    15. Senga T
    16. Amiya N
    17. Azuma T
    18. Amano M
    19. Abe H
    20. Yamamoto N
    21. Yoshimura T

    (2013) The saccus vasculosus of fish is a sensor of seasonal changes in day length

    Nature Communications 4:.

    https://doi.org/10.1038/ncomms3108

    • Google Scholar
73. 1. Nakao N
    2. Ono H
    3. Yamamura T
    4. Anraku T
    5. Takagi T
    6. Higashi K
    7. Yasuo S
    8. Katou Y
    9. Kageyama S
    10. Uno Y
    11. Kasukawa T
    12. Iigo M
    13. Sharp PJ
    14. Iwasawa A
    15. Suzuki Y
    16. Sugano S
    17. Niimi T
    18. Mizutani M
    19. Namikawa T
    20. Ebihara S
    21. Ueda HR
    22. Yoshimura T

    (2008) Thyrotrophin in the pars tuberalis triggers photoperiodic response

    Nature 452:317–322.

    https://doi.org/10.1038/nature06738

    • Google Scholar
74. 1. Near TJ
    2. Eytan RI
    3. Dornburg A
    4. Kuhn KL
    5. Moore JA
    6. Davis MP
    7. Wainwright PC
    8. Friedman M
    9. Smith WL

    (2012) Resolution of ray-finned fish phylogeny and timing of diversification

    Proceedings of the National Academy of Sciences of the United States of America 109:13698–13703.

    https://doi.org/10.1073/pnas.1206625109

    • Google Scholar
75. 1. Ohno S
    2. Muramoto J
    3. Klein J
    4. Atkin NB

    (1969)

    Chromosomes Today

    Diploid-tetraploid relationship in clupeoid and salmonoid fish, Chromosomes Today, Edinburgh, Oliver & Boyd.

    • Google Scholar
76. 1. Ono H
    2. Hoshino Y
    3. Yasuo S
    4. Watanabe M
    5. Nakane Y
    6. Murai A
    7. Ebihara S
    8. Korf HW
    9. Yoshimura T

    (2008) Involvement of thyrotropin in photoperiodic signal transduction in mice

    Proceedings of the National Academy of Sciences of the United States of America 105:18238–18242.

    https://doi.org/10.1073/pnas.0808952105

    • Google Scholar
77. 1. Parra G
    2. Bradnam K
    3. Korf I

    (2007) Cegma: A pipeline to accurately annotate core genes in eukaryotic genomes

    Bioinformatics 23:1061–1067.

    https://doi.org/10.1093/bioinformatics/btm071

    • Google Scholar
78. 1. Parra G
    2. Bradnam K
    3. Ning Z
    4. Keane T
    5. Korf I

    (2009) Assessing the gene space in draft genomes

    Nucleic Acids Research 37:289–297.

    https://doi.org/10.1093/nar/gkn916

    • Google Scholar
79. 1. Pedersen BS
    2. Schwartz DA
    3. Yang IV
    4. Kechris KJ

    (2012) Comb-p: Software for combining, analyzing, grouping and correcting spatially correlated p-values

    Bioinformatics 28:2986–2988.

    https://doi.org/10.1093/bioinformatics/bts545

    • Google Scholar
80. 1. Price MN
    2. Dehal PS
    3. Arkin AP

    (2010) Fasttree 2--approximately maximum-likelihood trees for large alignments

    PLoS One 5:e9490.

    https://doi.org/10.1371/journal.pone.0009490

    • Google Scholar
81. 1. Purcell S
    2. Neale B
    3. Todd-Brown K
    4. Thomas L
    5. Ferreira MA
    6. Bender D
    7. Maller J
    8. Sklar P
    9. de Bakker PI
    10. Daly MJ
    11. Sham PC

    (2007) Plink: A tool set for whole-genome association and population-based linkage analyses

    American Journal of Human Genetics 81:559–575.

    https://doi.org/10.1086/519795

    • Google Scholar
82. 1. Rausch T
    2. Zichner T
    3. Schlattl A
    4. Stütz AM
    5. Benes V
    6. Korbel JO

    (2012) Delly: Structural variant discovery by integrated paired-end and split-read analysis

    Bioinformatics 28:i333–i339.

    https://doi.org/10.1093/bioinformatics/bts378

    • Google Scholar
83. 1. Rubin CJ
    2. Zody MC
    3. Eriksson J
    4. Meadows JR
    5. Sherwood E
    6. Webster MT
    7. Jiang L
    8. Ingman M
    9. Sharpe T
    10. Ka S
    11. Hallböök F
    12. Besnier F
    13. Carlborg O
    14. Bed'hom B
    15. Tixier-Boichard M
    16. Jensen P
    17. Siegel P
    18. Lindblad-Toh K
    19. Andersson L

    (2010) Whole-genome resequencing reveals loci under selection during chicken domestication

    Nature 464:587–591.

    https://doi.org/10.1038/nature08832

    • Google Scholar
84. 1. Ryman N
    2. Lagercrantz U
    3. Andersson L
    4. Chakraborty R
    5. Rosenberg R

    (1984) Lack of correspondence between genetic and morphologic variability patterns in Atlantic herring (Clupea harengus)

    Heredity 53:687–704.

    https://doi.org/10.1038/hdy.1984.127

    • Google Scholar
85. 1. Schennink A
    2. Trott JF
    3. Manjarin R
    4. Lemay DG
    5. Freking BA
    6. Hovey RC

    (2015) Comparative genomics reveals tissue-specific regulation of prolactin receptor gene expression

    Journal of Molecular Endocrinology 54:1–15.

    https://doi.org/10.1530/JME-14-0212

    • Google Scholar
86. 1. Sheehan S
    2. Harris K
    3. Song YS

    (2013) Estimating variable effective population sizes from multiple genomes: A sequentially Markov conditional sampling distribution approach

    Genetics 194:647–662.

    https://doi.org/10.1534/genetics.112.149096

    • Google Scholar
87. Software

    1. Smit A
    2. Hubley R

    (2010) Repeatmodeler open-1.0

    Repeatmodeler open-1.0.

    http://www.repeatmasker.org/RepeatModeler.html
88. Software

    1. Smit AFA
    2. Hubley R
    3. Green P

    (2015) Repeatmasker

    Repeatmasker.

    http://repeatmasker.org
89. 1. Sperber GO
    2. Airola T
    3. Jern P
    4. Blomberg J

    (2007) Automated recognition of retroviral sequences in genomic data--retrotector

    Nucleic Acids Research 35:4964–4976.

    https://doi.org/10.1093/nar/gkm515

    • Google Scholar
90. 1. Stanke M
    2. Diekhans M
    3. Baertsch R
    4. Haussler D

    (2008) Using native and syntenically mapped cDNA alignments to improve de novo gene finding

    Bioinformatics 24:637–644.

    https://doi.org/10.1093/bioinformatics/btn013

    • Google Scholar
91. 1. Star B
    2. Nederbragt AJ
    3. Jentoft S
    4. Grimholt U
    5. Malmstrøm M
    6. Gregers TF
    7. Rounge TB
    8. Paulsen J
    9. Solbakken MH
    10. Sharma A
    11. Wetten OF
    12. Lanzén A
    13. Winer R
    14. Knight J
    15. Vogel JH
    16. Aken B
    17. Andersen O
    18. Lagesen K
    19. Tooming-Klunderud A
    20. Edvardsen RB
    21. Tina KG
    22. Espelund M
    23. Nepal C
    24. Previti C
    25. Karlsen BO
    26. Moum T
    27. Skage M
    28. Berg PR
    29. Gjøen T
    30. Kuhl H
    31. Thorsen J
    32. Malde K
    33. Reinhardt R
    34. Du L
    35. Johansen SD
    36. Searle S
    37. Lien S
    38. Nilsen F
    39. Jonassen I
    40. Omholt SW
    41. Stenseth NC
    42. Jakobsen KS

    (2011) The genome sequence of Atlantic cod reveals a unique immune system

    Nature 477:207–210.

    https://doi.org/10.1038/nature10342

    • Google Scholar
92. Software

    1. Tate R
    2. Hall B
    3. Derego T
    4. Geib S

    (2014) Annie: The annotation information extractor

    Annie: The annotation information extractor.

    http://genomeannotation.github.io/Annie
93. 1. Trapnell C
    2. Williams BA
    3. Pertea G
    4. Mortazavi A
    5. Kwan G
    6. van Baren MJ
    7. Salzberg SL
    8. Wold BJ
    9. Pachter L

    (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation

    Nature Biotechnology 28:511–515.

    https://doi.org/10.1038/nbt.1621

    • Google Scholar
94. 1. Vinogradov AE

    (1998) Genome size and gc-percent in vertebrates as determined by flow cytometry: The triangular relationship

    Cytometry 31:100–109.

    https://doi.org/10.1002/(SICI)1097-0320(19980201)31:2<100::AID-CYTO5>3.0.CO;2-Q

    • Google Scholar
95. 1. Visser ME
    2. Gienapp P
    3. Husby A
    4. Morrisey M
    5. de la Hera I
    6. Pulido F
    7. Both C

    (2015) Effects of spring temperatures on the strength of selection on timing of reproduction in a long-distance migratory bird

    PLoS Biology 13:e1002120.

    https://doi.org/10.1371/journal.pbio.1002120

    • Google Scholar
96. 1. Voskoboynik A
    2. Neff NF
    3. Sahoo D
    4. Newman AM
    5. Pushkarev D
    6. Koh W
    7. Passarelli B
    8. Fan HC
    9. Mantalas GL
    10. Palmeri KJ
    11. Ishizuka KJ
    12. Gissi C
    13. Griggio F
    14. Ben-Shlomo R
    15. Corey DM
    16. Penland L
    17. White RA
    18. Weissman IL
    19. Quake SR

    (2013) The genome sequence of the colonial chordate, Botryllus schlosseri

    eLife 2:e00569.

    https://doi.org/10.7554/eLife.00569

    • Google Scholar
97. 1. Wang G
    2. Yang E
    3. Smith KJ
    4. Zeng Y
    5. Ji G
    6. Connon R
    7. Fangue NA
    8. Cai JJ

    (2014) Gene expression responses of threespine stickleback to salinity: Implications for salt-sensitive hypertension

    Frontiers in Genetics 5:312.

    https://doi.org/10.3389/fgene.2014.00312

    • Google Scholar
98. 1. Wood AR
    2. Esko T
    3. Yang J
    4. Vedantam S
    5. Pers TH
    6. Gustafsson S
    7. Chu AY
    8. Estrada K
    9. Luan J
    10. Kutalik Z
    11. Amin N
    12. Buchkovich ML
    13. Croteau-Chonka DC
    14. Day FR
    15. Duan Y
    16. Fall T
    17. Fehrmann R
    18. Ferreira T
    19. Jackson AU
    20. Karjalainen J
    21. Lo KS
    22. Locke AE
    23. Mägi R
    24. Mihailov E
    25. Porcu E
    26. Randall JC
    27. Scherag A
    28. Vinkhuyzen AA
    29. Westra HJ
    30. Winkler TW
    31. Workalemahu T
    32. Zhao JH
    33. Absher D
    34. Albrecht E
    35. Anderson D
    36. Baron J
    37. Beekman M
    38. Demirkan A
    39. Ehret GB
    40. Feenstra B
    41. Feitosa MF
    42. Fischer K
    43. Fraser RM
    44. Goel A
    45. Gong J
    46. Justice AE
    47. Kanoni S
    48. Kleber ME
    49. Kristiansson K
    50. Lim U
    51. Lotay V
    52. Lui JC
    53. Mangino M
    54. Mateo Leach I
    55. Medina-Gomez C
    56. Nalls MA
    57. Nyholt DR
    58. Palmer CD
    59. Pasko D
    60. Pechlivanis S
    61. Prokopenko I
    62. Ried JS
    63. Ripke S
    64. Shungin D
    65. Stancáková A
    66. Strawbridge RJ
    67. Sung YJ
    68. Tanaka T
    69. Teumer A
    70. Trompet S
    71. van der Laan SW
    72. van Setten J
    73. Van Vliet-Ostaptchouk JV
    74. Wang Z
    75. Yengo L
    76. Zhang W
    77. Afzal U
    78. Arnlöv J
    79. Arscott GM
    80. Bandinelli S
    81. Barrett A
    82. Bellis C
    83. Bennett AJ
    84. Berne C
    85. Blüher M
    86. Bolton JL
    87. Böttcher Y
    88. Boyd HA
    89. Bruinenberg M
    90. Buckley BM
    91. Buyske S
    92. Caspersen IH
    93. Chines PS
    94. Clarke R
    95. Claudi-Boehm S
    96. Cooper M
    97. Daw EW
    98. De Jong PA
    99. Deelen J
    100. Delgado G
    101. Denny JC
    102. Dhonukshe-Rutten R
    103. Dimitriou M
    104. Doney AS
    105. Dörr M
    106. Eklund N
    107. Eury E
    108. Folkersen L
    109. Garcia ME
    110. Geller F
    111. Giedraitis V
    112. Go AS
    113. Grallert H
    114. Grammer TB
    115. Gräßler J
    116. Grönberg H
    117. de Groot LC
    118. Groves CJ
    119. Haessler J
    120. Hall P
    121. Haller T
    122. Hallmans G
    123. Hannemann A
    124. Hartman CA
    125. Hassinen M
    126. Hayward C
    127. Heard-Costa NL
    128. Helmer Q
    129. Hemani G
    130. Henders AK
    131. Hillege HL
    132. Hlatky MA
    133. Hoffmann W
    134. Hoffmann P
    135. Holmen O
    136. Houwing-Duistermaat JJ
    137. Illig T
    138. Isaacs A
    139. James AL
    140. Jeff J
    141. Johansen B
    142. Johansson Å
    143. Jolley J
    144. Juliusdottir T
    145. Junttila J
    146. Kho AN
    147. Kinnunen L
    148. Klopp N
    149. Kocher T
    150. Kratzer W
    151. Lichtner P
    152. Lind L
    153. Lindström J
    154. Lobbens S
    155. Lorentzon M
    156. Lu Y
    157. Lyssenko V
    158. Magnusson PK
    159. Mahajan A
    160. Maillard M
    161. McArdle WL
    162. McKenzie CA
    163. McLachlan S
    164. McLaren PJ
    165. Menni C
    166. Merger S
    167. Milani L
    168. Moayyeri A
    169. Monda KL
    170. Morken MA
    171. Müller G
    172. Müller-Nurasyid M
    173. Musk AW
    174. Narisu N
    175. Nauck M
    176. Nolte IM
    177. Nöthen MM
    178. Oozageer L
    179. Pilz S
    180. Rayner NW
    181. Renstrom F
    182. Robertson NR
    183. Rose LM
    184. Roussel R
    185. Sanna S
    186. Scharnagl H
    187. Scholtens S
    188. Schumacher FR
    189. Schunkert H
    190. Scott RA
    191. Sehmi J
    192. Seufferlein T
    193. Shi J
    194. Silventoinen K
    195. Smit JH
    196. Smith AV
    197. Smolonska J
    198. Stanton AV
    199. Stirrups K
    200. Stott DJ
    201. Stringham HM
    202. Sundström J
    203. Swertz MA
    204. Syvänen AC
    205. Tayo BO
    206. Thorleifsson G
    207. Tyrer JP
    208. van Dijk S
    209. van Schoor NM
    210. van der Velde N
    211. van Heemst D
    212. van Oort FV
    213. Vermeulen SH
    214. Verweij N
    215. Vonk JM
    216. Waite LL
    217. Waldenberger M
    218. Wennauer R
    219. Wilkens LR
    220. Willenborg C
    221. Wilsgaard T
    222. Wojczynski MK
    223. Wong A
    224. Wright AF
    225. Zhang Q
    226. Arveiler D
    227. Bakker SJ
    228. Beilby J
    229. Bergman RN
    230. Bergmann S
    231. Biffar R
    232. Blangero J
    233. Boomsma DI
    234. Bornstein SR
    235. Bovet P
    236. Brambilla P
    237. Brown MJ
    238. Campbell H
    239. Caulfield MJ
    240. Chakravarti A
    241. Collins R
    242. Collins FS
    243. Crawford DC
    244. Cupples LA
    245. Danesh J
    246. de Faire U
    247. den Ruijter HM
    248. Erbel R
    249. Erdmann J
    250. Eriksson JG
    251. Farrall M
    252. Ferrannini E
    253. Ferrières J
    254. Ford I
    255. Forouhi NG
    256. Forrester T
    257. Gansevoort RT
    258. Gejman PV
    259. Gieger C
    260. Golay A
    261. Gottesman O
    262. Gudnason V
    263. Gyllensten U
    264. Haas DW
    265. Hall AS
    266. Harris TB
    267. Hattersley AT
    268. Heath AC
    269. Hengstenberg C
    270. Hicks AA
    271. Hindorff LA
    272. Hingorani AD
    273. Hofman A
    274. Hovingh GK
    275. Humphries SE
    276. Hunt SC
    277. Hypponen E
    278. Jacobs KB
    279. Jarvelin MR
    280. Jousilahti P
    281. Jula AM
    282. Kaprio J
    283. Kastelein JJ
    284. Kayser M
    285. Kee F
    286. Keinanen-Kiukaanniemi SM
    287. Kiemeney LA
    288. Kooner JS
    289. Kooperberg C
    290. Koskinen S
    291. Kovacs P
    292. Kraja AT
    293. Kumari M
    294. Kuusisto J
    295. Lakka TA
    296. Langenberg C
    297. Le Marchand L
    298. Lehtimäki T
    299. Lupoli S
    300. Madden PA
    301. Männistö S
    302. Manunta P
    303. Marette A
    304. Matise TC
    305. McKnight B
    306. Meitinger T
    307. Moll FL
    308. Montgomery GW
    309. Morris AD
    310. Morris AP
    311. Murray JC
    312. Nelis M
    313. Ohlsson C
    314. Oldehinkel AJ
    315. Ong KK
    316. Ouwehand WH
    317. Pasterkamp G
    318. Peters A
    319. Pramstaller PP
    320. Price JF
    321. Qi L
    322. Raitakari OT
    323. Rankinen T
    324. Rao DC
    325. Rice TK
    326. Ritchie M
    327. Rudan I
    328. Salomaa V
    329. Samani NJ
    330. Saramies J
    331. Sarzynski MA
    332. Schwarz PE
    333. Sebert S
    334. Sever P
    335. Shuldiner AR
    336. Sinisalo J
    337. Steinthorsdottir V
    338. Stolk RP
    339. Tardif JC
    340. Tönjes A
    341. Tremblay A
    342. Tremoli E
    343. Virtamo J
    344. Vohl MC
    345. Amouyel P
    346. Asselbergs FW
    347. Assimes TL
    348. Bochud M
    349. Boehm BO
    350. Boerwinkle E
    351. Bottinger EP
    352. Bouchard C
    353. Cauchi S
    354. Chambers JC
    355. Chanock SJ
    356. Cooper RS
    357. de Bakker PI
    358. Dedoussis G
    359. Ferrucci L
    360. Franks PW
    361. Froguel P
    362. Groop LC
    363. Haiman CA
    364. Hamsten A
    365. Hayes MG
    366. Hui J
    367. Hunter DJ
    368. Hveem K
    369. Jukema JW
    370. Kaplan RC
    371. Kivimaki M
    372. Kuh D
    373. Laakso M
    374. Liu Y
    375. Martin NG
    376. März W
    377. Melbye M
    378. Moebus S
    379. Munroe PB
    380. Njølstad I
    381. Oostra BA
    382. Palmer CN
    383. Pedersen NL
    384. Perola M
    385. Pérusse L
    386. Peters U
    387. Powell JE
    388. Power C
    389. Quertermous T
    390. Rauramaa R
    391. Reinmaa E
    392. Ridker PM
    393. Rivadeneira F
    394. Rotter JI
    395. Saaristo TE
    396. Saleheen D
    397. Schlessinger D
    398. Slagboom PE
    399. Snieder H
    400. Spector TD
    401. Strauch K
    402. Stumvoll M
    403. Tuomilehto J
    404. Uusitupa M
    405. van der Harst P
    406. Völzke H
    407. Walker M
    408. Wareham NJ
    409. Watkins H
    410. Wichmann HE
    411. Wilson JF
    412. Zanen P
    413. Deloukas P
    414. Heid IM
    415. Lindgren CM
    416. Mohlke KL
    417. Speliotes EK
    418. Thorsteinsdottir U
    419. Barroso I
    420. Fox CS
    421. North KE
    422. Strachan DP
    423. Beckmann JS
    424. Berndt SI
    425. Boehnke M
    426. Borecki IB
    427. McCarthy MI
    428. Metspalu A
    429. Stefansson K
    430. Uitterlinden AG
    431. van Duijn CM
    432. Franke L
    433. Willer CJ
    434. Price AL
    435. Lettre G
    436. Loos RJ
    437. Weedon MN
    438. Ingelsson E
    439. O'Connell JR
    440. Abecasis GR
    441. Chasman DI
    442. Goddard ME
    443. Visscher PM
    444. Hirschhorn JN
    445. Frayling TM
    446. Electronic Medical Records and Genomics (eMEMERGEGE) Consortium
    447. MIGen Consortium
    448. PAGEGE Consortium
    449. LifeLines Cohort Study

    (2014) Defining the role of common variation in the genomic and biological architecture of adult human height

    Nature Genetics 46:1173–1186.

    https://doi.org/10.1038/ng.3097

    • Google Scholar
99. 1. Worm B
    2. Barbier EB
    3. Beaumont N
    4. Duffy JE
    5. Folke C
    6. Halpern BS
    7. Jackson JB
    8. Lotze HK
    9. Micheli F
    10. Palumbi SR
    11. Sala E
    12. Selkoe KA
    13. Stachowicz JJ
    14. Watson R

    (2006) Impacts of biodiversity loss on ocean ecosystem services

    Science 314:787–790.

    https://doi.org/10.1126/science.1132294

    • Google Scholar

Author details

1. Alvaro Martinez Barrio

   1. Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
   2. Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden

   Contribution

   AMB, Contributed to the assembly and genome annotation, Performed bioinformatic analysis of population data, Wrote the paper

   Contributed equally with

   Sangeet Lamichhaney, Guangyi Fan and Nima Rafati

   Competing interests

   The authors declare that no competing interests exist.

   "This ORCID iD identifies the author of this article:" 0000-0001-5064-2093

2. Sangeet Lamichhaney

   Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

   Contribution

   SL, Performed bioinformatic analysis of population data, Wrote the paper

   Contributed equally with

   Alvaro Martinez Barrio, Guangyi Fan and Nima Rafati

   Competing interests

   The authors declare that no competing interests exist.

   "This ORCID iD identifies the author of this article:" 0000-0003-4826-0349

3. Guangyi Fan

   1. State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau, China
   2. BGI-Shenzhen, Shenzen, China

   Contribution

   GF, Contributed to the assembly and genome annotation

   Contributed equally with

   Alvaro Martinez Barrio, Sangeet Lamichhaney and Nima Rafati

   Competing interests

   The authors declare that no competing interests exist.

4. Nima Rafati

   Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

   Contribution

   NRa, Performed bioinformatic analysis of population data

   Contributed equally with

   Alvaro Martinez Barrio, Sangeet Lamichhaney and Guangyi Fan

   Competing interests

   The authors declare that no competing interests exist.

   "This ORCID iD identifies the author of this article:" 0000-0002-3687-9745

5. Mats Pettersson

   Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

   Contribution

   MP, Performed bioinformatic analysis of population data

   Competing interests

   The authors declare that no competing interests exist.

6. He Zhang

   1. BGI-Shenzhen, Shenzen, China
   2. College of Physics, Qingdao University, Qingdao, China

   Contribution

   HZ, Contributed to the assembly and genome annotation

   Competing interests

   The authors declare that no competing interests exist.

7. Jacques Dainat

   1. Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
   2. Bioinformatics Infrastructure for Life Sciences, Uppsala University, Uppsala, Sweden

   Contribution

   JD, Contributed to the assembly and genome annotation

   Competing interests

   The authors declare that no competing interests exist.

8. Diana Ekman

   Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden

   Contribution

   DE, Contributed to the assembly and genome annotation

   Competing interests

   The authors declare that no competing interests exist.

9. Marc Höppner

   1. Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
   2. Bioinformatics Infrastructure for Life Sciences, Uppsala University, Uppsala, Sweden

   Contribution

   MH, Contributed to the assembly and genome annotation

   Competing interests

   The authors declare that no competing interests exist.

10. Patric Jern

    Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

    Contribution

    PJ, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

11. Marcel Martin

    Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden

    Contribution

    MM, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

    "This ORCID iD identifies the author of this article:" 0000-0002-0680-200X

12. Björn Nystedt

    Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden

    Contribution

    BN, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

13. Xin Liu

    BGI-Shenzhen, Shenzen, China

    Contribution

    XLi, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

14. Wenbin Chen

    BGI-Shenzhen, Shenzen, China

    Contribution

    WC, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

15. Xinming Liang

    BGI-Shenzhen, Shenzen, China

    Contribution

    XLia, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

16. Chengcheng Shi

    BGI-Shenzhen, Shenzen, China

    Contribution

    CS, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

17. Yuanyuan Fu

    1. BGI-Shenzhen, Shenzen, China
    2. School of Biological Science and Medical Engineering, Southeast University, Nanjing, China

    Contribution

    YF, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

18. Kailong Ma

    BGI-Shenzhen, Shenzen, China

    Contribution

    KM, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

19. Xiao Zhan

    BGI-Shenzhen, Shenzen, China

    Contribution

    XZ, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

20. Chungang Feng

    Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

    Contribution

    CF, Performed experimental work

    Competing interests

    The authors declare that no competing interests exist.

21. Ulla Gustafson

    Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden

    Contribution

    UG, Performed experimental work

    Competing interests

    The authors declare that no competing interests exist.

22. Carl-Johan Rubin

    Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

    Contribution

    C-JR, Performed bioinformatic analysis of population data

    Competing interests

    The authors declare that no competing interests exist.

23. Markus Sällman Almén

    Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden

    Contribution

    MSA, Analysis and interpretation of data, Contributed unpublished essential data or reagents

    Competing interests

    The authors declare that no competing interests exist.

24. Martina Blass

    Department of Aquatic Resources, Institute of Coastal Research, Swedish University of Agricultural Sciences, Öregrund, Sweden

    Contribution

    MB, Provided population samples

    Competing interests

    The authors declare that no competing interests exist.

25. Michele Casini

    Department of Aquatic Resources, Institute of Marine Research, Swedish University of Agricultural Sciences, Lysekil, Sweden

    Contribution

    MC, Provided population samples

    Competing interests

    The authors declare that no competing interests exist.

26. Arild Folkvord

    1. Department of Biology, University of Bergen, Bergen, Norway
    2. Hjort Center of Marine Ecosystem Dynamics, Bergen, Norway
    3. Institute of Marine Research, Bergen, Norway

    Contribution

    AF, Provided population samples

    Competing interests

    The authors declare that no competing interests exist.

27. Linda Laikre

    Department of Zoology, Stockholm University, Stockholm, Sweden

    Contribution

    LL, Provided population samples

    Competing interests

    The authors declare that no competing interests exist.

28. Nils Ryman

    Department of Zoology, Stockholm University, Stockholm, Sweden

    Contribution

    NRy, Provided population samples

    Competing interests

    The authors declare that no competing interests exist.

29. Simon Ming-Yuen Lee

    State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau, China

    Contribution

    SM-YL, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

30. Xun Xu

    BGI-Shenzhen, Shenzen, China

    Contribution

    XX, Contributed to the assembly and genome annotation

    Competing interests

    The authors declare that no competing interests exist.

31. Leif Andersson

    1. Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
    2. Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
    3. Department of Veterinary Integrative Biosciences, Texas A&M University, Texas, United States

    Contribution

    LA, Conceived the study, Led bioinformatic analysis of population data, Wrote the paper, Acquisition of data, Analysis and interpretation of data

    For correspondence

    leif.andersson@imbim.uu.se

    Competing interests

    The authors declare that no competing interests exist.

    "This ORCID iD identifies the author of this article:" 0000-0002-4085-6968

• 7,429

  views

• 1,212

  downloads

• 138

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.