The evolution of thermal performance curves in semi-aquatic newts: Thermal specialists on land and thermal generalists in water?

The position and shape of thermal performance curves (TPCs, the functions relating temperature to physiological performance) for ecologically relevant functions will directly affect the fitness of ectotherms and therefore should be under strong selection. However, thermodynamic considerations predict that relationships between the different components of the TPC will confound its evolutionary optimization. For instance, the “jack-of-all-temperatures” hypothesis predicts a trade-off between the breadth of the TPC and the maximal performance capacity; the “warmer is better” hypothesis suggests that low thermal optima will come with low absolute performances. Semi-aquatic organisms face the additional challenge of having to adjust their TPCs to two environments that are likely to differ in mean temperature and thermal variability. In this paper, we examine how parameters of the TPCs for maximal running and swimming speed have co-evolved in the semi-aquatic newt genus Triturus. We consider evolutionary relationships between the width and the height of the TPCs, the optimal temperatures and maximal performance. Phylogenetic comparative analyses reveal that in Triturus, swimming and running differ substantially in the (co-)variation of TPC parameters. Whereas evolutionary changes in the TPC for swimming primarily concern the shape of the curve (generalist versus specialist), most interspecific variation in running speed TPCs involves shifts in overall performance across temperatures.     

Thermal performance curves, phenotypic plasticity, and the time scales of temperature exposure

Thermal performance curves (TPCs) describe the effects of temperature on biological rate processes. Here, we use examples from our work on common killifish (Fundulus heteroclitus) to illustrate some important conceptual issues relating to TPCs in the context of using these curves to predict the responses of organisms to climate change. Phenotypic plasticity has the capacity to alter the shape and position of the TPCs for acute exposures, but these changes can be obscured when rate processes are measured only following chronic exposures. For example, the acute TPC for mitochondrial respiration in killifish is exponential in shape, but this shape changes with acclimation. If respiration rate is measured only at the acclimation temperature, the TPC is linear, concealing the underlying mechanistic complexity at an acute time scale. These issues are particularly problematic when attempting to use TPCs to predict the responses of organisms to temperature change in natural environments. Many TPCs are generated using laboratory exposures to constant temperatures, but temperature fluctuates in the natural environment, and the mechanisms influencing performance at acute and chronic time scales, and the responses of the performance traits at these time scales may be quite different. Unfortunately, our current understanding of the mechanisms underlying the responses of organisms to temperature change is incomplete, particularly with respect to integrating from processes occurring at the level of single proteins up to whole-organism functions across different time scales, which is a challenge for the development of strongly grounded mechanistic models of responses to global climate change.     

Plasticity of size and growth in fluctuating thermal environments: comparing reaction norms and performance curves

Ectothermic animals exhibit two distinct kinds of plasticityin response to temperature: Thermal performance curves (TPCs),in which an individual's performance (e.g., growth rate) variesin response to current temperature; and developmental reactionnorms (DRNs), in which the trait value (e.g., adult body sizeor development time) of a genotype varies in response to developmentaltemperatures experienced over some time period during development.Here we explore patterns of genetic variation and selectionon TPCs and DRNs for insects in fluctuating thermal environments.First, we describe two statistical methods for partitioningtotal genetic variation into variation for overall size or performanceand variation in plasticity, and apply these methods to availabledatasets on DRNs and TPCs for insect growth and size. Our resultsindicate that for the datasets we considered, genetic variationin plasticity represents a larger proportion of the total geneticvariation in TPCs compared to DRNs, for the available datasets.Simulations suggest that estimates of the genetic variationin plasticity are strongly affected by the number and rangeof temperatures considered, and by the degree of nonlinearityin the TPC or DRN. Second, we review a recent analysis of fieldselection studies which indicates that directional selectionfavoring increased overall size is common in many systems—thatbigger is frequently fitter. Third, we use a recent theoreticalmodel to examine how selection on thermal performance curvesrelates to environmental temperatures during selection. Themodel predicts that if selection acts primarily on adult sizeor development time, then selection on thermal performance curvesfor larval growth or development rates is directly related tothe frequency distribution of temperatures experienced duringlarval development. Using data on caterpillar temperatures inthe field, we show that the strength of directional selectionon growth rate is predicted to be greater at the modal (mostfrequent) temperatures, not at the mean temperature or at temperaturesat which growth rate is maximized. Our results illustrate someof the differences in genetic architecture and patterns of selectionbetween thermal performance curves and developmental reactionnorms.     

