Model of gene expression in extreme cold - reference transcriptome for the high-Antarctic cryopelagic notothenioid fish Pagothenia borchgrevinki

De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome

BACKGROUND

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited.

RESULTS

A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome.

CONCLUSIONS

We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Warm acclimation alters antioxidant defences but not metabolic capacities in the Antarctic fish, Notothenia coriiceps

The Southern Ocean surrounding the Western Antarctic Peninsula region is rapidly warming. Survival of members of the dominant suborder of Antarctic fishes, the Notothenioidei, will likely require thermal plasticity and adaptive capacity in key traits delimiting thermal tolerance. Herein, we have assessed the thermal plasticity of several cellular and biochemical pathways, many of which are known to be associated with thermal tolerance in notothenioids, including mitochondrial function, activities of aerobic and anaerobic enzymes, antioxidant defences, protein ubiquitination and degradation in cardiac, oxidative skeletal muscles and gill of Notothenia coriiceps warm acclimated to 4°C for 22 days or 5°C for 42 days. Levels of triacylglycerol (TAG) were measured in liver and oxidative and glycolytic skeletal muscles, and glycogen in liver and glycolytic muscle to assess changes in energy stores. Metabolic pathways displayed minimal thermal plasticity, yet antioxidant defences were lower in heart and oxidative skeletal muscles of warm-acclimated animals compared with animals held at ambient temperature. Despite higher metabolic rates at elevated temperature, energy storage depots of TAG and glycogen increase in liver and remain unchanged in muscle with warm acclimation. Overall, our studies reveal that N. coriiceps displays thermal plasticity in some key traits that may contribute to their survival as the Southern Ocean continues to warm.

Resilience of cardiac performance in Antarctic notothenioid fishes in a warming climate

Warming in the region of the Western Antarctic Peninsula is occurring at an unprecedented rate, which may threaten the survival of Antarctic notothenioid fishes. Herein, we review studies characterizing thermal tolerance and cardiac performance in notothenioids - a group that includes both red-blooded species and the white-blooded, haemoglobinless icefishes - as well as the relevant biochemistry associated with cardiac failure during an acute temperature ramp. Because icefishes do not feed in captivity, making long-term acclimation studies unfeasible, we focus only on the responses of red-blooded notothenioids to warm acclimation. With acute warming, hearts of the white-blooded icefish Chaenocephalus aceratus display persistent arrhythmia at a lower temperature (8°C) compared with those of the red-blooded Notothenia coriiceps (14°C). When compared with the icefish, the enhanced cardiac performance of N. coriiceps during warming is associated with greater aerobic capacity, higher ATP levels, less oxidative damage and enhanced membrane integrity. Cardiac performance can be improved in N. coriiceps with warm acclimation to 5°C for 6-9 weeks, accompanied by an increase in the temperature at which cardiac failure occurs. Also, both cardiac mitochondrial and microsomal membranes are remodelled in response to warm acclimation in N. coriiceps, displaying homeoviscous adaptation. Overall, cardiac performance in N. coriiceps is malleable and resilient to warming, yet thermal tolerance and plasticity vary among different species of notothenioid fishes; disruptions to the Antarctic ecosystem driven by climate warming and other anthropogenic activities endanger the survival of notothenioids, warranting greater protection afforded by an expansion of marine protected areas.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

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.

Transcriptome analysis of the growth performance of hybrid mandarin fish after food conversion

