Rhodopsin Gene Polymorphism Associated with Divergent Light Environments in Atlantic Cod

Genetic assignment predicts depth of benthic settlement for 0-group Atlantic cod

Atlantic cod is a keystone species that remains among the most economically important demersal fish in the North Atlantic. Throughout its distribution range, Atlantic cod is composed of populations with varying environmental preferences and migratory propensities. This life-history variation is likely to have contributed to the niche width and large population sizes of Atlantic cod, and its relative resilience to environmental change and exploitation. The Icelandic cod stock is currently managed as a single unit, but early research indicates population variation by depth and temperature and distinct offshore and inshore spawning components. Pelagic 0-group juveniles from different spawning grounds coexist in nursery areas around Iceland, but their genetic composition or habitat partitioning had not been examined post benthic settlement. In the current study we examine the genetic composition of Atlantic cod juvenile aggregations at nearshore nursery grounds in NW-Iceland and report distinct segregation by the depth of offshore and inshore juvenile cod. The physiological mechanism of this segregation is not known, but the pattern demonstrates the need to consider population structure at nursery grounds in the application of marine spatial planning and other area-based conservation tools.

Hidden but revealed: After years of genetic studies behavioural monitoring combined with genomics uncover new insight into the population dynamics of Atlantic cod in Icelandic waters

Stock structure is of paramount importance for sustainable management of exploited resources. In that context, genetic markers have been used for more than two decades to resolve spatial structure of marine exploited resources and to fully fathom stock dynamics and interactions. While genetic markers such as allozymes and RFLP dominated the debate in the early era of genetics, technology advances have provided scientists with new tools every decade to better assess stock discrimination and interactions (i.e. gene flow). Here, we provide a review of genetic studies performed to understand stock structure of Atlantic cod in Icelandic waters, from the early allozyme approaches to the genomic work currently carried out. We further highlight the importance of the generation of a chromosome-anchored genome assembly together with whole-genome population data, which drastically changed our perception of the possible management units to consider. After nearly 60 years of genetic investigation of Atlantic cod structure in Icelandic waters, genetic (and later genomic) data combined with behavioural monitoring using Data Storage Tags shifted the attention from geographical population structures to behavioural ecotypes. This review also demonstrates the need for future research to further disentangle the impact of these ecotypes (and gene flow among them) on the population structure of Atlantic cod in Icelandic waters. It also highlights the importance of whole-genome data to unravel unexpected within-species diversity related to chromosomal inversions and associated supergenes, which are important to consider for future development of sustainable management programmes of the species within the North Atlantic.

