Differential Expression of Genes that Control Respiration Contribute to Thermal Adaptation in Redband Trout (Oncorhynchus mykiss gairdneri)

The Constrained Disorder Principle Accounts for the Variability That Characterizes Breathing: A Method for Treating Chronic Respiratory Diseases and Improving Mechanical Ventilation

Variability characterizes breathing, cellular respiration, and the underlying quantum effects. Variability serves as a mechanism for coping with changing environments; however, this hypothesis does not explain why many of the variable phenomena of respiration manifest randomness. According to the constrained disorder principle (CDP), living organisms are defined by their inherent disorder bounded by variable boundaries. The present paper describes the mechanisms of breathing and cellular respiration, focusing on their inherent variability. It defines how the CDP accounts for the variability and randomness in breathing and respiration. It also provides a scheme for the potential role of respiration variability in the energy balance in biological systems. The paper describes the option of using CDP-based artificial intelligence platforms to augment the respiratory process's efficiency, correct malfunctions, and treat disorders associated with the respiratory system.

Intraspecific variation among Chinook Salmon populations indicates physiological adaptation to local environmental conditions

Understanding interpopulation variation is important to predicting species responses to climate change. Recent research has revealed interpopulation variation among several species of Pacific salmonids; however, the environmental drivers of population differences remain elusive. We tested for local adaptation and countergradient variation by assessing interpopulation variation among six populations of fall-run Chinook Salmon from the western United States. Juvenile fish were reared at three temperatures (11, 16 and 20°C), and five physiological metrics were measured (routine and maximum metabolic rate, aerobic scope, growth rate and critical thermal maximum). We then tested associations between these physiological metrics and 15 environmental characteristics (e.g. rearing temperature, latitude, migration distance, etc.). Statistical associations between the five physiological metrics and 15 environmental characteristics supported our hypotheses of local adaptation. Notably, latitude was a poor predictor of population physiology. Instead, our results demonstrate that populations from warmer habitats exhibit higher thermal tolerance (i.e. critical thermal maxima), faster growth when warm acclimated and greater aerobic capacity at high temperatures. Additionally, populations with longer migrations exhibit higher metabolic capacity. However, overall metabolic capacity declined with warm acclimation, indicating that future climate change may reduce metabolic capacity, negatively affecting long-migrating populations. Linking physiological traits to environmental characteristics enables flexible, population-specific management of disparate populations in response to local conditions.

Dreissena polymorpha responses under thermal and hypoxic stress: New insights in the tolerance of this freshwater sentinel species

Dreissena polymorpha is a sentinel freshwater mussel providing key functional ecosystemic services like nutrient recycling and suspended matter filtration. Global warming and especially extreme events imply rapid fluctuations of environmental parameters that sessile organisms could not escape. The increase occurrence of heat waves and the subsequent expansion of hypoxic areas could challenge the survival of mussels. This study provided a deeper knowledge of energy management and cellular function during thermal (+15 °C) or hypoxic (30% of dissolved oxygen saturation) stress for 7 days. A potential metabolic rate depression was highlighted in D. polymorpha under hypoxia through a decline in the mitochondrial activity and a constant AMP content over time. A contrasted pattern of response was observed in thermal-stressed mussels between 24 h and 7 days of exposure. A global increase of metabolic activity was noticed in mussels after 24 h while a return to control level was noticed at the end of the experiment. Although D. polymorpha is considered as a temperature tolerant species, a significant increase of ADP:ATP ratio, related to a decrease of mitochondrial activity and density, suggested an overwhelming of organisms. This study pointed to the importance of considering time of exposure to natural factor variations in tolerance window of organisms in a long-term changing environment. The apparent short-term tolerance of D. polymorpha could hide much more deleterious consequences, i.e. mortality, if abiotic stresses persist, as suggested by climate change models.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

Whole genome resequencing identifies local adaptation associated with environmental variation for redband trout

Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low-coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L, ROCK1, SIX6). Genotype-environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high-elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high-elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Mitogenomics of the endemic Ethiopian rats: looking for footprints of adaptive evolution in sky islands

