A baseline for the genetic stock identification of Atlantic herring, Clupea harengus, in ICES Divisions 6.a, 7.b-c.

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy.

Increasing temperature and prey availability affect the growth and swimming kinematics of Atlantic herring (Clupea harengus) larvae.

Climate change will increase the magnitude and duration of warming events and the variability in the phenology and abundance of available prey to the early life stages of fish. These factors influence physiological, behavioral and ecological processes, impacting growth, development and survival. Using a fully factorial design with two prey-availability treatments (1200 prey items L-1 (high prey abundance) or 40 prey items L-1 (low prey abundance)) under three temperature regimes (8, 10 and 12°C), the swimming kinematics of 6-week old spring-spawning Atlantic herring larvae were examined using silhouette video photography. Higher temperatures combined with food limitation significantly decreased the growth and swimming kinematics of larval herring, with the most negative effect observed in larvae reared at 12°C and exposed to low food abundances. Specifically, larvae displayed reduced locomotory behaviors and reduced vertical movements. By contrast, larvae reared at high prey abundance and at 12°C displayed more active swimming and exploratory behavior, as evidenced by an increase in both locomotory behavior and vertical and horizontal turn angles, suggesting increased motivation to search for food. This research highlights the importance of determining to what degree fish larvae are sensitive to changes in temperature and how these changes might be further influenced by food availability.

Tracking oocyte development and the timing of skipped spawning for north-east Arctic haddock (Melanogrammus aeglefinus).

The present study tracked oocyte development over 9 months and noted incidences of 'skipping', i.e., adults terminating their upcoming reproductive cycle, in field-caught north-east Arctic (NEA) haddock (Melanogrammus aeglefinus), currently the largest stock of this species. Applications of advanced image and histological techniques revealed the presence of cortical alveoli oocytes (CAO), which prevailed as the most advanced oocyte phase for 4-5 months. This new finding of an extended and early appearance of CAOs in this gadoid was supported by that vitellogenesis first started to appear 3 months later. The subsequent oocyte growth trajectories indicated that larger individuals [total length (TL) = 70 cm] typically spawn in the order of 3 weeks earlier than the smaller ones (TL = 40 cm). The spawning season appeared stretched over about 3 months. The majority of skipping females arrested oocyte growth at the CAO phase followed by atretic reabsorption. Compared to those individuals maturing for the spawning season, 'skippers' generally exhibited lower body condition, characterized also by relatively lower liver sizes at the time of the main spawning season. This study demonstrated well-developed skipping dynamics, but also that the CAO period, i.e., when skipping takes place, may be exceedingly long in this commercially valuable gadoid and that its reproductive cycle in many ways deviates from that of the data-rich, sympatric NEA cod (Gadus morhua).

Quantitative molecular detection of larval Atlantic herring (Clupea harengus) in stomach contents of Atlantic mackerel (Scomber scombrus) marks regions of predation pressure.

Mortality rates in the early life-history stages of fishes are generally high yet identifying the causes remain unclear. Faltering recruitment rates of Atlantic herring (Clupea harengus) in the Norwegian Sea indicate a need to identify which mortality factors influence larval herring survival. Previous research suggests that increased predation pressure by Atlantic mackerel (Scomber scombrus) may contribute to the disconnect between spawning stock biomass and recruitment. To quantify the contribution of predation pressure by Atlantic mackerel to herring larval mortality, two research cruises were conducted within a probable "hot spot" (67-72° N) for intensified mackerel predation based on particle drift simulations. Mackerel stomach contents were analysed for herring larvae content using droplet digital polymerase chain reaction (ddPCR) with a quantitative molecular detection assay specific for herring. The ddPCR results demonstrate clear predation by mackerel on herring larvae and also suggest that the alternative use of visual examination may give misleading results. Our results show that mackerel should be considered a potentially important predator on herring larvae. The quantitative molecular assay presented here shows great promise as an efficient and specific tool to correctly identify and quantify predation pressure on early life-history stages of fishes.

Development stage distribution as a proxy for feeding success and growth for first feeding Norwegian spring spawning herring larvae.

