Comparison of the beta-like globin gene families of rabbits and humans indicates that the gene cluster 5''-epsilon-gamma-delta-beta-3'' predates the mammalian radiation

The members of the rabbit and human beta-like globin gene families have been compared both by a computer-generated dot matrix graphical analysis of each entire gene and by calculating divergences in the coding regions. The rabbit-human gene pairs beta 4-epsilon, beta 3- gamma, psi beta 2-delta, and beta 1-beta were identified as orthologous on the basis of sequence similarities found in flanking and intervening sequences as well as by quantitative divergence calculations. The orthologous genes are in the same order on the chromosome in each species, which suggests that an ancestral family with the arrangement 5'-epsilon-gamma-delta-beta-3' preceded the mammalian radiation. Descendants of ancestral epsilon have diverged more slowly than other beta-like genes and are expressed only in embryonic life. Descendants of ancestral gamma and beta diverged at a higher rate and are expressed at wider range of developmental times. Descendants of delta have undergone nonreciprocal recombination at a high frequency and are often pseudogenes. Paralogous comparisons among the rabbit beta-like globin genes show that the beta 4-beta 3 and psi beta 2-beta 1 pairs are most similar and that beta 4 and beta 3 are more closely related to beta 1 than to psi beta 2. This fits with a branching pattern where the primordial beta split into ancestral epsilon/gamma and delta/beta genes, which later split into epsilon and gamma or delta and beta, respectively. Rabbit genes beta 4 and beta 1 acquired similar 3' untranslated regions after the epsilon/gamma split but prior to the mammalian radiation, presumably via a gene conversion event. The 5' end of beta 2 apparently converted with beta 1 after the radiation, and afterward it became a pseudogene.      

Concerted evolution led to high expression of a prosimian primate delta globin gene locus

The delta globin gene in simian primates is either weakly expressed (in hominoids and New World monkeys) or silent (in Old World monkeys). In prosimian primates, however, an unequal homologous crossover between the psi eta and delta loci of lemurs produced a hybrid psi eta delta pseudogene locus, whereas in tarsier the delta locus encodes a beta-type chain found in 18% of adult tarsier hemoglobin molecules. In the present study, the nucleotide and amino acid sequences of the galago delta and beta globin genes and their encoded peptides were determined, and evidence is provided showing that the galago delta locus encodes a beta-type chain (beta 2) found in 40% of the galago fetal and postnatal hemoglobin molecules, whereas the beta locus encodes the remaining 60% of the beta-type chain (beta 1). Galago beta 1 and beta 2 chains differ from each other by only one amino acid residue. The homology between the galago delta and beta loci extends from 800 base pairs 5' of the proximal CCAAT element to near the end of exon 3 as a result of a recombination event in which beta sequence replaced delta sequence. After this initial recombination event, concerted evolution between the loci continued over their conserved coding, intron 1, and promoter regions but failed to occur between the two loci in their intron 2 and distal 5'-flanking sequences where the two loci have now diverged by 20%. Calculations based on this divergence value and on a rate of noncoding sequence evolution of 4.2 x 10(-9) to 5.5 x 10(-9) substitutions/site/year for the lorisiform lineage to galago yielded a date of 18-24 million years ago for the initial recombination event. The fact that the promoter sequences of the galago delta locus are the same as that of the galago beta locus may account for the high level of expression of the galago delta gene.     

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.      

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.     

Structure and evolution of the horse zeta globin locus

The equine zeta globin gene locus consists of an intact 5' gene and a truncated 3' pseudogene (psi zeta) that has only 5' control sequences and a first exon and intron. Nevertheless, the psi zeta gene has retained almost perfect homology with its neighbour, presumably by gene conversion. The first introns of both zeta and psi zeta genes contain a number of degenerate tandem repeats of a 14 base-pair sequence that has been found in the zeta genes of goats and humans and that is related to a family of human minisatellite sequences. Comparisons of sequences flanking the zeta and psi zeta genes reveal areas of considerable interspecies homology, which can be explained by a zeta gene duplication that pre-dated the mammalian radiation.     

The primate psi beta 1 gene. An ancient beta-globin pseudogene

The human beta-globin gene cluster contains five functional genes plus a single pseudogene termed psi beta 1. Hybridization and comparative sequence analysis show that this pseudogene is not the product of a recent gene duplication, but is ancient and has been maintained in all major primate groups ranging from prosimians to anthropoids, at the same position as in man, between gamma- and delta-globin genes. In the lemur, a prosimian, the central exons of the psi beta 1 and delta-globin genes have undergone an unequal exchange, which has resulted in a contraction of the beta-globin gene cluster and the formation of a Lepore-type psi beta 1-delta globin pseudogene. Comparisons of defects shared by prosimian, New World monkey and human psi beta 1 sequences suggest that the ancestral primate gene was probably a pseudogene with an abnormal initiation codon but few if any additional defects, and that most contemporary pseudogene defects were accumulated relatively recently by slow neutral drift. We suggest that psi beta 1 arose early in primate evolution by silencing of a pre-existing discrete functional gene, and show that psi beta 1-related sequences are also present in other mammalian orders. In view of the antiquity of psi beta 1-related sequences, we propose that this gene be renamed the eta-globin gene.     

