The genome sequence of the Antarctic bullhead notothen reveals evolutionary adaptations to a cold environment

A conserved phenylalanine motif among teleost fish provides insight for improving electromagnetic perception

Magnetoreceptive biology as a field remains relatively obscure; compared with the breadth of species believed to sense magnetic fields, it remains under-studied. Here, we present grounds for the expansion of magnetoreception studies among teleosts. We begin with the electromagnetic perceptive gene (EPG) from Kryptopterus vitreolus and expand to identify 72 teleosts with homologous proteins containing a conserved three-phenylalanine (3F) motif. Phylogenetic analysis provides insight as to how EPG may have evolved over time and indicates that certain clades may have experienced a loss of function driven by different fitness pressures. One potential factor is water type with freshwater fish significantly more likely to possess the functional motif version (FFF), and saltwater fish to have the non-functional variant (FXF). It was also revealed that when the 3F motif from the homologue of Brachyhypopomus gauderio (B.g.) is inserted into EPG-EPG(B.g.)-the response (as indicated by increased intracellular calcium) is faster. This indicates that EPG has the potential to be engineered to improve upon its response and increase its utility to be used as a controller for specific outcomes.

A conserved phenylalanine motif among Teleost fish provides insight for improving electromagnetic perception

Magnetoreceptive biology as a field remains relatively obscure; compared to the breadth of species believed to sense magnetic fields, it remains under-studied. Here, we present grounds for the expansion of magnetoreception studies among Teleosts. We begin with the electromagnetic perceptive gene (EPG) from Kryptopterus vitreolus and expand to identify 72 Teleosts with homologous proteins containing a conserved three-phenylalanine (3F) motif. Phylogenetic analysis provides insight as to how EPG may have evolved over time, and indicates that certain clades may have experienced a loss of function driven by different fitness pressures. One potential factor is water type with freshwater fish significantly more likely to possess the functional motif version (FFF), and saltwater fish to have the non-functional variant (FXF). It was also revealed that when the 3F motif from the homolog of Brachyhypopomus gauderio (B.g.) is inserted into EPG - EPG(B.g.) - the response (as indicated by increased intracellular calcium) is faster. This indicates that EPG has the potential to be engineered to improve upon its response and increase its utility to be used as a controller for specific outcomes.

Chromosome-level genome assembly and annotation of the Antarctica whitefin plunderfish Pogonophryne albipinna

The Antarctic whitefin plunderfish Pogonophryne albipinna belongs to the family Artedidraconidae, a key component of Antarctic benthic ecosystems within the order Perciformes and the suborder Notothenioidei. While genome research on P. albipinna using short-read sequencing is available, high-quality genome assembly and annotation employing long-read sequencing have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for P. albipinna, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The resulting genome assembly spans approximately 1.07 Gb, with a longest scaffold measuring 59.39 Mb and an N50 length of 41.76 Mb. Of the 1,111 Hi-C scaffolds, 23 exceeded 10 Mb and were thus classified as chromosome-level. BUSCO completeness was assessed at 95.6%. The assembled genome comprises 50.68% repeat sequences, and a total of 31,128 protein-coding genes were predicted. This study will enhance our understanding of the genomic characteristics of cryonotothenioids and facilitate comparative analyses of their adaptation and evolution in extreme environments.

De novo transcriptome assembly of an Antarctic nematode for the study of thermal adaptation in marine parasites

Understanding the genomic underpinnings of thermal adaptation is a hot topic in eco-evolutionary studies of parasites. Marine heteroxenous parasites have complex life cycles encompassing a free-living larval stage, an ectothermic intermediate host and a homeothermic definitive host, thus representing compelling systems for the study of thermal adaptation. The Antarctic anisakid Contracaecum osculatum sp. D is a marine parasite able to survive and thrive both at very cold and warm temperatures within the environment and its hosts. Here, a de novo transcriptome of C. osculatum sp. D was generated for the first time, by performing RNA-Seq experiments on a set of individuals exposed to temperatures experienced by the nematode during its life cycle. The analysis generated 425,954,724 reads, which were assembled and then annotated. The high-quality assembly was validated, achieving over 88% mapping against the transcriptome. The transcriptome of this parasite will represent a valuable genomic resource for future studies aimed at disentangling the genomic architecture of thermal tolerance and metabolic pathways related to temperature stress.

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.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

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.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

