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.      

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.      

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.     

Isolation and nucleotide sequence analysis of the beta-type globin pseudogene from human, gorilla and chimpanzee

The beta-globin gene cluster of human, gorilla and chimpanzee contain the same number and organization of beta-type globin genes: 5'-epsilon (embryonic)-G gamma and A gamma (fetal)-psi beta (inactive)-delta and beta (adult)-3'. We have isolated the psi beta-globin gene regions from the three species and determined their nucleotide sequences. These three pseudogenes each share the same substitutions in the initiator codon (ATG----GTA), a substitution in codon 15 which generates a termination signal TGG----TGA, nucleotide deletion in codon 20 and the resulting frame shift which yields many termination signals in exons 2 and 3. The basic structure of these psi beta-globin genes, however, remains consistent with that found for functional beta-globin genes: their coding regions are split by two introns, IVS 1 (which splits codon 30, 121 base-pairs in length) and IVS 2 (which splits codon 104, 840 to 844 base-pairs in length). These introns retain the normal splice junctions found in other eukaryotic split genes. The three hominoid psi beta-globin genes show a high degree of sequence correspondence, with the number of differences found among them being only about one-third of that predicted for DNA sites evolving at the neutral rate (i.e. for sites evolving in the absence of purifying selection). Thus, there appears to be a deceleration in the rate of evolution of the psi beta-globin locus in higher primates.     

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.     

Structure of the goat psi beta y beta-globin pseudogene. Analysis of goat pseudogene evolutionary patterns

The 12-member beta-globin gene locus of the goat contains three beta(adult)-type pseudogenes, one in each of three four-gene subsets of the locus. We have determined the complete nucleotide sequence of psi beta y, the pseudogene present in the most downstream four-gene subset, which also contains the functional fetal gene, beta F. psi beta y contains, throughout its length, numerous incapacitating mutations in common with the previously sequenced goat psi beta x and psi beta z pseudogenes consistent with the model that all were descended from a common pseudogene ancestor which became defective prior to the expansion of the beta-globin locus in the goat lineage. Evolutionary analysis of the psi beta y sequence in comparison to psi beta x and psi beta z provides evidence that nucleotide substitutions were fixed in a random manner within these pseudogenes with respect to polarity, coding versus non-coding regions, and replacement sites versus silent sites. However, substitutions appear to have accumulated asymmetrically between different pseudogenes in a manner that provides evidence for partial gene conversion. Moreover, the presence of deletions in goat psi beta y, which are also observed in the cow pseudogene psi 2, but not in the cow psi 1 pseudogene, indicate that goat psi beta y and cow psi 2 are orthologous but cow psi 1 actually arose prior to the goat/cow divergence. The authentic goat orthologue to cow psi 1 temporarily existed in the goat lineage but was deleted, probably prior to the divergence of goats and sheep.     

Reexamination of the African hominoid trichotomy with additional sequences from the primate beta-globin gene cluster.

Additional DNA sequence information from a range of primates, including 13.7 kb from pygmy chimpanzee (Pan paniscus), was added to data sets of beta-globin gene cluster sequence alignments that span the gamma 1, gamma 2, and psi eta loci and their flanking and intergenic regions. This enlarged body of data was used to address the issue of whether the ancestral separations of gorilla, chimpanzee, and human lineages resulted from only one trichotomous branching or from two dichotomous branching events. The degree of divergence, corrected for superimposed substitutions, seen in the beta-globin gene cluster between human alleles is about a third to a half that observed between two species of chimpanzee and about a fourth that between human and chimpanzee. The divergence either between chimpanzee and gorilla or between human and gorilla is slightly greater than that between human and chimpanzee, suggesting that the ancestral separations resulted from two closely spaced dichotomous branchings. Maximum parsimony analysis further strengthened the evidence that humans and chimpanzees share the longest common ancestry. Support for this human-chimpanzee clade is statistically significant at P = 0.002 over a human-gorilla clade or a chimpanzee-gorilla clade. An analysis of expected and observed homoplasy revealed that the number of sequence changes uniquely shared by human and chimpanzee lineages is too large to be attributed to homoplasy. Molecular clock calculations that accommodated lineage variations in rates of molecular evolution yielded hominoid branching times that ranged from 17-19 million years ago (MYA) for the separation of gibbon from the other hominoids to 5-7 MYA for the separation of chimpanzees from humans. Based on the relatively late dates and mounting corroborative evidence from unlinked nuclear genes and mitochondrial DNA for the close sister grouping of humans and chimpanzees, a cladistic classification would place all apes and humans in the same family. Within this family, gibbons would be placed in one subfamily and all other extant hominoids in another subfamily. The later subfamily would be divided into a tribe for orangutans and another tribe for gorillas, chimpanzees, and humans. Finally, gorillas would be placed in one subtribe with chimpanzees and humans in another, although this last division is not as strongly supported as the other divisions.     

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.      