Trade‐offs between morphology and thermal niches mediate adaptation in response to competing selective pressures

The effects of climate change—such as increased temperature variability and novel predators—rarely happen in isolation, but it is unclear how organisms cope with multiple stressors simultaneously. To explore this, we grew replicate Paramecium caudatum populations in either constant or variable temperatures and exposed half to predation. We then fit thermal performance curves (TPCs) of intrinsic growth rate (r_max) for each replicate population (N = 12) across seven temperatures (10°C–38°C). TPCs of P. caudatum exposed to both temperature variability and predation responded only to one or the other (but not both), resulting in unpredictable outcomes. These changes in TPCs were accompanied by changes in cell morphology. Although cell volume was conserved across treatments, cells became narrower in response to temperature variability and rounder in response to predation. Our findings suggest that predation and temperature variability produce conflicting pressures on both thermal performance and cell morphology. Lastly, we found a strong correlation between changes in cell morphology and TPC parameters in response to predation, suggesting that responses to opposing selective pressures could be constrained by trade‐offs. Our results shed new light on how environmental and ecological pressures interact to elicit changes in characteristics at both the individual and population levels. We further suggest that morphological responses to interactive environmental forces may modulate population‐level responses, making prediction of long‐term responses to environmental change challenging.     

Inferring temperature adaptation from thermal performance curves of somatic growth rate: The importance of growth measurements and mortality

When comparing somatic growth thermal performance curves (TPCs), higher somatic growth across experimental temperatures is often observed for populations originating from colder environments. Such countergradient variation has been suggested to represent adaptation to seasonality, or shorter favourable seasons in colder climates. Alternatively, populations from cold climates may outgrow those from warmer climates at low temperature, and vice versa at high temperature, representing adaptation to temperature. Using modelling, we show that distinguishing between these two types of adaptation based on TPCs requires knowledge about (i) the relationship between somatic growth rate and population growth rate, which in turn depends on the scale of somatic growth (absolute or proportional), and (ii) the relationship between somatic growth rate and mortality rate in the wild. We illustrate this by quantifying somatic growth rate TPCs for three populations of Daphnia magna where population growth scales linearly with proportional somatic growth. For absolute somatic growth, the northern population outperformed the two more southern populations across temperatures, and more so at higher temperatures, consistent with adaptation to seasonality. In contrast, for the proportional somatic growth TPCs, and hence population growth rate, TPCs tended to converge towards the highest temperatures. Thus, if the northern population pays an ecological mortality cost of rapid growth in the wild, this may create crossing population growth TPCs consistent with adaptation to temperature. Future studies within this field should be more explicit in how they extrapolate from somatic growth in the lab to fitness in the wild.     

The combined effects of reactant kinetics and enzyme stability explain the temperature dependence of metabolic rates

A mechanistic understanding of the response of metabolic rate to temperature is essential for understanding thermal ecology and metabolic adaptation. Although the Arrhenius equation has been used to describe the effects of temperature on reaction rates and metabolic traits, it does not adequately describe two aspects of the thermal performance curve (TPC) for metabolic rate—that metabolic rate is a unimodal function of temperature often with maximal values in the biologically relevant temperature range and that activation energies are temperature dependent. We show that the temperature dependence of metabolic rate in ectotherms is well described by an enzyme‐assisted Arrhenius (EAAR) model that accounts for the temperature‐dependent contribution of enzymes to decreasing the activation energy required for reactions to occur. The model is mechanistically derived using the thermodynamic rules that govern protein stability. We contrast our model with other unimodal functions that also can be used to describe the temperature dependence of metabolic rate to show how the EAAR model provides an important advance over previous work. We fit the EAAR model to metabolic rate data for a variety of taxa to demonstrate the model's utility in describing metabolic rate TPCs while revealing significant differences in thermodynamic properties across species and acclimation temperatures. Our model advances our ability to understand the metabolic and ecological consequences of increases in the mean and variance of temperature associated with global climate change. In addition, the model suggests avenues by which organisms can acclimate and adapt to changing thermal environments. Furthermore, the parameters in the EAAR model generate links between organismal level performance and underlying molecular processes that can be tested for in future work.     

The complex drivers of thermal acclimation and breadth in ectotherms

