Current knowledge on the biological responses of freshwater fish under projected scenarios of climate change remains limited. Here, we examine differences in the protein configuration of two endemic Iberian freshwater fish species, Squalius carolitertii and the critically endangered S. torgalensis that inhabit in the Atlantic-type northern and in the Mediterranean-type southwestern regions, respectively. We performed protein structure modeling of fourteen genes linked to protein folding, energy metabolism, circadian rhythms and immune responses. Structural differences in proteins between the two species were found for HSC70, FKBP52, HIF1α and GPB1. For S. torgalensis, besides structural differences, we found higher thermostability for two proteins (HSP90 and GBP1), which can be advantageous in a warmer environment. Additionally, we investigated how these species might respond to projected scenarios of 3° climate change warming, acidification (ΔpH = -0.4), and their combined effects. Significant changes in gene expression were observed in response to all treatments, particularly under the combined warming and acidification. While S. carolitertii presented changes in gene expression for multiple proteins related to folding (hsp90aa1, hsc70, fkbp4 and stip1), only one such gene was altered in S. torgalensis (stip1). However, S. torgalensis showed a greater capacity for energy production under both the acidification and combined scenarios by increasing cs gene expression and maintaining ldha gene expression in muscle. Overall, these findings suggest that S. torgalensis is better prepared to cope with projected climate change. Worryingly, under the simulated scenarios, disturbances to circadian rhythm and immune system genes (cry1aa, per1a and gbp1) raise concerns for the persistence of both species, highlighting the need to consider multi-stressor effects when evaluating climate change impacts upon fish. This work also highlights that assessments of the potential of endangered freshwater species to cope with environmental change are crucial to help decision-makers adopt future conservation strategies.

Background The circadian clock is a biological timing system that improves the ability of organisms to deal with environmental fluctuations. At the molecular level it consists of a network of transcription-translation feedback loops, involving genes that activate (bmal and clock – positive loop) and repress expression (cryptochrome (cry) and period (per) – negative loop). This is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, depend on the environmental temperature and thus are good models to study its integration within the circadian system. Here, we studied the molecular evolution of circadian genes in four Squalius freshwater fish species, distributed across Western Iberian rivers affected by two climatic types with different environmental conditions (e.g., light and temperature). S. carolitertii and S. pyrenaicus inhabit the colder northern region under Atlantic climate type, while S. torgalensis, S. aradensis and some populations of S. pyrenaicus inhabit the warmer southern region affected by summer droughts, under Mediterranean climate type. Results We identified 16 circadian-core genes in the Squalius species using a comparative transcriptomics approach. We detected evidence of positive selection in 12 of these genes using methods based on dN/dS. Positive selection was mainly found in cry and per genes of the negative loop, with 55 putatively adaptive substitutions, 16 located on protein domains. Evidence for positive selection is predominant in southern populations affected by the Mediterranean climate type. By predicting protein features we found that changes at sites under positive selection can impact protein thermostability by changing their aliphatic index and isoelectric point. Additionally, in nine genes, the phylogenetic clustering of species that belong to different clades but inhabit southern basins with similar environmental conditions indicated evolutionary convergence. We found evidence for increased nonsynonymous substitution rate in convergent lineages, likely due to positive selection at 27 sites, mostly in cry genes. Conclusions Our results support that temperature may be a selective pressure driving the evolution of genes involved in the circadian system. By integrating sequence-based functional protein prediction with dN/dS-based methods to detect selection we uncovered adaptive convergence in the southern populations, probably related to their similar thermal conditions.

The circadian clock is a biological timing system that improves the inherent ability of organisms to deal with environmental fluctuations. It is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, have an increased dependence on the environmental temperature and thus are good models to study the integration of temperature within the circadian system. Here, we studied four species of freshwater fish of Squalius genus, distributed across a latitudinal gradient in Portugal with variable conditions of light and temperature. We identified and characterised the expected sixteen genes belonging to four main gene families (Cryptochromes, Period, CLOCK and BMAL) involved in the circadian system. We conducted phylogenetic analysis and tested these genes for the presence of signatures of selection based on dN/dS that revealed the presence of mutations under positive selection in nine genes. Furthermore, we analysed functional features of the proteins and found that these mutations alter protein isoelectric point and aliphatic index. Additionally, in 50% of the genes we found a pattern of convergence uncovered by the phylogenetic clustering of S. aradensis and S. pyrenaicus from Almargem, which inhabit basins with similar environmental conditions. Although this clustering can be explained by incomplete lineage sorting, evidence of positive selection in one gene in the two species and similar functional features in seven proteins are consistent with convergence at protein level. Our results highlight the importance of integrating sequence-based functional protein characterisation with dN/dS-based methods when studying the molecular adaptation of species.

