Transcriptomic and genomic evolution under constant cold in Antarctic notothenioid fish

Transcriptomic annotation of the Chungtien schizothoracin (Ptychobarbus chungtienensis) using Iso-seq and RNA-seq data

The Chungtien schizothoracin (Ptychobarbus chungtienensis), an endangered fish species endemic to the Zhongdian Plateau, remains underexplored in terms of transcriptomic sequencing. This investigation used tissues from five distinct organs (heart, liver, spleen, kidney, and brain) of the Chungtien schizothoracin for PacBio Iso-seq and RNA-seq analyses, yielding a repertoire of 16,598 full-length transcripts spanning lengths from 363 bp to 7,157 bp. Gene family clustering and phylogenetic analysis encompassed a comprehensive set of 13 fish species, all of which were cyprinids, including the zebrafish and the examined species Ptychobarbus chungtienensis. Moreover, the identification of long non-coding RNAs (lncRNAs) and coding sequences was accomplished across all five tissues. Comprehensive analyses of gene expression profiles and differentially expressed genes among the above five tissues were performed. In summary, the obtained full-length transcripts and detailed gene expression profiles of the Chungtien schizothoracin tissues furnish crucial expression data and genetic sequences, laying the groundwork for future investigations and fostering a holistic comprehension of the adaptive mechanisms inherent in the Chungtien schizothoracin under various conditions.

Effects of hibernation on two important contractile tissues in tibetan frogs, Nanorana parkeri: a perspective from transcriptomics and metabolomics approaches

BACKGROUND

In response to seasonal cold and food shortage, the Xizang plateau frogs, Nanorana parkeri (Anura: Dicroglossidae), enter a reversible hypometabolic state where heart rate and oxygen consumption in skeletal muscle are strongly suppressed. However, the effect of winter hibernation on gene expression and metabolic profiling in these two tissues remains unknown. In the present study, we conducted transcriptomic and metabolomic analyses of heart and skeletal muscle from summer- and winter-collected N. parkeri to explore mechanisms involved in seasonal hibernation.

RESULTS

We identified 2407 differentially expressed genes (DEGs) in heart and 2938 DEGs in skeletal muscle. Enrichment analysis showed that shared DEGs in both tissues were enriched mainly in translation and metabolic processes. Of these, the expression of genes functionally categorized as "response to stress", "defense mechanisms", or "muscle contraction" were particularly associated with hibernation. Metabolomic analysis identified 24 and 22 differentially expressed metabolites (DEMs) in myocardium and skeletal muscle, respectively. In particular, pathway analysis showed that DEMs in myocardium were involved in the pentose phosphate pathway, glycerolipid metabolism, pyruvate metabolism, citrate cycle (TCA cycle), and glycolysis/gluconeogenesis. By contrast, DEMs in skeletal muscle were mainly involved in amino acid metabolism.

CONCLUSIONS

In summary, natural adaptations of myocardium and skeletal muscle in hibernating N. parkeri involved transcriptional alterations in translation, stress response, protective mechanisms, and muscle contraction processes as well as metabolic remodeling. This study provides new insights into the transcriptional and metabolic adjustments that aid winter survival of high-altitude frogs N. parkeri.

Biomineralization in a cold environment: Insights from shield compositions and transcriptomics of polar sternaspids (Sternaspidae, Polychaeta)

The survival and physiological functions of polar marine organisms are impacted by global climate changes. Investigation of the adaptation mechanisms underlying biomineralization in polar organisms at low temperatures is important for understanding mineralized organismal sensitivity to climate change. Here, we performed electron probe analysis on the shields of Antarctic polychaete Sternaspis sendalli and Arctic polychaete Sternaspis buzhinskajae (Sternaspidae), and sequenced the transcriptomes of the tissues surrounding shields to examine biomineral characteristics and adaptive mechanisms in persistently cold environments. Compared to the temperate relative species, the relative abundance of iron, phosphorus, calcium, magnesium, nitrogen, sulfur and silicon in two polar sternaspid shields was similar to Sternaspis chinensis. However, the diversity and expression levels of biomineralization-related shell matrix proteins differed between the polar and temperate species, suggesting distinct molecular mechanisms underlying shield formation in cold environments. Tubulin and cyclophilin were upregulated compared to the temperate species. Furthermore, 42 positively selected genes were identified in Antarctic S. sendalli, with functions in cytoskeletal structure, DNA repair, immunity, transcription, translation, protein synthesis, and lipid metabolism. Highly expressed genes in both polar species were associated with cytoskeleton, macromolecular complexes and cellular component biosynthesis. Overall, this study reveals conserved elemental composition yet distinct biomineralization processes in the shields of polar sternaspids. The unique expression of biomineralization related genes and other cold-adaptation related genes provide molecular insights into biomineralization in cold marine environments.

MSTRG3207 promotes apoptosis in zebrafish ZF4 cells via sponging dre-miR-736/bbc3/LOC101885512 axis during cold acclimation

Long non-coding RNAs (lncRNAs) play essential roles in a variety of biological processes. It has been recently reported that lncRNAs can regulate mRNA expression by binding to microRNAs (miRNAs) as competing endogenous RNAs (ceRNAs). However, the involvement of this regulatory mechanism during cold acclimation in fish remains unclear. In this study, we constructed a ceRNA network mediated by lncRNAs in cold-acclimated zebrafish ZF4 cells through bioinformatic analysis of the mRNA, miRNA, and lncRNA profiles obtained from ZF4 cells cultured at 18 °C for 30 days. A previously uncharacterized lncRNA, MSTRG3207, was selected for further analysis. MSTRG3207 was upregulated and dre-miR-736 was downregulated during cold acclimation. MSTRG3207 was cloned by rapid amplification of cDNA ends (RACE) and functionally characterized. The binding of MSTRG3207 to dre-miR-736 was validated by dual-luciferase reporter assay. Under cold acclimation, MSTRG3207 promoted apoptosis by sponging dre-miR-736 and upregulating bbc3 and LOC101885512, two apoptotic genes targeted by dre-miR-736. Taken together, our findings indicate that MSTRG3207 upregulation promotes apoptosis by sponging dre-miR-736 during cold acclimation in fish.

Cold Stress Induces Apoptosis in Silver Pomfret via DUSP-JNK Pathway

