Whole-genome duplication and the functional diversification of teleost fish hemoglobins

Hemoglobin wonders: a fascinating gas transporter dive into molluscs

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

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.

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Functional exploration of SNP mutations in HIF2αb gene correlated with hypoxia tolerance in blunt snout bream (Megalobrama amblycephala)

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia environment. Hypoxia-inducible factor (HIF) is the most critical factor in the HIF pathway, which strictly regulates the hypoxia stress process of fish. In this study, we found six hifα genes in blunt snout bream that demonstrated different expressions under hypoxia conditions. In HEK293T cells, all six hifαs were detected to activate the HRE region by luciferase reporter assay. More importantly, we identified two linkage-disequilibrium SNP sites at exon 203 and 752 of the hif2αb gene in blunt snout bream. Haplotype II (A²⁰³A⁷⁵²) and its homozygous diplotype II (A²⁰³A²⁰³A⁷⁵²A⁷⁵²) appeared frequently in a selected strain of blunt snout bream with hypoxia tolerance. Diplotype II has a lower oxygen tension threshold for loss of equilibrium (LOE_crit) over a similar range of temperatures. Moreover, its erythrocyte number increased significantly (p < 0.05) than those in diplotype I and diplotype III strains at 48 h of hypoxia. The enzymes related with hypoxia tolerant traits, i.e., reduced glutathione, superoxide dismutase, and catalase, were also significantly (p < 0.05) induced in diplotype II than in diplotype I or III. In addition, the expression of epo in the liver of diplotype II was significantly (p < 0.01) higher than that in the diplotype I or III strains at 48 h of hypoxia. Taken together, our results found that the hypoxia-tolerant-related diplotype II of hif2αb has the potential to be used as a molecular marker in future genetic breeding of hypoxia-tolerant strain.

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.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as “MNˮ and “LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

Hb adaptation to hypoxia in high-altitude fishes: Fresh evidence from schizothoracinae fishes in the Qinghai-Tibetan Plateau

Uncovering the genetic basis of hypoxic adaptation is one of the most active research areas in evolutionary biology. Among air-breathing vertebrates, modifications of hemoglobin (Hb) play a pivotal role in mediating an adaptive response to high-altitude hypoxia. However, the relative contributions in water-breathing organisms are still unclear. Here, we tested the Hb concentration of fish at different altitudes. All species showed species-specific Hb concentration, which has a non-positive correlation with altitude. Moreover, we investigated the expression of Hb genes by the RNA-seq and quantitative real-time PCR (qRT-PCR), and Hb composition by two-dimensional electrophoresis (2-DE). The results showed that the multiple Hb genes and isoforms are co-expressed in schizothoracinae fishes endemic to the Qinghai-Tibetan Plateau (QTP). Phylogenetic analyses of Hb genes indicated that the evolutionary relationships are not easily reconciled with the organismal phylogeny. Furthermore, evidence of positive selection was found in the Hb genes of schizothoracinae fishes through the selection pressure analysis. We demonstrated that positively selected sites likely facilitated the functional divergence of Hb isoforms. Taken together, this study indicated that the long-term maintenance of high Hb concentration may be a disadvantage for physiologically acclimating to high altitude hypoxia. Meanwhile, the genetically based modification of Hb-O₂ affinity in schizothoracinae fishes might facilitate the evolutionary adaptation to Tibetan aqueous environments.

Zebrafish Model of Neuroblastoma Metastasis

Zebrafish has emerged as an important animal model to study human diseases, especially cancer. Along with the robust transgenic and genome editing technologies applied in zebrafish modeling, the ease of maintenance, high-yield productivity, and powerful live imaging altogether make the zebrafish a valuable model system to study metastasis and cellular and molecular bases underlying this process in vivo. The first zebrafish neuroblastoma (NB) model of metastasis was developed by overexpressing two oncogenes, MYCN and LMO1, under control of the dopamine-beta-hydroxylase (dβh) promoter. Co-overexpressed MYCN and LMO1 led to the reduced latency and increased penetrance of neuroblastomagenesis, as well as accelerated distant metastasis of tumor cells. This new model reliably reiterates many key features of human metastatic NB, including involvement of clinically relevant and metastasis-associated genetic alterations; natural and spontaneous development of metastasis in vivo; and conserved sites of metastases. Therefore, the zebrafish model possesses unique advantages to dissect the complex process of tumor metastasis in vivo.

Zebrafish as a Neuroblastoma Model: Progress Made, Promise for the Future