Signatures of selection in indigenous Chinese cattle genomes reveal adaptive genes and genetic variations to cold climate

Cold climate shapes the genome of animals and drives them to carry sufficient genetic variations to adapt to changes in temperature. However, limited information is available about the genome-wide pattern of adaptations to cold environments in cattle. In the present study, we used 777K SNP genotyping (15 Chinese cattle breeds, 198 individuals) and whole genome resequencing data (54 cattle breeds of the world, 432 individuals) to disentangle divergent selection signatures, especially between the cold-adapted (annual average temperature of habitat, 6.24 °C to 10.3 °C) and heat-adapted (20.2 °C to 24.73 °C) Chinese native cattle breeds. Genomic analyses revealed a set of candidate genes (e.g., UQCR11, DNAJC18, EGR1, and STING1) were functionally associated with thermogenesis and energy metabolism. We also characterized the adaptive loci of cattle exposed to cold temperatures. Our study finds new candidate genes and provides new insights into adaptations to cold climates in cattle.

Novel genome sequence of Chinese cavefish (Triplophysa rosa) reveals pervasive relaxation of natural selection in cavefish genomes

All cavefishes, living exclusively in caves across the globe, exhibit similar phenotypic traits, including the characteristic loss of eyes. To understand whether such phenotypic convergence shares similar genomic bases, here we investigated genome-wide evolutionary signatures of cavefish phenotypes by comparing whole-genome sequences of three pairs of cavefishes and their surface fish relatives. Notably, we newly sequenced and generated a whole-genome assembly of the Chinese cavefish Triplophysa rosa. Our comparative analyses revealed several shared features of cavefish genome evolution. Cavefishes had lower mutation rates than their surface fish relatives. In contrast, the ratio of nonsynonymous to synonymous substitutions (ω) was significantly elevated in cavefishes compared to in surface fishes, consistent with the relaxation of purifying selection. In addition, cavefish genomes had an increased mutational load, including mutations that alter protein hydrophobicity profiles, which were considered harmful. Interestingly, however, we found no overlap in positively selected genes among different cavefish lineages, indicating that the phenotypic convergence in cavefishes was not caused by positive selection of the same sets of genes. Analyses of previously identified candidate genes associated with cave phenotypes supported this conclusion. Genes belonging to the lipid metabolism functional ontology were under relaxed purifying selection in all cavefish genomes, which may be associated with the nutrient-poor habitat of cavefishes. Our work reveals previously uncharacterized patterns of cavefish genome evolution and provides comparative insights into the evolution of cave-associated phenotypic traits.

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.

Genomic Survey and Microsatellite Marker Investigation of Patagonian Moray Cod (Muraenolepis orangiensis)

Simple Summary

Patagonian moray cod is known to inhabit the cold waters near Antarctica, and it belongs to the Muraenolepis genus. This genus has seven species, and five of them are recently reported. The Muraenolepis genus has similar morphological characters, and this is a limitation of taxonomical classification. In this study, a genome survey and microsatellite marker analysis were conducted to characterize the genome profile for classification. As a result, genomic data such as genome size and microsatellite motifs were obtained.

Abstract

The Muraenolepididae family of fishes, known as eel cods, inhabits continental slopes and shelves in the Southern Hemisphere. This family belongs to the Gadiformes order, which constitutes one of the most important commercial fish resources worldwide, but the classification of the fish species in this order is ambiguous because it is only based on the morphological and habitat characteristics of the fishes. Here, the genome of Patagonian moray cod was sequenced using the Illumina HiSeq platform, and screened for microsatellite motifs. The genome was predicted to be 748.97 Mb, with a heterozygosity rate of 0.768%, via K-mer analysis (K = 25). The genome assembly showed that the total size of scaffolds was 711.92 Mb and the N50 scaffold length was 1522 bp. Additionally, 4,447,517 microsatellite motifs were identified from the genome survey assembly, and the most abundant motif type was found to be AC/GT. In summary, these data may facilitate the identification of molecular markers in Patagonian moray cod, which would be a good basis for further whole-genome sequencing with long read sequencing technology and chromosome conformation capture technology, as well as population genetics.