We cultured silver pomfret for 20 days, decreasing water temperature from 18 to 8 ℃, and sampled muscle every 5 days. Muscle fiber degeneration and apoptosis began to increase at 13 ℃ detected by HE and TUNEL staining. Further analysis of transcriptome revealed that several apoptosis-related pathways were highly enriched by differentially expressed genes (DEGs). We analyzed 10 DEGs from these pathways by RT-qPCR during the temperature-decreasing process. JNK1, PIDD, CytC, Casp 3, and GADD45 were up-regulated after 15 and 20 days, while DUSP3, JNK2, and PARP genes were down-regulated after 15 and 20 days. DUSP5 was up-regulated from 10 to 20 days, and C-JUN was up-regulated after 20 days. We analyzed apoptosis in PaM cells under different temperatures (26 ℃, 23 ℃, 20 ℃, 17 ℃, and 14 ℃). The cell viability significantly declined from 14 to 20 ℃; the TUNEL and IHC results showed that the apoptosis signal increased with the temperature dropping, especially in 17 ℃ and 14 ℃; DUSP5, JNK1, CytC, C-JUN, Casp 3, and GADD45 were up-regulated at 17 ℃ and 14 ℃, and PIDD was up-regulated at 20 ℃, 17 ℃, and 14 ℃. DUSP3 was up-regulated at 20 ℃ but down-regulated at 17 ℃ and 14 ℃, and PARP was down-regulated at 17 ℃ and 14 ℃. JNK2 was up-regulated at 20 ℃ but down-regulated at 17 ℃ and 14 ℃. Our results suggest that DUSP could help inhibit apoptosis in the initial stage of cold stress, but low temperature could down-regulate it and up-regulate JNK-C-JUN, inducing apoptosis in a later stage. These data provide a basis for the study of the response mechanism of fish to cold.

Genomics and transcriptomics of the Chinese mitten crabs (Eriocheir sinensis)

To gain a deeper understanding of the genetic factors influencing the growth and development of Eriocheir sinensis, a well-known species of hairy crab found in Yangcheng Lake, this study focused on the de novo genome and full-length transcriptome information of the selected subjects. Specifically, Yangcheng Lake hairy crabs were chosen as the experimental samples. Initially, a genome analysis was performed, resulting in the identification of gene fragments with a combined length of 1266,092,319 bp. Subsequently, a transcriptome analysis was conducted on a mixture of tissues from four different sites, namely muscle, brain, eye, and heart, to further investigate the genetic characteristics at the transcriptome level. The Pacific Biosciences (Pacio) single-molecule real-time sequencing system generated a total of 36.93 G sub-fragments and 175,90041 effective inserts. This research contributes to the indirect comprehension of genetic variations underlying individual traits. Furthermore, a comparison of the obtained data with relevant literature emphasizes the advantages of this study and establishes a basis for further investigations on the Chinese mitten crab.

Regulation of PGC-1α of the Mitochondrial Energy Metabolism Pathway in the Gills of Indian Medaka (Oryzias dancena) under Hypothermal Stress

Ectothermic fish exposure to hypothermal stress requires adjusting their metabolic molecular machinery, which was investigated using Indian medaka (Oryzias dancena; 10 weeks old, 2.5 ± 0.5 cm) cultured in fresh water (FW) and seawater (SW; 35‱) at room temperature (28 ± 1 °C). The fish were fed twice a day, once in the morning and once in the evening, and the photoperiod was 12 h:12 h light: dark. In this study, we applied two hypothermal treatments to reveal the mechanisms of energy metabolism via pgc-1α regulation in the gills of Indian medaka; cold-stress (18 °C) and cold-tolerance (extreme cold; 15 °C). The branchial ATP content was significantly higher in the cold-stress group, but not in the cold-tolerance group. In FW- and SW-acclimated medaka, the expression of genes related to mitochondrial energy metabolism, including pgc-1α, prc, Nrf2, tfam, and nd5, was analyzed to illustrate differential responses of mitochondrial energy metabolism to cold-stress and cold-tolerance environments. When exposed to cold-stress, the relative mRNA expression of pgc-1α, prc, and Nrf2 increased from 2 h, whereas that of tfam and nd5 increased significantly from 168 h. When exposed to a cold-tolerant environment, prc was significantly upregulated at 2 h post-cooling in the FW and SW groups, and pgc-1α was significantly upregulated at 2 and 12 h post-cooling in the FW group, while tfam and nd5 were downregulated in both FW and SW fish. Hierarchical clustering revealed gene interactions in the cold-stress group, which promoted diverse mitochondrial energy adaptations, causing an increase in ATP production. However, the cold-tolerant group demonstrated limitations in enhancing ATP levels through mitochondrial regulation via the PGC-1α energy metabolism pathway. These findings suggest that ectothermic fish may develop varying degrees of thermal tolerance over time in response to climate change. This study provides insights into the complex ways in which fish adjust their metabolism when exposed to cold stress, contributing to our knowledge of how they adapt.

Temperature activated transient receptor potential ion channels from Antarctic fishes

Antarctic notothenioid fishes (cryonotothenioids) live in waters that range between -1.86°C and an extreme maximum +4°C. Evidence suggests these fish sense temperature peripherally, but the molecular mechanism of temperature sensation in unknown. Previous work identified transient receptor potential (TRP) channels TRPA1b, TRPM4 and TRPV1a as the top candidates for temperature sensors. Here, cryonotothenioid TRPA1b and TRPV1a are characterized using Xenopus oocyte electrophysiology. TRPA1b and TRPV1a showed heat-evoked currents with Q10s of 11.1 ± 2.2 and 20.5 ± 2.4, respectively. Unexpectedly, heat activation occurred at a threshold of 22.9 ± 1.3°C for TRPA1b and 32.1 ± 0.6°C for TRPV1a. These fish have not experienced such temperatures for at least 15 Myr. Either (1) another molecular mechanism underlies temperature sensation, (2) these fishes do not sense temperatures below these thresholds despite having lethal limits as low as 5°C, or (3) native cellular conditions modify the TRP channels to function at relevant temperatures. The effects of osmolytes, pH, oxidation, phosphorylation, lipids and accessory proteins were tested. No conditions shifted the activity range of TRPV1a. Oxidation in combination with reduced cholesterol significantly dropped activation threshold of TRPA1b to 11.3 ± 2.3°C, it is hypothesized the effect may be due to lipid raft disruption.

