De novo assembly and characterization of tissue specific transcriptomes in the emerald notothen, Trematomus bernacchii

The hypoxia response pathway in the Antarctic fish Notothenia coriiceps is functional despite a poly Q/E insertion mutation in HIF-1α

Antarctic notothenioid fishes, inhabiting the oxygen-rich Southern Ocean, possess a polyglutamine and glutamic acid (poly Q/E) insertion mutation in the master transcriptional regulator of oxygen homeostasis, hypoxia- inducible factor-1α (HIF-1α). To determine if this mutation impairs the ability of HIF-1 to regulate gene expression in response to hypoxia, we exposed Notothenia coriiceps, with a poly Q/E insertion mutation in HIF-1α that is 9 amino acids long, to hypoxia (2.3 mg L-1 O2) or normoxia (10 mg L -1 O2) for 12 h. Heart ventricles, brain, liver, and gill tissue were harvested and changes in gene expression quantified using RNA sequencing. Levels of glycogen and lactate were also quantified to determine if anaerobic metabolism increases in response to hypoxia. Exposure to hypoxia resulted in 818 unique differentially expressed genes (DEGs) in liver tissue of N. coriiceps. Many hypoxic genes were induced, including ones involved in the MAP kinase and FoxO pathways, glycolytic metabolism, and vascular remodeling. In contrast, there were fewer than 104 unique DEGs in each of the other tissues sampled. Lactate levels significantly increased in liver in response to hypoxia, indicating that anaerobic metabolism increases in response to hypoxia in this tissue. Overall, our results indicate that the hypoxia response pathway is functional in N. coriiceps despite a poly Q/E mutation in HIF-1α, and confirm that Antarctic fishes are capable of altering gene expression in response to hypoxia.

Tissue-specific transcriptional response of post-larval clownfish to ocean warming

Anthropogenically driven climate change is predicted to increase average sea surface temperatures, as well as the frequency and intensity of marine heatwaves in the future. This increasing temperature is predicted to have a range of negative physiological impacts on multiple life-stages of coral reef fish. Nevertheless, studies of early-life stages remain limited, and tissue-specific transcriptomic studies of post-larval coral reef fish are yet to be conducted. Here, in an aquaria-based study we investigate the tissue-specific (brain, liver, muscle, and digestive tract) transcriptomic response of post-larval (20 dph) Amphiprion ocellaris to temperatures associated with future climate change (+3 °C). Additionally, we utilized metatranscriptomic sequencing to investigate how the microbiome of the digestive tract changes at +3 °C. Our results show that the transcriptional response to elevated temperatures is highly tissue-specific, as the number of differentially expressed genes (DEGs) and gene functions varied amongst the brain (102), liver (1785), digestive tract (380), and muscle (447). All tissues displayed DEGs associated with thermal stress, as 23 heat-shock protein genes were upregulated in all tissues. Our results indicate that post-larval clownfish may experience liver fibrosis-like symptoms at +3 °C as genes associated with extracellular matrix structure, oxidative stress, inflammation, glucose transport, and metabolism were all upregulated. We also observe a shift in the digestive tract microbiome community structure, as Vibrio sp. replace Escherichia coli as the dominant bacteria. This shift is coupled with the dysregulation of various genes involved in immune response in the digestive tract. Overall, this study highlights post-larval clownfish will display tissue-specific transcriptomic responses to future increases in temperature, with many potentially harmful pathways activated at +3 °C.

Antibacterial and anticancer activity of two NK-lysin-derived peptides from the Antarctic teleost Trematomus bernacchii

The NK-lysin antimicrobial peptide, first identified in mammals, possesses both antibacterial and cytotoxic activity against cancer cell lines. Homologue peptides isolated from different fish species have been examined for their functional characteristics in the last few years. In this study, a NK-lysin transcript was identified in silico from the head kidney transcriptome of the Antarctic teleost Trematomus bernacchii. The corresponding amino acid sequence, slightly longer than NK-lysins of other fish species, contains six cysteine residues that in mammalian counterparts form three disulphide bridges. Real time-PCR analysis indicated its predominant expression in T. bernacchii immune-related organs and tissues, with greatest mRNA abundance detected in gills and spleen. Instead of focusing on the full T. bernacchii derived NK-lysin mature molecule, we selected a 27 amino acid residue peptide (named NKL-WT), corresponding to the potent antibiotic NK-2 sequence found in human NK-lysin. Moreover, we designed a mutant peptide (named NKL-MUT) in which two alanine residues substitute the two cysteines found in the NKL-WT. The two peptides were obtained by solid phase organic synthesis to investigate their functional features. NKL-WT and NKL-MUT displayed antibacterial activity against the human pathogenic bacterium Enterococcus faecalis and the ESKAPE pathogen Acinetobacter baumannii, respectively. Moreover, at the determined Minimum Inhibitory Concentration and Minimum Bactericidal Concentration values against these pathogens, both peptides showed high selectivity as they did not exhibit any haemolytic activity on erythrocytes or cytotoxic activity against mammalian primary cell lines. Finally, the NKL-MUT selectively triggers the killing of the melanoma cell line B16F10 by means of a pro-apoptotic pathway at a concentration range in which no effects were found in normal mammalian cell lines. In conclusion, the two peptides could be considered as promising candidates in the fight against antibiotic resistance and tumour proliferative action, and also be used as innovative adjuvants, either to decrease chemotherapy side effects or to enhance anticancer drug activity.

Transcriptome analysis and gene expression analysis related to salinity-alkalinity and low temperature adaptation of Triplophysa yarkandensis

T. yarkandensis is a common species of Triplophysa, and it is distributed in Shule river of Hexi Corridor, of Gansu province in China. In order to enrich gene database resources and explore the environment adaptation of T. yarkandensis, fifteen tissues were collected from three adult T. yarkandensis for transcriptome sequencing and de novo assembly. Nine major international gene annotation databases (NR, COG, egg_NOG, TrEMBL, Pfam, KOG, Swiss prot, KEGG and Gene Ontology) were utilized to annotate unigenes. A detailed study was conducted to explore the gene expression and the differentially expressed genes among five tissues (brain, heart, kidney, liver and spleen). In addition, the current study showed that candidate genes involved in salinity-alkalinity and low temperature adaptation were differentially expressed in tissues of T. yarkandensis. Precisely, mapk1, abcc1, gpx1, gpx4, cat and aqp1 genes participated in the regulation process of salinity-alkalinity adaptation, and elovl4, acaca, fasn, acaa2, acox1 and acox3 genes were involved in fatty acid metabolism and closely associated with low temperature adaptation. On the one hand, it was found that the expression of these genes varied among different tissues, and the important pathways involved in these genes were mapped. Furthermore, we analyzed mapk1 and acox1 genes in depth to obtain the predicted gene structure and important amino acid sites. The transcriptome information in this study will be conducive to provide further understanding for the molecular level research and exploration of the environmental adaptation of T. yarkandensis.

Morphological analysis of erythrocytes of an Antarctic teleost under heat stress: Bias of the stabling effect

