Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “pseudogene”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68Linked to original sources
PubMed1996

Amylase gene structures in primates: retroposon insertions and promoter evolution.

Amylase transcription in the human salivary gland results from the evolutionary juxtaposition of two inserted elements, a gamma-actin pseudogene and an endogenous retrovirus, to create an unusual salivary-specific promoter. We utilized these structures as molecular tags to characterize the amylase genes in extant primates by polymerase chain reaction amplification of promoter fragments from genomic DNA. Six distinct amylase promoter structures were identified, which allowed us to infer the structures of common ancestors and trace the evolution of the modern human amylase promoters. Our data show that integration of the pseudogene and retrovirus were evolutionarily recent events. The gamma-actin pseudogene integrated after the divergence of the New World monkeys from the primate ancestral tree, and the retrovirus integrated later, after the divergence of the Old World monkeys. The New World monkey amylase promoter represents the mammalian amylase precursor structure before integration of the two retroposons. Two distinct amylase genes were identified in the Old World monkeys, one with a complete gamma-actin pseudogene insert and another novel structure with a truncation of the gamma-actin sequences. We demonstrated abundant amylase expression in the saliva of an Old World monkey, indicating that the endogenous retrovirus is not required for amylase transcription in the primate salivary gland.

Actins↗
PubMed2000

Frequent assimilation of mitochondrial DNA by grasshopper nuclear genomes.

Multiple copies of mitochondrial-like DNA were found in the brown mountain grasshopper, Podisma pedestris (Orthoptera: Acrididae), paralogous to COI and ND5 regions. The same was discovered using the ND5 regions of nine other grasshopper species from four separate subfamilies (Podisminae, Calliptaminae, Cyrtacanthacridinae, and Gomphocerinae). The extra ND5-like sequences were shown to be nuclear in the desert locust, Schistocerca gregaria (Cyrtacanthacridinae), and probably so in P. pedestris and an Italopodisma sp. (Podisminae). Eighty-seven different ND5-like nuclear mitochondrial pseudogenes (Numts) were sequenced from 12 grasshopper individuals. Different nuclear mitochondrial pseudogenes, if descended from the same mitochondrial immigrant, will have diverged from each other under no selective constraints because of their loss of functionality. Evidence of selective constraints in the differences between any two Numt sequences (e.g., if most differences are at third positions of codons) implies that they have separate mitochondrial origins. Through pairwise comparisons of pseudogene sequences, it was established that there have been at least 12 separate mtDNA integrations into P. pedestris nuclear genomes. This is the highest reported rate of horizontal transfer between organellar and nuclear genomes within a single animal species. The occurrence of numerous mitochondrial pseudogenes in nuclear genomes derived from separate integration events appears to be a common phenomenon among grasshoppers. More than one type of mechanism appears to have been involved in generating the observed grasshopper Numts.

Animals↗
PubMed2005

The vlhA loci of Mycoplasma synoviae are confined to a restricted region of the genome.

Mycoplasma synoviae, a major pathogen of poultry, contains a single expressed, full-length vlhA gene encoding its haemagglutinin, and a large number of vlhA pseudogenes that can be recruited by multiple site-specific recombination events to generate chimaeric variants of the expressed gene. The position and distribution of the vlhA pseudogene regions, and their relationship with the expressed gene, have not been investigated. To determine the relationship between these regions, a physical map of the M. synoviae genome was constructed using the restriction endonucleases SmaI, I-CeuI, BsiWI, ApaI and XhoI and radiolabelled probes for rrnA, recA and tufA. A cloned fragment encoding the unique portion of the expressed vlhA gene and two PCR products containing conserved regions of the ORF 3 and ORF 6 vlhA pseudogenes were used to locate the regions containing these genes on the map. The chromosome of M. synoviae was found to be 890.4 kb and the two rRNA operons were in the same orientation. Both the expressed vlhA gene and the vlhA pseudogenes were confined to the same 114 kb region of the chromosome. These findings indicate that, unlike Mycoplasma gallisepticum, in which the vlhA genes are located in several loci around the chromosome and in which antigenic variation is generated by alternating transcription of over 40 translationally competent genes, M. synoviae has all of the vlhA sequences clustered together, suggesting that close proximity is needed to facilitate the site-specific recombinations used to generate diversity in the expressed vlhA gene.

Animals↗
PubMed2006

Identification and targeted disruption of the mouse gene encoding ESG1 (PH34/ECAT2/DPPA5).