The phylogenetic history of New World monkey beta globin reveals a platyrrhine beta to delta gene conversion in the atelid ancestry

Orthologues of the beta globin gene locus from 10 New World monkey species were sequenced and aligned against available beta and delta globin sequences from rabbit and other primates. Where needed, additional primate sequencing was performed. Phylogenetic analysis identified a beta to delta conversion in the stem of the Anthropoidea, stretching from the 3' part of the proximal promotor to the 5' start of intron 2, consistent with earlier findings. No further conversion appeared to have occurred in the descent of the catarrhines. Within the New World monkey lineage that led to spider monkey and other atelids, another shorter gene conversion was found, spanning adjacent parts of exon 1 and intron 1. The analysis also confirmed that galago beta had replaced galago delta, that an earlier loriform-specific gene conversion extended over intron 2, and that gene conversion throughout the main gene conversion region occurred in the tarsiiform lineage. Platyrrhine phylogenetic relationships were investigated with beta sequences restricted to those that were not involved in gene conversions. This phylogeny generally agreed with results from other nuclear genes. The one exception was that the beta sequences did not place the callitrichine clade within the Cebidae but weakly joined the callitrichine and atelid clades.     

Isolation of the chicken beta-globin gene and a linked embryonic beta-like globin gene from a chicken DNA recombinant library.

A library of random chicken DNA fragments, 15-22 kb long, has been prepared in the vector lambda Charon 4A. This library was screened with combined adult and embryonic globin cDNA, and several independent globin gene-containing recombinants were isolated. One of these recombinants, lambda Chicken beta-globin 1 (lambda C beta G1), contains the adult chicken beta-globin gene and a closely linked embryonic beta-like globin gene. Both genes are transcribed in the same direction with the adult gene located 5' to the embryonic gene. Electron microscopic visualization of R loop structures generated by hybridization of globin RNA to lambda C beta G1 demonstrates that both globin genes contain major intervening sequences about 800 bp long, similar to those present in mammalian beta-globin genes. The adult beta-globin gene also contains a minor (approximately 100 bp long) intervening sequence analogous to the one observed in mammalian beta-globin genes. Restriction enzyme analysis of the adult beta-globin gene on lambda C beta G1 is consistent with the hypothesis that its two intervening sequences occur in the same positions with respect to the beta-globin amino acid sequence as do the corresponding mammalian intervening sequences.     

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.     

Molecular cloning and characterization of two sets of alpha-theta genes in the rat alpha-like globin gene cluster

The highly heterogeneous rat hemoglobin system was investigated at the gene level. Two regions of the alpha-like globin gene cluster from a Wistar rat were isolated. Four lambda Dash recombinant clones carrying rat alpha-like globin genes were localized on two distinct gene regions. A region of approximately 16kb was found to contain the 5'-IIalpha1-psi theta 1-3' loci, and another of approximately 24kb the 5'-IIalpha2-psi theta2-psiI alpha3-3' loci. Both IIalpha1 and IIalpha2 are considered to be active, coding the IIalpha-globin chain. The nt sequences of IIalpha1 and IIalpha2 are identical except for six nt in the non-coding region. The psiI alpha3 locus is a truncated pseudogene. The putative promoter region of an alpha-like globin gene is joined directly to the third exon, homologous to that of Ialpha-globin cDNA. psi theta1 and psi theta2 are also pseudogenes, as evidenced by several deletions located in the protein-coding regions of these loci. The psi theta1 and psi theta2 loci exhibit extensive homology, but the restriction maps of these genes and their flanking regions differ considerably. Genomic Southern blot analyses of the total liver DNA from six rats showed the existence of three theta-globin-related genes, including psi theta1 and psi theta2. These results indicate that the two gene regions investigated are not allelic variants, but may be generated by block duplication. This is the first report of the existence of rodent theta-globin genes.     

Assignment of orthologous relationships among mammalian alpha-globin genes by examining flanking regions reveals a rapid rate of evolution

In order to study the relationships among mammalian alpha-globin genes, we have determined the sequence of the 3' flanking region of the human alpha 1 globin gene and have made pairwise comparisons between sequenced alpha-globin genes. The flanking regions were examined in detail because sequence matches in these regions could be interpreted with the least complication from the gene duplications and conversions that have occurred frequently in mammalian alpha-like globin gene clusters. We found good matches between the flanking regions of human alpha 1 and rabbit alpha 1, human psi alpha 1 and goat I alpha, human alpha 2 and goat II alpha, and horse alpha 1 and goat II alpha. These matches were used to align the alpha-globin genes in gene clusters from different mammals. This alignment shows that genes at equivalent positions in the gene clusters of different mammals can be functional or nonfunctional, depending on whether they corrected against a functional alpha-globin gene in recent evolutionary history. The number of alpha-globin genes (including pseudogenes) appears to differ among species, although highly divergent pseudogenes may not have been detected in all species examined. Although matching sequences could be found in interspecies comparisons of the flanking regions of alpha- globin genes, these matches are not as extensive as those found in the flanking regions of mammalian beta-like globin genes. This observation suggests that the noncoding sequences in the mammalian alpha-globin gene clusters are evolving at a faster rate than those in the beta-like globin gene clusters. The proposed faster rate of evolution fits with the poor conservation of the genetic linkage map around alpha-globin gene clusters when compared to that of the beta-like globin gene clusters. Analysis of the 3' flanking regions of alpha-globin genes has revealed a conserved sequence approximately 100-150 bp 3' to the polyadenylation site; this sequence may be involved in the expression or regulation of alpha-globin genes.      

