Tolerance traits related to climate change resilience are independent and polygenic

Relatedness of hypoxia and hyperthermia tolerances in the Nile tilapia (Oreochromis niloticus) and their relationships with cardiac and gill traits

In fish, thermal and hypoxia tolerances may be functionally related, as suggested by the oxygen- and capacity-limited thermal tolerance (OCLTT) concept, which explains performance failure at high temperatures due to limitations in oxygen delivery. In this study the interrelatedness of hyperthermia and hypoxia tolerances in the Nile tilapia (Oreochromis niloticus), and their links to cardiorespiratory traits were examined. Different groups of O. niloticus (n = 51) were subjected to hypoxia and hyperthermia challenges and the O2 tension for aquatic surface respiration (ASR pO2) and critical thermal maximum (CTmax) were assessed as measurement endpoints. Gill filament length, total filament number, ventricle mass, length and width were also measured. Tolerance to hypoxia, as evidenced by ASR pO2 thresholds of the individual fish, was highly variable and varied between 0.26 and 3.39 kPa. ASR events increased more profoundly as O2 tensions decreased below 2 kPa. The CTmax values recorded for the O. niloticus individuals ranged from 43.1 to 44.8 °C (Mean: 44.2 ± 0.4 °C). Remarkably, there was a highly significant correlation between ASR pO2 and CTmax in O. niloticus (r = -0.76, p < 0.0001) with ASR pO2 increasing linearly with decreasing CTmax. There were, however, no discernible relationships between the measured cardiorespiratory properties and hypoxia or hyperthermia tolerances. The strong relationship between hypoxia and hyperthermia tolerances in this study may be related to the ability of the cardiorespiratory system to provide oxygen to respiring tissues under thermal stress, and thus provides some support for the OCLTT concept in this species, at least at the level of the entire organism.

Hypoxia alters the upper thermal limits and blood physiology in zebrafish, Danio rerio

Hypoxic aquatic environments occur more frequently as a result of climate change, thereby exerting challenges on the physiological and metabolic functions of aquatic animals. In this study, a model fish, zebrafish (Danio rerio) was used to observe the climate-induced hypoxic effect on the upper thermal limit (critical thermal maximum; CTmax), hemoglobin, and blood glucose levels, and abnormalities of erythrocytes at cellular and nuclear level. The value of CTmax decreased significantly under hypoxia (39.10 ± 0.96 °C) compared to normoxia (43.70 ± 0.91 °C). At CTmax, hemoglobin levels were much lower (9.33 ± 0.60 g/dL) and blood glucose levels were significantly higher (194.20 ± 11.33 mg/L) under hypoxia than they were under normoxia and at the beginning of the experiment. Increased frequencies of abnormalities in the erythrocytes at both cellular (fusion, twin, elongated, spindle and tear drop shaped) and nuclear (micronucleus, karyopyknosis, binuclei, nuclear degeneration and notched nuclei) levels were also found under hypoxia compared to normoxia. These results suggest that hypoxic conditions significantly alter the temperature tolerance and subsequent physiology in zebrafish. Our findings will aid in the development of effective management techniques for aquatic environments with minimum oxygen availability.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod

Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.

Population-specific vulnerability to ocean change in a multistressor environment

Variation in environmental conditions across a species' range can alter their responses to environmental change through local adaptation and acclimation. Evolutionary responses, however, may be challenged in ecosystems with tightly coupled environmental conditions, where changes in the covariance of environmental factors may make it more difficult for species to adapt to global change. Here, we conduct a 3-month-long mesocosm experiment and find evidence for local adaptation/acclimation in populations of red sea urchins, Mesocentrotus franciscanus, to multiple environmental drivers. Moreover, populations differ in their response to projected concurrent changes in pH, temperature, and dissolved oxygen. Our results highlight the potential for local adaptation/acclimation to multivariate environmental regimes but suggest that thresholds in responses to a single environmental variable, such as temperature, may be more important than changes to environmental covariance. Therefore, identifying physiological thresholds in key environmental drivers may be particularly useful for preserving biodiversity and ecosystem functioning.

Comparing mixed models and Random Forest association tests using naturalGWAS and a Striped Bass SNP dataset