For nearly a decade, researchers in the field of pediatric oncology have been using zebrafish as a model for understanding the contributions of genetic alternations to the pathogenesis of neuroblastoma (NB), and exploring the molecular and cellular mechanisms that underlie neuroblastoma initiation and metastasis. In this review, we will enumerate and illustrate the key advantages of using the zebrafish model in NB research, which allows researchers to: monitor tumor development in real-time; robustly manipulate gene expression (either transiently or stably); rapidly evaluate the cooperative interactions of multiple genetic alterations to disease pathogenesis; and provide a highly efficient and low-cost methodology to screen for effective pharmaceutical interventions (both alone and in combination with one another). This review will then list some of the common challenges of using the zebrafish model and provide strategies for overcoming these difficulties. We have also included visual diagram and figures to illustrate the workflow of cancer model development in zebrafish and provide a summary comparison of commonly used animal models in cancer research, as well as key findings of cooperative contributions between MYCN and diverse singling pathways in NB pathogenesis.

Reduced Oxygen as an Environmental Pressure in the Evolution of the Blind Mexican Cavefish

Extreme environmental features can drive the evolution of extreme phenotypes. Over the course of evolution, certain environmental changes may be so drastic that they lead to extinction. Conversely, if an organism adapts to harsh environmental changes, the adaptations may permit expansion of a novel niche. The interaction between environmental stressors and adaptive changes is well-illustrated by the blind Mexican cavefish, Astyanaxmexicanus, which has recurrently adapted to the stark subterranean environment. The transition from terrestrial rivers and streams (occupied by extant surface morphs of the same species) to the cave has been accompanied by the resorption of eyes, diminished pigmentation and reduced metabolism in cave-dwelling morphs. The principal features of caves most often associated with evolution of these common cave features are the absence of light and limited nutrition. However, a putatively essential cave feature that has received less attention is the frequently low concentration of oxygen within natural karst environments. Here, we review the potential role of limited oxygen as a critical environmental feature of caves in the Sierra de El Abra. Additionally, we review evidence that Astyanax cavefish may have evolved adaptive features enabling them to thrive in lower oxygen compared to their surface-dwelling counterparts.

The hemoglobin Gly16β1Asp polymorphism in turbot (Scophthalmus maximus) is differentially distributed across European populations

Turbot is an important flatfish widely distributed along the European coasts, whose fishery is centered in the North Sea. The commercial value of the species has boosted a successful aquaculture sector in Europe and China. Body growth is the main target of turbot breeding programs and is also a key trait related to local adaptation to temperature and salinity. Differences in growth rate and optimal growth temperature in turbot have been shown to be associated with a hemoglobin polymorphism reported more than 50 years ago. Here, we identified a Gly16Asp amino acid substitution in the β1 globin subunit by searching for genetic variation in the five functional globin genes within the whole annotated turbot genome. We predicted increased stability of the turbot hemoglobin by the replacement of the conserved Gly with the negative charged Asp residue that is consistent with the higher rate of αβ dimer assembly in the human J-Baltimore Gly16β->Asp mutant than in normal HbA. The turbot Hbβ1-Gly16 variant dominated in the northern populations examined, particularly in the Baltic Sea, while the Asp allele showed elevated frequencies in southern populations and was the prevalent variant in the Adriatic Sea. Body weight did not associate with the Hbβ1 genotypes at farming conditions (i.e., high oxygen levels, feeding ad libitum) after analyzing 90 fish with high growth dispersal from nine turbot families. Nevertheless, all data at hand suggest that the turbot hemoglobin polymorphism has an adaptive significance in the variable wild conditions regarding temperature and oxygen availability.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

The Atlantic salmon whole blood transcriptome and how it relates to major locus maturation genotypes and other tissues

The Atlantic salmon (Salmo salar) is important to many ecosystems and local economies and has therefore become the focus of a broad range of research questions that have benefited from the availability of high-quality genomic resources. Albeit gene expression studies have been extensive for this species, the transcriptome information for Atlantic salmon whole blood has been lacking. A transcriptome of Atlantic salmon blood would be a valuable resource for future studies, especially those wishing to take non-lethal samples. Here, we report a whole blood transcriptome for Atlantic salmon constructed from twelve 8-month old salmon parr using RNA-seq. We identify transcriptomic proxies for the genotype at the major maturation timing locus vestigial-like 3 (vgll3). Differentially expressed genes between the early and late maturing genotypes showed overrepresented Gene Ontology (GO) terms with the strongest result linked to 13 ribosomal subunit genes. To assess how the whole blood gene expression profile relates to other tissues, we compare the blood transcriptome to the reference transcriptome of fourteen other tissue types using both a common PCA method and a novel method. The novel method compares transcriptomes when gene expression is visualised as a layer using thin-plate spline smoothers. Both methods found similar patterns with the blood transcriptome being quite unique compared to the transcription profiles of other tissues.

Evolution of an Amniote-Specific Mechanism for Modulating Ubiquitin Signaling via Phosphoregulation of the E2 Enzyme UBE2D3