BACKGROUND: Embryonic stem cell-specific gene (ESG) 1, which encodes a KH-domain containing protein, is specifically expressed in early embryos, germ cells, and embryonic stem (ES) cells. Previous studies identified genomic clones containing the mouse ESG1 gene and five pseudogenes. However, their chromosomal localizations or physiological functions have not been determined. RESULTS: A Blast search of mouse genomic databases failed to locate the ESG1 gene. We identified several bacterial artificial clones containing the mouse ESG1 gene and an additional ESG1-like sequence with a similar gene structure from chromosome 9. The ESG1-like sequence contained a multiple critical mutations, indicating that it was a duplicated pseudogene. The 5' flanking region of the ESG1 gene, but not that of the pseudogene, exhibited strong enhancer and promoter activity in undifferentiated ES cells by luciferase reporter assay. To study the physiological functions of the ESG1 gene, we replaced this sequence in ES cells with a beta-geo cassette by homologous recombination. Despite specific expression in early embryos and germ cells, ESG1-/- mice developed normally and were fertile. We also generated ESG1-/- ES cells both by a second independent homologous recombination and directly from blastocysts derived from heterozygous intercrosses. Northern blot and western blot analyses confirmed the absence of ESG1 in these cells. These ES cells demonstrated normal morphology, proliferation, and differentiation. CONCLUSION: The mouse ESG1 gene, together with a duplicated pseudogene, is located on chromosome 9. Despite its specific expression in pluripotent cells and germ cells, ESG1 is dispensable for self-renewal of ES cells and establishment of germcells.

5' Flanking Region↗
PubMed2006

The pseudo-mitochondrial genome influences mistakes in heteroplasmy interpretation.

BACKGROUND: Nuclear mitochondrial pseudogenes (numts) are a potential source of contamination during mitochondrial DNA PCR amplification. This possibility warrants careful experimental design and cautious interpretation of heteroplasmic results. RESULTS: Here we report the cloning and sequencing of numts loci, amplified from human tissue and rho-zero (rho0) cells (control) with primers known to amplify the mitochondrial genome. This paper is the first to fully sequence 46 paralogous nuclear DNA fragments that represent the entire mitochondrial genome. This is a surprisingly small number due primarily to the primer sets used in this study, because prior to this, BLAST searches have suggested that nuclear DNA harbors between 400 to 1,500 paralogous mitochondrial DNA fragments. Our results indicate that multiple numts were amplified simultaneously with the mitochondrial genome and increased the load of pseudogene signal in PCR reactions. Further, the entire mitochondrial genome was represented by multiple copies of paralogous nuclear sequences. CONCLUSION: These findings suggest that mitochondrial genome disease-associated biomarkers must be rigorously authenticated to preclude any affiliation with paralogous nuclear pseudogenes. Importantly, the common perception that mitochondrial template "swamps" numts loci precluding detectable amplification, depends on the region of the mitochondrial genome targeted by the PCR reaction and the number of pseudogene loci that may co-amplify. Cloning and relevant sequencing data will facilitate the correct interpretation. This is the first complete, wet-lab characterization of numts that represent the entire mitochondrial genome.

Base Sequence↗
PubMed2004

'Frankenstein genes', or the Mad Magazine version of the human pseudogenome.

Annotation of the human genome is inching forward. Seven human chromosomes have now been fully annotated, covering 17 per cent of the genome, and at least one chromosome has been re-annotated. The enormity of the task forces a dependence on automated tools for detecting and assembling the genes, followed by hand curation to correct errors and polish the gene models. The accuracy of gene prediction algorithms is very good for internal exons from intact genes, but these programs do peculiar and exasperating things to pseudogenes. These programs can actually resurrect pseudogenes from the dead, making them into viable gene models for intact proteins, albeit science-fictional proteins. This process is demonstrated for four human pseudogenes from the cytochrome P450 family and one putatively functional P450 gene, CYP2U1, having a non-consensus intron boundary. These examples are offered as a call-to-arms to improve pseudogene prediction as an art in itself, and not as a by-product of gene annotation. Failure to do so will flood the databases with thousands of false-positive predictions. Indeed, they are already there.

Amino Acid Sequence↗
PubMed1998

Cloning and mapping of the mouse Gpx2 gene encoding gastrointestinal glutathione peroxidase.

