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.     

Molecular cloning and characterization of the human beta-like globin gene cluster

The genes encoding human embryonic (epsilon), fetal (G gamma, A gamma) and adult (delta, beta) beta-like globin polypeptides were isolated as a set of overlapping cloned DNA fragments from bacteriophage lambda libraries of high molecular weight (15-20 kb) chromosomal DNA. The 65 kb of DNA represented in these overlapping clones contains the genes for all five beta-like polypeptides, including the embryonic epsilon-globin gene, for which the chromosomal location was previously unknown. All five genes are transcribed from the same DNA strand and are arranged in the order 5'-epsilon-(13.3 kb)-G gamma-(3.5 kb)-A gamma-(13.9 kb)-delta-(5.4 kb)-beta-3'. Thus the genes are positioned on the chromosome in the order of their expression during development. In addition to the five known beta-like globin genes, we have detected two other beta-like globin sequences which do not correspond to known polypeptides. One of these sequences has been mapped to the A gamma-delta intergenic region while the other is located 6-9 kb 5' to the epsilon gene. Cross hybridization experiments between the intergenic sequences of the gene cluster have revealed a nonglobin repeat sequence (*) which is interspersed with the globin genes in the following manner: 5'-**epsilon-*G gamma-A gamma*-**delta-beta*-3'. Fine structure mapping of the region located 5' to the delta-globin gene revealed two repeats with a maximum size of 400 bp, which are separated by approximately 700 bp of DNA not repeated within the cluster. Preliminary experiments indicate that this repeat family is also repeated many times in the human genome.     

Structure and organization of the bovine beta-globin genes

Genomic clones spanning the entire cow beta-globin gene locus have been isolated and characterized. These clones demonstrate that the linkage of embryonic-like (epsilon) genes and pseudogenes (psi) to the previously described fetal (gamma) and adult (beta) genes is as follows: 5'-epsilon 3-epsilon 4-psi 3-beta-epsilon 1-epsilon 2-psi 1- psi 2-gamma-3'. Present data indicate that, like that of the goat, the fetal and adult genes arose via block duplication of an ancestral four- gene set: epsilon-epsilon-psi-beta. This duplication event preceded the divergence of cows and goats, which occurred greater than or equal to 18-20 Myr ago. However, cows do not have the additional four-gene block containing a preadult/stress globin gene (beta C). Furthermore, the cow fetal cluster contains an extra beta-like pseudogene, which apparently arose by a small-scale duplication. The fixation of this duplication may indicate a possible evolutionary role for pseudogenes.      

New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents

Species differences in the size or membership composition of multigene families can be attributed to lineage-specific additions of new genes via duplication, losses of genes via deletion or inactivation, and the creation of chimeric genes via domain shuffling or gene fusion. In principle, it should be possible to infer the recombinational pathways responsible for each of these different types of genomic change by conducting detailed comparative analyses of genomic sequence data. Here, we report an attempt to unravel the complex evolutionary history of the beta-globin gene family in a taxonomically diverse set of rodent species. The main objectives were: 1) to characterize the genomic structure of the beta-globin gene cluster of rodents; 2) to assign orthologous and paralogous relationships among duplicate copies of beta-like globin genes; and 3) to infer the specific recombinational pathways responsible for gene duplications, gene deletions, and the creation of chimeric fusion genes. Results of our comparative genomic analyses revealed that variation in gene family size among rodent species is mainly attributable to the differential gain and loss of later expressed beta-globin genes via unequal crossing-over. However, two distinct recombinational mechanisms were implicated in the creation of chimeric fusion genes. In muroid rodents, a chimeric gamma/epsilon fusion gene was created by unequal crossing-over between the embryonic epsilon- and gamma-globin genes. Interestingly, this gamma/epsilon fusion gene was generated in the same fashion as the "anti-Lepore" 5'-delta-(beta/delta)-beta-3' duplication mutant in humans (the reciprocal exchange product of the pathological hemoglobin Lepore deletion mutant). By contrast, in the house mouse, Mus musculus, a chimeric beta/delta fusion pseudogene was created by a beta-globin --> delta-globin gene conversion event. Although the gamma/epsilon and beta/delta fusion genes share a similar chimeric gene structure, they originated via completely different recombinational pathways.     

