Transcriptomic and genomic evolution under constant cold in Antarctic notothenioid fish

The Genomic Substrate for Adaptive Radiation: Copy Number Variation across 12 Tribes of African Cichlid Species

The initial sequencing of five cichlid genomes revealed an accumulation of genetic variation, including extensive copy number variation in cichlid lineages particularly those that have undergone dramatic evolutionary radiation. Gene duplication has the potential to generate substantial molecular substrate for the origin of evolutionary novelty. We use array-based comparative heterologous genomic hybridization to identify copy number variation events (CNVEs) for 168 samples representing 53 cichlid species including the 5 species for which full genome sequence is available. We identify an average of 50-100 CNVEs per individual. For those species represented by multiple samples, we identify 150-200 total CNVEs suggesting a substantial amount of intraspecific variation. For these species, only ∼10% of the detected CNVEs are fixed. Hierarchical clustering of species according to CNVE data recapitulates phylogenetic relationships fairly well at both the tribe and radiation level. Although CNVEs are detected on all linkage groups, they tend to cluster in "hotspots" and are likely to contain and be flanked by transposable elements. Furthermore, we show that CNVEs impact functional categories of genes with potential roles in adaptive phenotypes that could reasonably promote divergence and speciation in the cichlid clade. These data contribute to a more complete understanding of the molecular basis for adaptive natural selection, speciation, and evolutionary radiation.

Early Stage Adaptation of a Mesophilic Green Alga to Antarctica: Systematic Increases in Abundance of Enzymes and LEA Proteins

It is known that adaptive evolution in permanently cold environments drives cold adaptation in enzymes. However, how the relatively high enzyme activities were achieved in cold environments prior to cold adaptation of enzymes is unclear. Here we report that an Antarctic strain of Chlorella vulgaris, called NJ-7, acquired the capability to grow at near 0 °C temperatures and greatly enhanced freezing tolerance after systematic increases in abundance of enzymes/proteins and positive selection of certain genes. Having diverged from the temperate strain UTEX259 of the same species 2.5 (1.1-4.1) to 2.6 (1.0-4.5) Ma, NJ-7 retained the basic mesophilic characteristics and genome structures. Nitrate reductases in the two strains are highly similar in amino acid sequence and optimal temperature, but the NJ-7 one showed significantly higher abundance and activity. Quantitative proteomic analyses indicated that several cryoprotective proteins (LEA), many enzymes involved in carbon metabolism and a large number of other enzymes/proteins, were more abundant in NJ-7 than in UTEX259. Like nitrate reductase, most of these enzymes were not upregulated in response to cold stress. Thus, compensation of low specific activities by increased enzyme abundance appears to be an important strategy for early stage cold adaptation to Antarctica, but such enzymes are mostly not involved in cold acclimation upon transfer from favorable temperatures to near 0 °C temperatures.

Identification and characterization of miRNAs involved in cold acclimation of zebrafish ZF4 cells

MicroRNAs (miRNAs) play vital roles in various biological processes under multiple stress conditions by leading to mRNA cleavage or translational repression. However, the detailed roles of miRNAs in cold acclimation in fish are still unclear. In the present study, high-throughput sequencing was performed to identify miRNAs from 6 small RNA libraries from the zebrafish embryonic fibroblast ZF4 cells under control (28°C, 30 days) and cold-acclimation (18°C, 30 days) conditions. A total of 414 miRNAs, 349 known and 65 novel, were identified. Among those miRNAs, 24 (19 known and 5 novel) were up-regulated, and 23 (9 known and 14 novel) were down-regulated in cold acclimated cells. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses indicated that the target genes of known differentially expressed miRNAs (DE-miRNA) are involved in cold acclimation by regulation of phosphorylation, cell junction, intracellular signal transduction, ECM-receptor interaction and so on. Moreover, both miR-100-3p inhibitor and miR-16b mimics could protect ZF4 cells under cold stress, indicating the involvement of miRNA in cold acclimation. Further study showed that miR-100-3p and miR-16b could regulate inversely the expression of their target gene (atad5a, cyp2ae1, lamp1, rilp, atxn7, tnika, btbd9), and that overexpression of miR-100-3p disturbed the early embryonic development of zebrafish. In summary, the present data show that miRNAs are closely involved in cold acclimation in zebrafish ZF4 cells and provide information for further understanding of the roles of miRNAs in cold acclimation in fish.

Cytogenetic characterization of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902) through analysis of mitotic chromosomes from early larvae

This paper describes the cytogenetic features of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902), a keystone species of the Antarctic coastal marine ecosystem. Conventional cytogenetic analyses and physical mapping of repetitive DNA sequences were performed on metaphase plates obtained through direct chromosome preparation from P. antarctica early larvae. The Antarctic silverfish have a diploid number (2n) = 48, and a karyotype made up of a majority of two-armed chromosomes (karyotype formula36m/sm + 10st + 2a, fundamental number = 94). Major ribosomal gene repeats were detected on three chromosome pairs (20, 21, and 23), in correspondence of dim DAPI stained regions. Long Interspersed Nuclear Elements (LINEs) were abundant and wide spread over all chromosomes. Overall, the cytogenetic data presented herein are consistent with a long independent cytogenetic and evolutionary history for the species. The large number of two-armed chromosomes, indicative of highly-rearranged karyotype, coupled with a diploid number of 48, a presumed primitive character for this fish group, and the spread of the major ribosomal genes on three chromosome pairs, make the Antarctic silverfish distinct from all other notothenioid species.

The genomic basis for colonizing the freezing Southern Ocean revealed by Antarctic toothfish and Patagonian robalo genomes

Background

The Southern Ocean is the coldest ocean on Earth but a hot spot of evolution. The bottom-dwelling Eocene ancestor of Antarctic notothenioid fishes survived polar marine glaciation and underwent adaptive radiation, forming >120 species that fill all water column niches today. Genome-wide changes enabling physiological adaptations and the rapid expansion of the Antarctic notothenioids remain poorly understood.

Results