Positive and Relaxed Selective Pressures Have Both Strongly Influenced the Evolution of Cryonotothenioid Fishes during Their Radiation in the Freezing Southern Ocean

Evolution in the chronic cold of the Southern Ocean has had a profound influence on the physiology of cryonotothenioid fishes. However, the suite of genetic changes underlying the physiological gains and losses in these fishes is still poorly surveyed. By identifying the genomic signatures of selection, this study aims to identify the functional classes of genes that have been changed following two major physiological transitions: the onset of freezing temperatures and the loss of hemoproteins. Looking at the changes that followed the onset of freezing temperatures, positive selective pressure was found among a set of broadly acting gene regulatory factors, suggesting a route through which cryonotothenioid gene expression has been retooled for life in the cold. Further, genes related to the cell cycle and cellular adhesion were found under positive selection suggesting that both present key challenges to life in freezing waters. By contrast, genes showing signatures of the relaxation of selective pressure showed a narrower biological impact, acting on genes related to mitochondrial function. Finally, although chronic cold-water temperatures appear correlated with substantial genetic change, the loss of hemoproteins resulted in little observable change in protein-coding genes relative to their red-blooded relatives. Combined, the influence of positive and relaxed selection shows that long-term exposure to cold has led to profound changes in cryonotothenioid genomes that may make it challenging for them to adapt to a rapidly changing climate.

Transcriptome wide analyses reveal intraspecific diversity in thermal stress responses of a dominant habitat-forming species

The impact of climate change on biodiversity has stimulated the need to understand environmental stress responses, particularly for ecosystem engineers whose responses to climate affect large numbers of associated organisms. Distinct species differ substantially in their resilience to thermal stress but there are also within-species variations in thermal tolerance for which the molecular mechanisms underpinning such variation remain largely unclear. Intertidal mussels are well-known for their role as ecosystem engineers. First, we exposed two genetic lineages of the intertidal mussel Perna perna to heat stress treatments in air and water. Next, we ran a high throughput RNA sequencing experiment to identify differences in gene expression between the thermally resilient eastern lineage and the thermally sensitive western lineage. We highlight different thermal tolerances that concord with their distributional ranges. Critically, we also identified lineage-specific patterns of gene expression under heat stress and revealed intraspecific differences in the underlying transcriptional pathways in response to warmer temperatures that are potentially linked to the within-species differences in thermal tolerance. Beyond the species, we show how unravelling within-species variability in mechanistic responses to heat stress promotes a better understanding of global evolutionary trajectories of the species as a whole in response to changing climate.

Dusp1 regulates thermal tolerance limits in zebrafish by maintaining mitochondrial integrity

Temperature tolerance restricts the distribution of a species. However, the molecular and cellular mechanisms that set the thermal tolerance limits of an organism are poorly understood. Here, we report on the function of dual-specificity phosphatase 1 (DUSP1) in thermal tolerance regulation. Notably, we found that dusp1 ^-/- zebrafish grew normally but survived within a narrowed temperature range. The higher susceptibility of these mutant fish to both cold and heat challenges was attributed to accelerated cell death caused by aggravated mitochondrial dysfunction and over-production of reactive oxygen species in the gills. The DUSP1-MAPK-DRP1 axis was identified as a key pathway regulating these processes in both fish and human cells. These observations suggest that DUSP1 may play a role in maintaining mitochondrial integrity and redox homeostasis. We therefore propose that maintenance of cellular redox homeostasis may be a key mechanism for coping with cellular thermal stress and that the interplay between signaling pathways regulating redox homeostasis in the most thermosensitive tissue (i.e., gills) may play an important role in setting the thermal tolerance limit of zebrafish.

Population genomics of an icefish reveals mechanisms of glacier-driven adaptive radiation in Antarctic notothenioids

Background

Antarctica harbors the bulk of the species diversity of the dominant teleost fish suborder—Notothenioidei. However, the forces that shape their evolution are still under debate.

Results

We sequenced the genome of an icefish, Chionodraco hamatus, and used population genomics and demographic modelling of sequenced genomes of 52 C. hamatus individuals collected mainly from two East Antarctic regions to investigate the factors driving speciation. Results revealed four icefish populations with clear reproduction separation were established 15 to 50 kya (kilo years ago) during the last glacial maxima (LGM). Selection sweeps in genes involving immune responses, cardiovascular development, and photoperception occurred differentially among the populations and were correlated with population-specific microbial communities and acquisition of distinct morphological features in the icefish taxa. Population and species-specific antifreeze glycoprotein gene expansion and glacial cycle-paced duplication/degeneration of the zona pellucida protein gene families indicated fluctuating thermal environments and periodic influence of glacial cycles on notothenioid divergence.

Conclusions

We revealed a series of genomic evidence indicating differential adaptation of C. hamatus populations and notothenioid species divergence in the extreme and unique marine environment. We conclude that geographic separation and adaptation to heterogeneous pathogen, oxygen, and light conditions of local habitats, periodically shaped by the glacial cycles, were the key drivers propelling species diversity in Antarctica.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01432-x.

Neutral and adaptive loci reveal fine‐scale population structure in Eleginops maclovinus from north Patagonia

Patagonia is an understudied area, especially when it comes to population genomic studies with relevance to fishery management. However, the dynamic and heterogeneous landscape in this area can harbor an important but cryptic genetic population structure. Once such information is revealed, it can be integrated into the management of infrequently investigated species. Eleginops maclovinus is a protandrous hermaphrodite species with economic importance for local communities that are currently managed as a single genetic unit. In this study, we sampled five locations distributed across a salinity cline from Northern Patagonia to investigate the genetic population structure of E. maclovinus. We used restriction site‐associated DNA (RAD) sequencing and outlier tests to obtain neutral and adaptive loci, using F_ST and GEA approaches. We identified a spatial pattern of structuration with gene flow and spatial selection by environmental association. Neutral and adaptive loci showed two and three genetic groups, respectively. The effective population sizes estimated ranged from 572 (Chepu) to 14,454 (Chaitén) and were influenced more by locality than by salinity cline. We found loci putatively associated with salinity suggesting that salinity may act as a selective driver in E. maclovinus populations. These results suggest a complex interaction between genetic drift, gene flow, and natural selection in this area. Our findings also suggest several evolutionary significant units in this area, and the information should be integrated into the management of this species. We discussed the significance of these results for fishery management and suggest future directions to improve our understanding of how E. maclovinus has adapted to the dynamic waters of Northern Patagonia.