A molecular and evolutionary study of the beta-globin gene family of the Australian marsupial Sminthopsis crassicaudata

Beta-globin gene families in eutherians (placental mammals) consist of a set of four or more developmentally regulated genes which are closely linked and, in general, arranged in the order 5'-embryonic/fetal genes- adult genes-3'. This cluster of genes is proposed to have arisen by tandem duplication of ancestral beta-globin genes, with the first duplication occurring 200 to 155 MYBP just prior to a period in mammalian evolution when eutherians and marsupials diverged from a common ancestor. In this paper we trace the evolutionary history of the beta-globin gene family back to the origins of these mammals by molecular characterization of the beta-globin gene family of the Australian marsupial Sminthopsis crassicaudata. Using Southern and restriction analysis of total genomic DNA and bacteriophage clones of beta-like globin genes, we provide evidence that just two functional beta-like globin genes exist in this marsupial, including one embryonic- expressed gene (S.c-epsilon) and one adult-expressed gene (S.c-beta), linked in the order 5'-epsilon-beta-3'. The entire DNA sequence of the adult beta-globin gene is reported and shown to be orthologous to the adult beta-globin genes of the North American marsupial Didelphis virginiana and eutherian mammals. These results, together with results from a phylogenetic analysis of mammalian beta-like globin genes, confirm the hypothesis that a two-gene cluster, containing an embryonic- and an adult-expressed beta-like globin gene, existed in the most recent common ancester of marsupials and eutherians. Northern analysis of total RNA isolated from embryos and neonatals indicates that a switch from embryonic to adult gene expression occurs at the time of birth, coinciding with the transfer of the marsupial from a uterus to a pouch environment.      

The beta globin gene cluster of the prosimian primate Galago crassicaudatus: nucleotide sequence determination of the 41-kb cluster and comparative sequence analyses.

The nucleotide sequence of the beta globin gene cluster of the prosimian Galago crassicaudatus has been determined. A total sequence spanning 41,101 bp contains and links together previously published sequences of the five galago beta-like globin genes (5'-epsilon-gamma-psi eta-delta-beta-3'). A computer-aided search for middle interspersed repetitive sequences identified 10 LINE (L1) elements, including a 5' truncated repeat that is orthologous to the full-length L1 element found in the human epsilon-gamma intergenic region. SINE elements that were identified included one Alu type I repeat, four Alu type II repeats, and two methionine tRNA-derived Monomer (type III) elements. Alu type II and Monomer sequences are unique to the galago genome. Structural analyses of the cluster sequence reveals that it is relatively A+T rich (about 62%) and regions with high G+C content are associated primarily with globin coding regions. Comparative analyses with the beta globin cluster sequences of human, rabbit, and mouse reveal extensive sequence homologies in their genic regions, but only human, galago, and rabbit sequences share extensive intergenic sequence homologies. Divergence analyses of aligned intergenic and flanking sequences from orthologous human, galago, and rabbit sequences show a gradation in the rate of nucleotide sequence evolution along the cluster where sequences 5' of the epsilon globin gene region show the least sequence divergence and sequences just 5' of the beta globin gene region show the greatest sequence divergence.     

Rabbit globin pseudogene psi beta 2 is a hybrid of delta- and beta-globin gene sequences