We sequenced and compared 2 notothenioid genomes—the cold-adapted and neutrally buoyant Antarctic toothfish Dissostichus mawsoni and the basal Patagonian robalo Eleginops maclovinus, representing the temperate ancestor. We detected >200 protein gene families that had expanded and thousands of genes that had evolved faster in the toothfish, with diverse cold-relevant functions including stress response, lipid metabolism, protein homeostasis, and freeze resistance. Besides antifreeze glycoprotein, an eggshell protein had functionally diversified to aid in cellular freezing resistance. Genomic and transcriptomic comparisons revealed proliferation of selcys–transfer RNA genes and broad transcriptional upregulation across anti-oxidative selenoproteins, signifying their prominent role in mitigating oxidative stress in the oxygen-rich Southern Ocean. We found expansion of transposable elements, temporally correlated to Antarctic notothenioid diversification. Additionally, the toothfish exhibited remarkable shifts in genetic programs towards enhanced fat cell differentiation and lipid storage, and promotion of chondrogenesis while inhibiting osteogenesis in bone development, collectively contributing to the achievement of neutral buoyancy and pelagicism.

Conclusions

Our study revealed a comprehensive landscape of evolutionary changes essential for Antarctic notothenioid cold adaptation and ecological expansion. The 2 genomes are valuable resources for further exploration of mechanisms underlying the spectacular notothenioid radiation in the coldest marine environment.

Draft genome assembly and transcriptome data of the icefish Chionodraco myersi reveal the key role of mitochondria for a life without hemoglobin at subzero temperatures

Antarctic fish belonging to Notothenioidei represent an extraordinary example of radiation in the cold. In addition to the absence of hemoglobin, icefish show a number of other striking peculiarities including large-diameter blood vessels, high vascular densities, mitochondria-rich muscle cells, and unusual mitochondrial architecture. In order to investigate the bases of icefish adaptation to the extreme Southern Ocean conditions we sequenced the complete genome of the icefish Chionodraco myersi. Comparative analyses of the icefish genome with those of other teleost species, including two additional white-blooded and five red-blooded notothenioids, provided a new perspective on the evolutionary loss of globin genes. Muscle transcriptome comparative analyses against red-blooded notothenioids as well as temperate fish revealed the peculiar regulation of genes involved in mitochondrial function in icefish. Gene duplication and promoter sequence divergence were identified as genome-wide patterns that likely contributed to the broad transcriptional program underlying the unique features of icefish mitochondria.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

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.

Historical contingency shapes adaptive radiation in Antarctic fishes

Adaptive radiation illustrates links between ecological opportunity, natural selection and the generation of biodiversity. Central to adaptive radiation is the association between a diversifying lineage and the evolution of phenotypic variation that facilitates the use of new environments or resources. However, is not clear whether adaptive evolution or historical contingency is more important for the origin of key phenotypic traits in adaptive radiation. Here we use targeted sequencing of >250,000 loci across 46 species to examine hypotheses concerning the origin and diversification of key traits in the adaptive radiation of Antarctic notothenioid fishes. Contrary to expectations of adaptive evolution, we show that notothenioids experienced a punctuated burst of genomic diversification and evolved key skeletal modifications before the onset of polar conditions in the Southern Ocean. We show that diversifying selection in pathways associated with human skeletal dysplasias facilitates ecologically important variation in buoyancy among Antarctic notothenioid species, and demonstrate the sufficiency of altered trip11, col1a2 and col1a1a function in zebrafish (Danio rerio) to phenocopy skeletal reduction in Antarctic notothenioids. Rather than adaptation being driven by the cooling of the Antarctic, our results highlight the role of historical contingency in shaping the adaptive radiation of notothenioids. Understanding the historical and environmental context for the origin of key traits in adaptive radiations extends beyond reconstructing events that result in evolutionary innovation, as it also provides a context in forecasting the effects of climate change on the stability and evolvability of natural populations.

The impact of acute thermal stress on the metabolome of the black rockfish (Sebastes schlegelii)

Acute change in water temperature causes heavy economic losses in the aquaculture industry. The present study investigated the metabolic and molecular effects of acute thermal stress on black rockfish (Sebastes schlegelii). Gas chromatography time-of-flight mass spectrometry (GC-TOF-MS)-based metabolomics was used to investigate the global metabolic response of black rockfish at a high water temperature (27°C), low water temperature (5°C) and normal water temperature (16°C). Metabolites involved in energy metabolism and basic amino acids were significantly increased upon acute exposure to 27°C (P < 0.05), and no change in metabolite levels occurred in the low water temperature group. However, certain fatty acid levels were elevated after cold stress (P < 0.05), and this effect was not observed in the 27°C group, suggesting that acute high and low temperature exposures caused different physiological responses. Using quantitative real-time PCR, we analyzed the expression of ubiquitin (ub), hypoxia-inducible factor (hif), lactate dehydrogenase (ldh), and acetyl-CoA carboxylase (acac). Higher expression levels of ub, hif, and ldh (P < 0.05) were observed in the high water temperature group, but no changes in these expression levels occurred in the low water temperature group. Our findings provide a potential metabolic profile for black rockfish when exposed to acute temperature stress and provide some insights into host metabolic and molecular responses to thermal stress.

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.

Proteomic and structural differences in lumpfish skin among the dorsal, caudal and ventral regions

Fish skin is a vital organ that serves a multitude of functions including mechanical protection, homeostasis, osmoregulation and protection against diseases. The expression of skin proteins changes under different physiological conditions. However, little is known about differences in protein expression among various body sites in naïve fish. The objectives of this work is to study potential differences in protein and gene expression among dorsal, caudal and ventral regions of lumpfish skin employing 2D gel based proteomics and real-time PCR and to assess structural differences between these regions by using Alcian blue and Periodic acid Schiff stained skin sections. The proteins collagen alfa-1, collagen alfa-2, heat shock cognate 71 kDa, histone H4, parvalbumin, natterin-2, 40S ribosomal protein S12, topoisomerase A and topoisomerase B were differentially expressed among the three regions. mRNA expression of apoa1, hspa8 and hist1h2b showed significant differences between regions. Skin photomicrographs showed differences in epidermal thickness and goblet cell counts. The ventral region showed relatively high protein expression, goblet cell count and epidermal thickness compared to dorsal and caudal regions. Overall, this study provides an important benchmark for comparative analysis of skin proteins and structure between different parts of the lumpfish body.

Screening of the candidate genes related to low-temperature tolerance of Fenneropenaeus chinensis based on high-throughput transcriptome sequencing