Extremophiles in Earth's Deep Seas: A View Toward Life in Exo-Oceans

Humanity's search for extraterrestrial life is a modern manifestation of the exploratory and curious nature that has led us through millennia of scientific discoveries. With the ongoing exploration of extraterrestrial bodies, the potential for discovery of extraterrestrial life has expanded. We may better inform this search through an understanding of how life persists and flourishes on Earth in a myriad of environmental extremes. A significant proportion of our knowledge of extremophiles on Earth comes from studies on deep ocean life. Here, we review and synthesize the range of environmental extremes observed in the deep sea, the life that persists in these extreme conditions, and the biological adaptations utilized by these remarkable life-forms. We also review confirmed and predicted extraterrestrial oceans in our solar system and propose deep-sea sites that may serve as planetary field analog environments. We show that the clever ingenuity of evolution under deep-sea conditions suggests that the plausibility of extraterrestrial life is much greater than previously thought.

Insights Into miRNA-mRNA Regulatory Mechanisms of Cold Adaptation in Gymnocypris eckloni: Ubiquitin-Mediated Proteolysis Is Pivotal for Adaptive Energy Metabolism

This study aimed to understand cold stress adaptations mechanism in fish. Thus, the transcriptional response to cold conditions in Gymnocypris eckloni was evaluated using RNA-seq and microRNA (miRNA)-seq analyses. Low-temperature (LT) group G. eckloni was cultivated outdoors in waters cooled to 2–4°C for 3 weeks, while individuals in the control temperature (CT) group were exposed to 14–16°C. Significantly different responses were observed in both mRNA and miRNA expression profiles, with more mRNAs (1,833 and 1,869 mRNAs were up- and downregulated, respectively) and fewer miRNAs (15 and 6 were up- and downregulated, respectively) observed in the LT group individuals relative to the CT group individuals. A miRNA-mRNA network involved in the regulation of G. eckloni responses to cold stress was constructed; this network included ubiquitin-mediated proteolysis, protein processing, and oxidative phosphorylation. These results provided new insights into mechanisms of cold tolerance by fish, including decreased metabolic activity in addition to proteolysis.

Tissue-Specific and Differential Cold Responses in the Domesticated Cold Tolerant Fugu

Domestication can be defined as the artificial selection in animals to achieve morphological, physiological, and developmental conformity to human needs, with the aim of improving various limitations in species under a human feeding environment. The future sustainability of aquaculture may rely partly on the availability of numerous domesticated fish species. However, the underlying adaptive mechanisms that result in the domestication of fish are still unclear. Because they are poikilothermic, temperature is a key environmental element that affects the entire life of fish, so studying the association between physiological and behavioral changes in low-temperature domesticated fish can provide a model for understanding the response mechanisms of fish under cold stress. Through 5 generations and 10 years of artificial selection at low temperatures, we used cold-tolerant fugu as a biological model to compare transcriptome changes in brain and liver tissues to study the effects of cold stress on fish. It was found that the expression of genes such as apoptosis, p53, oxidative phosphorylation, and mitochondrial β-oxidation in the brain of cold-tolerant fugu was significantly lower than the wild type due to cold stress, while excessive energy metabolism would lead to the production of reactive oxygen species (ROS) and exacerbate the brain damage, thus causing rollover and coma. Meanwhile, under cold stress, the signaling pathways involved in glycogenolysis and lipid metabolism, such as insulin signaling, adipocytokines, and mTOR signaling pathways, were significantly up-regulated in the liver of cold-tolerant fugu. Although the mitochondrial β-oxidation pathway was increased in cold-tolerant fugu liver tissues, the transcriptome was not enriched in apoptotic. These phenomena predict that in response to low-temperature conditions, cold-tolerant fugu employs a dynamic inter-organ metabolic regulation strategy to cope with cold stress and reduce damage to brain tissues.

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

Toll-Like Receptor Evolution: Does Temperature Matter?

Toll-like receptors (TLRs) recognize conserved pathogen-associated molecular patterns (PAMPs) and are an ancient and well-conserved group of pattern recognition receptors (PRRs). The isolation of the Antarctic continent and its unique teleost fish and microbiota prompted the present investigation into Tlr evolution. Gene homologues of tlr members in teleosts from temperate regions were present in the genome of Antarctic Nototheniidae and the non-Antarctic sister lineage Bovichtidae. Overall, in Nototheniidae apart from D. mawsoni, no major tlr gene family expansion or contraction occurred. Instead, lineage and species-specific changes in the ectodomain and LRR of Tlrs occurred, particularly in the Tlr11 superfamily that is well represented in fish. Positive selective pressure and associated sequence modifications in the TLR ectodomain and within the leucine-rich repeats (LRR), important for pathogen recognition, occurred in Tlr5, Tlr8, Tlr13, Tlr21, Tlr22, and Tlr23 presumably associated with the unique Antarctic microbiota. Exposure to lipopolysaccharide (Escherichia coli O111:B4) Gram negative bacteria did not modify tlr gene expression in N. rossii head–kidney or anterior intestine, although increased water temperature (+4°C) had a significant effect.

OUP accepted manuscript

Animals rely on their sensory systems to inform them of ecologically relevant environmental variation. In the Southern Ocean, the thermal environment has remained between −1.9 and 5 °C for 15 Myr, yet we have no knowledge of how an Antarctic marine organism might sense their thermal habitat as we have yet to discover a thermosensitive ion channel that gates (opens/closes) below 10 °C. Here, we investigate the evolutionary dynamics of transient receptor potential (TRP) channels, which are the primary thermosensors in animals, within cryonotothenioid fishes—the dominant fish fauna of the Southern Ocean. We found cryonotothenioids have a similar complement of TRP channels as other teleosts (∼28 genes). Previous work has shown that thermosensitive gating in a given channel is species specific, and multiple channels act together to sense the thermal environment. Therefore, we combined evidence of changes in selective pressure, gene gain/loss dynamics, and the first sensory ganglion transcriptome in this clade to identify the best candidate TRP channels that might have a functional dynamic range relevant for frigid Antarctic temperatures. We concluded that TRPV1a, TRPA1b, and TRPM4 are the likeliest putative thermosensors, and found evidence of diversifying selection at sites across these proteins. We also put forward hypotheses for molecular mechanisms of other cryonotothenioid adaptations, such as reduced skeletal calcium deposition, sensing oxidative stress, and unusual magnesium homeostasis. By completing a comprehensive and unbiased survey of these genes, we lay the groundwork for functional characterization and answering long-standing thermodynamic questions of thermosensitive gating and protein adaptation to low temperatures.