The evolutionary history of the rabbit globin pseudogene psi beta 2 was studied by completing its nucleotide sequence and aligning the sequence with that of the rabbit adult globin gene beta 1 and the human minor adult globin gene delta. The 5' flanking region and exon 1 of psi beta 2 were most similar to rabbit beta 1, but the large intervening sequence and the 3' untranslated region were most similar to human delta. Intron 1 and exon 2 were equally similar to both delta and beta 1. This pattern indicates that psi beta 2 was originally a delta-like gene that acquired the 5' portion of gene beta 1 by intrachromosomal gene conversion. The presence of a delta-globin gene sequence in both rabbits and humans shows that it is an ancient gene, predating the mammalian radiation that occurred over 85 Myr ago. Delta has shown a pronounced tendency to be altered in its 5' end during the course of mammalian evolution. Quantitative divergence analysis shows that the ancestor to rabbit psi beta 2 was active until 20-30 Myr ago, during which time the lagomorph beta-globin gene family apparently functioned without a pseudogene.      

Use of long sequence alignments to study the evolution and regulation of mammalian globin gene clusters

The determination of long segments of DNA sequences encompassing the beta- and alpha-globin gene clusters has provided an unprecedented data base for analysis of genome evolution and regulation of gene clusters. A newly developed computer tool kit generates local alignments between such long sequences in a space-efficient manner, helps the user analyze the alignments effectively, and finds consistently aligning blocks of sequences in multiple pairwise comparisons. Such sequence analyses among the beta-like globin gene clusters of human, galago, rabbit, and mouse have revealed the general patterns of evolution of this gene cluster. Alignments in the flanking regions are very useful in assigning orthologous relationships. Investigation of such matches between the mouse and human beta-like globin gene clusters has led to a reassessment of some orthologous assignments in mouse and to a revision of the proposed pathway for evolution of this gene cluster. In general, the interspersed repetitive elements have inserted independently, presumably via a retrotransposition mechanism, in the different mammalian lineages. However, some examples of ancient L1 repeats are found, including one between the epsilon- and gamma-globin genes that appears to have been in the ancestral eutherian gene cluster. Prominent matching sequences are found in a long region 5' to the epsilon-globin gene, the locus control region (LCR) that is a positive regulator of the entire gene cluster. Three-way alignments among the human, goat, and rabbit sequences can extend for > or = 3 kb in part of the LCR (DNase hypersensitive site 3), indicating that the cis-acting components of this complex regulatory region cover a long segment of DNA. In contrast to the beta-like globin gene clusters, the alpha-like globin gene clusters of many mammals occur in very G+C-rich isochores and contain prominent CpG islands. The regions between the alpha-like globin genes are evolving faster than the intergenic regions of the beta-like globin gene clusters. The contrasts between the two gene clusters can be attributed to differences in DNA metabolism in the isochore. The proximal control elements of the rabbit alpha-globin gene are located both 5' to and within the gene. All of this region is part of a prominent CpG island that may be acting as an extended, enhancer- independent promoter. One can hypothesize that the analogue to the LCR in the alpha-globin gene cluster may interface with the distinctive alpha-globin promoter in ways different from the interaction between the beta LCR and the promoters of beta-like globin genes.(ABSTRACT TRUNCATED AT 400 WORDS)      

Haplotype analysis of the human beta-globin gene complex using multiple locus specific oligonucleotide probes

Three oligonucleotide probes complementary to specific DNA sequences of the six human globin genes (epsilon, G gamma, A gamma, psi beta, delta, beta) were synthesized. The oligonucleotides were used either singly or in combination as hybridization probes to determine the haplotype of the human beta-globin gene cluster employing the four conventionally used restriction endonucleases HincII, HindIII, AvaII, and BamHI, in addition to HpaI. Polymorphism in the epsilon- and psi beta-genes (HincII) can be simultaneously determined with a single probe mixture. One of the probes complementary to both the psi beta- and gamma-genes is useful for determining both HindIII and HincII polymorphisms. The advantages of these probes relative to conventional cDNA probes are discussed.     

Physical mapping of the globin gene deletion in hereditary persistence of foetal haemoglobin (HPFH).