Immunoglobulin CH gene family in hominoids and its evolutionary history.

The organization of the human immunoglobulin CH gene suggests that a gene duplication involving the C gamma-C gamma-C epsilon-C alpha region has occurred during evolution. We previously showed that both chimpanzee and gorilla have two 5'-C epsilon-C alpha-3', as in human, and that orangutan, gibbon, and Old World monkeys have one C epsilon gene and one, two, and one C alpha gene(s), respectively. In addition to these clustered CH genes, there is one processed C epsilon pseudogene in each species. The present study revealed that orangutan and crab-eating macaque (an Old World monkey) both have one 5'-C epsilon-C alpha-3' and that gibbon has two 5'-C epsilon-C alpha-3', one C epsilon gene of which is completely deleted. By Southern analysis, the number of C gamma genes in all the nonhuman hominoids was estimated to be four to five, as in human, in comparison with two for crab-eating macaque. The C mu and C delta genes were estimated to be present as single copies in both hominoids and crab-eating macaque. Furthermore, it was proved that there are two copies of the C epsilon 5'-flanking region in both the orangutan and the gibbon genomes. These results show that gene duplication including the C gamma-C gamma-C epsilon-C alpha genes occurred in the common ancestor of hominoids and that subsequent deletion of the C epsilon gene (in orangutan, including one of the C alpha genes) occurred independently in each hominoid species.     

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.     

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.     

Evolution and molecular characterization of a beta-globin gene from the Australian Echidna Tachyglossus aculeatus (Monotremata).

Coinciding with a period in evolution when monotremes, marsupials, and eutherians diverged from a common ancestor, a proto-beta-globin gene duplicated, producing the progenitors of mammalian embryonic and adult beta-like globin genes. To determine whether monotremes contain orthologues of these genes and to further investigate the evolutionary relationships of monotremes, marsupials, and eutherians, we have determined the complete DNA sequence of an echidna (Tachyglossus aculeatus) beta-like globin gene. Conceptual translation of the gene and sequence comparisons with eutherian and marsupial beta-like globin genes and echidna adult beta-globin indicate that the gene is adult expressed. Phylogenetic analyses do not clearly resolve the branching pattern of mammalian beta-like globin gene lineages and it is therefore uncertain whether monotremes have orthologues of the embryonic beta-like globin genes of marsupials and eutherians. Four models are proposed that provide a framework for interpreting further studies on the evolution of beta-like globin genes in the context of the evolution of monotremes, marsupials, and eutherians.     

Evolutionary rate of immunoglobulin alpha noncoding region is greater in hominoids than in Old World monkeys.

Recent studies on the molecular evolution of primates show that the evolutionary rate among hominoids is considerably slower than that among nonhominoid primates. However, this observation at the nucleotide-sequence level is restricted to the beta-globin family region. In this study, we sequenced orthologous immunoglobulin alpha (C alpha) genes of chimpanzee, gorilla, orangutan, and crab-eating macaque (an Old World monkey) and compared them with that of the human by using noncoding regions for analysis. Since significant differences in rates among hominoids were not found by using the relative rate test, we evaluated the ratio (R) of the evolutionary distance between Old World monkey and human to the distance between orangutan and human. The R value (1.12) for the C alpha gene was much smaller than the expected value (1.38-2.33), showing that the nucleotide substitution rate (= mutation rate per year under selective neutrality) of the C alpha gene is greater in the human lineage than in the Old World monkey lineage. We also did a similar analysis for the gamma 1-, gamma 2-, psi eta-, and delta-globin genes and found a considerable heterogeneity (1.12-2.37) among the R values, including that for the C alpha gene. This indicates that the hominoid slowdown of the evolutionary rate is not a universal phenomenon in primate evolution.     