The stenothermal Antarctic fish that live in the coastal waters of the Terra Nova Bay (Ross Sea) are rarely exposed to temperatures above zero during the year. We tested whether a slight temperature rise of 1.5 °C affects a sensitive biomarker such as erythrocytes morphology in sections of blood pellets of a small demersal notothen. The erythrocytes' shape descriptors showed significant or highly significant differences temporally from the capture of fish to the conclusion of the experiment. Surprisingly, the erythrocyte's morphology did not show significant differences between the two experimental conditions, returning similar results in control fish stabled at -0.9 °C and in the fish treated at +0.6 °C, although the values of the shape descriptors were often lower in the latter. This study demonstrates the critical issues of comparative physiology in the study of extremely sensitive organisms, such as the fish of the High Antarctic Zone. Moreover, the stabling effect inside the aquarium facilities appears to significantly obscure the effects of the experimental heat treatment.

Regulation of globin expression in Antarctic fish under thermal and hypoxic stress

In the freezing waters of the Southern Ocean, Antarctic teleost fish, the Notothenioidei, have developed unique adaptations to cope with cold, including, at the extreme, the loss of hemoglobin in icefish. As a consequence, icefish are thought to be the most vulnerable of the Antarctic fish species to ongoing ocean warming. Some icefish also fail to express myoglobin but all appear to retain neuroglobin, cytoglobin-1, cytoglobin-2, and globin-X. Despite the lack of the inducible heat shock response, Antarctic notothenioid fish are endowed with physiological plasticity to partially compensate for environmental changes, as shown by numerous physiological and genomic/transcriptomic studies over the last decade. However, the regulatory mechanisms that determine temperature/oxygen-induced changes in gene expression remain largely unexplored in these species. Proteins such as globins are susceptible to environmental changes in oxygen levels and temperature, thus playing important roles in mediating Antarctic fish adaptations. In this study, we sequenced the full-length transcripts of myoglobin, neuroglobin, cytoglobin-1, cytoglobin-2, and globin-X from the Antarctic red-blooded notothenioid Trematomus bernacchii and the white-blooded icefish Chionodraco hamatus and evaluated transcripts levels after exposure to high temperature and low oxygen levels. Basal levels of globins are similar in the two species and both stressors affect the expression of Antarctic fish globins in brain, retina and gills. Temperature up-regulates globin expression more effectively in white-blooded than in red-blooded fish while hypoxia strongly up-regulates globins in red-blooded fish, particularly in the gills. These results suggest globins function as regulators of temperature and hypoxia tolerance. This study provides the first insights into globin transcriptional changes in Antarctic fish.

Characterizing Gene Copy Number of Heat Shock Protein Gene Families in the Emerald Rockcod, Trematomus bernacchii

The suborder Notothenioidae is comprised of Antarctic fishes, several of which have lost their ability to rapidly upregulate heat shock proteins in response to thermal stress, instead adopting a pattern of expression resembling constitutive genes. Given the cold-denaturing effect that sub-zero waters have on proteins, evolution in the Southern Ocean has likely selected for increased expression of molecular chaperones. These selective pressures may have also enabled retention of gene duplicates, bolstering quantitative output of cytosolic heat shock proteins (HSPs). Given that newly duplicated genes are under more relaxed selection, it is plausible that gene duplication enabled altered regulation of such highly conserved genes. To test for evidence of gene duplication, copy number of various isoforms within major heat shock gene families were characterized via qPCR and compared between the Antarctic notothen, Trematomus bernacchii, which lost the inducible heat shock response, and the non-Antarctic notothen, Notothenia angustata, which maintains an inducible heat shock response. The results indicate duplication of isoforms within the hsp70 and hsp40 super families have occurred in the genome of T. bernacchii. The findings suggest gene duplications may have been critical in maintaining protein folding efficiency in the sub-zero waters and provided an evolutionary mechanism of alternative regulation of these conserved gene families.

Molecular characterization of novel mitochondrial peroxiredoxins from the Antarctic emerald rockcod and their gene expression in response to environmental warming

In the present study we describe the molecular characterization of the two paralogous mitochondrial peroxiredoxins from Trematomus bernacchii, a teleost that plays a pivotal role in the Antarctic food chain. The two putative amino acid sequences were compared with orthologs from other fish, highlighting a high percentage of identity and similarity with the respective variant, in particular for the residues that are essential for the characteristic peroxidase activity of these enzymes. The temporal expression of Prdx3 and Prdx5 mRNAs in response to short-term thermal stress showed a general upregulation of prdx3, suggesting that this isoform is the most affected by temperature increase. These data, together with the peculiar differences between the molecular structures of the two mitochondrial Prdxs in T. bernacchii as well as in the tropical species Stegastes partitus, suggest an adaptation that allowed these poikilothermic aquatic vertebrates to colonize very different environments, characterized by different temperature ranges.

Transcriptome analyses reveal alterations in muscle metabolism, immune responses and reproductive behavior of Japanese mantis shrimp (Oratosquilla oratoria) at different cold temperature

Low temperature reduction is thought to cause widespread effects on the physical and behavioral traits of marine organisms, which include metabolic processes, immune responses, and reproductive behavior. Crustaceans are generally considered sensitive to temperature reduction due to the lack of efficient regulators. To better understand the molecular regulatory mechanisms of crustacean exposure to cold stress, Japanese mantis shrimp (Oratosquilla oratoria) was chosen as a representative crustacean. Transcriptomic responses in O. oratoria from five temperatures (25 °C, 22 °C, 19 °C, 16 °C, and 13 °C) were studied using RNA-seq. A total of 64.91 Gb of clean transcriptomic data were generated in 10 libraries and then spliced into 52,107 unigenes with an average length of 1089 bp and an N50 length of 1872 bp. A total of 14,841 unigenes was annotated in at least one database using Blastx alignment. Compared with the control temperature (25 °C), 7, 21, 58, and 236 unigenes were significantly differentially expressed at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. GO analysis showed that 6, 20, 27, and 35 terms were significantly enriched at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. In addition, 2, 5, 2, and 10 significant pathways were presented at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. Combining NR, GO, and KEGG annotation information, many genes significantly differentially expressed at low temperatures may be associated with metabolic processes, immune response, and reproductive behavior. Additionally, we reconstructed the phylogenetic relationship based on 366 orthologous genes and the predicted differentiation time of O. oratoria and P. vannamei range from 212.82 to 365.30 Mya. Furthermore, 16 orthologous genes were identified as PSGs and 30 orthologous genes were identified as FEGs and these adaptive genes were associated with energy metabolism, stress response and immunity, and multiple cellular processing. These results provide fundamental information about molecular mechanisms regulating cold stress response of O. oratoria.

Characterization of thermally sensitive miRNAs reveals a central role of the FoxO signaling pathway in regulating the cellular stress response of an extreme stenotherm, Trematomus bernacchii

