Adult body growth and reproductive investment vary markedly within and across Atlantic and Pacific herring: a meta-analysis and review of 26 stocks

Life-history traits of Pacific (Clupea pallasii) and Atlantic (Clupea harengus) herring, comprising both local and oceanic stocks subdivided into summer-autumn and spring spawners, were extensively reviewed. The main parameters investigated were body growth, condition, and reproductive investment. Body size of Pacific herring increased with increasing latitude. This pattern was inconsistent for Atlantic herring. Pacific and local Norwegian herring showed comparable body conditions, whereas oceanic Atlantic herring generally appeared stouter. Among Atlantic herring, summer and autumn spawners produced many small eggs compared to spring spawners, which had fewer but larger eggs—findings agreeing with statements given several decades ago. The 26 herring stocks we analysed, when combined across distant waters, showed clear evidence of a trade-off between fecundity and egg size. The size-specific individual variation, often ignored, was substantial. Additional information on biometrics clarified that oceanic stocks were generally larger and had longer life spans than local herring stocks, probably related to their longer feeding migrations. Body condition was only weakly, positively related to assumingly in situ annual temperatures (0–30 m depth). Contrarily, body growth (cm × y−1), taken as an integrator of ambient environmental conditions, closely reflected the extent of investment in reproduction. Overall, Pacific and local Norwegian herring tended to cluster based on morphometric and reproductive features, whereas oceanic Atlantic herring clustered separately. Our work underlines that herring stocks are uniquely adapted to their habitats in terms of trade-offs between fecundity and egg size whereas reproductive investment mimics the productivity of the water in question.


Introduction
''Herrings'' consist of a complex group of phylogenetically primitive (in the sense of being ancient) fishes (Near et al. 2012) that have been further categorized at different hierarchical levels (Cushing 1967;Iles and Sinclair 1982;Rass and Wheeler 1991). At the highest (species) level, the Atlantic (Clupea harengus) and Pacific (C. pallasii) herring (Cushing 1967) are distinguished both genetically and morphologically (Libungan et al. 2016;Lamichhaney et al. 2017). In North European waters, however, a few herring populations are genetically related to Pacific herring (Jørstad et al. 1994). This distribution is attributed to inter-oceanic migration of Pacific herring to Atlantic waters during the Pleistocene (Laakkonen et al. 2013), which led to some degree of hybridization with local Atlantic herring (Laakkonen et al. 2015).
Plasticity in life-history parameters is common within all herring stocks (Iles and Sinclair 1982;Geffen 2009). Herring show group-synchronous oocyte development (Murua et al. 2003) resulting in both species in the production of demersal/benthic eggs-either shed at once (total spawner) or in smaller portions over hours or even days (Hay 1986;Kurita et al. 2003). This very short spawning period (at the individual level) speaks for that herring can be generally considered as true capital breeder (Kurita et al. 2003). Therefore, the number of oocytes estimated just prior to spawning reflects the total realized production, as further atresia (oocyte resorption) is unlikely (Ó skarsson et al. 2002). As a species, Atlantic herring spawn over the entire year, with each stock having a specific spawning season, and, as such, are frequently identified according to their specific spawning time (e.g., spring spawners and autumn spawners) and/or spawning grounds (Parrish and Saville 1965;Iles and Sinclair 1982;Sinclair and Tremblay 1984;Blaxter 1985). In contrast, Pacific herring are exclusively spring spawners, though the precise spawning time varies across their range, from mid-winter at the southern end to early summer at the extreme northern end (Blaxter 1985;Hay 1985). Egg size and fecundity differ between spawning seasons (Parrish and Saville 1965;Hempel and Blaxter 1967;van Damme et al. 2009), but egg production may also differ within the same spawning season across stocks and geographic regions (Silva et al. 2013;dos Santos Schmidt et al. 2017).
Herring stocks are also categorized by their spatial distribution and associated migratory behaviour (i.e., oceanic vs. local herring). The oceanic category contains large stock sizes that undergo long migrations, whereas the local category consists of relatively small stock sizes, which tend to be restricted to coastlines, fjords, or sea lochs (Parrish and Saville 1965;Holst et al. 2004). As the schools migrate between feeding, overwintering, and spawning areas on different spatial and seasonal scales (Lie et al. 1978;Varpe et al. 2005;Pampoulie et al. 2015), various stocks may overlap spatially Pampoulie et al. 2015;Berg et al. 2017). Furthermore, sympatric stocks that share common feeding and overwintering areas can have different spawning seasons (Jakobsson 1980;Husebø et al. 2005;van Damme et al. 2009;dos Santos Schmidt et al. 2017). So, as suggested by McQuinn (1997) and Ware and Schweigert (2001), herring fit well within the metapopulation concept due to the high complexity of stock dynamics and structure.
Based on the above outlines, the main objective of this article is to compare variation in observed adult life-history traits of different herring stocks distributed both in the North Atlantic and Northeast Pacific Oceans. The special focus was on reproduction and associated biometrics, key attributes for evaluating stock productivity and adaptation to environmental conditions, currently represented by downloaded annual temperature data sets. Such analyses are particularly relevant today because high-latitude waters have significantly been warmed (IPCC 2018) since the earliest herring investigations of this kind (see below). Further these analyses underscore the need for up to date and quality assured reference points considering future climate scenarios. In this metaanalysis we collated a modern database on 26 stocks, facilitated by international cooperation across marine laboratories. The data provided a basis for contrasting quantitative information of the various herring stocks studied, according to their reproductive biology including relationships between life-history, biometrics, reproductive trade-offs and environmental conditions. Because in situ external drivers (or possible stressors under climate change), such as relevant plankton abundance or ambient temperature could not be consistently reported across stocks (following a series of pilot analyses), we used body growth as an expression of experienced living conditions (Silva et al. 2013). We aimed to consider both fecundity and egg size jointly. If such pairs did not exist, we opted for ovary size as an indicator of reproductive investment. Further, we evaluated whether historic reproductive trade-offs in Atlantic herring persist today, including corresponding data on local Norwegian and Pacific herring. Any actual changes in egg size and fecundity over time (dos Santos Schmidt et al. 2017) were not pursued further due to examples of incomplete or lacking historical time series.

