Mariana Álvarez-Noriega,Craig R. White,Jan Kozłowski,Troy Day,Dustin J. Marshall
Life history optimisation drives latitudinal gradients and responses to global change in marine fishes
Within many species, and particularly fish, fecundity does not scale with mass linearly; instead, it scales disproportionately. Disproportionate intraspecific size–reproduction relationships contradict most theories of biological growth and present challenges for the management of biological systems. Yet the drivers of reproductive scaling remain obscure and systematic predictors of how and why reproduction scaling varies are lacking. Here, we parameterise life history optimisation model to predict global patterns in the life histories of marine fishes. Our model predict latitudinal trends in life histories: Polar fish should reproduce at a later age and show steeper reproductive scaling than tropical fish. We tested and confirmed these predictions using a new, global dataset of marine fish life histories, demonstrating that the risks of mortality shape maturation and reproductive scaling. Our model also predicts that global warming will profoundly reshape fish life histories, favouring earlier reproduction, smaller body sizes, and lower mass-specific reproductive outputs, with worrying consequences for population persistence.
Life history optimisation drives latitudinal gradients and responses to global change in marine fishes

Within many species, and particularly fish, fecundity does not scale with mass linearly; instead, it scales disproportionately. Disproportionate intraspecific size–reproduction relationships contradict most theories of biological growth and present challenges for the management of biological systems. Yet the drivers of reproductive scaling remain obscure and systematic predictors of how and why reproduction scaling varies are lacking. Here, we parameterise life history optimisation model to predict global patterns in the life histories of marine fishes. Our model predict latitudinal trends in life histories: Polar fish should reproduce at a later age and show steeper reproductive scaling than tropical fish. We tested and confirmed these predictions using a new, global dataset of marine fish life histories, demonstrating that the risks of mortality shape maturation and reproductive scaling. Our model also predicts that global warming will profoundly reshape fish life histories, favouring earlier reproduction, smaller body sizes, and lower mass-specific reproductive outputs, with worrying consequences for population persistence.

Reproductive allometry occurs when reproductive output increases disproportionately with body size. Reproductive allometry (as opposed to isometry) has recently been shown to be ubiquitous in marine fishes for batch fecundity [1], but it can also occur in other taxa ([2,3]; but see [4]), presenting challenges to our basic understanding of biological growth, how we manage disease vectors [3], and the way we use biological resources [5]. Most growth and life history models assumed that reproduction scales linearly with body size and failed to predict the existence of reproductive hyperallometry [1,2,6,7]. From a practical perspective, failing to account for reproductive hyperallometry risks underestimating the contribution of larger females to population replenishment [1,8] and results in massive, systematic errors in calculations of sustainable harvests [5]. One branch of theory–life history optimisation models [9,10] anticipated reproductive hyperallometry, but this prediction has largely been overlooked. In these models, the ideal resource allocation strategy is the one that maximises fitness for a given set of environmental conditions (e.g., for a given mortality regime). Like all models, life history optimisation models are a simplified representation of reality, and they use phenomenological descriptions of the mechanisms underlying the physiological constraints, leading to different energy intakes, usages, and allocation patterns. However, this simplification allows the models to focus on the determination of the best strategy of resource allocation for a given set of ecological conditions and physiological constraints [11]. While many studies investigating broad-scale biological patterns have focused on the effects of physiological constraints [12–14], far fewer have examined these patterns from an evolutionary optimisation perspective. Yet our capacity to predict how organisms will respond to human-mediated change depends critically on our knowledge of how organisms evolve and adapt to these conditions. Life history models predict that reproductive scaling should be a product of mortality regimes and growth rates. Here, we conduct such a test for marine fishes—a group for which there are important practical reasons for understanding and predicting reproductive scaling [1,5,8]. Because fish life histories are known to vary along a latitudinal gradient [15,16], we gathered information on biogeographical variation in the key parameters that shape predictions of classic life history optimisation models for fish across latitudes: offspring size (n = 52, 47 species), rates of growth (n = 52, 47 species), and adult mortality (n = 85, 77 species). We then used these inputs to parameterise Day and Taylor’s [9] life history evolution model to make predictions about how both reproductive schedules and reproductive scaling should covary with latitude. Finally, we tested these predictions by gathering data on age at 50% maturity (A50, n = 127, 59 species) and fecundity–mass relationships (n = 13,760, 153 species) across the globe.