Thermal acclimation capacity, the degree to which organisms can alter their optimal performance temperature and critical thermal limits with changing temperatures, reflects their ability to respond to temperature variability and thus might be important for coping with global climate change. Here, we combine simulation modelling with analysis of published data on thermal acclimation and breadth (range of temperatures over which organisms perform well) to develop a framework for predicting thermal plasticity across taxa, latitudes, body sizes, traits, habitats and methodological factors. Our synthesis includes > 2000 measures of acclimation capacities from > 500 species of ectotherms spanning fungi, invertebrates, and vertebrates from freshwater, marine and terrestrial habitats. We find that body size, latitude, and methodological factors often interact to shape acclimation responses and that acclimation rate scales negatively with body size, contributing to a general negative association between body size and thermal breadth across species. Additionally, we reveal that acclimation capacity increases with body size, increases with latitude (to mid‐latitudinal zones) and seasonality for smaller but not larger organisms, decreases with thermal safety margin (upper lethal temperature minus maximum environmental temperatures), and is regularly underestimated because of experimental artefacts. We then demonstrate that our framework can predict the contribution of acclimation plasticity to the IUCN threat status of amphibians globally, suggesting that phenotypic plasticity is already buffering some species from climate change.     

Predicting the growth of the amphibian chytrid fungus in varying temperature environments

Environmental temperature is a crucial abiotic factor that influences the success of ectothermic organisms, including hosts and pathogens in disease systems. One example is the amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), which has led to widespread amphibian population declines. Understanding its thermal ecology is essential to effectively predict outbreaks. Studies that examine the impact of temperature on hosts and pathogens often do so in controlled constant temperatures. Although varying temperature experiments are becoming increasingly common, it is unrealistic to test every temperature scenario. Thus, reliable methods that use constant temperature data to predict performance in varying temperatures are needed. In this study, we tested whether we could accurately predict Bd growth in three varying temperature regimes, using a Bayesian hierarchical model fit with constant temperature Bd growth data. We fit the Bayesian hierarchical model five times, each time changing the thermal performance curve (TPC) used to constrain the logistic growth rate to determine how TPCs influence the predictions. We then validated the model predictions using Bd growth data collected from the three tested varying temperature regimes. Although all TPCs overpredicted Bd growth in the varying temperature regimes, some functional forms performed better than others. Varying temperature impacts on disease systems are still not well understood and improving our understanding and methodologies to predict these effects could provide insights into disease systems and help conservation efforts.     

Transgenerational and within‐generation plasticity shape thermal performance curves

Thermal performance curves (TPCs) compute the effects of temperature on the performance of ectotherms and are frequently used to predict the effect of environmental conditions and currently, climate change, on organismal vulnerability and sensitivity. Using Drosophila melanogaster as an animal model, we examined how different thermal environments affected the shape of the performance curve and their parameters. We measured the climbing speed as a measure of locomotor performance in adult flies and tested the ontogenetic and transgenerational effects of thermal environment on TPC shape. Parents and offspring were reared at 28 ± 0ºC (28C), 28 ± 4ºC (28V), and 30 ± 0ºC (30C). We found that both, environmental thermal variability (28V) and high temperature (30C) experienced during early ontogeny shaped the fruit fly TPC sensitivity. Flies reared at variable thermal environments shifted the TPC to the right and increased heat tolerance. Flies held at high and constant temperature exhibited lower maximum performance than flies reared at the variable thermal environment. Furthermore, these effects were extended to the next generation. The parental thermal environment had a significative effect on TPC and its parameters. Indeed, flies reared at 28V whose parents were held at a high and constant temperature (30C) had a lower heat tolerance than F1 of flies reared at 28C or 28V. Also, offspring of flies reared at variable thermal environment (28V) reached the maximum performance at a higher temperature than offspring of flies reared at 28C or 30C. Consequently, since TPC parameters are not fixed, we suggest cautiousness when using TPCs to predict the impact of climate change on natural populations.     

Integrating patterns of thermal tolerance and phenotypic plasticity with population genetics to improve understanding of vulnerability to warming in a widespread copepod