In order to screen the candidate genes of Fenneropenaeus chinensis related to low-temperature tolerance, this research takes juvenile prawns of F. chinensis (P₄₀) in low temperature stress group (4°C) and normal temperature group (18°C) as experimental materials. The results showed that a total of 127,939 Unigenes with average length of 1,190 bp were obtained by assembly, of which 46% were annotated in the Nr database. A total of 1,698 differentially expressed genes were screened by differential gene expression analysis, of which 920 genes showed up-regulated expression and 778 genes showed down-regulated expression. Both GO and KEGG enrichment analysis revealed that differentially expressed genes were enriched in spliceosomes, ribosomes, bile secretion, ABC transport pathways, and cellular nitrogen compound synthesis. A further in-depth analysis obtained 8 genes that may be associated with low-temperature traits of F. chinensis. Five of them displayed up-regulated expression, including ATP-binding cassette protein C, acid ceramidase, glutathione transferase, C-type lectin and heat shock protein HSP70. The remaining three genes, γ-butyl betaine hydroxylase, β-hexosaminidase A and long chain fatty acid-CoA ligase displayed down-regulated expression. Eight differentially expressed genes were randomly selected and the real time RT-PCR verification showed that their expression levels were consistent with the sequencing results, demonstrating the accuracy of the sequencing results. The results of this study provide basic data for revealing the molecular mechanisms of F. chinensis in response to low temperature stress and the molecular assisted breeding of F. chinensis in low temperature.

iTRAQ-based quantitative proteomic analysis of Takifugu fasciatus liver in response to low-temperature stress

Low temperatures profoundly influence the physiological and behavioural processes of ectotherms, especially teleosts, which have made them the subjects of strong interest over time. However, the characteristics of fish cold-tolerance at the protein level remain unclear. Therefore, to shed further light on the molecular mechanisms of low temperature adaptation in fish, we conducted quantitative proteomics on the T. fasciatus liver using iTRAQ. Comparing the proteomic profiles of the T. fasciatus liver at 12 °C and 26 °C, a total of 3741 proteins were identified, and 160 were differentially abundant proteins (DAPs). Among the DAPs, the most significant changes were noted in proteins involved in oxidative stress (nine proteins), mitochondrial enzymes (eleven proteins) and signal transduction (thirteen proteins). The KEGG enrichment analysis indicated significant enhancement of D-arginine and D-ornithine metabolism, MAPK signalling, Wnt signalling and Gap junction pathway. Subsequently, three significantly up-regulated proteins (CIRB, HSP90 and GST) and two significantly down-regulated proteins (FLNB and A2ML1) were validated with parallel reaction monitoring (PRM) assays. Furthermore, the changes in abundance of proteins that are involved in oxidative stress, mitochondrial enzymes and signal transduction were validated at the transcriptional level with qPCR. These verification results show that the experimental data of iTRAQ are reliable. Our results not only deepen the understanding of the mechanisms underlying low-temperature tolerance in fish, but they also may contribute to the enhancement of cold tolerance during its breeding process. SIGNIFICANCE OF THE STUDY: The study focused on a comparative quantitative proteomics analysis of the T. fasciatus liver in response to low temperatures using iTRAQ, which has not yet been reported in the literatures. The results showed that the effect of low temperature on T. fasciatus is significant, including a detoxification of metabolic by-products and oxidative stress, an activation of the mitochondrial enzyme to strengthen energy metabolism, and a negative effect on signal transduction, which result in dysfunction or suboptimal performance. These low-temperature-related changes in the liver proteome of T. fasciatus can facilitate the understanding of the low temperature-related response that takes place in similar conditions in the liver and may contribute to the breeding of cold-resistant strains.

Network of microRNA-transcriptional factor-mRNA in cold response of turbot Scophthalmus maximus

The aim of this study is to understand fish cold-tolerant mechanism. We analyzed the transcriptional reactions to the cold condition in turbot Scophthalmus maximus by using RNA-seq and microRNA (miRNA)-seq. Meio-gynogenetic diploid turbots were treated at 0 °C to distinguish the cold-tolerant (CT) and cold-sensitive (CS) groups. The results showed that there were quite different responses at both mRNA and miRNA levels, with more up-regulated mRNAs (1069 vs. 194) and less down-regulated miRNAs (4 vs. 1) in CT versus CS relative to the control group. The network of miRNA-transcription factor-mRNA, regulating turbot different response to cold stress, was constructed, which involved in cell cycle, component of cell membrane, signal transduction, and circadian rhythm pathways. The above information demonstrates mechanisms by which cold tolerance is increased in fish.

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.

Effects of fulvic acid and fulvic ions on Escherichia coli survival in river under repeated freeze-thaw cycles

The effects of fulvic acid (FA) and ions on mesophilic pathogenic bacteria survival under freeze-thaw (FT) stress in natural water and its resistant mechanisms are rarely understood. Therefore, survival patterns of Escherichia coli in river water added with various concentrations of FA or FA-ion under FT stress were studied in this work. Meanwhile, cell surface hydrophobicity (CSH), unit activities of superoxide dismutase (SOD) and catalase (CAT) were determined and Escherichia coli morphologies were observed to explore the bacterial resistant mechanisms against FT stress. The results demonstrated that FT cycles significantly reduced bacterial quantities as sampling time, i.e. freeze-thaw cycle time increased. And the biggest reducing rate was observed after the first FT cycle in every system. Ttd values, time needed to reach detection limit under FT stress decreased under FT stress as FA was added into water, while the changes of ttd values were quite complicated when FA and various ions existed together. Generally, the ttd values of FA-cation systems exceeded that of FA system except FA-Ca²⁺ systems, but it was opposite for FA-anion systems. CSH was heightened after FT cycles and reached peak value at last sampling time in every system. Mechanical constraint from extracellular ice crystals and high CSH induced bacterial aggregation, which protect inner cells of aggregation from extracellular ice crystals. And the unit activities of SOD were significantly higher than those of CAT. Unit activities of SOD and CAT in large part of tested systems increased with sampling time under FT stress, which reduced reactive oxygen species produced from repeated FT cycles. Thus, these could improve the resistance of Escherichia coli to freeze-thaw stress and promote their survival. This work explored the survival pattern and strategy of Escherichia coli in natural water under FT stress.

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.

Transcriptomic responses to low temperature stress in the Nile tilapia, Oreochromis niloticus