Evolution of thermal physiology alters the projected range of threespine stickleback under climate change

Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of its range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographical range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data were either static ("no evolution" models), allowed to evolve at observed evolutionary rates ("evolution" models) or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL; "scaled evolution" models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.

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.

De novo transcriptome assemblies of red king crab (Paralithodes camtschaticus) and snow crab (Chionoecetes opilio) molting gland and eyestalk ganglia - Temperature effects on expression of molting and growth regulatory genes in adult red king crab

Red king crab (Paralithodes camtschaticus) and snow crab (Chionoecetes opilio) are deep-sea crustaceans widely distributed in the North Pacific and Northwest Atlantic Oceans. These giant predators have invaded the Barents Sea over the past decades, and climate-driven temperature changes may influence their distribution and abundance in the sub-Arctic region. Molting and growth in crustaceans are strongly affected by temperature, but the underlying molecular mechanisms are little known, particularly in cold-water species. Here, we describe multiple regulatory factors in the two high-latitude crabs by developing de novo transcriptomes from the molting gland (Y-organ or YO) and eye stalk ganglia (ESG), in addition to the hepatopancreas and claw muscle of red king crab. The Halloween genes encoding the ecdysteroidogenic enzymes were expressed in YO, and the ESG contained multiple neuropeptides, including molt-inhibiting hormone (MIH), crustacean hyperglycemic hormone (CHH), and ion-transport peptide (ITP). Both crabs expressed a diversity of growth-related factors, such as mTOR, AKT, Rheb and AMPKα, and stress-responsive factors, including multiple heat shock proteins (HSPs). Temperature effects on the expression of key regulatory genes were quantified by qPCR in adult red king crab males kept at 4 °C or 10 °C for two weeks during intermolt. The Halloween genes tended to be upregulated in YO at high temperature, while the ecdysteroid receptor and several growth regulators showed tissue-specific responses to elevated temperature. Constitutive and heat-inducible HSPs were expressed in an inverse temperature-dependent manner, suggesting that adult red king crabs can acclimate to increased water temperatures.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

The HSP/co-chaperone network in environmental cold adaptation of Chilo suppressalis

Winter cold is one of the major environmental stresses for ectotherm species. Chilo suppressalis, a notorious lepidopteran pest of rice, has a wide geographic region that includes temperate zones with severe environmental conditions. Although C. suppressalis exhibits remarkable cold tolerance, its cold-adaptation mechanisms remain unclear. Here, we used bioinformatics approaches to evaluate transcript levels of genes comprising the C. suppressalis heat shock protein (HSP)/co-chaperone network in response to cold-induced stress. Using all such genes identified in the C. suppressalis genome, we experimentally examined the corresponding transcript levels under cold-acclimation or intermittent cold-shock stresses in diapause and non-diapausing larvae. In total, we identified 19 HSPs and 8 HSP co-chaperones in the C. suppressalis genome. Nine (hsp90, hsp75, hsp70, hsp40, small hsp, activator of 90 kDa heat shock protein ATPase-like, heat shock factor, heat shock factor binding protein 1-like and HSPB1-associated protein 1) were highly cold-inducible and likely comprise the principal cold-response HSP/co-chaperone network in C. suppressalis. We also found that transcriptional regulation of the HSP/co-chaperone networks response differs between cold-acclimation and short-term cold-shock. Moreover, activation of the HSP/co-chaperone network depends on the diapause state of overwintering larvae and cold acclimation may further increase larval cold tolerance. These results provide key new insights in the cold-adaptation mechanisms in C. suppressalis.

Genome survey and microsatellite motif identification of Pogonophryne albipinna