In this study, we used the phenotype simulation package naturalGWAS to test the performance of Zhao's Random Forest method in comparison to an uncorrected Random Forest test, latent factor mixed models (LFMM), genome-wide efficient mixed models (GEMMA), and confounder adjusted linear regression (CATE). We created 400 sets of phenotypes, corresponding to five effect sizes and 2, 5, 15, or 30 causal loci, simulated from two empirical datasets containing SNPs from Striped Bass representing three and 13 populations. All association methods were evaluated for their ability to detect genotype-phenotype associations based on power, false discovery rates, and number of false positives. Genomic inflation was highest for uncorrected Random Forest and LFMM tests and lowest for Gemma and Zhao's Random Forest. All association tests had similar power to detect causal loci, and Zhao's Random Forest had the lowest false discovery rate in all scenarios. To measure the performance of association tests in small datasets with few loci surrounding a causal gene we also ran analyses again after removing causal loci from each dataset. All association tests were only able to find true positives, defined as loci located within 30k bp of a causal locus, in 3%-18% of simulations. In contrast, at least one false positive was found in 17%-44% of simulations. Zhao's Random Forest again identified the fewest false positives of all association tests studied. The ability to test the power of association tests for individual empirical datasets can be an extremely useful first step when designing a GWAS study.

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

Evolution of thermal physiology alters the projected range of threespine stickleback under climate change

Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of its range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographical range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data were either static ("no evolution" models), allowed to evolve at observed evolutionary rates ("evolution" models) or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL; "scaled evolution" models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.

Use of complex physiological traits as ecotoxicological biomarkers in tropical freshwater fishes

We review the use of complex physiological traits, of tolerance and performance, as biomarkers of the toxicological effects of contaminants in subtropical and tropical freshwater fishes. Such traits are growing in relevance due to climate change, as exposure to contaminants may influence the capacity of fishes to tolerate and perform in an increasingly stressful environment. We review the evidence that the critical oxygen level, a measure of hypoxia tolerance, provides a valuable biomarker of impacts of diverse classes of contaminants. When coupled with measures of cardiorespiratory variables, it can provide insight into mechanisms of toxicity. The critical thermal maximum, a simple measure of tolerance of acute warming, also provides a valuable biomarker despite a lack of understanding of its mechanistic basis. Its relative ease of application renders it useful in the rapid evaluation of multiple species, and in understanding how the severity of contaminant impacts depends upon prevailing environmental temperature. The critical swimming speed is a measure of exercise performance that is widely used as a biomarker in temperate species but very few studies have been performed on subtropical or tropical fishes. Overall, the review serves to highlight a critical lack of knowledge for subtropical and tropical freshwater fishes. There is a real need to expand the knowledge base and to use physiological biomarkers in support of decision making to manage tropical freshwater fish populations and their habitats, which sustain rich biodiversity but are under relentless anthropogenic pressure.

Interindividual plasticity in metabolic and thermal tolerance traits from populations subjected to recent anthropogenic heating

To better understand temperature's role in the interaction between local evolutionary adaptation and physiological plasticity, we investigated acclimation effects on metabolic performance and thermal tolerance among natural Fundulus heteroclitus (small estuarine fish) populations from different thermal environments. Fundulus heteroclitus populations experience large daily and seasonal temperature variations, as well as local mean temperature differences across their large geographical cline. In this study, we use three populations: one locally heated (32°C) by thermal effluence (TE) from the Oyster Creek Nuclear Generating Station, NJ, and two nearby reference populations that do not experience local heating (28°C). After acclimation to 12 or 28°C, we quantified whole-animal metabolic (WAM) rate, critical thermal maximum (CTmax) and substrate-specific cardiac metabolic rate (CaM, substrates: glucose, fatty acids, lactate plus ketones plus ethanol, and endogenous (i.e. no added substrates)) in approximately 160 individuals from these three populations. Populations showed few significant differences due to large interindividual variation within populations. In general, for WAM and CTmax, the interindividual variation in acclimation response (log2 ratio 28/12°C) was a function of performance at 12°C and order of acclimation (12-28°C versus 28-12°C). CTmax and WAM were greater at 28°C than 12°C, although WAM had a small change (2.32-fold) compared with the expectation for a 16°C increase in temperature (expect 3- to 4.4-fold). By contrast, for CaM, the rates when acclimatized and assayed at 12 or 28°C were nearly identical. The small differences in CaM between 12 and 28°C temperature were partially explained by cardiac remodeling where individuals acclimatized to 12°C had larger hearts than individuals acclimatized to 28°C. Correlation among physiological traits was dependent on acclimation temperature. For example, WAM was negatively correlated with CTmax at 12°C but positively correlated at 28°C. Additionally, glucose substrate supported higher CaM than fatty acid, and fatty acid supported higher CaM than lactate, ketones and alcohol (LKA) or endogenous. However, these responses were highly variable with some individuals using much more FA than glucose. These findings suggest interindividual variation in physiological responses to temperature acclimation and indicate that additional research investigating interindividual may be relevant for global climate change responses in many species.