Organisms living in high altitude must adapt to environmental conditions with hypoxia and low temperature, e.g. by changes in the structure and function of proteins associated with oxidative phosphorylation in mitochondria. Here we analysed the signs of adaptive evolution in 27 mitogenomes of endemic Ethiopian rats (Stenocephalemys), where individual species adapted to different elevation. Significant signals of positive selection were detected in 10 of the 13 mitochondrial protein-coding genes, with a majority of functional substitutions in the NADH dehydrogenase complex. Higher frequency of positively selected sites was found in phylogenetic lineages corresponding to Afroalpine specialists.

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole-genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency-based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%-75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.

Adjustments of cardiac mitochondrial phenotype in a warmer thermal habitat is associated with oxidative stress in European perch, Perca fluviatilis

Mitochondria are playing key roles in setting the thermal limits of fish, but how these organelles participate in selection mechanisms during extreme thermal events associated with climate warming in natural populations is unclear. Here, we investigated the thermal effects on mitochondrial metabolism, oxidative stress, and mitochondrial gene expression in cardiac tissues of European perch (Perca fluviatilis) collected from an artificially heated ecosystem, the "Biotest enclosure", and an adjacent reference area in the Baltic sea with normal temperatures (~ 23 °C and ~ 16 °C, respectively, at the time of capture in summer). Fish were sampled one month after a heat wave that caused the Biotest temperatures to peak at ~ 31.5 °C, causing significant mortality. When assayed at 23 °C, Biotest perch maintained high mitochondrial capacities, while reference perch displayed depressed mitochondrial functions relative to measurements at 16 °C. Moreover, mitochondrial gene expression of nd4 (mitochondrial subunit of complex I) was higher in Biotest fish, likely explaining the increased respiration rates observed in this population. Nonetheless, cardiac tissue from Biotest perch displayed higher levels of oxidative damage, which may have resulted from their chronically warm habitat, as well as the extreme temperatures encountered during the preceding summer heat wave. We conclude that eurythermal fish such as perch are able to adjust and maintain mitochondrial capacities of highly aerobic organs such as the heart when exposed to a warming environment as predicted with climate change. However, this might come at the expense of exacerbated oxidative stress, potentially threatening performance in nature.

Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.

Unravelling the changes during induced vitellogenesis in female European eel through RNA-Seq: What happens to the liver?

The life cycle of European eel (Anguilla anguilla), a catadromous species, is complex and enigmatic. In nature, during the silvering process prior to their long spawning migration, reproductive development is arrested, and they cease feeding. In studies of reproduction using hormonal induction, eels are equivalently not feed. Therefore, in female eels that undergo vitellogenesis, the liver plays different, essential roles being involved both in vitellogenins synthesis and in reallocating resources for the maintenance of vital functions, performing the transoceanic reproductive migration and completing reproductive development. The present work aimed at unravelling the major transcriptomic changes that occur in the liver during induced vitellogenesis in female eels. mRNA-Seq data from 16 animals (eight before induced vitellogenesis and eight after nine weeks of hormonal treatment) were generated and differential expression analysis was performed comparing the two groups. This analysis detected 1,328 upregulated and 1,490 downregulated transcripts. Overrepresentation analysis of the upregulated genes included biological processes related to biosynthesis, response to estrogens, mitochondrial activity and localization, while downregulated genes were enriched in processes related to morphogenesis and development of several organs and tissues, including liver and immune system. Among key genes, the upregulated ones included vitellogenin genes (VTG1 and VTG2) that are central in vitellogenesis, together with ESR1 and two novel genes not previously investigated in European eel (LMAN1 and NUPR1), which have been linked with reproduction in other species. Moreover, several upregulated genes, such as CYC1, ELOVL5, KARS and ACSS1, are involved in the management of the effect of fasting and NOTCH, VEGFA and NCOR are linked with development, autophagy and liver maintenance in other species. These results increase the understanding of the molecular changes that occur in the liver during vitellogenesis in this complex and distinctive fish species, bringing new insights on European eel reproduction and broodstock management.