The genus Pogonophryne is a speciose group that includes 28 species inhabiting the coastal or deep waters of the Antarctic Southern Ocean. The genus has been divided into five species groups, among which the P. albipinna group is the most deep-living group and is characterized by a lack of spots on the top of the head. Here, we carried out genome survey sequencing of P. albipinna using the Illumina HiSeq platform to estimate the genomic characteristics and identify genome-wide microsatellite motifs. The genome size was predicted to be ∼883.8 Mb by K-mer analysis (K = 25), and the heterozygosity and repeat ratio were 0.289 and 39.03%, respectively. The genome sequences were assembled into 571624 contigs, covering a total length of ∼819.3 Mb with an N50 of 2867 bp. A total of 2217422 simple sequence repeat (SSR) motifs were identified from the assembly data, and the number of repeats decreased as the length and number of repeats increased. These data will provide a useful foundation for the development of new molecular markers for the P. albipinna group as well as for further whole-genome sequencing of P. albipinna.

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.

Fish lysozyme gene family evolution and divergent function in early development

Lysozymes are an ancient group of antimicrobial enzymes of the innate immune system. Here we provide a comparative analysis of the evolution and function of lysozymes during early development in fish, the most speciose vertebrate group. In fishes, lineage and species-specific evolution of both C-type (chicken or conventional) and G-type (goose type) genes occurred. Phylogenetic analysis revealed that the teleost lysozyme G-type members group with the tetrapod homologues but the teleost C-type form three different clusters with the tetrapods. Most of the teleost C-type cluster with tetrapod Lyz but there are some that group with the mammalian Lyzl1/2 and LALBA. This suggests that early in gnathostome evolution these genes already existed and that lyzl1/2 and lalba genes are present in fish and tetrapods. Gene synteny analysis to confirm sequence orthologies failed to identify conserved genome regions between teleosts and other vertebrates lysozyme gene regions suggesting that in the ancestral bony fish genome lyz, lyzl1/2, lalba and lyg precursor genes were transposed to different chromosome regions. The homologue of the mammalian lactalbumin (LALBA) gene was identified for the first time in teleosts and was expressed in skin and during egg and larval development. Lysozyme activity was detected in teleost eggs and varied between species and in the gilthead sea bream lyg and lalba transcript abundance differed in eggs and larvae from different brood stock suggesting differences exist in maternal innate immune protection.

Antarctic ecosystems in transition - life between stresses and opportunities

Important findings from the second decade of the 21st century on the impact of environmental change on biological processes in the Antarctic were synthesised by 26 international experts. Ten key messages emerged that have stakeholder-relevance and/or a high impact for the scientific community. They address (i) altered biogeochemical cycles, (ii) ocean acidification, (iii) climate change hotspots, (iv) unexpected dynamism in seabed-dwelling populations, (v) spatial range shifts, (vi) adaptation and thermal resilience, (vii) sea ice related biological fluctuations, (viii) pollution, (ix) endangered terrestrial endemism and (x) the discovery of unknown habitats. Most Antarctic biotas are exposed to multiple stresses and considered vulnerable to environmental change due to narrow tolerance ranges, rapid change, projected circumpolar impacts, low potential for timely genetic adaptation, and migration barriers. Important ecosystem functions, such as primary production and energy transfer between trophic levels, have already changed, and biodiversity patterns have shifted. A confidence assessment of the degree of 'scientific understanding' revealed an intermediate level for most of the more detailed sub-messages, indicating that process-oriented research has been successful in the past decade. Additional efforts are necessary, however, to achieve the level of robustness in scientific knowledge that is required to inform protection measures of the unique Antarctic terrestrial and marine ecosystems, and their contributions to global biodiversity and ecosystem services.

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.

Evolution of chaperome gene expression and regulatory elements in the antarctic notothenioid fishes

Confined within the cold-stable Southern Ocean, Antarctic notothenioid fishes have undergone an evolutionary loss of the inducible heat shock response (HSR), while facing perpetual low-temperature challenges to cellular proteostasis. This study examines how evolution in chronic cold has affected the shared cellular apparatus that mediates proteostasis under normal and heat stressed states. To deduce Antarctic-specific changes, we compared native expression levels across the full suite of chaperome genes and assessed the structural integrity of two crucial HSR regulators - Heat Shock Factor 1 (HSF1) that activates HSR, and heat shock elements (HSEs), the binding sites for HSF1 - between Antarctic fishes and the basal temperate notothenioid Eleginops maclovinus. Native expression levels of Antarctic fish chaperomes showed very modest changes overall, contrary to the common view of constitutive upregulation in the cold. Only a few cytosolic HSP70 genes showed greater transcription, with only the ancestrally-inducible HSPA6 strongly upregulated across all Antarctic species. Additionally, the constant cold has apparently not relaxed the selective pressures on maintaining HSF1 and HSEs in Antarctic fish. Instead, we found HSF1 experienced intensified selective pressure, with conserved sequence changes in Antarctic species suggesting optimization for non-heat-stress functional roles. HSEs of the HSP70 gene family have largely remained conserved in canonical sequence motifs and copy numbers as in E. maclovinus, showing limited impact of relaxed selective pressure. This study shows that evolution in chronic cold has led to both subtle and distinctive changes in the cellular apparatus for proteostasis and HSR, with functional consequences amenable to experimental evaluation.