Specific immunity proteomic profile of the skin mucus of Antarctic fish Chionodraco hamatus and Notothenia coriiceps

The white-blooded Antarctic icefish is the only known vertebrate lacking oxygen-transporting haemoglobins. Fish skin mucus, as the first line of defence against pathogens, can reflect fish welfare. In this study, we analysed the skin mucus proteome profiles of the two Antarctic fish species, the white-blooded Antarctic icefish, Chionodraco hamatus, and the red-blooded Antarctic fish, Notothenia coriiceps, unfolding the different proteins by liquid chromatography coupled with tandem mass spectrometry isobaric tags for relative and absolute quantitation (iTRAQ) technology. Of the 4444 totally identified proteins, 227 differentially expressed proteins (DEPs) were found in the comparison between C. hamatus and N. coriiceps, of which 121 were upregulated and 106 were downregulated in the icefish. In the Kyoto Encyclopedia of Genes and Genomes pathway annotation, we found two pathways "Legionellosis" and "Complement and coagulation cascades" were significantly enriched, among of which innate immune candidate proteins such as C3, CASP1, ASC, F3 and C9 were significantly upregulated, suggesting their important roles in C. hamatus immune system. Additionally, the DEP protein-protein interaction network analysis and "Response to stress" GO category provided candidate biomarkers for deep understanding of the distinct immune response of the two Antarctic fish underlying the cold adaptation.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Effect of heat stress on the antioxidant defense system and erythrocyte morphology of Antarctic fishes

This study analyzed the effect of thermal stress on erythrocytes of Notothenia rossii and Notothenia coriiceps, abundant notothenioids in Admiralty Bay, Antarctic Peninsula. In both species, the antioxidant defense system enzymes, superoxide dismutase, catalase, glutathione peroxidase, glutathione-S transferase, glutathione reductase were punctually altered (8°C for 1, 3 and 6 days) in erythrocytes, indicating that these markers are not ideal for termal stress. However, under the influence of thermal stress, morphological changes in Notothenia coriiceps erythrocytes were observed at all exposure times (1, 3 and 6 days at 8°C), and in Notothenia rossii occurred in 6 days. These results suggest that Notothenia corriceps presents a lower tolerance to thermal stress at 8°C for up to 6 days, since the cellular and nuclear alterations recorded are pathological and may be deleterious to the cells. Among the morphological markers analyzed in this work, we believe that the shape change and nuclear bubble formation may be good stress biomarkers in erythrocytes of Notothenia rossii and Notothenia coriiceps.

Correspondence of aCGH and long-read genome assembly for detection of copy number differences: A proof-of-concept with cichlid genomes

Copy number variation is an important source of genetic variation, yet data are often lacking due to technical limitations for detection given the current genome assemblies. Our goal is to demonstrate the extent to which an array-based platform (aCGH) can identify genomic loci that are collapsed in genome assemblies that were built with short-read technology. Taking advantage of two cichlid species for which genome assemblies based on Illumina and PacBio are available, we show that inter-species aCGH log2 hybridization ratios correlate more strongly with inferred copy number differences based on PacBio-built genome assemblies than based on Illumina-built genome assemblies. With regard to inter-species copy number differences of specific genes identified by each platform, the set identified by aCGH intersects to a greater extent with the set identified by PacBio than with the set identified by Illumina. Gene function, according to Gene Ontology analysis, did not substantially differ among platforms, and platforms converged on functions associated with adaptive phenotypes. The results of the current study further demonstrate that aCGH is an effective platform for identifying copy number variable sequences, particularly those collapsed in short read genome assemblies.

The macronuclear genome of the Antarctic psychrophilic marine ciliate Euplotes focardii reveals new insights on molecular cold adaptation

The macronuclear (MAC) genomes of ciliates belonging to the genus Euplotes species are comprised of numerous small DNA molecules, nanochromosomes, each typically encoding a single gene. These genomes are responsible for all gene expression during vegetative cell growth. Here, we report the analysis of the MAC genome from the Antarctic psychrophile Euplotes focardii. Nanochromosomes containing bacterial sequences were not found, suggesting that phenomena of horizontal gene transfer did not occur recently, even though this ciliate species has a substantial associated bacterial consortium. As in other euplotid species, E. focardii MAC genes are characterized by a high frequency of translational frameshifting. Furthermore, in order to characterize differences that may be consequent to cold adaptation and defense to oxidative stress, the main constraints of the Antarctic marine microorganisms, we compared E. focardii MAC genome with those available from mesophilic Euplotes species. We focussed mainly on the comparison of tubulin, antioxidant enzymes and heat shock protein (HSP) 70 families, molecules which possess peculiar characteristic correlated with cold adaptation in E. focardii. We found that α-tubulin genes and those encoding SODs and CATs antioxidant enzymes are more numerous than in the mesophilic Euplotes species. Furthermore, the phylogenetic trees showed that these molecules are divergent in the Antarctic species. In contrast, there are fewer hsp70 genes in E. focardii compared to mesophilic Euplotes and these genes do not respond to thermal stress but only to oxidative stress. Our results suggest that molecular adaptation to cold and oxidative stress in the Antarctic environment may not only be due to particular amino acid substitutions but also due to duplication and divergence of paralogous genes.

Different ways to play it cool: Transcriptomic analysis sheds light on different activity patterns of three amphipod species under long-term cold exposure