Gastrointestinal glutathione peroxidase (GSHPx-GI) is an enzyme expressed in intestinal epithelial cells and may reduce hydroperoxides generated from the ingested diet. We isolated a genomic clone containing the mouse Gpx2 gene encoding 190 amino acids of GSHPx-GI. This gene is composed of two exons and an intron. The nucleotide sequence of the coding region was 89.9% identical with that of the human GPX2 gene. A TGA opal codon predicted to encode a selenocysteine was identified at codon 40. A genomic clone containing a pseudogene for the Gpx2 gene was also isolated. The nucleotide sequence of the pseudogene was 98.3% identical with that of the mouse Gpx2 gene and showed characteristics of a processed pseudogene. Linkage analysis using backcross mouse progeny indicated the mouse Cpx2 gene and its pseudogene to be located on mouse chromosomes 12 and 7, respectively.

Amino Acid Sequence↗
PubMed2004

Absence of the functional Myosin heavy chain 2b isoform in equine skeletal muscles.

Nucleotide sequences which included the full coding region for three types of myosin heavy chain (MyHC) isoforms were determined from equine skeletal muscles. The deduced amino acid sequences were 1937, 1938, and 1935 residues for the MyHC-2a, -2x, and -slow, respectively. No MyHC-2b isoform was amplified from the equine muscle cDNA except for one pseudogene fragment. One nucleotide was inserted in the coding region of the equine pseudogene product, a minute amount of which was expressed in the skeletal muscle. The 596 bp sequence of the equine MyHC pseudogene was categorized into the MyHC-2b genes on the phylogenetic tree of the mammalian MyHC genes. These results suggest that an ancestral MyHC-2b gene had lost its function and changed to a pseudogene during the course of horse history. The MyHC genes in some ungulates were analyzed through the PCR amplifications using the MyHC isoform-specific primers to confirm the presence of the MyHC-2b and -2x genes. The exon coding the 3' untranslated region of the MyHC-2x was successfully amplified from the all ungulates examined; however, that of the MyHC-2b gene was amplified only from horses, pigs and lesser mouse deer. The PCR analyses from rhinoceros, sika deer, moose, giraffes, water buffalo, bovine, Japanese serow and sheep genes implied the absence of the MyHC-2b-specific sequence in their genomes. These results suggest that the MyHC-2b gene independently lost its function in some ungulate species.

Amino Acid Sequence↗
PubMed2006

Complete sequence assembly and characterization of the C57BL/6 mouse Ig heavy chain V region.

The mechanisms that regulate variable (V) gene selection during the development of the mouse IgH repertoire are not fully understood, due in part to the absence of the complete locus sequence. To better understand these processes, we have assembled the entire 2.5-Mb mouse IgH (Igh) V region sequence of the C57BL/6 strain from public sequences and present the first complete annotated map of the region, including V genes, pseudogenes, repeats, and nonrepetitive intergenic sequences. In so doing, we have discovered a new V gene family, VH16. We have identified clusters of conserved region-specific intergenic sequences and have verified our assembly by genic and intergenic Southern blotting. We have observed that V pseudogenes are not evenly spread throughout the V region, but rather cluster together. The largest J558 family, which spans more than half of the locus, has two strikingly different domains, which suggest points of evolutionary divergence or duplication. The 5' end contains widely spaced J558 genes interspersed with 3609 genes and is pseudogene poor. The 3' end contains closely spaced J558 genes, no 3609 genes, and is pseudogene rich. Each occupies a different branch of the phylogenetic tree. Detailed analysis of 500-bp upstream of all functional genes has revealed several conserved binding sites, general and B cell-specific, as well as key differences between families. This complete and definitive assembly of the mouse Igh V region will facilitate detailed study of promoter function and large-scale mechanisms associated with V(D)J recombination including locus contraction and antisense intergenic transcription.

Animals↗
PubMed2002

Cloning and analysis of the genomic DNA sequence of augmenter of liver regeneration from rat.

OBJECTIVE: To search for genomic DNA sequence of the augmenter of liver regeneration (ALR) of rat. METHODS: Polymerase chain reaction (PCR) with specific primers was used to amplify the sequence from the rat genome. RESULTS: A piece of genomic DNA sequence and a piece of pseudogene of rat ALR were identified. The lengths of the gene and pseudogene are 1508 bp and 442 bp, respectively. The ALR gene of rat includes 3 exons and 2 introns. The 442 bp DNA sequence may represent a pseudogene or a ALR-related peptide. Predicted amino acid sequence analysis showed that there were 14 different amino acid residues between the gene and pseudogene. ALR-related peptide is 84 amino acid residues in length and relates closely to ALR protein. CONCLUSION: There might be a multigene family of ALR in rat.