The Nile tilapia, Oreochromis niloticus, is a species of high economic value and extensively cultured. The limited stress tolerance of this species to a low temperature usually leads to mass mortality and great loss. Nevertheless, there is limited information on the molecular mechanisms underlying the susceptibility to low temperature in the tilapia. In this study, tilapia was treated at 28 °C to a lethal temperature of 8 °C by a gradual decrement. Transcriptomic response of the immune organ, kidney, in tilapia was characterized using RNA-seq. In total, 2191 genes were annotated for significant expression, which were mainly associated with metabolism and immunity. Pathway analysis showed that immune-related pathways of phagosome and cell adhesion molecules (CAMs) pathway were significantly down-regulated under low temperature. Moreover, ferroptosis, a significantly changed pathway involved in tissue damage and acute renal failure, is reported here for the first time. The levels of serum parameters associated with kidney damage such as urea and uric acid (UA) increased significantly under low temperature. The immunofluorescence staining of the kidney showed that cell apoptosis occurred at low temperature. The results of the present study indicate that exposure to low temperature can cause kidney disfunction and down-regulate the immune-related pathway in the kidney of tilapia. This study provides new insight into the mechanism of kidney damage in fish under low temperature.

Liver transcriptome analysis and de novo annotation of the orange-spotted groupers (Epinephelus coioides) under cold stress

Cold stress has caused great economic loss in fish culture worldwide. Orange-spotted grouper (Epinephelus coioides) is one of the most serious lost aquatic animals in 2016 cold fronts in South China. However, the molecular mechanism of grouper's cold resistance has remained largely unknown. In the present study, HiSeq™2000 (Illumina) was used to analyze the transcriptomic profiles of the liver from grouper under control temperature (CT, 28 °C) and low temperature (LT, 13 °C). Two normalized liver cDNA libraries of CT and LT groups were created. We obtained 51,944,970 and 51,905,036 clean reads from CT and LT groups, respectively. Comparing the LT group to the CT group, a total of 5905 significantly differentially expressed genes (DEGs) were identified, including 2093 up-regulated unigenes and 3812 down-regulated unigenes. GO annotation and functional enrichment analysis indicated that all of the DEGs were classified into three categories: biological process (23 subclasses), cellular component (18 subclasses) and molecular function (13 subclasses). KEGG analysis of the DEGs showed that 2732 DEGs were annotated to 253 signaling pathways. The most highly enriched pathways were cell adhesion molecules, Staphylococcus aureus infection, PPAR signaling pathway, Vibrio cholerae infection, primary immunodeficiency, fatty acid elongation, and we found cold stress mainly affects immunity, metabolic and signal transduction. Thirteen of the DEGs were further validated by qRT-PCR. Our results provide valuable information for further analysis of the mechanisms of groupers response under cold stress.

Proteomic Responses Under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei), one of the most widely cultured shrimp species in the world, often suffers from cold stress. To understand the molecular mechanism of cold tolerance in Pacific white shrimp, we conducted a proteomic analysis on two contrasting shrimp cultivars, namely, cold-tolerant Guihai2 (GH2) and cold-sensitive Guihai1 (GH1), under normal temperature (28°C), under cold stress (16°C), and during recovery to 28°C. In total, 3,349 proteins were identified, among which 2,736 proteins were quantified. Based on gene ontology annotations, differentially expressed proteins largely belonged to biological processes, cellular components, and molecular functions. KEGG pathway annotations indicated that the main changes were observed in the lysosome, ribosomes, and oxidative phosphorylation. Subcellular localization analysis showed a significant increase in proteins present in cytosol, extracellular regions, and mitochondria. Combining enrichment-based clustering analysis and qRT-PCR analysis, we found that glutathione S-transferase, zinc proteinase, m7GpppX diphosphatase, AP2 transcription complex, and zinc-finger transcription factors played a major role in the cold stress response in Pacific white shrimp. Moreover, structure proteins, including different types of lectin and DAPPUDRAFT, were indispensable for cold stress tolerance of the Pacific white shrimp. Results indicate the molecular mechanisms of the Pacific white shrimp in response to cold stress and provide new insight into breeding new cultivars with increased cold tolerance.

Direct-to-consumer DNA testing of 6,000 dogs reveals 98.6-kb duplication associated with blue eyes and heterochromia in Siberian Huskies

Consumer genomics enables genetic discovery on an unprecedented scale by linking very large databases of personal genomic data with phenotype information voluntarily submitted via web-based surveys. These databases are having a transformative effect on human genomics research, yielding insights on increasingly complex traits, behaviors, and disease by including many thousands of individuals in genome-wide association studies (GWAS). The promise of consumer genomic data is not limited to human research, however. Genomic tools for dogs are readily available, with hundreds of causal Mendelian variants already characterized, because selection and breeding have led to dramatic phenotypic diversity underlain by a simple genetic structure. Here, we report the results of the first consumer genomics study ever conducted in a non-human model: a GWAS of blue eyes based on more than 3,000 customer dogs with validation panels including nearly 3,000 more, the largest canine GWAS to date. We discovered a novel association with blue eyes on chromosome 18 (P = 1.3x10-68) and used both sequence coverage and microarray probe intensity data to identify the putative causal variant: a 98.6-kb duplication directly upstream of the Homeobox gene ALX4, which plays an important role in mammalian eye development. This duplication is largely restricted to Siberian Huskies, is strongly associated with the blue-eyed phenotype (chi-square P = 5.2x10-290), and is highly, but not completely, penetrant. These results underscore the power of consumer-data-driven discovery in non-human species, especially dogs, where there is intense owner interest in the personal genomic information of their pets, a high level of engagement with web-based surveys, and an underlying genetic architecture ideal for mapping studies.

Evolution in chronic cold: varied loss of cellular response to heat in Antarctic notothenioid fish

Background

Confined within the freezing Southern Ocean, the Antarctic notothenioids have evolved to become both cold adapted and cold specialized. A marked signature of cold specialization is an apparent loss of the cellular heat shock response (HSR). As the HSR has been examined in very few notothenioid species to-date, it remains unknown whether HSR loss pervades the Antarctic radiation, or whether the broader cellular responses to heat stress has sustained similar loss. Understanding the evolutionary status of these responses in this stenothermal taxon is crucial for evaluating its adaptive potential to ocean warming under climate change.

Results

In this study, we used an acute heat stress protocol followed by RNA-Seq analyses to study the evolution of cellular-wide transcriptional responses to heat stress across three select notothenioid lineages - the basal temperate and nearest non-Antarctic sister species Eleginops maclovinus serving as ancestral proxy, the cryopelagic Pagothenia borchgrevinki and the icefish Chionodraco rastrospinosus representing cold-adapted red-blooded and hemoglobinless Antarctic notothenioids respectively. E. maclovinus displayed robust cellular stress responses including the ER Unfolded Protein Response and the cytosolic HSR, cementing the HSR as a plesiomorphy that preceded Antarctic notothenioid radiation. While the transcriptional response to heat stress was minimal in P. borchgrevinki, C. rastrospinosus exhibited robust responses in the broader cellular networks especially in inflammatory responses despite lacking the classic HSR and UPR.