Genetic variation in the enzymes that catalyze posttranslational modification of proteins is a potentially important source of phenotypic variation during evolution. Ubiquitination is one such modification that affects turnover of virtually all of the proteins in the cell in addition to roles in signaling and epigenetic regulation. UBE2D3 is a promiscuous E2 enzyme, which acts as an ubiquitin donor for E3 ligases that catalyze ubiquitination of developmentally important proteins. We have used protein sequence comparison of UBE2D3 orthologs to identify a position in the C-terminal α-helical region of UBE2D3 that is occupied by a conserved serine in amniotes and by alanine in anamniote vertebrate and invertebrate lineages. Acquisition of the serine (S138) in the common ancestor to modern amniotes created a phosphorylation site for Aurora B. Phosphorylation of S138 disrupts the structure of UBE2D3 and reduces the level of the protein in mouse embryonic stem cells (ESCs). Substitution of S138 with the anamniote alanine (S138A) increases the level of UBE2D3 in ESCs as well as being a gain of function early embryonic lethal mutation in mice. When mutant S138A ESCs were differentiated into extraembryonic primitive endoderm, levels of the PDGFRα and FGFR1 receptor tyrosine kinases were reduced and primitive endoderm differentiation was compromised. Proximity ligation analysis showed increased interaction between UBE2D3 and the E3 ligase CBL and between CBL and the receptor tyrosine kinases. Our results identify a sequence change that altered the ubiquitination landscape at the base of the amniote lineage with potential effects on amniote biology and evolution.

Genomic Survey of Tyrosine Kinases Repertoire in Electrophorus electricus With an Emphasis on Evolutionary Conservation and Diversification

Tyrosine kinases (TKs) play key roles in the regulation of multicellularity in organisms and involved primarily in cell growth, differentiation, and cell-to-cell communication. Genome-wide characterization of TKs has been conducted in many metazoans; however, systematic information regarding this superfamily in Electrophorus electricus (electric eel) is still lacking. In this study, we identified 114 TK genes in the E electricus genome and investigated their evolution, molecular features, and domain architecture using phylogenetic profiling to gain a better understanding of their similarities and specificity. Our results suggested that the electric eel TK (EeTK) repertoire was shaped by whole-genome duplications (WGDs) and tandem duplication events. Compared with other vertebrate TKs, gene members in Jak, Src, and EGFR subfamily duplicated specifically, but with members lost in Eph, Axl, and Ack subfamily in electric eel. We also conducted an exhaustive survey of TK genes in genomic databases, identifying 1674 TK proteins in 31 representative species covering all the main metazoan lineages. Extensive evolutionary analysis indicated that TK repertoire in vertebrates tended to be remarkably conserved, but the gene members in each subfamily were very variable. Comparative expression profile analysis showed that electric organ tissues and muscle shared a similar pattern with specific highly expressed TKs (ie, epha7, musk, jak1, and pdgfra), suggesting that regulation of TKs might play an important role in specifying an electric organ identity from its muscle precursor. We further identified TK genes exhibiting tissue-specific expression patterns, indicating that members in TKs participated in subfunctionalization representing an evolutionary divergence required for the performance of different tissues. This work generates valuable information for further gene function analysis and identifying candidate TK genes reflecting their unique tissue-function specializations in electric eel.

Analysis of Multiplicity of Hypoxia-Inducible Factors in the Evolution of Triplophysa Fish (Osteichthyes: Nemacheilinae) Reveals Hypoxic Environments Adaptation to Tibetan Plateau

HIF (Hypoxia-inducible factor) gene family members function as master regulators of cellular and systemic oxygen homeostasis during changes in oxygen availability. Qinghai-Tibet Plateau is a natural laboratory for for long-term hypoxia and cold adaptation. In this context, T. scleroptera that is restricted to >3500 m high-altitude freshwater rivers was selected as the model to compare with a representative species from the plain, P. dabryanus. We cloned different HIF-α and carried out a phylogenetic analysis from invertebrates to vertebrates for identifying HIF-α genes and analyzing their evolutionary history. Intriguingly, the HIF-α has undergone gene duplications might be due to whole-genome duplication (WGD) events during evolution. PAML analysis indicated that HIF-1αA was subjected to positive selection acted on specific sites in Triplophysa lineages. To investigate the relationship between hypoxia adaptation and the regulation of HIF-α stability by pVHL in plateau and plain fish, a series of experiments were carried out. Comparison the luciferase transcriptional activity and protein levels of HIF-αs and the differing interactions of HIF-αs with pVHL, show clear differences between plateau and plain fish. T. scleroptera pVHL could enhance HIF-α transcriptional activity under hypoxia, and functional validation through pVHL protein mutagenesis showed that these mutations increased the stability of HIF-α and its hetero dimerization affinity to ARNT. Our research shows that missense mutations of pVHL induced evolutionary molecular adaptation in Triplophysa fishes living in high altitude hypoxic environments.