Characterization of the hypoxia-inducible factor-1 pathway in hearts of Antarctic notothenioid fishes

The ability of Antarctic notothenioid fishes to mount a robust molecular response to hypoxia is largely unknown. The transcription factor, hypoxia-inducible factor-1 (HIF-1), a heterodimer of HIF-1α and HIF-1β subunits, is the master regulator of oxygen homeostasis in most metazoans. We sought to determine if, in the hearts of Antarctic notothenioids, HIF-1 is activated and functional in response to either an acute heat stress or hypoxia. The red-blooded Notothenia coriiceps and the hemoglobinless icefish, Chaenocephalus aceratus, were exposed to their critical thermal maximum (CT_MAX) or hypoxia (5.0 ± 0.3 mg of O₂ L^-1) for 2 h. Additionally, N. coriiceps was exposed to 2.3 ± 0.3 mg of O₂ L^-1 for 12 h, and red-blooded Gobionotothen gibberifrons was exposed to both levels of hypoxia. Levels of HIF-1α were quantified in nuclei isolated from heart ventricles using western blotting. Transcript levels of genes involved in anaerobic metabolism, and known to be regulated by HIF-1, were quantified by real-time PCR, and lactate levels were measured in heart ventricles. Protein levels of HIF-1α increase in nuclei of hearts of N. coriiceps and C. aceratus in response to exposure to CT_MAX and in hearts of N. coriiceps exposed to severe hypoxia, yet mRNA levels of anaerobic metabolic genes do not increase in any species, nor do lactate levels increase, suggesting that HIF-1 does not stimulate metabolic remodeling in hearts of notothenioids under these conditions. Together, these data suggest that Antarctic notothenioids may be vulnerable to hypoxic events, which are likely to increase with climate warming.

The complete mitochondrial genome of the Antarctic marbled rockcod, Notothenia rossii (Perciformes, Nototheniidae)

The complete mitochondrial genome of Notothenia rossii was obtained using PacBio Sequel long-read sequencing platform. The mitogenome of N. rossii was circular form and 18,274 bp long, which consists of 13 protein-coding genes, 24 tRNAs, 2 rRNAs, and non-coding control region. Particularly, we found duplicated tRNA^Thr and tRNA^Pro in addition to the typical 22 tRNAs. The phylogenetic tree revealed that N. rossii was most closely related to N. coriiceps among species in the Nototheniidae clade within the suborder Notothenioidei.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae)

Background

Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades.

Results

We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia.

Conclusions

The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.

LPS Modulates the Expression of Iron-Related Immune Genes in Two Antarctic Notothenoids

The non-specific immunity can induce iron deprivation as a defense mechanism against potential bacterial pathogens, but little information is available as to its role in Antarctic fish. In this study the response of iron metabolism related genes was evaluated in liver and head kidney of the Antarctic notothenoids Notothenia coriiceps and Notothenia rossii 7 days after lipopolysaccharide (LPS) injection. Average plasma Fe²⁺ concentration was unaffected by treatment in any of the species. The gene expression response to LPS varied between tissues and species, being stronger in N. coriiceps and more prominent in the head kidney than liver. The reaction to LPS was marked by increased individual variability in most genes analyzed, even when the change in expression was not statistically significant, suggesting different individual sensitivity and coping responses in these wild fish. We found that iron related genes had an attenuated and homogenous response to LPS but there was no detectable relationship between plasma Fe²⁺ and gene expression. However, overall in both tissues and species LPS exposure set a multilevel response that concur to promote intracellular accumulation of iron, an indication that Antarctic Notothenoids use innate nutritional immunity as a resistance mechanism against pathogens.

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.