Despite the lack of an inducible heat shock response (HSR), the Antarctic notothenioid fish, Trematomus bernacchii, has retained a level of physiological plasticity that can at least partially compensate for the effects of acute heat stress. Over the last decade, both physiological and transcriptomic studies have signaled these fish can mitigate the effects of acute heat stress by employing other aspects of the cellular stress response (CSR) that help confer thermotolerance as well as drive homeostatic mechanisms during long-term thermal acclimations. However, the regulatory mechanisms that determine temperature-induced changes in gene expression remain largely unexplored in this species. Therefore, this study utilized next generation sequencing coupled with an in silico approach to explore the regulatory role of microRNAs in governing the transcriptomic level response observed in this Antarctic notothenioid with respect to the CSR. Using RNAseq, we characterized the expression of 125 distinct miRNA orthologues in T. bernacchii gill tissue. Additionally, we identified 12 miRNAs that appear to be thermally responsive based on differential expression (DE) analyses performed between fish acclimated to control (^-1.5 °C) and an acute heat stress (⁺4 °C). We further characterized the functional role of these DE miRNAs using bioinformatics pipelines to identify putative gene targets of the DE miRNAs and subsequent gene set enrichment analyses, which together suggest these miRNAs are involved in regulating diverse aspects of the CSR in T. bernacchii.

Metabolic responses in Antarctic Nototheniidae brains subjected to thermal stress

Antarctic Nototheniidae is an attractive group for studying metabolic and physiological responses at high temperatures. The present work investigated the metabolic responses of the carbohydrate metabolism and antioxidant system to thermal stress at 8 °C (for 2-144 h) in the brains of Notothenia rossii and Notothenia coriiceps. In N. coriiceps, glycogenolysis was essential in the first hours of exposure (2 h) at 8 °C and, in addition to inhibiting glucose-6-phosphatase activity, was important for activating the pentose phosphate pathway. In N. rossii, anaerobic metabolism was reduced in the first hours of exposure (2 and 6 h) at 8 °C, followed by reduced hexokinase activity, suggesting energy regulation between neurons and astrocytes. The antioxidant system results indicated the importance of the actions of the glutathione-dependent antioxidant enzymes glutathione-S-transferase and glutathione peroxidase as well as those of catalase in N. coriiceps and the action of glutathione-S-transferase, glutathione peroxidase and glutathione reductase in N. rossii, especially during the first 12 h of thermal stress exposure. These results indicate tissue-specific patterns and species-specific responses to this stress.

Adaptation of Proteins to the Cold in Antarctic Fish: A Role for Methionine?

The evolution of antifreeze glycoproteins has enabled notothenioid fish to flourish in the freezing waters of the Southern Ocean. Whereas successful at the biodiversity level to life in the cold, paradoxically at the cellular level these stenothermal animals have problems producing, folding, and degrading proteins at their ambient temperatures of -1.86 °C. In this first multi-species transcriptome comparison of the amino acid composition of notothenioid proteins with temperate teleost proteins, we show that, unlike psychrophilic bacteria, Antarctic fish provide little evidence for the mass alteration of protein amino acid composition to enhance protein folding and reduce protein denaturation in the cold. The exception was the significant overrepresentation of positions where leucine in temperate fish proteins was replaced by methionine in the notothenioid orthologues. We hypothesize that these extra methionines have been preferentially assimilated into the genome to act as redox sensors in the highly oxygenated waters of the Southern Ocean. This redox hypothesis is supported by analyses of notothenioids showing enrichment of genes associated with responses to environmental stress, particularly reactive oxygen species. So overall, although notothenioid fish show cold-associated problems with protein homeostasis, they may have modified only a selected number of biochemical pathways to work efficiently below 0 °C. Even a slight warming of the Southern Ocean might disrupt the critical functions of this handful of key pathways with considerable impacts for the functioning of this ecosystem in the future.

Hypoxia-Inducible Factor-1α in Antarctic notothenioids contains a polyglutamine and glutamic acid insert that varies in length with phylogeny

The long evolution of the Antarctic perciform suborder of Notothenioidei in the icy, oxygen-rich waters of the Southern Ocean may have reduced selective pressure to maintain a hypoxic response. To test this hypothesis, cDNA of the key transcriptional regulator of hypoxic genes, hypoxia-inducible factor-1α (HIF-1α), was sequenced in heart ventricles of the red-blooded notothenioid, Notothenia coriiceps, and the hemoglobinless icefish, Chaenocephalus aceratus. HIF-1α cDNA is 4500 base pairs (bp) long and encodes 755 amino acids in N. coriiceps, and in C. aceratus, HIF-1α is 3576 bp long and encodes 779 amino acids. All functional domains of HIF-1α are highly conserved compared to other teleosts, but HIF-1α contains a polyglutamine/glutamic acid (polyQ/E) insert 9 amino acids long in N. coriiceps and 34 amino acids long in C. aceratus. Sequencing of this region in four additional species, representing three families of notothenioids, revealed that the length of the polyQ/E insert varies with phylogeny. Icefishes, the crown family of notothenioids, contain the longest polyQ/E inserts, ranging between16 and 34 amino acids long, whereas the basal, cold-temperate notothenioid, Eleginops maclovinus, contains a polyQ/E insert only 4 amino acids long. PolyQ/E inserts may affect dimerization of HIF-1α and HIF-1β, HIF-1 translocation into the nucleus and/or DNA binding.

Antarctic emerald rockcod have the capacity to compensate for warming when uncoupled from CO2 -acidification

Increases in atmospheric CO₂ levels and associated ocean changes are expected to have dramatic impacts on marine ecosystems. Although the Southern Ocean is experiencing some of the fastest rates of change, few studies have explored how Antarctic fishes may be affected by co-occurring ocean changes, and even fewer have examined early life stages. To date, no studies have characterized potential trade-offs in physiology and behavior in response to projected multiple climate change stressors (ocean acidification and warming) on Antarctic fishes. We exposed juvenile emerald rockcod Trematomus bernacchii to three PCO₂ treatments (~450, ~850, and ~1,200 μatm PCO₂ ) at two temperatures (-1 or 2°C). After 2, 7, 14, and 28 days, metrics of physiological performance including cardiorespiratory function (heart rate [f_H ] and ventilation rate [f_V ]), metabolic rate (M˙O2), and cellular enzyme activity were measured. Behavioral responses, including scototaxis, activity, exploration, and escape response were assessed after 7 and 14 days. Elevated PCO₂ independently had little impact on either physiology or behavior in juvenile rockcod, whereas warming resulted in significant changes across acclimation time. After 14 days, f_H , f_V and M˙O2 significantly increased with warming, but not with elevated PCO₂ . Increased physiological costs were accompanied by behavioral alterations including increased dark zone preference up to 14%, reduced activity by 12%, as well as reduced escape time suggesting potential trade-offs in energetics. After 28 days, juvenile rockcod demonstrated a degree of temperature compensation as f_V , M˙O2, and cellular metabolism significantly decreased following the peak at 14 days; however, temperature compensation was only evident in the absence of elevated PCO₂ . Sustained increases in f_V and M˙O2 after 28 days exposure to elevated PCO₂ indicate additive (f_V ) and synergistic (M˙O2) interactions occurred in combination with warming. Stressor-induced energetic trade-offs in physiology and behavior may be an important mechanism leading to vulnerability of Antarctic fishes to future ocean change.

Cold Fusion: Massive Karyotype Evolution in the Antarctic Bullhead Notothen Notothenia coriiceps