During recent years, China has become a hotspot for the domestication of mandarin fish, and this is of great commercial value. Although the food preference of domesticated mandarin fish has been studied, little is known about genes regulating their growth. We raised hybrid mandarin fish on artificial feed for 3 months, the results showed that the survival rate of hybrid mandarin fish was 60.00%. Their total length and body weight were 18.34 ±0.43 cm and 100.44 ±4.87 g. The absolute length and weight gain rates were 0.14 cm/d and 1.08 g/d, respectively. Finally, RNA sequencing (RNA-Seq) was performed to identify potential genes and pathways activated in response to growth performance. The transcriptome analysis generated 68, 197 transcripts and 45,871 unigenes. Among them, 1025 genes were up-regulated and 593 genes were down-regulated between the fast- and slow-growth fish. Finally, we obtained 32 differentially expressed genes, which were mainly related to fatty acid biosynthesis (e.g. FASN and ACACB), collecting duct acid secretion (e.g. ATP6E and KCC4), cell cycle (e.g. CDC20 and CCNB), and the insulin-like growth factor (IGF) system (IGFBP1). These pathways might be related to the growth of hybrid mandarin fish. In addition, more potential single nucleotide polymorphisms (SNPs) were detected in the fast-growth fish than in the slow-growth fish. The results suggest that the interaction of metabolism and abundant alleles might determine the growth of hybrid mandarin fish after food conversion.

Metabolic responses in Scophthalmus maximus kidney subjected to thermal stress

Turbot (Scophthalmus maximus) is an economically important marine fish cultured in China. In this study, fish in the experimental group were exposed to four temperatures: 15, 20, 25 and 28 °C. Metabolomics analysis and quantitative real-time PCR were used to assess changes in metabolic profiling and gene expression associated with thermal stress. The results showed the levels of heat shock protein 70 (HSP70), heat shock protein 90 (HSP90), blood creatinine and cortisol in S. maximus were all significantly upregulated (P < 0.05), indicating a stress response at 25 °C or higher. Challenge with thermal stress significantly increased expression levels of succinate dehydrogenase (SDH), fructose-1, 6-bisphosphatase (FBPase), malate dehydrogenase (MDH), cytosolic phosphoenolpyruvate carboxykinase (cPEPCK), glucose-6-phosphatase (G6Pase) and aspartate aminotransferase (AST) (P < 0.05). However, there was no effect on the expression levels of lactate dehydrogenase (LDH), alanine aminotransferase (ALT) and mitochondrial phosphoenolpyruvate carboxykinase (mPEPCK). Moreover, high temperature decreased levels of glycogenic amino acids, including histidine, threonine, glutamine, phenylalanine, arginine, serine, tyrosine, methionine and isoleucine. These findings suggest a significant correlation between gene expression and regulation of carbohydrate and amino acid metabolism in heat-stressed S. maximus kidney. In addition, the maintenance of aerobic metabolism and activation of gluconeogenesis appeared to be a critical metabolic strategy in combating heat stress in turbot kidney.

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.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

A tale of two genes: divergent evolutionary fate of haptoglobin and hemopexin in hemoglobinless antarctic icefishes

Evolution of Antarctic notothenioid fishes in the isolated freezing Southern Ocean have led to remarkable trait gains and losses. One of the most extraordinary was the loss of the major oxygen carrier hemoglobin (Hb) in the icefishes (family Channichthyidae). While the mechanisms of this loss and the resulting compensatory changes have been well studied, the impact of Hb loss on the network of genes that once supported its recycling and disposal has remained unexplored. Here we report the functional fate and underlying molecular changes of two such key Hb-supporting proteins across the icefish family - haptoglobin (Hp) and hemopexin (Hx), crucial in removing cytotoxic free Hb and heme respectively. Hp plays a critical role in binding free Hb for intracellular recycling and absent its primary client, icefish Hp transcription is now vanishingly little and translation into a functional protein is nearly silenced. Hp genotype degeneration has manifested in separate lineages of the icefish phylogeny with three distinct nonsense mutations and a deletion-frameshift, as well as mutated polyadenylation signal sequences. Thus, Hb loss appears to have diminished selective constraint on Hp maintenance, resulting in its stochastic, co-evolutionary drift towards extinction. Hx binds free heme for iron recycling in hepatocytes. In contrast to Hp, Hx genotype integrity is preserved in the icefishes and transcription occurs at comparable levels to the red-blooded notothenioids. The persistence of Hx likely owes to continued selective pressure for its function from mitochondrial and non-Hb cellular hemoproteins.

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.

De Novo Transcriptome Characterization and Growth-Related Gene Expression Profiling of Diploid and Triploid Bighead Catfish (Clarias macrocephalus Günther, 1864)