Conclusion

The disparate patterns observed in these two archetypal Antarctic species indicate the evolutionary status in cellular ability to mitigate acute heat stress varies even among Antarctic lineages, which may affect their adaptive potential in coping with a warming world.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1254-6) contains supplementary material, which is available to authorized users.

Transcriptional fates of human-specific segmental duplications in brain

Despite the importance of duplicate genes for evolutionary adaptation, accurate gene annotation is often incomplete, incorrect, or lacking in regions of segmental duplication. We developed an approach combining long-read sequencing and hybridization capture to yield full-length transcript information and confidently distinguish between nearly identical genes/paralogs. We used biotinylated probes to enrich for full-length cDNA from duplicated regions, which were then amplified, size-fractionated, and sequenced using single-molecule, long-read sequencing technology, permitting us to distinguish between highly identical genes by virtue of multiple paralogous sequence variants. We examined 19 gene families as expressed in developing and adult human brain, selected for their high sequence identity (average >99%) and overlap with human-specific segmental duplications (SDs). We characterized the transcriptional differences between related paralogs to better understand the birth-death process of duplicate genes and particularly how the process leads to gene innovation. In 48% of the cases, we find that the expressed duplicates have changed substantially from their ancestral models due to novel sites of transcription initiation, splicing, and polyadenylation, as well as fusion transcripts that connect duplication-derived exons with neighboring genes. We detect unannotated open reading frames in genes currently annotated as pseudogenes, while relegating other duplicates to nonfunctional status. Our method significantly improves gene annotation, specifically defining full-length transcripts, isoforms, and open reading frames for new genes in highly identical SDs. The approach will be more broadly applicable to genes in structurally complex regions of other genomes where the duplication process creates novel genes important for adaptive traits.

Full-length transcriptome of Misgurnus anguillicaudatus provides insights into evolution of genus Misgurnus

Reconstruction and annotation of transcripts, particularly for a species without reference genome, plays a critical role in gene discovery, investigation of genomic signatures, and genome annotation in the pre-genomic era. This study generated 33,330 full-length transcripts of diploid M. anguillicaudatus using PacBio SMRT Sequencing. A total of 6,918 gene families were identified with two or more isoforms, and 26,683 complete ORFs with an average length of 1,497 bp were detected. Totally, 1,208 high-confidence lncRNAs were identified, and most of these appeared to be precursor transcripts of miRNAs or snoRNAs. Phylogenetic tree of the Misgurnus species was inferred based on the 1,905 single copy orthologous genes. The tetraploid and diploid M. anguillicaudatus grouped into a clade, and M. bipartitus showed a closer relationship with the M. anguillicaudatus. The overall evolutionary rates of tetraploid M. anguillicaudatus were significantly higher than those of other Misgurnus species. Meanwhile, 28 positively selected genes were identified in M. anguillicaudatus clade. These positively selected genes may play critical roles in the adaptation to various habitat environments for M. anguillicaudatus. This study could facilitate further exploration of the genomic signatures of M. anguillicaudatus and provide potential insights into unveiling the evolutionary history of tetraploid loach.

Characterization of calmodulin in the clam Anodonta woodiana: differential expressions in response to environmental Ca2+ and Cd2+

Aims

To explore effect of Ca2+ and Cd2+ on the calmodulin (CaM), one complete cDNA sequence (AwCaM1) was cloned and characterized from the freshwater mussel Anodonta woodiana and its expressions were analyzed.

Materials and methods

The AwCaM1 was cloned from the A. woodiana using the rapid amplification of cDNA ends methods and its expression was determined by real-time PCR.

Results

In the hepatopancreas, AwCaM1 expression was up-regulated with a time and dose dependent pattern in the Ca2+ treated groups (0.01, 0.02, 0.04 and 0.08 mg/L) during experiment observed, and increased more than 56.15% (p<0.05) compared with that of control group. AwCaM1 mRNA level increased more 65.04% (p<0.05) in the Cd2+ treated groups (8 and 16 mg/L). In the gill, AwCaM1 expression increased more than 79.41% (p<0.05) compared with that of control group in all the Ca2+ treated groups, and more than 88.23% (p<0.05) in all the Cd2+ treated groups.

Conclusion

These results indicated that up-regulations of AwCaM1 expression in bivalve A. woodiana are associated with Ca2+ absorb and environmental adaption derived from Ca2+ and Cd2+ treatment.

Characterization of photosynthetic ferredoxin from the Antarctic alga Chlamydomonas sp. UWO241 reveals novel features of cold adaptation

The objective of this work was to characterize photosynthetic ferredoxin from the Antarctic green alga Chlamydomonas sp. UWO241, a key enzyme involved in distributing photosynthetic reducing power. We hypothesize that ferredoxin possesses characteristics typical of cold-adapted enzymes, namely increased structural flexibility and high activity at low temperatures, accompanied by low stability at moderate temperatures. To address this objective, we purified ferredoxin from UWO241 and characterized the temperature dependence of its enzymatic activity and protein conformation. The UWO241 ferredoxin protein, RNA, and DNA sequences were compared with homologous sequences from related organisms. We provide evidence for the duplication of the main ferredoxin gene in the UWO241 nuclear genome and the presence of two highly similar proteins. Ferredoxin from UWO241 has both high activity at low temperatures and high stability at moderate temperatures, representing a novel class of cold-adapted enzymes. Our study reveals novel insights into how photosynthesis functions in the cold. The presence of two distinct ferredoxin proteins in UWO241 could provide an adaptive advantage for survival at cold temperatures. The primary amino acid sequence of ferredoxin is highly conserved among photosynthetic species, and we suggest that subtle differences in sequence can lead to significant changes in activity at low temperatures.

Molecular Evolutionary Analysis of the HCRTR Gene Family in Vertebrates