Oxygenation properties of hemoglobin and the evolutionary origins of isoform multiplicity in an amphibious air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris)

Among the numerous lineages of teleost fish that have independently transitioned from obligate water breathing to facultative air breathing, evolved properties of hemoglobin (Hb)-O₂ transport may have been shaped by the prevalence and severity of aquatic hypoxia (which influences the extent to which fish are compelled to switch to aerial respiration) as well as the anatomical design of air-breathing structures and the cardiovascular system. Here, we examined the structure and function of Hbs in an amphibious, facultative air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris). We also characterized the genomic organization of the globin gene clusters of the species and we integrated phylogenetic and comparative genomic analyses to unravel the duplicative history of the genes that encode the subunits of structurally distinct mudskipper Hb isoforms (isoHbs). The B. pectinirostris isoHbs exhibit high intrinsic O₂ affinities, similar to those of hypoxia-tolerant, water-breathing teleosts, and remarkably large Bohr effects. Genomic analysis of conserved synteny revealed that the genes that encode the α-type subunits of the two main adult isoHbs are members of paralogous gene clusters that represent products of the teleost-specific whole-genome duplication. Experiments revealed no appreciable difference in the oxygenation properties of co-expressed isoHbs in spite of extensive amino acid divergence between the alternative α-chain subunit isoforms. It therefore appears that the ability to switch between aquatic and aerial respiration does not necessarily require a division of labor between functionally distinct isoHbs with specialized oxygenation properties.

Multiplicity and Polymorphism of Fish Hemoglobins

The diversity of fish hemoglobins and the association with oxygen availability and physiological requirements during the life cycle has attracted scientists since the first report on multiple hemoglobin in fishes (Buhler and Shanks 1959). The functional heterogeneity of the fish hemoglobins enables many species to tolerate hypoxic conditions and exhausting swimming, but also to maintain the gas pressure in the swim bladder at large depths. The hemoglobin repertoire has further increased in various species displaying polymorphic hemoglobin variants differing in oxygen binding properties. The multiplicity of fish hemoglobins as particularly found in the tetraploid salmonids strongly contrasts with the complete loss of hemoglobins in Antarctic icefishes and illustrates the adaptive radiation in the oxygen transport of this successful vertebrate group.

Tandem repeats lead to sequence assembly errors and impose multi-level challenges for genome and protein databases

The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with 'ready-to-use' deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others.

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.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

Duplication of toll-like receptor 22 in teleost fishes

The TLRs of teleost fishes have distinct features and are highly diverse, but the duplication characteristics and expression patterns of the tlr22 gene remain unclear. Here, we identified paralogous tlr22 genes in 13 teleost fishes by screening available fish genomic resources and using molecular cloning. We then conducted comprehensive bioinformatics analyses and investigated spatiotemporal differences in the expression patterns of the tlr22 genes in G. eckloni. The results indicated that more than three paralogous tlr22 genes were possessed by some teleost fishes. Of these, tlr22c is specific to some subfamilies of the Cyprinidae (e.g., Barbinae, Cyprininae, Schizothoracinae, and Leuciscinae). Phylogenetic and syntenic analyses showed that the paralogous tlr22 genes originated from two single-gene duplication events. Molecular clock calculations dated the two gene duplication events at 49.5 and 39.3 MYA, which is before the common carp-specific genome duplication event and well after the fish-specific genome duplication. Gene duplication of tlr22 was followed by gene loss or pseudogene events in certain lineages. Spatiotemporal expression differences between the three duplicated tlr22 genes from G. eckloni suggested that these genes diverged functionally after gene duplication.

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.

Long-read sequence capture of the haemoglobin gene clusters across codfish species

Combining high-throughput sequencing with targeted sequence capture has become an attractive tool to study specific genomic regions of interest. Most studies have so far focused on the exome using short-read technology. These approaches are not designed to capture intergenic regions needed to reconstruct genomic organization, including regulatory regions and gene synteny. Here, we demonstrate the power of combining targeted sequence capture with long-read sequencing technology for comparative genomic analyses of the haemoglobin (Hb) gene clusters across eight species separated by up to 70 million years. Guided by the reference genome assembly of the Atlantic cod (Gadus morhua) together with genome information from draft assemblies of selected codfishes, we designed probes covering the two Hb gene clusters. Use of custom-made barcodes combined with PacBio RSII sequencing led to highly continuous assemblies of the LA (~100 kb) and MN (~200 kb) clusters, which include syntenic regions of coding and intergenic sequences. Our results revealed an overall conserved genomic organization of the Hb genes within this lineage, yet with several, lineage-specific gene duplications. Moreover, for some of the species examined, we identified amino acid substitutions at two sites in the Hbb1 gene as well as length polymorphisms in its regulatory region, which has previously been linked to temperature adaptation in Atlantic cod populations. This study highlights the use of targeted long-read capture as a versatile approach for comparative genomic studies by generation of a cross-species genomic resource elucidating the evolutionary history of the Hb gene family across the highly divergent group of codfishes.

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1-0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation.