To enhance understanding of triploid gene expression, the transcriptome information from bighead catfish (Clarias macrocephalus Günther, 1864) was studied using the paired-end Illumina HiSeq™ 2000 sequencing platform. In total, 68,227,832 raw reads were generated from liver tissues and 53,149 unigenes were assembled, with an average length of 765 bp and N50 length of 1283 bp. Of these unigenes, 33,428 (62.89%) could be annotated according to their homology with matches in the NCBI non-redundant (Nr), NCBI nucleotide (Nt), Swiss-Prot, Clusters of Orthologous Groups (COG), gene ontology (GO), or Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Relative expression of liver genes between diploid and triploid bighead catfish revealed more than 90% of the annotated unigenes similarly expressed, regardless of ploidy, whereas 362 upregulated and 83 downregulated with at least a twofold change in triploid relative to diploid. Quantitative real-time PCR of 15 differentially expressed growth-related genes showed consistency between the expression profiles of those genes with the results from RNA-seq analysis. Our results showed that genes in C. macrocephalus liver responded independently to triploidy with the majority showing similar expression levels between diploid and triploid (a dosage compensation phenomenon). The underlying mechanism of the varying gene expression patterns was discussed. Notably, 5 of the top 20 upregulated genes associated with stress response and thus may reflect stress caused by triploidy. The present study adds a substantial contribution to the sequence data available for C. macrocephalus and hence provides valuable resources for further studies. Furthermore, it gives information that may enhance understanding of triploid physiology.

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.

Transcriptome comparison reveals a genetic network regulating the lower temperature limit in fish

Transcriptional plasticity is a major driver of phenotypic differences between species. The lower temperature limit (LTL), namely the lower end of survival temperature, is an important trait delimiting the geographical distribution of a species, however, the genetic mechanisms are poorly understood. We investigated the inter-species transcriptional diversification in cold responses between zebrafish Danio rerio and tilapia Oreochromis niloticus, which were reared at a common temperature (28 °C) but have distinct LTLs. We identified significant expressional divergence between the two species in the orthologous genes from gills when the temperature cooled to the LTL of tilapia (8 °C). Five KEGG pathways were found sequentially over-represented in the zebrafish/tilapia divergently expressed genes in the duration (12 hour) of 8 °C exposure, forming a signaling cascade from metabolic regulation to apoptosis via FoxO signaling. Consistently, we found differential progression of apoptosis in the gills of the two species in which zebrafish manifested a delayed and milder apoptotic phenotype than tilapia, corresponding with a lower LTL of zebrafish. We identified diverged expression in 25 apoptosis-related transcription factors between the two species which forms an interacting network with diverged factors involving the FoxO signaling and metabolic regulation. We propose a genetic network which regulates LTL in fishes.

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.

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.

Genomic incompatibilities in the diploid and tetraploid offspring of the goldfish × common carp cross

Polyploidy is much rarer in animals than in plants but it is not known why. The outcome of combining two genomes in vertebrates remains unpredictable, especially because polyploidization seldom shows positive effects and more often results in lethal consequences because viable gametes fail to form during meiosis. Fortunately, the goldfish (maternal) × common carp (paternal) hybrids have reproduced successfully up to generation 22, and this hybrid lineage permits an investigation into the genomics of hybridization and tetraploidization. The first two generations of these hybrids are diploids, and subsequent generations are tetraploids. Liver transcriptomes from four generations and their progenitors reveal chimeric genes (>9%) and mutations of orthologous genes. Characterizations of 18 randomly chosen genes from genomic DNA and cDNA confirm the chimera. Some of the chimeric and differentially expressed genes relate to mutagenesis, repair, and cancer-related pathways in 2nF1. Erroneous DNA excision between homologous parental genes may drive the high percentage of chimeric genes, or even more potential mechanisms may result in this phenomenon. Meanwhile, diploid offspring show paternal-biased expression, yet tetraploids show maternal-biased expression. These discoveries reveal that fast and unstable changes are mainly deleterious at the level of transcriptomes although some offspring still survive their genomic abnormalities. In addition, the synthetic effect of genome shock might have resulted in greatly reduced viability of 2nF2 hybrid offspring. The goldfish × common carp hybrids constitute an ideal system for unveiling the consequences of intergenomic interactions in hybrid vertebrate genomes and their fertility.

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.