Complete nucleotide sequence of the rabbit beta-like globin gene cluster. Analysis of intergenic sequences and comparison with the human beta-like globin gene cluster

The nucleotide sequence of the entire beta-like globin gene cluster of rabbits has been determined. This sequence of a continuous stretch of 44.5 x 10(3) base-pairs (bp) starts about 6 x 10(3) bp upstream from epsilon (the 5'-most gene) and ends about 12 x 10(3) bp downstream from beta (the 3'-most gene). Analysis of the sequence reveals that: (1) the sequence is relatively A + T rich (about 60%); (2) regions with high G + C content are associated with OcC repeats, a short interspersed repeated DNA in rabbits; (3) the distribution of polypurines, polypyrimidines and alternating purine/pyrimidine tracts is not random within the cluster; (4) most open reading frames are associated with known globin coding regions, OcC repeats or long interspersed repeats (L1 repeats); (5) the most prominent open reading frames are found in the L1 repeats; (6) different strand asymmetries in base composition are associated with embyronic and adult genes as well as the tandem L1 repeats at the 3' end of the cluster; and (7) essentially all the repeats appear to have been inserted by a transposon mechanism. A comparison of the sequence with itself by a dot-plot analysis has revealed nine new members of the OcC family of repeats in addition to the six previously reported. The OcC repeats tend to be clustered, particularly in the epsilon-gamma and gamma-psi delta intergenic regions. Dot-plot comparisons between the rabbit and the human clusters have revealed extensive sequence matches. Homology starts about 6 x 10(3) bp 5' to epsilon or as far upstream as the rabbit sequence is available. It continues throughout the entire cluster and stops about 0.7 x 10(3) bp 3' to beta, at which point several repeats have inserted in both rabbits and humans. Throughout the gene cluster, the homology is interrupted mainly by insertions or deletions in either the rabbit or the human genome. Almost all of the insertions are of known short or long repeated DNAs. The positions of the insertions are different in the two gene clusters, which indicates that both short and long repeats have been transposing throughout the genome for the time since the mammalian radiation. An alignment of rabbit and human sequences allows the calculation of the substitution rate around epsilon. Sequences far removed from the gene are evolving at a rate equivalent to the pseudogene rate, although some short regions show an apparently higher rate.(ABSTRACT TRUNCATED AT 400 WORDS)     

Sequence and comparative analysis of the rabbit alpha-like globin gene cluster reveals a rapid mode of evolution in a G + C-rich region of mammalian genomes

A sequence of 10,621 base-pairs from the alpha-like globin gene cluster of rabbit has been determined. It includes the sequence of gene zeta 1 (a pseudogene for the rabbit embryonic zeta-globin), the functional rabbit alpha-globin gene, and the theta 1 pseudogene, along with the sequences of eight C repeats (short interspersed repeats in rabbit) and a J sequence implicated in recombination. The region is quite G + C-rich (62%) and contains two CpG islands. As expected for a very G + C-rich region, it has an abundance of open reading frames, but few of the long open reading frames are associated with the coding regions of genes. Alignments between the sequences of the rabbit and human alpha-like globin gene clusters reveal matches primarily in the immediate vicinity of genes and CpG islands, while the intergenic regions of these gene clusters have many fewer matches than are seen between the beta-like globin gene clusters of these two species. Furthermore, the non-coding sequences in this portion of the rabbit alpha-like globin gene cluster are shorter than in human, indicating a strong tendency either for sequence contraction in the rabbit gene cluster or for expansion in the human gene cluster. Thus, the intergenic regions of the alpha-like globin gene clusters have evolved in a relatively fast mode since the mammalian radiation, but not exclusively by nucleotide substitution. Despite this rapid mode of evolution, some strong matches are found 5' to the start sites of the human and rabbit alpha genes, perhaps indicating conservation of a regulatory element. The rabbit J sequence is over 1000 base-pairs long; it contains a C repeat at its 5' end and an internal region of homology to the 3'-untranslated region of the alpha-globin gene. Part of the rabbit J sequence matches with sequences within the X homology block in human. Both of these regions have been implicated as hot-spots for recombination, hence the matching sequences are good candidates for such a function. All the interspersed repeats within both gene clusters are retroposon SINEs that appear to have inserted independently in the rabbit and human lineages.