Intraspecific variation in tolerance of warming in fishes

Intraspecific variation in key traits such as tolerance of warming can have profound effects on ecological and evolutionary processes, notably responses to climate change. The empirical evidence for three primary elements of intraspecific variation in tolerance of warming in fishes is reviewed. The first is purely mechanistic that tolerance varies across life stages and as fishes become mature. The limited evidence indicates strongly that this is the case, possibly because of universal physiological principles. The second is intraspecific variation that is because of phenotypic plasticity, also a mechanistic phenomenon that buffers individuals' sensitivity to negative impacts of global warming in their lifetime, or to some extent through epigenetic effects over successive generations. Although the evidence for plasticity in tolerance to warming is extensive, more work is required to understand underlying mechanisms and to reveal whether there are general patterns. The third element is intraspecific variation based on heritable genetic differences in tolerance, which underlies local adaptation and may define long-term adaptability of a species in the face of ongoing global change. There is clear evidence of local adaptation and some evidence of heritability of tolerance to warming, but the knowledge base is limited with detailed information for only a few model or emblematic species. There is also strong evidence of structured variation in tolerance of warming within species, which may have ecological and evolutionary significance irrespective of whether it reflects plasticity or adaptation. Although the overwhelming consensus is that having broader intraspecific variation in tolerance should reduce species vulnerability to impacts of global warming, there are no sufficient data on fishes to provide insights into particular mechanisms by which this may occur.

The role of mechanistic physiology in investigating impacts of global warming on fishes

Warming of aquatic environments as a result of climate change is already having measurable impacts on fishes, manifested as changes in phenology, range shifts and reductions in body size. Understanding the physiological mechanisms underlying these seemingly universal patterns is crucial if we are to reliably predict the fate of fish populations with future warming. This includes an understanding of mechanisms for acute thermal tolerance, as extreme heatwaves may be a major driver of observed effects. The hypothesis of gill oxygen limitation (GOL) is claimed to explain asymptotic fish growth, and why some fish species are decreasing in size with warming; but its underlying assumptions conflict with established knowledge and direct mechanistic evidence is lacking. The hypothesis of oxygen- and capacity-limited thermal tolerance (OCLTT) has stimulated a wave of research into the role of oxygen supply capacity and thermal performance curves for aerobic scope, but results vary greatly between species, indicating that it is unlikely to be a universal mechanism. As thermal performance curves remain important for incorporating physiological tolerance into models, we discuss potentially fruitful alternatives to aerobic scope, notably specific dynamic action and growth rate. We consider the limitations of estimating acute thermal tolerance by a single rapid measure whose mechanism of action is not known. We emphasise the continued importance of experimental physiology, particularly in advancing our understanding of underlying mechanisms, but also the challenge of making this knowledge relevant to the more complex reality.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Mitochondria and the thermal limits of ectotherms

Temperature is a critical abiotic factor shaping the distribution and abundance of species, but the mechanisms that underpin organismal thermal limits remain poorly understood. One possible mechanism underlying these limits is the failure of mitochondrial processes, as mitochondria play a crucial role in animals as the primary site of ATP production. Conventional measures of mitochondrial performance suggest that these organelles can function at temperatures much higher than those that limit whole-organism function, suggesting that they are unlikely to set organismal thermal limits. However, this conclusion is challenged by recent data connecting sequence variation in mitochondrial genes to whole-organism thermal tolerance. Here, we review the current state of knowledge of mitochondrial responses to thermal extremes and ask whether they are consistent with a role for mitochondrial function in shaping whole-organism thermal limits. The available data are fragmentary, but it is possible to draw some conclusions. There is little evidence that failure of maximal mitochondrial oxidative capacity as assessed in vitro sets thermal limits, but there is some evidence to suggest that temperature effects on ATP synthetic capacity may be important. Several studies suggest that loss of mitochondrial coupling is associated with the thermal limits for organismal growth, although this needs to be rigorously tested. Most studies have utilized isolated mitochondrial preparations to assess the effects of temperature on these organelles, and there remain many untapped opportunities to address these questions using preparations that retain more of their biological context to better connect these subcellular processes with whole-organism thermal limits.

Mountain stoneflies may tolerate warming streams: Evidence from organismal physiology and gene expression