Gene Turnover and Diversification of the α- and β-Globin Gene Families in Sauropsid Vertebrates

The genes that encode the α- and β-chain subunits of vertebrate hemoglobin have served as a model system for elucidating general principles of gene family evolution, but little is known about patterns of evolution in amniotes other than mammals and birds. Here, we report a comparative genomic analysis of the α- and β-globin gene clusters in sauropsids (archosaurs and nonavian reptiles). The objectives were to characterize changes in the size and membership composition of the α- and β-globin gene families within and among the major sauropsid lineages, to reconstruct the evolutionary history of the sauropsid α- and β-globin genes, to resolve orthologous relationships, and to reconstruct evolutionary changes in the developmental regulation of gene expression. Our comparisons revealed contrasting patterns of evolution in the unlinked α- and β-globin gene clusters. In the α-globin gene cluster, which has remained in the ancestral chromosomal location, evolutionary changes in gene content are attributable to the differential retention of paralogous gene copies that were present in the common ancestor of tetrapods. In the β-globin gene cluster, which was translocated to a new chromosomal location, evolutionary changes in gene content are attributable to differential gene gains (via lineage-specific duplication events) and gene losses (via lineage-specific deletions and inactivations). Consequently, all major groups of amniotes possess unique repertoires of embryonic and postnatally expressed β-type globin genes that diversified independently in each lineage. These independently derived β-type globins descend from a pair of tandemly linked paralogs in the most recent common ancestor of sauropsids.

Hypoxia-induced changes in hemoglobins of Lake Victoria cichlids

In a previous study, broods of the Lake Victoria cichlid Haplochromis ishmaeli raised under hypoxic or normoxic conditions showed striking differences in isohemoglobin (isoHb) pattern that were not observed in two other cichlids that do not belong to the Lake Victoria species flock. We therefore hypothesized that the adaptive mechanism seen in H. ishmaeli in response to hypoxia constitutes a trait that the Lake Victoria species flock inherited from ancestors that lived in hypoxic environments. We tested this hypothesis by designing split-brood experiments with three other representative species from the same species flock: the insectivorous Haplochromis thereuterion, the mollusk-shelling Platytaeniodus degeni and the zooplanktivorous Haplochromis piceatus, while keeping H. ishmaeli as a reference. Split broods were raised, under either normoxia or hypoxia. All hypoxia-raised (HR) individuals of each of the four species exhibited a distinctly different isoHb pattern compared with their normoxia-raised (NR) siblings. The hemoglobin of HR H. thereuterion showed higher O₂ affinity compared with NR siblings particularly in the presence of ATP and GTP, indicating that blood of HR juveniles has significantly improved O₂-binding affinity under hypoxic conditions. We also tested the capacity to acclimate at greater age in two species by reversing the O₂ condition after 7 (H. thereuterion) and 4 (H. ishmaeli) months. After reacclimation for 1 and 2 months, respectively, we found incomplete reversal with intermediate isoHb patterns. As three of the four species do not encounter hypoxic conditions in their environment, this unique trait seems to be a relic inherited from predecessors that lived in hypoxic environments.

Role of elasmobranchs and holocephalans in understanding peptide evolution in the vertebrates: Lessons learned from gonadotropin releasing hormone (GnRH) and corticotropin releasing factor (CRF) phylogenies

The cartilaginous fishes (Class Chondrichthyes) comprise two morphologically distinct subclasses; Elasmobranchii and Holocephali. Evidence indicates early divergence of these subclasses, suggesting monophyly of their lineage. However, such a phylogenetic understanding is not yet developed within two highly conserved peptide lineages, GnRH and CRF. Various GnRH forms exist across the Chondrichthyes. Although 4-7 immunoreactive forms have been described in Elasmobranchii, only one has been elucidated in Holocephali. In contrast, Chondrichthyan CRF phylogeny follows a pattern more consistent with vertebrate evolution. For example, three forms are expressed within the lamprey, with similar peptides present within the genome of the Callorhinchus milii, a holocephalan. Although these findings are consistent with recent evidence regarding the phylogenetic age of Chondrichthyan lineages, CRF evolution in vertebrates remains elusive. Assuming that the Elasmobranchii and Holocephali are part of a monocladistic clade within the Chondrichthyes, we interpret the findings of GnRH and CRF to be products of their respective lineages.