Half of all vertebrate species share a series of chromosome fusions that preceded the teleost genome duplication (TGD), but we do not understand the causative evolutionary mechanisms. The "Robertsonian-translocation hypothesis" suggests a regular fusion of each ancestral acro- or telocentric chromosome to just one other by centromere fusions, thus halving the karyotype. An alternative "genome-stirring hypothesis" posits haphazard and repeated fusions, inversions, and reciprocal and nonreciprocal translocations. To study large-scale karyotype reduction, we investigated the decrease of chromosome numbers in Antarctic notothenioid fish. Most notothenioids have 24 haploid chromosomes, but bullhead notothen (Notothenia coriiceps) has 11. To understand mechanisms, we made a RAD-tag meiotic map with ∼10,000 polymorphic markers. Comparative genomics aligned about a thousand orthologs of platyfish and stickleback genes along bullhead chromosomes. Results revealed that 9 of 11 bullhead chromosomes arose by fusion of just two ancestral chromosomes and two others by fusion of three ancestral chromosomes. All markers from each ancestral chromosome remained contiguous, implying no inversions across fusion borders. Karyotype comparisons support a history of: (1) Robertsonian fusions of 22 ancestral chromosomes in pairs to yield 11 fused plus two small unfused chromosomes, like N. angustata; (2) fusion of one of the remaining two ancestral chromosomes to a preexisting fused pair, giving 12 chromosomes like N. rossii; and (3) fusion of the remaining ancestral chromosome to another fused pair, giving 11 chromosomes in N. coriiceps These results raise the question of what selective forces promoted the systematic fusion of chromosomes in pairs and the suppression of pericentric inversions in this lineage, and provide a model for chromosome fusions in stem teleosts.

Loggerhead sea turtle embryos (Caretta caretta) regulate expression of stress response and developmental genes when exposed to a biologically realistic heat stress

Oviparous reptile embryos are expected to breach their critical thermal maxima if temperatures reach those predicted under current climate change models due to the lack of the maternal buffering processes and parental care. Heat-shock proteins (HSPs) are integral in the molecular response to thermal stress, and their expression is heritable, but the roles of other candidate families such as the heat-shock factors (HSFs) have not been determined in reptiles. Here, we subject embryonic sea turtles (Caretta caretta) to a biologically realistic thermal stress and employ de novo transcriptomic profiling of brain tissue to investigate the underlying molecular response. From a reference transcriptome of 302 293 transcripts, 179 were identified as differentially expressed between treatments. As anticipated, genes enriched in the heat-shock treatment were primarily associated with the Hsp families, or were genes whose products play similar protein editing and chaperone functions (e.g. bag3, MYOC and serpinh1). Unexpectedly, genes encoding the HSFs were not significantly upregulated under thermal stress, indicating their presence in unstressed cells in an inactive state. Genes that were downregulated under thermal stress were less well functionally defined but were associated with stress response, development and cellular organization, suggesting that developmental processes may be compromised at realistically high temperatures. These results confirm that genes from the Hsp families play vital roles in the thermal tolerance of developing reptile embryos and, in addition with a number of other genes, should be targets for evaluating the capacity of oviparous reptiles to respond adaptively to the effects of climate change.

The effects of elevated temperature and ocean acidification on the metabolic pathways of notothenioid fish

The adaptations used by notothenioid fish to combat extreme cold may have left these fish poorly poised to deal with a changing environment. As such, the expected environmental perturbations brought on by global climate change have the potential to significantly affect the energetic demands and subsequent cellular processes necessary for survival. Despite recent lines of evidence demonstrating that notothenioid fish retain the ability to acclimate to elevated temperatures, the underlying mechanisms responsible for temperature acclimation in these fish remain largely unknown. Furthermore, little information exists on the capacity of Antarctic fish to respond to changes in multiple environmental variables. We have examined the effects of increased temperature and pCO₂ on the rate of oxygen consumption in three notothenioid species, Trematomus bernacchii, Pagothenia borchgrevinki, and Trematomus newnesi. We combined these measurements with analysis of changes in aerobic and anaerobic capacity, lipid reserves, fish condition, and growth rates to gain insight into the metabolic cost associated with acclimation to this dual stress. Our findings indicated that temperature is the major driver of the metabolic responses observed in these fish and that increased pCO₂ plays a small, contributing role to the energetic costs of the acclimation response. All three species displayed varying levels of energetic compensation in response to the combination of elevated temperature and pCO₂. While P. borchgrevinki showed nearly complete compensation of whole animal oxygen consumption rates and aerobic capacity, T. newnesi and T. bernacchii displayed only partial compensation in these metrics, suggesting that at least some notothenioids may require physiological trade-offs to fully offset the energetic costs of long-term acclimation to climate change related stressors.

Comparative transcriptome analysis reveals candidate genes for the biosynthesis of natural insecticide in Tanacetum cinerariifolium

Background

Pyrethrins are monoterpenoids and consist of either a chrysanthemic acid or pyrethric acid with a rethrolone moiety. Natural pyrethrins are safe and eco-friendly while possessing strong insecticidal properties. Despite such advantages of commercial value coming with the eco-friendly tag, most enzymes/genes involved in the pyrethrin biosynthesis pathway remain unidentified and uncharacterized. Since the flowers of Tanacetum cinerariifolium are rich in major pyrethrins, next generation transcriptome sequencing was undertaken to compare the flowers and the leaves of the plant de novo to identify differentially expressed transcripts and ascertain which among them might be involved in and responsible for the differential accumulation of pyrethrins in T. cinerariifolium flowers.

Results

In this first tissue specific transcriptome analysis of the non-model plant T. cinerariifolium, a total of 23,200,000 and 28,500,110 high quality Illumina next generation sequence reads, with a length of 101 bp, were generated for the flower and leaf tissue respectively. After functional enrichment analysis and GO based annotation using public protein databases such as UniRef, PFAM, SMART, KEGG and NR, 4443 and 8901 unigenes were identified in the flower and leaf tissue respectively. These could be assigned to 13344 KEGG pathways and the pyrethrin biosynthesis contextualized. The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was involved in the biosynthesis of acid moiety of pyrethrin and this pathway predominated in the flowers as compared to the leaves. However, enzymes related to oxylipin biosynthesis were found predominantly in the leaf tissue, which suggested that major steps of pyrethrin biosynthesis occurred in the flowers.

Conclusions

Transcriptome comparison between the flower and leaf tissue of T. cinerariifolium provided an elaborate list of tissue specific transcripts that was useful in elucidating the differences in the expression of the biosynthetic pathways leading to differential presence of pyrethrin in the flowers. The information generated on genes, pathways and markers related to pyrethrin biosynthesis in this study will be helpful in enhancing the production of these useful compounds for value added breeding programs. Related proteome comparison to overlay our transcriptome comparison can generate more relevant information to better understand flower specific accumulation of secondary metabolites in general and pyrethrin accumulation in particular.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3409-4) contains supplementary material, which is available to authorized users.

Immune response of the Antarctic teleost Trematomus bernacchii to immunization with Psychrobacter sp. (TAD1)