Transcriptomic Analysis of Metabolic Pathways in Milkfish That Respond to Salinity and Temperature Changes

Milkfish (Chanos chanos), an important marine aquaculture species in southern Taiwan, show considerable euryhalinity but have low tolerance to sudden drops in water temperatures in winter. Here, we used high throughput next-generation sequencing (NGS) to identify molecular and biological processes involved in the responses to environmental changes. Preliminary tests revealed that seawater (SW)-acclimated milkfish tolerated lower temperatures than the fresh water (FW)-acclimated group. Although FW- and SW-acclimated milkfish have different levels of tolerance for hypothermal stress, to date, the molecular physiological basis of this difference has not been elucidated. Here, we performed a next-generation sequence analysis of mRNAs from four groups of milkfish. We obtained 29669 unigenes with an average length of approximately 1936 base pairs. Gene ontology (GO) analysis was performed after gene annotation. A large number of genes for molecular regulation were identified through a transcriptomic comparison in a KEGG analysis. Basal metabolic pathways involved in hypothermal tolerance, such as glycolysis, fatty acid metabolism, amino acid catabolism and oxidative phosphorylation, were analyzed using PathVisio and Cytoscape software. Our results indicate that in response to hypothermal stress, genes for oxidative phosphorylation, e.g., succinate dehydrogenase, were more highly up-regulated in SW than FW fish. Moreover, SW and FW milkfish used different strategies when exposed to hypothermal stress: SW milkfish up-regulated oxidative phosphorylation and catabolism genes to produce more energy budget, whereas FW milkfish down-regulated genes related to basal metabolism to reduce energy loss.

Omics for aquatic ecotoxicology: control of extraneous variability to enhance the analysis of environmental effects

There are multiple sources of biological and technical variation in a typical ecotoxicology study that may not be revealed by traditional endpoints but that become apparent in an omics dataset. As researchers increasingly apply omics technologies to environmental studies, it will be necessary to understand and control the main source(s) of variability to facilitate meaningful interpretation of such data. For instance, can variability in omics studies be addressed by changing the approach to study design and data analysis? Are there statistical methods that can be employed to correctly interpret omics data and make use of unattributed, inherent variability? The present study presents a review of experimental design and statistical considerations applicable to the use of omics methods in systems toxicology studies. In addition to highlighting potential sources that contribute to experimental variability, this review suggests strategies with which to reduce and/or control such variability so as to improve reliability, reproducibility, and ultimately the application of omics data for systems toxicology.

Global identification of the genetic networks and cis-regulatory elements of the cold response in zebrafish

The transcriptional programs of ectothermic teleosts are directly influenced by water temperature. However, the cis- and trans-factors governing cold responses are not well characterized. We profiled transcriptional changes in eight zebrafish tissues exposed to mildly and severely cold temperatures using RNA-Seq. A total of 1943 differentially expressed genes (DEGs) were identified, from which 34 clusters representing distinct tissue and temperature response expression patterns were derived using the k-means fuzzy clustering algorithm. The promoter regions of the clustered DEGs that demonstrated strong co-regulation were analysed for enriched cis-regulatory elements with a motif discovery program, DREME. Seventeen motifs, ten known and seven novel, were identified, which covered 23% of the DEGs. Two motifs predicted to be the binding sites for the transcription factors Bcl6 and Jun, respectively, were chosen for experimental verification, and they demonstrated the expected cold-induced and cold-repressed patterns of gene regulation. Protein interaction modeling of the network components followed by experimental validation suggested that Jun physically interacts with Bcl6 and might be a hub factor that orchestrates the cold response in zebrafish. Thus, the methodology used and the regulatory networks uncovered in this study provide a foundation for exploring the mechanisms of cold adaptation in teleosts.

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.