Species of littoral freshwater environments in regions with continental climate experience pronounced seasonal temperature changes. Coping with long cold winters and hot summers requires specific physiological and behavioural adaptations. Endemic amphipods of Lake Baikal, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, show high metabolic activity throughout the year; E. verrucosus even reproduces in winter. In contrast, the widespread Holarctic amphipod Gammarus lacustris overwinters in torpor. This study investigated the transcriptomic hallmarks of E. verrucosus, E. cyaneus and G. lacustris exposed to low water temperatures. Amphipods were exposed to 1.5°C and 12°C (corresponding to the mean winter and summer water temperatures, respectively, in the Baikal littoral) for one month. At 1.5°C, G. lacustris showed upregulation of ribosome biogenesis and mRNA processing genes, as well as downregulation of genes related to growth, reproduction and locomotor activity, indicating enhanced energy allocation to somatic maintenance. Our results suggest that the mitogen-activated protein kinase (MAPK) signalling pathway is involved in the preparation for hibernation; downregulation of the actin cytoskeleton pathway genes could relate to the observed low locomotor activity of G. lacustris at 1.5°C. The differences between the transcriptomes of E. verrucosus and E. cyaneus from the 1.5°C and 12°C exposures were considerably smaller than for G. lacustris. In E. verrucosus, cold-exposure triggered reproductive activity was indicated by upregulation of respective genes, whereas in E. cyaneus, genes related to mitochondria functioning were upregulated, indicating cold compensation in this species. Our data elucidate the molecular characteristics behind the different adaptations of amphipod species from the Lake Baikal area to winter conditions.

Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog

Temperature is a critical driver of ectotherm life‐history strategies, whereby a warmer environment is associated with increased growth, reduced longevity and accelerated senescence. Increasing evidence indicates that thermal adaptation may underlie such life‐history shifts in wild populations. Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) can help uncover the molecular mechanisms of temperature‐driven variation in growth, longevity and senescence. However, our understanding of these mechanisms is still limited, which reduces our ability to predict the response of non‐model ectotherms to global temperature change.

In this study, we examined the potential role of thermal adaptation in clinal shifts of life‐history traits (i.e. life span, senescence rate and recruitment) in the Columbia spotted frog Rana luteiventris along a broad temperature gradient in the western United States.

We took advantage of extensive capture–recapture datasets of 20,033 marked individuals from eight populations surveyed annually for 14–18 years to examine how mean annual temperature and precipitation influenced demographic parameters (i.e. adult survival, life span, senescence rate, recruitment and population growth). After showing that temperature was the main climatic predictor influencing demography, we used RAD‐seq data (50,829 SNPs and 6,599 putative CNVs) generated for 352 individuals from 31 breeding sites to identify the genomic signatures of thermal adaptation.

Our results showed that temperature was negatively associated with annual adult survival and reproductive life span and positively associated with senescence rate. By contrast, recruitment increased with temperature, promoting the long‐term viability of most populations. These temperature‐dependent demographic changes were associated with strong genomic signatures of thermal adaptation. We identified 148 SNP candidates associated with temperature including three SNPs located within protein‐coding genes regulating resistance to cold and hypoxia, immunity and reproduction in ranids. We also identified 39 CNV candidates (including within 38 transposable elements) for which normalized read depth was associated with temperature.

Our study indicates that both SNPs and structural variants are associated with temperature and could eventually be found to play a functional role in clinal shifts in senescence rate and life‐history strategies in R. luteiventris. These results highlight the potential role of different sources of molecular variation in the response of ectotherms to environmental temperature variation in the context of global warming.

Induced Torpor as a Countermeasure for Low Dose Radiation Exposure in a Zebrafish Model

The development of the Artemis programme with the goal of returning to the moon is spurring technology advances that will eventually take humans to Mars and herald a new era of interplanetary space travel. However, long-term space travel poses unique challenges including exposure to ionising radiation from galactic cosmic rays and potential solar particle events, exposure to microgravity and specific nutritional challenges arising from earth independent exploration. Ionising radiation is one of the major obstacles facing future space travel as it can generate oxidative stress and directly damage cellular structures such as DNA, in turn causing genomic instability, telomere shortening, extracellular-matrix remodelling and persistent inflammation. In the gastrointestinal tract (GIT) this can lead to leaky gut syndrome, perforations and motility issues, which impact GIT functionality and affect nutritional status. While current countermeasures such as shielding from the spacecraft can attenuate harmful biological effects, they produce harmful secondary particles that contribute to radiation exposure. We hypothesised that induction of a torpor-like state would confer a radioprotective effect given the evidence that hibernation extends survival times in irradiated squirrels compared to active controls. To test this hypothesis, a torpor-like state was induced in zebrafish using melatonin treatment and reduced temperature, and radiation exposure was administered twice over the course of 10 days. The protective effects of induced-torpor were assessed via RNA sequencing and qPCR of mRNA extracted from the GIT. Pathway and network analysis were performed on the transcriptomic data to characterise the genomic signatures in radiation, torpor and torpor + radiation groups. Phenotypic analyses revealed that melatonin and reduced temperature successfully induced a torpor-like state in zebrafish as shown by decreased metabolism and activity levels. Genomic analyses indicated that low dose radiation caused DNA damage and oxidative stress triggering a stress response, including steroidal signalling and changes to metabolism, and cell cycle arrest. Torpor attenuated the stress response through an increase in pro-survival signals, reduced oxidative stress via the oxygen effect and detection and removal of misfolded proteins. This proof-of-concept model provides compelling initial evidence for utilizing an induced torpor-like state as a potential countermeasure for radiation exposure.

Chromosome-level assembly of the mangrove plant Aegiceras corniculatum genome generated through Illumina, PacBio and Hi-C sequencing technologies

Aegiceras corniculatum is a major mangrove plant species adapted to waterlogging and saline conditions, grows in the coastal intertidal zone of tropical and subtropical regions. Here, we present a chromosome-level genome assembly of A. corniculatum by incorporating PacBio long-read sequencing and Hi-C technology. The results showed that the PacBio draft genome size is 906.63 Mb. Hi-C scaffolding anchored 885.06 Mb contigs (97.62% of draft assembly) onto 24 pseudochromosomes. The contig N50 and scaffold N50 were 7.1 Mb and 37.74 Mb, respectively. Out of 40,727 protein-coding genes predicted in the study, 89% have functional annotations in public databases. We also showed that of the 603.93 Mb repetitive sequences predicted in the assembled genome, long terminal repeat retrotransposons constitute 41.52%. The genome evolution analysis showed that the A. corniculatum genome experienced two whole-genome duplication events and shared the ancient γ whole-genome triplication event. A comparative genomic analysis revealed an incidence of expansion in 1,488 gene families associated with essential metabolism and biosynthetic pathways, including photosynthesis, oxidative phosphorylation, phenylalanine, glyoxylate, dicarboxylate metabolism, and DNA replication, which probably constitute adaptation traits that allow the A. corniculatum to survive in the intertidal zone. Also, the systematic characterization of genes associated with flavonoid biosynthesis pathway and the AcNHX gene family conducted in this study will provide insight into the adaptation mechanism of A. corniculatum to intertidal environments. The high-quality genome reported here can provide historical insights into genomic transformations that support the survival of A. corniculatum under harsh intertidal habitats.