Adult Trematomus bernacchii have been immunized intraperitoneally with heat-killed cells of the Antarctic marine bacterium Psychrobacter sp. (TAD1) up to 60 days. After immunizations and sampling at various times, fish sera were tested for specific IgM by ELISA, and different tissues (head kidney and spleen) were investigated for transcription of master genes of the acquired immune response (IgM, IgT, TRβ, TRγ). Results from ELISA assays showed a time-dependent induction of specific serum anti-TAD1 IgM, and western blot analysis of TAD1 lysates probed with fish sera revealed enhanced immunoreactivity in immunized animals compared to controls. Quantitative PCR analysis of transcripts coding for IgM, IgT, TRβ, TRγ was performed in T. bernacchii tissues to assess basal expression, and then on cDNAs of cells from head kidney and spleen of fish injected for 8, 24, and 72 h with inactivated TAD1. The results showed a differential basal expression of transcripts in the examined tissues, and a time-dependent strong up-regulation of IgT, TRβ, TRγ genes upon in vivo stimulation with TAD1. These results represent a first in vivo study on the mounting of a specific immune response in an Antarctic teleost species.

Life on the Edge-the Biology of Organisms Inhabiting Extreme Environments: An Introduction to the Symposium

Life persists, even under extremely harsh conditions. While the existence of extremophiles is well known, the mechanisms by which these organisms evolve, perform basic metabolic functions, reproduce, and survive under extreme physical stress are often entirely unknown. Recent technological advances in terms of both sampling and studying extremophiles have yielded new insight into their evolution, physiology and behavior, from microbes and viruses to plants to eukaryotes. The goal of the "Life on the Edge-the Biology of Organisms Inhabiting Extreme Environments" symposium was to unite researchers from taxonomically and methodologically diverse backgrounds to highlight new advances in extremophile biology. Common themes and new insight that emerged from the symposium included the important role of symbiotic associations, the continued challenges associated with sampling and studying extremophiles and the important role these organisms play in terms of studying climate change. As we continue to explore our planet, especially in difficult to reach areas from the poles to the deep sea, we expect to continue to discover new and extreme circumstances under which life can persist.

Integrated studies of organismal plasticity through physiological and transcriptomic approaches: examples from marine polar regions

Transcriptomic methods are now widely used in functional genomic research. The vast amount of information received from these studies comes along with the challenge of developing a precise picture of the functional consequences and the characteristic regulatory mechanisms. Here we assess recent studies in marine species and their adaptation to polar (and seasonal) cold and explore how they have been able to draw reliable conclusions from transcriptomic patterns on functional consequences in the organisms. Our analysis indicates that the interpretation of transcriptomic data suffers from insufficient understanding of the consequences for whole organism performance and fitness and comes with the risk of supporting only preliminary and superficial statements.We propose that the functional understanding of transcriptomic data may be improved by their tighter integration into overarching physiological concepts that support the more specific interpretation of the 'omics' data and, at the same time, can be developed further through embedding the transcriptomic phenomena observed. Such possibilities have not been fully exploited.In the context of thermal adaptation and limitation, we explore preliminary evidence that the concept of oxygen and capacity limited thermal tolerance (OCLTT) may provide sufficient complexity to guide the integration of such data and the development of associated functional hypotheses. At the same time, we identify a lack of methodological approaches linking genes and function to higher levels of integration, in terms of organism and ecosystem functioning, at temporal and geographical scales, to support more reliable conclusions and be predictive with respect to the effects of global changes.

Antarctic notothenioid fish: what are the future consequences of 'losses' and 'gains' acquired during long-term evolution at cold and stable temperatures?

Antarctic notothenioids dominate the fish fauna of the Southern Ocean. Evolution for millions of years at cold and stable temperatures has led to the acquisition of numerous biochemical traits that allow these fishes to thrive in sub-zero waters. The gain of antifreeze glycoproteins has afforded notothenioids the ability to avert freezing and survive at temperatures often hovering near the freezing point of seawater. Additionally, possession of cold-adapted proteins and membranes permits them to sustain appropriate metabolic rates at exceptionally low body temperatures. The notothenioid genome is also distinguished by the disappearance of traits in some species, losses that might prove costly in a warmer environment. Perhaps the best-illustrated example is the lack of expression of hemoglobin in white-blooded icefishes from the family Channichthyidae. Loss of key elements of the cellular stress response, notably the heat shock response, has also been observed. Along with their attainment of cold tolerance, notothenioids have developed an extreme stenothermy and many species perish at temperatures only a few degrees above their habitat temperatures. Thus, in light of today's rapidly changing climate, it is critical to evaluate how these extreme stenotherms will respond to rising ocean temperatures. It is conceivable that the remarkable cold specialization of notothenioids may ultimately leave them vulnerable to future thermal increases and threaten their fitness and survival. Within this context, our review provides a current summary of the biochemical losses and gains that are known for notothenioids and examines these cold-adapted traits with a focus on processes underlying thermal tolerance and acclimation capacity.

Transcriptomic responses to environmental temperature in eurythermal and stenothermal fishes

Ectothermic species like fishes differ greatly in the thermal ranges they tolerate; some eurythermal species may encounter temperature ranges in excess of 25°C, whereas stenothermal species in polar and tropical waters live at essentially constant temperatures. Thermal specialization comes with fitness trade-offs and as temperature increases due to global warming, the physiological basis of specialization and thermal plasticity has become of great interest. Over the past 50 years, comparative physiologists have studied the physiological and molecular differences between stenothermal and eurythermal fishes. It is now well known that many stenothermal fishes have lost an inducible heat shock response (HSR). Recent advances in transcriptomics have now made it possible to examine genome-wide changes in gene expression (GE) in non-model ecologically important fish, broadening our view beyond the HSR to regulation of genes involved in hundreds of other cellular processes. Here, we review the major findings from transcriptomic studies of extreme eurythermal and stenothermal fishes in response to acute and long-term exposure to temperature, both time scales being critically important for predicting climate change responses. We consider possible molecular adaptations that underlie eurythermy and stenothermy in teleosts. Furthermore, we highlight the challenges that still face the field of comparative environmental genomics and suggest fruitful paths of future investigation.

Comparison of De Novo Transcriptome Assemblers and k-mer Strategies Using the Killifish, Fundulus heteroclitus

Background

De novo assembly of non-model organism’s transcriptomes has recently been on the rise in concert with the number of de novo transcriptome assembly software programs. There is a knowledge gap as to what assembler software or k-mer strategy is best for construction of an optimal de novo assembly. Additionally, there is a lack of consensus on which evaluation metrics should be used to assess the quality of de novo transcriptome assemblies.

Result

Six different assembly strategies were evaluated from four different assemblers. The Trinity assembly was used in its default 25 single k-mer value while Bridger, Oases, and SOAPdenovo-Trans were performed with multiple k-mer strategies. Bridger, Oases, and SOAPdenovo-Trans used a small multiple k-mer (SMK) strategy consisting of the k-mer lengths of 21, 25, 27, 29, 31, and 33. Additionally, Oases and SOAPdenovo-Trans were performed using a large multiple k-mer (LMK) strategy consisting of k-mer lengths of 25, 35, 45, 55, 65, 75, and 85. Eleven metrics were used to evaluate each assembly strategy including three genome related evaluation metrics (contig number, N50 length, Contigs >1 kb, reads) and eight transcriptome evaluation metrics (mapped back to transcripts (RMBT), number of full length transcripts, number of open reading frames, Detonate RSEM-EVAL score, and percent alignment to the southern platyfish, Amazon molly, BUSCO and CEGMA databases). The assembly strategy that performed the best, that is it was within the top three of each evaluation metric, was the Bridger assembly (10 of 11) followed by the Oases SMK assembly (8 of 11), the Oases LMK assembly (6 of 11), the Trinity assembly (4 of 11), the SOAP LMK assembly (4 of 11), and the SOAP SMK assembly (3 of 11).