Adaptation of plasticity to projected maximum temperatures and across climatically defined bioregions

Resilience to environmental stressors due to climate warming is influenced by local adaptations, including plastic responses. The recent literature has focused on genomic signatures of climatic adaptation, but little is known about how plastic capacity may be influenced by biogeographic and evolutionary processes. We investigate phenotypic plasticity as a target of climatic selection, hypothesizing that lineages that evolved in warmer climates will exhibit greater plastic adaptive resilience to upper thermal stress. This was experimentally tested by comparing transcriptomic responses within and among temperate, subtropical, and desert ecotypes of Australian rainbowfish subjected to contemporary and projected summer temperatures. Critical thermal maxima were estimated, and ecological niches delineated using bioclimatic modeling. A comparative phylogenetic expression variance and evolution model was used to assess plastic and evolved changes in gene expression. Although 82% of all expressed genes were found in the three ecotypes, they shared expression patterns in only 5 out of 236 genes that responded to the climate change experiment. A total of 532 genes showed signals of adaptive (i.e., genetic-based) plasticity due to ecotype-specific directional selection, and 23 of those responded to projected summer temperatures. Network analyses demonstrated centrality of these genes in thermal response pathways. The greatest adaptive resilience to upper thermal stress was shown by the subtropical ecotype, followed by the desert and temperate ecotypes. Our findings indicate that vulnerability to climate change will be highly influenced by biogeographic factors, emphasizing the value of integrative assessments of climatic adaptive traits for accurate estimation of population and ecosystem responses.

Transcriptomic responses to pollution in natural riverine environment in Rita rita

Aquatic pollution is one of the most common threats to the ecological health of aquatic ecosystems and its biota. Fish as lower vertebrates are excellent model to study the impact and responses of aquatic pollution. In fish, gill is the main organ indicator of whole animal health as it comes in contact with the surrounding water and absorbs many pollutants and contaminants; therefore, investigations on alterations in fish gill at transcriptome level could provide newer insights to the stress response mechanism(s) and pathways. For comprehensive evaluation of the impacts of pollutants (joint toxicity) prevalent in the riverine environment, comparative transcriptome analysis, by Next Generation Sequencing under Illumina HiSeq 2500 platform, was carried out in gill tissues of Rita rita collected from two stretches of river Ganga (Kanpur and Farakka) and results were validated by RT-qPCR. Out of 154,077 unigenes (Accession SRR548008), a total of 2024 differentially expressed genes (DEGs) including 942 up-regulated and 1082 down-regulated genes were identified by DESeq program. Further, Gene Ontology (GO) of DEGs showed that ribosomal large subunit biogenesis, mitochondrial ribosome and box H/ACA SnoRNA binding categories are highly affected by pollution. Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis showed the involvement of the DEGs in energy metabolism, translational and transcriptional machinery, protein folding and degradation suggesting that these signalling pathways are highly affected by aquatic pollution. Among the DEGs, up-regulation of cytochrome c oxidase subunit (cox) 7a2 (69.47 fold), hsp70 subunit 14 (hsp70-14, 5.27 fold), muscle related coiled-coil protein (MURC, 21.55 fold), lysozyme G (40.14 fold), cox17 (29.36 fold) were the conspicuous ones which showed similar trends in expression when analysed by RT-qPCR. Based on fold change, perturbation values, correlation analysis by PCA and RT-qPCR validation, up-regulation of cox7a2, MURC and hsp70-14 appeared to be the most promising biomarker responses and could be useful in the evaluation of gill health and possibly be extended towards aquatic ecosystem health assessment.

Signals of selection in the mitogenome provide insights into adaptation mechanisms in heterogeneous habitats in a widely distributed pelagic fish

