Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts

Abstract

Polar fishes have evolved antifreeze proteins (AFPs) that allow them to survive in subzero temperatures. We performed deep transcriptomic sequencing on a postlarval/juvenile variegated snailfish, Liparis gibbus (Actinopterygii: Scorpaeniformes: Cottoidei: Liparidae), living in an iceberg habitat (−2°C) in Eastern Greenland and report detection of highly expressed transcripts that code for putative AFPs from 2 gene families, Type I and LS-12-like proteins (putative Type IV AFPs). The transcripts encoding both proteins have expression levels among the top <1% of expressed genes in the fish. The Type I AFP sequence is different from a reported Type I AFP from the same species, possibly expressed from a different genetic locus. While prior findings from related adult sculpins suggest that LS-12-like/Type IV AFPs may not have a role in antifreeze protection, our finding of very high relative gene expression of the LS-12-like gene suggests that highly active transcription of the gene is important to the fish in the iceberg habitat and raises the possibility that weak or combinatorial antifreeze activity could be beneficial. These findings highlight the physiological importance of antifreeze proteins to the survival of fishes living in polar habitats.

Introduction

Fish antifreeze proteins

The icy waters of polar oceans are an extreme environment for marine life that limits inhabitants to those with mechanisms to cope with freezing (−1.9°C) temperatures. The body fluids of fishes are hypoosmotic to seawater, and hence, have a higher freezing point than seawater.^1,2 In addition, fishes are incapable of surviving even partial freezing of their body fluids, unlike some species of reptiles and insects.^3,4 Antifreeze proteins (AFPs) were first discovered in Antarctic notothenioid fishes almost 5 decades ago ⁵ and many polar and subpolar fishes are known to possess ice-binding AFPs that allow them to successfully inhabit icy habitats. ⁶ Proteins with ice-binding properties have evolved independently in many different organisms, spanning bacteria, lichens, plants, insects, and vertebrates, to enable them to survive in sub-zero environments.^7,8 AFPs offer an advantage to animals that dwell in places below the freezing point by preventing the formation of large ice grains inside cells and body fluids, thus preserving body fluids in a liquid state. AFPs limit the growth of ice crystals to manageable sizes, rather than fully inhibiting the growth of ice crystals inside fishes. Thus, they are also referred to as ice-structuring proteins. ⁹ However, for the purpose of this article, we use the term AFP.

Fish AFPs, that are primarily synthesized in the liver and exported to the blood, can be seasonally regulated (eg, by water temperature or photoperiod) or maintained at high concentrations year-round. ¹⁰ Five classes of AFPs (Types I, II, III, and IV, as well as antifreeze glycoproteins (AFGPs)), that are differentiated based on their origin, amino acid composition, and the resulting 3D structure of the polypeptides, have been identified and characterized from fishes living in mid- to high-latitude seawater.^11,12 According to Graham et al, ¹³ Type I AFPs, that consist of single alanine-rich alpha-helices, are found in cunner, winter flounder, sculpins and snailfishes. Type II AFPs, derived from lectins, are found in smelt, herring, and sea raven. More specifically, Type II AFPs share high sequence similarity with the carbohydrate recognition domain of Ca²⁺-dependent type (C-type) lectins. ¹⁴ However, Nishimiya et al ¹⁴ identified the AFP of the longsnout poacher, Brachyopsis segaliensis, as a Ca²⁺-independent Type II AFP. Graham et al ¹⁵ hypothesized that herring acquired their Type II AFP via horizontal gene transfer of foreign extracellular DNA into herring sperm during external fertilization. The presence of a multitude of transposable elements (TEs) within the herring genome ¹⁶ allowed Graham and Davies ¹⁷ to demonstrate that the Type II AFP was transferred from herring to smelt because the smelt AFP is bordered by 3 herring-specific TEs. Type III AFPs, derived from the C-terminal domain of the enzyme sialic acid synthase, ¹⁸ have been found in eelpouts and wolffish. LS-12 proteins, exhibiting weak antifreeze activity, were described as a new class of AFP, Type IV AFPs, and were first isolated and characterized in the longhorn sculpin, Myoxocephalus octodecimspinosis. ¹⁹ AFGPs have been found in Arctic codfishes (evolved from non-coding DNA) ²⁰ and Antarctic notothenioids (derived from a trypsinogen-like protease).^21,22