Differences in population vulnerability to warming are defined by spatial patterns in thermal adaptation. These patterns may be driven by natural selection over spatial environmental gradients, but can also be shaped by gene flow, especially in marine taxa with high dispersal potential. Understanding and predicting organismal responses to warming requires disentangling the opposing effects of selection and gene flow. We begin by documenting genetic divergence of thermal tolerance and developmental phenotypic plasticity. Ten populations of the widespread copepod Acartia tonsa were collected from sites across a large thermal gradient, ranging from the Florida Keys to Northern New Brunswick, Canada (spanning over 20° latitude). Thermal performance curves (TPCs) from common garden experiments revealed local adaptation at the sampling range extremes, with thermal tolerance increasing at low latitudes and decreasing at high latitudes. The opposite pattern was observed in phenotypic plasticity, which was strongest at high latitudes. No relationship was observed between phenotypic plasticity and environmental variables. Instead, the results are consistent with the hypothesis of a trade‐off between thermal tolerance and the strength of phenotypic plasticity. Over a large portion of the sampled range, however, we observed a remarkable lack of differentiation of TPCs. To examine whether this lack of divergence is the result of selection for a generalist performance curve or constraint by gene flow, we analyzed cytochrome oxidase I mtDNA sequences, which revealed four distinct genetic clades, abundant genetic diversity, and widely distributed haplotypes. Strong divergence in thermal performance within genetic clades, however, suggests that the pace of thermal adaptation can be relatively rapid. The combined insight from the laboratory physiological experiments and genetic data indicate that gene flow constrains differentiation of TPCs. This balance between gene flow and selection has implications for patterns of vulnerability to warming. Taking both genetic differentiation and phenotypic plasticity into account, our results suggest that local adaptation does not increase vulnerability to warming, and that low‐latitude populations in general may be more vulnerable to predicted temperature change over the next century.     

Connecting thermal performance curve variation to the genotype: a multivariate QTL approach

Thermal performance curves (TPCs) are continuous reaction norms that describe the relationship between organismal performance and temperature and are useful for understanding trade‐offs involved in thermal adaptation. Although thermal trade‐offs such as those between generalists and specialists or between hot‐ and cold‐adapted phenotypes are known to be genetically variable and evolve during thermal adaptation, little is known of the genetic basis to TPCs – specifically, the loci involved and the directionality of their effects across different temperatures. To address this, we took a multivariate approach, mapping quantitative trait loci (QTL) for locomotor activity TPCs in the fly, Drosophila serrata, using a panel of 76 recombinant inbred lines. The distribution of additive genetic (co)variance in the mapping population was remarkably similar to the distribution of mutational (co)variance for these traits. We detected 11 TPC QTL in females and 4 in males. Multivariate QTL effects were closely aligned with the major axes genetic (co)variation between temperatures; most QTL effects corresponded to variation for either overall increases or decreases in activity with a smaller number indicating possible trade‐offs between activity at high and low temperatures. QTL representing changes in curve shape such as the ‘generalist–specialist’ trade‐off, thought key to thermal adaptation, were poorly represented in the data. We discuss these results in the light of genetic constraints on thermal adaptation.     

Predictions of response to temperature are contingent on model choice and data quality

The equations used to account for the temperature dependence of biological processes, including growth and metabolic rates, are the foundations of our predictions of how global biogeochemistry and biogeography change in response to global climate change. We review and test the use of 12 equations used to model the temperature dependence of biological processes across the full range of their temperature response, including supra‐ and suboptimal temperatures. We focus on fitting these equations to thermal response curves for phytoplankton growth but also tested the equations on a variety of traits across a wide diversity of organisms. We found that many of the surveyed equations have comparable abilities to fit data and equally high requirements for data quality (number of test temperatures and range of response captured) but lead to different estimates of cardinal temperatures and of the biological rates at these temperatures. When these rate estimates are used for biogeographic predictions, differences between the estimates of even the best‐fitting models can exceed the global biological change predicted for a decade of global warming. As a result, studies of the biological response to global changes in temperature must make careful consideration of model selection and of the quality of the data used for parametrizing these models.     

Variation,selection and evolution of function-valued traits

We describe an emerging framework for understanding variation, selection and evolution of phenotypic traits that are mathematical functions. We use one specific empirical example – thermal performance curves (TPCs) for growth rates of caterpillars – to demonstrate how models for function-valued traits are natural extensions of more familiar, multivariate models for correlated, quantitative traits. We emphasize three main points. First, because function-valued traits are continuous functions, there are important constraints on their patterns of variation that are not captured by multivariate models. Phenotypic and genetic variation in function-valued traits can be quantified in terms of variance-covariance functions and their associated eigenfunctions: we illustrate how these are estimated as well as their biological interpretations for TPCs. Second, selection on a function-valued trait is itself a function, defined in terms of selection gradient functions. For TPCs, the selection gradient describes how the relationship between an organism's performance and its fitness varies as a function of its temperature. We show how the form of the selection gradient function for TPCs relates to the frequency distribution of environmental states (caterpillar temperatures) during selection. Third, we can predict evolutionary responses of function-valued traits in terms of the genetic variance-covariance and the selection gradient functions. We illustrate how non-linear evolutionary responses of TPCs may occur even when the mean phenotype and the selection gradient are themselves linear functions of temperature. Finally, we discuss some of the methodological and empirical challenges for future studies of the evolution of function-valued traits.     