Oceans are vast, dynamic, and complex ecosystems characterized by fluctuations in environmental parameters like sea surface temperature (SST), salinity, oxygen availability, and productivity. Environmental variability acts as the driver of organismal evolution and speciation as organisms strive to cope with the challenges. We investigated the evolutionary consequences of heterogeneous environmental conditions on the mitogenome of a widely distributed small pelagic fish of Indian ocean, Indian oil sardine, Sardinella longiceps. Sardines were collected from different eco-regions of the Indian Ocean and selection patterns analyzed in coding and non-coding regions. Signals of diversifying selection were observed in key functional regions involved in OXPHOS indicating OXPHOS gene regulation as the critical factor to meet enhanced energetic demands. A characteristic control region with 38–40 bp tandem repeat units under strong selective pressure as evidenced by sequence conservation and low free energy values was also observed. These changes were prevalent in fishes from the South Eastern Arabian Sea (SEAS) followed by the Northern Arabian Sea (NAS) and rare in Bay of Bengal (BoB) populations. Fishes belonging to SEAS exhibited accelerated substitution rate mainly due to the selective pressures to survive in a highly variable oceanic environment characterized by seasonal hypoxia, variable SST, and food availability.

Divergence history and hydrothermal vent adaptation of decapod crustaceans: A mitogenomic perspective

Decapod crustaceans, such as alvinocaridid shrimps, bythograeid crabs and galatheid squat lobsters are important fauna in the hydrothermal vents and have well adapted to hydrothermal vent environments. In this study, eighteen mitochondrial genomes (mitogenomes) of hydrothermal vent decapods were used to explore the evolutionary history and their adaptation to the hydrothermal vent habitats. BI and ML algorithms produced consistent phylogeny for Decapoda. The phylogenetic relationship revealed more evolved positions for all the hydrothermal vent groups, indicating they migrated from non-vent environments, instead of the remnants of ancient hydrothermal vent species, which support the extinction/repopulation hypothesis. The divergence time estimation on the Alvinocarididae, Bythograeidae and Galatheoidea nodes are located at 75.20, 56.44 and 47.41–50.43 Ma, respectively, which refers to the Late Cretaceous origin of alvinocaridid shrimps and the Early Tertiary origin of bythograeid crabs and galatheid squat lobsters. These origin stories are thought to associate with the global deep-water anoxic/dysoxic events. Total eleven positively selected sites were detected in the mitochondrial OXPHOS genes of three lineages of hydrothermal vent decapods, suggesting a link between hydrothermal vent adaption and OXPHOS molecular biology in decapods. This study adds to the understanding of the link between mitogenome evolution and ecological adaptation to hydrothermal vent habitats in decapods.

Salmonid gene expression biomarkers indicative of physiological responses to changes in salinity and temperature, but not dissolved oxygen

An organism's ability to respond effectively to environmental change is critical to its survival. Yet, life stage and overall condition can dictate tolerance thresholds to heightened environmental stressors, such that stress may not be equally felt across individuals and at all times. Also, the transcriptional responses induced by environmental changes can reflect both generalized responses as well as others that are highly specific to the type of change being experienced. Thus, if transcriptional biomarkers specific to a stressor, even under multi-stressor conditions, can be identified, the biomarkers could then be applied in natural environments to determine when and where an individual experiences such a stressor. Here, we experimentally challenged juvenile Chinook salmon (Oncorhynchus tshawytscha) to validate candidate gill gene expression biomarkers. A sophisticated experimental design manipulated salinity (freshwater, brackish water and seawater), temperature (10, 14 and 18°C) and dissolved oxygen (normoxia and hypoxia) in all 18 possible combinations for 6 days using separate trials for three smolt statuses (pre-smolt, smolt and de-smolt). In addition, changes in juvenile behaviour, plasma variables, gill Na+/K+-ATPase activity, body size, body morphology and skin pigmentation supplemented the gene expression responses. We identified biomarkers specific to salinity and temperature that transcended the multiple stressors, smolt status and mortality (live, dead and moribund). Similar biomarkers for dissolved oxygen were not identified. This work demonstrates the unique power of gene expression biomarkers to identify a specific stressor even under multi-stressor conditions, and we discuss our next steps for hypoxia biomarkers using an RNA-seq study.