Natural mortality (M) (expressed as a size-independent mortality rate that is not a product of fishing) is highest at the tropics and decreases by 80% closer to the poles (from 0.85 y−1 at 0° to 0.17 y−1 at 60°), which constitutes nearly a 2-fold increase in yearly survival probabilities (from 0.43 at 0° to 0.84 at 60°, 95% CI 1.33- to 3.69-fold) (R2 = 0.16) (S1A Fig). The effect of latitude on mass-specific growth rate (k) could not be differentiated from zero (the slope of ln(k) with latitude had an estimate of −0.01, 95% CI: −0.03 to 0.01; S1 Table), although the tendency was for k to decrease with increasing latitude (from 4.36 g1/4 y−1 at 0° to 2.11 g1/4 y−1 at 60°) (R2 = 0.13) (S1B Fig). We did not detect an effect of latitude on the theoretical initial offspring weight (w0) (similar to t0 in the von Bertalanffy equation) (the slope of ln(w0) with latitude had an estimate of 0.00, 95% CI: −0.5 to 0.2), with the median expectation for w0 ranging between 12.9 and 14.4g (R2 = 0.14) (S1C Fig). Using these estimated relationships (with their attendant uncertainty), we calibrated the life history model to predict optimal reproductive schedules (see Methods).

Age at maturity was predicted to increase from the tropics to the poles (reproduce in the first year at 0° and 12.1 y at 60°; Fig 1A). Reproductive scaling was predicted to triple across the same latitudes (2.8 at 0° to 8.7 at 60°; Fig 2A). Consequently, reproductive scaling and age at maturity were predicted to positively covary with each other. Similar predictions arise under the assumption of size-dependent mortality (S2 Fig).