Material and methods
Herring stocks and the overall approach This work examined a series of life-history traits of adult herring females (Table 1) from 26 stocks inhabiting different regions in the Pacific and Atlantic Oceans (Fig. 1). In every case, the data were from samples taken prior to spawning but referred generally to a single season due to restricted access to reproductive time series. Examples of missing parameters are identified in Table 1. The material was grouped by geography: (1) Northeast Pacific, seven stocks off California to Alaska; (2) Northwest Atlantic, four stocks on the east coast of Canada and the US; (3) Northeast Atlantic, three stocks in central to eastern parts of the Nordic Seas, and two of the North Sea; (4) local Norwegian, eight stocks along the Norwegian coast; (5) Baltic Sea, one stock from the Central Baltic Sea; and (6) White Sea, the White Sea stock (WSH) ( Table 1). The data-rich Norwegian spring-spawning herring (NSSH) was used to exemplify the extent of annual variation in trait expression, highlighting seasons with known higher, typical, and lower fecundity levels and associated biometrics in response to varying abiotic and biotic conditions: 2007(N = 117), 1999(N = 38), and 2014dos Santos Schmidt et al. 2017) (Table 1). We also searched other databases for additional information on yearly variation in body condition because this metric impacts reproductive performance (Stearns 1992). As a result, data sets for 9 of the 26 evaluated stocks were extended by 1-2 additional years to give an impression of body condition variability. Length-at-age is considered a more resilient factor than body condition (Sande et al. 2019), especially in slow-growing fishes, such as herring (dos Santos Schmidt et al. 2017). Consequently, any interannual dynamics in body growth were only touched upon for NSSH. This decision was supported by a wealth of existing published information on growth dynamics for both Atlantic (Dragesund et al. 1980;Nash et al. 2005;Ó skarsson 2008) and Pacific herring (Tanasichuk et al. 1993;Hay et al. 2008;Ito et al. 2015), though less so for local Norwegian herring stocks (Hognestad 1994;da Silva et al. 2013;Johannessen et al. 2014) and across regions (Berg et al. 2017).

Individual metrics
Each herring was measured for length, either total length (TL, in cm) or, in the case of Pacific herring, standard length (SL, in mm). Gulf of Maine autumnspawning herring (GMH) measurements were based on fork length (FL, in mm). SL and FL in centimetres were converted to TL according to Karpov and Kwiecien (1988): TL = 0.766 ? 1.223 9 SL (r 2 = 0.994), and TL = 1.702 ? 1.109 9 FL (r 2 = 0.996). Whole body weight (W) was measured in grams. Age (in years) was determined from either scales or otoliths. No age data were available for NSSH in 1999. Other stocks that lacked age information included Scotia-Fundy autumn-spawning herring (SFH), Gulf of Saint Lawrence autumn-spawning herring (aGSLH), and WSH. Reproductive parameters that were indirectly calculated are marked in Table 1, but such parameters did not exist for Prince William Sound herring (PWSH). Ovary weight (OW) was determined to 0.1 g for Northeast Atlantic, local Norwegian herring stocks, and SFH, where fresh ovarian subsamples were preserved in 3.6% phosphate buffered formaldehyde (Bancroft and Stevens 1996;Ó skarsson et al. 2002;dos Santos Schmidt et al. 2017). Whole Baltic Sea herring (BSH) ovaries were preserved, and for North Sea winter-spawning herring (NSWH) oocytes were subsampled using a pipette, then preserved in buffered formaldehyde (van Damme et al. 2009;Bucholtz et al. 2013). GMH ovary samples were also preserved in buffered formaldehyde, but samples were either fresh or frozen prior to preserving (Wuenschel and Deroba 2019). Pacific herring stocks (except PWSH) and Gulf of Saint Lawrence spring-spawning herring (sGSLH) as well as aGSLH were preserved in Gilson's fluid (Messieh 1976;Bagenal 1978;. Ovaries of Sitka Sound herring (SSH) were boiled (S.C. Dressel, personal communication). The gonadosomatic index (GSI S , in %) was calculated as: GSI S = 100 9 OW/(W -OW).
Body growth analyses considered TL-at-age, fitted by the von Bertalanffy equation: TL(t) = TL ?--(1 -e -k(t-t 0 ) ), where TL(t) refers to size at age t, TL ? to the asymptotic value, t 0 to size at age zero, and k to the growth coefficient (year -1 ). Fitting of the von Bertalanffy equation and a power function were performed with the aid of R package FSA (Ogle 2016) and by standard routines, respectively. Body condition was reported as W-at-TL, where the curve was fitted using a power function of the type W = a 9 TL b , or by Fulton's condition factor K, K = 100 9 W/TL 3 . Any evidence of allometric patterns (b = 3) was specified with the listed equation (Table 2); associated plots are presented in the Supplementary Information (SI). Length dependency in K was added as additional information to the box plots (positive or negative correlation). Relative condition (K n ) was also calculated based on the W-at-TL curve for all stock combined (a = 0.003, b = 3.28, r 2 = 0.905), by the formula K n = W observed /W expected .

Oocyte size
The auto-diametric method (Thorsen and Kjesbu 2001) was used to determine the mean oocyte diameter (OD, in lm) of at least 100 developing oocytes per individual for Gulf of Maine, local Norwegian, and Northeast Atlantic herring samples, except for Icelandic summer-spawning herring (ISSH), for which 50 of these oocytes were measured manually from digital images (Ó skarsson and Taggart 2006). The difference in sample size did not influence measurement accuracy, which were tested on 10 NSSH specimens (0.416 B p B 0.976; N = 50 vs. N = 100; t-test). For BSH, oocyte diameter was estimated by the oocyte area of histological sections, corrected afterwards for shrinkage (Bucholtz et al. 2013).
Where OD data were missing, but OW and potential fecundity (F P , number of developing oocytes) data were available (Table 1), OD was given via oocyte packing density (OPD, in g -1 ). Two approaches were used to calculate the OPD. If possible, the preferred option was by using ovary subsample weights (OW S ) and associated oocyte counts (N S ), i.e., OPD = N S /OW S , if not by using the F P value and OW, i.e., OPD = F P /OW. Applying the preferred option on NSSH (dos Santos Schmidt et al. 2017) resulted in the following relationship: The validity of this relationship was then tested with BFH and Norwegian summer-autumn spawning herring (NASH), but included a stock-specific tuning factor (9.15 lm). BFH was selected due to its genetic similarity to Pacific herring (Jørstad et al. 1991), whereas NASH was selected because it spawns a few months later in the year (summer spawners) than NSSH (dos Santos Schmidt et al. 2017). Therefore, for BFH, the OD equation became: This revised BFH relationship was then, successfully applied to all Pacific herring stocks and WSH for the reasons abovementioned (Table 1). The validity of this relationship was also applied for SFH, a (summer)-autumn spawner, after including a stock-specific constant factor (15.40 lm): When information on ovarian subsample weights and accompanying oocyte counts were unavailable, OPD was then determined from the F P /OW relationship; this was done for California herring (CAH), sGSLH, and aGSLH. This approach represents an indirect estimate and was calibrated using Craig  Table 1) herring (CRH), where each individual fecundity estimate was based on five replicates. Thereafter, the OD of sGSLH, CAH, and aGSLH was determined by Eqs. 1, 2, or 3, respectively. As OPD and OD are strongly related, all ODs were considered accurate and used in all further analyses.