We have mapped the globin gene region in the DNA of two HPFH patients. In a patient homozygous for the G gamma A gamma type of HPFH at least 24 kb of DNA in the globin gene region has been deleted to remove most of the gamma-delta intergenic region and the delta and beta globin genes. The 5' break point of the deletion is located about 9 kb upstream from the delta globin gene. The 3' break point has not been precisely located but is at least 7 kb past the beta globin gene. DNA from an individual heterozygous for the Greek (A gamma) type of HPFH, however, shows no detectable deletion in the entire gamma delta beta-globin gene region. HPFH, therefore, appears to occur in different molecular forms. These results are discussed in terms of a model for the regulation of globin gene expression in man.     

Nucleotide sequence analysis of the lemur beta-globin gene family: evidence for major rate fluctuations in globin polypeptide evolution

Lemur beta-related globin genes have been isolated and sequenced. Orthology of prosimian and human epsilon-, gamma-, and beta-related globin genes was established by dot-matrix analysis. All of these lemur globin genes potentially encode functional beta-related globin polypeptides, though precisely when the gamma-globin gene is expressed remains unknown. The organization of the 18-kb brown lemur beta-globin gene cluster (5' epsilon-gamma-[psi eta-delta]-beta 3') is consistent with its evolution by contraction via unequal crossing-over from the putative ancestral mammalian beta-globin gene cluster (5' epsilon-gamma- eta-delta-beta 3'). The dwarf lemur nonadult globin genes are arranged as in the brown lemur. Similar levels of synonymous (silent) nucleotide substitutions and noncoding DNA sequence differences have accumulated between species in all of these genes, suggesting a uniform rate of noncoding DNA divergence throughout primate beta-globin gene clusters. These differences are comparable with those observed in the nonfunctional psi eta pseudogene and have therefore accumulated at the presumably maximal neutral rate. In contrast, nonsynonymous (replacement) nucleotide substitutions show a significant heterogeneity in distribution for both the same gene in different lineages and different genes in the same lineage. These major fluctuations in replacement but not silent substitution rates cannot be attributed to changes in mutation rate, suggesting that changes in the rate of globin polypeptide evolution in primates is not governed solely by variable mutation rates.      

Duplication of a four-gene set during the evolution of the goat beta-globin locus produced genes now expressed differentially in development

Genomic clones which link the goat preadult (beta C) and adult (beta A) beta-globin genes have been isolated. These overlapping clones contain a previously unidentified embryonic like globin gene (epsilon III) and establish the following linkage map of eight genes in the goat beta-globin locus: epsilon I-epsilon II-psi beta X-beta C-epsilon III-epsilon IV-psi beta Z-beta A. This linkage map and the nucleotide sequence of the eight genes document a relatively recent duplication of a four-gene set: epsilon-epsilon-psi beta-beta. This duplication produced two genes (beta C and beta A) which are now expressed differentially during development. An embryonic like globin gene located downstream from beta A has also been isolated. The embryonic nature of this gene as well as the adult beta-like sequence of the goat fetal globin gene (gamma) suggest that a duplication of the four-gene set also produced the globin gene now expressed during fetal development.     

Identification by nucleotide sequence analysis of a goat pseudoglobin gene

We have recently prepared a recombinant library of goat genomic DNA and have isolated clones containing th known goat globin genes. These include the alpha, gamma, beta C and beta A genes. In addition to these, another beta-like sequence has been observed. In this communication we report the complete nucleotide sequence of this gene, excluding a portion of the large intervening sequence. Several features suggest that this is a non-functional or pseudoglobin gene. The alterations include a frameshift mutation, substitution of the heme-binding histidines, a mutated termination codon, a change in the GT/AG excision sequence of the 5' end of the first intervening sequence, an AT rich sequence in the 3' untranslated region, and a mutated Hogness-Goldberg box. We conclude that this gene cannot function in the synthesis of globin.