We found qualitative support for the predictions of our model when we compiled the available data from the literature. Age at maturity is difficult to measure, but a proxy for maturity, age at which 50% of individuals become mature (A50), increased 7.7-fold across the same range of latitudes (0.76 y at 0° to 5.87 at 60°; R2 = 0.84, 95% CI: 0.78 to 0.88) (Fig 1B). When we examined the available batch fecundity data for 1,000s of individuals in over 100 species, we find support for the model predictions. Intraspecific reproductive scaling—the slope of the relationship between fecundity and mass (in log–log scale) for individuals of the same species differing in size—is always greater than 1 and increases with latitude (13% increase from 0° to 60°; Fig 2B). For example, a 10-kg fish is predicted to produce approximately 2.2 million eggs at 0° latitude (0.5 to 9.2 million, 95% CI) and approximately 3.5 million at 60° latitude (0.9 to 14.4 million, 95% CI), a 20-kg fish is predicted to produce approximately 4.7 million eggs at 0° latitude (1.2 to 19.1 million, 95% CI) and approximately 8.1 million at 60° latitude (2.2 to 14.4 million, 95% CI). The number of eggs per female is higher for high-latitude species relative to low-latitude species across all female sizes except for the smallest females, where fecundities are essentially equivalent (Fig 2C). Because only a few species had data across a range of latitudes, different species were analysed at different latitudes, and, therefore, the dependence of within-species reproductive scaling on latitude must be treated as interspecific in Fig 2. For those species where mass and fecundity data were available across at least 20 degrees latitude, reproductive scaling also increased with latitude within species (S3 Fig). In contrast, when the relationship between mass and fecundity is estimated by taking a snapshot of the within-species relationship (e.g., only mass and fecundity at maturity instead of fecundity across all mature sizes within a population) and comparing this snapshot across species, such interspecific relationship is typically hypoallometric [6,17]. Although the direction of the relationship between reproductive scaling and latitude is successfully predicted by life history optimisation, the values of reproductive scaling estimated from the data are markedly lower than those predicted by the life history model. We suspect this discrepancy arises for 2 reasons: (i) how reproductive output is estimated; and (ii) how reproductive costs are modelled. First, our model makes predictions about the scaling of total reproductive output while we test these predictions based on the available data, which includes only 1 component of reproduction (batch fecundity). As well as having higher fecundity, larger females also tend to spawn more frequently and produce larger, heavier eggs than smaller females [18], so the scaling of total reproductive output (which incorporates the effects of spawning frequency, egg size, and batch fecundity) will be steeper than estimates based on batch fecundity alone. For example, Sardiops sagax shows a batch fecundity scaling of only approximately 1, but the scaling of total reproductive output in this species is approximately 3.5 because larger females spawn much more frequently than smaller females [18]. Hence, more comprehensive estimates of the scaling of total reproductive output are likely to be closer to those predicted by our model. Second, our model assumes all of the energy that goes into reproduction manifests as reproductive output, but in reality, a nontrivial proportion of reproductive allocation is likely to be expended as “overhead.” While the overheads costs of reproduction are poorly resolved empirically, initial theoretical explorations show that including such costs would lower our predicted estimates of reproductive scaling and bring them closer to those observed in nature (S4 Fig). Empirically estimating how energy allocated into reproduction translates to reproductive output across different body sizes is necessary for better quantitative predictions of size–fecundity relationships. While we acknowledge that our predictions differed quantitatively from the observed patterns, it is striking that a simple life history optimisation model [9], parameterised only with data on mortality, somatic growth, and offspring size across latitudes, can successfully predict qualitative global patterns in maturation and reproductive scaling for marine fish. Indeed, even when the effect of offspring size is neglected, the model still recovers the direction of change in age at maturity and reproductive scaling across latitudes. Broadly, and in relative terms, tropical fish suffer high mortality and maximise their fitness by diverting energy into reproduction earlier in life, reaching smaller sizes, having only shallow reproductive scaling exponents (though they are still greater than 1), and producing fewer eggs per unit mass. In contrast, polar fish experience lower mortality, and, therefore, their optimal strategy is to continue growing until later in life, reaching a larger size, and having steep reproductive scaling. In other words, in a fast–slow continuum, low-latitude fish tend to have fast life histories, whereas high-latitude fish tend towards the slow-paced end of the spectrum. These results contradict the expectation that highly seasonal environments should select for faster life history strategies, while in stable environments—such as the tropics—it would be beneficial to tend towards a slow-paced demography [19]. For example, low-latitude birds tend to show higher parental care and delayed maturation than temperate birds [20,21]. This discrepancy with our findings may be explained by patterns in mortality: in birds, mortality increases with latitude, as many die when temperatures drop in winter [22]. Our findings, in conjunction with others [11], imply that many of the life history patterns we observe in nature are driven, not only by mechanistic constraints, but also by the evolutionary optimisation of life histories. At the very least, life history optimisation approaches may be a useful tool for explaining patterns in life history based only on a few key parameters. Hyperallometric scaling occurs because, as an increasing proportion of energy is diverted from growth to reproduction with time, the increments of female size decrease with age while the increments of reproductive output increase. Our model makes the novel prediction that fishes delaying reproduction because of low mortality produces a higher reproductive scaling. This pattern arises because the delay in the onset of reproduction yields larger sizes at maturity, and a shallower rate of increase in the size–production relationship whenever the scaling of production is <1 (Fig 3A and 3B). Hence, when a fish that has delayed reproduction finally starts allocating resources into reproduction, proportional changes in size between time-steps are small. If energy allocation from growth to reproduction changes with age rather than size, the amount of production allocated to reproduction increases more rapidly with size than for a fish that matures early. For example, for a fish in the tropics, the total amount of production allocated into reproduction increases slowly across sizes (Fig 3A). In contrast, for a fish at high latitudes, while the postmaturity change in energy allocation with time is the same, the change with respect to size occurs much more abruptly (Fig 3B). Hence, the rate of increase in reproduction is much greater in a relative sense than the rate of growth compared to fish that mature early (Fig 3C). This means that reproductive scaling should be steeper in organisms that delay reproduction for longer—the compiled data support this prediction: higher-latitude species have later reproduction and higher reproductive scaling. Importantly, our model makes no predictions about how reproductive allocation should be divided in terms of offspring size and number though other models make such predictions [23]. Given that higher-latitude species tend to produce larger offspring, it would be interesting to see whether both fecundity and offspring size show steeper scaling at higher latitudes. A steeper fecundity scaling in cold temperatures also occurs in taxa other than fish (e.g., flies [24], water fleas [25], and water snails [26]), but with exceptions (e.g., mosquitoes [27] and water striders [28]). A previous life history optimisation model from Arendt [29] argues that the size–temperature relationship, where organisms living in colder environments reach larger sizes than organisms in warm areas, can be explained by a steeper size–fecundity relationship in cold environments. If large organisms have a greater advantage in terms of fecundity, it may be optimal to delay maturation and reach a larger size. Here, we make the inverse argument: that delaying maturity is optimal when mortality is low, and the steeper reproductive scaling is the product of delaying reproduction. It is important to note that, although the changes in mortality and growth with latitude are certainly affected by a temperature gradient, other factors, such as light and predation [30], are likely influencing these patterns too. Experimental evolution studies seem ideally suited to teasing apart these 2 hypotheses: for example, multigenerational, orthogonal manipulations of mortality rate and temperature, and measuring life history changes (e.g., maturation schedules, size, and fecundity scaling) should provide a test for whether mortality drives life history evolution or temperature effects drive fecundity and size evolution. Interestingly, these data may already exist for some species but have not been analysed in this context [31]. In terms of fecundity, the average tropical fish produces fewer offspring than the average polar fish, and this disparity increases with female size such that larger, colder fish produce far more offspring for a given female mass than smaller, warmer fish. These results are also consistent with life history patterns of freshwater species [32,33] as well as other marine studies [34]. This is not to say that higher-latitude fish necessarily have higher fitness, as contrasting life history strategies may yield similar population growth in different environments [35]. The congruence between the model predictions and the compiled reproductive data suggests that reproductive schedules might well be a product of life history optimisation given local demographic or environmental constraints, particularly those affecting mortality. It is worth noting that in our model, we consider mortality as the driver of resource allocation, yet it is equally likely that different life histories will result in different levels of mortality-determining causality is difficult. A useful illustration of this point comes from a recent synthesis whereby mortality rates are the product of life history strategies [23]. Certainly, there is evidence that mortality can shape life history evolution [36] and that life history strategies affect mortality (e.g., in semelparous species). We suspect that causality is bidirectional—there is an interplay between mortality and life history, but reconciling these effects into the same modelling framework represents a formidable challenge. Hyperallometric reproduction has important consequences for the replenishment of fish populations and how they are managed [37]. Fishing practices tend to disproportionally remove older, larger fish (i.e., “age-truncation”; [38,39]), shifting population size structure towards smaller sizes. Therefore, reductions in spawning biomass through harvesting can have disproportionate consequences on reproductive output [40,41]. Our model predicts that accelerated mortality schedules associated with harvesting should lead to the evolution of lower mass-specific reproductive outputs—with concomitant losses in replenishment. Similarly, reproductive scaling is likely to strongly affect the spillover benefits of marine protected areas (MPAs) [8]. Both our model and data compilations show that high-latitude fisheries are likely to show steep reproductive scaling relative to the tropics. The degree to which reproductive hyperallometry matters for both fisheries models and MPAs is likely to systematically vary with latitude—while hyperallometry is ubiquitous, it is particularly pronounced near the poles. We would therefore predict that MPAs in polar regions, in particular, are likely to yield the greatest benefit for exploited populations. If the latitudinal patterns observed here are mostly driven by temperature, global warming will reshape fish life histories and their demography. The general expectation is that fish will shrink in response to warming oceans [42–45] and that the productivity of most fisheries will decline with rising temperatures [46]. However, patterns in fish sizes associated to climate change are variable in nature [47]; tropical species have shrunk with increasing temperature, while temperate species show much more variable responses [48]. Assuming that temperature drives latitudinal changes in mortality, growth, and offspring size, our results suggest that regardless of how fish sizes change with increasing temperatures, mass-specific reproductive outputs are still likely to decline. We find that fish in warmer waters have lower mass-specific reproductive output than fish in cooler waters, and this discrepancy is particularly pronounced in larger fish. If future populations show the same temperature relationships as we observe here, a given standing spawner biomass of fish in the future will have much lower reproductive outputs than that same biomass today. For example, if CO2 emissions remain high (i.e., an increase of 2.58°C in sea surface temperature; [49]), our model predicts that a 25-kg fish at 60° latitude will produce about 300,000 fewer eggs by the end of the century, a difference of 5% of its mass-specific fecundity (Figs 4, S5 and S6). Such changes would have worrying consequences for replenishment and, thus, the productivity of global fisheries and may impact the food web in ways that are difficult to anticipate [50]. However, these predictions critically depend on all the inputs (offspring size, growth, and mortality) changing with temperature in the same ways as seen in a latitude gradient—an untested, yet critical assumption. Certainly, temperature is not the only factor that affects fecundity and changes across latitudes, and the effects of climate change on fish reproduction will depend on the joint effects of climate change on all drivers. Predictions of marine fish reproductive schedules using life history optimisation prove similar to patterns observed in the compiled data. A valuable extension of the approach taken here would be to resolve the other parameters that likely change with latitude. For example, processes affecting resource acquisition and transformation to biomass and how they change with temperature were not considered here. Particularly, we expect the scaling of production to be temperature dependent, hence affecting the optimal life history strategy at each latitude. As an essential next step, we need a better understanding of how production scales with size, and how this relationship is affected by temperature.