Hints on T cell responses in a fish-parasite model: Enteromyxum leei induces differential expression of T cell signature molecules depending on the organ and the infection status

BACKGROUD

Enteromyxum leei is a myxozoan parasite that produces a slow-progressing intestinal disease. This parasite invades the paracellular space of the intestinal epithelium and progresses from the posterior to the anterior intestine. The aim of the present study was to gain insights into fish T cell responses in the gilthead sea bream-E. leei infection model using a PCR-array with 30 signature molecules for different leukocyte responses in head kidney, spleen, anterior and posterior intestine.

RESULTS

The PCR-array results suggest that E. leei induced migration of T cells from head kidney to intestines where T_H1, CTL and T_H17 profiles were activated and kept in balance by the upregulation of regulatory cytokines. These results were partially validated by the use of cross-reacting antibodies and BrdU immunostaining to monitor proliferation. Zap70 immunostaining supported the increased number of T cells in the anterior intestine detected by gene expression, but double staining with BrdU did not show active proliferation of this cell type at a local level, supporting the migration from lymphohaematopoietic tissues to the site of infection. Global analyses of the expression profiles revealed a clear separation between infected and exposed, but non-infected fish, more evident in the target organ. Exposed, non-infected animals showed an intermediate phenotype closer to the control fish.

CONCLUSIONS

These results evidence a clear modulation of the T cell response of gilthead sea bream upon E. leei infection. The effects occurred both at local and systemic levels, but the response was stronger and more specific at the site of infection, the intestine. Altogether, this research poses a promising basis to understand the response against this important parasite and establish effective preventive or palliative measures.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Molecular evolution of myoglobin in the Tibetan Plateau endemic schizothoracine fish (Cyprinidae, Teleostei) and tissue-specific expression changes under hypoxia

Myoglobin (Mb) is an oxygen-binding hemoprotein that was once thought to be exclusively expressed in oxidative myocytes of skeletal and cardiac muscle where it serves in oxygen storage and facilitates intracellular oxygen diffusion. In this study, we cloned the coding sequence of the Mb gene from four species, representing three groups, of the schizothoracine fish endemic to the Qinghai-Tibetan Plateau (QTP), then conducted molecular evolution analyses. We also investigated tissue expression patterns of Mb and the expression response to moderate and severe hypoxia at the mRNA and protein levels in a representative of the highly specialized schizothoracine fish species, Schizopygopsis pylzovi. Molecular evolution analyses showed that Mb from the highly specialized schizothoracine fish have undergone positive selection and one positively selected residue (81L) was identified, which is located in the F helix, close to or in contact with the heme. We present tentative evidence that the Mb duplication event occurred in the ancestor of the schizothoracine and Cyprininae fish (common carp and goldfish), and that the Mb2 paralog was subsequently lost in the schizothoracine fish. In S. pylzovi, Mb mRNA is expressed in various tissues with the exception of the intestine and gill, but all such tissues, including the liver, muscle, kidney, brain, eye, and skin, expressed very low levels of Mb mRNA (< 8.0%) relative to that of the heart. The trace levels of Mb expression in non-muscle tissues are perhaps the major reason why non-muscle Mb remained undiscovered for so long. The expression response of the Mb gene to hypoxia at the mRNA and protein levels was strikingly different in S. pylzovi compared to that found in the common carp, medaka, zebrafish, and goldfish, suggesting that the hypoxia response of Mb in fish may be species and tissue-specific. Notably, severe hypoxia induced significant expression of Mb at the mRNA and protein levels in the S. pylzovi heart, which suggests Mb has a major role in the supply of oxygen to the heart of Tibetan Plateau fish.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical “non-blood” globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Identification, molecular evolution of toll-like receptors in a Tibetan schizothoracine fish (Gymnocypris eckloni) and their expression profiles in response to acute hypoxia

Hypoxia plays an important role in regulating a variety of physiological responses as well as in pathological situations, but to date the roles of Toll-like receptors (TLRs) in fish in response to hypoxia are still poorly understood. Here, we sequenced the transcriptome of G. eckloni and identified the members of TLR family by scanning transcriptome, and then investigated the expression profiles of a complete set of TLRs in G. eckloni in response to acute hypoxia (4 h at DO = 0.3 ± 0.1 mg/L). The de novo-assembled transcriptome consisted of a total of 162,235 transcripts, further clustered into 110,231 unigenes. Based on the transcriptome, a total of 18 TLRs were identified in G. eckloni, and of them three TLRs (TLR5, TLR8 and TLR22) possessed two distinct paralogous genes. The duplicated genes of TLR22 were discovered for the first time in cyprinid fish, but did not origin from a recent duplication event. Of them TLR22b may be specific for schizothoracine fish, at least for G. eckloni. Phylogenetic analysis supported the classification of TLRs into six families as in other vertebrates but was partly different from the previous study. The sliding window analysis showed strong signals of positive selection in TLR2, TLR 4, TLR 5a, TLR 7, TLR 19, TLR 20, TLR 21, TLR 22a and TLR 22b, but most codons under positive selection were located in the putative LRR regions. The mRNA expression of most TLRs in head kidney, spleen and gill decreased significantly or remained unchanged under acute hypoxia, whereas acute hypoxia increased expressions of TLR2 and TLR3 in head kidney, of TLR8a, TLR12 and TLR19 in spleen, and of TLR1 in gill, suggesting tissues-specific expressions of TLRs play important roles in mediating innate immune responses for host defense against tissue damages or physiological changes induced by hypoxia.

