Concurrent warming and browning eliminate cold-water fish habitat in many temperate lakes.

Cold-water species in temperate lakes face two simultaneous climate-driven ecosystem changes: warming and browning of their waters. Browning refers to reduced transparency arising from increased dissolved organic carbon (DOC), which absorbs solar energy near the surface. It is unclear whether the net effect is mitigation or amplification of climate warming impacts on suitable oxythermal habitat (<20 °C, >5 mgO/L) for cold-loving species because browning expands the vertical distribution of both cool water and oxygen depletion. We analyzed long-term trends and high-frequency sensor data from browning lakes in New York's Adirondack region to assess the contemporary status of summertime habitat for lacustrine brook trout. Across two decades, surface temperatures increased twice as fast and bottom dissolved oxygen declined >180% faster than average trends for temperate lakes. We identify four lake categories based on oxythermal habitat metrics: constrained, squeezed, overheated, and buffered. In most of our study lakes, trout face either seasonal loss (7 of 15) or dramatic restriction (12 to 21% of the water column; 5 of 15) of suitable habitat. These sobering statistics reflect rapid upward expansion of oxygen depletion in lakes with moderate or high DOC relative to compression of heat penetration. Only in very clear lakes has browning potentially mitigated climate warming. Applying our findings to extensive survey data suggests that decades of browning have reduced oxythermal refugia in most Adirondack lakes. We conclude that joint warming and browning may preclude self-sustaining cold-water fisheries in many temperate lakes; hence, oxythermal categorization is essential to guide triage strategies and management interventions.

CHARR efficiently estimates contamination from DNA sequencing data.

DNA sample contamination is a major issue in clinical and research applications of whole-genome and -exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a metric to estimate DNA sample contamination from variant-level whole-genome and -exome sequence data called CHARR, contamination from homozygous alternate reference reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VariantDataset format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole-genome and exome sequencing datasets.

Sex-specific effects of inbreeding in juvenile brown trout.

Inbreeding depression, that is, the reduction of health and vigour in individuals with high inbreeding coefficients, is expected to increase with environmental, social, or physiological stress. It has therefore been predicted that sexual selection and the associated stress usually lead to higher inbreeding depression in males than in females. However, sex-specific differences in life history may reverse that pattern during certain developmental stages. In some salmonids, for example, female juveniles start developing their gonads earlier than males who instead grow faster. We tested whether the sexes are differently affected by inbreeding during that time. To study the effects of inbreeding coefficients that may be typical for natural populations of brown trout (Salmo trutta), and also to control for potentially confounding maternal or paternal effects, we sampled males and females from the wild, used their gametes in a block-wise full-factorial breeding design to produce 60 full-sib families, released the offspring as yolk-sac larvae into the wild, sampled them 6 months later, identified their genetic sex, and used microsatellites to assign them to their parents. We used whole-genome resequencing to calculate the kinship coefficients for each breeding pair and hence the expected average inbreeding coefficient per family. Juvenile growth could be predicted from these expected inbreeding coefficients and the genetic sex: Females reached lower body sizes with increasing inbreeding coefficient, while no such link could be found in males. This sex-specific inbreeding depression led to the overall pattern that females were on average smaller than males by the end of their first summer.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs.

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

New indicators for monitoring genetic diversity applied to alpine brown trout populations using whole genome sequence data.

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's Ï´ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (FROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

Climate change vulnerability of Arctic char across Scandinavia.

Climate change is anticipated to cause species to shift their ranges upward and poleward, yet space for tracking suitable habitat conditions may be limited for range-restricted species at the highest elevations and latitudes of the globe. Consequently, range-restricted species inhabiting Arctic freshwater ecosystems, where global warming is most pronounced, face the challenge of coping with changing abiotic and biotic conditions or risk extinction. Here, we use an extensive fish community and environmental dataset for 1762 lakes sampled across Scandinavia (mid-1990s) to evaluate the climate vulnerability of Arctic char (Salvelinus alpinus), the world's most cold-adapted and northernly distributed freshwater fish. Machine learning models show that abiotic and biotic factors strongly predict the occurrence of Arctic char across the region with an overall accuracy of 89 percent. Arctic char is less likely to occur in lakes with warm summer temperatures, high dissolved organic carbon levels (i.e., browning), and presence of northern pike (Esox lucius). Importantly, climate warming impacts are moderated by habitat (i.e., lake area) and amplified by the presence of competitors and/or predators (i.e., northern pike). Climate warming projections under the RCP8.5 emission scenario indicate that 81% of extant populations are at high risk of extirpation by 2080. Highly vulnerable populations occur across their range, particularly near the southern range limit and at lower elevations, with potential refugia found in some mountainous and coastal regions. Our findings highlight that range shifts may give way to range contractions for this cold-water specialist, indicating the need for pro-active conservation and mitigation efforts to avoid the loss of Arctic freshwater biodiversity.