Optimizing de novo transcriptome assembly and extending genomic resources for striped catfish (Pangasianodon hypophthalmus)

Striped catfish (Pangasianodon hypophthalmus) is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The culture industry is facing a significant challenge however from saltwater intrusion into many low topographical coastal provinces across the Mekong Delta as a result of predicted climate change impacts. Developing genomic resources for this species can facilitate the production of improved culture lines that can withstand raised salinity conditions, and so we have applied high-throughput Ion Torrent sequencing of transcriptome libraries from six target osmoregulatory organs from striped catfish as a genomic resource for use in future selection strategies. We obtained 12,177,770 reads after trimming and processing with an average length of 97bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 66,451 contigs with an average length of 478bp and N50 length of 506bp. A total of 37,969 contigs (57%) possessed significant similarity with proteins in the non-redundant database. Comparative analyses revealed that a significant number of contigs matched sequences reported in other teleost fishes, ranging in similarity from 45.2% with Atlantic cod to 52% with zebrafish. In addition, 28,879 simple sequence repeats (SSRs) and 55,721 single nucleotide polymorphisms (SNPs) were detected in the striped catfish transcriptome. The sequence collection generated in the current study represents the most comprehensive genomic resource for P. hypophthalmus available to date. Our results illustrate the utility of next-generation sequencing as an efficient tool for constructing a large genomic database for marker development in non-model species.

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.

An experimental and bioinformatics protocol for RNA-Seq analyses of photoperiodic diapause in the Asian tiger mosquito, Aedes albopictus

Photoperiodic diapause is an important adaptation that allows individuals to escape harsh seasonal environments via a series of physiological changes, most notably developmental arrest and reduced metabolism. Global gene expression profiling via RNA-Seq can provide important insights into the transcriptional mechanisms of photoperiodic diapause. The Asian tiger mosquito, Aedes albopictus, is an outstanding organism for studying the transcriptional bases of diapause due to its ease of rearing, easily induced diapause, and the genomic resources available. This manuscript presents a general experimental workflow for identifying diapause-induced transcriptional differences in A. albopictus. Rearing techniques, conditions necessary to induce diapause and non-diapause development, methods to estimate percent diapause in a population, and RNA extraction and integrity assessment for mosquitoes are documented. A workflow to process RNA-Seq data from Illumina sequencers culminates in a list of differentially expressed genes. The representative results demonstrate that this protocol can be used to effectively identify genes differentially regulated at the transcriptional level in A. albopictus due to photoperiodic differences. With modest adjustments, this workflow can be readily adapted to study the transcriptional bases of diapause or other important life history traits in other mosquitoes.

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.

RNA-seq profiles from grass carp tissues after reovirus (GCRV) infection based on singular and modular enrichment analyses

Hemorrhagic disease of the grass carp, Ctenopharyngodon idella, is a fatal disease in fingerlings and yearlings caused by a reovirus, GCRV. RNA-seq data from four diseased grass carp tissues (gill, intestine, liver and spleen) were obtained at 2h before and six times after (2h, 24h, 48h, 72h, 96h and 120h) GCRV challenge. A total of 7.25±0.18 million (M) clean reads and 3.53±0.37M unique reads were obtained per RNA-seq analysis. Compared with controls, there were 9060 unique differentially expressed genes (DEGs) in the four tissues at the six time points post-GCRV challenge. Hierarchical clustering analysis of the DEGs showed that the data from the six time points fell into three branches: 2h, 24h/48h, and 72h/96h/120h. Singular (SEA) and modular enrichment analyses of DEGs per RNA-seq dataset were performed based on gene ontology. The results showed that immune responses occurred in all four tissues, indicating that GCRV probably does not target any tissue specifically. Moreover, during the course of disease, disturbances were observed in lipid and carbohydrate metabolism in each of the organs. SEA of DEGs based on the Kyoto Encyclopedia of Genes and Genomes database was also performed, and this indicated that the complement system and cellular immunity played an important role during the course of hemorrhagic disease. The qPCR of pooled samples of duplicate challenge experiment were used to confirm our RNA-seq approach.