Organization of a β and α globin gene set in the teleost Atlantic cod, Gadus morhua

Structural and functional characterization of haemoglobin genes in Labeo catla: Insights into hypoxic adaptation and survival

Haemoglobin (HB) protein comprises four subunits: two identical α-subunits (HBA) and two identical β-subunits (HBB), encoded by the HBA and HBB genes. In this investigation, 5'/3' RACE PCR (Rapid Amplification of cDNA Ends) was used to obtain complete coding sequences (CDSs) of both the genes from farmed Labeo catla. The resulting CDSs were 432 base pairs and 447 base pairs for HBA and HBB, respectively, corresponding to 143 and 148 amino acids. Phylogenetic analysis revealed close relationships with other cyprinids, with Labeo rohita being the closest relative. Functional analysis and protein structure prediction were conducted using bioinformatics tools. Expression profiling of both genes was checked in various tissues under control (C) and hypoxic (H) conditions. Notably, under hypoxia, HBA and HBB genes were significantly upregulated (P < 0.05) initially, followed by a return to normal expression levels. Similar trends were observed for Hif1α (Hypoxia-inducible factor one alpha) and EPO (Erythropoietin) genes. Additionally, haematological indices also significantly increased corresponding to the gene expressions. However, with the decrease in the expression of these genes an onset of mortality was observed in the hypoxia (H) treated groups. The results of the current study explored the role of haemoglobin genes in adaptation to the hypoxic condition.

Nucleotide variation and balancing selection at the Ckma gene in Atlantic cod: analysis with multiple merger coalescent models

High-fecundity organisms, such as Atlantic cod, can withstand substantial natural selection and the entailing genetic load of replacing alleles at a number of loci due to their excess reproductive capacity. High-fecundity organisms may reproduce by sweepstakes leading to highly skewed heavy-tailed offspring distribution. Under such reproduction the Kingman coalescent of binary mergers breaks down and models of multiple merger coalescent are more appropriate. Here we study nucleotide variation at the Ckma (Creatine Kinase Muscle type A) gene in Atlantic cod. The gene shows extreme differentiation between the North (Canada, Greenland, Iceland, Norway, Barents Sea) and the South (Faroe Islands, North-, Baltic-, Celtic-, and Irish Seas) with F_ST > 0.8 between regions whereas neutral loci show no differentiation. This is evidence of natural selection. The protein sequence is conserved by purifying selection whereas silent and non-coding sites show extreme differentiation. The unfolded site-frequency spectrum has three modes, a mode at singleton sites and two high frequency modes at opposite frequencies representing divergent branches of the gene genealogy that is evidence for balancing selection. Analysis with multiple-merger coalescent models can account for the high frequency of singleton sites and indicate reproductive sweepstakes. Coalescent time scales vary with population size and with the inverse of variance in offspring number. Parameter estimates using multiple-merger coalescent models show that times scales are faster than under the Kingman coalescent.

Evolution and function of the globin intergenic regulatory regions of the antarctic dragonfishes (Notothenioidei: Bathydraconidae)

As the Southern Ocean cooled to -1.8 °C over the past 40 My, the teleostean clade Notothenioidei diversified and, under reduced selection pressure for an oxygen-transporting apparatus, became less reliant on hemoglobin and red blood cells. At the extreme of this trend, the crown group of Antarctic icefishes (Channichthyidae) lost both components of oxygen transport. Under the decreased selection scenario, we hypothesized that the Antarctic dragonfishes (Bathydraconidae, the red-blooded sister clade to the icefishes) evolved lower blood hemoglobin concentrations because their globin gene complexes (α- and β-globin gene pairs linked by a regulatory intergene) transcribe globin mRNAs less effectively than those of basal notothenioids (e.g., the Nototheniidae [notothens]). To test our hypothesis, we 1) sequenced the α/β-intergenes of the adult globin complexes of three notothen and eight dragonfish species and 2) measured globin transcript levels in representative species from each group. The typical nototheniid intergene was ∼3-4 kb in length. The bathydraconid intergenes resolved into three subclasses (long [3.8 kb], intermediate [3.0 kb], and short [1.5-2.3 kb]) that corresponded to the three subclades proposed for the taxon. Although they varied in length due to indels, the three notothen and eight dragonfish intergenes contained a conserved ∼90-nt element that we have previously shown to be required for globin gene transcription. Using the quantitative polymerase chain reaction, we found that globin mRNA levels in red cells from one notothen species and from one species of each dragonfish subclade were equivalent statistically. Thus, our results indicate that the bathydraconids have evolved adult globin loci whose regulatory intergenes tend to be shorter than those of the more basal nototheniids yet are equivalent in transcriptional efficacy. Their low blood hemoglobin concentrations are most likely due to reduction in hematocrit.

Genomic organization and gene expression of the multiple globins in Atlantic cod: conservation of globin-flanking genes in chordates infers the origin of the vertebrate globin clusters

Background