We use a life history model proposed by Day and Taylor [9] to investigate the reproductive scaling arising from optimising energy allocation between growth and reproduction across ontogeny. In the model, growth follows different trajectories before and after maturity, and the change in energy allocation is determined by a continuous function. Investment into reproduction is initially zero and increases through time after maturity as determined by 2 parameters: α (age at maturity) and h (the rate at which the proportion of surplus energy invested into reproduction increases). Weight w(t) is a function of time (t) and initially only depends on a growth parameter k (not equivalent to the growth coefficient K in the von Bertalanffy growth equation) and an offspring size parameter (w0; theoretical weight at time = 0 as a proxy). We chose this model because it depends only on a few parameters for which we can get relatively reliable data from many fish species, yet it is able to accurately predict fish life histories. While more complex models might predict reproductive patterns for species or locations more accurately, increasing model complexity would have limited the number of species and locations for which we could find reliable estimates of all model parameters and thus would have come at the cost of generality. Before maturity, the growth trajectory is described by the differential equation dw / dt = P(w), where the production function P(w) is assumed to be After maturity, if u(t) denotes the proportion of resources invested in growth f(t), a proxy for fecundity, is given by If we further assume that the function u(t) is given by we can solve for weight as a function of time t, knowing age at maturity α, initial size w0, and the rate at which resources diverted to reproduction change over time h: These equations assume that production scales at 3/4 with body size, but results are qualitatively consistent when a scaling equals to 2/3 (i.e., as assumed in Day and Taylor [9]; it is enough to change 4 with 3, and 1/4 with 1/3). In fact, Eq 4 can be adjusted to any scaling parameter of production. To find the optimal energy allocation strategy for fish across latitudes, we predicted how demographic parameters change across latitudes using data from the literature. From the data (see data compilation and data analysis below), we obtained the relationships for natural mortality (M), k, and w0 with latitude. Using these relationships, we looked at each latitude for the values of α that would maximise lifetime reproductive output (R0) under a median value of h: where p is the initial burst of mortality (here arbitrarily fixed at 0.8; this parameter does not change optimization results), and T is the longevity and was set to the age where only 1% of the individuals are expected to survive (). We generated mass and fecundity relationships with the optimal age at maturity and optimal surplus energy allocation for each latitude. The reproductive scaling was then determined as the change in fecundity (log-scale) over the change in mass (log-scale) from maturity to T. Upon first glance, it would appear that fecundity should scale with mass at ¾, like production. However, changes in energy allocation from growth to fecundity also affect the scaling by diminishing the change in growth, and increasing the change in fecundity, allowing for hyperallometry. To test the robustness of the model predictions, we implemented a second life history model proposed by Kozłowski [10] (see below). In this model, resource allocation also changes throughout ontogeny. The main difference with the Day and Taylor model is that Kozłowski’s model assumes a seasonal environment and that changes in resource allocation are abrupt instead of gradual within each season starting from maturation; such bang-bang switching is optimal as proved by Kozłowski and Teriokhin [51]. Thus, each season is optimally divided between growth and reproduction. Numerical studies showed that results from Kozłowski’s model can be well approximated by Day and Taylor equation (see S1 Supporting Information), and, therefore, simpler Day and Taylor’s model was applied in the main body of the paper.