According to Deng et al, ¹⁹ Cheng, ²³ and Cheng and Zhuang, ²⁴ the evolutionary precursor of LS-12 proteins is a blood plasma apolipoprotein. The subsequent discovery of homologs of the LS-12 genes across a broad range of tropical and subtropical teleost fishes called into question whether the LS-12 protein, and therefore the new class of Type IV AFPs, are actually apolipoproteins that coincidentally have negligible antifreeze properties in polar species, or whether they have different functions altogether (eg, gastrulation; lipid metabolism).^25-28 Li et al ²⁹ analyzed the expression of the LS-12-like gene homolog in the gibel carp (Carassius auratus gibelio) and found that the gene is predominantly expressed during embryogenesis; and that the protein was no longer detectable in larvae older than 13 days. Deng and Laursen ³⁰ showed that the concentration of the LS-12 gene in the longhorn sculpin, M. octodecimspinosis, is so low that without concentration of the protein, antifreeze activity could not be detected. A protein homologlous to the M. octodecimspinosis LS-12/TypeIV AFP from the subtropical olive flounder (Paralichthys olivaceus), ³¹ which does not require AFPs for survival, has been cloned, heterologously expressed (in E. coli) and assayed for antifreeze activity. Based on its comparatively “very weak” thermal hysteresis activity (0.07 ± 0.01°C at 0.5 mg/mL), the P. olivaceus LS-12 homolog would not prevent freezing at subzero temperatures. Additionally, the thermal hysteresis activity of the LS-12 homolog in the subtropical P. olivaceus is comparable to that of LS-12 proteins in cold-water fishes. The authors concluded that the activity, inheritance pattern, and sequence conservation of the LS-12 gene family suggests that its antifreeze activity was not gained due to selective pressure imparted by a cold habitat, but is rather a chance activity of a protein needed for other purposes.

Fish often have more than one copy of AFP genes in their genetic repertoire. The shorthorn sculpin Myoxocephalus scorpius and longhorn sculpin Myoxocephalus octodecemspinosus have both Type I AFPs and LS-12 gene homologs. In those fishes, the Type I AFP that is present in peripheral tissues (skin) was hypothesized to provide sufficient protection from freezing,^25,32 while their LS-12-like protein was found in insufficient concentrations in fish tissues to provide freezing protection. ²⁵ In fishes that have either Type I, II or III AFPs, the fish often carry multiple copies of genes encoding proteins from the same AFP family in their genomes. ³³ The expanded gene families may account for high expression of AFPs, while LS-12 gene homologs are hypothesized to be less highly expressed because they are found in a single locus. ²⁵ Yet, in this study we report exceptionally high expression of transcripts coding for Type I AFPs and an LS-12-like gene homolog in the same juvenile individual of the variegated snailfish, Liparis gibbus.

Classification and general biology of snailfishes (Liparidae)

Liparidae is a diverse family of cottoid fishes with approximately 420 species arrayed in 32 genera. Liparids are characterized by their asquamate, elongate bodies, and by a sucking disk formed by modified pelvic fins. Liparidae is bipolar in distribution, with members of the widespread genus Liparis (~60 spp.) restricted to the northern hemisphere. Liparids inhabit the intertidal zone to the extreme depths of the Mariana Trench, and exhibit a greater depth range than any other family of fishes. ³⁴ They are among the most common hadal vertebrate species, and Pseudoliparis swirei and the undescribed “ethereal snailfish” are currently recognized as the deepest dwelling vertebrates (8076 and 8178 m, respectively).^34,35