The estimation of growth rates in young herring larvae (Clupea harengus) in the field can be difficult because the primary increments in the otoliths may not be discernible or formed at a daily level. Likewise, the estimation of mortality rates of fish larvae in the field is very difficult to achieve, especially in a rigorous quantitative manner. In this study, the authors suggest the use of a stage-based proxy of feeding success, growth and potential survival or mortality risk of field-caught larvae. The stage-based proxy is derived based on observations from previous laboratory studies where larvae successfully completing start-feeding on external food sources will advance through the early development stages, whereas those that do not (unsuccessful larvae) remain and accumulate in the development stage preceding first feeding. The relative occurrence of larvae in the early development stages is therefore expected to reflect feeding conditions of the larvae, with higher ratios of unsuccessful larvae indicative of poor feeding success and higher mortality risk. Using field data on Norwegian spring spawning herring, the authors document that the relative occurrence of larvae in the late non-feeding stage is significantly higher at lower average zooplankton concentrations, in line with the predictions of the authors that this novel approach of using a stage-based proxy could be a useful indication of feeding success, growth and mortality in the field. Further, there was a significant interaction effect with ambient temperature, with the ratio being higher at low zooplankton concentrations at higher temperatures. This study also suggests that these findings are not population specific as the same accumulation of non-feeding larvae in the late non-feeding stage was observed in laboratory-reared larvae of both autumn and spring spawning herring populations.

Ecological adaptation in Atlantic herring is associated with large shifts in allele frequencies at hundreds of loci.

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection.

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.

Reconstruction of the birth of a male sex chromosome present in Atlantic herring.

The mechanisms underlying sex determination are astonishingly plastic. Particularly the triggers for the molecular machinery, which recalls either the male or female developmental program, are highly variable and have evolved independently and repeatedly. Fish show a huge variety of sex determination systems, including both genetic and environmental triggers. The advent of sex chromosomes is assumed to stabilize genetic sex determination. However, because sex chromosomes are notoriously cluttered with repetitive DNA and pseudogenes, the study of their evolution is hampered. Here we reconstruct the birth of a Y chromosome present in the Atlantic herring. The region is tiny (230 kb) and contains only three intact genes. The candidate male-determining gene BMPR1BBY encodes a truncated form of a BMP1B receptor, which originated by gene duplication and translocation and underwent rapid protein evolution. BMPR1BBY phosphorylates SMADs in the absence of ligand and thus has the potential to induce testis formation. The Y region also contains two genes encoding subunits of the sperm-specific Ca2+ channel CatSper required for male fertility. The herring Y chromosome conforms with a characteristic feature of many sex chromosomes, namely, suppressed recombination between a sex-determining factor and genes that are beneficial for the given sex. However, the herring Y differs from other sex chromosomes in that suppression of recombination is restricted to an ∼500-kb region harboring the male-specific and sex-associated regions. As a consequence, any degeneration on the herring Y chromosome is restricted to those genes located in the small region affected by suppressed recombination.

First estimates of metabolic rate in Atlantic bluefin tuna larvae.

Atlantic bluefin tuna is an iconic scombrid species with a high commercial and ecological value. Despite their importance, many physiological aspects, especially during the larval stages, are still unknown. Metabolic rates are one of the understudied aspects in scombrid larvae, likely due to challenges associated to larval handling before and during respirometry trials. Gaining reliable estimates of metabolic rates is essential to understand how larvae balance their high growth needs and activity and other physiological functions, which can be very useful for fisheries ecology and aquaculture. This is the first study to (a) estimate the relationship between routine metabolic rate (RMR) and larval dry weight (DW) (mass scaling exponent) at a constant temperature of 26°C, (b) measure the RMR under light and darkness and (c) test whether the interindividual differences in the RMR are related to larval nutritional status (RNA/DNA and DNA/DW). The RMR scaled nearly isometrically with body size (b = 0.99, 0.60-31.56 mg DW) in contrast to the allometric relationship observed in most fish larvae (average b = 0.87). The results show no significant differences in larval RMR under light and darkness, suggesting similar larval activity levels in both conditions. The size explained most of the variability in RMR (97%), and nutritional condition was unrelated to the interindividual differences in routine metabolism. This is the first study to report the metabolic rates of Atlantic bluefin tuna larvae and discuss the challenges of performing bioenergetic studies with early life stages of scombrids.