The unique structural features of carbonmonoxy hemoglobin from the sub-Antarctic fish Eleginops maclovinus

Tetrameric hemoglobins (Hbs) are prototypical systems for the investigations of fundamental properties of proteins. Although the structure of these proteins has been known for nearly sixty years, there are many aspects related to their function/structure that are still obscure. Here, we report the crystal structure of a carbonmonoxy form of the Hb isolated from the sub-Antarctic notothenioid fish Eleginops maclovinus characterised by either rare or unique features. In particular, the distal site of the α chain results to be very unusual since the distal His is displaced from its canonical position. This displacement is coupled with a shortening of the highly conserved E helix and the formation of novel interactions at tertiary structure level. Interestingly, the quaternary structure is closer to the T-deoxy state of Hbs than to the R-state despite the full coordination of all chains. Notably, these peculiar structural features provide a rationale for some spectroscopic properties exhibited by the protein in solution. Finally, this unexpected structural plasticity of the heme distal side has been associated with specific sequence signatures of various Hbs.

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.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

The ancient origins and functional versatility of globins make them ideal subjects for studying physiological adaptation to environmental change. Our goals in this review are to describe the evolution of the vertebrate globin gene superfamily and to explore the structure/function relationships of hemoglobin, myoglobin, neuroglobin and cytoglobin in teleost fishes. We focus on the globins of Antarctic notothenioids, emphasizing their adaptive features as inferred from comparisons with human proteins. We dedicate this review to Guido di Prisco, our co-author, colleague, friend, and husband of C.V. Ever thoughtful, creative, and enthusiastic, Guido spearheaded study of the structure, function, and evolution of the hemoglobins of polar fishes - this review is testimony to his wide-ranging contributions. Throughout his career, Guido inspired younger scientists to embrace polar biological research, and he challenged researchers of all ages to explore evolutionary adaptation in the context of global climate change. Beyond his scientific contributions, we will miss his warmth, his culture, and his great intellect. Guido has left an outstanding legacy, one that will continue to inspire us and our research.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Phylogenomics investigation of sparids (Teleostei: Spariformes) using high-quality proteomes highlights the importance of taxon sampling

Sparidae (Teleostei: Spariformes) are a family of fish constituted by approximately 150 species with high popularity and commercial value, such as porgies and seabreams. Although the phylogeny of this family has been investigated multiple times, its position among other teleost groups remains ambiguous. Most studies have used a single or few genes to decipher the phylogenetic relationships of sparids. Here, we conducted a thorough phylogenomic analysis using five recently available Sparidae gene-sets and 26 high-quality, genome-predicted teleost proteomes. Our analysis suggested that Tetraodontiformes (puffer fish, sunfish) are the closest relatives to sparids than all other groups used. By analytically comparing this result to our own previous contradicting finding, we show that this discordance is not due to different orthology assignment algorithms; on the contrary, we prove that it is caused by the increased taxon sampling of the present study, outlining the great importance of this aspect in phylogenomic analyses in general.

The sequencing and de novo assembly of the Larimichthys crocea genome using PacBio and Hi-C technologies

Larimichthys crocea is an endemic marine fish in East Asia that belongs to Sciaenidae in Perciformes. L. crocea has now been recognized as an "iconic" marine fish species in China because not only is it a popular food fish in China, it is a representative victim of overfishing and still provides high value fish products supported by the modern large-scale mariculture industry. Here, we report a chromosome-level reference genome of L. crocea generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The genome sequences were assembled into 1,591 contigs with a total length of 723.86 Mb and a contig N50 length of 2.83 Mb. After chromosome-level scaffolding, 24 scaffolds were constructed with a total length of 668.67 Mb (92.48% of the total length). Genome annotation identified 23,657 protein-coding genes and 7262 ncRNAs. This highly accurate, chromosome-level reference genome of L. crocea provides an essential genome resource to support the development of genome-scale selective breeding and restocking strategies of L. crocea.

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.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

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.