Do ectotherms partition thermal resources? We still do not know

Partitioning of the niche space is a mechanism used to explain the coexistence of similar species. Ectotherms have variable body temperatures and their body temperatures influence performance and, ultimately, fitness. Therefore, many ectotherms use behavioral thermoregulation to avoid reduced capacities associated with body temperatures far from the optimal temperature for performance. Several authors have proposed that thermal niche partitioning in response to interspecific competition is a mechanism that allows the coexistence of similar species of ectotherms. We reviewed studies on thermal resource partitioning to evaluate the evidence for this hypothesis. In almost all studies, there was insufficient evidence to conclude unequivocally that thermal resource partitioning allowed species coexistence. Future studies should include sites where species are sympatric and sites where they are allopatric to rule out alternative mechanisms that cause differences in thermal traits between coexisting species. There is evidence of conservatism in the evolution of most thermal traits across a wide range of taxa, but thermal performance curves and preferred temperatures do respond to strong selection under laboratory conditions. Thus, there is potential for selection to act on thermal traits in response to interspecific competition. Nevertheless, more stringent tests of the thermal resource partitioning hypothesis are required before we can assess whether it is widespread in communities of ectotherms in nature.     

Variation in Thermal Sensitivity and Thermal Tolerances in an Invasive Species across a Climatic Gradient: Lessons from the Land Snail Cornu aspersum

The ability of organisms to perform at different temperatures could be described by a continuous nonlinear reaction norm (i.e., thermal performance curve, TPC), in which the phenotypic trait value varies as a function of temperature. Almost any shift in the parameters of this performance curve could highlight the direct effect of temperature on organism fitness, providing a powerful framework for testing thermal adaptation hypotheses. Inter-and intraspecific differences in this performance curve are also reflected in thermal tolerances limits (e.g., critical and lethal limits), influencing the biogeographic patterns of species’ distribution. Within this context, here we investigated the intraspecific variation in thermal sensitivities and thermal tolerances in three populations of the invasive snail Cornu aspersum across a geographical gradient, characterized by different climatic conditions. Thus, we examined population differentiation in the TPCs, thermal-coma recovery times, expression of heat-shock proteins and standard metabolic rate (i.e., energetic costs of physiological differentiation). We tested two competing hypotheses regarding thermal adaptation (the “hotter is better” and the generalist-specialist trade-offs). Our results show that the differences in thermal sensitivity among populations of C. aspersum follow a latitudinal pattern, which is likely the result of a combination of thermodynamic constraints (“hotter is better”) and thermal adaptations to their local environments (generalist-specialist trade-offs). This finding is also consistent with some thermal tolerance indices such as the Heat-Shock Protein Response and the recovery time from chill-coma. However, mixed responses in the evaluated traits suggest that thermal adaptation in this species is not complete, as we were not able to detect any differences in neither energetic costs of physiological differentiation among populations, nor in the heat-coma recovery.     

Can we predict ectotherm responses to climate change using thermal performance curves and body temperatures?

Thermal performance curves (TPCs), which quantify how an ectotherm's body temperature (T_b) affects its performance or fitness, are often used in an attempt to predict organismal responses to climate change. Here, we examine the key – but often biologically unreasonable – assumptions underlying this approach; for example, that physiology and thermal regimes are invariant over ontogeny, space and time, and also that TPCs are independent of previously experienced T_b. We show how a critical consideration of these assumptions can lead to biologically useful hypotheses and experimental designs. For example, rather than assuming that TPCs are fixed during ontogeny, one can measure TPCs for each major life stage and incorporate these into stage‐specific ecological models to reveal the life stage most likely to be vulnerable to climate change. Our overall goal is to explicitly examine the assumptions underlying the integration of TPCs with T_b, to develop a framework within which empiricists can place their work within these limitations, and to facilitate the application of thermal physiology to understanding the biological implications of climate change.     

Comparative experimental evolution reveals species-specific idiosyncrasies in marine phytoplankton adaptation to warming