Spatial asynchrony and cross-scale climate interactions in populations of a coldwater stream fish.

Climate change affects populations over broad geographic ranges due to spatially autocorrelated abiotic conditions known as the Moran effect. However, populations do not always respond to broad-scale environmental changes synchronously across a landscape. We combined multiple datasets for a retrospective analysis of time-series count data (5-28 annual samples per segment) at 144 stream segments dispersed over nearly 1,000 linear kilometers of range to characterize the population structure and scale of spatial synchrony across the southern native range of a coldwater stream fish (brook trout, Salvelinus fontinalis), which is sensitive to stream temperature and flow variations. Spatial synchrony differed by life stage and geographic region: it was stronger in the juvenile life stage than in the adult life stage and in the northern sub-region than in the southern sub-region. Spatial synchrony of trout populations extended to 100-200 km but was much weaker than that of climate variables such as temperature, precipitation, and stream flow. Early life stage abundance changed over time due to annual variation in summer temperature and winter and spring stream flow conditions. Climate effects on abundance differed between sub-regions and among local populations within sub-regions, indicating multiple cross-scale interactions where climate interacted with local habitat to generate only a modest pattern of population synchrony over space. Overall, our analysis showed higher degrees of response heterogeneity of local populations to climate variation and consequently population asynchrony than previously shown based on analysis of individual, geographically restricted datasets. This response heterogeneity indicates that certain local segments characterized by population asynchrony and resistance to climate variation could represent unique populations of this iconic native coldwater fish that warrant targeted conservation. Advancing the conservation of this species can include actions that identify such priority populations and incorporate them into landscape-level conservation planning. Our approach is applicable to other widespread aquatic species sensitive to climate change.

Parental kinship coefficient but not paternal coloration predicts early offspring growth in lake char.

The 'good genes' hypotheses of sexual selection predict that females prefer males with strong ornaments because they are in good health and vigor and can afford the costs of the ornaments. A key assumption of this concept is that male health and vigor are useful predictors of genetic quality and hence offspring performance. We tested this prediction in wild-caught lake char (Salvelinus umbla) whose breeding coloration is known to reveal aspects of male health. We first reanalyzed results from sperm competition trials in which embryos of known parenthood had been raised singly in either a stress- or non-stress environment. Paternal coloration did not correlate with any measures of offspring performance. However, offspring growth was reduced with higher kinship coefficients between the parents. To test the robustness of these first observations, we collected a new sample of wild males and females, used their gametes in a full-factorial in vitro breeding experiment, and singly raised about 3000 embryos in either a stress- or non-stress environment (stress induced by microbes). Again, paternal coloration did not predict offspring performance, while offspring growth was reduced with higher kinship between the parents. We conclude that, in lake char, the genetic benefits of mate choice would be strongest if females could recognize and avoid genetically related males, while male breeding colors may be more relevant in intra-sexual selection.

Male lake char release taurocholic acid as part of a mating pheromone.

The evolutionary origins of sexual preferences for chemical signals remain poorly understood, due, in part, to scant information on the molecules involved. In the current study, we identified a male pheromone in lake char (Salvelinus namaycush) to evaluate the hypothesis that it exploits a non-sexual preference for juvenile odour. In anadromous char species, the odour of stream-resident juveniles guides migratory adults into spawning streams. Lake char are also attracted to juvenile odour but have lost the anadromous phenotype and spawn on nearshore reefs, where juvenile odour does not persist long enough to act as a cue for spawning site selection by adults. Previous behavioural data raised the possibility that males release a pheromone that includes components of juvenile odour. Using metabolomics, we found that the most abundant molecule released by males was also released by juveniles but not females. Tandem mass spectrometry and nuclear magnetic resonance were used to identify the molecule as taurocholic acid (TCA), which was previously implicated as a component of juvenile odour. Additional chemical analyses revealed that males release TCA at high rates via their urine during the spawning season. Finally, picomolar concentrations of TCA attracted pre-spawning and spawning females but not males. Taken together, our results indicate that male lake char release TCA as a mating pheromone and support the hypothesis that the pheromone is a partial match of juvenile odour.