Developing specific molecular biomarkers for thermal stress in salmonids

Background

Pacific salmon (Oncorhynchus spp.) serve as good biological indicators of the breadth of climate warming effects on fish because their anadromous life cycle exposes them to environmental challenges in both marine and freshwater environments. Our study sought to mine the extensive functional genomic studies in fishes to identify robust thermally-responsive biomarkers that could monitor molecular physiological signatures of chronic thermal stress in fish using non-lethal sampling of gill tissue.

Results

Candidate thermal stress biomarkers for gill tissue were identified using comparisons among microarray datasets produced in the Molecular Genetics Laboratory, Pacific Biological Station, Nanaimo, BC, six external, published microarray studies on chronic and acute temperature stress in salmon, and a comparison of significant genes across published studies in multiple fishes using deep literature mining. Eighty-two microarray features related to 39 unique gene IDs were selected as candidate chronic thermal stress biomarkers. Most of these genes were identified both in the meta-analysis of salmon microarray data and in the literature mining for thermal stress markers in salmonids and other fishes. Quantitative reverse transcription PCR (qRT-PCR) assays for 32 unique genes with good efficiencies across salmon species were developed, and their activity in response to thermally challenged sockeye salmon (O. nerka) and Chinook salmon (O. tshawytscha) (cool, 13–14 °C and warm temperatures 18–19 °C) over 5–7 days was assessed. Eight genes, including two transcripts of each SERPINH1 and HSP90AA1, FKBP10, MAP3K14, SFRS2, and EEF2 showed strong and robust chronic temperature stress response consistently in the discovery analysis and both sockeye and Chinook salmon validation studies.

Conclusions

The results of both discovery analysis and gene expression showed that a panel of genes involved in chaperoning and protein rescue, oxidative stress, and protein biosynthesis were differentially activated in gill tissue of Pacific salmon in response to elevated temperatures. While individually, some of these biomarkers may also respond to other stressors or biological processes, when expressed in concert, we argue that a biomarker panel comprised of some or all of these genes could provide a reliable means to specifically detect thermal stress in field-caught salmon.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5108-9) contains supplementary material, which is available to authorized users.

Stoichiometric responses to nano ZnO under warming are modified by thermal evolution in Daphnia magna

Effects of stressors on body stoichiometry are important as these may cascade through food webs. Contamination and global warming are two key anthropogenic stressors, yet their effects on body stoichiometry have been rarely tested. Further, while thermal evolution may increase the ability to deal with warming, it is unknown how thermal evolution modifies the effect of contaminants under warming. Using resurrection ecology, we studied two Daphnia magna subpopulations (old/recent) of which the recent subpopulation evolved a higher heat tolerance. We exposed both subpopulations to a sublethal concentration of nano zinc oxide (nZnO) and 4 °C warming and quantified their effects on body stoichiometry: carbon (C), nitrogen (N), phosphorus (P) contents and their ratios (C:N, C:P, N:P). In the old subpopulation, nZnO only marginally decreased the C content and had no effect on N and P contents and their ratios. In contrast, in the recent subpopulation nZnO strongly increased the body P content (+51%) and reduced the C:P (-34%) and N:P (-34%) ratios at 24 °C but not at 20 °C. Moreover, these stoichiometric changes were not explained by changes of corresponding macromolecules as assumed by theory. Our results indicate that the stoichiometric responses to nZnO in Daphnia are temperature-dependent and modified by rapid evolution. The observed changes in body stoichiometry may affect the food quality of this important prey and have the potential to cascade through food webs and shape nutrients cycling.

Thermal evolution offsets the elevated toxicity of a contaminant under warming: A resurrection study in Daphnia magna