Rapid glacier recession is altering the physical conditions of headwater streams. Stream temperatures are predicted to rise and become increasingly variable, putting entire meltwater-associated biological communities at risk of extinction. Thus, there is a pressing need to understand how thermal stress affects mountain stream insects, particularly where glaciers are likely to vanish on contemporary timescales. In this study, we measured the critical thermal maximum (CT_MAX ) of stonefly nymphs representing multiple species and a range of thermal regimes in the high Rocky Mountains, USA. We then collected RNA-sequencing data to assess how organismal thermal stress translated to the cellular level. Our focal species included the meltwater stonefly, Lednia tumana, which was recently listed under the U.S. Endangered Species Act due to climate-induced habitat loss. For all study species, critical thermal maxima (CT_MAX  > 20°C) far exceeded the stream temperatures mountain stoneflies experience (<10°C). Moreover, while evidence for a cellular stress response was present, we also observed constitutive expression of genes encoding proteins known to underlie thermal stress (i.e., heat shock proteins) even at low temperatures that reflected natural conditions. We show that high-elevation aquatic insects may not be physiologically threatened by short-term exposure to warm temperatures and that longer-term physiological responses or biotic factors (e.g., competition) may better explain their extreme distributions.

Thermal tolerance and hypoxia tolerance are associated in blacktip reef shark (Carcharhinus melanopterus) neonates

Thermal dependence of growth and metabolism can influence thermal preference and tolerance in marine ectotherms, including threatened and data-deficient species. Here, we quantified the thermal dependence of physiological performance in neonates of a tropical shark species (blacktip reef shark, Carcharhinus melanopterus) from shallow, nearshore habitats. We measured minimum and maximum oxygen uptake rates (ṀO2), calculated aerobic scope, excess post-exercise oxygen consumption and recovery from exercise, and measured critical thermal maxima (CTmax), thermal safety margins, hypoxia tolerance, specific growth rates, body condition and food conversion efficiencies at two ecologically relevant acclimation temperatures (28 and 31°C). Owing to high post-exercise mortality, a third acclimation temperature (33°C) was not investigated further. Acclimation temperature did not affect ṀO2 or growth, but CTmax and hypoxia tolerance were greatest at 31°C and positively associated. We also quantified in vitro temperature (25, 30 and 35°C) and pH effects on haemoglobin–oxygen (Hb–O2) affinity of wild-caught, non-acclimated sharks. As expected, Hb–O2 affinity decreased with increasing temperatures, but pH effects observed at 30°C were absent at 25 and 35°C. Finally, we logged body temperatures of free-ranging sharks and determined that C. melanopterus neonates avoided 31°C in situ. We conclude that C. melanopterus neonates demonstrate minimal thermal dependence of whole-organism physiological performance across a seasonal temperature range and may use behaviour to avoid unfavourable environmental temperatures. The association between thermal tolerance and hypoxia tolerance suggests a common mechanism warranting further investigation. Future research should explore the consequences of ocean warming, especially in nearshore, tropical species.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Mitochondria, sex and variation in routine metabolic rate

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (ṀO_{2 routine} ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with ṀO_{2 routine} , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with ṀO_{2 routine} involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in ṀO_{2 routine} . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on ṀO_{2 routine} are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

Variation in Thermal Tolerance and Its Relationship to Mitochondrial Function Across Populations of Tigriopus californicus

Variation in thermal tolerance plays a key role in determining the biogeographic distribution of organisms. Consequently, identifying the mechanistic basis for thermal tolerance is necessary for understanding not only current species range limits but also the capacity for range limits to shift in response to climate change. Although variation in mitochondrial function likely contributes to variation in thermal tolerance, the extent to which mitochondrial function underlies local thermal adaptation is not fully understood. In the current study, we examine variation in thermal tolerance and mitochondrial function among three populations of the intertidal copepod Tigriopus californicus found across a latitudinal thermal gradient along the coast of California, USA. We tested (1) acute thermal tolerance using survivorship and knockdown assays, (2) chronic thermal tolerance using survivorship of nauplii and developmental rate, and (3) mitochondrial performance at a range of temperatures using ATP synthesis fueled by complexes I, II, and I&II, as well as respiration of permeabilized fibers. We find evidence for latitudinal thermal adaptation: the southernmost San Diego population outperforms the northernmost Santa Cruz in measures of survivorship, knockdown temperature, and ATP synthesis rates during acute thermal exposures. However, under a chronic thermal regime, survivorship and developmental rate are more similar in the southernmost and northernmost population than in the mid-range population (Abalone Cove). Though this pattern is unexpected, it aligns well with population-specific rates of ATP synthesis at these chronic temperatures. Combined with the tight correlation of ATP synthesis decline and knockdown temperature, these data suggest a role for mitochondria in setting thermal range limits and indicate that divergence in mitochondrial function is likely a component of adaptation across latitudinal thermal gradients.