A number of experimental studies have demonstrated that phytoplankton can display rapid thermal adaptation in response to warmed environments. While these studies provide insight into the evolutionary responses of single species, they tend to employ different experimental techniques. Consequently, our ability to compare the potential for thermal adaptation across different, ecologically relevant, species remains limited. Here, we address this limitation by conducting simultaneous long-term warming experiments with the same experimental design on clonal isolates of three phylogenetically diverse species of marine phytoplankton; the cyanobacterium Synechococcus sp., the prasinophyte Ostreococcus tauri and the diatom Phaeodoactylum tricornutum. Over the same experimental time period, we observed differing levels of thermal adaptation in response to stressful supra-optimal temperatures. Synechococcus sp. displayed the greatest improvement in fitness (i.e., growth rate) and thermal tolerance (i.e., temperature limits of growth). Ostreococcus tauri was able to improve fitness and thermal tolerance, but to a lesser extent. Finally, Phaeodoactylum tricornutum showed no signs of adaptation. These findings could help us understand how the structure of phytoplankton communities may change in response to warming, and possible biogeochemical implications, as some species show relatively more rapid adaptive shifts in their thermal tolerance.     

Two-pore channels (TPCs): Novel voltage-gated ion channels with pleiotropic functions

Two-pore channels (TPC1-3) comprise a subfamily of the eukaryotic voltage-gated ion channels (VGICs) superfamily that are mainly expressed in acidic stores in plants and animals. TPCS are widespread across the animal kingdom, with primates, mice and rats lacking TPC3, and mainly act as Ca⁺ and Na⁺ channels, although it was also suggested that they could be permeable to other ions. Nowadays, TPCs have been related to the development of different diseases, including Parkinson´s disease, obesity or myocardial ischemia. Due to this, their study has raised the interest of the scientific community to try to understand their mechanism of action in order to be able to develop an efficient drug that could regulate TPCs activity. In this review, we will provide an updated view regarding TPCs structure, function and activation, as well as their role in different pathophysiological processes.     

Two-pore channels form homo- and heterodimers

Two-pore channels (TPCs) have been recently identified as NAADP-regulated Ca(2+) release channels, which are localized on the endolysosomal system. TPCs have a 12-transmembrane domain (TMD) structure and are evolutionary intermediates between the 24-TMD α-subunits of Na(+) or Ca(2+) channels and the transient receptor potential channel superfamily, which have six TMDs in a single subunit and form tetramers with 24 TMDs as active channels. Based on this relationship, it is predicted that TPCs dimerize to form functional channels, but the dimerization of human TPCs has so far not been studied. Using co-immunoprecipitation studies and a mass spectroscopic analysis of the immunocomplex, we show the presence of homo- and heteromeric complexes for human TPC1 and TPC2. Despite their largely distinct localization, we identified a discrete number of endosomes that coexpressed TPC1 and TPC2. Homo- and heteromerization were confirmed by a FRET study, showing that both proteins interacted in a rotational (N- to C-terminal/head-to-tail) symmetry. This is the first report describing the presence of homomultimeric TPC1 channels and the first study showing that TPCs are capable of forming heteromers.     

The importance of parentage in assessing temperature effects on fish early life history: a review of the experimental literature

Climate change and human development are altering aquatic thermal regimes, highlighting the need to understand how fish fitness may be impacted across a generational boundary. We reviewed experimental temperature studies investigating the links between parents and progeny, asking questions regarding the taxa studied, broodfish used, offspring traits examined, experimental durations and research motivations. We identified forty-one peer-reviewed articles examining the effects of pre-spawning adult temperature holding on offspring. These studies showed a strong focus on the order Salmoniformes (46% of studies) and aquaculturally driven research (66%). The use of wild broodfish was rare (12%) and the majority of experiments (83%) did not examine offspring consequences beyond hatch. We also identified 56 articles investigating how incubation temperature and parental influences affect embryonic and larval development. We demonstrate that these studies are not common in comparison to the majority of incubation thermal experiments that do not employ controlled parental breeding designs. However, 52 out of 56 studies we reviewed reported maternal, paternal or family identity influenced offspring responses to temperature. In characterizing these studies, Salmoniformes were the most studied order (52%), wild broodfish were more commonly used (55%), aquaculture motivations were less evident (23%), and few studies investigated offspring performance or traits beyond endogenous yolk stages. Overall, we suggest it is beneficial to experimentally examine temperature with consideration to parent-progeny relationships. To broaden our current understanding of intergenerational temperature effects, we recommend an increased focus on wild populations, offspring physiological and performance measures, later offspring development stages, and expanding research in non-salmonid species.