Synergistic interactions between temperature and contaminants are a major challenge for ecological risk assessment, especially under global warming. While thermal evolution may increase the ability to deal with warming, it is unknown whether it may also affect the ability to deal with the many contaminants that are more toxic at higher temperatures. We investigated how evolution of genetic adaptation to warming affected the interactions between warming and a novel stressor: zinc oxide nanoparticles (nZnO) in a natural population of Daphnia magna using resurrection ecology. We hatched resting eggs from two D. magna subpopulations (old: 1955-1965, recent: 1995-2005) from the sediment of a lake that experienced an increase in average temperature and in recurrence of heat waves but was never exposed to industrial waste. In the old "ancestral" subpopulation, exposure to a sublethal concentration of nZnO decreased the intrinsic growth rate, metabolic activity, and energy reserves at 24°C but not at 20°C, indicating a synergism between warming and nZnO. In contrast, these synergistic effects disappeared in the recent "derived" subpopulation that evolved a lower sensitivity to nZnO at 24°C, which indicates that thermal evolution could offset the elevated toxicity of nZnO under warming. This evolution of reduced sensitivity to nZnO under warming could not be explained by changes in the total internal zinc accumulation but was partially associated with the evolution of the expression of a key metal detoxification gene under warming. Our results suggest that the increased sensitivity to the sublethal concentration of nZnO under the predicted 4°C warming by the end of this century may be counteracted by thermal evolution in this D. magna population. Our results illustrate the importance of evolution to warming in shaping the responses to another anthropogenic stressor, here a contaminant. More general, genetic adaptation to an environmental stressor may ensure that synergistic effects between contaminants and this environmental stressor will not be present anymore.

Mechanisms of thermal adaptation and evolutionary potential of conspecific populations to changing environments

Heterogeneous and ever-changing thermal environments drive the evolution of populations and species, especially when extreme conditions increase selection pressure for traits influencing fitness. However, projections of biological diversity under scenarios of climate change rarely consider evolutionary adaptive potential of natural species. In this study, we tested for mechanistic evidence of evolutionary thermal adaptation among ecologically divergent redband trout populations (Oncorhynchus mykiss gairdneri) in cardiorespiratory function, cellular response and genomic variation. In a common garden environment, fish from an extreme desert climate had significantly higher critical thermal maximum (p < .05) and broader optimum thermal window for aerobic scope (>3°C) than fish from cooler montane climate. In addition, the desert population had the highest maximum heart rate during warming (20% greater than montane populations), indicating improved capacity to deliver oxygen to internal tissues. In response to acute heat stress, distinct sets of cardiac genes were induced among ecotypes, which helps to explain the differences in cardiorespiratory function. Candidate genomic markers and genes underlying these physiological adaptations were also pinpointed, such as genes involved in stress response and metabolic activity (hsp40, ldh-b and camkk2). These markers were developed into a multivariate model that not only accurately predicted critical thermal maxima, but also evolutionary limit of thermal adaptation in these specific redband trout populations relative to the expected limit for the species. This study demonstrates mechanisms and limitations of an aquatic species to evolve under changing environments that can be incorporated into advanced models to predict ecological consequences of climate change for natural organisms.

Adaptive Patterns of Mitogenome Evolution Are Associated with the Loss of Shell Scutes in Turtles

The mitochondrial genome encodes several protein components of the oxidative phosphorylation (OXPHOS) pathway and is critical for aerobic respiration. These proteins have evolved adaptively in many taxa, but linking molecular-level patterns with higher-level attributes (e.g., morphology, physiology) remains a challenge. Turtles are a promising system for exploring mitochondrial genome evolution as different species face distinct respiratory challenges and employ multiple strategies for ensuring efficient respiration. One prominent adaptation to a highly aquatic lifestyle in turtles is the secondary loss of keratenized shell scutes (i.e., soft-shells), which is associated with enhanced swimming ability and, in some species, cutaneous respiration. We used codon models to examine patterns of selection on mitochondrial protein-coding genes along the three turtle lineages that independently evolved soft-shells. We found strong evidence for positive selection along the branches leading to the pig-nosed turtle (Carettochelys insculpta) and the softshells clade (Trionychidae), but only weak evidence for the leatherback (Dermochelys coriacea) branch. Positively selected sites were found to be particularly prevalent in OXPHOS Complex I proteins, especially subunit ND2, along both positively selected lineages, consistent with convergent adaptive evolution. Structural analysis showed that many of the identified sites are within key regions or near residues involved in proton transport, indicating that positive selection may have precipitated substantial changes in mitochondrial function. Overall, our study provides evidence that physiological challenges associated with adaptation to a highly aquatic lifestyle have shaped the evolution of the turtle mitochondrial genome in a lineage-specific manner.