Amino Acid Sequence↗
PubMed2000

Genomic structure of the human p47-phox (NCF1) gene.

The cytosolic factor p47-phox, encoded by the NCF1 gene, is an essential component of the phagocyte NADPH-oxidase system. Upon activation of this multicomponent system, p47-phox translocates to the membrane and participates in the electron transfer from NADPH to molecular oxygen. A deficiency or absence of p47-phox is the most common autosomal form of chronic granulomatous disease (CGD). We have cloned and characterized the NCF1 gene from four bacteriophage clones, a P1 clone and genomic DNA from normal individuals. The gene is 15,236 base pairs long and includes 11 exons. It is 98.6% homologous in sequence to at least one pseudogene that maps to the same region of chromosome 7q11.23. Slightly more than half (50.37%) of the wild-type NCF1 gene consists of repetitive elements. In particular, the density of Alu sequences is high (1.4 Alu/kb); there are 21 Alu repeats interspersed through 10 introns. These findings are consistent with the observation that recombination events between the wild-type gene and its highly homologous pseudogenes account for the majority of potentially lethal mutations in p47-phox-deficient chronic granulomatous disease. Analysis of 1.96 kb of sequence 5' of the start of translation revealed a high homology (99.6%) between wild-type and pseudogene clones. Characterization of NCF1 establishes a foundation for detailed molecular analysis of p47-phox-deficient CGD patients as well as for the study of the regulation of the NCF1 gene and pseudogenes, both of which are present as full-length transcripts in normal individuals.

5' Untranslated Regions↗
PubMed1994

Studies on the antigenicity and nucleotide sequence of the rabies virus Nishigahara strain, a current seed strain used for dog vaccine production in Japan.

The Nishigahara strain of rabies virus, a current seed strain used for animal vaccine production in Japan, is believed to derived from the original Pasteur strain obtained from Paris in or before 1915. In Japan, the virus was serially passaged through several kinds of animals and cell cultures. Reactions with anti-nucleocapsid protein monoclonal antibodies (MAb-N) indicated the Nishigahara strain had maintained the antigenic profile of the Pasteur virus. Reactions with monoclonal antibodies to the glycoprotein (MAb-G) revealed differences between the Nishigahara strain and the Pasteur strain; however, the Nishigahara strain maintained a closer resemblance to the Pasteur virus than to other Pasteur-related viruses or to rabies strains unrelated to the Pasteur strain. Comparative amino acid sequence analysis of cloned cDNA encoding the G gene confirmed the antigenic differences among these strains and the resemblance of the Nishigahara strain to the original Pasteur strain. Comparative nucleotide sequence analysis of the noncoding pseudogene region (Tordo et al., Proc Natl Acad Sci USA 83, 3914-3918, 1986) revealed different relationships. Unlike the Pasteur strain, which encodes a transcription-terminating signal at the end of the G gene (marking the beginning of the pseudogene), a long G-L intergenic sequence in the Nishigahara strain was connected to the 3' end of the cDNA, and the transcription-terminating signal was present only at the end of, but not before, the pseudogene. These results are not inconsistent with the documented origin of the Nishigahara strain, but the genome structure around the pseudogene region suggests divergence from the Pasteur strain and a closer resemblance to other strains of rabies virus.

Amino Acid Sequence↗
PubMed2004

Polymorphic changes in the KAL1 gene: not all of them should be classified as polymorphisms.

The KAL1 gene has a closely related nonfunctional pseudogene on the Y chromosome; a high degree of X-Y sequence similarity is observed. Some individuals present a T to C substitution at position 1833 (exon 12). Because this nucleotide differs in the X (thymine) and in the Y (cytosine) chromosome, we investigated if this was truly a polymorphism, or if in some cases the Y sequence had been amplified. The complete sequence of exon 12 of KAL1 was analyzed in 11 Kallmann Syndrome (KS) males, in 50 normal males, in 50 normal females, and in 16 patients with Ullrich-Turner Syndrome (UTS). Nucleotide 1833 was found in a heterozygous or a homozygous state in KS, normal males and normal females; UTS patients were always homozygous. Of the 61 males, 17 were heterozygous, while 11 were TT and 33 were CC. With these observations we can not assure whether these patients present a "real" polymorphism. Besides, all males were heterozygous in nucleotides 1678, 1694, 1699, 1708 and 1825, whilst females were homozygous; and in these positions, KAL1 also differs from its pseudogene. These results indicate that we are identifying the X and the Y nucleotide and these variants are not polymorphisms. Sequence variations may be pseudogene products rather than true polymorphisms, so we should always determine if the position where the variation is located differs between KAL1 and its pseudogene, because it has been suggested that the presence of various polymorphisms in affected individuals could be the cause of KS.