Thermal adaptation rather than demographic history drives genetic structure inferred by copy number variants in a marine fish

Increasing evidence shows that structural variants represent an overlooked aspect of genetic variation with consequential evolutionary roles. Among those, copy number variants (CNVs), including duplicated genomic regions and transposable elements (TEs), may contribute to local adaptation and/or reproductive isolation among divergent populations. Those mechanisms suppose that CNVs could be used to infer neutral and/or adaptive population genetic structure, whose study has been restricted to microsatellites, mitochondrial DNA and Amplified fragment length polymorphism markers in the past and more recently the use of single nucleotide polymorphisms (SNPs). Taking advantage of recent developments allowing CNV analysis from RAD-seq data, we investigated how variation in fitness-related traits, local environmental conditions and demographic history are associated with CNVs, and how subsequent copy number variation drives population genetic structure in a marine fish, the capelin (Mallotus villosus). We collected 1538 DNA samples from 35 sampling sites in the north Atlantic Ocean and identified 6620 putative CNVs. We found associations between CNVs and the gonadosomatic index, suggesting that six duplicated regions could affect female fitness by modulating oocyte production. We also detected 105 CNV candidates associated with water temperature, among which 20% corresponded to genomic regions located within the sequence of protein-coding genes, suggesting local adaptation to cold water by means of gene sequence amplification. We also identified 175 CNVs associated with the divergence of three previously defined parapatric glacial lineages, of which 24% were located within protein-coding genes, making those loci potential candidates for reproductive isolation. Lastly, our analyses unveiled a hierarchical, complex CNV population structure determined by temperature and local geography, which was in stark contrast to that inferred based on SNPs in a previous study. Our findings underline the complementarity of those two types of genomic variation in population genomics studies.

Overexpression of the V-ATPase c subunit gene from Antarctic notothenioid fishes enhances freezing tolerance in transgenic Arabidopsis plants

Theoretical and experimental studies have demonstrated that temperature is an important environmental factor that affects the regional distribution of plants. However, how to modify the distribution pattern of plants in different regions is a focus of current research. Obtain the information of cold tolerance genes from cold tolerance species, cloning genes with real cold tolerance effects is one of the most important ways to find the genes related to cold tolerance. In this study, we investigated whether transferring the VHA-c gene from Antarctic notothenioid fishes into Arabidopsis enhances freezing tolerance of Arabidopsis. The physiological response and molecular changes of VHA-c overexpressing pedigree and wildtype Arabidopsis were studied at -20 °C. The results showed that the malondialdehyde (MDA) and membrane leakage rates of WT plants were significantly higher than those of VHA-c8 and VHA-c11 plants, but the soluble sugar, soluble protein, proline and ATP contents of WT plants were significantly lower than those of VHA-c8 and VHA-c11 plants under -20 °C freezing treatment. The survival rate, VHA-c gene expression level and VHA-c protein contents of WT plants were significantly lower than those of VHA-c8 and VHA-c11 plants under -20 °C freezing treatment. Correlation analysis showed that ATP content was significantly negatively correlated with MDA and membrane leakage rate, and positively correlated with soluble sugar, soluble protein and proline content under -20 °C freezing treatment. These results demonstrated that overexpression of the VHA-c gene provided strong freezing tolerance to Arabidopsis by increasing the synthesis of ATP and improved the adaptability of plants in low temperature environment.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

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.

Gene Copy Number Variation Does Not Reflect Structure or Environmental Selection in Two Recently Diverged California Populations of Suillus brevipes

Gene copy number variation across individuals has been shown to track population structure and be a source of adaptive genetic variation with significant fitness impacts. In this study, we report opposite results for both predictions based on the analysis of gene copy number variants (CNVs) of Suillus brevipes, a mycorrhizal fungus adapted to coastal and montane habitats in California. In order to assess whether gene copy number variation mirrored population structure and selection in this species, we investigated two previously studied locally adapted populations showing a highly differentiated genomic region encompassing a gene predicted to confer salt tolerance. In addition, we examined whether copy number in the genes related to salt homeostasis was differentiated between the two populations. Although we found many instances of CNV regions across the genomes of S. brevipes individuals, we also found CNVs did not recover population structure and known salt-tolerance-related genes were not under selection across the coastal population. Our results contrast with predictions of CNVs matching single-nucleotide polymorphism divergence and showed CNVs of genes for salt homeostasis are not under selection in S. brevipes.

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 reveals molecular mechanisms responsive to acute cold stress in the tropical stenothermal fish tiger barb (Puntius tetrazona)

BACKGROUND

Tropical stenothermal fish exhibit special tolerance and response to cold stress. However current knowledge of the molecular mechanisms response to cold stress in aquatic ectotherms is largely drawn from eurythermal or extreme stenothermal species. The tiger barb Puntius tetrazona is a tropical stenothermal fish, with great popularity in aquarium trade and research.

RESULTS

To investigate the response mechanism of P. tetrazona to low temperature, fish were exposed to increasing levels of acute cold stress. Histopathological analysis showed that the brain, gill, liver and muscle tissues appeared serious damage after cold stress (13 °C). Brain, gill, liver and muscle tissues from control (CTRL) groups (27 °C) and COLD stress groups (13 °C) of eight-month fish (gender-neutral) were sampled and assessed for transcriptomic profiling by high-throughput sequencing. 83.0 Gb of raw data were generated, filtered and assembled for de novo transcriptome assembly. According to the transcriptome reference, we obtained 392,878 transcripts and 238,878 unigenes, of which 89.29% of the latter were annotated. There were 23,743 differently expressed genes (DEGs) been filtered from four pairs of tissues (brain, gill, liver and muscle) between these cold stress and control groups. These DEGs were mainly involved in circadian entrainment, circadian rhythm, biosynthesis of steroid and fatty acid. There were 64 shared DEGs between the four pairs of groups, and five were related to ubiquitylation/deubiquitylation. Our results suggested that ubiquitin-mediated protein degradation might be necessary for tropical stenothermal fish coping with acute cold stress. Also, the significant cold-induced expression of heat shock 70 kDa protein (HSP70) and cold-induced RNA-binding protein (CIRBP) was verified. These results suggested that the expression of the molecular chaperones HSP70 and CIRBP in P. tetrazona might play a critical role in coping with acute cold stress.

