Integrated organismal responses induced by projected levels of CO2 and temperature exposures in the early life stages of lake sturgeon.

Atmospheric CO2 and temperature are rising concurrently, and may have profound impacts on the transcriptional, physiological and behavioural responses of aquatic organisms. Further, spring snowmelt may cause transient increases of pCO2 in many freshwater systems. We examined the behavioural, physiological and transcriptomic responses of an ancient fish, the lake sturgeon (Acipenser fulvescens) to projected levels of warming and pCO2 during its most vulnerable period of life, the first year. Specifically, larval fish were raised in either low (16°C) or high (22°C) temperature, and/or low (1000 µatm) or high (2500 µatm) pCO2 in a crossed experimental design over approximately 8 months. Following overwintering, lake sturgeon were exposed to a transient increase in pCO2 of 10,000 µatm, simulating a spring melt based on data in freshwater systems. Transcriptional analyses revealed potential connections to otolith formation and reduced growth in fish exposed to high pCO2 and temperature in combination. Network analyses of differential gene expression revealed different biological processes among the different treatments on the edges of transcriptional networks. Na+/K+-ATPase activity increased in fish not exposed to elevated pCO2 during development, and mRNA abundance of the ß subunit was most strongly predictive of enzyme activity. Behavioural assays revealed a decrease in total activity following an acute CO2 exposure. These results demonstrate compensatory and compounding mechanisms of pCO2 and warming dependent on developmental conditions in lake sturgeon. Conserved elements of the cellular stress response across all organisms provide key information for how other freshwater organisms may respond to future climate change.

Climate change introduces threatened killer whale populations and conservation challenges to the Arctic.

The Arctic is the fastest-warming region on the planet, and the lengthening ice-free season is opening Arctic waters to sub-Arctic species such as the killer whale (Orcinus orca). As apex predators, killer whales can cause significant ecosystem-scale changes. Setting conservation priorities for killer whales and their Arctic prey species requires knowledge of their evolutionary history and demographic trajectory. Using whole-genome resequencing of 24 killer whales sampled in the northwest Atlantic, we first explored the population structure and demographic history of Arctic killer whales. To better understand the broader geographic relationship of these Arctic killer whales to other populations, we compared them to a globally sampled dataset. Finally, we assessed threats to Arctic killer whales due to anthropogenic harvest by reviewing the peer-reviewed and gray literature. We found that there are two highly genetically distinct, non-interbreeding populations of killer whales using the eastern Canadian Arctic. These populations appear to be as genetically different from each other as are ecotypes described elsewhere in the killer whale range; however, our data cannot speak to ecological differences between these populations. One population is newly identified as globally genetically distinct, and the second is genetically similar to individuals sampled from Greenland. The effective sizes of both populations recently declined, and both appear vulnerable to inbreeding and reduced adaptive potential. Our survey of human-caused mortalities suggests that harvest poses an ongoing threat to both populations. The dynamic Arctic environment complicates conservation and management efforts, with killer whales adding top-down pressure on Arctic food webs crucial to northern communities' social and economic well-being. While killer whales represent a conservation priority, they also complicate decisions surrounding wildlife conservation and resource management in the Arctic amid the effects of climate change.

Comparison of metabolic rate between two genetically distinct populations of lake sturgeon.

Environmental temperatures differ across latitudes in the temperate zone, with relatively lower summer and fall temperatures in the north leading to a shorter growing season prior to winter. As an adaptive response, during early life stages, fish in northern latitudes may grow faster than their conspecifics in southern latitudes, which potentially manifests as different allometric relationships between body mass and metabolic rate. In the present study, we examined if population or year class had an effect on the variation of metabolic rate and metabolic scaling of age-0 lake sturgeon (Acipenser fulvescens) by examining these traits in both a northern (Nelson River) and a southern (Winnipeg River) population. We compiled 6 years of data that used intermittent flow respirometry to measure metabolic rate within the first year of life for developing sturgeon that were raised in the same environment at 16°C. We then used a Bayesian modeling approach to examine the impacts of population and year class on metabolic rate and mass-scaling of metabolic rate. Despite previous reports of genetic differences between populations, our results showed that there were no significant differences in standard metabolic rate, routine metabolic rate, maximum metabolic rate, and metabolic scaling between the two geographically separated populations at a temperature of 16°C. Our analysis implied that the lack of metabolic differences between populations could be due to family effects/parental contribution, or the rearing temperature used in the study. The present research provided insights for conservation and reintroduction strategies for these populations of lake sturgeon, which are endangered or threatened across most of their natural range.