Liparis gibbus is a polar benthic species occurring in the Arctic, North Atlantic, and the Northeast Pacific oceans (Greenland, Canada, Russia, Svalbard, and southeastern Alaska). ³⁶ It generally occurs between 100 and 200 m depth, but is frequently encountered at much shallower depths, has an average length of 11 cm (max: 52 cm), and typically feeds on amphipods and crabs.^37-39 L. gibbus is currently the only polar fish reported to have biofluorescence. ⁴⁰ Prior research of Liparis atlanticus and L. gibbus AFPs discovered 1 Type I AFP in L. atlanticus and 2 Type I AFP proteins with confirmed antifreeze activity in blood plasma and an identical Type I AFP in the skin of L. gibbus. Additionally, Southern blot analysis of L. gibbus genomic DNA indicated that it codes for multiple copies of Type I AFP coding sequences. ⁴¹

Results

Transcriptome statistics

In this study, we sequenced and assembled the transcriptome of a post-larval juvenile variegated snailfish, L. gibbus (Figure 1), collected on an iceberg in Eastern Greenland (Figure 2). The data comes from the whole animal minus the gut. After filtering based on read depth (retaining transcripts with >100 reads mapped), the transcriptome assembly consisted of 93 433 contigs with 68 031 of those with a length greater than 1000 bps. The average transcript length was 2190 bp. BUSCO completeness of the filtered assembly shows that it contained complete sequences of 84.4% of the BUSCOS from the vertebrate set, and is only missing 7.4% of them, suggesting the assembly captured the majority of expected conserved genes in vertebrates (Supplemental Table S1).

Figure 1.

Post-larval/juvenile Liparis gibbus: (A) dorsal view and (B) ventral view. (AMNH 277096, 18 mm Standard Length (SL), measured from the tip of the snout (upper jaw) to the hypural flexure (at base of the caudal fin)).

Figure 2.

Study site showing the iceberg habitat in Greenland where Liparis gibbus was collected. The dive boat, a zodiac inflatable craft, can been seen at the bottom left, and the divers are visible near the center of the image. The image was shot using a small drone.

Liparis gibbus AFPs

The transcriptome was screened for putative AFPs from all 5 AFP families. We detected AFPs from Type I and Type IV AFP (LS-12-like) families, but not Type II, Type III, or AFGPs. We also investigated Type I AFPs and LS-12 like proteins in available whole genome assemblies of 2 fish from the Liparidae, the hadal zone snailfish Pseudoliparis swirei and Tanaka’s snailfish Liparis tanakae.