Environmental variation associated with overwintering elicits marked metabolic plasticity in a temperate salmonid, Salvelinus fontinalis.

Poleward winters commonly expose animals, including fish, to frigid temperatures and low food availability. Fishes that remain active over winter must therefore balance trade-offs between conserving energy and maintaining physiological performance in the cold, yet the extent and underlying mechanisms of these trade-offs are not well understood. We investigated the metabolic plasticity of brook char (Salvelinus fontinalis), a temperate salmonid, from the biochemical to whole-animal level in response to cold and food deprivation. Acute cooling (1°C day-1) from 14°C to 2°C had no effect on food consumption but reduced activity by 77%. We then assessed metabolic performance and demand over 90 days with exposure to warm (8°C) or cold winter (2°C) temperatures while fish were fed or starved. Resting metabolic rate (RMR) decreased substantially during initial cooling from 8°C to 2°C (Q10=4.2-4.5) but brook char exhibited remarkable thermal compensation during acclimation (Q10=1.4-1.6). Conversely, RMR was substantially lower (40-48%) in starved fish, conserving energy. Thus, the absolute magnitude of thermal plasticity may be masked or modified under food restriction. This reduction in RMR was associated with atrophy and decreases in in vivo protein synthesis rates, primarily in non-essential tissues. Remarkably, food deprivation had no effect on maximum oxygen uptake rates and thus aerobic capacity, supporting the notion that metabolic capacity can be decoupled from RMR in certain contexts. Overall, our study highlights the multi-faceted energetic flexibility of Salvelinus spp. that likely contributes to their success in harsh and variable environments and may be emblematic of winter-active fishes more broadly.

Smaller body size under warming is not due to gill-oxygen limitation in a cold-water salmonid.

Declining body size in fishes and other aquatic ectotherms associated with anthropogenic climate warming has significant implications for future fisheries yields, stock assessments and aquatic ecosystem stability. One proposed mechanism seeking to explain such body-size reductions, known as the gill oxygen limitation (GOL) hypothesis, has recently been used to model future impacts of climate warming on fisheries but has not been robustly empirically tested. We used brook trout (Salvelinus fontinalis), a fast-growing, cold-water salmonid species of broad economic, conservation and ecological value, to examine the GOL hypothesis in a long-term experiment quantifying effects of temperature on growth, resting metabolic rate (RMR), maximum metabolic rate (MMR) and gill surface area (GSA). Despite significantly reduced growth and body size at an elevated temperature, allometric slopes of GSA were not significantly different than 1.0 and were above those for RMR and MMR at both temperature treatments (15°C and 20°C), contrary to GOL expectations. We also found that the effect of temperature on RMR was time-dependent, contradicting the prediction that heightened temperatures increase metabolic rates and reinforcing the importance of longer-term exposures (e.g. >6 months) to fully understand the influence of acclimation on temperature-metabolic rate relationships. Our results indicate that although oxygen limitation may be important in some aspects of temperature-body size relationships and constraints on metabolic supply may contribute to reduced growth in some cases, it is unlikely that GOL is a universal mechanism explaining temperature-body size relationships in aquatic ectotherms. We suggest future research focus on alternative mechanisms underlying temperature-body size relationships, and that projections of climate change impacts on fisheries yields using models based on GOL assumptions be interpreted with caution.

Species Difference? Bovine, Trout, and Human Plasma Protein Binding of Per- and Polyfluoroalkyl Substances.

Per- and polyfluoroalkyl substances (PFAS) strongly bind to proteins and lipids in blood, which govern their accumulation and distribution in organisms. Understanding the plasma binding mechanism and species differences will facilitate the quantitative in vitro-to-in vivo extrapolation and improve risk assessment of PFAS. We studied the binding mechanism of 16 PFAS to bovine serum albumin (BSA), trout, and human plasma using solid-phase microextraction. Binding of anionic PFAS to BSA and human plasma was found to be highly concentration-dependent, while trout plasma binding was linear for the majority of the tested PFAS. At a molar ratio of PFAS to protein ν < 0.1 molPFAS/molprotein, the specific protein binding of anionic PFAS dominated their human plasma binding. This would be the scenario for physiological conditions (ν < 0.01), whereas in in vitro assays, PFAS are often dosed in excess (ν > 1) and nonspecific binding becomes dominant. BSA was shown to serve as a good surrogate for human plasma. As trout plasma contains more lipids, the nonspecific binding to lipids affected the affinities of PFAS for trout plasma. Mass balance models that are parameterized with the protein-water and lipid-water partitioning constants (chemical characteristics), as well as the protein and lipid contents of the plasma (species characteristics), were successfully used to predict the binding to human and trout plasma.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Distinct dynamics in mountain watersheds: Exploring mercury and microplastic pollution-Unraveling the influence of atmospheric deposition, human activities, and hydrology.