Egg dry weight
Egg dry weight (EDW, in lg) was reported based on analyses of hydrated oocytes of specimens belonging to Northwest Atlantic, Northeast Atlantic, and local Norwegian herring stocks (Table 1). The use of hydrated instead of postovulated oocytes (eggs) was considered adequate (Kurita et al. 2003) and reduced the risk of losing any eggs (and underestimating fecundity). For SFH and ISSH, triplicates of 10 fresh eggs were dried at 62°C for 48 h and placed in a desiccator for a further 24 h (Ó skarsson et al. 2019). The rest of the herring stocks samples had been fixed in formalin; prefixation has been shown to not influence EDW (Hempel and Blaxter 1967;dos Santos Schmidt et al. 2017). A single sample of 50 formalinfixed eggs per female were dried at 60°C for 3 days (dos Santos Schmidt et al. 2017), except for GMH, where two or three replicates contained 50 eggs were used, when the eggs dissociation was difficult a total of 25 eggs were used instead. Differences in egg dry weight procedures (Hislop and Bell 1987;Kjesbu et al. 1996) were considered but found to be insignificant for the present analyses.

Oocyte wet weight and egg water content
Mean oocyte wet weight (OoW, in mg) was given as: OoW = 1000 9 OW/F P . This estimate included ovarian stroma. The corresponding water content (in %) was: 100 -100 9 EDW/OoW, marked as STEP 1. The following estimations were restricted to NSSH, verified from earlier proximate analyses to show an egg water content of typically 70% (Kurita et al. 2003). OoW was multiplied with IF (interaction factor; IF \ 1.0); this resulting correction was applied to the 2014 NSSH data set, marked as STEP 2.