In the juvenile L. gibbus transcriptome assembly we identified 2 distinct full length, and 3 truncated, transcripts encoding putative Type I AFPs. The truncation in assembled contigs may be an artifact of a strong predicted secondary structure found in the 3′ end of the full-length transcripts (Supplemental Figure 1). In one of 2 available whole genome assemblies for L. tanakae (NCBI BioProject# PRJNA523297; https://doi.org/10.6084/m9.figshare.7819418.v1), we found 3 distinct copies of genes encoding full length Type I AFP proteins. In the other independent assembly of the L. tanakae genome, ³⁴ the regions containing the Type I AFP were not found. Likewise, the P. swirei genome assembly ³⁴ did not contain evidence of Type I AFP genes, however we cannot say whether that is because they are wholly absent or whether the Type I AFP genes are recalcitrant to assembly, as seems to be the case in L. tanakae.

Protein alignments of translations of the 5 complete or partial Type I AFPs discovered in our data with previously reported L. gibbus and L. atlanticus sequences, as well as the Type I AFP sequences discovered in the L. tanakae geneome, reveals that our whole transcriptome assembly uncovered the sequences of previously unreported L. gibbus Type I AFPs. The protein alignments show that the 2 novel full length AFPs share 66% and 62% identity, respectively, with AY455863 at the amino acid level (Figure 3). We did not detect the previously reported L. gibbus Type I AFP, AY455863, in our de novo transcriptome assembly. Additionally, the 3 Type I AFP sequences found in the L. tanakae genome vary from one another and are different from other reported Type I AFP sequences from fishes in the genus Liparis (Figure 3).

Figure 3.

Protein alignments between observed juvenile L. gibbus Type I AFP and LS-12-like proteins and homologous proteins from the same or related species. (A) Type I AFPs. (B) Type IV AFPs. Colors correspond to amino acid identity where light and dark green indicates primarily hydrophobic residues, red indicates positively charged residues, light blue indicates polar-alcohol residues, purple indicates polar-carboxamide residues, and dark blue indicates negatively charged residues. No highlighting indicates the indicated residue is different from the reference sequence which is the top sequence in each alignment. The column “cov” indicates percent coverage of the reference sequence and “pid” indicates overall percent identity to the reference sequence.

We identified 2 full length transcripts encoding LS-12-like proteins (putative Type IV AFPs) in the juvenile L. gibbus transcriptome assembly. Alignment of the 2 juvenile L. gibbus LS-12-like transcripts to the closely related L. tanakae genome suggests that they are splice variants. The longer contig from our juvenile L. gibbus assembly, contig 445 731, retains the first intron of the gene, producing an additional ATG start codon that adds a predicted secretory signal peptide on that isoform (Supplemental Figure 2). The shorter isoform, contig 445 757, has that intron spliced out and does not have any predicted signal peptides. Protein alignments showed that the LS-12-like protein from L. gibbus is nearly identical to the LS-12-like protein found in L. tanakae (Figure 3, 98.4% similarity). Additionally, unlike Type I AFPs, LS-12 proteins appear to occur at a single locus and a phylogenetic tree based on LS-12-like protein sequences suggests that the gene follows a pattern of divergence in accordance with speciation in the group (Figure 4).

Figure 4.

Phylogenetic tree depicting evolutionary relationships of the L. gibbus LS-12-like protein. The maximum likelihood tree was created in IQ-TREE using homologs of the longhorn sculpin, M. octodecemspinous, LS-12 protein (blue highlight). Support values are from 1000 ultrarapid bootstraps. L. gibbus sequences are highlighted in pink.

We performed 3D homology modeling using Alpha-Fold ⁴² to infer the tertiary structure of the putative antifreeze proteins. The Type I AFP exhibits a primarily alpha-helical structure, consistent with prior examinations of Type I AFP structure ⁴³ (Figure 5). The LS-12-like protein also folded into a primarily alpha helical structure, consistent with predicted structures of LS-12 proteins from fish in the literature. ⁴⁴

Figure 5.

Structural models of the Liparis gibbus Type I AFP and LS-12-like protein. Models were built using Alphafold and rendered in ChimeraX. (A) Putative Type I AFP, consisting of 2 alpha helices joined by an unstructured region. (B) Putative LS-12-like protein, consisting of 3 alpha helices joined by unstructured regions. The color scale is based on the predicted Local Distance Difference Test (pLDDT) score for each residue. A yellow-blue color, as seen for the bulk of the 2 predictions, indicates a pLDDT > 70% which can be interpreted as confidence in an accurate representation of the backbone structure of the protein model. Position of individual atoms (not shown) are not predicted with confidence, and in the case of these 2 proteins likely depend on interactions with their respective substrates of ice crystals and potentially ice crystals and lipids.

The putative full-length Type I AFP and LS-12-like transcript sequences of the snailfish L. gibbus were deposited under GenBank Accessions numbers MT678484 and MT678485, respectively. All raw sequence reads were deposited in NCBI’s publicly available Short Read Archive under BioProject ID PRJNA624035.

Relative expression levels

We found that one of the full-length Type I AFP transcripts and the shorter contig of the LS-12-like gene transcripts were among the top <1% most highly expressed contigs in our transcriptome (Figure 6). Two of the truncated Type I AFP transcripts were also highly expressed (Supplemental Figure 3).

Figure 6.

Expression plot indicating very high relative expression of a Type I AFP and the LS-12-like protein in the collected juvenile L. gibbus specimen. The Type I AFP and LS-12-like protein are among the top <1% most highly expressed genes in the L. gibbus transcriptome. Expression level, shown as the transcripts per million (TPM), the expected number of molecules of each transcript per million total transcripts, is plotted on a log₁₀ scale. The width of the violin plot is representative of the number of genes at a given expression level.

Discussion

The whole animal transcriptome from a juvenile variegated snailfish (L. gibbus) collected in Eastern Greenland revealed transcripts encoding expressed Type I AFPs and an LS-12-like protein that were among the most highly expressed transcripts in the animal (Figure 6). The Type I AFP sequences reported here are different from a previously reported Type I AFP from the same species (Figure 3), and are therefore likely expressed from a different genetic locus.

In an L. tanakae genome assembly (https://doi.org/10.6084/m9.figshare.7819418.v1), we found evidence of 3 genetic loci encoding full length Type I AFP genes and a single locus encoding an LS-12-like protein. The gene structure of the 2 types of genes is telling regarding their mode of inheritance and evolution. In the L. tanakae genome, all 3 detected Type I AFP genes lack introns and contain multiple non-synonymous differences indicative of selection on the genes. The multiple loci of intron-less genes are indicative of recent gene duplications, possibly by retrotransposition. A BLAST search of the leading genomic sequence, 5′ to the start codon of the Type I AFP coding sequence, reveals that the leading sequence is found in multiple copies in other fish genomes (Supplemental Figure 4), but not followed by an AFP, further evidence that the Type I AFP is associated with a mobile genetic element. Active gene duplication is often correlated with gene loss of one or more copies of the gene ⁴⁵ leading to multiple hypotheses regarding why we did not recover the exact previously reported L. gibbus Type I AFP sequence in our transcriptome. The gene could have been lost in the individual that we sequenced, or it could be repressed at this particular developmental stage, or it could be unassembled due to difficulties such as DNA repeats and RNA secondary structure.

The structure of the multiple copies of the Type I AFP gene stand in contrast to the genomic structure of the LS-12-like gene in L. tanakae. The LS-12-like gene is found at a single locus and contains introns. The intron/exon structure is also well conserved with related fish, like the hadal snailfish P. swirei, where the conserved exons are the same length with the same breakpoint in the 2 fish, but the introns, which are presumably under relaxed selection, vary in length and sequence (Figure 7). While some authors have hypothesized that the rapid expansion and multiple copies of Type I AFP genes in fish genomes explain high expression levels, it is worth noting that the LS-12 gene is present at a single, mature locus and is still among the <1% of expressed genes in the juvenile L. gibbus transcriptome. Transcriptional (eg, high transcription rate) and possibly post-transcriptional (eg, low degradation rate) regulation account for high expression levels in the absence of multiple gene copies.

Figure 7.

Genomic structure of Type I AFPs and LS-12 like genes. Three copies of a Type I AFP gene were found in 1 of 2 available L. tanakae whole genome assemblies. The genes were on separate contigs and all had the same intronless structure indicative of recent gene duplications, possibly by retrotransposition. In contrast, the LS-12-like gene contains introns, has a conserved structure in 2 fish genomes, and is single copy, indicative of a stable, conserved gene product.

High expression of a LS-12-like mRNA in a juvenile L. gibbus provides evidence that the protein product might be present at a high concentration. While a prior study concluded that LS-12-like Type IV AFPs do not exhibit strong enough antifreeze activity to be physiologically relevant, we argue that it is worth considering whether highly expressed Type I AFPs with strong antifreeze properties and LS-12-like proteins with weak antifreeze activity act synergistically to lower the freezing point of fish blood and tissues. To our knowledge a controlled experiment with a mixture of AFPs from different families (or putative families) has not been done. An alternative explanation is that LS-12 proteins do not contribute to antifreeze activity and are highly expressed for a reason other than promoting freeze tolerance. For example, in a freshwater trout that experiences seasonal variation in temperature, the proportion of unsaturated lipids carried by high-density lipoproteins (HDLs) increases with a decrease in temperature. ⁴⁶ Since LS-12-like proteins are structural analogs of Apo-A1 lipoproteins ⁴⁴ and Apo-A1 lipoproteins are the major component of HDLs, ⁴⁷ the very high expression of the L. gibbus LS-12-like gene on an iceberg may be related to a temperature dependent transport of unsaturated lipids rather than, or in addition to the reported weak antifreeze activity of the protein.

Since the mid-20th century, temperatures have increased twice as fast in the Arctic as in mid-latitudes and it is predicted that if Arctic sea ice decline continues at this current rate, it will be ice-free by the next 3 decades.^48,49 Arctic seas are inhabited by only ~270 fish species,^50,51 with 15 endemic marine fish species, mostly sculpins, snailfishes and eelpouts. ⁵² We hypothesize that with increasingly warming oceanic temperatures in the Arctic, ⁵² ice-dwelling specialists such as L. gibbus and the few other teleosts that use substantial metabolic resources to generate AFPs and regulate lipid composition may encounter increased competition by more temperate species that were previously unable to survive in sub-zero seawater temperatures. ⁵³

Materials and Methods

Collection

The juvenile Liparis gibbus specimen (18 mm SL) was collected on August 26, 2019 in an iceberg habitat, near Qaattu (65.788, −37.879) in Eastern Greenland, at 15 m depth and was preserved in RNA Stabilization Solution (RNAlater, Qiagen) for transport back to the laboratory (Figures 1 and 2). The specimen was collected via SCUBA using a hand-net, via the application of quinaldine, an anesthetic used to temporarily immobilize fishes for capture, to a targeted habitat. This study was approved and carried out in strict accordance with the recommendations in the Guidelines for the Use of Fishes in Research of the American Fisheries Society and the American Museum of Natural History’s Institutional Animal Care and Use Committee (IACUC).

Analysis

Given the small size of the snailfish, total RNA was extracted from the whole body (including muscle and skin) except the gut, which was dissected out to limit contamination with DNA from food. Illumina NovaSeq 2 × 150 bp sequencing and initial bioinformatic analyses were performed by Genewiz (South Plainfield, NJ). The NGS run yielded 281 412 170 reads (84 424 megabases) with a mean Q score of 35.71 (92.79% of bases ⩾Q30). Trimmomatic v0.36 ⁵⁴ was used to trim adapters and quality trim the raw data. Trinity v2.5 ⁵⁵ was used to assemble the transcriptome de-novo with a minimum contig length of 200 bp. SeqKit v0.11.0 ⁵⁶ was used to generate statistics for the transcriptome (Supplemental Table S1). Transcript quantification was accomplished by alignment and quantification of sequencing reads against the transcriptome assembly using Salmon (v1.2.0). ⁵⁷ For plotting the transcript abundance distribution, contigs that recruited low numbers of reads (⩽100 reads) were removed. The threshold was determined empirically by plotting the number of contigs remaining after removing contigs below threshold. The curve flattened out at 100 reads mapped (Supplemental Figure 5). Relative abundance of transcriptome contigs was expressed using transcripts per million (TPM). The abundance distribution was visualized in R (v3.6.3) ⁵⁸ using the ggplot2 package. ⁵⁹

The EMBOSS tool getorf was used to find open reading frames. The transcriptome was annotated using Diamond BLASTX. Any transcripts that were identified by Diamond BLASTX as AFPs were reciprocally blasted to nr and UniProtKB/Swiss-Prot using BLASTP. The results of the reciprocal blasts are presented. To discover additional putative AFPs, representative sequences of Type I (NCBI nucleotide: AY455862.2), Type II (UniProtKB: P84493), Type III (UniProtKB: P19614), Type IV (UniProtKB: P80961) AFPs, and AGFP (UniProtKB: P24856) were searched against the transcriptome assembly using BLASTN (nucleotide query, nucleotide database) and tBLASTN (protein query, nucleotide database) searches. The dual search strategy was needed because Type I AFPs contain repetitive amino acid motifs of 2 to 4 alanine residues in a row repeated throughout the length of the peptide sequence (65 out of 115 residues are alanine). The repetitive nature of the protein sequence leads to poor BLAST search scores. Due to the degenerate nature of the protein code, however, the nucleotide sequence that codes for the Type I AFP is not as repetitive and specific, and significant hits can be identified where present. In addition to the general BLAST search using known AFPs as a query, the transcripts in the top 1% of expression levels were searched against nr using BLASTX (nucleotide query, protein database) to discover hits to highly expressed transcripts containing c-type lectin or apolipoprotein domains, from which AFPs have evolved in other species.

Alignments to known AFP sequences were completed using MUSCLE v3.8 ⁶⁰ and visualized with MView v1.63 ⁶¹ accessed via EBI Web Services. ⁶²

Phylogenetic trees

Homologs of select L. gibbus transcripts and relevant antifreeze proteins were collected by BLASTP searches of potential L. gibbus AFPs against the NCBI database nr ⁶³ for Type IV and Type II AFPs. The top 50 BLAST hits were recovered for phylogenetic analysis. For an outgroup, the seed proteins were searched against nr limited to Chondrichthyes sequences (cartilaginous fishes, NCBI taxon ID 7777). Canonical antifreeze protein sequences were included in the BLAST hits by searching nr using antifreeze proteins as a seed: Uniprot protein P80961 was used as a seed for LS-12-like/Type IV AFPs. For Type I AFPs, nucleotide searches using BLAST were conducted using AY455863 as a seed, and Type I AFP sequences were collected manually from the NCBI nucleotide collection by searching with the terms: “type I” antifreeze protein.

Sequences were aligned using MAFFT v7.310, ⁶⁴ the alignments were trimmed with trimAl v1.2rev59, ⁶⁵ ModelTest-NG v0.1.7 ⁶⁶ was used to find best-fit models for each alignment (Type IV: LG), and the trimmed alignments were used to build a maximum likelihood phylogenetic tree in IQ-TREE 2 ⁶⁷ using 1000 ultrafast bootstraps. Sequence alignments are available in respective Supplemental Files (typeIV_apoCI.trimal).

Molecular modeling

An Alphafold2 colab notebook ⁶⁸ was used to generate the protein models of all putative AFPs from Liparis gibbus. ChimeraX ⁶⁹ was used to render models.

Supplemental Material

sj-trimal-1-evb-10.1177_11769343221118347 – Supplemental material for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts

Supplemental material, sj-trimal-1-evb-10.1177_11769343221118347 for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts by John A Burns, David F Gruber, Jean P Gaffney, John S Sparks and Mercer R Brugler in Evolutionary Bioinformatics

Supplemental Material

sj-trimal-2-evb-10.1177_11769343221118347 – Supplemental material for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts

Supplemental material, sj-trimal-2-evb-10.1177_11769343221118347 for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts by John A Burns, David F Gruber, Jean P Gaffney, John S Sparks and Mercer R Brugler in Evolutionary Bioinformatics

Supplemental Material

sj-trimal-3-evb-10.1177_11769343221118347 – Supplemental material for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts

Supplemental material, sj-trimal-3-evb-10.1177_11769343221118347 for Transcriptomics of a Greenlandic Snailfish Reveals Exceptionally High Expression of Antifreeze Protein Transcripts by John A Burns, David F Gruber, Jean P Gaffney, John S Sparks and Mercer R Brugler in Evolutionary Bioinformatics