The intensification of human activities all around the globe has led to the spread of micropollutants in high-mountain freshwater environments. We therefore aimed to assess the geospatial distribution and determine the potential sources of (total-) mercury (THg) and microplastics (MPs) in mountain freshwater ecosystems. To do so, we analyzed THg and MP concentrations in brown trout, biofilm, and sediments from lotic and lentic ecosystems in the Pyrenees - all subjected to different types of human pressure. Additionally, we assessed the potential impacts of these pollutants on fish, and explored the bioindication capacity of brown trout (Salmo trutta fario) and biofilm regarding THg and MP pollution. For the first time, we measured concentrations of MPs trapped in the matrix of freshwater biofilm. Our results suggest that THg in the Pyrenees might be explained by both legacy (regional) and distant sources, in combination with environmental characteristics such as the presence of peatlands or streamwater physicochemistry, while MPs in fish are linked to recent local pollution sources such as single-use plastics. In contrast, MPs in biofilm matrix and sediments indicate a combination of distant (i.e., atmospheric deposition) and recent local pollution sources. Moreover, hydrodynamics and plastic density likely control MP distribution in rivers. Based on Fulton's condition factor, we also found that higher THg concentrations caused a negative impact on fish health (K < 1), while no impact of MPs could be seen. Therefore, we suggest that brown trout and biofilm can serve as bioindicators of atmospheric deposition of THg in high-altitude lakes and that biofilm is a reliable bioindicator to assess MP pollution in remote environments. Brown trout may also act as a bioindicator of MP pollution, but only efficiently in more polluted areas.

Consistent Negative Correlations between Parasite Infection and Host Body Condition Across Seasons Suggest Potential Harmful Impacts of Salmincola markewitschi on Wild White-Spotted Charr, Salvelinus leucomaenis.

Assessing the impacts of parasites on wild fish populations is a fundamental and challenging aspect of the study of host-parasite relationships. Salmincola, a genus of ectoparasitic copepods, mainly infects salmonid species. This genus, which is notorious in aquaculture, damages host fishes, but its impacts under natural conditions remain largely unknown or are often considered negligible. In this study, we investigated the potential impacts of mouth-attaching Salmincola markewitschi on white-spotted charr (Salvelinus leucomaenis) through intensive field surveys across four seasons using host body condition as an indicator of harmful effects. The prevalence and parasite abundance were highest in winter and gradually decreased in summer and autumn, which might be due to host breeding and/or wintering aggregations that help parasite transmissions. Despite seasonal differences in prevalence and parasite abundance, consistent negative correlations between parasite abundance and host body condition were observed across all seasons, indicating that the mouth-attaching copepods could reduce the body condition of the host fish. This provides field evidence suggesting that S. markewitschi has a potential negative impact on wild white-spotted charr.

Association of ectoparasite Lepeophtheirus salmonis counts on farmed Atlantic salmon and wild sea trout in Scotland.

Parasitic sea lice (Copepoda: Caligidae) colonising marine salmonid (Salmoniformes: Salmonidae) aquaculture production facilities have been implicated as a possible pressure on wild salmon and sea trout populations. This investigation uses monitoring data from the mainland west coast and Western Isles of Scotland to estimate the association of the abundance of adult female Lepeophtheirus salmonis (Krøyer) colonising farmed Atlantic salmon Salmo salar L. with the occurrence of juvenile and mobile L. salmonis on wild sea trout, anadromous S. trutta L. The associations were evaluated using generalised linear mixed models incorporating farmed adult female salmon louse abundances which are temporally lagged relative to dependent wild trout values. The pattern of lags, which is consistent with time for L. salmonis development between egg and infective stage, was evaluated using model deviances. A significant positive association is identified between adult female L. salmonis abundance on farms and juvenile L. salmonis on wild trout. This association is consistent with a causal relationship in which increases in the number of L. salmonis copepodids originating from lice colonising farmed Atlantic salmon cause an increase of L. salmonis abundance on wild sea trout.