Conclusion

This study provides an in-depth multi k-mer strategy investigation concluding that the assembler itself had a greater impact than k-mer size regardless of the strategy employed. Additionally, the comprehensive performance transcriptome evaluation metrics utilized in this study identified the need for choosing metrics centered on user defined research goals. Based on the evaluation metrics performed, the Bridger assembly was able to construct the best assembly of the testis transcriptome in Fundulus heteroclitus.

Integrated analysis of mRNA-seq and miRNA-seq in the liver of Pelteobagrus vachelli in response to hypoxia

Pelteobagrus vachelli is a well-known commercial species in Asia. However, a sudden lack of oxygen will result in mortality and eventually to pond turnover. Studying the molecular mechanisms of hypoxia adaptation in fishes will not only help us to understand fish speciation and the evolution of the hypoxia-signaling pathway, but will also guide us in the breeding of hypoxia-tolerant fish strains. Despite this, the genetic regulatory network for miRNA-mRNA and the signaling pathways involved in hypoxia responses in fish have remained unexamined. In the present study, we used next-generation sequencing technology to characterise mRNA-seq and miRNA-seq of control- and hypoxia-treated P. vachelli livers to elucidate the molecular mechanisms of hypoxia adaptation. We were able to find miRNA-mRNA pairs using bioinformatics analysis and miRNA prediction algorithms. Furthermore, we compared several key pathways which were identified as involved in the hypoxia response of P. vachelli. Our study is the first report on integrated analysis of mRNA-seq and miRNA-seq in fishes and offers a deeper insight into the molecular mechanisms of hypoxia adaptation. qRT-PCR analysis further confirmed the results of mRNA-Seq and miRNA-Seq analysis. We provide a good case study for analyzing mRNA/miRNA expression and profiling a non-model fish species using next-generation sequencing technology.

RNA-seq reveals a diminished acclimation response to the combined effects of ocean acidification and elevated seawater temperature in Pagothenia borchgrevinki

PURPOSE

The IPCC has reasserted the strong influence of anthropogenic CO2 contributions on global climate change and highlighted the polar-regions as highly vulnerable. With these predictions the cold adapted fauna endemic to the Southern Ocean, which is dominated by fishes of the sub-order Notothenioidei, will face considerable challenges in the near future. Recent physiological studies have demonstrated that the synergistic stressors of elevated temperature and ocean acidification have a considerable, although variable, impact on notothenioid fishes. The present study explored the transcriptomic response of Pagothenia borchgrevinki to increased temperatures and pCO2 after 7, 28 and 56days of acclimation. We compared this response to short term studies assessing heat stress alone and foretell the potential impacts of these stressors on P. borchgrevinki's ability to survive a changing Southern Ocean.

RESULTS

P. borchgrevinki did demonstrate a coordinated stress response to the dual-stressor condition, and even indicated that some level of inducible heat shock response may be conserved in this notothenioid species. However, the stress response of P. borchgrevinki was considerably less robust than that observed previously in the closely related notothenioid, Trematomus bernacchii, and varied considerably when compared across different acclimation time-points. Furthermore, the molecular response of these fish under multiple stressors displayed distinct differences compared to their response to short term heat stress alone.

CONCLUSIONS

When exposed to increased sea surface temperatures, combined with ocean acidification, P. borchgrevinki demonstrated a coordinated stress response that has already peaked by 7days of acclimation and quickly diminished over time. However, this response is less dramatic than other closely related notothenioids under identical conditions, supporting previous research suggesting that this notothenioid species is less sensitive to environmental variation.

Transcriptome wide analyses reveal a sustained cellular stress response in the gill tissue of Trematomus bernacchii after acclimation to multiple stressors

BACKGROUND

As global climate change progresses, the Southern Ocean surrounding Antarctica is poised to undergo potentially rapid and substantial changes in temperature and pCO2. To survive in this challenging environment, the highly cold adapted endemic fauna of these waters must demonstrate sufficient plasticity to accommodate these changing conditions or face inexorable decline. Previous studies of notothenioids have focused upon the short-term response to heat stress; and more recently the longer-term physiological response to the combined stress of increasing temperatures and pCO2. This inquiry explores the transcriptomic response of Trematomus bernacchii to increased temperatures and pCO2 at 7, 28 and 56 days, in an attempt to discern the innate plasticity of T. bernacchii available to cope with a changing Southern Ocean.

RESULTS

Differential gene expression analysis supported previous research in that T. bernacchii exhibits no inducible heat shock response to stress conditions. However, T. bernacchii did demonstrate a strong stress response to the multi-stressor condition in the form of metabolic shifts, DNA damage repair, immune system processes, and activation of apoptotic pathways combined with negative regulation of cell proliferation. This response declined in magnitude over time, but aspects of this response remained detectable throughout the acclimation period.

CONCLUSIONS

When exposed to the multi-stressor condition, T. bernacchii demonstrates a cellular stress response that persists for a minimum of 7 days before returning to near basal levels of expression at longer acclimation times. However, subtle changes in expression persist in fish acclimated for 56 days that may significantly affect the fitness T. bernacchii over time.

Transcriptional response to heat shock in liver of snow trout (Schizothorax richardsonii)--a vulnerable Himalayan Cyprinid fish

The snow trout (Schizothorax richardsonii) belonging to family Cyprinidae, is an endemic fish of the Himalayan region. The species is tagged as vulnerable species in the IUCN red list of threatened species. The fish thrives well in snowmelt water of several streams and rivers in the region but are occasionally exposed to more than 20 °C during the summer season. Therefore, we have used deep RNA sequencing to decipher the transcriptome of snow trout and characterize the genes and molecular pathways involved in heat shock response. In this study 72,601,298 and 65,428,283 raw reads for heat-shocked and control, respectively, were obtained by Illumina paired-end sequencing technology. The de novo assembled transcriptome was tested for differential gene expression across the treatment groups. The quality of assembly was high with N75 and N50 lengths of 461 and 1274 bases, respectively. A total of 65 unique transcripts were differentially expressed in liver under heat shock and control. Annotated blast matches reveal that differentially expressed transcripts correspond to critical chaperones and molecular pathways, previously shown to be important for thermal stress in other fish species. Eight randomly selected heat-stressed responsive transcripts were also observed to be upregulated during qRT-PCR analysis. This study is the preliminary step to understanding the responses during sudden environmental changes like heat shock. The reference transcriptome database would also aid further studies on biological and physiological aspects of the snow trout under abiotic stresses.

Non-Antarctic notothenioids: Past phylogenetic history and contemporary phylogeographic implications in the face of environmental changes