Reproductive investment
Reproductive investment (RI, in mg) was derived from multiplying egg mass by the number of eggs: RI = EDW 9 F P . This was followed by standardization by body mass, creating an index for reproductive investment (RII S ): RII S = RI/W S (mg g -1 ), where W S is ovary-free (somatic) weight (dos Santos Schmidt et al. 2017). Finally, as EDW was missing for several stocks (Table 1), OW was used as a proxy for reproductive investment.
In situ temperatures Temperature data were extracted from the NOAA World Ocean Atlas (WOA2018) (https://www.nodc. noaa.gov/access/index.html) to represent, as far as possible, the encountered thermal conditions for the various herring stocks. WOA2018 contains objectively analysed climatological fields of in situ temperatures and other physical and chemical variables at standard depths for various periods . Decadal averaged annual temperatures, from two discrete depths (0 and 30 m) were extracted for each of the 26 areas where the herring stocks are found ( Fig. 1; Table S1). The time of sampling/observation for the 26 herring stocks (Table 1) varied from the 1970s (CAH, sGSLH, aGSLH) to the recent decade for the Northeast Atlantic herring, except WSH, which was sampled in 1984. Due to low resolution of temperature data in WOA2018 from the White Sea we selected data from https://www.nodc.noaa.gov/OC5/ WH_SEA/index1.html that provide annual temperature data for the specific sampling year 1984. We assigned the decadal averaged temperature that overlaps with the biological data of the herring stocks (Table S1).
The habitat extents of herring stocks vary interannually and seasonally (e.g. Dragesund et al. 1997) but are typically not exactly known, or, in cases, indicative only. In regions with small spatial temperature gradients the differences between the tabulated temperature values and the actual ambient herring temperature can be assumed small. However, in regions where vertical and/or horizontal temperature gradients are large, the actual ambient temperatures might be quite different from the tabulated ones. Therefore, we excluded the two herring stocks in the Gulf of Saint Lawrence, aGSLH and sGSLH. Not only because of the strong vertical temperature gradients seen in this region (Table S1), causing estimates of ambient temperature to be sensitive to the depth position of the herring, but also because these stocks seasonally leave the cool Gulf of Saint Lawrence entering the warmer open ocean (Comeau et al. 2002) requiring detailed knowledge about migration routes as well. Note here that also BSH and WSH stay in waters with strong temperature gradients (Table S1).

Statistical analysis
All plots and statistical analyses were performed in R (version 3.5.2; R Core Team 2017). Biometric and reproductive parameters were tested for normality (Shapiro-Wilk test). As this assumption often did not hold, the non-parametric Kruskal-Wallis test was then used instead of parametric tests (ANOVA, t-test) when contrasting stocks. Analysis of covariance (ANCOVA) was applied to body growth (TL-at-age) and body condition (W-at-TL) data following log transformation. ANOVA was used to test difference between temperature and condition (Fulton's K and relative condition), due the differences in number of observations for herring stocks (Table 1) the weighting factor was included in the model. Post-hoc pairwise comparisons were made using Tukey's honest significant-differences (HSD) test (p [ 0.05). As the spatial configuration was known a priori (Zuur et al. 2007) (Table 1), discriminant analysis (package vegan; Oksanen et al. 2016) was used to differentiate between herring stocks based on either their biometric (TL, W, and, age) or reproductive (OD and RF P,TL 3 ) parameters, but also, as supplementary information, a combination of these two trait categories (TL, W, age, OD, and RF P,TL 3 ). A MANOVA (the Wilks lambda) test was applied to verify overall group effect significance (Zuur et al. 2007). A subset of 150 of 250 collected specimens of ISSH was randomly selected to maintain a comparable amount of data across stocks (Table 1).

Results
The physical setting of herring stocks The Northeast Pacific and North Atlantic herring are largely confined to boreal ecosystems with CAH at the upper thermal habitat range of about 13°C (with the annual mean as reference) associated with Californian upwelling ecosystem, and WSH of the lower thermal habitat range of less than 4°C associated with ecosystems close the Arctic ( Fig. 1; Table S1). The spatial extent of the stocks reflects the spatial structure of the ocean climate in the various regions: (1) the Northeast Pacific herring stocks were distributed along a 3000 km coastline from around Cape Mendocino at the Californian coast to Shelikof Strait in Alaska over which distance the annual mean ambient temperature changes from 13 to 6°C (Table S1); (2) the Northwest Atlantic herring stocks differ in spatial distribution from the other groups of herring stocks, because of the extraordinary strong latitudinal thermal gradient along the Canadian east coast where the cold southwardflowing Labrador Current encounters the warm northward-flowing Gulf Stream (Sundby 2000;Sundby and Drinkwater 2007). Here, the distributions of the herring stocks are limited from the Bay of Fundy to the northern Gulf of Saint Lawrence, a distance of only 1200 km, but where annual mean ambient temperature changes from 10 to 4°C in the upper layers (Table S1), and 3) The Northeast Atlantic stocks, including the Icelandic stock, are distributed along a 4000 km distance from the English Channel to the White Sea over which distance the ambient temperature changes from 12 to less than 4°C (Table S1).
The herring stocks are latitudinally distributed from 40°N (CAH) to 70°N (BFH) (cf. Fig. 1) implying that they all are confined to spring-bloom ecosystems with increasing seasonality in the primary production with increasing latitude, and where the increasing seasonality has a particular influence on life-cycle dynamics of planktivorous species when exceeding ''critical latitude'' near the Arctic Circle (Sundby et al. 2016). Here, primary production during winter is too low to sustain active feeding leading to dominating overwintering strategy for herbivorous zooplankton with cascading effect on planktivorous fish. These processes are relevant for the northernmost stocks of the Northeast Atlantic.

Analysis of body growth
Prespawning females of Northeast Atlantic herring, except for NSAH (North Sea autumn-spawning herring) and NSWH, were larger than local Norwegian and Northeast Pacific herring (C 26 cm vs. C 18-19 cm) (Fig. 2a, c, d). The TL-at-age curves overlapped frequently across regions (Fig. 2). Northeast Pacific herring grew from 19 to 35 cm between 2 and 13 years of age, with Kamishak Bay herring (KBH) exhibiting the fastest growth and British Columbia (BCH) and CAH the slowest (Fig. 2). The regional TL ? ranged from 28.16 (CAH) to 36.72 cm (SSH) ( Table 3). sGSLH apparently grew relatively faster than GMH (Fig. 2b), though TL ? of sGSLH was uncertain whereas settled at 29.92 cm for GMH (Table 3). Northeast Atlantic herring were 26-37 cm at 2-15 years of age, except for NSAH, in which the curve flattened out relatively quickly, noting 10 years and 31 cm as the maxima (Fig. 2c). The corresponding von Bertalanffy parameters indicated a latitudinal trend in the Northeast Atlantic herring stocks, with the northern herring stocks (ISSH, NSSH, and NASH) showing the largest TL ? and smaller growth coefficient compared to the southern herring stocks (NSAH and NWSH) ( Table 3) . Local Norwegian herring had comparable values to Northeast Pacific herring: 18-33 cm for 2-12 years of age. The examined material for Lindåspollene herring (LPH) included older fish (13-18 years) and none of the Lake Rossfjord herring (LRH) females exceeded 7 years (Fig. 2d); these two stocks span extreme situations. LRH showed the highest growth coefficient compared to other local Norwegian herring stocks, except possibly GLH (Table 3). BSH had similar length-atage as LRH at younger ages (3-6 years) but varied considerably at older ages (7 to 13 years) (Fig. 2), this large variation may be related to the lower number of observations (Table 1).
In line with the above results, both TL and age varied significantly across stocks (Kruskal-Wallis test; p \ 0.001 for both variables). In addition, the corresponding TL-at-age was different (ANCOVA; slope and intercept: p \ 0.001). For local Norwegian herring, southern stocks (Landvik herring [LVH], Coastal Skagerrak [CSH], and LPH) grew faster and became larger than those located either partly or more clearly northwards (Sognefjord [SGH], GLH, Trondheimsfjord [TRH], and BFH) (Fig. 2d). LRH, found geographically near BFH (Fig. 1), formed a third, much slower growing category (Fig. 2d). The sensitivity analysis on NSSH (2007 vs. 2014; Fig. 2c) showed interannual variation in growth (ANCOVA; slope and intercept: p \ 0.001) but was minor compared to the variation among stocks, both within and between regions (Fig. 2).

Analysis of body condition
Body condition (W-at-TL) varied among the 26 herring stocks (ANCOVA: slope and intercept: p \ 0.001) (Fig. S1), considering here also detailed a posteriori comparison (Table S2). Any detected nuance in curve appearance was also significant, as tested in NSSH (ANCOVA: slope and intercept: p \ 0.001) (Fig. S1c). W and TL were undoubtedly related (r 2 C 0.559, p \ 0.001) ( Table 2). For Northeast Pacific herring, higher W-at-TL values were found further south (Fig. S1a), indicating a latitudinal trend in body condition in an opposite direction as the one for body growth (Fig. 2a). Thus, CAH appeared to be in the best condition, followed by BCH and the investigated stocks of Alaska herring (PWSH, KBH, SSH, CRH, and Revilla Channel herring [RCH]), with the latter five showing overlapping patterns (Fig. S1a).
Overlapping W-at-TL curves were also seen among Northeast Atlantic herring stocks, except for NSSH, which had lower values (Fig. S1c). For local Norwegian herring, LPH stood out positively (p \ 0.05; Fig. S1d). The curve of WSH fell at the midpoint for local Norwegian herring (Fig. S1d). In the Northwest Atlantic, autumn spawners in the Gulf of Saint Lawrence (aGSLH) appeared stouter than the complementary spring spawners (sGSLH) (Fig. S1b). SFH, the other examined category of autumn spawners in this region, overlapped with aGSLH at similar TLs (Fig. S1b). Studying all functional relationships presently established between W and TL, the exponent b ranged from 1.61 for LRH to 3.68 for KBH (Table 2). Local Norwegian herring had the highest variation in b across stocks (Table 2). Grand mean b for all 26 stocks was 2.93. Thus, the use of Fulton's condition factor K, in which b is set at 3 (isometric growth), seemed reasonable, though many stocks had b values either significantly below 3 (8 stocks) or above 3 (8 stocks; Table 2).
The following analyses of K strengthened the above findings based on W-at-TL, but also clarified that herring condition often exhibits length dependency, either positive or negative (correlation analysis; p [ 0.05) (Fig. 3). This dependency could switch sign between seasons (cf. NSSH). Herring females from the Northeast Pacific and along the Norwegian coast up to the White Sea typically had K \ 0.8, whereas K [ 0.8 for those from the Northwest and Northeast Atlantic (Fig. 3). Notable exceptions to these regional patterns were CAH (fatter), NSSH (thinner), and LPH (fatter) (Fig. 3). Interannual variation in K was, as expected, observed; K for NSSH and NSAH fluctuated significantly (ANOVA; p \ 0.001 and p = 0.015, respectively). The other stocks tested also indicated that K is highly dynamic (p B 0.039; 9 stocks) (Fig. 3). However, these fluctuations in K for each stock were restricted and did not change the general regional and stock-specific patterns. To cancel out any bias attributed to length dependency, K n outputs were also examined finding comparable pattern as outlined for K: a negative latitudinal trend for Pacific herring and Northwest Atlantic herring (provided restricted to autumn spawners), a positive trend for Northeast Atlantic herring, and a slightly flat average for local Norwegian herring stocks (Fig. S2).

Biometric discriminant analysis
The separational effect of TL, W, and age resulted in a tilted U-shaped plot in which smaller-sized Northeast Pacific, local Norwegian and Baltic Sea herring stocks clustered, whereas larger-sized local Norwegian herring clustered with the remaining pool of oceanic stocks (Fig. 4). The sensitivity test of NSSH (2007 vs. 2014) did not alter this impression (Fig. 4). The first principle component explained 69.1%, and the second trace, 23.5%. A significant group effect was also observed (MANOVA-the Wilks lambda test; p \ 0.001).

3-
[ 200 cm -3 and OD \ 1100 lm, i.e., were located in Q4 (Fig. 5). Deviations from this pattern included NSAH (in Q2 instead of Q4) and both BSH and sGSLH (in Q4 instead of Q2). CAH and WSH were outside of this pattern, in Q1 and Q3, respectively (Fig. 5). For WSH, this may be due to high uncertainty because of few observations (Table 1). Higher withinstock variability (± 95% CI) was generally observed for OD than RF P,TL 3 , especially for spring spawners (Fig. 5). The existence of interannual variability between OD and RF P,TL 3 was exemplified for NSSH, but these three sets all fell within the same quadrant, Q2 (Fig. 5). Finally, trade-offs by season were tested at the individual level (OD vs. W and RF P,TL 3 vs. W), often with a large spread in data points within a given stock, including crossings of the above-mentioned threshold values (Figs. S3-S4), even in cases when the overall trade-off was within the expected quadrant (e.g., ISSH; Fig. 5).

Reproductive discriminant analysis
The selected set of reproductive variables, OD and RF P,TL 3 , resulted in a U-shaped scatterplot (Fig. 6). Overall, Northeast Pacific and local Norwegian herring stocks, jointly labelled as ''local'' (Table 1), formed one cluster, whereas the oceanic stocks in the Northwest and Northeast Atlantic formed another. However, NSSH and NSAH, both oceanic stocks, fell into the former category, whereas CAH, a local stock, appeared in-between (Fig. 6). BSH was isolated from all the other herring stocks. For NSSH, all three plotted points fell within the same area of the graph (Fig. 6). The first trace explained 87.5% of the variation, and the second trace, 12.5%, a significant group effect was recorded (MANOVA-the Wilks lambda test; p \ 0.001).

Biometric-reproductive discriminant analysis
This combined analysis (Fig. S5) largely supported the above two separate analyses of a similar kind: Pacific and local Norwegian herring clustered, except their larger body-size members, KBH, LPH, CSH, and LVH (Fig. 2), which joined with NSAH and NSSH (2007). GMH, sGSLH, ISSH, and NASH formed a third cluster, whereas NSSH (2014) appeared as an Fig. 3 Fulton's condition factor K for prespawning females across the examined herring stocks (Table 1), sorted by geographic region. For each box plot (main study), the thick line is the median value, top and bottom lines indicate the 75th and 25th percentiles, respectively, and whiskers indicate ± 2 SD. Seasonal variation is indicated for nine stocks (two additional years for BCH, NSAH, and NSSH instead of one), with the filled circle indicating the median value (with ± 2 SD whiskers). TL dependency, positive (?) or negative (-), is indicated as *p \ 0.05, **p \ 0.01, and ***p \ 0.001. K = 0.8 is defined as a general reference line, whereas K = 0.7 represents the threshold between poor and good body condition of NSSH; K \ 0.7 is associated with intensified atresia (Ó skarsson et al. 2002) isolated spot in the plot, as well as BSH (Fig. S5). The first trace explained 49.7%, and the second trace, 28.3%, a significant group effect was recorded (MANOVA -the Wilks lambda test; p \ 0.001).

Oocyte and egg characteristics
OoW (oocyte wet weight) and OD were closely related, with NSAH possibly slightly off compared to the other stocks (Fig. S6a). Although EDW and OD also showed a strong positive relationship, at least for OD [ 1100 lm, SFH with its the smaller ODs (700-1100 lm) apparently followed a different trajectory for EDW (Fig. S7). Furthermore, EDW of NSAH hardly increased with OD. Due to many instances of little spread along the x axis regression analyses were not performed (Fig. S7). Water content estimates steadily increased with OD for SFH, whereas these were established at a higher plateau for the other stocks (Fig. S6b). All these presented water content figures were, however, inflated by the inclusion of ovarian stroma. When the ovarian stroma were excluded the estimates fell from an average of 81.6 (STEP 1) to 70% (STEP 2, published observations) in the case of NSSH, or, as determined by iteration, to be &15% of OoW, i.e., IF = 0.15 (see Material and methods), in this stock at this prespawning stage.

Reproductive investment
In terms of absolute values, oceanic herring invested generally more in reproduction than local herring, regardless of spawning season (Fig. 7a). Following standardization by somatic body weight, this difference decreased but was still significant (Kruskal-Wallis test, p \ 0.001) (Fig. 7b).

Relationship between fecundity, ovary size and body growth
Overall, both F P and OW increased in synchrony with body growth, but the similarity in this across-stock Fig. 4 Discriminant analysis scatterplot combining biometric (total length and whole body weight) and age data for herring stocks (Table 1), grouped by region: 1 = Northeast Pacific; 2 = Northwest Atlantic; 3 = Northeast Atlantic; 4 = Local Norwegian; and 5 = Baltic Sea. The ellipsoid indicates the 90% tolerance interval, and the symbols represent the average value per stock. Missing age information for SFH and aGSLH (both in Region 2), and WSH (Region 6) precluded inclusion of these stocks. For NSSH, the output for 2014 (default) is compared to the output for 2007. Within each region, stocks are sorted by general herring stock location, from north (top) to south (except for the supplementary year for NSSH) Fig. 5 Trade-off between mean oocyte diameter (OD) (± 95% CI) and relative length-based fecundity (RF P,TL 3 ) (± 95% CI) in prespawning individuals from different herring stocks (Table 1) located in the following geographic regions: 1 = Northeast Pacific; 2 = Northwest Atlantic; 3 = Northeast Atlantic; 4 = Local Norwegian; 5 = Baltic Sea; and 6 = White Sea. Lines at 1100 lm and 200 g -1 were used to divide the plot area into four quadrants (Q1-4), with most spring spawners occurring in Q2 and summer/autumn spawners in Q4. No reproductive data were available for PWSH (Region 1). For NSSH, the default analysis (2014) was extended to include two additional years, 1999 and 2007. Within each region, stocks are sorted by main place attachment, from north (top) to south (except for the supplementary years for NSSH) trend was more evident for OW (Fig. 8). As the herring grew from age 4 to 6 years, the relationship turned steeper, again more clearly seen with OW as the response variable. Focusing on the latter, TRH exemplified a ''catching up'' trend, while the observation for the CAH set became progressively more extreme.

Relationship between body condition and in situ temperature
Mean prespawning K and annual mean temperature showed a weak positive relationship (ANOVA, r = 0.298; p \ 0.001), where CAH was isolated, with the highest condition at the highest annual mean temperature, followed by most of Northwest and Northeast Atlantic herring stocks (except NSSH) but also LPH (cf. Figure 3). K was typically around 0.9 in these cases referring to a temperature range of 7.5 to 12.5°C (Fig. 9a). The remaining herring stocks showed mean K B 0.8 for 4-10.5°C. Mean K n and annual mean temperature showed a flatter relationship but still significant (ANOVA, r = 0.274; p \ 0.001), except for CAH that showed an extremely high value compared to the other herring stocks (Fig. 9b). The two GSLH stocks were not included in these relationships for reasons outlined in the Material and methods.

Discussion
This analysis of 26 herring stocks had three main functions: (i) the demonstration of possibly several unique stock-specific traits, (ii) illustrate the potential utility of the present work as a platform for further studies, and (iii) as a useful reference for monitoring effects of on-going climate change (IPCC 2018) on herring trait expressions. Due to the rich variety of reproductively isolated stocks for both Clupea pallasii and C. harengus (Hay et al. 2001;Geffen 2009), this review and meta-analysis required extensive co-operation across laboratories to adequately contrast reproductive performance and corresponding body metrics of herring located in different and distant waters. Thus, this work was made possible by a compilation of available national data sets and by acquiring new information on several herring stocks. Our work is special in the herring literature because of the comprehensive list of variables considered along the axes of both reproduction and biometrics, replacing historic laboratory techniques within the field of reproductive analyses with modern techniques, and, not at least, examination of the material for individual variability rather than addressing only pooled data. To the best of our knowledge, this study is the first thorough comparison of Pacific and Atlantic herring reproductive traits. Most of the earlier, relevant information seems to be restricted to either Pacific (e.g. Tanasichuk and Ware 1987;Hay et al. 2008) or Atlantic herring stocks (Parrish and Saville 1965;Blaxter 1985). The life-history of the 26 herring stocks refer to different areas and migratory styles, i.e., from a relatively stationary (local) to highly migratory (oceanic) behaviour. We also considered herring stocks distributed over extensive latitudes (from approximately 40 to 70°N), with different day lengths (Sundby et al. 2016). Furthermore, both unproductive habitats (e.g., the high-latitude sea loch of LRH) (Hognestad 1994;Mikkelsen et al. 2018) and highly productive habitats (e.g., the temperate eastern  (Table 1) grouped by region: 1 = Northeast Pacific; 2 = Northwest Atlantic; 3 = Northeast Atlantic; 4 = Local Norwegian; 5 = Baltic Sea; and 6 = White Sea. The ellipsoid indicates a 90% tolerance interval, and symbols represent the average value per stock. PWSH (Region 1) is not represented. For NSSH, the output for 2014 (default) is compared to the outputs for 1999 and 2007. Within each region, stocks are sorted by main place attachment, from north (top) to south (except for the supplementary years for NSSH) boundary upwelling zone (Field et al. 1998) inhabited by CAH) were included. So, our results should provide a representative overview of the scope of reproductive trade-offs for these ecologically important clupeids, including several stocks of high economic value (Hay et al. 2001). Though it is beyond the scope of this work to elaborate on phylogenetics, both LRH and BFH were identified as ''C. pallasii peripheral populations'' in Libungan et al. (2016). In Mikkelsen et al. (2018), LRH and BFH were identified as C. pallasii rather than C. harengus. However, because the other local Norwegian herring stocks considered should all be considered C. harengus (Martinez Barrio et al. 2016), we opted for geographic region as one of the separation criteria instead of genetics.
A limitation of our work is the general lack of time series data, especially regarding complete sets of reproductive traits of current interests, although we applied sensitivity analyses attempting to address this concern. However, we understand that significant lag effects might be important requiring additional complex data sets, but they were inaccessible across most of stocks in our pilot studies. For example, the fecundity of NSSH is not only a function of prey availability and thereby body condition during the current season, but also during previous seasons (i.e., at oogonial proliferation) (dos Santos Schmidt et al. 2017). Moreover, environmental conditions among different regions may not fluctuate synchronously, e.g., a positive North Atlantic Oscillation Index (NAOI) creates a relatively warmer situation in the northeast Atlantic but a relatively cooler situation in the northwest Atlantic (Bjerknes 1962;Sundby and Drinkwater 2007). A negative NAOI will reverse this situation. Over the last 50 years, Northeast Atlantic herring stocks have been exposed to warmer temperatures caused by positive phase of another mode of natural climate variability, the larger-scale Atlantic  (Table 1), sorted by spawning season. Stocks are further distinguished based on their distribution: oceanic versus local (grey box). For each box plot, the thick line is the median value, top and bottom lines indicate the 75th and 25th percentiles, respectively, and whiskers indicate ± 2 SD. Stocks in the Northeast Pacific are not represented due to lack of egg dry weight data Multidecadal Oscillation (AMO; Drinkwater et al. 2014;Hoegh-Guldberg et al. 2014) that has influenced the entire North Atlantic synchronously. Temperature also varies in these regions based on different water currents and/or local habitat condition where each herring stock lives (see Material and Methods). Therefore, herring collected on each side of the Atlantic should not be uncritically compared regarding seasonal attributes, such as body condition (Sande et al. 2019). Along the same lines, Pacific herring can be affected by the El Niño Southern Oscillation (Bailey et al. 1995;Lehodey 2001), and the largerscale and interannual to multidecadal time scale, Pacific Decadal Oscillation (PDO) (Mantua and Hare 2002) with particular impact in the North Pacific north of 20°N with opposite phases in the Northwest and Northeast Pacific. PDO has been ascribed to influence Northeast Pacific herring stocks, such as stock size of BCH (Landis et al. 2003). Additionally, within the herring habitats of the North Pacific and North Atlantic there are recent differences in temperature trends. In the California Current (habitat of CAH) and in the Gulf of Saint Lawrence (habitat of sGSLH and aGSLH) the long-term temperature has remained nearly constant, while in the northern parts of the Northeast Pacific (habitat of BCH, RCH, CRH, SSH, KBH, and PWSH) the temperature has increased, as noticed in the above-mentioned climate literature. All of the Northeast Atlantic herring stocks have experienced a substantial temperature increase since 1970s. This is caused by amplification of the combined effects of the positive phase of AMO and the global anthropogenic climate change (Hoegh-Guldberg et al. 2014). These changes are clearly reflected in the three Fig. 8 Relationship between mean body growth (± 95% CI) and mean potential fecundity (± 95% CI) (a, c, and e) and mean ovary weight (± 95% CI) (b, d, and f), split by age class. The herring stocks are organized by region: 1 = Northeast Pacific; 2 = Northwest Atlantic; 3 = Northeast Atlantic; 4 = Local Norwegian; 5 = Baltic Sea; and 6 = White Sea. PWSH, aGSLH, SFH, and WSH were excluded due to lack of fecundity or age data temperature time intervals presented in Table S1. The largest sea surface temperature (SST) increase has been experienced by BSH and WSH. This is linked to the larger thermal stratification in these two regions. Additionally, sGSLH and aGSLH are also living in highly stratified waters; the two Gulf of Saint Lawrence stocks are subjected to special thermal conditions caused by inflow of the so-called Cool Intermediate Layer (CIL) from the Labrador Current with its core vertically extending between 30 and 100 m depth (Gilbert and Pettigrew 1997). Here, temperature at 30 m depth was about B 3°C during the sampling of the stocks in the 1970s and at 100 m depth only 2°C as an annual average. The CIL is more dominant in the northernmost part of the Gulf of Saint Lawrence. Hence, the spring-spawning stock (sGSLH) is most influenced of the two. It is uncertain if the herring would avoid such depths due to the low temperature. Hopefully, our work will encourage lengthy sampling programs in both oceans to assess differences in environmental impacts on herring reproductive ecology in a more comprehensive manner. Body growth was currently used as a metric to reflect environmental influences (see Sande et al. 2019 and references therein), as addressed below.
Distinct differences in biometric and reproductive features were identified when compared across all herring stocks. Northeast Pacific and local Norwegian herring stocks often exhibited similarities, whereas oceanic Atlantic herring stocks tended to form another separate subset. Thus, the noted patterns often referred to geographic location (or region) rather than species (C. harengus vs. C. pallasii), emphasizing the importance of local adaptations. Northeast Pacific herring showed a positive effect of latitude on body size, which was well-documented previously (Tanasichuk et al. 1993;Hay et al. 2008;Ito et al. 2015). This northward trend seems to also occur among Northeast Atlantic herring stocks, but it was inversed for local Norwegian herring. Hay et al. (2008) speculated that this latitudinal gradient in body size may be related to environmental conditions in each living area. Generally, body growth (length-at-age) and condition (weight-at-length, K or K n ) varied between oceanic and local herring stocks; oceanic herring typically had faster growth rates and stouter body shape. Parrish and Saville (1965) extensively compared the morphological features of oceanic and ''shelf'' stocks, with the latter represented by different North Sea herring stocks. They summarized oceanic herring as having a long-life span, large body size, and late maturation, whereas coastal/local herring were completely opposite. The generally larger body size of oceanic herring is probably attributable to their longer migration Fig. 9 Relationship between body condition factor, either Fulton's condition (a) or relative condition (b), and annual mean temperature from 1955 to 2017. aGSLH and sGSLH are not included due to strong vertical thermal stratification (Table S1) along with complex migration routes (see Material and methods). The herring stocks are organized by region: 1 = Northeast Pacific; 2 = Northwest Atlantic; 3 = Northeast Atlantic; 4 = Local Norwegian; 5 = Baltic Sea; and 6 = White Sea routes (Dragesund et al. 1997). For example, in recent years, NSSH have migrated from feeding grounds in the Norwegian Sea to overwintering areas outside northern Norway before moving on again to the main spawning grounds further south along the coast off western Norway, an annual migration distance [ 2700 km (Huse et al. 2010). In contrast, local Norwegian stocks are mainly restricted to the fjords (Lie et al. 1978;Hognestad, 1994). LRH migrate in an irregular fashion, without forming any school, inside the small 12 km long and 1 km wide lake (Hognestad 1994). Pacific herring also exhibit variability in migration distance. The prime example of long migration in these waters is Togiak herring in the Alaskan Bering Sea, which migrate into inshore areas (Bristol Bay) during the spawning season but stay offshore (Unimak Pass and between Pribilof Islands and St. Matthew Island) during the overwintering and feeding period, migrating * 2100 km over this time period (Tojo et al. 2007). At the other end of the scale in these waters are BCH, a herring of smaller body size, which spawn in the Strait of Georgia and feed in shelf regions off the west coast of Vancouver Island, covering a total distance of 400 km (Hay et al. 2001). Not surprisingly, herring migration distance is related to body size; larger fish will spend less metabolic energy during migration than smaller fish (Slotte 1999;Nøttestad et al. 1999). Thus, larger herring can utilize distant productive areas, which should further support faster growth and greater reproductive investment. There are exceptions to these generalizations; the local Norwegian herring stocks LPH, CSH, and LVH grew roughly in the same way as the oceanic stocks. Their large size has been attributed to cooccurrence with NSSH (Silva et al. 2013) and other coastal herring stocks (e.g., Skagerrak, and Kattegat herring) (Johannessen et al. 2009Eggers et al. 2014;Berg et al. 2017). More specifically, LPH have interacted with NSSH over prolonged periods, but mainly during the NSSH stock collapse in the late 1960s when its distribution was highly restricted over a couple of decades to near the Norwegian coast ). Here we found that autumn spawners tend to have higher body condition, but lower body growth compared to spring spawners, agreeing with findings in Berg et al. (2020) studying a herring fjord population complex near Bergen, Norway. The higher prespawning body condition of autumn spawners may be explained by that the feeding season to a larger extent overlap with the main part of gametogenesis reducing the depletion of surplus energy compared to spring spawners displaying fast gonad growth during late autumn and early winter (Kurita et al. 2003;dos Santos Schmidt et al. 2017) (see below). Recently, it has been demonstrated that the combination of three approaches, such as maturity stage, otoliths microstructure and, genetics analysis based on single nucleotide polymorphisms can be a better tool to distinguish between spring and autumn spawners (Berg et al. 2020). Although, most of the herring stocks studied here were spring spawners (n = 19) and samples were collected prior to spawning which unequivocally allows for the correct assessment of spawning seasons, the applicability of these three approaches combined could be an effective tool on sympatric stocks, such as Gulf of Saint Lawrence herring.
Apart from intrinsic factors, extrinsic factors may also influence the body condition of herring. Although the relationship with ambient temperature across the entire thermal range appears highly scattered, there is, statistically speaking, a positive correlation. One may speculate what is the mechanistic link between temperature and the consulted measures of body condition (K and K n ). Most probably, temperature is a proxy for other extrinsic factors such as integrated primary production, which might have some latitudinal variation and which in turn is affected by ocean temperature. Moreover, spatial changes in primary production will cascade into the abundance of zooplankton food abundance.
Earlier studies indicate that a trade-off between fecundity and egg size is to be expected among herring stocks, at least for Atlantic herring (Blaxter and Holliday 1963;Hempel and Blaxter 1967;Blaxter 1985). The most evident explanation is the timing of egg production, spring versus autumn spawners, in relation to peaks in prey abundance for the startfeeding larvae, which has been studied extensively (Cushing 1967;Hempel and Blaxter 1967;van Damme et al. 2009;dos Santos Schmidt et al. 2017). We systematized the available reproductive information by splitting the RF P,TL 3 -OD diagram into four quadrants using RF P,TL 3 = 200 oocytes cm -3 and OD = 1100 lm as thresholds. In general, summer and autumn spawners produce many smaller eggs (Q4), whereas spring and winter spawners produce relatively fewer, but large eggs (Q2) (Parrish and Saville 1965;Hempel and Blaxter 1967;dos Santos Schmidt et al. 2017). We found that Pacific herring (spring spawners) fit into this scheme, though not the southernmost stock located in an upwelling area, the CAH, which is apparently energetically able to produce many large eggs. Local Norwegian herring (spring spawners) also fit, though not WSH, perhaps due to the extreme cold fjord habitat (Pesciaroli et al. 2012). Fundamentally speaking, these differences in reproductive allocation patterns (i.e., egg size vs. fecundity) are regulated by the length of the reproductive cycle: winter and spring spawners with a longer period of vitellogenesis (extended yolk uptake), whereas autumn spawners have a shorter period (Parrish and Saville 1965;van Damme et al. 2009;dos Santos Schmidt et al. 2017) (see also above). However, a lengthened reproductive cycle comes at the cost of increased energy expenditure and higher prevalence of atresia (Kurita et al. 2003;van Damme et al. 2009van Damme et al. ). Ó skarsson et al. (2002 found that EDW and OD are tightly coupled in NSSH. However, this link varies among stocks; for example, for NSAH we found indications that oocytes may grow without any further increase in EDW (arrested yolk uptake). The underlying reason is evasive, as NSAH fell in Q2 rather than Q4 as expected in the RF P,TL 3 -OD diagram. Another pattern in the associated EDW-OD plot was that SFH stayed in the so-called atretic window (800 \ OD \ 1000 lm) towards spawning, i.e., past this window, energetic costs markedly increase due to the exponential increase in OD, and therefore EDW (Kurita et al. 2003). Intensified atresia is thought to provide energy (or vitellogenin directly) to neighbouring developing oocytes (Kurita et al. 2003). The fecundity of SFH was extremely high, i.e., at the same level as the much larger aGSLH adults, which are also found in the same region, indicating that atresia was not a major issue in SFH females. Apparently, both SFH and NSAH have opted for different reproductive tactics to lessen the extra cost associated with accelerated vitellogenic sequestration. Complementing estimations of the overall reproductive investment standardized for somatic body size (RII S ) clarified that these figures were much higher in SFH than NSAH. Unfortunately, EDW data were not available for Pacific herring, but we found that local herring clearly invest less in reproduction than their oceanic counterparts, a finding that is in line with the above elaborations on body size-related issues in less and more productive habitats.
The use of ovary weight (OW) as a proxy for reproductive investment turned out successful when related to body growth, the latter taken as a spatiotemporal integrator of ambient environmental conditions (Sande et al. 2019). First, OW was reported for all stocks except one. Likewise, body growth at age 4, 5, and 6 years could be provided for most stocks. An unexpectedly tight relationship between OW and body growth occurred across stocks, underlining similar fundamental energy allocation patterns for both C. harengus and C. pallasii, and some of the variation is likely explained by spawning season: spring versus summer-autumn spawners. These plots also clarified that older specimens progressively invest relatively more in reproduction. Logically, the relationship between F P and body growth was less clear due to the outlined phenomenon of well-developed trade-offs between F P and OD (EDW) implying that F P only partly explains OW dynamics. Further to this, the related exercise of estimating oocyte wet weight (OoW) resulted in an exceedingly strong physical relationship between OoW and OD, not only telling that F P and OW were recorded correctly across laboratories, but also that the relative amount of ovarian stroma appears consistent. The indicated level of 15% stroma for NSSH is, however, higher than expected;  state this to be\ 5%. This topic should be followed up by measurements of spent ovary weights  but also stereological examinations (Serrat et al. 2019).
Our approach was extensive, studying 26 herring stocks in different waters in the Northern hemisphere. Newer data were compiled and analysed, split by stock. Prespawning Pacific and Atlantic herring, as well as local Norwegian herring, were compared methodically in terms of reproductive traits, with local Norwegian herring largely exhibiting a similar pattern as Pacific herring, whereas Atlantic herring deviated. Plots on EDW as a function of OD revealed unexpected patterns (e.g., less or arrested accumulation of egg mass, likely to save energy). Overall, the analysis demonstrated a significant degree of plasticity in reproductive traits across stocks in association with great variability both in body condition and growth. This kind of plasticity is an important mechanism for herring as a species to respond to changing environmental conditions, including seasonal pulses of energy that may or may not coincide with oocyte development. Nevertheless, body growth and reproductive investment seem to be closely coupled for all herring stocks examined speaking for general energy allocation principles. However higher efficiency of direct energy allocation to ovary (offsetting metabolic costs of energy transfer and storage/maintenance) may allow summer-autumn spawners to allocate more energy to reproduction than spring spawners. Our study does not provide any evidence that the present extent of climate change (IPCC 2018) has affected herring reproductive trade-offs. Therefore, we would recommend the continuation or establishment of dedicated reproductive ecology time series, in parallel with data collection as part of stock assessment because the reproductive potential (status) is an integral component of the following stock advice provided.