Evolution of Hemoglobin Genes in Codfishes Influenced by Ocean Depth

Understanding the genetic basis of adaptation is one of the main enigmas of evolutionary biology. Among vertebrates, hemoglobin has been well documented as a key trait for adaptation to different environments. Here, we investigate the role of hemoglobins in adaptation to ocean depth in the diverse teleost order Gadiformes, with species distributed at a wide range of depths varying in temperature, hydrostatic pressure and oxygen levels. Using genomic data we characterized the full hemoglobin (Hb) gene repertoire for subset of species within this lineage. We discovered a correlation between expanded numbers of Hb genes and ocean depth, with the highest numbers in species occupying shallower, epipelagic regions. Moreover, we demonstrate that the Hb genes have functionally diverged through diversifying selection. Our results suggest that the more variable environment in shallower water has led to selection for a larger Hb gene repertoire and that Hbs have a key role in adaptive processes in marine environments.

Acclimation to prolonged hypoxia alters hemoglobin isoform expression and increases hemoglobin oxygen affinity and aerobic performance in a marine fish

Hemoglobin (Hb) multiplicity is common in fish, yet despite its ubiquitous nature, the functional significance is unclear. Here we explore the hypothesis that Hb multiplicity plays a role in hypoxia tolerance using the red drum (Sciaenops ocellatus). Red drum is an economically and ecologically important species native to coastal regions and estuaries of the Gulf of Mexico – habitats that routinely experience pronounced hypoxic events. Using a transcriptomic approach, we demonstrate that red drum red blood cells express 7 and 5 Hbα and Hbβ isoforms, respectively. Phylogenetic analysis grouped these isoforms into distinct isoHb clades, and provided evidence of lineage specific expression of particular isoHbs. In normoxia, three isoHbs predominated (Hbα-3.1, -3.2, and Hbβ-3.1). A three-week hypoxia acclimation (48 mmHg) resulted in significant up-regulation of Hbα-2, Hbα-3.2, and Hbβ-3.1, effectively switching the predominantly expressed isoforms. Changes in subunit expression were correlated with a decrease in non-stripped hemolysate P₅₀. Similarly, hypoxia acclimation resulted in a 20% reduction in whole animal critical oxygen threshold (P_crit). Hypoxia acclimation was not associated with changes in gill morphology, hematocrit, or relative ventricular mass. Overall, these data provide support for the hypothesis that Hb isoform switching can provide a physiological benefit to counteract environmental stress in fishes.

Expansion by whole genome duplication and evolution of the sox gene family in teleost fish

It is now recognized that several rounds of whole genome duplication (WGD) have occurred during the evolution of vertebrates, but the link between WGDs and phenotypic diversification remains unsolved. We have investigated in this study the impact of the teleost-specific WGD on the evolution of the sox gene family in teleostean fishes. The sox gene family, which encodes for transcription factors, has essential role in morphology, physiology and behavior of vertebrates and teleosts, the current largest group of vertebrates. We have first redrawn the evolution of all sox genes identified in eleven teleost genomes using a comparative genomic approach including phylogenetic and synteny analyses. We noticed, compared to tetrapods, an important expansion of the sox family: 58% (11/19) of sox genes are duplicated in teleost genomes. Furthermore, all duplicated sox genes, except sox17 paralogs, are derived from the teleost-specific WGD. Then, focusing on five sox genes, analyzing the evolution of coding and non-coding sequences, as well as the expression patterns in fish embryos and adult tissues, we demonstrated that these paralogs followed lineage-specific evolutionary trajectories in teleost genomes. This work, based on whole genome data from multiple teleostean species, supports the contribution of WGDs to the expansion of gene families, as well as to the emergence of genomic differences between lineages that might promote genetic and phenotypic diversity in teleosts.

Early exposure to chronic hypoxia induces short- and long-term regulation of hemoglobin gene expression in European sea bass (Dicentrarchus labrax)