Tissue-specific transcriptomes reveal potential mechanisms of microbiome heterogeneity in an ancient fish.

The lake sturgeon (Acipenser fulvescens) is an ancient, octoploid fish faced with conservation challenges across its range in North America, but a lack of genomic resources has hindered molecular research in the species. To support such research, we created a transcriptomic database from 13 tissues: brain, esophagus, gill, head kidney, heart, white muscle, liver, glandular stomach, muscular stomach, anterior intestine, pyloric cecum, spiral valve and rectum. The transcriptomes for each tissue were sequenced and assembled individually from a mean of 98.3 million (±38.9 million SD) reads each. In addition, an overall transcriptome was assembled and annotated with all data used for each tissue-specific transcriptome. All assembled transcriptomes and their annotations were made publicly available as a scientific resource. The non-gut transcriptomes provide important resources for many research avenues. However, we focused our analysis on messenger ribonucleic acid (mRNA) observations in the gut because the gut represents a compartmentalized organ system with compartmentalized functions, and seven of the sequenced tissues were from each of these portions. These gut-specific analyses were used to probe evidence of microbiome regulation by studying heterogeneity in microbial genes and genera identified from mRNA annotations. Gene set enrichment analyses were used to reveal the presence of photoperiod and circadian-related transcripts in the pyloric cecum, which may support periodicity in lake sturgeon digestion. Similar analyses were used to identify different types of innate immune regulation across the gut, while analyses of unique transcripts annotated to microbes revealed heterogeneous genera and genes among different gut tissues. The present results provide a scientific resource and information about the mechanisms of compartmentalized function across gut tissues in a phylogenetically ancient vertebrate. Database URL: https://figshare.com/projects/Lake_Sturgeon_Transcriptomes/133143.

Elevated temperatures reduce population-specific transcriptional plasticity in developing lake sturgeon (Acipenser fulvescens).

Rising mean and variance in temperatures elevates threats to endangered freshwater species such as lake sturgeon, Acipenser fulvescens. Previous research demonstrated that higher temperatures during development result in physiological consequences for lake sturgeon populations throughout Manitoba, Canada, with alteration of metabolic rate, thermal tolerance, transcriptional responses, growth and mortality. We acclimated lake sturgeon (30-60 days post fertilization, a period of high mortality) from northern and southern populations (56°02'46.5″N, 96°54'18.6″W and 50°17'52″N, 95°32'51″W, respectively, separated by approximately 650 km) within Manitoba to current (summer highs of 20-23°C) and future projected (+2-3°C) environmental temperatures of 16, 20 and 24°C for 30 days, and we measured gill transcriptional responses using RNAseq. Transcripts revealed SNPs consistent with genetically distinct populations and transcriptional responses altered by acclimation temperature. There were a higher number of differentially expressed transcripts observed in the southern, compared to the northern, population as temperatures increased, indicating enhanced transcriptional plasticity. Both lake sturgeon populations responded to elevated acclimation temperatures by downregulating the transcription of genes involved in protein synthesis and energy production. Furthermore, there were population-specific thresholds for the downregulation of processes promoting transcriptional plasticity as well as mitochondrial function as the northern population showed decreases at 20°C, while this capacity was not diminished until 24°C in the southern population. These transcriptional responses highlight the molecular impacts of increasing temperatures for divergent lake sturgeon populations during vulnerable developmental periods and the critical influence of transcriptome plasticity on acclimation capacity.

Genomic architecture of migration timing in a long-distance migratory songbird.

The impact of climate change on spring phenology poses risks to migratory birds, as migration timing is controlled predominantly by endogenous mechanisms. Despite recent advances in our understanding of the underlying genetic basis of migration timing, the ways that migration timing phenotypes in wild individuals may map to specific genomic regions requires further investigation. We examined the genetic architecture of migration timing in a long-distance migratory songbird (purple martin, Progne subis subis) by integrating genomic data with an extensive dataset of direct migratory tracks. A moderate to large amount of variance in spring migration arrival timing was explained by genomics (proportion of phenotypic variation explained by genomics = 0.74; polygenic score R2 = 0.24). On chromosome 1, a region that was differentiated between migration timing phenotypes contained genes that could facilitate nocturnal flights and act as epigenetic modifiers. Overall, these results advance our understanding of the genomic underpinnings of migration timing.