Genomic differentiation in Pacific cod using Pool-Seq

Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (F ST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (d XY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (d XY) but roughly double the relative divergence (F ST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.

The pantophysin gene and its relationship with survival in early life stages of Atlantic cod

Genetic markers are widely used in fisheries management around the world. While the genetic structure and markers selected are usually based on samples from the wild, very few controlled experiments have been carried out to investigate possible differences in influence on traits between markers. Here we examine the bi-allelic gene pantophysin (Pan I), widely used in the management of Atlantic cod, in a series of in vitro crosses under a range of temperatures. It has been proposed that this gene, or another tightly linked gene, may be under strong divergent selection. Resolving this issue is essential in order to interpret results when using this gene marker for stock management. We found no evidence of departure from the expected 1 : 2 : 1 Mendelian ratio for any of the three genotypes during the egg stage, while both the 6 and 12°C temperature regimes in tank experiments favoured the survival of the Pan IAA genotype. No difference in genotype survival was, however, found in a more natural mesocosm environment. Collectively, these results suggest that for the early life stages of Atlantic cod, and under the current experimental conditions, there is no strong consistent influence of Pan I genotype on survival. The results also emphasize the importance of varied experimental studies to verify the importance of environmental factors influencing genotype selection.

Identifying loci under positive selection in complex population histories

Detailed modeling of a species' history is of prime importance for understanding how natural selection operates over time. Most methods designed to detect positive selection along sequenced genomes, however, use simplified representations of past histories as null models of genetic drift. Here, we present the first method that can detect signatures of strong local adaptation across the genome using arbitrarily complex admixture graphs, which are typically used to describe the history of past divergence and admixture events among any number of populations. The method-called graph-aware retrieval of selective sweeps (GRoSS)-has good power to detect loci in the genome with strong evidence for past selective sweeps and can also identify which branch of the graph was most affected by the sweep. As evidence of its utility, we apply the method to bovine, codfish, and human population genomic data containing panels of multiple populations related in complex ways. We find new candidate genes for important adaptive functions, including immunity and metabolism in understudied human populations, as well as muscle mass, milk production, and tameness in specific bovine breeds. We are also able to pinpoint the emergence of large regions of differentiation owing to inversions in the history of Atlantic codfish.

Environmental, population and life-stage plasticity in the visual system of Atlantic cod

The visual system is for many fishes essential in guiding behaviors, such as foraging, predator avoidance and mate choice. The marine environment is characterized by large spatio-temporal fluctuations in light intensity and spectral composition. However, visual capabilities are restricted by both space limitations set by eye size and by the genomic content of light-absorbing opsin genes. The rich array of visual opsins in teleosts may be used differentially to tune vision towards specific needs during ontogeny and to changing light. Yet, to what extent visual plasticity is a pre-programmed developmental event, or is triggered by photic environment, is unclear. Our previous studies on Atlantic cod revealed an evolutionary genomic loss of UV-sensitive sws1 and red-sensitive lws opsin families, while blue-sensitive sws2 and green-sensitive rh2 opsins had duplicated. The current study has taken an opsin expression approach to characterize visual plasticity in cod towards different spectral light during the larval stage, to maturation and extreme seasonal changes in the Barents Sea. Our data suggest that opsin plasticity in cod larvae is controlled by developmental programme rather than immediate light environment. The lack of expressional changes during maturation suggests a less important role for visual modulation related to mate choice. Although no seasonal effects on visual opsins were detected in migratory Northeast Arctic cod, the expressed opsin subset differed from the more stationary Norwegian coastal cod described in previous studies. Interestingly, these data provide the first indications of a population difference in actively used visual opsins associated with cod ecotypes.

Trans-oceanic genomic divergence of Atlantic cod ecotypes is associated with large inversions

Chromosomal rearrangements such as inversions can play a crucial role in maintaining polymorphism underlying complex traits and contribute to the process of speciation. In Atlantic cod (Gadus morhua), inversions of several megabases have been identified that dominate genomic differentiation between migratory and nonmigratory ecotypes in the Northeast Atlantic. Here, we show that the same genomic regions display elevated divergence and contribute to ecotype divergence in the Northwest Atlantic as well. The occurrence of these inversions on both sides of the Atlantic Ocean reveals a common evolutionary origin, predating the >100 000-year-old trans-Atlantic separation of Atlantic cod. The long-term persistence of these inversions indicates that they are maintained by selection, possibly facilitated by coevolution of genes underlying complex traits. Our data suggest that migratory behaviour is derived from more stationary, ancestral ecotypes. Overall, we identify several large genomic regions—each containing hundreds of genes—likely involved in the maintenance of genomic divergence in Atlantic cod on both sides of the Atlantic Ocean.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

Three chromosomal rearrangements promote genomic divergence between migratory and stationary ecotypes of Atlantic cod

Identification of genome-wide patterns of divergence provides insight on how genomes are influenced by selection and can reveal the potential for local adaptation in spatially structured populations. In Atlantic cod – historically a major marine resource – Northeast-Arctic- and Norwegian coastal cod are recognized by fundamental differences in migratory and non-migratory behavior, respectively. However, the genomic architecture underlying such behavioral ecotypes is unclear. Here, we have analyzed more than 8.000 polymorphic SNPs distributed throughout all 23 linkage groups and show that loci putatively under selection are localized within three distinct genomic regions, each of several megabases long, covering approximately 4% of the Atlantic cod genome. These regions likely represent genomic inversions. The frequency of these distinct regions differ markedly between the ecotypes, spawning in the vicinity of each other, which contrasts with the low level of divergence in the rest of the genome. The observed patterns strongly suggest that these chromosomal rearrangements are instrumental in local adaptation and separation of Atlantic cod populations, leaving footprints of large genomic regions under selection. Our findings demonstrate the power of using genomic information in further understanding the population dynamics and defining management units in one of the world’s most economically important marine resources.