"A cleaner break": Genetic divergence between geographic groups and sympatric phenotypes revealed in ballan wrasse (Labrus bergylta).

Capture and long-distance translocation of cleaner fish to control lice infestations on marine salmonid farms has the potential to influence wild populations via overexploitation in source regions, and introgression in recipient regions. Knowledge of population genetic structure is therefore required. We studied the genetic structure of ballan wrasse, a phenotypically diverse and extensively used cleaner fish, from 18 locations in Norway and Sweden, and from Galicia, Spain, using 82 SNP markers. We detected two very distinct genetic groups in Scandinavia, northwestern and southeastern. These groups were split by a stretch of sandy beaches in southwest Norway, representing a habitat discontinuity for this rocky shore associated benthic egg-laying species. Wrasse from Galicia were highly differentiated from all Scandinavian locations, but more similar to northwestern than southeastern locations. Distinct genetic differences were observed between sympatric spotty and plain phenotypes in Galicia, but not in Scandinavia. The mechanisms underlying the geographic patterns between phenotypes are discussed, but not identified. We conclude that extensive aquaculture-mediated translocation of ballan wrasse from Sweden and southern Norway to western and middle Norway has the potential to mix genetically distinct populations. These results question the sustainability of the current cleaner fish practice.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection.

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.

Recurrent convergent evolution at amino acid residue 261 in fish rhodopsin.

The evolutionary process that occurs when a species colonizes a new environment provides an opportunity to explore the mechanisms underlying genetic adaptation, which is essential knowledge for understanding evolution and the maintenance of biodiversity. Atlantic herring has an estimated total breeding stock of about 1 trillion (1012) and has colonized the brackish Baltic Sea within the last 10,000 y. Minute genetic differentiation between Atlantic and Baltic herring populations at selectively neutral loci combined with this rapid adaptation to a new environment facilitated the identification of hundreds of loci underlying ecological adaptation. A major question in the field of evolutionary biology is to what extent such an adaptive process involves selection of novel mutations with large effects or genetic changes at many loci, each with a small effect on phenotype (i.e., selection on standing genetic variation). Here we show that a missense mutation in rhodopsin (Phe261Tyr) is an adaptation to the red-shifted Baltic Sea light environment. The transition from phenylalanine to tyrosine differs only by the presence of a hydroxyl moiety in the latter, but this results in an up to 10-nm red-shifted light absorbance of the receptor. Remarkably, an examination of the rhodopsin sequences from 2,056 species of fish revealed that the same missense mutation has occurred independently and been selected for during at least 20 transitions between light environments across all fish. Our results provide a spectacular example of convergent evolution and how a single amino acid change can have a major effect on ecological adaptation.

Genetic factors have a major effect on growth, number of vertebrae and otolith shape in Atlantic herring (Clupea harengus).

Atlantic herring, Clupea harengus, have complex population structures. Mixing of populations is known, but the extent of connectivity is still unclear. Phenotypic plasticity results in divergent phenotypes in response to environmental factors. A marked salinity gradient occurs from Atlantic Ocean (salinity 35) into the Baltic Sea (salinity range 2-12). Herring from both habitats display phenotypic and genetic variability. To explore how genetic factors and salinity influence phenotypic traits like growth, number of vertebrae and otolith shape an experimental population consisting of Atlantic purebreds and Atlantic/Baltic F1 hybrids were incubated and co-reared at two different salinities, 16 and 35, for three years. The F1-generation was repeatedly sampled to evaluate temporal variation. A von Bertalanffy growth model indicated that reared Atlantic purebreds had a higher maximum length (26.2 cm) than Atlantic/Baltic hybrids (24.8 cm) at salinity 35, but not at salinity 16 (25.0 and 24.8 cm, respectively). In contrast, Atlantic/Baltic hybrids achieved larger size-at-age than the wild caught Baltic parental group. Mean vertebral counts and otolith aspect ratios were higher for reared Atlantic purebreds than Atlantic/Baltic hybrids, consistent with the differences between parental groups. There were no significant differences in vertebral counts and otolith aspect ratios between herring with the same genotype but raised in different salinities. A Canonical Analysis of Principal Coordinates was applied to analyze the variation in wavelet coefficients that described otolith shape. The first discriminating axis identified the differences between Atlantic purebreds and Atlantic/Baltic hybrids, while the second axis represented salinity differences. Assigning otoliths based on genetic groups (Atlantic purebreds vs. Atlantic/Baltic hybrids) yielded higher classification success (~90%) than based on salinities (16 vs. 35; ~60%). Our results demonstrate that otolith shape and vertebral counts have a significant genetic component and are therefore useful for studies on population dynamics and connectivity.