Acute thermal stress elicits interactions between gene expression and alternative splicing in a fish of conservation concern.

Transcriptomic research provides a mechanistic understanding of an organism's response to environmental challenges such as increasing temperatures, which can provide key insights into the threats posed by thermal challenges associated with urbanization and climate change. Differential gene expression and alternative splicing are two elements of the transcriptomic stress response that may work in tandem, but relatively few studies have investigated these interactions in fishes of conservation concern. We studied the imperilled redside dace (Clinostomus elongatus) as thermal stress is hypothesized to be an important cause of population declines. We tested the hypothesis that gene expression-splicing interactions contribute to the thermal stress response. Wild fish exposed to acute thermal stress were compared with both handling controls and fish sampled directly from a river. Liver tissue was sampled to study the transcriptomic stress response. With a gene set enrichment analysis, we found that thermally stressed fish showed a transcriptional response related to transcription regulation and responses to unfolded proteins, and alternatively spliced genes related to gene expression regulation and metabolism. One splicing factor, prpf38b, was upregulated in the thermally stressed group compared with the other treatments. This splicing factor may have a role in the Jun/AP-1 cellular stress response, a pathway with wide-ranging and context-dependent effects. Given large gene interaction networks and the context-dependent nature of transcriptional responses, our results highlight the importance of understanding interactions between gene expression and splicing for understanding transcriptomic responses to thermal stress. Our results also reveal transcriptional pathways that can inform conservation breeding, translocation and reintroduction programs for redside dace and other imperilled species by identifying appropriate source populations.

A chromosomal inversion may facilitate adaptation despite periodic gene flow in a freshwater fish.

Differences in genomic architecture between populations, such as chromosomal inversions, may play an important role in facilitating adaptation despite opportunities for gene flow. One system where chromosomal inversions may be important for eco-evolutionary dynamics is in freshwater fishes, which often live in heterogenous environments characterized by varying levels of connectivity and varying opportunities for gene flow. In the present study, reduced representation sequencing was used to study possible adaptation in n = 345 walleye (Sander vitreus) from three North American waterbodies: Cedar Bluff Reservoir (Kansas, USA), Lake Manitoba (Manitoba, Canada), and Lake Winnipeg (Manitoba, Canada). Haplotype and outlier-based tests revealed a putative chromosomal inversion that contained three expressed genes and was nearly fixed in walleye assigned to Lake Winnipeg. These patterns exist despite the potential for high gene flow between these proximate Canadian lakes, suggesting that the inversion may be important for facilitating adaptive divergence between the two lakes despite gene flow. However, a specific adaptive role for the putative inversion could not be tested with the present data. Our study illuminates the importance of genomic architecture consistent with local adaptation in freshwater fishes. Furthermore, our results provide additional evidence that inversions may facilitate local adaptation in many organisms that inhabit connected but heterogenous environments.

Non-Lethal Sampling Supports Integrative Movement Research in Freshwater Fish.

Freshwater ecosystems and fishes are enormous resources for human uses and biodiversity worldwide. However, anthropogenic climate change and factors such as dams and environmental contaminants threaten these freshwater systems. One way that researchers can address conservation issues in freshwater fishes is via integrative non-lethal movement research. We review different methods for studying movement, such as with acoustic telemetry. Methods for connecting movement and physiology are then reviewed, by using non-lethal tissue biopsies to assay environmental contaminants, isotope composition, protein metabolism, and gene expression. Methods for connecting movement and genetics are reviewed as well, such as by using population genetics or quantitative genetics and genome-wide association studies. We present further considerations for collecting molecular data, the ethical foundations of non-lethal sampling, integrative approaches to research, and management decisions. Ultimately, we argue that non-lethal sampling is effective for conducting integrative, movement-oriented research in freshwater fishes. This research has the potential for addressing critical issues in freshwater systems in the future.

RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish.

BACKGROUND: Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well-including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis. RESULTS: We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13). CONCLUSIONS: These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem.

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.