The non-Antarctic Notothenioidei families, Bovichtidae, Pseudaphritidae and Eleginopsidae, diverged early from the main notothenioid lineage. They are important in clarifying the early evolutionary processes that triggered notothenioid evolution in the Antarctic. The early-diverged group represents 8% of all notothenioid species and never established themselves on the Antarctic shelf. Most attention has been paid to the Antarctic notothenioids and their limited physiological tolerance to climate change and increased temperatures. In this review, we discuss key life history traits that are characteristic of the non-Antarctic early-diverged notothenioid taxa as well as the genetic resources and population differentiation information available for this group. We emphasise the population fitness and dynamics of these species and indicate how resource management and conservation of the group can be strengthened through an integrative approach. Both Antarctic waters and the non-Antarctic regions face rapid temperature rises combined with strong anthropogenic exploitation. While it is expected that early-diverged notothenioid species may have physiological advantages over high Antarctic species, it is difficult to predict how climate changes might alter the geographic range, behaviour, phenology and ultimately genetic variability of these species. It is possible, however, that their high degree of endemism and dependence on local environmental specificities to complete their life cycles might enhance their vulnerability.

De novo assembly and characterization of tissue-specific transcriptome in the endangered golden mahseer, Tor putitora

The golden mahseer (Tor putitora) graces most of the Himalayan Rivers of India and neighboring South Asian countries. Despite its several importance as a research model, as food, and in sport fishing, knowledge on transcriptome database is nil. Therefore, it was targeted to develop reference transcriptome databases of the species using next-generation sequencing. In the present study, 100,540,130 high-quality paired-end reads were obtained from six cDNA libraries of spleen, liver, gill, kidney, muscle, and brain with 28.4 GB data using Illumina paired-end sequencing technology. Tissue-specific transcriptomes as well as complete transcriptome assembly were analyzed for concise representation of the study. In brief, the de novo assembly of individual tissue resulted in an average of 31,829 (18,512–46,348) contigs per sample, while combined transcriptome comprised 77,907 unique transcript fragments (unigenes) assembled from reads of six tissues. Approximately 75,407 (96.8%) unigenes could be annotated according to their homology matches in the nr, SwisseProt, GO, or KEGG databases. Comparative analysis showed that 84% of the unigenes have significant similarity to zebra fish RefSeq proteins. Tissue-specific-dominated genes were also identified to hypothesize their localization and expression in individual tissue. In addition, 2485 simple sequence repeats (SSRs) were detected from 77,907 transcripts in the combined transcriptome of the golden mahseer. This study has generated organ-specific transcriptome profiles, which will be helpful to understand the local adaptation, genome evolution, and also future functional studies on immune system of the golden mahseer.

•

Organ specific and concatenated transcriptome of 6 organs of golden mahseer was generated using illumina sequencing.

•

A total of 100,540,130 high-quality paired-ends reads with 28.4 GB data that assembled into 77,907 contigs were obtained.

•

Contigs were annotated using GO and KEGG pathway and organ-specific few genes were identified.

•

Thirty four usable SSR markers were extracted from 2485 identified.

Peroxiredoxin 6 from the Antarctic emerald rockcod: molecular characterization of its response to warming

In the present study, we describe the purification and molecular characterization of two peroxiredoxins (Prdxs), referred to as Prdx6A and Prdx6B, from Trematomus bernacchii, a teleost widely distributed in many areas of Antarctica, that plays a pivotal role in the Antarctic food chain. The two putative amino acid sequences were compared with Prdx6 orthologs from other fish, highlighting a high percentage of identity and similarity with the respective variant, in particular for the residues that are essential for the characteristic peroxidase and phospholipase activities of these enzymes. Phylogenetic analyses suggest the appearance of the two prdx6 genes through a duplication event before the speciation that led to the differentiation of fish families and that the evolution of the two gene variants seems to proceed together with the evolution of fish orders and families. The temporal expression of Prdx6 mRNA in response to short-term thermal stress showed a general upregulation of prdx6b and inhibition of prdx6a, suggesting that the latter is the variant most affected by temperature increase. The variations of mRNA accumulation are more conspicuous in heart than the liver, probably related to behavioral changes of the specimens in response to elevated temperature. These data, together with the peculiar differences between the molecular structures of the two Prdx6s in T. bernacchii as well as in the tropical species Stegastes partitus, suggest an adaptation that allowed these poikilothermic aquatic vertebrates to colonize very different environments, characterized by different temperature ranges.

Elevated temperature causes metabolic trade-offs at the whole-organism level in the Antarctic fish Trematomus bernacchii

As a response to ocean warming, shifts in fish species distribution and changes in production have been reported that have been partly attributed to temperature effects on the physiology of animals. The Southern Ocean hosts some of the most rapidly warming regions on earth and Antarctic organisms are reported to be especially temperature sensitive. While cellular and molecular organismic levels appear, at least partially, to compensate for elevated temperatures, the consequences of acclimation to elevated temperature for the whole organism are often less clear. Growth and reproduction are the driving factors for population structure and abundance. The aim of this study was to assess the effect of long-term acclimation to elevated temperature on energy budget parameters in the high-Antarctic fish Trematomus bernacchii. Our results show a complete temperature compensation for routine metabolic costs after 9 weeks of acclimation to 4°C. However, an up to 84% reduction in mass growth was measured at 2 and 4°C compared with the control group at 0°C, which is best explained by reduced food assimilation rates at warmer temperatures. With regard to a predicted temperature increase of up to 1.4°C in the Ross Sea by 2200, such a significant reduction in growth is likely to affect population structures in nature, for example by delaying sexual maturity and reducing production, with severe impacts on Antarctic fish communities and ecosystems.

Neural transcriptome reveals molecular mechanisms for temporal control of vocalization across multiple timescales

Background

Vocalization is a prominent social behavior among vertebrates, including in the midshipman fish, an established model for elucidating the neural basis of acoustic communication. Courtship vocalizations produced by territorial males are essential for reproductive success, vary over daily and seasonal cycles, and last up to hours per call. Vocalizations rely upon extreme synchrony and millisecond precision in the firing of a homogeneous population of motoneurons, the vocal motor nucleus (VMN). Although studies have identified neural mechanisms driving rapid, precise, and stable neuronal firing over long periods of calling, little is known about underlying genetic/molecular mechanisms.

Results

We used RNA sequencing-based transcriptome analyses to compare patterns of gene expression in VMN to the surrounding hindbrain across three daily and seasonal time points of high and low sound production to identify candidate genes that underlie VMN’s intrinsic and network neuronal properties. Results from gene ontology enrichment, enzyme pathway mapping, and gene category-wide expression levels highlighted the importance of cellular respiration in VMN function, consistent with the high energetic demands of sustained vocal behavior. Functionally important candidate genes upregulated in the VMN, including at time points corresponding to high natural vocal activity, encode ion channels and neurotransmitter receptors, hormone receptors and biosynthetic enzymes, neuromodulators, aerobic respiration enzymes, and antioxidants. Quantitative PCR and RNA-seq expression levels for 28 genes were significantly correlated. Many candidate gene products regulate mechanisms of neuronal excitability, including those previously identified in VMN motoneurons, as well as novel ones that remain to be investigated. Supporting evidence from previous studies in midshipman strongly validate the value of transcriptomic analyses for linking genes to neural characters that drive behavior.

Conclusions

Transcriptome analyses highlighted a suite of molecular mechanisms that regulate vocalization over behaviorally relevant timescales, spanning milliseconds to hours and seasons. To our knowledge, this is the first comprehensive characterization of gene expression in a dedicated vocal motor nucleus. Candidate genes identified here may belong to a conserved genetic toolkit for vocal motoneurons facing similar energetic and neurophysiological demands.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1577-2) contains supplementary material, which is available to authorized users.