Species invasion progressively disrupts the trophic structure of native food webs.

Species invasions can have substantial impacts on native species and ecosystems, with important consequences for biodiversity. How these disturbances drive changes in the trophic structure of native food webs through time is poorly understood. Here, we quantify trophic disruption in freshwater food webs to invasion by an apex fish predator, lake trout, using an extensive stable isotope dataset across a natural gradient of uninvaded and invaded lakes in the northern Rocky Mountains, USA. Lake trout invasion increased fish diet variability (trophic dispersion), displaced native fishes from their reference diets (trophic displacement), and reorganized macroinvertebrate communities, indicating strong food web disruption. Trophic dispersion was greatest 25 to 50 y after colonization and dissipated as food webs stabilized in later stages of invasion (>50 y). For the native apex predator, bull trout, trophic dispersion preceded trophic displacement, leading to their functional loss in late-invasion food webs. Our results demonstrate how invasive species progressively disrupt native food webs via trophic dispersion and displacement, ultimately yielding biological communities strongly divergent from those in uninvaded ecosystems.

Non-lethal detection of Eubothrium crassum (Cestoda) in farmed Atlantic salmon, Salmo salar, using anal swabs and real-time PCR.

Detection of intestinal parasites in fish typically requires autopsy, resulting in the sacrifice of the fish. Here, we describe a non-lethal method for detecting the tapeworm Eubothrium crassum in fish using anal swabs and real-time PCR detection. Two assays were developed to detect cytochrome oxidase I (COI) mitochondrial DNA and 18S ribosomal DNA sequences of E. crassum, respectively. The assays were tested on swab samples from confirmed pathogen free Atlantic salmon (Salmo salar L.) and on samples from farmed Atlantic salmon, where the presence and intensity of parasites had been established through autopsy. The COI assay was shown to be specific to E. crassum, while the 18S assay also amplified the closely related E. salvelini, a species infecting Arctic charr (Salvelinus alpinus L.) in freshwater. The COI assay detected E. crassum in all field samples regardless of parasite load while the 18S assay failed to detect the parasite in two samples. The results thus demonstrates that this non-lethal approach can effectively detect E. crassum and can be a valuable tool in assessing the prevalence of infection in farmed salmon, aiding in treatment decisions and evaluating treatment effectiveness.

Thermal tolerance of cultured and wild Icelandic arctic charr (Salvelinus alpinus) at self-selected flow rates.

Climate change is predicted to change not only the temperature of many freshwater systems but also flow dynamics. Understanding how fishes will fare in the future requires knowing how they will respond to both extended variations of temperature and flow. Arctic charr have had their thermal tolerance measured, but never with respect to flow. Additionally, this circumpolar species has multiple populations exhibiting dramatic phenotypic plasticity which may mean that regional differences in thermal tolerance are unaccounted for. In Iceland, Arctic charr populations have experienced highly variable flow and temperature conditions over the past 10,000 years. The Icelandic climate, topography and geothermal activity have created a mosaic of freshwater habitats inhabited by charr that vary substantially in both temperature and flow. Our purpose was to test whether populations from these varied environments had altered thermal tolerance and whether phenotypic plasticity of thermal tolerance in charr depends on flow. We raised cultured Icelandic charr from hatch under a 2 X 2 matrix of flow and temperature and compared them to wild charr captured from matching flow and temperature environments. Wild fish were more thermally tolerant than cultured fish at both acclimation temperatures and were more thermally plastic. Icelandic Arctic charr were more thermally tolerant than comparison charr populations across Europe and North America, but only when acclimated to 13 °C; fish acclimated to 5 °C compared equably with comparison charr populations. Icelandic Arctic charr were also more thermally plastic than all but one other salmonine species. Neither flow of rearing or the flow selected during a thermal tolerance (CTmax) test factored into thermal tolerance. Thermal tolerance was also independent of body size, condition factor, heart and gill size. In summary, wild Icelandic Arctic charr have greater thermal tolerance and plasticity than predicted from the literature and their latitude, but artificial selection for properties like growth rate or fecundity may be breeding that increased tolerance out of cultured fish. As the world moves toward a warmer climate and increased dependence on cultured fish, this is a noteworthy result and merits further study.