Moderate nucleotide diversity in the Atlantic herring is associated with a low mutation rate.

The Atlantic herring is one of the most abundant vertebrates on earth but its nucleotide diversity is moderate (π = 0.3%), only three-fold higher than in human. Here, we present a pedigree-based estimation of the mutation rate in this species. Based on whole-genome sequencing of four parents and 12 offspring, the estimated mutation rate is 2.0 × 10-9 per base per generation. We observed a high degree of parental mosaicism indicating that a large fraction of these de novo mutations occurred during early germ cell development. The estimated mutation rate - the lowest among vertebrates analyzed to date - partially explains the discrepancy between the rather low nucleotide diversity in herring and its huge census population size. But a species like the herring will never reach its expected nucleotide diversity because of fluctuations in population size over the millions of years it takes to build up high nucleotide diversity.

The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing.

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.

Organ damage in Atlantic herring larvae as a result of ocean acidification.

The dissolution of anthropogenically emitted excess carbon dioxide lowers the pH of the world's ocean water. The larvae of mass spawning marine fishes may be particularly vulnerable to such ocean acidification (OA), yet the generality of earlier results is unclear. Here we show the detrimental effects of OA on the development of a commercially important fish species, the Atlantic herring (Clupea harengus). Larvae were reared at three levels of CO2: today (0.0385 kPa), end of next century (0.183 kPa), and a coastal upwelling scenario (0.426 kPa), under near-natural conditions in large outdoor tanks. Exposure to elevated CO2 levels resulted in stunted growth and development, decreased condition, and severe tissue damage in many organs, with the degree of damage increasing with CO2 concentration. This complements earlier studies of OA on Atlantic cod larvae that revealed similar organ damage but at increased growth rates and no effect on condition.

Experimental parameterisation of principal physics in buoyancy variations of marine teleost eggs.

It is generally accepted that the high buoyancy of pelagic marine eggs is due to substantial influx of water across the cell membrane just before ovulation. Here we further develop the theoretical basis by applying laboratory observations of the various components of the fertilized egg in first-principle equations for egg specific gravity (ρ(egg)) followed by statistical validation. We selected Atlantic cod as a model animal due to the affluent amount of literature on this species, but also undertook additional dedicated experimental works. We found that specific gravity of yolk plus embryo is central in influencing ρ(egg) and thereby the buoyancy. However, our established framework documents the effect on ρ(egg) of the initial deposition of the heavy chorion material in the gonad prior to spawning. Thereafter, we describe the temporal changes in ρ(egg) during incubation: Generally, the eggs showed a slight rise in ρ(egg) from fertilization to mid-gastrulation followed by a gradual decrease until full development of main embryonic organs just before hatching. Ontogenetic changes in ρ(egg) were significantly associated with volume and mass changes of yolk plus embryo. The initial ρ(egg) at fertilization appeared significantly influenced by the chorion volume fraction which is determined by the combination of the final chorion volume of the oocyte and of the degree of swelling (hydrolyzation) prior to spawning. The outlined principles and algorithms are universal in nature and should therefore be applicable to fish eggs in general.

Genetic ontogeny of pancreatic enzymes in Labrus bergylta larvae and the effect of feed type on enzyme activities and gene expression.