To test whether the optimal reproductive schedules change across latitudes, we compiled data on location-specific mortality, growth, and initial size to calibrate the life history model (S7 Fig). For mortality, we compiled data on the rate of natural mortality (M) from literature, including datasets in [52–54]. We used only estimates of M estimated from age-based catch curves or tag and release data (i.e., direct estimates) where the reference specified that the population was unfished or that fishing was insignificant, giving us a total of 131 observations for 74 species from 59 references. However, some references had multiple estimates of M for the same species across a small range of latitudes. For those references, we grouped estimates of the same species and used the mean M and the mean latitude instead, giving a total of 82 observations across latitudes 0° to 72.5°. To estimate the growth parameter k and initial weight w0, we compiled length at age data before maturity for 44 species from 32 references (n = 49), with latitudes ranging from 7.5° to 63°. Length was converted into mass using length–weight relationships in FishBase that were either species specific or based on closely related species with similar body shape [55,56]. For each species at each location, we fitted Eq 4 (with t<0) to the data. Similarly, to restrict the range of the possible values of h when optimising R0 using Day and Taylor’s model, we used the subset of these references where species approached their asymptotic mass to estimate h. In these instances, all mass at age data was used, and both Eq 2 and Eq 3 were fitted to the data. The subset of references for which h was estimated included 28 species across latitudes ranging from 15.4° to 61°. When optimizing maturity schedules in Day and Taylor’s model, h was fixed at its median value (h = 0.26; S8 Fig), but α was allowed to vary. Results are robust to fixing h at its median observed value because variation in R0 is consistently higher across values of α than across values of h (S9 and S10 Figs). Data on age at 50% maturity (A50) were also compiled from published literature. Age at 50% maturity is the age at which half of the individuals are expected to be mature. Although often estimates of A50 could be obtained indirectly when both length at 50% maturity (L50) and the von Bertalanffy growth equation were known, we only used estimates of A50 that were directly reported in the literature. The R2 decreased by 42% if indirect estimates of A50 were included (R2: 0.33, CI: 0.21 to 0.43 if indirect estimates of A50 were excluded versus R2: 0.14, CI: 0.07 to 0.20 if indirect estimates of A50 were included). We compiled a total of 127 observations of A50 from 34 references and 59 species, with latitudes ranging from 9° to 73°. We used data on batch fecundity–female size relationships for marine fish from the published literature. Most of the data had been previously compiled by Barneche and colleagues [1], but we searched for additional studies, particularly at high latitudes and in underrepresented areas to have a broader geographical scope. Only marine species were considered for the analysis; species that migrate to or from freshwater were not included. Sharks and rays were also excluded since their fecundity was much lower than for the rest of the species because of their viviparity or large egg envelops, and tuna were excluded because their sizes were much larger than for the rest of the species. We used the number of eggs per reproductive pulse (i.e., batch fecundity) as the measurement for fecundity. While the literature also includes other measurements of fecundity, such as annual fecundity, such data are much less common and therefore were excluded from the analysis. Importantly, when the model included changes in seasonality across latitudes (i.e., longer periods favourable for reproduction at low versus high latitudes), age at maturity and reproductive scaling were still predicted to increase with latitude (S13 Fig). If a species had fecundity–female size data at one location for more than 1 year, only the year with the highest number of data points was included. To accurately estimate mass–fecundity relationships for each species, we excluded species for which less than 10 data points per location were available. We compiled a total of 13,739 observations from 157 references and 153 species, with latitudes ranging from 0.9° to 74.5°.