The environmental genomics of metazoan thermal adaptation

Continued and accelerating change in the thermal environment places an ever-greater priority on understanding how organisms are going to respond. The paradigm of ‘move, adapt or die', regarding ways in which organisms can respond to environmental stressors, stimulates intense efforts to predict the future of biodiversity. Assuming that extinction is an unpalatable outcome, researchers have focussed attention on how organisms can shift in their distribution to stay in the same thermal conditions or can stay in the same place by adapting to a changing thermal environment. How likely these respective outcomes might be depends on the answer to a fundamental evolutionary question, namely what genetic changes underpin adaptation to the thermal environment. The increasing access to and decreasing costs of next-generation sequencing (NGS) technologies, which can be applied to both model and non-model systems, provide a much-needed tool for understanding thermal adaptation. Here we consider broadly what is already known from non-NGS studies about thermal adaptation, then discuss the benefits and challenges of different NGS methodologies to add to this knowledge base. We then review published NGS genomics and transcriptomics studies of thermal adaptation to heat stress in metazoans and compare these results with previous non-NGS patterns. We conclude by summarising emerging patterns of genetic response and discussing future directions using these increasingly common techniques.

Analysis and characterization of the head kidney transcriptome from the Antarctic fish Trematomus bernacchii (Teleostea, Notothenioidea): a source for immune relevant genes

In this study we describe the de novo assembled head kidney transcriptome of the Antarctic notothenioid fish Trematomus bernacchii, an important model species for biochemical, environmental and immunological studies. RNA-seq data was generated using Illumina paired-end sequencing, obtaining ~7 Gbp of sequence data, which were assembled into 96,641 contigs and annotated with the Trinotate pipeline. Since this sequence collection is expected to contain a relevant number of immunity-related transcripts, it will be used as a reference for future immunological studies in this species.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

Is warmer better? Decreased oxidative damage in notothenioid fish after long-term acclimation to multiple stressors

Antarctic fish of the suborder Notothenioidei have evolved several unique adaptations to deal with subzero temperatures. However, these adaptations may come with physiological trade-offs, such as an increased susceptibility to oxidative damage. As such, the expected environmental perturbations brought on by global climate change have the potential to significantly increase the level of oxidative stress and cellular damage in these endemic fish. Previous single stressor studies of the notothenioids have shown they possess the capacity to acclimate to increased temperatures, but the cellular-level effects remain largely unknown. Additionally, there is little information on the ability of Antarctic fish to respond to ecologically relevant environmental changes where multiple variables change concomitantly. We have examined the potential synergistic effects that increased temperature and Ṗ(CO2) have on the level of protein damage in Trematomus bernacchii, Pagothenia borchgrevinki and Trematomus newnesi, and combined these measurements with changes in total enzymatic activity of catalase (CAT) and superoxide dismutase (SOD) in order to gauge tissue-specific changes in antioxidant capacity. Our findings indicate that total SOD and CAT activity levels displayed only small changes across treatments and tissues. Short-term acclimation to decreased seawater pH and increased temperature resulted in significant increases in oxidative damage. Surprisingly, despite no significant change in antioxidant capacity, cellular damage returned to near-basal levels, and significantly decreased in T. bernacchii, after long-term acclimation. Overall, these data suggest that notothenioid fish currently maintain the antioxidant capacity necessary to offset predicted future ocean conditions, but it remains unclear whether this capacity comes with physiological trade-offs.

RNA-seq analyses of cellular responses to elevated body temperature in the high Antarctic cryopelagic nototheniid fish Pagothenia borchgrevinki

Through evolution in the isolated, freezing (-1.9°C) Southern Ocean, Antarctic notothenioid fish have become cold-adapted as well as cold-specialized. Notothenioid cold specialization is most evident in their limited tolerance to heat challenge, and an apparent loss of the near universal inducible heat shock (HSP70) response. Beyond these it remains unclear how broadly cold specialization pervades the underlying tissue-wide cellular responses. We report the first analysis of massively parallel RNA sequencing (RNA-seq) to identify gene expression changes in the liver in response to elevated body temperature of a high-latitude Antarctic nototheniid, the highly cold-adapted and cold-specialized cryopelagic bald notothen, Pagothenia borchgrevinki. From a large (14,873) mapped set of qualified, annotated liver transcripts, we identified hundreds of significantly differentially expressed genes following two and four days of 4°C exposure, suggesting substantial transcriptional reorganization in the liver when body temperature was raised 5°C above native water temperature. Most notably, and in sharp contrast to heat stressed non-polar fish species, was a widespread down-regulation of nearly all classes of molecular chaperones including HSP70, as well as polyubiquitins that are associated with proteosomal degradation of damaged proteins. In parallel, genes involved in the cell cycle were down-regulated by day two of 4°C exposure, signifying slowing cellular proliferation; by day four, genes associated with transcriptional and translational machineries were down-regulated, signifying general slowing of protein biosynthesis. The log2 fold differential transcriptional changes are generally of small magnitudes but significant, and in total portray a broad down turn of cellular activities in response to four days of elevated body temperature in the cold-specialized bald notothen.

Resistance to Botrytis cinerea in Solanum lycopersicoides involves widespread transcriptional reprogramming

BACKGROUND

Tomato (Solanum lycopersicum), one of the world's most important vegetable crops, is highly susceptible to necrotrophic fungal pathogens such as Botrytis cinerea and Alternaria solani. Improving resistance through conventional breeding has been hampered by a shortage of resistant germplasm and difficulties in introgressing resistance into elite germplasm without linkage drag. The goal of this study was to explore natural variation among wild Solanum species to identify new sources of resistance to necrotrophic fungi and dissect mechanisms underlying resistance against B. cinerea.

RESULTS

Among eight wild species evaluated for resistance against B. cinerea and A. solani, S. lycopersicoides expressed the highest levels of resistance against both pathogens. Resistance against B. cinerea manifested as containment of pathogen growth. Through next-generation RNA sequencing and de novo assembly of the S. lycopersicoides transcriptome, changes in gene expression were analyzed during pathogen infection. In response to B. cinerea, differentially expressed transcripts grouped into four categories: genes whose expression rapidly increased then rapidly decreased, genes whose expression rapidly increased and plateaued, genes whose expression continually increased, and genes with decreased expression. Homology-based searches also identified a limited number of highly expressed B. cinerea genes. Almost immediately after infection by B. cinerea, S. lycopersicoides suppressed photosynthesis and metabolic processes involved in growth, energy generation, and response to stimuli, and simultaneously induced various defense-related genes, including pathogenesis-related protein 1 (PR1), a beta-1,3-glucanase (glucanase), and a subtilisin-like protease, indicating a shift in priority towards defense. Moreover, cluster analysis revealed novel, uncharacterized genes that may play roles in defense against necrotrophic fungal pathogens in S. lycopersicoides. The expression of orthologous defense-related genes in S. lycopersicum after infection with B. cinerea revealed differences in the onset and intensity of induction, thus illuminating a potential mechanism explaining the increased susceptibility. Additionally, metabolic pathway analyses identified putative defense-related categories of secondary metabolites.

CONCLUSIONS

In sum, this study provided insight into resistance against necrotrophic fungal pathogens in the Solanaceae, as well as novel sequence resources for S. lycopersicoides.