A newly cultivated wrasse species, Labrus bergylta, have shown great potential for use in Atlantic salmon (Salmo salar) farms in the battle against sea lice (Lepeoptheirus salmonis) infections. Hatchery reared L. bergylta were studied from 2 to 55 DPH to examine the molecular basis of digestive ontogeny related to the pancreas. An isolated feeding trial was performed on 27-34 DPH larvae to compare the effect of diet on enzyme activity and the possible exogenous contribution by live feed. The following genes coding for key pancreatic enzymes were analyzed by qPCR: trypsin, Cyp7 A1, BAL, sPLA(2) 1B, amylase and pancreatic chitinase. Enzyme activity was measured on trypsin, neutral lipase, sPLA(2), amylase and chitinase in fed and unfed larvae. We did not observe any effects of the formulated diet v.s. rotifers on enzyme activities of neutral lipase, chitinase and sPLA(2). However, a probable feed-dependency was observed at a transcriptional level, where rotifers seem to stimulate upregulation. The regulation of BAL was the only exception, where an upregulation was observed after weaning both in the ontogeny series and the experimental part. Our data on pancreatic chitinase and amylase mRNA levels suggest the importance of carbohydrates in the diet of early larval and juvenile L. bergylta.

Egg buoyancy variability in local populations of Atlantic cod (Gadus morhua).

Previous studies have found strong evidences for Atlantic cod (Gadus morhua) egg retention in fjords, which are caused by the combination of vertical salinity structure, estuarine circulation, and egg specific gravity, supporting small-scaled geographical differentiations of local populations. Here, we assess the variability in egg specific gravity for selected local populations of this species, that is, two fjord-spawning populations and one coastal-spawning population from Northern Norway (66-71°N/10-25°E). Eggs were naturally spawned by raised broodstocks (March to April 2009), and egg specific gravity was measured by a density-gradient column. The phenotype of egg specific gravity was similar among the three local populations. However, the associated variability was greater at the individual level than at the population level. The noted gradual decrease in specific gravity from gastrulation to hatching with an increase just before hatching could be a generic pattern in pelagic marine fish eggs. This study provides needed input to adequately understand and model fish egg dispersal.

Development of Atlantic cod (Gadus morhua) exposed to produced water during early life stages: Effects on embryos, larvae, and juvenile fish.

Produced water (PW) contains numerous toxic compounds of natural origin, such as dispersed oil, metals, alkylphenols (APs), and polycyclic aromatic hydrocarbons (PAHs). In addition, PW also contains many different chemicals which have been added during the oil production process. In the study described here, cod were exposed to real PW collected from an oil production platform in the North Sea. This was done in order to best recreate the most realistic field-exposure regime in which fish will be affected by a wide range of chemicals. The biological effects found in this study therefore cannot be assigned to one group of chemicals alone, but are the result of exposure to the complex chemical mixture found in real PW. Since APs are well known to cause endocrine disruption in marine organisms, we focused our chemical analysis on APs in an attempt to better understand the long-term effects of APs from PW on the biology of fish. In this study, cod were exposed to several concentrations of real PW and 17ß-oestradiol (E(2)), a natural oestrogen, at different developmental stages. Cod were exposed to PW either during the embryo and early larvae stage (up to 3 months of age) or during the early juvenile stage (from 3 to 6 months of age). Results showed that, in general, APs bioconcentrate in fish tissue in a dose and developmental stage dependent manner during PW exposure. However, juveniles appeared able to effectively metabolise the short chain APs. Importantly, PW exposure had no effect on embryo survival or hatching success. However, 1% PW clearly interfered with the development of normal larval pigmentation. After hatching most of the larvae exposed to 1% PW failed to begin feeding and died of starvation. This inability to feed may be linked to the increased incidence of jaw deformities seen in these larvae. In addition, cod exposed to 1% PW, had significantly higher levels of the biomarkers vitellogenin and CYP1A in plasma and liver, respectively. No similar effects were seen in cod exposed to either 0.1% or 0.01% PW.