Embryonic epsilon and gamma globin genes of a prosimian primate (Galago crassicaudatus). Nucleotide and amino acid sequences, developmental regulation and phylogenetic footprints

Sequence analysis of epsilon and gamma genes and encoded globins and high-pressure liquid chromatography analysis of globin compositions in blood hemolysates obtained from embryos, fetuses and adults show that the prosimian primate Galago crassicaudatus expresses its epsilon and gamma genes only embryonically. Since rabbit, mouse and galago all have embryonic gamma genes but simian primates have fetal gamma genes, we conclude that gamma E evolved into gamma F in stem-simians. An elevated non-synonymous substitution rate characterizes this transition. The alignment of epsilon and gamma nucleotide sequences and the parsimoniously reconstructed evolutionary history of these sequences identify several anciently conserved cis-regulatory elements (phylogenetic footprints) important for gamma expression in primates and also cis-mutations which may have been involved in the recruitment of the gamma gene to a fetal program in simian primates.     

Concerted evolution of the cow epsilon 2 and epsilon 4 beta-globin genes

The nucleotide sequences of the cow epsilon 2 and epsilon 4 globin genes were determined. The sequences were 95% identical. These genes arose via a four-gene block duplication that also gave rise to the bovine fetal (gamma) and adult (beta) genes. Their deduced amino acid sequences are unlike any previously reported fetal or adult globins; rather, comparison to other mammalian globin genes indicates that they are embryonic in nature. The sequence data indicate that these two genes have converted each other during evolution. Pairwise comparison to the corresponding goat genes shows greater similarity between paralogues than between more directly related orthologues. This is in direct contrast to the situation between the cow and goat fetal and adult genes. These observations suggest that the frequency of DNA conversion or the fixation of conversion events may vary in different locations of the cow beta-globin cluster.      

Structural and functional analysis of the goat epsilon-globin genes

Since none of the vertebrate beta-globin loci studied to date has more than two functional embryonic beta-like globin genes, it would be unique if all six goat embryonic beta-globin genes were required for its survival. In this study we have asked whether all six embryonic genes in the goat are functional. This question has been addressed by examining the transient expression of these genes in HeLa cells and correlating these results with the sequence information obtained to date. Our studies show that only epsilon I and epsilon II are functional while the remaining four epsilon-globin genes are nonfunctional, i.e., pseudogenes. Interestingly, the two active epsilon-globin genes are located at the 5' end of the locus. While this unusual inactivation pattern may be the result of chance, it could also have resulted because the two duplication events, of the ancestral gene set epsilon-epsilon-psi beta-beta, did not include distally located regulatory region(s) essential for epsilon-globin gene expression. Once separated from the 5'-regulatory sequences the remaining four embryonic genes (epsilon III, epsilon IV, epsilon V and epsilon VI) accumulated mutations and became pseudogenes.     

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.     

Genetic diversity at class II DRB loci of the primate MHC

The evolution of polymorphism at loci encoding the beta-chains of the MHC class II DR Ag was studied in primates by DNA amplification (polymerase chain reaction). Phylogenetic analysis of 63 DRB sequences from the polymorphic second exon (first domain) of nonhuman primates and 53 human sequences indicates the presence of five DRB loci in primates, derived from a DRB1-like ancestral locus over 20 million yr ago. Many of the allelic types at the DRB1 locus predate the divergence of hominoids (5 million yr ago) and some (DR4, DR3, 5, 6) predate the divergence of Old world monkeys and hominoids (20 million yr ago). The DRB3 locus appears to have arisen before the divergence of hominoids on an ancestral DRB1 lineage. The DRB2 and DRB5 loci were generated more than 20 million yr ago and the DRB4 locus more than 5 million yr ago. The DRB2 locus, a pseudogene in humans, is polymorphic in the nonhuman primates.