Adult↗
PubMed1984

The eta-globin gene. Its long evolutionary history in the beta-globin gene family of mammals.

In phylogenetic reconstructions by the parsimony method, utilizing 62 sequenced globin genes and pseudogenes (including 34 of the beta-globin gene family from eutherian orders Primates, Lagomorpha, Artiodactyla and Rodentia), the branch of primate psi beta pseudogenes and the goat embryonically expressed epsilon II gene group monophyletically together as orthologues of a common ancestral gene (labelled eta) distinct from orthologues of epsilon, gamma, delta and beta. This primate psi eta-goat eta branch is cladistically closer to epsilon and gamma than to delta and beta branches. In each eutherian order gene conversions replaced portions of delta by beta sequences, whereas in descent of Primates epsilon, gamma and eta mostly retained their separate ancient identities predating the radiation of Eutheria in all their exons and non-coding regions. The loci of the ancestral beta-globin gene cluster in basal eutherians and proto-primates, as deduced from beta-clusters representing the four eutherian orders, were linked 5'-epsilon-gamma-eta-delta-beta-3' with epsilon, gamma and eta being embryonically expressed genes, and delta and beta ontogenetically later expressed genes. Through deletions gamma was lost in artiodactyl evolution, eta in lagomorph and rodent evolution, and all DNA between exon 2 3' boundaries of eta and delta in prosimian lemuriform evolution (lemur having the hybrid pseudogene psi eta delta). Simian primates retained intact the five loci of the ancestral cluster. Not only did eta, after it became a pseudogene in the basal primates, persist intact in descent to present-day simians but in the line to hominoids it evolved during the last 40 million years at the decelerated rate of 1 X 10(-9) substitutions/site per year which is one-fifth the expected neutral rate. The possibility is suggested that the psi eta locus situated between fetal and adult chromosomal domains of the simian beta-globin gene cluster might play some role in a mechanism for ontogenetic switches of globin gene expression. However, not enough sequence data on genes and intergenic regions in DNA of species of primates and other mammals as yet exist to know if the slow rate of 1 X 10(-9) reflects the rate of a conserved functional gene or primarily reflects a decelerated neutral rate of hominoid DNA evolution, conceivably from enhanced DNA repair and longer generation times in hominoids. The further possibility is raised that gene correction (repair of damaged DNA that prevents emergence of new alleles) and gene conversion both more often involve strand copying of conserved than of rapidly evolving DNA.

Animals↗
PubMed2002

Amplification of multiple copies of mitochondrial Cytochrome b gene fragments in the Australian freshwater crayfish, Cherax destructor Clark (Parastacidae: Decapoda).

Non-coding copies of fragments of the mitochondrial genome translocated to the nucleus or pseudogenes are being found with increasing frequency in a diversity of organisms. As part of a study to evaluate the utility of a range of mitochondrial gene regions for population genetic and systematic studies of the Australian freshwater crayfish, Cherax destructor (the yabby), we report the first detection of Cytochrome b (Cyt b) pseudogenes in crustaceans. We amplified and sequenced fragments of the mitochondrial Cyt b gene from 14 individuals of C. destructor using polymerase chain reaction (PCR) with primers designed from conserved regions of Penaeus monodon and Drosophila melanogaster mitochondrial genomes. The phylogenetic tree produced from the amplified fragments using these primers showed a very different topology to the trees obtained from sequences from three other mitochondrial genes, suggesting one or more nuclear pseudogenes have been amplified. Supporting this conclusion, two highly divergent sequences were isolated from each of two single individuals, and a 2 base pair (bp) deletion in one sequence was observed. There was no evidence to support inadvertent amplification of parasite DNA or contamination of samples from other sources. These results add to other recent observations of pseudogenes suggesting the frequent transfer of mitochondrial DNA (mtDNA) genes to the nucleus and reinforces the necessity of great care in interpreting PCR-generated Cyt b sequences used in population or evolutionary studies in freshwater crayfish and crustaceans more generally.

Animals↗
PubMed1984

Chromosomal organization of the human dihydrofolate reductase genes: dispersion, selective amplification, and a novel form of polymorphism.