Hypocretin system is composed of hypocretins (hcrts) and their receptors (hcrtrs), which has multiple vital functions. Hypocretins work via hypocretin receptors and it is reported that functional differentiation occurred in hcrtrs. It is necessary to figure out the evolution process of hypocretin receptors. In our study, we adopt a comprehensive approach and various bioinformatics tools to analyse the evolution process of HCRTR gene family. It turns out that the second round of whole genome duplication in early vertebrate ancestry and the independent round in fish ancestry may contribute to the diversity of HCRTR gene family. HCRTR1 of fishes and mammals are not the same receptor, which means that there are three members in the family. HCRTR2 is proved to be the most ancient one in HCRTR gene family. After duplication events, the structure of HCRTR1 diverged from HCRTR2 owing to relaxed selective pressure. Negative selection is the predominant evolutionary force acting on the HCRTR gene family but HCRTR1 of mammals is found to be subjected to positive selection. Our study gains insight into the molecular evolution process of HCRTR gene family, which contributes to the further study of the system.

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

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

Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells

Long non-coding RNAs (lncRNAs) are increasingly regarded as a key role in regulating diverse biological processes in various tissues and species. Although the cold responsive lncRNAs have been reported in plants, no data is available on screening and functional prediction of lncRNAs in cold acclimation in fish so far. Here we compared the expression profile of lncRNAs in cold acclimated zebrafish embryonic fibroblast cells (ZF4) cultured at 18°C for 30 days with that of cells cultured at 28°C as control by high-throughput sequencing. Totally 8,363 novel lncRNAs were identified. Including known and novel lncRNAs, there are 347 lncRNAs up-regulated and 342 lncRNAs down-regulated in cold acclimated cells. Among the differentially expressed lncRNAs, 74 and 61 were detected only in control cells or cold-acclimated cells, respectively. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses of adjacent genes to the differentially expressed lncRNAs showed that the enriched genes are involved in electron transport, cell adhesion, oxidation-reduction process, and so on. We also predicted the target genes of the differentially expressed lncRNAs by looking for interactions between lncRNAs and mRNAs, and constructed an interaction network. In summary, our genome-wide systematic identification and functional prediction of cold responsive lncRNAs in zebrafish cells suggests a crucial role of lincRNAs in cold acclimation in fish.

Gene duplication and dosage effects during the early emergence of C4 photosynthesis in the grass genus Alloteropsis

The importance of gene duplication for evolutionary diversification has been mainly discussed in terms of genetic redundancy allowing neofunctionalization. In the case of C4 photosynthesis, which evolved via the co-option of multiple enzymes to boost carbon fixation in tropical conditions, the importance of genetic redundancy has not been consistently supported by genomic studies. Here, we test for a different role for gene duplication in the early evolution of C4 photosynthesis, via dosage effects creating rapid step changes in expression levels. Using genome-wide data for accessions of the grass genus Alloteropsis that recently diversified into different photosynthetic types, we estimate gene copy numbers and demonstrate that recurrent duplications in two important families of C4 genes coincided with increases in transcript abundance along the phylogeny, in some cases via a pure dosage effect. While increased gene copy number during the initial emergence of C4 photosynthesis probably offered a rapid route to enhanced expression, we also find losses of duplicates following the acquisition of genes encoding better-suited isoforms. The dosage effect of gene duplication might therefore act as a transient process during the evolution of a C4 biochemistry, rendered obsolete by the fixation of regulatory mutations increasing expression levels.

Characterization of Selenoprotein P cDNA of the Antarctic toothfish Dissostichus mawsoni and its role under cold pressure

Our previous study using comparative genome analysis revealed a significant gene copy number gain of Dissostichus mawsoni selenoprotein P (Dm-SEPP) during the evolutionary radiation of Antarctic notothenioids, suggesting that Dm-SEPP contribute to this process, but the detailed structure and function of this gene product remain unclear. In the present study, the Dm-SEPP cDNA was cloned and characterized. The Dm-SEPP cDNA contains 17 selenocysteines (Sec) encoded by TGA codons and 2 typical SECIS elements located in the 3'-UTR. Evolutionary analysis of the Dm-SEPP gene revealed that it's closely related to the SEPP gene of zebrafish (Danio rerio), showing 51% amino acid similarity. Over-expression of Dm-SEPP could protect mammalian cells under cold pressure, probably via eliminating ROS. Further study showed an increase of endogenous SEPP in zebrafish ZF4 cells under cold pressure, and knockdown of SEPP decreased cell viability, accompanied with increased ROS. Our results suggested a protective role of Dm-SEPP in cold adaptation in Antarctic notothenioids.

The loss of hemoglobin and myoglobin does not minimize oxidative stress in Antarctic icefishes

The unusual pattern of expression of hemoglobin (Hb) and myoglobin (Mb) among Antarctic notothenioid fishes provides an exceptional model system for assessing the impact of these proteins on oxidative stress. We tested the hypothesis that the lack of oxygen-binding proteins may reduce oxidative stress. Levels and activity of pro-oxidants and small-molecule and enzymatic antioxidants, and levels of oxidized lipids and proteins in the liver, oxidative skeletal muscle and heart ventricle were quantified in five species of notothenioid fishes differing in the expression of Hb and Mb. Levels of ubiquitinated proteins and rates of protein degradation by the 20S proteasome were also quantified. Although levels of oxidized proteins and lipids, ubiquitinated proteins, and antioxidants were higher in red-blooded fishes than in Hb-less icefishes in some tissues, this pattern did not persist across all tissues. Expression of Mb was not associated with oxidative damage in the heart ventricle, whereas the activity of citrate synthase and the contents of heme were positively correlated with oxidative damage in most tissues. Despite some tissue differences in levels of protein carbonyls among species, rates of degradation by the 20S proteasome were not markedly different, suggesting either alternative pathways for eliminating oxidized proteins or that redox tone varies among species. Together, our data indicate that the loss of Hb and Mb does not correspond with a clear pattern of either reduced oxidative defense or oxidative damage.

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.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Cold-induced retrotransposition of fish LINEs

Classes of retrotransposons constitute a large portion of metazoan genome. There have been cases reported that genomic abundance of retrotransposons is correlated with the severity of low environmental temperatures. However, the molecular mechanisms underlying such correlation are unknown. We show here by cell transfection assays that retrotransposition (RTP) of a long interspersed nuclear element (LINE) from an Antarctic notothenioid fish Dissostichus mawsoni (dmL1) could be activated by low temperature exposure, causing increased dmL1 copies in the host cell genome. The cold-induced dmL1 propagation was demonstrated to be mediated by the mitogen-activated protein kinases (MAPK)/p38 signaling pathway, which is activated by accumulation of reactive oxygen species (ROS) in cold-stressed conditions. Surprisingly, dmL1 transfected cells showed an increase in the number of viable cells after prolonged cold exposures than non-transfected cells. Features of cold inducibility of dmL1 were recapitulated in LINEs of zebrafish origin both in cultured cell lines and tissues, suggesting existence of a common cold-induced LINE amplification in fishes. The findings reveal an important function of LINEs in temperature adaptation and provid insights into the MAPK/p38 stress responsive pathway that shapes LINE composition in fishes facing cold stresses.