All of the analyses were performed in R version 3.5.2 [57]. We used linear mixed effects models to estimate M, k, w0, fecundity (i.e., number of eggs per batch), and A50 using the package “brms” version 2.7.0 [58] assuming a Gaussian distribution. We fitted a linear model to predict M (in log-scale) as a function of absolute latitude (fixed effect) using only direct measurements of M. To predict k and w0 (both in log-scale), we fitted linear models with absolute latitude as a fixed effect. To estimate k and w0, we corrected for phylogenetic relatedness between species by including a random intercept of phylogeny with a covariance matrix reflecting the phylogenetic tree using the R package “rotl” [59]. For M, we chose not to correct for phylogenetic relatedness because our goal here was to understand the spatial distribution of different mortality rates or, in other words, what is the likely mortality rate of the “average fish” for a given latitude. We use predictions from these 3 models (S1 Table) to calibrate the life history model across latitudes. Life history models calibrated using the phylogenetically controlled estimates of M yielded qualitatively similar predictions (i.e., an increase in age at maturity and reproductive scaling with increasing latitude). The model estimating fecundity (number of eggs or offspring per batch, in log-scale) included female weight (g, in log-scale), latitude (absolute value), and their interactions as fixed effects, assuming a Gaussian distribution. Here, the reproductive scaling is the slope of the regression between natural logarithm of fecundity and the natural logarithm of female weight. Since each species had multiple observations at each location, species was included as a random effect that could modify the intercept and the slope of relationship between fecundity and female weight. In addition, a random intercept was included for phylogeny, associated to a variance–covariance that corrected for phylogenetic relatedness. The model estimating A50 (years, in log-scale) included latitude as a fixed effect and the random intercept correcting for phylogenetic relatedness. The model was fitted assuming a Gaussian distribution.