The human dihydrofolate reductase (DHFR; tetrahydrofolate dehydrogenase; 5,6,7,8-tetrahydrofolate: NADP+ oxidoreductase, EC 1.5.1.3) gene family includes a functional gene (hDHFR) and at least four intronless genes. Three intronless genes (hDHFR-psi 2, hDHFR-psi 3, and hDHFR-psi 4) are identifiable as pseudogenes because of DNA sequence divergence from the functional gene with introns, while one intronless gene (hDHFR-psi 1) is completely homologous to the coding sequences of the functional gene. Analysis of genomic DNA from two panels of somatic human-rodent cell hybrids with specific molecular probes provide insight into the chromosomal organization and assignment of these genes. The five genes are dispersed in that each one is found on a different chromosome. The functional gene hDHFR has been assigned to chromosome 5, and one pseudogene (hDHFR-psi 4), to chromosome 3. In a human cell line (HeLa) that was selected for methotrexate resistance, the functional locus became amplified, while there was no amplification of the four intronless pseudogenes. hDHFR-psi 1 was found to be present in DNA of some individuals and absent from DNA of others, consistent with a recent evolutionary origin of this gene originally suggested by its sequence identity to the coding portions of the functional gene. The presence or absence of this intronless pseudogene represents a previously unreported form of DNA polymorphism.

Animals↗
PubMed1998

Characterization of a 190-kilobase pair domain of human type I hair keratin genes.

Polymerase chain reaction-based screening of an arrayed human P1 artificial chromosome (PAC) library using primer pairs specific for the human type I hair keratins hHa3-II or hHa6, led to the isolation of two PAC clones, which covered 190 kilobase pairs (kbp) of genomic DNA and contained nine human type I hair keratin genes, one transcribed hair keratin pseudogene, as well as one orphan exon. The hair keratin genes are 4-7 kbp in size, exhibit intergenic distances of 5-8 kbp, and display the same direction of transcription. With one exception, all hair keratin genes are organized into 7 exons and 6 positionally conserved introns. On the basis of sequence homologies, the genes can be grouped into three subclusters of tandemly arranged genes. One subcluster harbors the highly related genes hHa1, hHa3-I, hHa3-II, and hHa4. A second subcluster of highly related genes comprises the novel genes hHa7 and hHa8, as well as pseudogene PsihHaA, while the structurally less related genes hHa6, hHa5, and hHa2 are constituents of the third subcluster. As shown by reverse transcription-polymerase chain reaction, all hair keratin genes, including the pseudogene, are expressed in the human hair follicle. The transcribed pseudogene PsihHaA contains a premature stop codon in exon 4 and exhibits aberrant pre-mRNA splicing. Evolutionary tree construction reveals an early divergence of hair keratin genes from cytokeratin genes, followed by the segregation of the genes into the three subclusters. We suspect that the 190-kbp domain contains the entire complement of human type I hair keratin genes.

Amino Acid Sequence↗
PubMed2006

A genomic fossil reveals key steps in hemoglobin loss by the antarctic icefishes.

Antarctic icefishes are the only vertebrates that do not have hemoglobin and erythrocytes in their blood. These startling phenotypes are associated in several icefish species with deletions of juvenile and adult globin loci, which in red-blooded teleosts are typically composed of tightly linked pairs of alpha- and beta-globin genes. It is unknown if the loss of hemoglobin expression in icefishes was the direct result of such deletions or if other mutational events compromised globin chain synthesis prior to globin gene loss. In this study, we show that 15 of the 16 icefish species have lost the adult beta-globin gene but retain a truncated alpha-globin pseudogene. Surprisingly, a phylogenetically derived icefish species, Neopagetopsis ionah, possesses a complete, but nonfunctional, adult alphabeta-globin complex. This cluster contains 2 distinct beta-globin pseudogenes whose phylogenetic origins span the entire Antarctic notothenioid radiation, consistent with an origin via introgression. Maximum likelihood ancestral state reconstruction supports a scenario of icefish globin gene evolution that involves a single loss of the transcriptionally active adult alphabeta-globin cluster prior to the diversification of the extant species in the clade. Through lineage sorting of ancestral polymorphism, 2 types of alleles became fixed in the clade: 1) the alpha-globin pseudogene of the majority of species and 2) the inactive alphabeta-globin complex of N. ionah. We conclude that the globin pseudogene complex of N. ionah is a "genomic fossil" that reveals key intermediate steps on the pathway to loss of hemoglobin expression by all icefish species.

Animals↗