Growth arrest specific gene 2 in tilapia (Oreochromis niloticus): molecular characterization and functional analysis under low-temperature stress

BACKGROUND

Growth arrest specific 2 (gas2) gene is a component of the microfilament system that plays a major role in the cell cycle, regulation of microfilaments, and cell morphology during apoptotic processes. However, little information is available on fish gas2. In this study, the tilapia (Oreochromis niloticus) gas2 gene was cloned and characterized for the first time.

RESULTS

The open reading frame was 1020 bp, encoding 340 amino acids; the 5'-untranslated region (UTR) was 140 bp and the 3'-UTR was 70 bp, with a poly (A) tail. The highest promoter activity occurred in the regulatory region (-3000 to -2400 bp). The Gas2-GFP fusion protein was distributed within the cytoplasm. Quantitative reverse transcription-polymerase chain reaction and western blot analyses revealed that gas2 gene expression levels in the liver, muscle, and brain were clearly affected by low temperature stress. The results of gas2 RNAi showed decreased expression of the gas2 and P53 genes.

CONCLUSION

These results suggest that the tilapia gas2 gene may be involved in low temperature stress-induced apoptosis.

Antarctic fish versus human cytoglobins - The same but yet so different

The cytoglobins of the Antarctic fish Chaenocephalus aceratus and Dissostichus mawsoni have many features in common with human cytoglobin. These cytoglobins are heme proteins in which the ferric and ferrous forms have a characteristic hexacoordination of the heme iron, i.e. axial ligation of two endogenous histidine residues, as confirmed by electron paramagnetic resonance, resonance Raman and optical absorption spectroscopy. The combined spectroscopic analysis revealed only small variations in the heme-pocket structure, in line with the small variations observed for the redox potential. Nevertheless, some striking differences were also discovered. Resonance Raman spectroscopy showed that the stabilization of an exogenous heme ligand, such as CO, occurs differently in human cytoglobin in comparison with Antarctic fish cytoglobins. Furthermore, while it has been extensively reported that human cytoglobin is essentially monomeric and can form an intramolecular disulfide bridge that can influence the ligand binding kinetics, 3D modeling of the Antarctic fish cytoglobins indicates that the cysteine residues are too far apart to form such an intramolecular bridge. Moreover, gel filtration and mass spectrometry reveal the occurrence of non-covalent multimers (up to pentamers) in the Antarctic fish cytoglobins that are formed at low concentrations. Stabilization of these oligomers by disulfide-bridge formation is possible, but not essential. If intermolecular disulfide bridges are formed, they influence the heme-pocket structure, as is shown by EPR measurements.

Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Polynoid scale worms (Polynoidae, Annelida) invaded deep-sea chemosynthesis-based ecosystems approximately 60 million years ago, but little is known about their genetic adaptation to the extreme deep-sea environment. In this study, we reported the first two transcriptomes of deep-sea polynoids (Branchipolynoe pettiboneae, Lepidonotopodium sp.) and compared them with the transcriptome of a shallow-water polynoid (Harmothoe imbricata). We determined codon and amino acid usage, positive selected genes, highly expressed genes and putative duplicated genes. Transcriptome assembly produced 98,806 to 225,709 contigs in the three species. There were more positively charged amino acids (i.e., histidine and arginine) and less negatively charged amino acids (i.e., aspartic acid and glutamic acid) in the deep-sea species. There were 120 genes showing clear evidence of positive selection. Among the 10% most highly expressed genes, there were more hemoglobin genes with high expression levels in both deep-sea species. The duplicated genes related to DNA recombination and metabolism, and gene expression were only enriched in deep-sea species. Deep-sea scale worms adopted two strategies of adaptation to hypoxia in the chemosynthesis-based habitats (i.e., rapid evolution of tetra-domain hemoglobin in Branchipolynoe or high expression of single-domain hemoglobin in Lepidonotopodium sp.).

Genetic Adaptation of Schizothoracine Fish to the Phased Uplifting of the Qinghai-Tibetan Plateau

Many species of Schizothoracine, a subfamily of Cyprinidae, are highly endemic to the Qinghai–Tibetan Plateau (QTP). To characterize the adaptive changes associated with the Schizothoracine expansion at high altitudes, we sequenced tissue transcriptomes of two highland and two subhighland Schizothoracines and analyzed gene evolution patterns by comparing with lowland cyprinids. Phylogenetic tree reconstruction and divergence time estimation indicated that the common ancestor of Schizothoracine fish lived ∼32.7 million years ago (MYA), coinciding with the timing of the first phase of QTP uplifting. Both high- and subhigh-Schizothoracines demonstrated elevated dN/dS ratios in the protein-coding genes compared to lowland cyprinids, from which some biological processes implicated in altitude adaptation were commonly identified. On the other hand, the highland and subhighland lineages presented drastically divergent landscapes of positively selected genes (PSGs), enriched with very different gene ontology (GO) profiles, including those in “sensory organ morphogenesis,” “regulation of protein ubiquitination,” “blood circulation,” and “blood vessel development.” These results indicated different selection pressures imposed on the highland and subhighland lineages of the Schizothoracine subfamily, with a higher number of genes in the high-altitude species involved in adaptations such as sensory perception, blood circulation, and protein metabolism. Our study indicated divergent genetic adaptations in the aquatic species facing the phased uplifting of QTP.

Antioxidant defense system and oxidative status in Antarctic fishes: The sluggish rockcod Notothenia coriiceps versus the active marbled notothen Notothenia rossii