To estimate the expected change in fecundity for a given female mass with climate change, we used predicted increases in water temperature under different climate change scenarios: a low emissions scenario (RCP2.6) and a high emissions scenario (RCP8.5) [49]. We estimated mean sea surface water temperature (SST) (1981 to 2010) across the globe from the COBE-SST2 data provided by the NOAA/OAR/ESRL PSL (Boulder, Colorado, USA), from their website at https://psl.noaa.gov/. Using these data, we fitted a general additive model (GAM) to predict SST as a function of absolute latitude (S14 Fig). At each latitude, we calculated the predicted SST mid-century (2031 to 2050) and late-century (2081 to 2100) for the low and high emissions scenarios predicted by the IPCC (an increase of 0.64°C and 0.73°C for the mid- and late-century for the low emissions scenario, and a 0.95°C and a 2.58°C for the high emissions scenario). The IPCC predictions were relative to temperature in 1986 to 2005, which is within the time interval when most of the fish reproductive data were collected (S15 Fig). For each scenario at each latitude, we calculated the predicted mid- and late-century SST, and we found the latitude that currently experienced that mean temperature. Then, assuming fish reproductive scaling responds to changes in temperature, we estimated future reproduction by using the fitted mass–fecundity relationship for the latitude that matched the predicted SST. At low latitudes, an increase in SST was unmatched (e.g., for any future temperature >27°C, there were no latitudes that were currently experiencing that mean temperature), and, therefore, changes in reproduction could not be predicted.

Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file. Click here for additional data file.