BackgroundThe circadian clock is a biological timing system that improves the ability of organisms to deal with environmental fluctuations. At the molecular level it consists of a network of transcription-translation feedback loops, involving genes that activate (bmal and clock – positive loop) and repress expression (cryptochrome (cry) and period (per) – negative loop). This is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, depend on the environmental temperature and thus are good models to study its integration within the circadian system. Here, we studied the molecular evolution of circadian genes in four Squalius freshwater fish species, distributed across Western Iberian rivers affected by two climatic types with different environmental conditions (e.g., light and temperature). S. carolitertii and S. pyrenaicus inhabit the colder northern region under Atlantic climate type, while S. torgalensis, S. aradensis and some populations of S. pyrenaicus inhabit the warmer southern region affected by summer droughts, under Mediterranean climate type. ResultsWe identified 16 circadian-core genes in the Squalius species using a comparative transcriptomics approach. We detected evidence of positive selection in 12 of these genes using methods based on dN/dS. Positive selection was mainly found in cry and per genes of the negative loop of the cycle, with 55 putatively adaptive substitutions, 16 located on protein domains. Evidence for positive selection is predominant in southern populations affected by the Mediterranean climate type. By predicting protein features we found that changes at sites under positive selection can impact protein thermostability by changing their aliphatic index and isoelectric point. Additionally, in nine genes, the phylogenetic clustering of species that belong to different clades but inhabit southern basins with similar environmental conditions indicated evolutionary convergence. We found evidence for increased nonsynonymous substitution rate in convergent lineages, likely due to positive selection at 27 sites, mostly in cry genes. ConclusionsOur results support that temperature may be a selective pressure driving the evolution of genes involved in the circadian system. By integrating sequence-based functional protein prediction with dN/dS-based methods to detect selection we uncovered adaptive convergence in the southern populations, probably related to their similar thermal conditions.

BackgroundThe circadian clock is a biological timing system that improves the ability of organisms to deal with environmental fluctuations. At the molecular level it consists of a network of transcription-translation feedback loops, involving genes that activate (bmal and clock – positive loop) and repress expression (cryptochrome (cry) and period (per) – negative loop). This is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, depend on the environmental temperature and thus are good models to study its integration within the circadian system. Here, we studied the molecular evolution of circadian genes in four Squalius freshwater fish species, distributed across Western Iberian rivers affected by two climatic types with different environmental conditions (e.g. light and temperature). S. carolitertii and S. pyrenaicus inhabit the colder northern region under Atlantic climate type, while S. torgalensis, S. aradensis and some populations of S. pyrenaicus inhabit the warmer southern region affected by summer droughts, under Mediterranean climate type. ResultsWe identified 16 circadian-core genes in the Squalius species using a comparative transcriptomics approach. We detected evidence of positive selection in nine of these genes using methods based on dN/dS. Positive selection was mainly found in cry and per genes of the negative loop of the cycle, with 11 putatively adaptive substitutions mostly located on protein domains. Evidence for positive selection is predominant in southern populations affected by the Mediterranean climate type. By predicting protein features we found that changes at sites under positive selection can impact protein thermostability by changing their aliphatic index and isoelectric point. Additionally, in nine genes, the phylogenetic clustering of species that belong to different clades but inhabit southern basins with similar environmental conditions indicated evolutionary convergence.ConclusionsOur results support that temperature may be a strong selective pressure driving the evolution of genes involved in the circadian system. By integrating sequence-based functional protein prediction with dN/dS-based methods to detect selection we also uncovered adaptive convergence in the southern populations, probably related to their similar thermal conditions.