European sea bass (Dicentrarchus labrax) inhabits coastal waters and may be exposed to hypoxia at different life stages, requiring physiological and behavioral adaptation. In the present study, we attempted to determine whether regulation of hemoglobin (Hb) gene expression plays a role in the physiological response to chronic moderate hypoxia in whole larvae and hematopoietic tissues (head kidney and spleen) of juveniles. We also tested the hypothesis that hypoxia exposure at the larval stage could induce a long-term effect on the regulation of Hb gene expression. For this purpose, D. labrax were exposed to a non-lethal hypoxic condition (40% air saturation) at the larval stage from 28 to 50 days post-hatching (dph) and/or at the juvenile stage from 196 to 296 dph. Data obtained from larvae indicate that hypoxia induced a subtype-specific regulation of Hb gene expression, with a significant decrease of MN-Hbα3, MN-Hbβ4 and MN-Hbβ5 and increase of MN-Hbα2, LA-Hbα1 and LA-Hbβ1 transcript levels. Hypoxia did not induce regulation of Hb gene expression in juveniles, except in the head kidney for those that experienced hypoxia at the larval stage. The latter exhibited a significant hypoxia-induced stimulation of MN-Hbα2, LA-Hbα1 and LA-Hbβ1 gene expression, associated with stimulation of the PHD-3 gene involved in the hypoxia-inducible factor oxygen-sensing pathway. We conclude that subtype- and stage-specific regulation of Hb gene expression plays a role in the physiological response of D. labrax to cope with hypoxia and that early exposure to low oxygen concentration has a long-term effect on this response.

Reprimo tissue-specific expression pattern is conserved between zebrafish and human

Reprimo (RPRM), a member of the RPRM gene family, is a tumor-suppressor gene involved in the regulation of the p53-mediated cell cycle arrest at G2/M. RPRM has been associated with malignant tumor progression and proposed as a potential biomarker for early cancer detection. However, the expression and role of RPRM, as well as its family, are poorly understood and their physiology is as yet unstudied. In this scenario, a model system like the zebrafish could serve to dissect the role of the RPRM family members in vivo. Phylogenetic analysis reveals that RPRM and RPRML have been differentially retained by most species throughout vertebrate evolution, yet RPRM3 has been retained only in a small group of distantly related species, including zebrafish. Herein, we characterized the spatiotemporal expression of RPRM (present in zebrafish as an infraclass duplication rprma/rprmb), RPRML and RPRM3 in the zebrafish. By whole-mount in situ hybridization (WISH) and fluorescent in situ hybridization (FISH), we demonstrate that rprm (rprma/rprmb) and rprml show a similar spatiotemporal expression profile during zebrafish development. At early developmental stages rprmb is expressed in somites. After one day post-fertilization, rprm (rprma/rprmb) and rprml are expressed in the notochord, brain, blood vessels and digestive tube. On the other hand, rprm3 shows the most unique expression profile, being expressed only in the central nervous system (CNS). We assessed the expression patterns of RPRM gene transcripts in adult zebrafish and human RPRM protein product in tissue samples by RT-qPCR and immunohistochemistry (IHC) staining, respectively. Strikingly, tissue-specific expression patterns of the RPRM transcripts and protein are conserved between zebrafish and humans. We propose the zebrafish as a powerful tool to elucidate the both physiological and pathological roles of the RPRM gene family.

Gene Duplication and Evolutionary Innovations in Hemoglobin-Oxygen Transport

During vertebrate evolution, duplicated hemoglobin (Hb) genes diverged with respect to functional properties as well as the developmental timing of expression. For example, the subfamilies of genes that encode the different subunit chains of Hb are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different developmental stages. In some vertebrate taxa, functional differentiation between co-expressed Hb isoforms may also contribute to physiologically important divisions of labor.

Carp Il10a and Il10b exert identical biological activities in vitro, but are differentially regulated in vivo

We recently reported on the functional characterization of carp Il10. We showed that carp Il10 is able to downregulate proinflammatory activities by carp phagocytes and promote B cell proliferation, differentiation and antibody production as well as proliferation of memory T cells. Taking advantage of the recent annotation of the carp genome, we completed the sequence of a second il10 paralogue, named il10b, the presence of which was expected owing to the recent (8 million years ago) fourth round of whole genome duplication that occurred in common carp. In the present study we closely compared the two Il10 paralogues and show that Il10a and Il10b have almost identical gene structure, synteny, protein sequence as well as bioactivity on phagocytes. Although the two il10 paralogues show a large overlap in tissue expression, il10b has a low constitutive expression and is highly upregulated upon infection, whereas il10a is higher expressed under basal conditions but its gene expression remains constant during viral and parasitic infections. This differential regulation is most likely due to the observed differences in their promoter regions. Altogether our results demonstrate that gene duplication in carp, although recent, led to sub-functionalization and expression divergence rather than functional redundancy of the Il10 paralogues, yet with very similar protein sequences.