Antarctic blackfin icefish genome reveals adaptations to extreme environments

Icefishes (suborder Notothenioidei; family Channichthyidae) are the only vertebrates that lack functional haemoglobin genes and red blood cells. Here, we report a high-quality genome assembly and linkage map for the Antarctic blackfin icefish Chaenocephalus aceratus, highlighting evolved genomic features for its unique physiology. Phylogenomic analysis revealed that Antarctic fish of the teleost suborder Notothenioidei, including icefishes, diverged from the stickleback lineage about 77 million years ago and subsequently evolved cold-adapted phenotypes as the Southern Ocean cooled to sub-zero temperatures. Our results show that genes involved in protection from ice damage, including genes encoding antifreeze glycoprotein and zona pellucida proteins, are highly expanded in the icefish genome. Furthermore, genes that encode enzymes that help to control cellular redox state, including members of the sod3 and nqo1 gene families, are expanded, probably as evolutionary adaptations to the relatively high concentration of oxygen dissolved in cold Antarctic waters. In contrast, some crucial regulators of circadian homeostasis (cry and per genes) are absent from the icefish genome, suggesting compromised control of biological rhythms in the polar light environment. The availability of the icefish genome sequence will accelerate our understanding of adaptation to extreme Antarctic environments.

Insertion Hot Spots of DIRS1 Retrotransposon and Chromosomal Diversifications among the Antarctic Teleosts Nototheniidae

By their faculty to transpose, transposable elements are known to play a key role in eukaryote genomes, impacting both their structuration and remodeling. Their integration in targeted sites may lead to recombination mechanisms involved in chromosomal rearrangements. The Antarctic fish family Nototheniidae went through several waves of species radiations. It is a suitable model to study transposable element (TE)-mediated mechanisms associated to genome and chromosomal diversifications. After the characterization of Gypsy (GyNoto), Copia (CoNoto), and DIRS1 (YNoto) retrotransposons in the genomes of Nototheniidae (diversity, distribution, conservation), we focused on their chromosome location with an emphasis on the three identified nototheniid radiations (the Trematomus, the plunderfishes, and the icefishes). The strong intrafamily TE conservation and wide distribution across species of the whole family suggest an ancestral acquisition with potential secondary losses in some lineages. GyNoto and CoNoto (including Hydra and GalEa clades) mostly produced interspersed signals along chromosomal arms. On the contrary, insertion hot spots accumulating in localized regions (mainly next to centromeric and pericentromeric regions) highlighted the potential role of YNoto in chromosomal diversifications as facilitator of the fusions which occurred in many nototheniid lineages, but not of the fissions.

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.

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.

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.

De Novo Gene Evolution of Antifreeze Glycoproteins in Codfishes Revealed by Whole Genome Sequence Data

New genes can arise through duplication of a pre-existing gene or de novo from non-coding DNA, providing raw material for evolution of new functions in response to a changing environment. A prime example is the independent evolution of antifreeze glycoprotein genes (afgps) in the Arctic codfishes and Antarctic notothenioids to prevent freezing. However, the highly repetitive nature of these genes complicates studies of their organization. In notothenioids, afgps evolved from an extant gene, yet the evolutionary origin of afgps in codfishes is unknown. Here, we demonstrate that afgps in codfishes have evolved de novo from non-coding DNA 13-18 Ma, coinciding with the cooling of the Northern Hemisphere. Using whole-genome sequence data from several codfishes and notothenioids, we find higher copy number of afgp in species exposed to more severe freezing suggesting a gene dosage effect. Notably, antifreeze function is lost in one lineage of codfishes analogous to the afgp losses in non-Antarctic notothenioids. This indicates that selection can eliminate the antifreeze function when freezing is no longer imminent. In addition, we show that evolution of afgp-assisting antifreeze potentiating protein genes (afpps) in notothenioids coincides with origin and lineage-specific losses of afgp. The origin of afgps in codfishes is one of the first examples of an essential gene born from non-coding DNA in a non-model species. Our study underlines the power of comparative genomics to uncover past molecular signatures of genome evolution, and further highlights the impact of de novo gene origin in response to a changing selection regime.