Adaptive responses of antioxidant defense systems (ADS) to changes in increased levels of activity are critical, especially in Antarctic fishes. The benthopelagic marbled notothen (Notothenia rossii) shows higher spontaneous activity than the benthonic and sluggish rockcod (N. coriiceps). Therefore, we hypothesize that species-related responses of ADS would occur to counteract different rates of reactive oxygen species formation in these two Antarctic fish. Here we evaluated ADS and oxidative damage in tissues (brain, gills, liver and white muscle) of the two Antarctic fish. Despite no significant differences in lipid and protein oxidative damage were observed, we actually found species- and tissue-specific differences in ADS. Gill metallothionein-like proteins (MTLP) and liver reduced glutathione (GSH) concentrations were higher in N. coriiceps than in N. rossii. Brain and gill antioxidant capacity against peroxyl radicals (ACAP); gill enzyme [glutamate-cysteine ligase (GSL), superoxide dismutase (SOD) and catalase (CAT)] activity; liver GCL and SOD activity; and white muscle CAT activity were higher in N. rossii than in N. coriiceps. Therefore, the more active fish (N. rossii) maintains higher activities of enzymes involved in superoxide ions (O₂^.-) detoxification and GSH production in peripheral tissues (gills, liver and white muscle). This allows the more active fish (N. rossii) to keep levels of lipid and protein oxidative damage similar to those observed in the sluggish fish (N. coriiceps). It is worth noting that the more active fish also shows a higher brain antioxidant capacity, which could involve other non-enzymatic antioxidants like vitamins C and E. In contrast, N. coriiceps shows lower consumption of non-enzymatic antioxidants in peripheral tissues than N. coriiceps. As hypothesized, our results indicate that differences in ADS profiles between fish species are likely related to their habits and metabolic rates. This would imply in different fish abilities to deal with oxidative stress associated with increasing seawater temperature.

Skin mucus proteins of lumpsucker (Cyclopterus lumpus)

Fish skin mucus serves as a first line of defense against pathogens and external stressors. In this study the proteomic profile of lumpsucker skin mucus was characterized using 2D gels coupled with tandem mass spectrometry. Mucosal proteins were identified by homology searches across the databases SwissProt, NCBInr and vertebrate EST. The identified proteins were clustered into ten groups based on their gene ontology biological process in PANTHER (www.patherdb.org). Calmodulin, cystatin-B, histone H2B, peroxiredoxin1, apolipoprotein A1, natterin-2, 14-3-3 protein, alfa enolase, pentraxin, warm temperature acclimation 65 kDa (WAP65kDa) and heat shock proteins were identified. Several of the proteins are known to be involved in immune and/or stress responses. Proteomic profile established in this study could be a benchmark for differential proteomics studies.

•

A proteome reference map of lumpsucker skin mucus was established.

•

Proteins involved in immune and stress responses were identified in skin mucus of Cyclopterus lumpus.

•

Mucosal proteins identified could be potential biomarkers.

Biodiversity in marine invertebrate responses to acute warming revealed by a comparative multi-omics approach

Understanding species' responses to environmental change underpins our abilities to make predictions on future biodiversity under any range of scenarios. In spite of the huge biodiversity in most ecosystems, a model species approach is often taken in environmental studies. To date, we still do not know how many species we need to study to input into models and inform on ecosystem-level responses to change. In this study, we tested current paradigms on factors setting thermal limits by investigating the acute warming response of six Antarctic marine invertebrates: a crustacean Paraceradocus miersi, a brachiopod Liothyrella uva, two bivalve molluscs, Laternula elliptica, Aequiyoldia eightsii, a gastropod mollusc Marseniopsis mollis and an echinoderm Cucumaria georgiana. Each species was warmed at the rate of 1 °C h^-1 and taken to the same physiological end point (just prior to heat coma). Their molecular responses were evaluated using complementary metabolomics and transcriptomics approaches with the aim of discovering the underlying mechanisms of their resilience or sensitivity to warming. The responses were species-specific; only two showed accumulation of anaerobic end products and three exhibited the classical heat shock response with expression of HSP70 transcripts. These diverse cellular measures did not directly correlate with resilience to heat stress and suggested that each species may have a different critical point of failure. Thus, one unifying molecular mechanism underpinning response to warming could not be assigned, and no overarching paradigm was supported. This biodiversity in response makes future ecosystems predictions extremely challenging, as we clearly need to develop a macrophysiology-type approach to cellular evaluations of the environmental stress response, studying a range of well-rationalized members from different community levels and of different phylogenetic origins rather than extrapolating from one or two arbitrary model species.

Genomic Basis of Adaptive Evolution: The Survival of Amur Ide (Leuciscus waleckii) in an Extremely Alkaline Environment

The Amur ide (Leuciscus waleckii) is a cyprinid fish that is widely distributed in Northeast Asia. The Lake Dali Nur population inhabits one of the most extreme aquatic environments on Earth, with an alkalinity up to 50 mmol/L (pH 9.6), thus providing an exceptional model with which to characterize the mechanisms of genomic evolution underlying adaptation to extreme environments. Here, we developed the reference genome assembly for L. waleckii from Lake Dali Nur. Intriguingly, we identified unusual expanded long terminal repeats (LTRs) with higher nucleotide substitution rates than in many other teleosts, suggesting their more recent insertion into the L. waleckii genome. We also identified expansions in genes encoding egg coat proteins and natriuretic peptide receptors, possibly underlying the adaptation to extreme environmental stress. We further sequenced the genomes of 10 additional individuals from freshwater and 18 from Lake Dali Nur populations, and we detected a total of 7.6 million SNPs from both populations. In a genome scan and comparison of these two populations, we identified a set of genomic regions under selective sweeps that harbor genes involved in ion homoeostasis, acid-base regulation, unfolded protein response, reactive oxygen species elimination, and urea excretion. Our findings provide comprehensive insight into the genomic mechanisms of teleost fish that underlie their adaptation to extreme alkaline environments.

Neofunctionalization of zona pellucida proteins enhances freeze-prevention in the eggs of Antarctic notothenioids

The mechanisms by which the eggs of the Antarctic notothenioid fishes avoid freezing are not fully understood. Zona pellucida proteins (ZPs) are constituents of the chorion which forms a protective matrix surrounding the egg. Here we report occurrence of freezing temperature-related gene expansion and acquisition of unusual ice melting-promoting (IMP) activity in a family of Antarctic notothenioid ZPs (AnnotoZPs). Members of AnnotoZPs are shown to bind with ice and non-colligatively depress the melting point of a solution in a range of 0.26 to 0.65 °C at a moderate concentration. Eggs of zebrafishes expressing an AnnotoZP transgene show improved melting point depression and enhanced survival in freezing conditions. Mutational analyses in a representative AnnotoZP indicate the ZP domain and patches of acidic residues are essential structures for the IMP activity. AnnotoZPs, therefore, represent a group of macromolecules that prevent freezing by a unique ZP-ice interaction mechanism distinct from the known antifreeze proteins.