Sex-specific influences of mtDNA mitotype and diet on mitochondrial functions and physiological traits in Drosophila melanogaster

Here we determine the sex-specific influence of mtDNA type (mitotype) and diet on mitochondrial functions and physiology in two Drosophila melanogaster lines. In many species, males and females differ in aspects of their energy production. These sex-specific influences may be caused by differences in evolutionary history and physiological functions. We predicted the influence of mtDNA mutations should be stronger in males than females as a result of the organelle's maternal mode of inheritance in the majority of metazoans. In contrast, we predicted the influence of diet would be greater in females due to higher metabolic flexibility. We included four diets that differed in their protein: carbohydrate (P:C) ratios as they are the two-major energy-yielding macronutrients in the fly diet. We assayed four mitochondrial function traits (Complex I oxidative phosphorylation, reactive oxygen species production, superoxide dismutase activity, and mtDNA copy number) and four physiological traits (fecundity, longevity, lipid content, and starvation resistance). Traits were assayed at 11 d and 25 d of age. Consistent with predictions we observe that the mitotype influenced males more than females supporting the hypothesis of a sex-specific selective sieve in the mitochondrial genome caused by the maternal inheritance of mitochondria. Also, consistent with predictions, we found that the diet influenced females more than males.

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Thermal plasticity in farmed, wild and hybrid Atlantic salmon during early development: has domestication caused divergence in low temperature tolerance?

BACKGROUND

In the past three decades, millions of domesticated Atlantic salmon Salmo salar L. have escaped from farms into the wild. Their offspring display reduced survival in the natural environment, which demonstrates that gene-flow is likely to have a negative effect on wild populations. However, inter-population differences in introgression of farmed salmon have been observed, and the underlying ecological mechanisms remain enigmatic. We hypothesised that domestication-driven divergence in tolerance to low temperatures during early development may contribute to lower survival of farmed salmon offspring in the wild, which in turn, may influence patterns of introgression among populations exposed to different temperature regimes. We reared the offspring of 35 families of wild, farmed and hybrid origin at three temperatures (3.9, 5.6 and 12°C) from the onset of exogenous feeding and throughout their first summer. Thermal reaction norms for growth and survival were investigated along the gradient.

RESULTS

The main results of this study, which is based upon the analysis of juvenile salmon from five wild strains, two farmed strains and two hybrid strains, can be summarised as; (i) salmon of all origins were able to successfully initiate feeding at all temperatures and similar survival reaction norms were detected in all strains across the temperature gradient; (ii) deviating growth reaction norms were detected between strains, although this result was most likely due to an overall lack of growth in the lower temperature treatments.

CONCLUSIONS

This study revealed no evidence of domesticated-driven divergence in low temperature tolerance in Atlantic salmon during early development. Although the potential interaction between low temperature and other river-specific factors cannot be excluded, our results indicate that the reduced survival of farmed offspring in the wild is not explained by farmed salmon displaying impaired abilities to initiate feeding at low temperatures. We therefore suggest that the observed inter-population patterns of introgression are not low-temperature driven and that other ecological or biological factors may explain why detection of farmed salmon in wild rivers is not synonymous with introgression. In general, our results support the literature indicating that phenotypic plasticity instead of thermal adaption has been selected for in Atlantic salmon.