Comparative genome analysis of 52 fish species suggests differential associations of repetitive elements with their living aquatic environments

BACKGROUND

Repetitive elements make up significant proportions of genomes. However, their roles in evolution remain largely unknown. To provide insights into the roles of repetitive elements in fish genomes, we conducted a comparative analysis of repetitive elements of 52 fish species in 22 orders in relation to their living aquatic environments.

RESULTS

The proportions of repetitive elements in various genomes were found to be positively correlated with genome sizes, with a few exceptions. More importantly, there appeared to be specific enrichment between some repetitive element categories with species habitat. Specifically, class II transposons appear to be more abundant in freshwater bony fish than in marine bony fish when phylogenetic relationship is not considered. In contrast, marine bony fish harbor more tandem repeats than freshwater species. In addition, class I transposons appear to be more abundant in primitive species such as cartilaginous fish and lamprey than in bony fish.

CONCLUSIONS

The enriched association of specific categories of repetitive elements with fish habitats suggests the importance of repetitive elements in genome evolution and their potential roles in fish adaptation to their living environments. However, due to the restriction of the limited sequenced species, further analysis needs to be done to alleviate the phylogenetic biases.

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.

Structural analyses of adenylate kinases from Antarctic and tropical fishes for understanding cold adaptation of enzymes

Psychrophiles are extremophilic organisms capable of thriving in cold environments. Proteins from these cold-adapted organisms can remain physiologically functional at low temperatures, but are structurally unstable even at moderate temperatures. Here, we report the crystal structure of adenylate kinase (AK) from the Antarctic fish Notothenia coriiceps, and identify the structural basis of cold adaptation by comparison with homologues from tropical fishes including Danio rerio. The structure of N. coriiceps AK (AKNc) revealed suboptimal hydrophobic packing around three Val residues in its central CORE domain, which are replaced with Ile residues in D. rerio AK (AKDr). The Val-to-Ile mutations that improve hydrophobic CORE packing in AKNc increased stability at high temperatures but decreased activity at low temperatures, suggesting that the suboptimal hydrophobic CORE packing is important for cold adaptation. Such linkage between stability and activity was also observed in AKDr. Ile-to-Val mutations that destabilized the tropical AK resulted in increased activity at low temperatures. Our results provide the structural basis of cold adaptation of a psychrophilic enzyme from a multicellular, eukaryotic organism, and highlight the similarities and differences in the structural adjustment of vertebrate and bacterial psychrophilic AKs during cold adaptation.

Draft genome of the Antarctic dragonfish, Parachaenichthys charcoti

The Antarctic bathydraconid dragonfish, Parachaenichthys charcoti, is an Antarctic notothenioid teleost endemic to the Southern Ocean. The Southern Ocean has cooled to -1.8ºC over the past 30 million years, and the seawater had retained this cold temperature and isolated oceanic environment because of the Antarctic Circumpolar Current. Notothenioids dominate Antarctic fish, making up 90% of the biomass, and all notothenioids have undergone molecular and ecological diversification to survive in this cold environment. Therefore, they are considered an attractive Antarctic fish model for evolutionary and ancestral genomic studies. Bathydraconidae is a speciose family of the Notothenioidei, the dominant taxonomic component of Antarctic teleosts. To understand the process of evolution of Antarctic fish, we select a typical Antarctic bathydraconid dragonfish, P. charcoti. Here, we have sequenced, de novo assembled, and annotated a comprehensive genome from P. charcoti. The draft genome of P. charcoti is 709 Mb in size. The N50 contig length is 6145 bp, and its N50 scaffold length 178 362 kb. The genome of P. charcoti is predicted to contain 32 712 genes, 18 455 of which have been assigned preliminary functions. A total of 8951 orthologous groups common to 7 species of fish were identified, while 333 genes were identified in P. charcoti only; 2519 orthologous groups were also identified in both P. charcoti and N. coriiceps, another Antarctic fish. Four gene ontology terms were statistically overrepresented among the 333 genes unique to P. charcoti, according to gene ontology enrichment analysis. The draft P. charcoti genome will broaden our understanding of the evolution of Antarctic fish in their extreme environment. It will provide a basis for further investigating the unusual characteristics of Antarctic fishes.