CONCLUSIONS

This is the first transcriptome analysis of P. tetrazona using RNA-Seq technology. Novel findings about tropical stenothermal fish under cold stress (such as HSP70 and CIRBP genes) are presented here. This study contributes new insights into the molecular mechanisms of tropical stenothermal species response to acute cold stress.

Insights into the Function and Evolution of Taste 1 Receptor Gene Family in the Carnivore Fish Gilthead Seabream (Sparus aurata)

A plethora of molecular and functional studies in tetrapods has led to the discovery of multiple taste 1 receptor (T1R) genes encoding G-protein coupled receptors (GPCRs) responsible for sweet (T1R2 + T1R3) and umami (T1R1 + T1R3) taste. In fish, the T1R gene family repertoires greatly expanded because of several T1R2 gene duplications, and recent studies have shown T1R2 functional divergence from canonical mammalian sweet taste perceptions, putatively as an adaptive mechanism to develop distinct feeding strategies in highly diverse aquatic habitats. We addressed this question in the carnivore fish gilthead seabream (Sparus aurata), a model species of aquaculture interest, and found that the saT1R gene repertoire consists of eight members including saT1R1, saT1R3 and six saT1R2a-f gene duplicates, adding further evidence to the evolutionary complexity of fishT1Rs families. To analyze saT1R taste functions, we first developed a stable gene reporter system based on Ca²⁺-dependent calcineurin/NFAT signaling to examine specifically in vitro the responses of a subset of saT1R heterodimers to L-amino acids (L-AAs) and sweet ligands. We show that although differentially tuned in sensitivity and magnitude of responses, saT1R1/R3, saT1R2a/R3 and saT1R2b/R3 may equally serve to transduce amino acid taste sensations. Furthermore, we present preliminary information on the potential involvement of the Gi protein alpha subunits saGαi1 and saGαi2 in taste signal transduction.

Global gene expression profile under low-temperature conditions in the brain of the grass carp (Ctenopharyngodon idellus)

Grass carp is an important commercial fish widely cultivated in China. A wide range of temperatures, particularly extremely low temperatures, have dramatic effects on the aquaculture of this teleost. However, relatively few studies have characterized the molecular responses of grass carp exposed to acute cooling in natural environment. Here, we investigated the transcriptome profiles of the grass carp brain in response to cooling. Through regulation pattern analyses, we identified 2,513 differentially expressed genes (DEGs) that responded to moderate cold stress (12°C), while 99 DEGs were induced by severe low temperature (4°C).The pathway analyses revealed that the DEGs sensitive to moderate cold were largely enriched in steroid biosynthesis, spliceosome, translation, protein metabolism, phagosome, gap junction and estrogen signaling pathways. Additionally, we discerned genes most likely involved in low temperature tolerance, of which the MAPK signaling pathway was dominantly enriched. Further examination and characterization of the candidate genes may help to elucidate the mechanisms underpinning extreme plasticity to severe cold stress in grass carp.

Transcriptomic evidences of local thermal adaptation for the native fish Colossoma macropomum (Cuvier, 1818)

Brazil has five climatically distinct regions, with an annual average temperature difference up to 14 ºC between the northern and southern extremes. Environmental variation of this magnitude can lead to new genetic patterns among farmed fish populations. Genetically differentiated populations of tambaqui (Colossoma macropomum Cuvier, 1818), an important freshwater fish for Brazilian continental aquaculture, may be associated with regional adaptation. In this study, we selected tambaquis raised in two thermally distinct regions, belonging to different latitudes, to test this hypothesis. De novo transcriptome analysis was performed to compare the significant differences of genes expressed in the liver of juvenile tambaqui from a northern population (Balbina) and a southeastern population (Brumado). In total, 2,410 genes were differentially expressed (1,196 in Balbina and 1,214 in Brumado). Many of the genes are involved in a multitude of biological functions such as biosynthetic processes, homeostasis, biorhythm, immunity, cell signaling, ribosome biogenesis, modification of proteins, intracellular transport, structure/cytoskeleton, and catalytic activity. Enrichment analysis based on biological networks showed a different protein interaction profile for each population, whose encoding genes may play potential functions in local thermal adaptation of fish to their respective farming environments.

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

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

Effects of di-(2-ethylhexyl) phthalate on Transcriptional Expression of Cellular Protection-Related HSP60 and HSP67B2 Genes in the Mud Crab Macrophthalmus japonicus

Di-2-ethylhexyl phthalate (DEHP) has attracted attention as an emerging dominant phthalate contaminant in marine sediments. Macrophthalmus japonicus, an intertidal mud crab, is capable of tolerating variations in water temperature and sudden exposure to toxic substances. To evaluate the potential effects of DEHP toxicity on cellular protection, we characterized the partial open reading frames of the stress-related heat shock protein 60 (HSP60) and small heat shock protein 67B2 (HSP67B2) genes of M. japonicus and further investigated the molecular effects on their expression levels after exposure to DEHP. Putative HSP60 and small HSP67B2 proteins had conserved HSP-family protein sequences with different C-terminus motifs. Phylogenetic analysis indicated that M. japonicus HSP60 (Mj-HSP60) and M. Japonicus HSP67B2 (Mj-HSP67B2) clustered closely with Eriocheir sinensis HSP60 and Penaeus vannamei HSP67B2, respectively. The tissue distribution of Heat shock proteins (HSPs) was the highest in the gonad for Mj-HSP60 and in the hepatopancreas for Mj-HSP67B2. The expression of Mj-HSP60 Messenger Ribonucleic Acid (mRNA) increased significantly at day 1 after exposure to all doses of DEHP, and then decreased in a dose-dependent and exposure time-dependent manner in the gills and hepatopancreas. Mj-HSP67B2 transcripts were significantly upregulated in both tissues at all doses of DEHP and at all exposure times. These results suggest that cellular immune protection could be disrupted by DEHP toxicity through transcriptional changes to HSPs in crustaceans. Small and large HSPs might be differentially involved in responses against environmental stressors and in detoxification in M. japonicus crabs.

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.