The vertebrate globin genes encoding the α- and β-subunits of the tetrameric hemoglobins are clustered at two unlinked loci. The highly conserved linear order of the genes flanking the hemoglobins provides a strong anchor for inferring common ancestry of the globin clusters. In fish, the number of α-β-linked globin genes varies considerably between different sublineages and seems to be related to prevailing physico-chemical conditions. Draft sequences of the Atlantic cod genome enabled us to determine the genomic organization of the globin repertoire in this marine species that copes with fluctuating environments of the temperate and Arctic regions.

Results

The Atlantic cod genome was shown to contain 14 globin genes, including nine hemoglobin genes organized in two unlinked clusters designated β5-α1-β1-α4 and β3-β4-α2-α3-β2. The diverged cod hemoglobin genes displayed different expression levels in adult fish, and tetrameric hemoglobins with or without a Root effect were predicted. The novel finding of maternally inherited hemoglobin mRNAs is consistent with a potential role played by fish hemoglobins in the non-specific immune response. In silico analysis of the six teleost genomes available showed that the two α-β globin clusters are flanked by paralogs of five duplicated genes, in agreement with the proposed teleost-specific duplication of the ancestral vertebrate globin cluster. Screening the genome of extant urochordate and cephalochordate species for conserved globin-flanking genes revealed linkage of RHBDF1, MPG and ARHGAP17 to globin genes in the tunicate Ciona intestinalis, while these genes together with LCMT are closely positioned in amphioxus (Branchiostoma floridae), but seem to be unlinked to the multiple globin genes identified in this species.

Conclusion

The plasticity of Atlantic cod to variable environmental conditions probably involves the expression of multiple globins with potentially different properties. The interspecific difference in number of fish hemoglobin genes contrasts with the highly conserved synteny of the flanking genes. The proximity of globin-flanking genes in the tunicate and amphioxus genomes resembles the RHBDF1-MPG-α-globin-ARHGAP17-LCMT linked genes in man and chicken. We hypothesize that the fusion of the three chordate linkage groups 3, 15 and 17 more than 800 MYA led to the ancestral vertebrate globin cluster during a geological period of increased atmospheric oxygen content.

Integrating the markers Pan I and haemoglobin with the genetic linkage map of Atlantic cod (Gadus morhua)

Background

Haemoglobin (Hb) and pantophysin (Pan I) markers have been used intensively in population studies of Atlantic cod (Gadus morhua) and in the analysis of traits such as temperature tolerance, growth characteristics and sexual maturation. We used an Illumina GoldenGate panel and the KASPar SNP genotyping system to analyse SNPs in three Atlantic cod families, one of which was polymorphic at the Hb β1 locus, and to generate a genetic linkage map integrating Pan I and multiple Hb loci.

Findings

Data generated allowed the mapping of nine Hb loci, the Pan I locus, and other 122 SNPs onto an existing linkage genetic map for Atlantic cod. Four Hb genes (i.e. α1, α4, β1 and β5) have been mapped on linkage group (LG) 2 while the other five (i.e. α2, α3, β2, β3 and β4) were placed on LG18. Pan I was mapped on LG 1 using a newly developed KASPar assay for a SNP variable only in Pan I^A allelic variants. The new linkage genetic map presented here comprises 1046 SNPs distributed between 23 linkage groups, with a length of 1145.6 cM. A map produced by forcing additional loci, resulting in a reduced goodness-of-fit for mapped markers, allowed the mapping of a total of 1300 SNPs. Finally, we compared our genetic linkage map data with the genetic linkage map data produced by a different group and identified 29 shared SNPs distributed on 10 different linkage groups.

Conclusions

The genetic linkage map presented here incorporates the marker Pan I, together with multiple Hb loci, and integrates genetic linkage data produced by two different research groups. This represents a useful resource to further explore if Pan I and Hbs or other genes underlie quantitative trait loci (QTL) for temperature sensitivity/tolerance or other phenotypes.

Multiple linked β and α globin genes in Atlantic cod: A PCR based strategy of genomic exploration

Allozyme variation in Atlantic cod hemoglobins shows various signs of natural selection. We report a genomic exploration of globin genes in this non-model organism. Applying a PCR based strategy with a strict criterion of phylogenetically informative sites we estimate the number of linked β and α globin genes. We estimate PCR error rate by PCR of cloned DNA and recloning and by analysis of singleton variable sites among clones. Based on the error rate we exclude variable sites so that the remaining variation meets successively stricter criteria of doubleton and triplet variable site. Applying these criteria we find ten clusters of linked β/α globin genes in the genome of Atlantic cod. Six variable amino acid changes in both genes were found in linkage disequilibrium with silent nucleotide substitutions. A phylogenetic tree, based on our strictly phylogenetically informative sites among 57 clones from 19 individuals, is split into two major branches by an amino acid change in a β gene. This change is supported by extensive linkage disequilibrium between the amino acid change and numerous other phylogenetically informative silent nucleotide sites. The different gene sets in the genome may represent different loci encoding different globins and/or allelic variation at some loci.